Figures
Abstract
Viruses belonging to the species Wallal virus and Warrego virus of the genus Orbivirus were identified as causative agents of blindness in marsupials in Australia during 1994/5. Recent comparisons of nucleotide (nt) and amino acid (aa) sequences have provided a basis for the grouping and classification of orbivirus isolates. However, full-genome sequence data are not available for representatives of all Orbivirus species. We report full-genome sequence data for three additional orbiviruses: Wallal virus (WALV); Mudjinabarry virus (MUDV) and Warrego virus (WARV). Comparisons of conserved polymerase (Pol), sub-core-shell ‘T2’ and core-surface ‘T13’ proteins show that these viruses group with other Culicoides borne orbiviruses, clustering with Eubenangee virus (EUBV), another orbivirus infecting marsupials. WARV shares <70% aa identity in all three conserved proteins (Pol, T2 and T13) with other orbiviruses, consistent with its classification within a distinct Orbivirus species. Although WALV and MUDV share <72.86%/67.93% aa/nt identity with other orbiviruses in Pol, T2 and T13, they share >99%/90% aa/nt identities with each other (consistent with membership of the same virus species - Wallal virus). However, WALV and MUDV share <68% aa identity in their larger outer capsid protein VP2(OC1), consistent with membership of different serotypes within the species - WALV-1 and WALV-2 respectively.
Citation: Belaganahalli MN, Maan S, Maan NS, Pritchard I, Kirkland PD, Brownlie J, et al. (2014) Full Genome Characterization of the Culicoides-Borne Marsupial Orbiviruses: Wallal Virus, Mudjinbarry Virus and Warrego Viruses. PLoS ONE 9(10): e108379. https://doi.org/10.1371/journal.pone.0108379
Editor: Robert Belshaw, Plymouth University, United Kingdom
Received: May 21, 2014; Accepted: August 19, 2014; Published: October 9, 2014
Copyright: © 2014 Belaganahalli et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. Full length sequences of Seg-1 to Seg-10 of WALV (AUS1978/09), MUDV (AUS1982/03) and WARV (AUS1969/01) have been submitted to GenBank with accession numbers KJ495745 to KJ495754; KJ495765 to KJ495774 and KJ495755 to KJ495764 respectively. These data were compared to other orbiviruses (accession numbers for the data used are given in Table S1 in File S1).
Funding: This work was partly funded by the Commonwealth Scholarship Commission (http://cscuk.dfid.gov.uk/) awarded to MB, partly by the Department for Environment, Food and Rural Affairs (Defra)(SE2617) awarded to SM. NSM, PCM and HA were funded by the Biotechnology and Biological Sciences Research Council (BBSRC), the European Union (EU) (SANCO/ 940/2002; EMIDA grant OrbiNet - K1303206) and WildTech - Grant no. 222633-2. IP was part of this work, funded by CSIRO & Dept of Agriculture, Fisheries and Forestry, Australian Government. PK was funded by NSW Department of Primary Industries. No additional external funding was received for this study. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: Houssam Attoui is a member of the PLOS ONE editorial board. This does not alter the authors' adherence to PLOS ONE Editorial policies and criteria.
Introduction
The orbiviruses are classified within the genus Orbivirus, family Reoviridae. They have icosahedral and non-enveloped virus-particles, containing a genome composed of ten linear pieces of dsRNA. The orbiviruses currently include 22 virus species (which also represent distinct ‘serogroups’) that have been recognised by the International Committee for the Taxonomy of Viruses (ICTV), as well as 14 unclassified isolates (some of which may represent additional species) [1]–[6]. The orbiviruses are primarily transmitted by ticks or haematophagous-insect-vectors (including mosquitoes, Culicoides [biting midges] and sand flies) and collectively have a wide host-range that includes domestic and wild ruminants, equines, marsupials, sloths, bats, birds and humans [2], [7], [8].
Orbiviruses, including Wallal virus (WALV), Warrego virus (WARV) and Eubenangee virus (EUBV) have previously been isolated from marsupials, and there is serological evidence for infections by Wongorr virus (WGRV) [7]. Evidence for a causal link between the orbiviruses and clinical disease in marsupials was provided by the isolation of WALV and WARV from animals affected by ‘kangaroo blindness’ in Australia during 1994/95 [9], [10]. Both viruses were isolated from Western Grey kangaroos (Macropus fuliginosis), red kangaroos (Macropus rufa), Eastern grey kangaroos (Macropus giganteus) and euros or wallaroo (Macropus robustus). Many thousands of kangaroos were affected in the states of New South Wales, South Australia and Victoria, suffering from non-suppurative choroiditis, retinal degeneration with inflammation and secondary degeneration of the optic nerve [9], [10]. A histopathology slide from a kangaroo with chorio-retinitis seen in the early 1970s, was found to have the same ‘structure’ as cases in 1994/95, suggesting periodic outbreaks of the disease [11]. Later in 1998–99, viruses belonging to the species Eubenangee virus were also isolated from wallabies affected by ‘Tammar Wallaby Sudden Death Syndrome’ (TSDS) [12], [13].
The species Wallal virus contains three distinct serotypes/strains: Wallal virus (WALV), Wallal K virus (WALKV) and Mudjinabarry virus (MUDV) [14]. WALV was first isolated in 1969 from adult Culicoides collected in Charleville, Queensland, Australia [15]. MUDV was isolated in September 1971 from Culicoides marksi collected in the Alligator Rivers area of the Northern Territory of Australia [16].
The species Warrego virus contains three distinct serotypes/strains: Warrego virus (WARV), Warrego K virus (WARKV) and Mitchell River virus (MRV). WARV virus was first isolated in 1970 from adult Culicoides collected in Charleville, Queensland, Australia [7], [15]. Mitchell River virus was isolated in April 1969 from Culicoides spp. trapped in the Mitchell River region in Australia [15].
The species Eubenangee virus contains four distinct viruses: Eubenangee virus (EUBV), Tilligerry virus (TILV), Pata virus (PATAV) and Ngoupe virus (NGOV) [2]. EUBV and TILV were isolated in Australia in 1963 and 1971 respectively from pools of mosquitoes [7]. They were also isolated from Culicoides spp. and antibodies to these viruses have been detected in cattle and marsupials [7], [17]. PATAV and NGOV were isolated in 1968 and 1974 respectively in the Central African Republic, also from pools of mosquitoes [7]. Nucleotide (nt) and amino-acid (aa) sequence comparisons indicate that PATAV represents a novel and distinct Orbivirus species [18].
The different species of Culicoides borne orbivirus (CBO) include Bluetongue virus (BTV) (the type species of the genus Orbivirus), African horse sickness virus (AHSV) and Epizootic haemorrhagic disease virus (EHDV), which are thought to include the most economically important orbiviruses. Complete genomes have been sequenced for representative isolates of BTV, EHDV, AHSV as well as some other CBO species including Palyam virus (PALV) and Equine encephalosis virus (EEV) (the members of which infect domesticated livestock) [19]–[23]. Although WALV, WARV and EUBV are also transmitted by adult Culicoides, [9], [12], [24] the full genome sequence of only one of these viruses – EUBV – has been reported [18].
Recent advances in RT-PCR and high throughput DNA sequencing methods have generated full-genome sequence data for multiple orbivirus isolates, helping to clarify the genetic relatedness and taxonomic status within and between different orbivirus species, including BTV, AHSV and EHDV [18], [19], [25]–[28]. Comparisons of nucleotide (nt) and amino acid (aa) sequence data now provide a basis for molecular epidemiology studies and for grouping distinct orbivirus isolates into serotypes and topotypes [19], [26], [29], [30].
As part of a programme to generate representative data for all Orbivirus species, we report full genome nucleotide sequences for WALV, MUDV and WARV, and present a comparison to isolates from other Orbivirus species as an aid to help identify novel orbivirus isolates. The full-genome sequences of WALV and MUDV indicate that they represent different serotypes within the Wallal virus species and that WARV represents a distinct species within the genus Orbivirus.
Materials and Methods
Viruses
The viruses used in this study [AUS1978/09 (WALV), AUS1982/03 (MUDV) and AUS1969/01 (WARV)] were obtained from the Orbivirus Reference Collection (ORC) at The Pirbright Institute. These viruses were originally taken from naturally infected animals by qualified veterinarians, as part of normal diagnostic testing procedures in the respective countries. MUDV was propagated in BHK-21 cells (clone 13 obtained from European Collection of Animal cell Cultures (ECACC – 84100501), while WALV and WARV were grown in BSR cells (a clone of BHK) [35] or BHK cells, in Dulbecco's minimum essential medium (DMEM) supplemented with antibiotics (100 units/ml penicillin and 100 µg/ml streptomycin) and 2 mM glutamine. Infected cell cultures were incubated until they showed widespread (100%) cytopathic effects (CPE). Viruses were harvested, aliquoted and used for the dsRNA extraction, or stored in the orbivirus reference collection (ORC) at −80°C.
Preparation of viral dsRNA
Intact genomic dsRNA was extracted from WALV, MUDV and WARV infected cell cultures using a guanidinium isothiocyanate extraction procedure as described by Attoui et al [36]. Briefly, the infected cell pellet was lysed in 1 ml of commercially available TRIZOL reagent (Invitrogen), 0.2 volume of chloroform was added, mixed by vortexing and the mixture incubated on ice for 10 min. The supernatant containing total RNA was separated from cellular debris and DNA by centrifuging at 10,000 xg for 10 min at 4°C. Single stranded RNA (ssRNA) was removed by 2 M LiCl precipitation at 4°C overnight, followed by centrifugation at 10,000 xg for 5 min. An equal volume of isopropanol containing 750 mM ammonium acetate was added to the supernatant, then mixed and the viral dsRNA was allowed to precipitate for a minimum of 2 hours at −20°C. The dsRNA was pelleted by centrifugation at 10,000 xg for 10 min, washed with 70% ethanol, air dried and suspended in nuclease free water (NFW). The RNA was either used immediately or stored at −20°C.
Reverse transcription of dsRNA and PCR amplification of cDNAs
The genome segments of WALV, MUDV and WARV were reverse-transcribed using a ‘full-length amplification of cDNA' (FLAC) technique described by Maan et al [37]. Briefly, a 35 base self-priming oligonucleotide ‘anchor-primer’ with a phosphorylated 5′ terminus was ligated to the 3′ ends of the viral dsRNAs using the T4 RNA ligase, followed by reverse transcription using RT system (Promega). The resulting cDNAs were amplified using complementary primers to the anchor primer. The resulting cDNA amplicons were analyzed by agarose gel electrophoresis. For cloning purposes, a high fidelity KOD polymerase enzyme (Novagen) was used in the PCR.
Cloning and sequencing of cDNAs
Amplified cDNAs were purified and cloned. MUDV genes were cloned into the Strataclone blunt-end PCR cloning vector ‘pSC-B-amp/kan’ supplied in the StrataClone Blunt PCR cloning kit (Stratagene). WALV (only Seg-2) and WARV amplicons were cloned into the ‘pCR-Blunt’ vector supplied with the Zero Blunt PCR Cloning Kit (Invitrogen). Recombinant plasmid-vectors containing MUDV inserts were transformed into Solopack competent cells (Agilent Technologies), while WALV and WARV plasmids were transformed into One Shot TOP10 competent cells supplied with the respective cloning kits. Clones containing relevant inserts were identified by touch PCR using M13 universal primers. Plasmids were extracted from the clones identified using the QIAprep Spin MiniPrep Kit (Qiagen). The plasmids and PCR products were sequenced using an automated ABI 3730 DNA sequencer (Applied Biosystems). Some of the MUDV sequencing primers detected other segments of WALV and using the sequences generated, running primers were designed to obtain full sequence for the remaining nine segments of WALV.
Sequence analysis and phylogenetic tree construction
‘Raw’ ABI sequence data were assembled into ‘contigs’ using the SeqManII sequence analysis package (DNAstar version 5.0). The ORFs of WALV, MUDV and WARV were identified and translated into aa sequences for further analysis using EditSeq (DNAstar version 5.0). The putative function of each protein was identified by BlastX comparisons to homologous orbivirus (e.g. BTV) proteins in GenBank (http://blast.ncbi.nlm.nih.gov/Blast.cgi?CMD=Web&PAGE_TYPE=BlastHome). Multiple alignments of consensus sequences were performed using ClustalX (Version 2.0) [38], Clustal Omega (http://www.ebi.ac.uk/Tools/msa/clustalo/) and MAFFT [39] to ensure proper alignment. Aligned protein sequences were back translated to nucleotide sequences using DAMBE [40]) or RevTrans 1.4 server available online (http://www.cbs.dtu.dk/services/RevTrans/) for further nucleotide analysis. The best fit amino acid (aa) and nucleotide (nt) models for Maximum likelihood (ML) analysis were determined using ProtTest 3.0 and jModeltest respectively [41], [42]. The models were also determined using MEGA 5 software. The consensus or simplest model given by Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC) was selected for ML tree construction. The aa model rtREV (G+F) was selected for T2 tree and rtREV (I+G+F) model was selected for both VP1 and T13 proteins. One nt model, GTR (+G) was used to construct ML trees for all three genes. All phylogenetic trees construction and pairwise distance calculations using p-distance parameter were performed using MEGA 5 [43], [44]. GenBank nucleotide accession numbers for the sequences used for analysis and phylogenetic studies are listed in the Table S1 in File S1.
Results
Virus propagation and genomic-segment migration-pattern (electropherotype)
WALV (isolate - AUS1978/09), MUDV (isolate AUS1982/03) and WARV (isolate AUS1969/01) all induced severe cytopathic effects (CPE) in BHK and BSR (a clone of BHK cells) cell monolayers at 48–72 hours post infection. When analysed by 1% agarose gel electrophoresis (AGE), the genomic dsRNA of these viruses migrated as distinct groups, consisting of 3 large segments, 3 medium size segments, 3 small segments, and 1 very small segment (giving a 3-3-3-1 pattern) (Figure 1). This pattern is similar to that observed in other CBOs (including BTV, EHDV and EUBV) [18], [37]. WALV and MUDV generated identical genome segment migration patterns (‘electropherotype’) during AGE, as expected for viruses belonging to the same virus species (Figure 1) [2], [37]. In contrast, WARV showed clear differences in the migration of Seg-3, -7, -8 and -9 (Lane 3 in Figure 1) confirming that it belongs to a different virus species.
1 = WALV (AUS1978/09); 2 = MUDV (AUS1982/03); 3 = WARV (AUS1969/01).
Genome sequence analyses of WALV (AUS1978/09), MUDV (AUS1982/03) and WARV (AUS1969/01)
Full length sequences of Seg-1 to Seg-10 of WALV (AUS1978/09), MUDV (AUS1982/03) and WARV (AUS1969/01) have been submitted to GenBank with accession numbers KJ495745 to KJ495754; KJ495765 to KJ495774 and KJ495755 to KJ495764 respectively. These data were compared to other orbiviruses (accession numbers for the data used are given in Table S1 in File S1). The total genome size of WALV, MUDV and WARV is 19,387 bp, 19,382 bp and 19,428 bp respectively, with segments ranging from 3959 bp to 862 bp (3964 bp to 838 bp in WARV), encoding proteins of 1306 aa to 87 aa (1303 aa to 78 aa in WARV) (Table 1 and 2). As suggested by the electropherotype analysis (Figure 1), most of the individual genome segments of WALV and MUDV, with the exception of Seg-2 (VP2) and Seg-8 (NS2(ViP)) have identical sizes with similar ORFs, identical start and stop positions generating proteins of the same length (Table 1).
The different genome segments of WALV and MUDV have the same conserved termini, sharing five fully conserved nucleotides at their 5′ and 3′ends (+ve: 5′-GUAU/AA………..AC/AA/GCU/A/CUAC-3′). WARV shares six fully conserved nucleotides at the 5′ ends and four at the 3′ ends (+ve: 5′-GUAAAA………..AC/AA/U/CUAC-3′). Two terminal nucleotides at the 5′ end and three nucleotides at the 3′ end (5′-GU……UAC-3′) of all three viruses are characteristic of the genus Orbivirus, and the first and last two nucleotides represent inverted complements (http://www.reoviridae.org/dsRNA_virus_proteins/CPV-RNA-Termin.htm).
The G+C content of the WALV, MUDV and WARV genomes are 39.14%, 39.41% and 38.99% respectively, well within the range previously identified for insect-borne orbiviruses (IBOs) (36.66% in PHSV to 45.86% in EEV) but less than tick-borne orbiviruses (TBOs) (51.93% in St Croix River virus (SCRV) to 57.29% in Great Island virus (GIV)) [5].
Like other orbiviruses, most genome segments of WALV, MUDV and WARV have shorter 5' than 3' terminal non-coding regions (NCRs). However, Seg-8 (encoding the viral inclusion body protein NS2) of WALV and MUDV, and Seg-5 (encoding the viral tubule protein NS1) of WARV have longer 5′ NCRs (Table 1). Exceptions were also observed in Umatilla virus (UMAV), Corriparta virus (CORV) and GIV, although the significance of these variations is uncertain [5]. Collectively the upstream and downstream terminal NCRs represent 4.29%, 4.3% and 4.08% of the WALV, MUDV and WARV genomes respectively. Comparisons with other CBOs showed similar or higher percentages, while the mosquito-borne orbiviruses (MBOs) and TBOs show higher values [5].
The coding assignments of the genome segments of WALV and MUDV (Table 1) are similar to those of BTV [45], [46]. The coding assignments of WARV are also similar, with the exception of Seg-7 and Seg-8 (coding for ViP and VP7 respectively), which have reversed their relative order of migration/size (Table 2). Most of the genome segments are monocistronic, encoding proteins from single open reading frames (ORF). However, Seg-10 of all three viruses has two in-frame initiation sites, at or near the start of the same ORF leading to translation of NS3 and NS3a (Table 1). Recently, a novel non-structural protein ‘NS4’ (which varies between 10 kDa to 22.5 kDa, depending on Orbivirus species), encoded by a downstream ORF (+2 frame) on Seg-9 (which also encodes VP6-Helicase) was identified in BTV and several other orbiviruses [33], [47], [48]. Similarly WALV, MUDV and WARV have a downstream ORF (+2 frame) on Seg-9 that encodes NS4.
Most of the WALV and MUDV ORF's have a ‘strong’ Kozak sequence (RNNAUGG), except for the first ORFs of Seg-9 and Seg-10, which have ‘weak’ (UAGAUGU) or ‘moderate’ (GUCAUGU) contexts respectively, potentially allowing read-through of the upstream initiation codon, enhancing expression of NS4 and NS3a respectively. Similarly, most ORF's of the WARV genome-segments have a strong (RNNAUGG) to moderate Kozak sequence [Seg-4 (AUCAUGC) and Seg-5 (UCUAUGA)], except Seg-9 (UGGAUGU) and Seg-10 (GUCAUGC), both of which have a weak ‘context’.
Weak or moderate Kozak sequences have previously been observed in several of the genome segments of other orbiviruses, although they still appear to be translated effectively [18], [33]. The down-stream NS4 ORF of Seg-9 in WALV and MUDV has a strong Kozak context (AAGAUGG) encoding an 87 aa protein, and in WARV has a moderate Kozak sequence (GGGAUGC) encoding 78 aa protein. In all three viruses WALV, MUDV and WARV, the second initiation site (for NS3a) in Seg-10 has the strong Kozak context (AA/CGAUGG).
Phylogenetic relationships of Seg-1/VP1(Pol)
The conserved nature of the orbivirus RNA dependent RNA polymerase (Pol) (encoded by Seg-1) has previously been used to help classify orbivirus isolates at the genus and species level based on sequence comparisons and identity levels [18], [25], [27], [31]. Phylogenetic comparisons of Seg-1/VP1(Pol) from WALV, MUDV and WARV with other orbiviruses (Figure 2), showed three groups containing the CBOs, MBOs and TBOs. In each case WALV, MUDV and WARV grouped with the CBOs.
One thousand bootstrap replicates were used for construction of the phylogenetic trees. The WALV, MUDV and WARV isolates characterised in this study are marked with green circle. BTV = Bluetongue virus; AHSV = African Horse sickness virus; EHDV = Epizootic Haemorrhagic Disease Virus; CHUV = Chuzan virus; EUBV = Eubenangee virus; TILV = Tilligerry virus; PATAV = Pata virus; EEV = Equine encephalosis virus; WARV = Warrego virus; WALV = Wallal virus; MUDV = Mudjinbarry virus; CORV = Corriparta virus; UMAV = Umatilla virus; SLOV = Stretch Lagoon orbivirus; SVIV = Sathuvachari virus; YUOV = Yunnan orbivirus; PHSV = Peruvian horse sickness virus; TRBV = Tribec virus; TRBVh = another isolate of Tribec virus; LIPV = Lipovnik virus; LIPV = Lipovnik virus isolate CzArLip 91; KEMV = Kemerovo virus; KEMVh = Kemerovo virus isolate EgAn 1169-61; GIV = Great Island virus; SCRV = St Croix river virus. Virus isolates and accession numbers of polymerase sequences used for comparative analysis are listed in Table S1 in File S1. ‘e’ and ‘w’ after serotype number indicate eastern and western strains, respectively.
Seg-1/VP1(Pol) of WALV and MUDV share 99.16/93.57 aa/nt identity, confirming that they belong to the same Orbivirus species. In comparisons to other orbiviruses they show a maximum of 70.96%/66.97% aa/nt identity with EUBV and a minimum of 36.08%/43.55% with SCRV (Table 3). VP1(Pol) of WARV shared up to a maximum of 65.41%/63.72% aa/nt identities with EHDV-2w, indicating its distinctiveness from other Orbivirus species within the CBOs.
Phylogenetic relationships of the subcore-shell ‘T2’ protein
BlastX analyses confirmed that the subcore shell T2 protein - VP3 - is encoded by Seg-3 of WALV, MUDV and WARV. Phylogenetic trees constructed for the T2 protein of different orbiviruses (Figure 3 and Figure S1 in File S1) showed that the CBOs (including BTV, AHSV, EHDV, WALV, EUBV, WARV and PALV) form a single group, in which T2(VP3) is encoded by Seg-3. The TBOs (GIV) and the MBOs (CORV, WGRV, Peruvian horse sickness virus [PHSV] and Yunnan orbivirus [YUOV]), form two further groups, in both of which the subcore T2(VP2) is encoded by Seg-2. The T2 protein (encoded by Seg-2) of SCRV branches separately, representing a more distantly related orbivirus from ticks.
One thousand bootstrap replicates were used for construction of the phylogenetic trees. The WALV, MUDV and WARV isolates characterised in this study are marked with green circle. WGRV = Wongorr virus; DAGV = D'Aguilar virus; PICV = Picola virus; PARV = Paroo river virus; CMPV = California mosquito pool virus; MPOV = Middle point orbivirus; BRDV = Broadhaven virus. Further abbreviations used in this figure are same as mentioned in Figure 2. Virus isolates and accession numbers of polymerase sequences used for comparative analysis are listed in Table S1 in File S1. ‘e’ and ‘w’ after serotype number indicate eastern and western strains, respectively.
The T2 protein/gene of WALV and MUDV clustered together, sharing 99.78%/90.36% aa/nt identities (Table 3), confirming that they belong to the same Orbivirus species (Figure 3). Comparisons to other orbiviruses showed a maximum of 72.86%/67.89% aa/nt identities with EUBV, down to 23.49%/37.44% with SCRV supporting classification of WALV and MUDV within a separate virus species. In contrast, WARV branched separately from other orbiviruses, sharing maximum aa/nt identities of 67.33%/65.55% with EUBV/WALV and a minimum of 24.4%/37.17% with SCRV (Table 3) confirming its classification within a distinct species.
Phylogenetic comparisons of Seg-7/VP7(T13)
The core surface protein VP7(T13) is the most abundant of the orbivirus structural proteins and is a strongly immuno-dominant serogroup/species-specific antigen [49]. The ML trees constructed for VP7(T13) protein and gene sequences of WALV, MUDV, WARV and other orbiviruses (Figure 4) exhibited similar topology to the T2 and VP1(Pol) trees (Figures 2 and 3) with three distinct ‘vector-groups’ (CBOs, MBOs and TBOs). WALV and MUDV again grouped together, while WARV branched separately within the CBOs. When compared to other orbiviruses, the VP7(T13) proteins/genes of WALV and MUDV (which share 100%/97.71% aa/nt identity) showed a maximum of 58.57%/60.67% aa/nt identity with EUBV and TILV. In contrast WARV again showed maximum aa/nt identity with EHDV (56.45%/59.45%). WALV and WARV share 39.94% to 58.57% aa sequence identity with the other CBOs; 20.06% to 25.21% aa identity with MBOs; 21.04% to 25.21% aa identity with the TBOs; but only 19.53% to 20.41% aa identity with the more distantly related SCRV (Table 3).
One thousand bootstrap replicates were used for construction of the phylogenetic trees. The WALV, MUDV and WARV isolates characterised in this study are marked with green circle. Abbreviations used in this figure are same as mentioned in Figure 2 and 3. Virus isolates and accession numbers of polymerase sequences used for comparative analysis are listed in Table S1 in File S1. ‘e’ and ‘w’ after serotype number indicate eastern and western strains, respectively.
Phylogenetic comparisons of Seg-2/VP2(OC1)
Orbivirus outer capsid protein (OC1) determines virus serotype and is highly variable, both in sequence identity levels and in size. VP2(OC1) of the CBOs is encoded by Seg-2 (VP2), by Seg-3 (VP3) of the MBOs and by Seg-4 (VP4) of the TBOs [5], [18], [33]. In WALV, MUDV and WARV viruses, VP2(OC1) is encoded by Seg-2, consistent with Culicoides transmission. An unrooted NJ phylogenetic tree constructed for OC1 of different orbiviruses (Figure 5) identified three major clusters, again corresponding to their ‘vector groups’ in a manner similar to the trees for VP1, T2 and T13 proteins (Figures 2, 3 and 4). These groupings suggest that despite a high level of aa sequence variation in response to immunological pressure from the vertebrate host (different serotypes), the orbiviruses have maintained a significant level of conservation in OC1 that is related to vector group. Despite much higher identity levels in their conserved genes and proteins, WALV and MUDV share only 67.89%/69.22% aa/nt identity in Seg-2/VP2(OC1), confirming that they represent different serotypes (WALV-1 and WALV-2 respectively) within the same Orbivirus species.
One thousand bootstrap replicates were used for construction of the phylogenetic trees. The VP2(OC1) protein of WALV, MUDV and WARV isolates characterised in this study are marked with green circle. CBOs = Culicoides borne orbiviruses; MBOs = Mosquito borne orbiviruses; TBOs = Tick borne orbiviruses. Abbreviations used for viruses in this figure are same as mentioned in Figure 2 and 3. Virus isolates and accession numbers of polymerase sequences used for comparative analysis are listed in Table S1 in File S1. ‘e’ and ‘w’ after serotype number indicate eastern and western strains, respectively.
Discussion
The majority of the existing Orbivirus species were initially recognised as distinct ‘serogroups’, sharing common antigens detectable by complement-fixation (CF), agar gel immuno-diffusion (AGID), fluorescent antibody (FA) tests, or enzyme-linked immune-sorbent assays (ELISA) [50], [51]. However, these serological assays are relatively insensitive, can lack specificity and can show significant levels of cross-reactivity between closely related Orbivirus species [52]–[54]. Identification of new orbivirus isolates (by species and serotype) was also dependent on the availability of reference antisera for the virus species that have already been recognised.
The advent of more rapid and reliable sequencing methods for dsRNA virus genomes has generated data for reference strains of the individual virus species, that can be easily disseminated for phylogenetic comparisons [5], [19], [34], [37]. These data sets now represent a primary tool for the identification and classification of novel orbivirus isolates [25]–[27], [31], [55].
The orbivirus sub-core-shell ‘T2’, polymerase ‘Pol’ and outer-core ‘T13’ proteins are all highly conserved, with >83%, >73% and >73% aa identity respectively within a single Orbivirus species. Consequently, these proteins and their comparisons have been used as ‘markers’ for identification and classification of existing and novel orbivirus isolates [5], [18], [26], [31], [32]. The highest inter-species identity levels that have been detected in these conserved proteins are between the closely related and intensively studied species: Bluetongue virus and Epizootic haemorrhagic disease virus (with maximum aa identities of 80%, 73%, and 66% in T2, Pol and T13 protein respectively). The genes encoding these proteins have also served as targets for development of orbivirus genus-specific and species-specific RT-PCR assays that can be used for diagnosis and virus discovery [18], [29], [56]–[58].
Nucleotide variation of up to 31.6% (equivalent to >68.4% nt identity) has been observed in the serotype controlling VP2(OC1) gene between different BTV isolates within a single serotype, with a maximum of 71.5% nt identity between different types providing useful ‘bench marks’ for typing of isolates in other orbivirus species [32], [34]
Phylogenetic analysis of the three most conserved VP1, T2 and T13 proteins/genes and the highly variable OC1 protein of the orbiviruses, that are presented here show three distinct ‘clusters’ that correspond to the arthropod vectors used by each virus (Figure 2, 3 and 4), indicating that the orbiviruses have evolved through ‘co-speciation’ with their arthropod vectors [18], [33]. The TBOs and MBOs form distinct phylogenetic clusters originating from a common branch but are more closely related to each other than to the CBOs (Figure 4 and 5). The TBOs appear to provide an ancestral ‘root’ for the insect transmitted orbiviruses. The identification of a smaller OC1 proteins in all of the MBOs (when compared to the CBOs) and a much smaller protein (∼50%) in the TBOs suggests that significant evolutionary changes in OC1 may have played a significant role in the initial vector preference and evolutionary separation of these groups of viruses. The changes in size of the OC1 protein of the IBOs compared to the TBOs could be related to an historic concatemerization event [59] doubling the size of the protein during the evolution of the MBOs and CBOs.
Phylogenetic trees for the different genome segments/proteins of WALV, MUDV and WARV confirm that members of the same Orbivirus species cluster more closely together, while members of distinct species branch separately. The trees constructed for orbivirus Pol, T2 and T13 genes are almost entirely congruent, also showing similarities to the OC1 tree, although with higher levels of identity within each virus species. This suggests that functional interactions between the different viral proteins/RNAs may restrict the exchange of genome segments (reassortment) between virus species.
WALV and MUDV are indistinguishable in complement fixation tests and can reassort with high frequency [60], demonstrating that they belong to the same Orbivirus species - Wallal virus. The data presented here, showing that WALV and MUDV share >99%/90% aa/nt identities in their conserved Pol, T2 and T13 proteins support this conclusion (Table 3). However, WALV and MUDV can be distinguished in serum neutralization tests [60], [61] and are therefore regarded as distinct serotypes - WALV-1 and WALV-2 respectively. This is supported by the data presented here, showing that they share <68%/70% aa/nt identity in outer capsid protein VP2(OC1).
The utility of sequence data for phylogenetic analysis of viral pathogens, their classification, as well as diagnostic test and vaccine development, suggests that full-genome sequencing is becoming an accepted standard for molecular epidemiological studies. It will therefore be useful to generate a full-genome sequence database for all orbivirus isolates, to help elucidate relationships, determine the frequency and significance of gene exchange (reassortment), the mechanisms involved in virus movement and emergence, and their impact on disease in vertebrate populations. The generation of full-genome sequences for reference strains of Wallal virus and Warrego virus, completes a genome data set for reference strains of the Culicoides borne Orbivirus species, including BTV, AHSV, EHDV, EUBV, PALV and EEV. The sequences generated in this study will not only help to identify novel Wallal and Warrego viruses, but will also provide a useful tool for identification and study of other orbiviruses, particularly other CBOs and those affecting marsupials.
Full genome references sequences are now available for fourteen of the twenty two Orbivirus species recognised by ICTV. These data will provide a basis for development of relevant diagnostic assays, and support molecular epidemiology/evolutionary studies, to promote our understanding of orbivirus disease in marsupials and other host groups.
Supporting Information
File S1.
Table S1 and Figure S1. Table S1. Nucleotide accession numbers for sequences used in phylogenetic analysis. Figure S1. Unrooted Neighbour-joining phylogenetic tree comparing orbivirus T2 proteins.
https://doi.org/10.1371/journal.pone.0108379.s001
(DOCX)
Acknowledgments
The authors would like to thank all collaborators for providing virus isolates for tests and analyses. P. Mertens is a Jenner Investigator.
Author Contributions
Conceived and designed the experiments: MB SM NSM HA PCM. Performed the experiments: MB SM NSM. Analyzed the data: MB SM HA NSM PCM. Contributed reagents/materials/analysis tools: MB SM HA PCM IP PK JB. Wrote the paper: MB SM HA NSM PCM.
References
- 1. Brown SE, Gorman BM, Tesh RB, Knudson DL (1992) Isolation of bluetongue and epizootic hemorrhagic disease viruses from mosquitoes collected in Indonesia. Vet Microbiol 32: 241–252.
- 2.
Mertens PPC, Maan S, Samuel A, Attoui H (2005) Orbiviruses, Reoviridae. In: Fauquet CM, M.AMayo, JManiloff, UDesselberger and L.ABall, editor. Virus Taxonomy Eighth Report of the International Committee on Taxonomy of Viruses. London: Elsevier/Academic Press. pp. 466–483.
- 3. Martins LC, Diniz JA, Silva EV, Barros VL, Monteiro HA, et al. (2007) Characterization of Minacu virus (Reoviridae: Orbivirus) and pathological changes in experimentally infected newborn mice. Int J Exp Pathol 88: 63–73.
- 4. Vieira Cde M, Nunes MR, da Silva EV, Carvalho VL, Nunes Neto JP, et al. (2009) Full-length sequencing and genetic characterization of Breu Branco virus (Reoviridae, Orbivirus) and two related strains isolated from Anopheles mosquitoes. J Gen Virol 90: 2183–2190.
- 5. Belaganahalli MN, Maan S, Maan NS, Nomikou K, Guimera M, et al. (2013) Full genome sequencing of Corriparta virus, identifies California mosquito pool virus as a member of the Corriparta virus species. PLoS One 8: e70779.
- 6. Kapoor A, Tesh RB, Duraisamy R, Popov VL, Travassos da Rosa AP, et al. (2013) A novel mosquito-borne Orbivirus species found in South-east Asia. J Gen Virol 94: 1051–1057.
- 7.
Karabatsos N (1985) International catalogue of arboviruses including certain other viruses of vertebrates. San Antonio: American Society of Tropical Medicine and Hygiene.
- 8.
Attoui H, Maan SS, Anthony SJ, Mertens PPC (2009) Bluetongue virus, other orbiviruses and other reoviruses: Their relationships and taxonomy. In: Mellor PS, Baylis M, Mertens PPC, editors. Bluetongue monograph. 1 ed. London: Elsevier/Academic Press. pp. 23–552.
- 9. Hooper PT, Lunt RA, Gould AR, Hyatt AD, Russell GM, et al. (1999) Epidemic of blindness in kangaroos–evidence of a viral aetiology. Aust Vet J 77: 529–536.
- 10. Reddacliff L, Kirkland P, Philbey A, Davis R, Vogelnest L, et al. (1999) Experimental reproduction of viral chorioretinitis in kangaroos. Aust Vet J 77: 522–528.
- 11. Curran G, Reddacliff L, Menke L, Witte I, Segboer L, et al. (1997) Investigating the epidemiology of a new disease: viral blindness in kangaroos. Epidemiol Sante Anim 31–32: 1–3.
- 12. Rose K, Kirkland P, Davis R, Cooper D, Blumstein D, et al. (2012) Epizootics of sudden death in tammar wallabies (Macropus eugenii) associated with an orbivirus infection. Aust Vet J 90
- 13.
Rose KA, Kirkland PD, Davies RJ, Cooper DW, Blumstein DT (2000) An epizootic of sudden death in Tammar Wallabies (Macropus eugenii). Proceedings AAZV and IAAAM joint conference. pp. 166–170.
- 14.
Attoui H, Mertens PPC, Becnel J, Belaganahalli M, Bergoin M, et al., editors (2012) Reoviridae. London, UK: Elsevier Academic Press. 497–650 p.
- 15. Doherty RL, Carley JG, Standfast HA, Dyce AL, Kay BH, et al. (1973) Isolation of arboviruses from mosquitoes, biting midges, sandflies and vertebrates collected in Queensland, 1969 and 1970. Trans R Soc Trop Med Hyg 67: 536–543.
- 16. Doherty RL, Standfast HA, Dyce AL, Carley JG, Gorman BM, et al. (1978) Mudjinbarry virus, an orbivirus related to Wallal virus isolated from midges from the northern territory of Australia. Aust J Biol Sci 31: 97–103.
- 17. Doherty RL, Whitehead RH, Wetters EJ, Gorman BM, Carley JG (1970) A survey of antibody to 10 arboviruses (Koongol Group, Mapputta group and ungrouped) isolated in Queensland. Trans R Soc Trop Med Hyg 64: 748–753.
- 18. Belaganahalli MN, Maan S, Maan NS, Nomikou K, Pritchard I, et al. (2012) Full genome sequencing and genetic characterization of Eubenangee viruses identify Pata virus as a distinct species within the genus Orbivirus. PLoS One 7: e31911.
- 19. Maan S, Maan NS, Ross-smith N, Batten CA, Shaw AE, et al. (2008) Sequence analysis of bluetongue virus serotype 8 from the Netherlands 2006 and comparison to other European strains. Virology 377: 308–318.
- 20. Anthony SJ, Maan N, Maan S, Sutton G, Attoui H, et al. (2009) Genetic and phylogenetic analysis of the non-structural proteins NS1, NS2 and NS3 of epizootic haemorrhagic disease virus (EHDV). Virus Res 145: 211–219.
- 21. Anthony SJ, Maan N, Maan S, Sutton G, Attoui H, et al. (2009) Genetic and phylogenetic analysis of the core proteins VP1, VP3, VP4, VP6 and VP7 of epizootic haemorrhagic disease virus (EHDV). Virus Res 145: 187–199.
- 22. Potgieter AC, Page NA, Liebenberg J, Wright IM, Landt O, et al. (2009) Improved strategies for sequence-independent amplification and sequencing of viral double-stranded RNA genomes. J Gen Virol 90: 1423–1432.
- 23. Yamakawa M, Kubo M, Furuuchi S (1999) Molecular analysis of the genome of Chuzan virus, a member of the Palyam serogroup viruses, and its phylogenetic relationships to other orbiviruses. J Gen Virol 80 (Pt 4): 937–941.
- 24. Hooper PT, Russell GM, Selleck PW, Lunt RA, Morrissy CJ, et al. (1999) Immunohistochemistry in the identification of a number of new diseases in Australia. Vet Microbiol 68: 89–93.
- 25. Belaganahalli MN, Maan S, Maan NS, Tesh R, Attoui H, et al. (2011) Umatilla virus genome sequencing and phylogenetic analysis: identification of stretch lagoon orbivirus as a new member of the Umatilla virus species. PLoS One 6: e23605.
- 26. Maan S, Maan NS, Nomikou K, Batten C, Antony F, et al. (2011) Novel bluetongue virus serotype from Kuwait. Emerg Infect Dis 17: 886–889.
- 27. Attoui H, Mohd Jaafar F, Belhouchet M, Aldrovandi N, Tao S, et al. (2005) Yunnan orbivirus, a new orbivirus species isolated from Culex tritaeniorhynchus mosquitoes in China. J Gen Virol 86: 3409–3417.
- 28. Anthony SJ, Maan S, Maan N, Kgosana L, Bachanek-Bankowska K, et al. (2009) Genetic and phylogenetic analysis of the outer-coat proteins VP2 and VP5 of epizootic haemorrhagic disease virus (EHDV): comparison of genetic and serological data to characterise the EHDV serogroup. Virus Res 145: 200–210.
- 29. Pritchard LI, Gould AR, Wilson WC, Thompson L, Mertens PP, et al. (1995) Complete nucleotide sequence of RNA segment 3 of bluetongue virus serotype 2 (Ona-A). Phylogenetic analyses reveal the probable origin and relationship with other orbiviruses. Virus Res 35: 247–261.
- 30. Nomikou K, Dovas CI, Maan S, Anthony SJ, Samuel AR, et al. (2009) Evolution and phylogenetic analysis of full-length VP3 genes of Eastern Mediterranean bluetongue virus isolates. PLoS One 4: e6437.
- 31. Attoui H, Stirling JM, Munderloh UG, Billoir F, Brookes SM, et al. (2001) Complete sequence characterization of the genome of the St Croix River virus, a new orbivirus isolated from cells of Ixodes scapularis. J Gen Virol 82: 795–804.
- 32. Maan S, Maan NS, van Rijn PA, van Gennip RG, Sanders A, et al. (2010) Full genome characterisation of bluetongue virus serotype 6 from the Netherlands 2008 and comparison to other field and vaccine strains. PLoS One 5: e10323.
- 33. Belhouchet M, Mohd Jaafar F, Tesh R, Grimes J, Maan S, et al. (2010) Complete sequence of Great Island virus and comparison with the T2 and outer-capsid proteins of Kemerovo, Lipovnik and Tribec viruses (genus Orbivirus, family Reoviridae). J Gen Virol 91: 2985–2993.
- 34. Maan S, Maan NS, Samuel AR, Rao S, Attoui H, et al. (2007) Analysis and phylogenetic comparisons of full-length VP2 genes of the 24 bluetongue virus serotypes. J Gen Virol 88: 621–630.
- 35. Sato M, Maeda N, Yoshida H, Urade M, Saito S (1977) Plaque formation of herpes virus hominis type 2 and rubella virus in variants isolated from the colonies of BHK21/WI-2 cells formed in soft agar. Arch Virol 53: 269–273.
- 36. Attoui H, Billoir F, Cantaloube JF, Biagini P, de Micco P, et al. (2000) Strategies for the sequence determination of viral dsRNA genomes. J Virol Methods 89: 147–158.
- 37. Maan S, Rao S, Maan NS, Anthony SJ, Attoui H, et al. (2007) Rapid cDNA synthesis and sequencing techniques for the genetic study of bluetongue and other dsRNA viruses. J Virol Methods 143: 132–139.
- 38. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, et al. (2007) Clustal W and Clustal X version 2.0. Bioinformatics 23: 2947–2948.
- 39. Katoh K, Asimenos G, Toh H (2009) Multiple alignment of DNA sequences with MAFFT. Methods Mol Biol 537: 39–64.
- 40. Xia X, Xie Z (2001) DAMBE: software package for data analysis in molecular biology and evolution. J Hered 92: 371–373.
- 41. Abascal F, Zardoya R, Posada D (2005) ProtTest: selection of best-fit models of protein evolution. Bioinformatics 21: 2104–2105.
- 42. Posada D (2008) jModelTest: phylogenetic model averaging. Mol Biol Evol 25: 1253–1256.
- 43. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, et al. (2011) MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 28: 2731–2739.
- 44. Saitou N, Nei M (1987) The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 4: 406–425.
- 45. Pedley S, Mohamed ME, Mertens PP (1988) Analysis of genome segments from six different isolates of bluetongue virus using RNA-RNA hybridisation: a generalised coding assignment for bluetongue viruses. Virus Res 10: 381–390.
- 46. Mertens PP, Brown F, Sangar DV (1984) Assignment of the genome segments of bluetongue virus type 1 to the proteins which they encode. Virology 135: 207–217.
- 47. Firth AE (2008) Bioinformatic analysis suggests that the Orbivirus VP6 cistron encodes an overlapping gene. Virol J 5: 48.
- 48. Belhouchet M, Mohd Jaafar F, Firth AE, Grimes JM, Mertens PP, et al. (2011) Detection of a fourth orbivirus non-structural protein. PLoS One 6: e25697.
- 49. Gumm ID, Newman JF (1982) The preparation of purified bluetongue virus group antigen for use as a diagnostic reagent. Arch Virol 72: 83–93.
- 50.
Verwoerd DW, Huismans H, Erasmus BJ (1979) Orbiviruses. In: Fraenkel-Conrat H, Wagner RR, editors. Comprehensive Virology. New York: Plenum Press. pp. 285–345.
- 51.
Gorman BM, Taylor J, Walker PJ (1983) Orbiviruses. In: Joklik WK, editor. The Reoviridae. New York and London: Plenum Press. pp. 287–357.
- 52. Moore DL (1974) Bluetongue and related viruses in Ibadan, Nigeria: serologic comparison of bluetongue, epizootic hemorrhagic disease of deer, and Abadina (Palyam) viral isolates. Am J Vet Res 35: 1109–1113.
- 53. Della-Porta AJ, Parsonson IM, McPhee DA (1985) Problems in the interpretation of diagnostic tests due to cross-reactions between orbiviruses and broad serological responses in animals. Prog Clin Biol Res 178: 445–453.
- 54. Borden EC, Shope RE, Murphy FA (1971) Physicochemical and morphological relationships of some arthropod-borne viruses to bluetongue virus–a new taxonomic group. Physiocochemical and serological studies. J Gen Virol 13: 261–271.
- 55. Attoui H, Mendez-Lopez MR, Rao S, Hurtado-Alendes A, Lizaraso-Caparo F, et al. (2009) Peruvian horse sickness virus and Yunnan orbivirus, isolated from vertebrates and mosquitoes in Peru and Australia. Virology 394: 298–310.
- 56. Shaw AE, Monaghan P, Alpar HO, Anthony S, Darpel KE, et al. (2007) Development and initial evaluation of a real-time RT-PCR assay to detect bluetongue virus genome segment 1. J Virol Methods 145: 115–126.
- 57. Mertens PP, Maan NS, Prasad G, Samuel AR, Shaw AE, et al. (2007) Design of primers and use of RT-PCR assays for typing European bluetongue virus isolates: differentiation of field and vaccine strains. J Gen Virol 88: 2811–2823.
- 58. Sakamoto K, Punyahotra R, Mizukoshi N, Ueda S, Imagawa H, et al. (1994) Rapid detection of African horsesickness virus by the reverse transcriptase polymerase chain reaction (RT-PCR) using the amplimer for segment 3 (VP3 gene). Arch Virol 136: 87–97.
- 59. Anthony SJ, Darpel KE, Belaganahalli MN, Maan N, Nomikou K, et al. (2011) RNA segment 9 exists as a duplex concatemer in an Australian strain of epizootic haemorrhagic disease virus (EHDV): Genetic analysis and evidence for the presence of concatemers as a normal feature of orbivirus replication. Virology 420: 164–171.
- 60. Gorman BM, Taylor J, Walker PJ, Young PR (1978) The isolation of recombinants between related orbiviruses. J Gen Virol 41: 333–342.
- 61. Walker PJ, Mansbridge JN, Gorman BM (1980) Genetic analysis of orbiviruses by using RNase T(1) oligonucleotide fingerprints. J Virol 34: 583–591.