Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

The taxonomy of an Australian nodavirus isolated from mosquitoes

  • David Warrilow ,

    Roles Formal analysis, Investigation, Methodology, Writing – original draft

    David.Warrilow@health.qld.gov.au (DW); j.peters2@uq.edu.au (JHP)

    Affiliation Public Health Virology Laboratory, Queensland Health Forensic and Scientific Services, Archerfield, Queensland

  • Bixing Huang,

    Roles Data curation, Investigation, Methodology

    Affiliation Public Health Virology Laboratory, Queensland Health Forensic and Scientific Services, Archerfield, Queensland

  • Natalee D. Newton,

    Roles Investigation, Methodology

    Affiliation Australian Infectious Diseases Research Centre, School of Chemistry and Molecular Biosciences, The University of Queensland, St Lucia, Queensland, Australia

  • Jessica J. Harrison,

    Roles Investigation, Methodology

    Affiliation Australian Infectious Diseases Research Centre, School of Chemistry and Molecular Biosciences, The University of Queensland, St Lucia, Queensland, Australia

  • Karyn N. Johnson,

    Roles Formal analysis, Writing – review & editing

    Affiliation School of Biological Sciences, The University of Queensland, St Lucia, Queensland, Australia

  • Weng Kong Chow,

    Roles Investigation, Methodology, Resources

    Affiliation Australian Defence Force Malaria and Infectious Disease Institute, Queensland, Australia

  • Roy A. Hall,

    Roles Conceptualization, Funding acquisition, Resources, Supervision, Writing – review & editing

    Affiliation Australian Infectious Diseases Research Centre, School of Chemistry and Molecular Biosciences, The University of Queensland, St Lucia, Queensland, Australia

  • Jody Hobson-Peters

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Resources, Supervision, Writing – original draft

    David.Warrilow@health.qld.gov.au (DW); j.peters2@uq.edu.au (JHP)

    Affiliation Australian Infectious Diseases Research Centre, School of Chemistry and Molecular Biosciences, The University of Queensland, St Lucia, Queensland, Australia

Abstract

We describe a virus isolated from Culex annulirostris mosquitoes in Australia. Phylogenetic analysis of its RNA-dependent RNA polymerase sequence and that of other related viruses revealed 6 clades, two of which corresponded wholly or partly with existing genera in the family Nodaviridae. There was greater genetic diversity within the family than previously recognized prompting us to suggest that additional genera should be considered within the family.

Introduction

Nodaviruses are positive-sense RNA viruses with bipartite genomes which are capped but not polyadenylated [1]. There are currently two genera recognized: Alphanodavirus (5 species) and Betanodavirus (4 species). The alphanodaviruses primarily infect insects, and betanodaviruses infect fish. Historically, the earliest isolations were of the alphanodaviruses in Australasia, initially with Nodamura virus (NoV) in Japan [2]; Flock House virus (FHV), and then Black beetle virus (BBV), and Boolarra virus (BoV) [3]. More recently, Pariacoto virus (PaV) [4] was isolated from a sample from Peru; the first alphanodavirus identified outside Australasia. The four recognized species in the betanodaviruses are Barfin flounder nervous necrosis virus (BFNNV), Redspotted grouper nervous necrosis virus (RGNNV), Striped jack nervous necrosis virus (SJNNV), and Tiger puffer nervous necrosis virus (TPNNV) [3]. In 2014, a nodavirus named mosinovirus (MoNV) was isolated from mosquito species of the Culex genus [5]. Currently, MoNV is not yet recognized as a species and has not been assigned to a genus. Most recently, a large number of nodavirus-like sequences were identified in a large-scale sequencing study of invertebrates [6]. This study has greatly enlarged our knowledge of the genetic diversity of this group of viruses, and will be complemented by future research on their biology.

The two genomic RNAs of nodaviruses, a ~3.1 kb RNA1 and ~1.4 kb RNA2 respectively, encode for two open reading frames (ORFs) [1]. RNA1 ORF1 (protein A) encodes the RNA-dependent RNA polymerase (RdRp) and the RNA2 ORF1 encodes the capsid protein. A third RNA identified in infected cells is sub-genomic and derived from RNA1, and potentially encodes two small proteins: B1 and B2. B2 has been linked to suppression of virus-inhibitory RNAi activity [7]. Interestingly, a fourth RNA of MoNV is also sub-genomic and encodes B2-like RNAi suppressor activity [5]. Also of interest is the observation that the MoNV RNA2 capsid seems to have originated from a distantly-related virus family as indicated by the genetic distance separating this genome segment from other nodaviruses [5].

In this work we describe a new species of nodavirus which was isolated from a sample of Culex annulirostris mosquitoes and which we have tentatively named culannivirus (CulV). Phylogenetic analysis of the RdRp sequence placed CulV distantly from any of the currently recognized genera. By contrast, phylogenetic analysis of the capsid protein placed CulV in a clade with the currently recognized alphanodaviruses. We suggest a re-organization of the family based on the RdRp sequence.

Materials and methods

Ethics statement

No specific permits were required for the described field studies, however permission was sought and obtained from the traditional land owners of Bradshaw Field Training Area to conduct mosquito trapping activities in the area. These field studies did not involve endangered or protected species.

Mosquito trapping and virus isolation

Adult mosquitoes were collected using CO2-baited light traps from the Bradshaw Field Training Area (BFTA) in the Northern Territory of Australia in 2014, using previously described methods [8]. Screening of mosquitoes for the presence of RNA viruses by mosquito homogenization, inoculation onto C6/36 cell monolayers and subsequent ELISA assessment was performed as previously described [9]. Successive passaging of the virus was performed by inoculating onto monolayers of C6/36 cells for 5–7 days before harvesting. Inoculations and virus passaging in C6/36 cells was performed in RPMI 1640 supplemented with 5% fetal bovine serum, 50 U/mL penicillin, 50 mg/mL streptomycin and 2 mM L-glutamine and incubation at 28 °C.

Sequencing and bioinformatics.

RNA extraction, first and second strand cDNA synthesis, library preparation and Illumina sequencing were as previously described [10]. A virus consensus sequence of RNA1 was assembled from the data using GeneiousPro v8.1 software and Hubei noda-like virus 5 (GenBank accession number KX883080.1) as a reference sequence. The assembly was then improved with iterations of the assembly using the consensus as a reference (GenBank accession number MH794142). RNA2 was obtained using a de novo assembly (GenBank accession number MH794143). To construct phylogenetic trees, global multiple protein sequence alignments of RNA1 and RNA2 ORFs were performed with the Geneious Alignment feature (Blosum62 cost matrix; gap open penalty, 12; gap extension penalty, 3) [11]. A mid-point rooted phylogenetic tree was then generated using the FastTree 2.1.5 (Jones-Taylor-Thornton model) optimizing for the gamma20 likelihood and using the Shimodaira-Hasegawa test to calculate the support values [12]. Within and between mean amino acid p-distances were calculated using MEGA7 [13] from the protein alignments using with complete deletion of gaps and a shape parameter of 5.

Results

An unidentified RNA virus was isolated during the screening of Culex annulirostris mosquitoes collected from the Northern Territory of Australia in 2014. We tentatively named this virus culannivirus (CulV). Preliminary analysis of deep sequencing of extracted RNA identified the virus as belonging to the Nodaviridae family. We performed more detailed analysis using RdRp protein sequence as it has been previously suggested that it would be preferable for taxonomic determinations [5]. Our alignments also seemed to support this, and phylogenetic analysis revealed at least 6 separate clades (Fig 1) all of which were well-supported (support values >90). The members of these clades were separated internally by genetic distances (p-distance) ranging from 0.36–0.55, whilst the groups themselves were separated by genetic distances of 0.65–0.80 (Table 1). Clade 1 corresponded to the alphanodaviruses with the notable omission of PaV. However, the grouping of PaV separate from the other alphanodaviruses has been previously observed when the RdRp was used for phylogenetic analysis [5]. This may indicate some re-assortment event in the history of this virus. Clade 2 included the betanodaviruses, but also contained additional viruses detected in a variety of arthropod hosts from recent sequencing studies [6, 14, 15]. The fish group, corresponding to the betanodavirues, composed a separate and much less diverse sub-clade (internal p-distance of 0.05). Clade 3 was a genetically diverse clade which contained PaV, MoNV and CulV in separate sub-clades. Notably, the latter two were both isolated from mosquitoes of the genus Culex. This clade has been greatly expanded by recent sequencing studies, with material again originating from diverse arthropod hosts [6]. Clade 4 consisted of several viruses sequenced in recent studies, but no isolate has been reported for this group. Hence, this group is one of the least biologically characterized but, interestingly, does contain two sequences identified in metagenomics studies of gut samples from a wolf and an otter [16, 17], perhaps originating from arthropods ingested by those animals. Whilst clades 3 and 4 were separate with strong support, together they also formed a larger clade (support value 99). Clade 5 includes one biologically characterized isolate obtained from the moth Helicoverpa zea as well as a number of nodavirus sequences from diverse insect hosts [18]. Clade 6 included sequences from virus isolates obtained from nematodes [19, 20]. These 6 clades indicate known taxonomic groupings, such as the recognized genera Alphanodavirus and Betanodavirus, or they may potentially reveal relationships based on biological characteristics which will be the basis of some future taxonomic delineation.

thumbnail
Fig 1. Phylogenetic analysis of the nodavirus RNA-dependent RNA polymerase sequence from RNA1.

A multiple sequence alignment was performed and a tree was generated. The 6 clades are numbered and highlighted. The alphanodaviruses correspond to clade 1 and the betanodaviruses are shown in a box with a dashed outline. Biologically characterized viruses are indicated in purple and sequences from recent metagenomics studies are given in black text. Branch support values are shown.

https://doi.org/10.1371/journal.pone.0210029.g001

The capsid protein sequence has been used as the major taxonomic determinant when assigning members of the nodaviruses [21]. Analysis of the capsid structure of several nodaviruses shows that members of this virus family have diverse capsid topologies [2227], suggesting that the capsid evolution may not be linear within the group. Available sequence was more limited for capsid. An alignment of capsid protein sequence was found to be of poor quality at the family level, making it difficult to determine the taxonomic relationships among more distantly-related nodaviruses. Phylogenetic trees constructed using this alignment revealed at least 4 clades corresponding to those identified in the RdRp tree (clades 1, 2, 3 and 6), but with weaker support values (Fig 2). Individual viruses grouped with clades that were identified for RdRp with exceptions such as CulV, which grouped with the alphanodaviruses. In addition, HzAM derived virus, Macrobrachium rosenbergii nodavirus and Penaeus vannamei nodavirus were also placed in different clades. Similar phylogenetic assignment differences with these viruses between RdRp and capsid were noted previously [5]. These data further suggest reassortment of genome RNAs during evolution.

thumbnail
Fig 2. Phylogenetic analysis of the nodavirus capsid sequence from RNA2.

A multiple sequence alignment was performed and a tree was generated. Viruses in the genus Alphanodavirus are highlighted in color. Viruses where there is some biological characterization are indicated in purple and sequences from recent metagenomics studies are given in black text. Branch support values are shown.

https://doi.org/10.1371/journal.pone.0210029.g002

Discussion

The outcome of nodavirus infection is influenced by both host and virus genetic determinants. Entry of the virus into the cell involves interaction between the virus capsid and cell membrane receptors. Studies on FHV have indicated that capsid protein gamma facilitates membrane penetration during host cell entry (reviewed in [22]). Once inside the cell, the viral capsid and RdRp proteins are directly translated from the genomic RNAs using the host translation machinery. The host protein Hsp90 has been indicated to be important for efficient translation of Protein A [23]. Nodavirus replication is driven by Protein A (RdRp) which binds to the outer mitochondrial membrane and induces spherules where replication of the viral proteins takes place [2426]. Following replication of the genomic RNAs and accumulation of the capsid proteins the virus self assembles and is released from the cell (reviewed in [27]). Many steps in the virus infection cycle involve interactions between the host and virus; hence, the genetic determinants of these interactions influence evolution of the virus thus impacting the phylogenetic resolution of the group.

The phylogenetic analysis of RdRp sequences has been reliably used for decades to determine the relationships of RNA viruses up to the family, and potentially higher, taxonomic level [2831]. Our work also supports the use of the RdRp as being a more reliable indicator of the taxonomic relationships among the nodaviruses than capsid, as has been recently suggested [5]. The RdRp analysis placed PaV distantly from clade 1, in agreement with previous analysis, suggesting that this virus may be more accurately assigned outside this genus. It was more closely related to MoNV and CulV, both of which were isolated from mosquitoes of the genus Culex; the latter virus being the subject of this work. The phylogenetic analysis revealed 6 monophyletic clades, two of which were constituted either wholly or partly of the existing genera. Clade 1 corresponded to the genus Alphanodavirus. Clade 2 included members of the betanodaviruses and additional viruses that had been isolated from organisms other than fish. The genetic diversity of this genus is relatively small (internal p distance of 0.05) and this possibly reflects either limited sampling or a genetic bottleneck in adapting to their common fish hosts. The internal diversity of clade 2 (internal p distance of 0.44) was larger than the alphanodaviruses (internal p distance of 0.36), but smaller than the internal diversity of three of the other clades identified (i.e. clades 3, 4 and 6). Clade 2 includes viruses collected from a variety of arthropod hosts, similarly to most of the other nodaviruses, as well as fish. Virus genetic groupings generally reflect their biological characteristics such as replication strategy, genome structure, host tropism and others. The nodaviruses have a common genome structure and share the ability to infect a variety of mostly arthropod hosts, which makes it difficult to determine taxonomic groupings below the family level. Future research may reveal biological characteristics which reflect the genetic groupings identified in this study.

In consideration of the greater genetic diversity revealed by this phylogenetic analysis, which includes metagenomics and other studies, there seems to be impetus for greater taxonomic division below the family level. On the basis of a characterized virus type strain, clear monophyletic groupings, which were well-separated from other clades and had good support values, were determined. Given the level of diversity identified, future biological characterization may reveal further genera within this family Nodaviridae.

Acknowledgments

The opinions expressed herein are those of the authors and do not necessarily reflect those of the Australian Defence Force and/or extant Defence Force Policy.

References

  1. 1. Friesen PD. Insect virses. In: Knipe DMaH, P M., editor. Fields Virology. 1. Philadelphia, USA: Lippincott Williams And Wilkins; 2007. p. 707–36.
  2. 2. Scherer WF, Hurlbut HS. Nodamura virus from Japan: a new and unusual arbovirus resistant to diethyl ether and chloroform. Am J Epidemiol. 1967;86(2):271–85. pmid:4383437.
  3. 3. Virus Taxonomy: The Classification and Nomenclature of Viruses. 2017 10.
  4. 4. Johnson KN, Zeddam JL, Ball LA. Characterization and construction of functional cDNA clones of Pariacoto virus, the first Alphanodavirus isolated outside Australasia. J Virol. 2000;74(11):5123–32. pmid:10799587.
  5. 5. Schuster S, Zirkel F, Kurth A, van Cleef KW, Drosten C, van Rij RP, et al. A unique nodavirus with novel features: mosinovirus expresses two subgenomic RNAs, a capsid gene of unknown origin, and a suppressor of the antiviral RNA interference pathway. J Virol. 2014;88(22):13447–59. pmid:25210176.
  6. 6. Shi M, Lin XD, Tian JH, Chen LJ, Chen X, Li CX, et al. Redefining the invertebrate RNA virosphere. Nature. 2016. pmid:27880757.
  7. 7. Aliyari R, Wu Q, Li HW, Wang XH, Li F, Green LD, et al. Mechanism of induction and suppression of antiviral immunity directed by virus-derived small RNAs in Drosophila. Cell Host Microbe. 2008;4(4):387–97. pmid:18854242.
  8. 8. Frances SP, Cooper RD, Rowcliffe KL, Chen N, Cheng Q. Occurrence of Ross River virus and Barmah Forest virus in mosquitoes at Shoalwater Bay military training area, Queensland, Australia. J Med Entomol. 2004;41(1):115–20. pmid:14989354.
  9. 9. McLean BJ, Hobson-Peters J, Webb CE, Watterson D, Prow NA, Nguyen HD, et al. A novel insect-specific flavivirus replicates only in Aedes-derived cells and persists at high prevalence in wild Aedes vigilax populations in Sydney, Australia. Virology. 2015;486:272–83. pmid:26519596.
  10. 10. Huang B, Pyke AT, McMahon J, Warrilow D. Complete Coding Sequence of a Case of Chikungunya Virus Imported into Australia. Genome Announc. 2017;5(19). pmid:28495775.
  11. 11. Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28(12):1647–9. Epub 2012/05/01. pmid:22543367.
  12. 12. Price MN, Dehal PS, Arkin AP. FastTree 2—approximately maximum-likelihood trees for large alignments. PLoS One. 2010;5(3):e9490. pmid:20224823.
  13. 13. Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol Biol Evol. 2016. pmid:27004904.
  14. 14. Li L, Victoria JG, Wang C, Jones M, Fellers GM, Kunz TH, et al. Bat guano virome: predominance of dietary viruses from insects and plants plus novel mammalian viruses. J Virol. 2010;84(14):6955–65. Epub 2010/05/14. pmid:20463061.
  15. 15. Webster CL, Waldron FM, Robertson S, Crowson D, Ferrari G, Quintana JF, et al. The Discovery, Distribution, and Evolution of Viruses Associated with Drosophila melanogaster. PLoS Biol. 2015;13(7):e1002210. pmid:26172158.
  16. 16. Conceicao-Neto N, Godinho R, Alvares F, Yinda CK, Deboutte W, Zeller M, et al. Viral gut metagenomics of sympatric wild and domestic canids, and monitoring of viruses: Insights from an endangered wolf population. Ecol Evol. 2017;7(12):4135–46. pmid:28649326.
  17. 17. Conceicao-Neto N, Zeller M, Heylen E, Lefrere H, Mesquita JR, Matthijnssens J. Fecal virome analysis of three carnivores reveals a novel nodavirus and multiple gemycircularviruses. Virol J. 2015;12:79. pmid:25986582.
  18. 18. Bai H, Wang Y, Li X, Mao H, Li Y, Han S, et al. Isolation and characterization of a novel alphanodavirus. Virol J. 2011;8:311. pmid:21682922.
  19. 19. Felix MA, Ashe A, Piffaretti J, Wu G, Nuez I, Belicard T, et al. Natural and experimental infection of Caenorhabditis nematodes by novel viruses related to nodaviruses. PLoS Biol. 2011;9(1):e1000586. pmid:21283608 accordance with the PLoS policy on competing interests she has been excluded from all stages of the review process for this article.
  20. 20. Franz CJ, Zhao G, Felix MA, Wang D. Complete genome sequence of Le Blanc virus, a third Caenorhabditis nematode-infecting virus. J Virol. 2012;86(21):11940. pmid:23043172.
  21. 21. Sahul Hameed AS, Ninawe AS, Nakai T, Chi SC, Johnson KL, Ictv Report C. ICTV Virus Taxonomy Profile: Nodaviridae. J Gen Virol. 2018. pmid:30431412.
  22. 22. Odegard A, Banerjee M, Johnson JE. Flock house virus: a model system for understanding non-enveloped virus entry and membrane penetration. Curr Top Microbiol Immunol. 2010;343:1–22. pmid:20407886.
  23. 23. Castorena KM, Weeks SA, Stapleford KA, Cadwallader AM, Miller DJ. A functional heat shock protein 90 chaperone is essential for efficient flock house virus RNA polymerase synthesis in Drosophila cells. J Virol. 2007;81(16):8412–20. pmid:17522196.
  24. 24. Miller DJ, Ahlquist P. Flock house virus RNA polymerase is a transmembrane protein with amino-terminal sequences sufficient for mitochondrial localization and membrane insertion. J Virol. 2002;76(19):9856–67. pmid:12208963.
  25. 25. Miller DJ, Schwartz MD, Ahlquist P. Flock house virus RNA replicates on outer mitochondrial membranes in Drosophila cells. J Virol. 2001;75(23):11664–76. pmid:11689648.
  26. 26. Van Wynsberghe PM, Chen HR, Ahlquist P. Nodavirus RNA replication protein a induces membrane association of genomic RNA. J Virol. 2007;81(9):4633–44. pmid:17301137.
  27. 27. Venter PA, Schneemann A. Recent insights into the biology and biomedical applications of Flock House virus. Cell Mol Life Sci. 2008;65(17):2675–87. pmid:18516498.
  28. 28. Cerny J, Cerna Bolfikova B, Valdes JJ, Grubhoffer L, Ruzek D. Evolution of tertiary structure of viral RNA dependent polymerases. PLoS One. 2014;9(5):e96070. pmid:24816789.
  29. 29. Koonin EV. The phylogeny of RNA-dependent RNA polymerases of positive-strand RNA viruses. J Gen Virol. 1991;72 (Pt 9):2197–206. pmid:1895057.
  30. 30. Simmonds P. Methods for virus classification and the challenge of incorporating metagenomic sequence data. J Gen Virol. 2015;96(Pt 6):1193–206. pmid:26068186.
  31. 31. Ward CW. Progress towards a higher taxonomy of viruses. Res Virol. 1993;144(6):419–53. pmid:8140287.