Next Article in Journal
High-Resolution Structure Analysis of Antibody V5 and U4 Conformational Epitopes on Human Papillomavirus 16
Previous Article in Journal
Herpesvirus and Autophagy: “All Right, Everybody Be Cool, This Is a Robbery!”
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 9
Issue 12
10.3390/v9120373
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genomic Characterisation of Vinegar Hill Virus, An Australian Nairovirus Isolated in 1983 from Argas Robertsi Ticks Collected from Cattle Egrets

by
Penelope J. Gauci
1,*,
Jane McAllister
1,
Ian R. Mitchell
1,
Daisy Cybinski
2,
Toby St George
2 and
Aneta J. Gubala
1
1
Land Division, Defence Science & Technology Group, Fishermans Bend, Victoria 3207, Australia
2
Formerly: Long Pocket Laboratories, Commonwealth Scientific and Industrial Research Organisation, Indooroopilly, Queensland 4068, Australia
*
Author to whom correspondence should be addressed.
Viruses 2017, 9(12), 373; https://doi.org/10.3390/v9120373
Submission received: 8 November 2017 / Revised: 27 November 2017 / Accepted: 28 November 2017 / Published: 5 December 2017
(This article belongs to the Section Animal Viruses)
Browse Figures
Versions Notes

Abstract

:
This report describes the near complete genomic sequence and subsequent analysis of Vinegar Hill virus (VINHV; tentative member of the genus Orthonairovirus, family Nairoviridae, order Bunyavirales). VINHV is the second nairovirus reported to be isolated on mainland Australia and the first to be sequenced and analysed. Our genetic analysis shows that VINHV belongs to the Dera Ghazi Khan genogroup, a group of viruses previously isolated in other parts of the world including Asia, South Africa, and the USA. We discuss possible routes of entry for nairoviruses into Australia and the need to understand the virome of Australian ticks in the context of new and emerging disease.

1. Introduction

The genus Orthonairovirus (family Nairoviridae, order Bunyavirales) comprises 12 species to which more than 60 predominantly tick-borne viruses have been assigned, including several associated with severe human and livestock disease such as Crimean-Congo haemorrhagic fever (CCHF) and Nairobi sheep disease (NSD) viruses, respectively [1,2,3]. The nairovirus genome consists of three negative-sense, single-stranded RNA (-ssRNA) segments; small (S), medium (M), and large (L), that encode the nucleoprotein (N protein), glycoprotein precursor (GPC), and the RNA-dependant RNA polymerase (L protein), respectively. Until recently, there was little genetic data available for viruses of this genus, and the available information was mostly restricted to viruses of the CCHF and NSD serogroups. However, with the recent publication of several full-length nairoviruses genomes, our knowledge is increasing and changing the way we relate members within this important genus [3,4,5]. Consequently, Walker et al. [5] proposed the assignment of nairoviruses into nine distinct genogroups. For the most part, the genogroups represent the corresponding established serogroups, with the exception of the NSD and CCHF serogroups, which are combined into a single genogroup, NSD.
Nairoviruses have rarely been isolated in Australia or its territories. The isolation of three Australian nairoviruses has previously been documented; two Sakhalin genogroup viruses (Taggert virus, TAGV; and Finch Creek virus, FCV) isolated from Ixodes uriae ticks on Macquarie Island 1500 km south-southeast of Tasmania [6,7], and a tentative Dera Ghazi Khan (DGK) genogroup virus (isolate NT15470) isolated from Argas robertsi ticks in the Northern Territory [8]. Whilst TAGV is the only Australian nairovirus that has been fully sequenced [5], a partial sequence of FCV is also available in GenBank. An additional nairovirus, isolate CSIRO1499, which we describe in this manuscript and propose the name Vinegar Hill virus (VINHV), was isolated on the Australian mainland in 1983 [9]. VINHV was isolated from a pool of five female A. robertsi ticks collected from a rookery of the cattle egret (Bulbulcus ibis) at Gatton, Queensland in December 1981. The ticks were collected during an investigation of a heavy mortality in nestling chicks in the rookery. Following collection, the ticks were held alive until their blood meal was fully digested, and were subsequently stored frozen for future virus isolation. In contrast to the high isolation rate usually experienced with bird tick viruses, VINHV was isolated from only one of several hundred pools of ticks processed. Although results are unavailable, historical information suggests that serological testing found VINHV to be unrelated to any other arboviruses in Australia [9]. It is not evident from the literature if this testing included isolate NT15470, which had been demonstrated to be a strain of, or closely related to, Kao Shuan virus of the DGK genogroup [8]. VINHV was subsequently sent for further testing to the former World Reference Centre for Arboviruses at Yale (now the World Reference Centre for Emerging Viruses and Arboviruses at Galveston, TX, USA). Antibodies to VINHV have been detected in 13 of 401 (3.2%) tested seabird sera and 1 of 101 (1%) tested human sera. VINHV has been shown to produce an antibody response in experimentally infected cattle egrets and is capable of causing death in birds; however, no data on the pathogenicity of this virus for mammalian hosts is available [10].
The threat of new and emerging disease, together with the ease and availability of next generation sequencing technologies, has seen an increase in the characterisation of new and historic viruses. Here we report the near complete genome sequence of VINHV and subsequent predictive genetic and amino acid analyses, and demonstrate that this virus is a tentative member of the family Nairoviridae, genus Orthonairovirus. This is the first report describing the sequence of a nairovirus from the Australian mainland.

2. Materials and Methods

2.1. Virus Culture and Genomic Sequencing

Vinegar Hill virus (VINHV; isolate CSIRO1499) was obtained from the Berrimah Veterinary Laboratories, Darwin, NT, Australia. Viruses were propagated in baby hamster kidney BHK-BSR cells (a derivative of the BHK-21 cell line) grown in supplemented Basal Medium Eagle (Gibco, Grand Island, NY, USA) at 37 °C, following which they were harvested, and total RNA was extracted and converted to double stranded cDNA as previously described [11]. The cDNA material was prepared for high-throughput sequencing using the TruSeq CHIP-seq (Illumina, San Diego, CA, USA) protocols and standard multiplex adaptors. A paired-end, 250-base-read protocol was used for sequencing on an Illumina MiSeq instrument at The Ramaciotti Centre for Genomics, University of New South Wales, Sydney, Australia. Primary assembly of raw data and generation of consensus sequences were performed as described previously [11].

2.2. Bioinformatic Analysis

Alignments and predictive open reading frame (ORF) analysis was performed using Clone Manager (SciEd, Denver, CO, USA). Analysis of deduced proteins and targeting predictions were generated using the Centre for Biological Sequence Analysis portal tools including SignalP 4.1, NetOGlyc 4.0, NetNGlyc 1.0, ProP 1.0, and TMHMM (http://www.cbs.dtu.dk/services/). Pairwise identities were determined using p-distance analysis in MEGA7 [12].

2.3. Phylogenetic Analysis

Phylogenetic trees were constructed using 38 complete L protein, GPC, and N protein sequences, of nairoviruses accessed from GenBank (listed in Table S1) and the cognate VINHV protein sequences. Sequences were aligned using the MUSCLE 3.6 algorithm [13]. Bayesian analyses of amino acid (aa) sequence alignments were performed with BEAST software [14], using a Whelan and Goldman (WAG) model of aa substitution with gamma + invariant site heterogeneity. A lognormal relaxed clock model was also used, with a tree prior set to coalescent:exponential growth. The model was run with a Markov-chain Monte Carlo (MCMC) chain length of 10,000,000 with the output logged every 1000 steps, producing 10,000 trees. The maximum clade credibility tree was chosen using Tree Annotator (1000 tree burn-in) and trees were created using FigTree v1.4 (http://tree.bio.ed.ac.uk/software/figtree/). Maximum likelihood (ML) trees were constructed using MEGA5 [15], employing the WAG model of aa substitution with a gamma distribution of rate variation and 1000 bootstrap replications.

3. Results and Discussion

3.1. VINHV Genome and Terminal Sequences

The complete S and L genomic segments of VINHV and near complete M genomic segment were sequenced using high throughput sequencing (GenBank Accession numbers MF176883, MF176881, and MF176882, respectively). The lengths of the S, M, and L segments are 1729 nucleotides (nt), 4473 nt (lacking the 3′-terminal region, genome sense), and 12,133 nt, respectively. The organisation of the VINHV genome is consistent with those found in other orthonairoviruses, each containing a single ORF encoding the N protein, GPC, and L protein (Table 1). BlastX analysis of the GenBank databases indicates that VINHV is most similar to Dera Ghazi Khan virus (DGKV) sharing 78%, 69% and 71% amino acid identity with the translated protein products from the L, M, and S segments, respectively.
A common feature of bunyaviruses is the conservation of genus-specific genome termini. In nairoviruses, the consensus terminal nt sequences are 3′ AGAGUUUCU- and 5′ UCUCAAAGA-. The genome termini of VINHV S and L segments are consistent with this, with the exception of a single nucleotide change (U→A) at position 9 of the 3′ terminus of the S segment. The terminal sequences of Dera Ghazi Khan genogroup viruses have been observed to differ at position 9 of both terminal ends for each segment, with the exception of DGKV (Figure S1) [3]. DGKV has a deviation from the consensus at position 9 of the 3′ termini only of the M and S segments, similar to the observation in VINHV, further supporting a close relationship between VINHV and DGKV. Attempts to obtain the 3′ terminal non-coding sequence of the VINHV M segment were unsuccessful, despite several attempts. As the coding sequence for the GPC was complete, allowing comparative analyses with other GPCs to be performed, further attempts to obtain the non-coding portion at the 3′ terminus were abandoned. It is anticipated, however, that given the overall similarity of VINHV to DGKV, a similar sequence would be present at the 3′ end, though this would need to be confirmed.

3.2. L Protein

The single ORF on the L segment of VINHV encodes a 3948-aa viral L polymerase protein. The L polymerases of -ssRNA viruses contain four conserved regions reflective of the universal functions of this protein [16,17,18]. Region I has a presumed cap-snatching endonuclease activity [19], whilst the function of region II is unknown. Region III, also called the polymerase module, contains six conserved motifs (pre-motif A and motifs A–E) and is predicted to be involved in catalytic functions of the polymerase, and in template and/or primer positioning [17]. Region IV is suggested to have a role in capped primer-cleavage and 5′ viral RNA binding [16]. All of these regions and motifs are highly conserved in the putative VINHV L protein (Figure 1a–d). Although zinc finger and leucine zipper sequence motifs have previously been identified in the CCHFV L protein [20], these are not always apparent in all nairoviruses [3]. Likewise, they are not apparent in the VINHV L protein.
In addition to these regions, an ovarian tumour (OTU)-like domain (pfam02338) has been identified in proximity to the N termini of the L protein of all nairoviruses with the possible exception of the “nairo-like” viruses including South Bay virus, which has a divergent sequence that shows some homology to the OTU-like domain [3,5,20,21]. Similarly, an OTU-like domain is predicted in the VINHV L protein (Figure 2). The observed functionality differences between the OTU domains of virulent CCHFV and less virulent DUGV lead some to speculate that this domain may be a virulence factor [3,22].

3.3. GPC

The VINHV M segment contains a single ORF that putatively encodes a 1414-aa polyprotein, which, like other bunyavirus M segment polyproteins, is predicted to be co- and post-translationally processed into mature viral glycoproteins [23]. The VINHV polyprotein shares similar sequence organisation to other nairoviruses and contains various conserved post-translational modification sites and structural features (Figure 3). The study of CCHFV provides much of our understanding of nairovirus GPC structure and processing [24,25,26,27,28]. The analysis of CCHFV shows that the GPC has an N-terminal mucin-like domain containing a large number of predicted O-glycosylation sites, followed by a protein of unknown function (GP38), an envelope glycoprotein (Gn), a non-structural protein (NSm), and a second envelope glycoprotein (Gc).
Similar to other nairoviruses, the VINHV GPC is predicted to contain an N-terminal signal peptide (at VLA30-NT) followed by a highly O-glycoslylated (16 sites) mucin-like domain, but one which is considerably shorter with less predicted O-glycosylation sites than in CCHFV. Although the function of the mucin-like domain of CCHFV GPC remains undetermined, a similar mucin-like domain in the Ebola virus glycoprotein GP1 is known to play a major role in pathogenesis [29]. The M segment is the most variable of the three segments, and this is particularly notable in the hypervariable N-terminal region that precedes the Gn protein. The characteristics of this region are generally genogroup-specific in relation to variation in the number of O-glycosylation sites and the length of the predicted mucin-like domain [5]. Viruses of the DGK genogroup are generally known to have one of the smallest mucin-like domains amongst the nairoviruses ranging between 56 to 124 aa in length, containing between seven and 22 O-glycosylation sites. The CCHFV mucin-like domain is cleaved by a furin or furin-like protease (site RSKR), generating a 247-aa protein [30]. There does not appear to be an equivalent furin-like protease cleavage site in any of the analysed DGK group GPCs, including in VINHV; however, there are a number of possible alternate protease cleavage sites in the vicinity of the domain (Figure 3). The VINHV M segment does not appear to encode an NSm protein, which is consistent with all other nairoviruses except for the NSD genogroup viruses, which do encode this protein [5].
The Gn and Gc glycoproteins of nairoviruses are relatively well conserved in size and structural characteristics [5]. Similar to predictions in other nairoviruses, the VINHV Gn and Gc proteins are predicted to be cleaved by the subtilisin/kexin-isozyme-1 (SKI-1) protease at sites RHLL383↓ and RRLL775↓, respectively. The VINHV Gn and Gc proteins are of similar size to those of other DGK group viruses, and they contain numerous conserved cysteine residues which have a functional role in protein folding, transmembrane domains, and zinc finger domains (Figure 3; Figures S2 and S3).
VINHV Gn contains three predicted glycosylation sites (Figure S2). The location of the first glycosylation site (NGTK432) is universally conserved amongst all nairoviruses. The second glycosylation site (NGSG498) appears to be conserved with DGKV and Sapphire II virus (SAPV), whilst the third site (NHTS509) appears to be unique to VINHV. Similarly, there are three predicted glycosylation sites in the VINHV Gc protein (Figure S3). The first (NNSV795) is conserved amongst all the analysed DGK viruses with the exception of SAPV, and the second (NGSI1151) is conserved in all the analysed DGK group viruses. The third site (NCTG1309) is generally conserved with most of the viruses in the orthonairovirus genus [5].

3.4. N Protein

The N protein of -ssRNA viruses binds to genomic RNA to form ribonucleoprotein complexes that associate with the polymerase for viral RNA synthesis (transcription and replication) and form the structural core of the virion [31]. The length of the VINHV N protein is 499 aa, which is similar to those of other nairoviruses. Crystal structure studies of the N protein of CCHFV demonstrated two major domains, a globular head and an extended stalk, with RNA/DNA binding-associated sites predominantly found on the head domain [31,32,33]. In comparative sequence analysis, the VINHV N protein exhibits conservation of these binding sites, either fully (K132, R134, K222, Q300, K343, R384, H453, and Q457), or with a conservative change (H197N, Y374H, E387D, and K411R) [31,32] (Figure S4). The caspase-3 cleavage site motif previously identified in some nairoviruses (CCHFV, Hazara virus (HAZV) and Thiafora genogroup viruses) is not apparent in the VINHV N protein [5].
Pairwise alignments show that the VINHV N protein shares 54.4 to 72% identity with the N proteins of other viruses within the DGK genogroup (Table 2) and 31.7 to 42.1% identity with the N proteins of representative viruses from the other genogroups (Table S2). Walker et al. [5] suggested a sequence identity cut-off of 52% to place viruses into genogroups. Using this criterion, the placement of VINHV into the Dera Ghazi Khan genogroup is well supported.

3.5. Phylogenetic Analysis

Until recently, only nairoviruses associated with hard ticks had been fully sequenced. This posed challenges for the phylogenetic analysis of nairoviruses associated with soft ticks, which were clearly different. Whilst partial L protein data for some soft tick nairoviruses has existed for some time [34], the recent work of Walker et al. [4,5] and Kuhn et al. [3] has produced sequence data enabling the comprehensive genomic analysis of numerous soft-tick viruses from this genus. It is evident from these studies that the phylogenetic relationship of nairoviruses broadly reflects vector preferences, genome organisation, and serological relationships.
Bayesian phylogenetic analyses were performed using the N protein, GPC, and L protein of VINHV and other representative nairoviruses (Figure 4a–c, respectively.) Maximum likelihood analyses were also performed for comparison and these produced trees with similar topologies (data not shown). The phylogenetic analyses demonstrate strong support for the formation of nine distinct clades representing the nine proposed genogroups [5]. Lower support present at some of the deeper nodes is most likely reflective of the divergence of some viruses, and will only be strengthened by the sequencing of additional viruses in this genus.
The inclusion of VINHV within the DGK genogroup is strongly supported and the relationships inferred within the group are uniform with all three segments. Viruses of the DGK genogroup are widespread throughout the world (Pakistan, Taiwan, Thailand, South Africa, USA, and Australia) and, in most instances, have been isolated from ticks feeding on birds. Thus, it is feasible that the distribution of this group of nairoviruses, including the introduction of VINHV and other nairoviruses to the Australian mainland, may be via avian migration. The cattle egret, which VINHV is associated with, has populated Australia only since the late 1940s.
Phylogenetic analyses also demonstrate that VINHV and DGKV share a common ancestor. DGKV was isolated from ticks feeding on camels in Pakistan in 1966. This may indicate another entry route for VINHV, or its ancestor, into Australia via the 10–20,000 camels that were brought into the country from India and Pakistan in the period of 1860–1907.

3.6. Ticks and Emerging Viruses in Australia

Though tick-borne diseases do not contribute greatly to the overall communicable disease burden in Australia, an increase in incidence may be seen in the future with climatic, population, and lifestyle changes [35]. Also, it is possible that a proportion of unknown or undiagnosed illnesses could be attributed to tick vectors. It is essential that we gain an understanding of the biome of the native ticks, particularly those that are known to bite humans. The A. robertsi tick, from which VINHV was isolated, is one of five soft tick species in Australia that possibly feed on humans and domestic animals [36]. However, the most important tick in Australia, from both a medical and veterinary perspective, is Ixodes holocyclus, a hard tick species. Hence, much of the tick research in Australia is focused on hard ticks and associated diseases, particularly of bacterial origin. I. holocyclus is the vector for Rickettsia australis and R. honei, the aetiological agents of the only two recognised tick-borne diseases in Australia—Queensland tick typhus and Flinders Island spotted fever, respectively [35,36]. It is speculated that this species of tick may also have a role in Hendra virus transmission [37]. Furthermore, amid debate regarding the presence of tick-borne Lyme disease in Australia, a Borrelia sp. related to the Lyme disease agent has been isolated from this species of tick [38]. However, despite patients presenting with Lyme-like disease, no aetiological agent has been linked to disease locally, and therefore the presence of Lyme disease in Australia is not confirmed.
Advances in sequencing technology have allowed us to investigate the biome of arthropods that are known vectors of disease. Although an ongoing study into tick-borne diseases in Australia has developed strategies to successfully identify low abundant bacteria in hard ticks [38], it is not evident whether this study will expand to include the identification of viral agents. Analysis of the viromes of three American ticks revealed a diverse array of viruses, including several novel viruses with genetic similarities to pathogens of humans and livestock [21]. Likewise, a similar metagenomics analysis of Australian mosquitoes detected the presence of viruses from families Flaviridae, Rhabdoviridae, Reoviridae, Togaviridae, and Bunyaviridae [39]. It is evident that there is great potential for novel and emerging viruses circulating in Australian arthropods.
Although some viruses have previously been isolated from Australian ticks [40,41], none as yet have been associated with human disease. However, it is important to note that antibodies to VINHV have been found in human sera and, as such, the potential threat to human health demands further investigation. With the ongoing sequencing of historic Australian virus isolates, our understanding of viruses circulating within Australia will also increase [4,5,11,41,42,43,44]. In other parts of the world, tick-borne infectious diseases are on the rise and becoming a serious world health problem affecting both human and animal health. For example, there has been a marked increase in the range and incidence of CCHF since 2000, and tick-borne encephalitis is a growing concern in Europe and Asia [45]. Similarly, the incursion of African swine fever into the Caucasus, and potentially from there into Europe, is of deep concern and requires preventative strategies to avoid the spread of this disease [46].
Bird migration plays an important part in the spread of tick-borne disease. For example, the massive expansion of the cattle egret range began after cattle became established in newly created, and expanding, cattle pastures in continents additional to Africa and Asia. They flew to the Americas in 1933, Australia 1948, and Europe in 1958, and have subsequently spread widely from there. Cattle were introduced to each of these new territories following European exploration and removal of forests. It is presumed that the associated ticks and tick-borne viruses have spread more slowly, as tick- and virus-free colonies exist within flight ranges of an infected colony. Thus, the pasture habitats are created for cattle, and the migrant egrets, ticks, and viruses follow, in that order. Investigations into the virome of Australian ticks will provide valuable information on the potential for the emergence of new viruses associated with this vector within the Australian landscape. Adequate biosurveillance in this area should be prioritised to mitigate any potential future emerging diseases.

Supplementary Materials

The following are available online at www.mdpi.com/1999-4915/9/12/373/s1. Figure S1: Genome terminal sequences of VINHV and other Dera Ghazi Khan genogroup viruses, Figure S2: Amino acid sequences alignment of glycoprotein Gn of the DGK group viruses and CCHFV, Figure S3: Amino acid sequence alignment of glycoprotein Gc of the DGK group viruses and CCHFV, Figure S4: Amino acid sequence alignment of the nucleoproteins of DGK group viruses and CCHFV, Table S1: Viruses and GenBank accession numbers of sequences used in phylogenetic and other comparative analyses.Table S2: Amino acid sequence identities (%) of nairovirus N proteins as determined by p-distance estimation in MEGA7.

Acknowledgments

We gratefully acknowledge the work of the scientists who studied and collected the ticks and sera from the infected rookeries, with specific thanks to Helen Zakrzewski and Steven Davis for their excellent technical assistance, and to Evan Harris for the great care he provided in looking after the cattle egrets experimentally infected with VINHV.

Author Contributions

Penelope J. Gauci, Aneta J. Gubala, and Toby St George conceived and designed the experiments; Penelope J. Gauci, Jane McAllister, and Daisy Cybinski performed the experiments; Penelope J. Gauci and Ian R. Mitchell analyzed the data; Penelope J. Gauci, Jane McAllister, Ian R. Mitchell, Daisy Cybinski, Toby St George, and Aneta J. Gubala contributed reagents/materials/analysis tools; Penelope J. Gauci and Aneta J. Gubala wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Plyusnin, A.; Beaty, B.J.; Elliott, R.M.; Goldbach, R.; Kormelink, R.; Lundkvist, A.; Schmaljohn, C.S.; Tesh, R.B. Bunyaviridae. In Virus Taxonomy: Ninth Report of the International Committee on Taxonomy of Viruses; King, A.M.Q., Adams, M.J., Carstens, E.B., Lefkowitz, E.J., Eds.; Elsevier Academic Press: San Diego, CA, USA, 2012; pp. 725–741. [Google Scholar]
  2. Adams, M.J.; Lefkowitz, E.J.; King, A.M.Q.; Harrach, B.; Harrison, R.L.; Knowles, N.J.; Kropinski, A.M.; Krupovic, M.; Kuhn, J.H.; Mushegian, A.R.; et al. Changes to taxonomy and the International Code of Virus Classification and Nomenclature ratified by the International Committee on Taxonomy of Viruses (2017). Arch. Virol. 2017, 162, 2505–2538. [Google Scholar] [CrossRef] [PubMed]
  3. Kuhn, J.H.; Wiley, M.R.; Rodriguez, S.E.; Yīmíng, B.; Prieto, K.; Travassos da Rosa, A.P.A.; Guzman, H.; Savji, N.; Ladner, J.T.; Tesh, R.B.; et al. Genomic Characterization of the Genus Nairovirus (Family Bunyaviridae). Viruses 2016, 8, 164. [Google Scholar] [CrossRef] [PubMed]
  4. Walker, P.J.; Widen, S.G.; Firth, C.; Blasdell, K.R.; Wood, T.G.; Travassos Da Rosa, A.P.A.; Guzman, H.; Tesh, R.B.; Vasilakis, N. Genomic characterization of yogue, kasokero, issyk-kul, keterah, gossas, and thiafora viruses: Nairoviruses naturally infecting bats, shrews, and ticks. Am. J. Trop. Med. Hyg. 2015, 93, 1041–1051. [Google Scholar] [CrossRef] [PubMed]
  5. Walker, P.J.; Widen, S.G.; Wood, T.G.; Guzman, H.; Tesh, R.B.; Vasilakis, N. A Global genomic characterization of nairoviruses identifies nine discrete genogroups with distinctive structural characteristics and host-vector associations. Am. J. Trop. Med. Hyg. 2016, 94, 1107–1122. [Google Scholar] [CrossRef] [PubMed]
  6. Doherty, R.L.; Carley, J.G.; Murray, M.D.; Main, A.J., Jr.; Kay, B.H.; Domrow, R. Isolation of arboviruses (Kemerovo group, Sakhalin group) from Ixodes uriae collected at Macquarie Island, Southern ocean. Am. J. Trop. Med. Hyg. 1975, 24, 521–526. [Google Scholar] [CrossRef] [PubMed]
  7. Major, L.; La Linn, M.; Slade, R.W.; Schroder, W.A.; Hyatt, A.D.; Gardner, J.; Cowley, J.; Suhrbier, A. Ticks Associated with macquarie island penguins carry arboviruses from four genera. PLoS ONE 2009, 4, e4375. [Google Scholar] [CrossRef] [PubMed]
  8. Doherty, R.L.; Kay, B.H.; Carley, J.G.; Filippich, C. Isolation of virus strains related to kao shuan virus from argas robertsi in northern territory, Australia. Search 1976, 7, 484. [Google Scholar]
  9. St George, T.D. Research on Milk Fever, Bluetongue, Akabane, Ephemeral Fever and Other Arboviruses at CSIRO Long Pocket Laboratories, Brisbane, Australia. Ph.D. Thesis, University of Queensland, Brisbane, Australia, 2011. [Google Scholar]
  10. Humphery-Smith, I.; Cybinski, D.H.; Byrnes, K.A.; St George, T.D. Seroepidemiology of arboviruses among seabirds and island residents of the Great Barrier Reef and Coral Sea. Epidemiol. Infect. 1991, 107, 435–440. [Google Scholar] [CrossRef] [PubMed]
  11. Gauci, P.J.; McAllister, J.; Mitchell, I.R.; Boyle, D.B.; Bulach, D.M.; Weir, R.P.; Melville, L.F.; Gubala, A.J. Genomic characterisation of three mapputta group viruses, a serogroup of Australian and papua new guinean bunyaviruses associated with human disease. PLoS ONE 2015, 10, e0116561. [Google Scholar] [CrossRef] [PubMed]
  12. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed]
  13. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [PubMed]
  14. Drummond, A.J.; Suchard, M.A.; Xie, D.; Rambaut, A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol. Biol. Evol. 2012, 29, 1969–1973. [Google Scholar] [CrossRef] [PubMed]
  15. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [PubMed]
  16. Aquino, V.H.; Moreli, M.L.; Moraes Figueiredo, L.T. Analysis of oropouche virus L protein amino acid sequence showed the presence of an additional conserved region that could harbour an important role for the polymerase activity. Arch. Virol. 2003, 148, 19–28. [Google Scholar] [CrossRef] [PubMed]
  17. Muller, R.; Poch, O.; Delarue, M.; Bishop, D.H.L.; Bouloy, M. Rift valley fever virus L segment: Correction of the sequence and possible functional role of newly identified regions conserved in RNA-dependent polymerases. J. Gen. Virol. 1994, 75, 1345–1352. [Google Scholar] [CrossRef] [PubMed]
  18. Poch, O.; Sauvaget, I.; Delarue, M.; Tordo, N. Identification of four conserved motifs among the RNA-dependent polymerase encoding elements. EMBO J. 1989, 8, 3867–3874. [Google Scholar] [PubMed]
  19. Reguera, J.; Weber, F.; Cusack, S. Bunyaviridae RNA polymerases (L-protein) have an N-terminal, influenza-like endonuclease domain, essential for viral cap-dependent transcription. PLoS Pathog. 2010, 6, e1001101. [Google Scholar] [CrossRef] [PubMed]
  20. Honig, J.E.; Osborne, J.C.; Nichol, S.T. Crimean–Congo hemorrhagic fever virus genome L RNA segment and encoded protein. Virology 2004, 321, 29–35. [Google Scholar] [CrossRef] [PubMed]
  21. Tokarz, R.; Williams, S.H.; Sameroff, S.; Sanchez Leon, M.; Jain, K.; Lipkin, W.I. Virome analysis of Amblyomma americanum, Dermacentor variabilis, and Ixodes scapularis ticks reveals novel highly divergent vertebrate and invertebrate viruses. J. Virol. 2014, 88, 11480–11492. [Google Scholar] [CrossRef] [PubMed]
  22. Bakshi, S.; Holzer, B.; Bridgen, A.; McMullan, G.; Quinn, D.G.; Baron, M.D. Dugbe virus ovarian tumour domain interferes with ubiquitin/ISG15-regulated innate immune cell signalling. J. Gen. Virol. 2013, 94, 298–307. [Google Scholar] [CrossRef] [PubMed]
  23. Clerx, J.P.M.; Casals, J.; Bishop, D.H.L. Structural characteristics of Nairoviruses (Genus Nairovirus, Bunyaviridae). J. Gen. Virol. 1981, 55, 165–178. [Google Scholar] [CrossRef] [PubMed]
  24. Lasecka, L.; Baron, M.D. The molecular biology of nairoviruses, an emerging group of tick-borne arboviruses. Arch. Virol. 2014, 159, 1249–1265. [Google Scholar] [CrossRef] [PubMed]
  25. Sanchez, A.J.; Vincent, M.J.; Nichol, S.T. Characterization of the glycoproteins of Crimean-Congo hemorrhagic fever virus. J. Virol. 2002, 76, 7263–7275. [Google Scholar] [CrossRef] [PubMed]
  26. Bergeron, É.; Vincent, M.J.; Nichol, S.T. Crimean-Congo hemorrhagic fever virus glycoprotein processing by the endoprotease SKI-1/S1P is critical for virus infectivity. J. Virol. 2007, 81, 13271–13276. [Google Scholar] [CrossRef] [PubMed]
  27. Vincent, M.J.; Sanchez, A.J.; Erickson, B.R.; Basak, A.; Chretien, M.; Seidah, N.G.; Nichol, S.T. Crimean-congo hemorrhagic fever virus glycoprotein proteolytic processing by subtilase SKI-1. J. Virol. 2003, 77, 8640–8649. [Google Scholar] [CrossRef] [PubMed]
  28. Altamura, L.A.; Bertolotti-Ciarlet, A.; Teigler, J.; Paragas, J.; Schmaljohn, C.S.; Doms, R.W. Identification of a novel C-terminal cleavage of Crimean-Congo hemorrhagic fever virus PreG(N) that leads to generation of an NSM protein. J. Virol. 2007, 81, 6632–6642. [Google Scholar] [CrossRef] [PubMed]
  29. Yang, Z.-Y.; Duckers, H.J.; Sullivan, N.J.; Sanchez, A.; Nabel, E.G.; Nabel, G.J. Identification of the Ebola virus glycoprotein as the main viral determinant of vascular cell cytotoxicity and injury. Nat. Med. 2000, 6, 886. [Google Scholar] [PubMed]
  30. Sanchez, A.J.; Vincent, M.J.; Erickson, B.R.; Nichol, S.T. Crimean-congo hemorrhagic fever virus glycoprotein precursor is cleaved by Furin-like and SKI-1 proteases to generate a novel 38-kilodalton glycoprotein. J. Virol. 2006, 80, 514–525. [Google Scholar] [CrossRef] [PubMed]
  31. Carter, S.D.; Surtees, R.; Walter, C.T.; Ariza, A.; Bergeron, É.; Nichol, S.T.; Hiscox, J.A.; Edwards, T.A.; Barr, J.N. Structure, function, and evolution of the Crimean-Congo hemorrhagic fever virus nucleocapsid protein. J. Virol. 2012, 86, 10914–10923. [Google Scholar] [CrossRef] [PubMed]
  32. Guo, Y.; Lou, Z.; Rao, Z.; Wang, W.; Ji, W.; Zhou, H.; Yang, C.; Deng, M.; Deng, F.; Wang, H.; et al. Crimean-Congo hemorrhagic fever virus nucleoprotein reveals endonuclease activity in bunyaviruses. Proc. Natl. Acad. Sci. USA 2012, 109, 5046–5051. [Google Scholar] [CrossRef] [PubMed]
  33. Wang, Y.; Dutta, S.; Karlberg, H.; Devignot, S.; Weber, F.; Hao, Q.; Tan, Y.J.; Mirazimi, A.; Kotaka, M. Structure of Crimean-Congo hemorrhagic fever virus nucleoprotein: Superhelical homo-oligomers and the role of caspase-3 cleavage. J. Virol. 2012, 86, 12294–12303. [Google Scholar] [CrossRef] [PubMed]
  34. Honig, J.E.; Osborne, J.C.; Nichol, S.T. The high genetic variation of viruses of the genus Nairovirus reflects the diversity of their predominant tick hosts. Virology 2004, 318, 10–16. [Google Scholar] [CrossRef] [PubMed]
  35. Lowbridge, C.P.; Doggett, S.L.; Graves, S. Bug Breakfast in the Bulletin. Tickborne diseases. N. S. W. Public Health Bull. 2011, 22, 237. [Google Scholar] [CrossRef] [PubMed]
  36. Barker, S.C.; Walker, A.R.; Campelo, D. A list of the 70 species of Australian ticks; diagnostic guides to and species accounts of Ixodes holocyclus (paralysis tick), Ixodes cornuatus (southern paralysis tick) and Rhipicephalus australis (Australian cattle tick); and consideration of the place of Australia in the evolution of ticks with comments on four controversial ideas. Int. J. Parasitol. 2014, 44, 941–953. [Google Scholar] [PubMed]
  37. Barker, S.C. The Australian paralysis tick may be the missing link in the transmission of Hendra virus from bats to horses to humans. Med. Hypotheses 2003, 60, 481–483. [Google Scholar] [CrossRef]
  38. Gofton, A.W.; Oskam, C.L.; Paparini, A.; Greay, T.L.; Ryan, U.; Irwin, P.; Lo, N.; Wei, H.; McCarl, V.; Beninati, T.; et al. Inhibition of the endosymbiont “Candidatus Midichloria mitochondrii” during 16S rRNA gene profiling reveals potential pathogens in Ixodes ticks from Australia. Parasites Vectors 2015, 8, 345. [Google Scholar] [CrossRef] [PubMed]
  39. Coffey, L.L.; Page, B.L.; Greninger, A.L.; Herring, B.L.; Russell, R.C.; Doggett, S.L.; Haniotis, J.; Wang, C.; Deng, X.; Delwart, E.L. Enhanced arbovirus surveillance with deep sequencing: Identification of novel rhabdoviruses and bunyaviruses in Australian mosquitoes. Virology 2014, 448, 146–158. [Google Scholar] [CrossRef] [PubMed]
  40. St. George, T.D.; Standeast, H.A.; Doherty, R.L.; Carley, J.G.; Fillipich, C.; Brandsma, J. The isolation of saumarez reef virus, a new flavivirus, from bird ticks ornithodoros capensis and ixodes eudyptidis in Australia. Aust. J. Exp. Biol. Med. 1977, 55, 493. [Google Scholar] [CrossRef]
  41. Gauci, P.J.; McAllister, J.; Mitchell, I.R.; George, T.D.S.; Cybinski, D.H.; Davis, S.S.; Gubala, A.J.; St George, T.D. Hunter island group phlebovirus in ticks, Australia. Emerg. Infect. Dis. 2015, 21, 2246–2248. [Google Scholar] [CrossRef] [PubMed]
  42. Bixing, H.; Firth, C.; Watterson, D.; Allcock, R.; Colmant, A.M.G.; Hobson-Peters, J.; Kirkland, P.; Hewitson, G.; McMahon, J.; Hall-Mendelin, S.; et al. Genetic characterization of archived bunyaviruses and their potential for emergence in Australia. Emerg. Infect. Dis. 2016, 22, 833–840. [Google Scholar]
  43. Huang, B.; Allcock, R.; Warrilow, D. Newly characterized arboviruses of northern Australia. Virol. Rep. 2016, 6, 11–17. [Google Scholar] [CrossRef]
  44. McAllister, J.; Gauci, P.J.; Mitchell, I.R.; Gubala, A.J.; Boyle, D.B.; Bulach, D.M.; Weir, R.P.; Melville, L.F.; Davis, S.S. Genomic characterisation of Almpiwar virus, Harrison Dam virus and Walkabout Creek virus; three novel rhabdoviruses from northern Australia. Virol. Rep. 2014, 3, 1–17. [Google Scholar] [CrossRef]
  45. Estrada-Peña, A.; de la Fuente, J. Review: The ecology of ticks and epidemiology of tick-borne viral diseases. Antivir. Res. 2014, 108, 104–128. [Google Scholar] [CrossRef] [PubMed]
  46. Manzano-Román, R.; Díaz-Martín, V.; de la Fuente, J.; Pérez-Sánchez, R. Soft Ticks as Pathogen Vectors: Distribution, Surveillance and Control, Parasitology; Shah, M.M., Ed.; InTech: Rijeka, Croatia, 2012. [Google Scholar]
Viruses 09 00373 g001 550
Figure 1. Conserved motifs on Vinegar Hill virus (VINHV) L protein. Amino acid sequence alignments of RNA-dependent RNA polymerase modules of VINHV and other Dera Ghazi Khan group viruses, with Crimean Congo haemorrhagic virus included as a reference. (a) Region I and (b) Region II are highly conserved in bunya- and arenaviruses; (c) Region III, including pre-motif A and motifs A–E, is highly conserved in negative sense RNA viruses, and (d) Region IV is highly conserved in bunyaviruses and other segmented negative sense RNA viruses. Conserved amino acid residues amongst the DGK group of viruses and CCHFV are shaded green.
Figure 1. Conserved motifs on Vinegar Hill virus (VINHV) L protein. Amino acid sequence alignments of RNA-dependent RNA polymerase modules of VINHV and other Dera Ghazi Khan group viruses, with Crimean Congo haemorrhagic virus included as a reference. (a) Region I and (b) Region II are highly conserved in bunya- and arenaviruses; (c) Region III, including pre-motif A and motifs A–E, is highly conserved in negative sense RNA viruses, and (d) Region IV is highly conserved in bunyaviruses and other segmented negative sense RNA viruses. Conserved amino acid residues amongst the DGK group of viruses and CCHFV are shaded green.
Viruses 09 00373 g001
Viruses 09 00373 g002 550
Figure 2. Ovarian tumour (OTU)-like domain on VINHV L protein. Amino acid sequence alignment of a fragment of the L protein of VINHV and other DGK group viruses depicting a predicted (OTU)-like protease domain close to the N-terminal region. Amino acids that constitute the domain are marked with * below the sequence. Amino acids conserved between the DGK group viruses and CCHFV are shaded green.
Figure 2. Ovarian tumour (OTU)-like domain on VINHV L protein. Amino acid sequence alignment of a fragment of the L protein of VINHV and other DGK group viruses depicting a predicted (OTU)-like protease domain close to the N-terminal region. Amino acids that constitute the domain are marked with * below the sequence. Amino acids conserved between the DGK group viruses and CCHFV are shaded green.
Viruses 09 00373 g002
Viruses 09 00373 g003 550
Figure 3. Illustration of the glycoprotein precursors (GPC) of VINHV and other DGK group nairoviruses. The CCHFV GPC is included as a reference for comparison. Regions corresponding to the mucin-like domain (light blue), GP38 (tan), Gn (orange), NSm (grey), and Gc (dark blue) are depicted. The numbers of predicted O-linked glycosylation sites in the mucin-like domain are shown. Signal peptidase cleavage sites (SP), potential furin-like (for CCHFV only) and subtilisin/kexin-isozyme-1 (SKI-I) cleavage sites that might be utilised to generate the individual proteins, were predicted using SignalP 4.1 and comparative studies and are shown above each schematic. A potential alternate cleavage site (rrll) identified in the VINHV protein is depicted in lower case. Transmembrane domains are indicated by yellow boxes. Additional cleavage sites predicted by ProP1.0 are indicated below the schematic with a green stripe (scores > 0.5, higher confidence) or a pink stripe (scores 0.3–0.49, lower confidence). DGKV: Dera Ghazi Khan virus; SAPV: Sapphire II virus; AMV: Abu Mina virus; TUNV: Tunis virus; AHV: Abu Hammad virus.
Figure 3. Illustration of the glycoprotein precursors (GPC) of VINHV and other DGK group nairoviruses. The CCHFV GPC is included as a reference for comparison. Regions corresponding to the mucin-like domain (light blue), GP38 (tan), Gn (orange), NSm (grey), and Gc (dark blue) are depicted. The numbers of predicted O-linked glycosylation sites in the mucin-like domain are shown. Signal peptidase cleavage sites (SP), potential furin-like (for CCHFV only) and subtilisin/kexin-isozyme-1 (SKI-I) cleavage sites that might be utilised to generate the individual proteins, were predicted using SignalP 4.1 and comparative studies and are shown above each schematic. A potential alternate cleavage site (rrll) identified in the VINHV protein is depicted in lower case. Transmembrane domains are indicated by yellow boxes. Additional cleavage sites predicted by ProP1.0 are indicated below the schematic with a green stripe (scores > 0.5, higher confidence) or a pink stripe (scores 0.3–0.49, lower confidence). DGKV: Dera Ghazi Khan virus; SAPV: Sapphire II virus; AMV: Abu Mina virus; TUNV: Tunis virus; AHV: Abu Hammad virus.
Viruses 09 00373 g003
Viruses 09 00373 g004 550
Figure 4. Phylogenetic analyses of the (a) nucleoprotein; (b) glycoprotein precursor, and (c) L protein of VINHV and other representative nairoviruses. Relationships were inferred by Bayesian analysis of the protein sequence alignments. A Whelan and Goldman (WAG) model of amino acid substitution with gamma + invariant site heterogeneity was used. Numbers represent Bayesian posterior probabilities. SSV1 was included as an outgroup. Trees are drawn to scale measured in substitutions per site as indicated by the scale bar.
Figure 4. Phylogenetic analyses of the (a) nucleoprotein; (b) glycoprotein precursor, and (c) L protein of VINHV and other representative nairoviruses. Relationships were inferred by Bayesian analysis of the protein sequence alignments. A Whelan and Goldman (WAG) model of amino acid substitution with gamma + invariant site heterogeneity was used. Numbers represent Bayesian posterior probabilities. SSV1 was included as an outgroup. Trees are drawn to scale measured in substitutions per site as indicated by the scale bar.
Viruses 09 00373 g004
Table
Table 1. Genome size comparisons of Vinegar Hill Virus (VINHV) and other Dera Ghazi Khan (DGK) group viruses (anti-genome sense).
Table 1. Genome size comparisons of Vinegar Hill Virus (VINHV) and other Dera Ghazi Khan (DGK) group viruses (anti-genome sense).
SegmentRegionLength (nt/aa) for Indicated Viruses
VINHVDGKVAMVAHVTUNVSAPVa CCHFV
L5′UTR53574135359676
L ORF11,847/394811,847/394811,925/397411,868/395511,868/395511,871/395611,837/3945
3′UTR233106249214400218194
segment total12,13312,01012,21512,11712,30312,18512,108
M5′UTR41 b421430301492
GPC ORF4245/14144239/14124203/14004200/13994167/13884341/14465055/1684
3′UTR187159224472609434219
segment total4473444044414702480647895366
S5′UTR58505154565455
NP ORF1500/4991500/4991497/4981497/4981497/4981476/4911448/482
3′UTR171221212226315137168
segment total1729177117601777186816681672
a CCHFV is included as a reference nairovirus. b Incomplete. VINHV: Vinegar Hill virus; DGK: Dera Ghazi Khan; UTR: Untranslated region; ORF: open reading frame; L: RNA dependent RNA polymerase; GPC: Glycoprotein precursor; NP: Nucleoprotein; DGKV: Dera Ghazi Khan virus; AMV: Abu Mina virus; AHV: Abu Hammad virus; TUNV: Tunis virus; SAPV: Sapphire II virus; CCHFV: Crimean Congo haemorrhagic virus.
Table
Table 2. Amino acid sequence identities (%) of the nucleoproteins of Dera Ghazi Khan genogroup viruses. Sequence identities were determined by p-distance estimation using MEGA7.
Table 2. Amino acid sequence identities (%) of the nucleoproteins of Dera Ghazi Khan genogroup viruses. Sequence identities were determined by p-distance estimation using MEGA7.
VirusVINHVDGKVAHVTUNVAMV
VINHV
DGKV72.0
AHV62.159.3
TUNV62.159.995.4
AMV60.663.476.075.8
SAPV54.452.952.452.452.4

Share and Cite

MDPI and ACS Style

Gauci, P.J.; McAllister, J.; Mitchell, I.R.; Cybinski, D.; St George, T.; Gubala, A.J. Genomic Characterisation of Vinegar Hill Virus, An Australian Nairovirus Isolated in 1983 from Argas Robertsi Ticks Collected from Cattle Egrets. Viruses 2017, 9, 373. https://doi.org/10.3390/v9120373

AMA Style

Gauci PJ, McAllister J, Mitchell IR, Cybinski D, St George T, Gubala AJ. Genomic Characterisation of Vinegar Hill Virus, An Australian Nairovirus Isolated in 1983 from Argas Robertsi Ticks Collected from Cattle Egrets. Viruses. 2017; 9(12):373. https://doi.org/10.3390/v9120373

Chicago/Turabian Style

Gauci, Penelope J., Jane McAllister, Ian R. Mitchell, Daisy Cybinski, Toby St George, and Aneta J. Gubala. 2017. "Genomic Characterisation of Vinegar Hill Virus, An Australian Nairovirus Isolated in 1983 from Argas Robertsi Ticks Collected from Cattle Egrets" Viruses 9, no. 12: 373. https://doi.org/10.3390/v9120373

APA Style

Gauci, P. J., McAllister, J., Mitchell, I. R., Cybinski, D., St George, T., & Gubala, A. J. (2017). Genomic Characterisation of Vinegar Hill Virus, An Australian Nairovirus Isolated in 1983 from Argas Robertsi Ticks Collected from Cattle Egrets. Viruses, 9(12), 373. https://doi.org/10.3390/v9120373

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop