Friday, November 22, 2019

Distantly Related Rotaviruses in Common Shrews, Germany, 2004–2014 - Volume 25, Number 12—December 2019 - Emerging Infectious Diseases journal - CDC

Distantly Related Rotaviruses in Common Shrews, Germany, 2004–2014 - Volume 25, Number 12—December 2019 - Emerging Infectious Diseases journal - CDC

Issue Cover for Volume 25, Number 12—December 2019

Volume 25, Number 12—December 2019
Dispatch

Distantly Related Rotaviruses in Common Shrews, Germany, 2004–2014

Reimar JohneComments to Author , Simon H. Tausch, Josephine Grützke, Alexander Falkenhagen, Corinna Patzina-Mehling, Martin Beer, Dirk Höper, and Rainer G. Ulrich
Author affiliations: German Federal Institute for Risk Assessment, Berlin, Germany (R. Johne, S.H. Tausch, J. Grützke, A. Falkenhagen, C. Patzina-Mehling)Friedrich-Loeffler-Institut, Greifswald-Insel Riems, Germany (M. Beer, D. Höper, R.G. Ulrich)Deutsches Zentrum für Infektionsforschung, partner site Hamburg-Lübeck-Borstel-Insel Riems, Germany (R.G. Ulrich)

Abstract

We screened samples from common shrews (Sorex araneus) collected in Germany during 2004–2014 and identified 3 genetically divergent rotaviruses. Virus protein 6 sequence similarities to prototype rotaviruses were low (64.5% rotavirus A, 50.1% rotavirus C [tentative species K], 48.2% rotavirus H [tentative species L]). Shrew-associated rotaviruses might have zoonotic potential.
Rotaviruses are a major cause of diarrhea in young children, causing an estimated 215,000 deaths worldwide every year (1). These viruses are nonenveloped and have a genome consisting of 11 segments of double-stranded RNA (2); each segment codes for either 1 of the structural proteins, virus protein (VP) 1–7, or 1 or 2 of the nonstructural proteins (NSPs), NSP1–6. Rotaviruses are classified into species A–I or the tentative species J on the basis of the amino acid sequence similarity of the conserved structural protein VP6 and the conserved nucleotide sequence of the genome segment ends (35). For rotavirus A, further classification into genome segment–specific genotypes has been established (6). Rotaviruses can infect a wide diversity of animals, and zoonotic transmission of rotaviruses has been reported (7).
Shrews are small insectivorous mammals that have been previously identified as reservoirs for other pathogens (e.g., hantaviruses and Leptospira spp.) (810). In this investigation, we aimed to determine whether common shrews (Sorex araneus, order Eulipotyphla) are also a reservoir for rotaviruses and, if so, assess the genetic variability of the viruses found in this species.

The Study


Thumbnail of Distribution of common shrews (Sorex araneus) collected at monitoring sites (9) and additional sites (10) in Germany, 2004–2014, positive and negative for RVA, RVC-like, and RVH-like species by reverse transcription PCR. Numbers in white circles indicate the number of negative samples at that collection site; white circles without numbers indicate 1 negative sample at that site. Circles with multiple colors indicate animals with co-infections. The collection sites of the 2 samples a
Figure 1. Distribution of common shrews (Sorex araneus) collected at monitoring sites (9) and additional sites (10) in Germany, 2004–2014, positive and negative for RVA, RVC-like, and RVH-like species...
During 2004–2014, small mammals were caught in different regions of Germany as part of local monitoring or pest control measures (9,10). From these animal collections, we acquired samples (intestine contents) collected from 49 common shrews (Figure 1). We combined these samples almost equally into 2 pools and performed RNA extraction followed by next-generation sequencing (NGS) using the Ion Torrent Personal Genome Machine system (ThermoFisher Scientific, https://www.thermofisher.comExternal LinkAppendix). By applying the RIEMS data analysis pipeline (11), we identified 3 short contigs with low sequence similarities to rotavirus H in both pools. To identify the positive animals, we extracted RNA from individual samples and screened for rotavirus RNA using reverse transcription PCR (RT-PCR) with primers specific to 1 of the 3 rotavirus H contigs we previously obtained (Appendix Table 1). In total, 7 (15.2%) of 46 samples turned out to be positive for species H–like rotavirus (Table 1); 2 samples, KS/12/0644 and KS/11/2281, generated the strongest signal on ethidium bromide staining. We subjected these 2 samples to RNase and DNase treatment followed by RNA extraction and NGS using the NextSeq 500 sequencing system (Illumina, https://www.illumina.comExternal Link); 8,576,782 read pairs for KS/12/0644 and 6,168,437 for KS/11/2281 were generated. After a RAMBO-K analysis (12) suggested a low abundance of highly deviant rotavirus sequences, we performed data analysis and contig assembly using a newly generated pipeline (Appendix). By this method, contig lengths were 164–3,017 nt, and we obtained 48 contigs with sequence similarities to rotavirus A, 17 with low sequence similarities to rotavirus C, and 23 with low sequence similarities to rotavirus H (Appendix Table 2). Because contigs of homologous genes from each of the 3 viruses were detected in these samples, we concluded 3 different rotaviruses were present in both.

We then performed RT-PCR with all samples using primers specific to the species A and species C–like rotavirus contigs from the previous analysis, and 21.7% (10/46) were positive for species A rotavirus and 10.9% (5/46) for species C–like rotavirus; rotavirus co-infections were also identified (Table 1). An analysis of the geographic distribution of shrew rotaviruses in Germany shows that species C–like rotaviruses were mainly located in the northeast and southwest, species H–like rotaviruses mainly in the south, and species A rotaviruses broadly throughout (Figure 1). At the monitoring site in Baden-Württemberg (southwest Germany), frequent detections of different rotaviruses and multiple co-infections were observed.
Despite several efforts, we could delineate only partial genomic sequences of rotaviruses from the NGS data. By application of primer ligation, rapid amplification of cDNA ends, and degenerated primer RT-PCR strategies, we acquired the complete open reading frames of VP1, VP6, and NSP5 of most viruses (Table 2). In addition, we reamplified and sequenced the VP6 genes of all viruses by dideoxy chain–termination sequencing and confirmed the VP6 sequences obtained. Sequence analysis of these genes and in silico translation indicated 14.1%–65.6% amino acid sequence similarity to the respective proteins of other rotavirus species (Table 2). The Rotavirus Classification Working Group reviewed the sequences of the shrew rotavirus A genes in sample KS/11/2281 and designated the new genotypes R23 for VP1, I27 for VP6, and H23 for NSP5. The maximum amino acid sequence similarities to established rotavirus type species of 50.1% for VP6 of species C–like rotavirus and 48.2% (species H) or 48.3% (species J) for VP6 of species H–like rotavirus suggest that these viruses should be classified as novel (tentative) rotavirus species (Table 2).

Thumbnail of Phylogenetic relationship of shrew rotaviruses (bold), Germany, 2011–2012, with RVA–RVJ determined by using the deduced amino acid sequences of virus protein 1 (A), virus protein 6 (B), and nonstructural protein 5 (C). Trees were constructed by using a neighbor-joining method implemented in the MegAlign module of DNASTAR (https://www.dnastar.com) and a bootstrap analysis with 1,000 trials and 111 random seeds. Bootstrap values of >50% are shown. The rotavirus species, host, strai
Figure 2. Phylogenetic relationship of shrew rotaviruses (bold), Germany, 2011–2012, with RVA–RVJ determined by using the deduced amino acid sequences of virus protein 1 (A), virus protein 6 (B), and nonstructural protein 5...
Phylogenetic analyses of the VP1, VP6, and NSP5 proteins indicate a consistent branching of shrew rotavirus A with other rotavirus A species and shrew species C–like rotavirus with other rotavirus C species. However, the species H–like rotavirus branches more variably within the rotavirus B-G-H-I-J cluster (Figure 2). A more detailed phylogenetic analysis of complete and additional partial genome segment nucleotide sequences of the shrew rotavirus A showed a basal branching at the cluster of other species A rotavirus sequences for most genes (Appendix Figure 1). In addition, phylogenetic analyses of partial amino acid sequences deduced from other genes of the shrew species C–like and H–like rotaviruses confirmed the relationship evident from analyses of the 3 completely sequenced open reading frames (Appendix Figure 2–4).

Shrews have been analyzed infrequently for rotaviruses. In 1 study, rotavirus antigen was detected in wild Chinese tree shrews (Tupaia chinensis, order Scandentia) (13), and in another study, species A rotaviruses not identical to those of our study (Appendix Figure 1) were identified in house shrews (Suncus murinus, order Eulipotyphla) from China (14). Here, a broader rotavirus screening of common shrew samples resulted in the identification of novel rotaviruses. The rotavirus detection rate of 10.9%–21.7% in the analyzed samples from animals from different regions of Germany suggests a wide circulation of rotaviruses in shrews, although more samples should be analyzed in the future to clarify the association of rotaviruses with these animals. We also identified co-infections with >1 rotavirus, a regular finding in other animal host species (15).
The shrew rotavirus A sequences showed low similarities with other species A rotaviruses, resulting in the assignment of novel genotypes and suggesting a long-term separate evolution of these viruses in this shrew species. The 2 other rotaviruses identified showed even lower sequence similarities to the known rotavirus species. According to the cutoff value of 53% suggested for the differentiation of rotavirus species on the basis of the encoded VP6 amino acid sequence (5), both viruses should be considered new rotavirus species, which we tentatively designate rotavirus species K (for the rotavirus C–like species) and L (for the rotavirus H–like species). However, because their complete genome sequences have not been determined, a final classification of these viruses remains to be accomplished. At least the 5′ and 3′ termini of these rotavirus genome segments, which are conserved within known rotavirus species (2), should be determined. The low virus amounts in samples, restricted available sample volumes, presence of multiple viruses in single samples, and low sequence similarities for some virus genes might help explain the failure to generate complete genome sequences in our study.

Conclusions

We identified multiple, genetically divergent rotavirus species in common shrews in Germany. These animals should be further investigated as a potential reservoir for rotaviruses capable of infecting humans.
Dr. Johne is a senior scientist at the German Federal Institute for Risk Assessment, Berlin, Germany. His research interests are the epidemiology, molecular evolution, pathogenesis, and transmission pathways of foodborne viruses.

Acknowledgments


The excellent technical assistance of Anke Mandelkow and Patrick Zitzow and the generation of Figure 1 by Patrick Wysocki are kindly acknowledged.
Sample collection and distribution were supported in part by the German Federal Ministry of Education and Research through the German Research Platform for Zoonoses (network rodent-borne pathogens, FKZ 01KI1018 and 01KI1303 to R.G.U.). The initial NGS screening was supported in part by the contract research project for the Bundeswehr Medical Service FV E/U2AD/CF512/DF557 META-InfRisk and by European Union Horizon 2020 Research and Innovation Program COMPARE (grant no. 643476). Further NGS and rotavirus sequence analyses were supported by grants of the German Federal Ministry of Education and Research (project ESS-B.A.R., FKZ 13N13982) and the Deutsche Forschungsgemeinschaft, Germany (grant nos. JO369/4-3 and JO369/5-1).

References

  1. Tate  JEBurton  AHBoschi-Pinto  CParashar  UDWorld Health Organization–Coordinated Global Rotavirus Surveillance NetworkGlobal, regional, and national estimates of rotavirus mortality in children <5 years of age, 2000–2013. Clin Infect Dis2016;62(Suppl 2):S96105DOIExternal LinkPubMedExternal Link
  2. Attoui  HMertens  PPCBecnel  JBelaganahalli  SBergoin  MBrussaard  CPet al. Family: Reoviridae. In: King AMQ, Adams MJ, Carstens EB, Lefkowitz EJ, editors. Virus taxonomy: ninth report of the International Committee on Taxonomy of Viruses. Amsterdam: Elsevier Academic Press; 2012. p. 541–637.
  3. Bányai  KKemenesi  GBudinski  IFöldes  FZana  BMarton  Set al. Candidate new rotavirus species in Schreiber’s bats, Serbia. Infect Genet Evol2017;48:1926DOIExternal LinkPubMedExternal Link
  4. International Committee on Taxonomy of Viruses. Taxonomic information. Virus taxonomy: 2018b release. 2018 Jul [cited 2019 Aug 21]. https://talk.ictvonline.org/taxonomyExternal Link
  5. Matthijnssens  JOtto  PHCiarlet  MDesselberger  UVan Ranst  MJohne  RVP6-sequence-based cutoff values as a criterion for rotavirus species demarcation. Arch Virol2012;157:117782DOIExternal LinkPubMedExternal Link
  6. Matthijnssens  JCiarlet  MRahman  MAttoui  HBányai  KEstes  MKet al. Recommendations for the classification of group A rotaviruses using all 11 genomic RNA segments. Arch Virol2008;153:16219DOIExternal LinkPubMedExternal Link
  7. Martella  VBányai  KMatthijnssens  JBuonavoglia  CCiarlet  MZoonotic aspects of rotaviruses. Vet Microbiol2010;140:24655DOIExternal LinkPubMedExternal Link
  8. Schlegel  MRadosa  LRosenfeld  UMSchmidt  STriebenbacher  CLöhr  PWet al. Broad geographical distribution and high genetic diversity of shrew-borne Seewis hantavirus in Central Europe. Virus Genes2012;45:4855DOIExternal LinkPubMedExternal Link
  9. Fischer  SMayer-Scholl  AImholt  CSpierling  NGHeuser  ESchmidt  Set al. Leptospira genomospecies and sequence type prevalence in small mammal populations in Germany. Vector Borne Zoonotic Dis2018;18:18899DOIExternal LinkPubMedExternal Link
  10. Mayer-Scholl  AHammerl  JASchmidt  SUlrich  RGPfeffer  MWoll  Det al. Leptospira spp. in rodents and shrews in Germany. Int J Environ Res Public Health2014;11:756274DOIExternal LinkPubMedExternal Link
  11. Scheuch  MHöper  DBeer  MRIEMS: a software pipeline for sensitive and comprehensive taxonomic classification of reads from metagenomics datasets. BMC Bioinformatics2015;16:69DOIExternal LinkPubMedExternal Link
  12. Tausch  SHRenard  BYNitsche  ADabrowski  PWRAMBO-K: rapid and sensitive removal of background sequences from next generation sequencing data. PLoS One2015;10:e0137896DOIExternal LinkPubMedExternal Link
  13. Wang  XXLi  JXWang  WGSun  XMHe  CYDai  JJ[Preliminary investigation of viruses to the wild tree shrews (Tupaia belangeri Chinese)]Dongwuxue Yanjiu2011;32:669.PubMedExternal Link
  14. Li  KLin  XDHuang  KYZhang  BShi  MGuo  WPet al. Identification of novel and diverse rotaviruses in rodents and insectivores, and evidence of cross-species transmission into humans. Virology2016;494:16877DOIExternal LinkPubMedExternal Link
  15. Otto  PHRosenhain  SElschner  MCHotzel  HMachnowska  PTrojnar  Eet al. Detection of rotavirus species A, B and C in domestic mammalian animals with diarrhoea and genotyping of bovine species A rotavirus strains. Vet Microbiol2015;179:16876DOIExternal LinkPubMedExternal Link
Figures
Tables
Cite This Article

DOI: 10.3201/eid2512.191225
Original Publication Date: 11/6/2019

No comments:

Post a Comment