Research | Open | Published:
Fecal virome analysis of three carnivores reveals a novel nodavirus and multiple gemycircularviruses
Virology Journalvolume 12, Article number: 79 (2015)
More knowledge about viral populations in wild animals is needed in order to better understand and assess the risk of zoonotic diseases. In this study we performed viral metagenomic analysis of fecal samples from three healthy carnivores: a badger (Meles meles), a mongoose (Herpestes ichneumon) and an otter (Lutra lutra) from Portugal.
We detected the presence of novel highly divergent viruses in the fecal material of the carnivores analyzed, such as five gemycircularviruses. Four of these gemycircularviruses were found in the mongoose and one in the badger. In addition we also identified an RNA-dependent RNA polymerase gene from a putative novel member of the Nodaviridae family in the fecal material of the otter.
Together these results underline that many novel viruses are yet to be discovered and that fecal associated viruses are not always related to disease. Our study expands the knowledge of viral species present in the gut, although the interpretation of the true host species of such novel viruses needs to be reviewed with great caution.
With the advent of next generation sequencing techniques, samples from a wide range of animal species have been screened to identify novel viruses and this approach has become the most important tool for early detection and characterization of possible emerging zoonotic agents [1–3]. It is important to monitor these emerging zoonotic agents as they can be responsible for minor or major epidemics worldwide . For example, zoonoses can range from the Middle East respiratory syndrome coronavirus (MERS-CoV), which recently drew a lot of attention worldwide , to the 2009 pandemic H1N1 influenza A virus or to the less publicized Hepatitis E virus . Not only are humans at risk, but also animals can be infected with viruses from other host species, resulting in disease, or further transmission to humans. The severe acute respiratory syndrome (SARS) coronavirus pandemic originated from wildlife, where bats where identified as the reservoir and civets as an intermediate host [7, 8]. A similar example comes from Hendra virus, which also originated from bats, but transmission to humans occurred with horses as an intermediate host, causing severe pneumonia in horses and humans . As such, zoonotic infections can have significant consequences for animal and public welfare. However, little is known about these pathogens before they emerge from unrecognized zoonotic sources and therefore a deeper understanding of the virome of wild animals will allow us to more rapidly identify the host of particular novel zoonotic viruses, and act appropriately to prevent further spread of such viruses.
Fourteen species of wild carnivores can be found in Portugal, often in relative close contact with humans. Only red fox, stone marten, badger, common genet and, more recently, the Egyptian mongoose, are known to have a generalized distribution throughout the country . From these animals we sampled two species, a badger (Meles meles) and an Egyptian mongoose (Herpestes ichneumon). In addition, we sampled a Eurasian otter (Lutra lutra) raised in a zoo, since otters are also widely distributed in Portugal . Currently little is known about the viral communities that populate the gut of these animals. Thus far, Bodewes and colleagues investigated the fecal virome of badgers and an otter from Spain, and found a fecal phlebovirus in an otter . A study by van den Brand and colleagues investigated the virome of badgers in the Netherlands, identifying two novel circularviruses . In Portugal, a study performed by Oliveira and colleagues screened otters for parvoviruses, adenoviruses and parainfluenza virus, however none of the samples had detectable levels of virus . Also in Portugal, another study screened for and identified parvoviruses in genets, badgers and mongooses . As these findings are probably only the tip of the iceberg, we were interested to further explore the viral communities of the gut in widely spread wild species in Portugal and a zoo specimen.
Results and discussion
Identification of five novel gemycircularviruses in the fecal samples from the badger and mongoose
In recent years, ssDNA viruses have been frequently found in fecal samples of a variety of animal species, such as badgers , bats , cows , turkey , rodents [17, 18], chimpanzees , pigs [20–22], New Zealand fur seal  foxes  and recently also in ancient caribou feces . ssDNA viruses are small viruses (1.4-8.5 kb) and can encode as little as two proteins: a capsid and a replication associated protein. Thus far, seven major ssDNA families have been identified based on the host range and type of ssDNA genome. Only the Parvoviridae family encodes for a linear genome, the remaining are circular. The Nanoviridae and Geminiviridae families infect plants, whereas Circoviridae, Parvoviridae and Anelloviridae are known to infect animals. Viruses belonging to the Inoviridae and Microviridae infect bacteria .
In our study, four complete novel gemycircularviruses were found in the mongoose and one in the badger feces. The first gemycircularvirus was discovered in 2010 in fungi . Since then, these novel ssDNA circular viruses have been found in a wide range of hosts, including the cassava plant , badger feces , mosquitoes , dragonflies , Hypericum japonicum, a flowering plant in the family Hypericaceae , as well as in the gut of a wide variety of mammals and birds from New Zealand . Although these viruses were initially named gemini-like viruses, the novel genus gemycircularvirus was recently proposed by Rosario and colleagues . Gemycircularviruses encode for a highly variable capsid protein and a more conserved replication associated protein (Rep). Figure 1 shows the phylogenetic analysis of known gemycircularviruses based on the amino acid level of the Rep gene. The viruses identified in this study showed to be distantly related to the currently known gemycircularviruses. The first gemycircularvirus found in the mongoose, named ‘Mongoose feces-associated gemycircularvirus a’, was found to be most closely related to the Badger feces-associated gemycircularvirus, sharing 57.2 % similarity on the amino-acid level. On the amino-acid level, ‘Mongoose feces-associated gemycircularvirus b’ shares its highest similarity (61.0 %) with ‘Mongoose feces-associated gemycircularvirus c’. The most closely related virus to the ‘Mongoose feces-associated gemycircularvirus d’ shares 71.7 % amino-acid similarity, and was found in the New Zealand bird Chatham gerygone (Gemycircularvirus 6 isolate P24a). All five gemycircularviruses contained nonanucleotide stem loop motifs, conserved Rolling Circle Replication (RCR) motifs (I, II, III and geminivirus-like Rep sequence motif (GRS)) and helicase motifs (Walker-A and Walker-B) (Table 1). Only the Mongoose feces-associated gemycircularvirus d had a nonanucleotide motif identical to the one found in most known gemycircularviruses (TAATATTAT), but the Rep motifs of the newly discovered viruses also showed a variable degree of similarity among each other (Table 1) and with previously characterized gemycircularviruses [28, 20]. The discovery of these novel viruses expands the knowledge on gemycircularviruses genetic diversity and their putative host range. Initially, it was hypothesized that these viruses infect fungi, since one of these viruses is known to infect the pathogenic fungi Sclerotinia sclerotiorum , however it is not known whether the remaining viruses do. Most recently, clones from a novel gemycircularvirus found in caribou feces were inoculated in a plant, resulting in a successful infection . Nevertheless, the true host(s) of these viruses remain to be determined. In our case, these viruses might have infected the mongoose and the badger, or alternatively, might have arisen from fungi inhabiting their intestines, or they could also be derived from insects or plants as part of the diet of the badger and the mongoose. Therefore, since the true host cannot be determined yet, the nomenclature of these novel viruses should be addressed cautiously. We tentatively named them feces-associated gemycircularviruses, preceded by the common name of the animal where they were isolated from.
Identification a partial RNA-dependent RNA polymerase of a nodavirus in the fecal material of the otter
The Nodaviridae family comprises two genera Alphanodavirus and the Betanodavirus of bipartite single stranded RNA viruses . Alphanodaviruses are usually insect viruses whereas betanodaviruses infect fish and are responsible for viral nervous necrosis in numerous fish species [32, 33]. Nodaviruses have two segments, RNA1 (3.2 kb) encodes for a RNA-dependent RNA polymerase (RdRp) responsible for its RNA replication and RNA2 (1.2 kb) encodes for a capsid protein . Nodaviruses are classified by the ICTV (International Committee on Taxonomy of Viruses) according to the genetic diversity of the RNA2 segment . In the fecal material of the otter we identified a partial RNA1 of a novel putative Nodavirus (1.7 kb), which was most closely related to the recently discovered Mosinovirus, isolated from mosquitoes, sharing 43 % similarity on the amino-acid level (Fig. 2). Adopting the convention of naming based on Schuster and colleagues for Mosinovirus virus (mosquito nodavirus), we have tentatively named the virus Lunovirus (Lutra lutra nodavirus) . As the RNA2 of Lunovirus was not found, most likely due to the fact that it is highly divergent and could not be detected by similarity searches in current viral databases, we should be reluctant to suggest a final classification. However, based on the large divergence of the Lunovirus RNA1 with the RNA1 of other nodaviruses, it seems likely that the Lunovirus is a novel member of the Nodaviridae family.
In the three healthy carnivores analyzed, viral sequences belonging to the Caudovirales order of bacteriophages were also detected, as previously reported . Furthermore, our study showed that even healthy wildlife seems to harbor many divergent viral communities that deserve to be explored further to expand our current knowledge and databases. From the limited data available from the fecal virome studies from an otter (Lutra lutra) and badger (Meles meles) of Bodewes and van der Brand [12, 13], the virus families discovered in our study are completely different, possibly resulting from different diets, or indicating a large unexplored area of the ‘virome space’. Also the fact that the viruses discovered are highly different from viruses available in databases, might explain the difficulty in finding viruses by regular PCR screening, as previously reported . These novel viruses reported in this study are likely derive from the diet, as the Nodavirus from the otter is likely to be from fish and the gemycircularviruses from insects, which are part of the animals diets. Viral discovery can be challenging because novel viruses, as seen in this study, can be quite divergent and their classification and true host determination can be difficult. In this regard, replication associated proteins have shown to be conserved, and the best strategy to create alignments. It is very interesting to see the ubiquity of different circular virus species found nowadays due to the availability of next-generation sequencing. Screening of larger groups of animals and species will help to increase our knowledge of viruses circulating in wild animals.
Fecal samples were collected from a badger (Meles meles), a mongoose (Herpestes ichneumon) from a rescue center and an Eurasian otter (Lutra lutra) from a zoo upon their arrival in 2011 at the wildlife center “Parque Biológico da Serra da Lousã” in Coimbra district, Portugal. This collection was part of the quarantine assessment program applied to newly introduced animals in the center. Samples were kept at −80 °C until further processing.
Ten percent fecal suspensions were homogenized for 1 min at 3000 rpm with a MINILYS homogenizer (Bertin Technologies) and filtered consecutively through 100 μm, 10 μm and 0.8 μm membrane filters (Millipore) for 30 s at 1250 g. The filtrate was then treated with a homemade buffer (1 M Tris, 100 mM CaCl2 and 30 mM MgCl2) and a cocktail of Benzonase (Novagen) and Micrococcal Nuclease (New England Biolabs) at 37 °C for 2 h to digest free-floating nucleic acids. RNA and DNA was extracted using the QIAamp Viral RNA Mini Kit (Qiagen) according to the manufacturer’s instructions but without addition of carrier RNA to the lysis buffer. First and second strand synthesis and random PCR amplification for 25 cycles were performed using a slightly modified Whole Transcriptome Amplification (WTA) Kit procedure (Sigma-Aldrich), with a denaturation temperature of 95 °C instead of 72 °C to allow for the denaturation of dsDNA and dsRNA. This modification leads to the amplification of both RNA and DNA. A size selection after library synthesis was performed using a 0.7 ratio of Agencourt AMPure XP beads (Beckman Coulter, Inc.). WTA products were purified with MSB Spin PCRapace spin columns (Stratec) and were prepared for Illumina sequencing using the KAPA Library Preparation Kit (Kapa Biosystems). Fragments ranging from 350–600 bp were selected using the BluePippin (Sage Science) according to the manufacturer’s instructions. Libraries were quantified with the KAPA Library Quantification kit (Kapa Biosystems) and sequencing of the samples was performed on a HiSeq 2500 platform (Illumina) for 301 cycles (150 bp paired ends). Each sample was attributed a total of 2 million paired end reads. Mongoose feces associated gemycircularvirus a, b, c and d yielded 262, 295, 37,644, and 356 reads respectively. The Badger feces associated gemycircularvirus yielded 88 reads and the RNA2 of the Nodavirus 20 reads.
Genomic and phylogenetic analysis
Raw reads were trimmed for quality and adapters using Skewer  and were de novo assembled into contigs using SPAdes . Scaffolds were classified using a tBLASTx search against all complete viral genomes in GenBank using an e-value cut-off of 10−10. Scaffolds with a significant tBLASTx hit were retained and used for a second tBlastx search against the GenBank nucleotide database using an e-value of 10−4 . Open reading frames (ORF) were identified with ORF Finder analysis tools and the conserved motifs in the amino acid sequences were identified with HMMER . Amino acid alignments of the viral sequences were performed with MAFFT version 7  using the --auto option. Maximum likelihood phylogenetic trees were constructed in MEGA6.0 , using JCC (best substitution model) with 500 bootstrap replicates. Potential intron acceptor and donor sites from the novel gemycircularviruses were identified manually. Using the method above we were able to retrieve the five complete gemycircularvirus and partial sequence of the RNA2 of the Lunovirus. Presence of the discovered novel viruses was then confirmed by PCR and Sanger sequencing using the original extracts. Gemycircularviruses’ primers were designed in the capsid gene, covering the complete circular genome and the nodavirus primer pair targeted the RNA2 sequence found. All sequences from the novel viruses were submitted to GenBank [KP263543, KP263544, KP263543-KP263548, KP263546, KP263547, KP263548].
Morse SS, Mazet JA, Woolhouse M, Parrish CR, Carroll D, Karesh WB, et al. Prediction and prevention of the next pandemic zoonosis. Lancet. 2012;380(9857):1956–65. doi:10.1016/S0140-6736(12)61684-5.
Wang LF. Discovering novel zoonotic viruses. N S W Public Health Bull. 2011;22(5–6):113–7. doi:10.1071/NB10078.
Temmam S, Davoust B, Berenger JM, Raoult D, Desnues C. Viral metagenomics on animals as a tool for the detection of zoonoses prior to human infection? Int J Mol Sci. 2014;15(6):10377–97. doi:10.3390/ijms150610377.
Cutler SJ, Fooks AR, van der Poel WH. Public health threat of new, reemerging, and neglected zoonoses in the industrialized world. Emerg Infect Dis. 2010;16(1):1–7. doi:10.3201/eid1601.081467.
Al-Tawfiq JA, Zumla A, Memish ZA. Travel implications of emerging coronaviruses: SARS and MERS-CoV. Trav Med Infec Dis. 2014. doi:10.1016/j.tmaid.2014.06.007.
Yugo DM, Meng XJ. Hepatitis E virus: foodborne, waterborne and zoonotic transmission. Int J Environ Res Public Health. 2013;10(10):4507–33. doi:10.3390/ijerph10104507.
Wang LF, Eaton BT. Bats, civets and the emergence of SARS. Curr Top Microbiol Immunol. 2007;315:325–44.
Li W, Shi Z, Yu M, Ren W, Smith C, Epstein JH, et al. Bats are natural reservoirs of SARS-like coronaviruses. Science. 2005;310(5748):676–9. doi:10.1126/science.1118391.
Middleton D. Hendra Virus. Vet Clin North Am Equine Pract. 2014;30(3):579–89. doi:10.1016/j.cveq.2014.08.004.
Duarte MD, Henriques AM, Barros SC, Fagulha T, Mendonca P, Carvalho P, et al. Snapshot of viral infections in wild carnivores reveals ubiquity of parvovirus and susceptibility of Egyptian mongoose to feline panleukopenia virus. PLoS One. 2013;8(3):e59399. doi:10.1371/journal.pone.0059399.
Oliveira M, Sales-Luís T, Duarte A, Nunes SF, Carneiro C, Tenreiro T, et al. First assessment of microbial diversity in faecal microflora of Eurasian otter (Lutra lutra Linnaeus, 1758) in Portugal. Eur J Wildl Res. 2008;54:245–52. doi:10.1007/s10344-007-0137-8.
Bodewes R, Ruiz-Gonzalez A, Schapendonk CM, van den Brand JM, Osterhaus AD, Smits SL. Viral metagenomic analysis of feces of wild small carnivores. Virol J. 2014;11:89. doi:10.1186/1743-422X-11-89.
van den Brand JM, van Leeuwen M, Schapendonk CM, Simon JH, Haagmans BL, Osterhaus AD, et al. Metagenomic analysis of the viral flora of pine marten and European badger feces. J Virol. 2012;86(4):2360–5. doi:10.1128/JVI.06373-11.
Ge X, Li J, Peng C, Wu L, Yang X, Wu Y, et al. Genetic diversity of novel circular ssDNA viruses in bats in China. J Gen Virol. 2011;92(Pt 11):2646–53. doi:10.1099/vir.0.034108-0.
Kim HK, Park SJ, Nguyen VG, Song DS, Moon HJ, Kang BK, et al. Identification of a novel single-stranded, circular DNA virus from bovine stool. J Gen Virol. 2012;93(Pt 3):635–9. doi:10.1099/vir.0.037838-0.
Reuter G, Boros A, Delwart E, Pankovics P. Novel circular single-stranded DNA virus from turkey faeces. Arch Virol. 2014. doi:10.1007/s00705-014-2025-3.
Phan TG, Kapusinszky B, Wang C, Rose RK, Lipton HL, Delwart EL. The fecal viral flora of wild rodents. PLoS Pathog. 2011;7(9):e1002218. doi:10.1371/journal.ppat.1002218.
Nishiyama S, Dutia BM, Stewart JP, Meredith AL, Shaw DJ, Simmonds P, et al. Identification of novel anelloviruses with broad diversity in UK rodents. J Gen Virol. 2014;95(Pt 7):1544–53. doi:10.1099/vir.0.065219-0.
Blinkova O, Victoria J, Li Y, Keele BF, Sanz C, Ndjango JB, et al. Novel circular DNA viruses in stool samples of wild-living chimpanzees. J Gen Virol. 2010;91(Pt 1):74–86. doi:10.1099/vir.0.015446-0.
Sikorski A, Dayaram A, Varsani A. Identification of a novel circular DNA virus in New Zealand fur seal (Arctocephalus forsteri) fecal matter. Genome Annc. 2013;1(4). doi:10.1128/genomeA.00558-13.
Cheung AK, Ng TF, Lager KM, Bayles DO, Alt DP, Delwart EL, et al. A divergent clade of circular single-stranded DNA viruses from pig feces. Arch Virol. 2013;158(10):2157–62. doi:10.1007/s00705-013-1701-z.
Lager KM, Ng TF, Bayles DO, Alt DP, Delwart EL, Cheung AK. Diversity of viruses detected by deep sequencing in pigs from a common background. J Vet Diagn Invest. 2012;24(6):1177–9. doi:10.1177/1040638712463212.
Bodewes R, van der Giessen J, Haagmans BL, Osterhaus AD, Smits SL. Identification of multiple novel viruses, including a parvovirus and a hepevirus, in feces of red foxes. J Virol. 2013;87(13):7758–64. doi:10.1128/JVI.00568-13.
Ng TF, Chen LF, Zhou Y, Shapiro B, Stiller M, Heintzman PD et al. Preservation of viral genomes in 700-y-old caribou feces from a subarctic ice patch. Proc Natl Acad Sci U S A. 2014. doi:10.1073/pnas.1410429111.
King AM, Adams MJ, Carstens EB, Lefkowitz EJ. Virus taxonomy: ninth report of the international committee on taxonomy of viruses. San Diego, CA: Academic; 2011.
Yu X, Li B, Fu Y, Jiang D, Ghabrial SA, Li G, et al. A geminivirus-related DNA mycovirus that confers hypovirulence to a plant pathogenic fungus. Proc Natl Acad Sci U S A. 2010;107(18):8387–92. doi:10.1073/pnas.0913535107.
Dayaram A, Opong A, Jaschke A, Hadfield J, Baschiera M, Dobson RC, et al. Molecular characterisation of a novel cassava associated circular ssDNA virus. Virus Res. 2012;166(1–2):130–5. doi:10.1016/j.virusres.2012.03.009.
Rosario K, Dayaram A, Marinov M, Ware J, Kraberger S, Stainton D, et al. Diverse circular ssDNA viruses discovered in dragonflies (Odonata: Epiprocta). J Gen Virol. 2012;93(Pt 12):2668–81. doi:10.1099/vir.0.045948-0.
Du Z, Tang Y, Zhang S, She X, Lan G, Varsani A, et al. Identification and molecular characterization of a single-stranded circular DNA virus with similarities to Sclerotinia sclerotiorum hypovirulence-associated DNA virus 1. Arch Virol. 2014;159(6):1527–31. doi:10.1007/s00705-013-1890-5.
Sikorski A, Massaro M, Kraberger S, Young LM, Smalley D, Martin DP, et al. Novel myco-like DNA viruses discovered in the faecal matter of various animals. Virus Res. 2013;177(2):209–16. doi:10.1016/j.virusres.2013.08.008.
Ball LA, Wohlrab B, Li Y. Nodavirus RNA replication: mechanism and harnessing to vaccinia virus recombinants. Arch Virol Suppl. 1994;9:407–16.
Shetty M, Maiti B, Shivakumar Santhosh K, Venugopal MN, Karunasagar I. Betanodavirus of marine and freshwater fish: distribution, genomic organization, diagnosis and control measures. Indian J Virol. 2012;23(2):114–23. doi:10.1007/s13337-012-0088-x.
Bailey L, Scott HA. The pathogenicity of Nodamura virus for insects. Nature. 1973;241(5391):545.
Thiéry R, Johnson KL, Nakai T, Schneemann A, Bonami JR, Lightner DV. Family Nodaviridae. In: King AMQ, Lefkowitz E, Adams MJ, Carstens EB, editors. Virus Taxonomy: Ninth Report of the International Committee on Taxonomy of Viruses. London: Elsevier Academic Press; 2011. p. 1061–7.
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. doi:10.1128/JVI.02144-14.
Jiang H, Lei R, Ding SW, Zhu S. Skewer: a fast and accurate adapter trimmer for next-generation sequencing paired-end reads. BMC Bioinf. 2014;15:182. doi:10.1186/1471-2105-15-182.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–77. doi:10.1089/cmb.2012.0021.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10. doi:10.1016/S0022-2836(05)80360-2.
Finn RD, Clements J, Eddy SR. HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 2011;39(Web Server issue):W29–37. doi:10.1093/nar/gkr367.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80. doi:10.1093/molbev/mst010.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9. doi:10.1093/molbev/mst197.
NCN and MZ were supported by the Institute for the Promotion of Innovation through Science and Technology in Flanders (IWT Vlaanderen). We would like to thank Prof. Maria São José Nascimento and Joren Raymenants who assisted in the proof-reading of the manuscript.
The authors declare that they have no competing interests.
NCN conceived, designed and performed the experiments, analyzed the data, drafted the manuscript. MZ and EH performed data analysis and contributed to manuscript drafting. HL was involved in the experiment performance and data analysis. JRM collected the samples and drafted the manuscript. JM conceived and designed the study and contributed to data analysis and manuscript drafting. All authors read and approved the final manuscript.