key: cord-0837414-idawavsd authors: Li, Kun; Lin, Xian-Dan; Huang, Kai-Yu; Zhang, Bing; Shi, Mang; Guo, Wen-Ping; Wang, Miao-Ruo; Wang, Wen; Xing, Jian-Guang; Li, Ming-Hui; Hong, Wang-Sheng; Holmes, Edward C.; Zhang, Yong-Zhen title: Identification of novel and diverse rotaviruses in rodents and insectivores, and evidence of cross-species transmission into humans date: 2016-07-31 journal: Virology DOI: 10.1016/j.virol.2016.04.017 sha: f39f4373a486587aaa540db5f18223e3db3780b8 doc_id: 837414 cord_uid: idawavsd Abstract Rotaviruses are an important cause of severe diarrheal illness in children globally. We characterized rotaviruses sampled in humans, insectivores (shrews) and rodents from urban and rural regions of Zhejiang province, China. Phylogenetic analyses revealed seven genotypic constellations of human rotaviruses with six different combinations of G and P genotypes – G3P[8] (50.06%), G9P[8] (36.16%), G1P[8] (8.92%), G2P[4] (4.63%), G3P[3] (0.12%), and G3P[9] (0.12%). In rodents and shrews sampled from the same locality we identified a novel genotype constellation (G32-P[46]-I24-R18-C17-M17-A28-N17-T19-E24-H19), a novel P genotype (P[45]), and two different AU-1-like rotaviruses associated with a G3P[3] genotype combination. Of particular note was a novel rotavirus from a human patient that was closely related to viruses sampled from rodents in the same region, indicative of a local species jump. In sum, these data are suggestive of the cross-species transmission of rodent rotaviruses into humans and for reassortment among human and animal rotaviruses. Rotaviruses are an important cause of severe diarrheal illness in infants and young children globally, causing 453,000 deaths per year in those o5 years of age Greenberg, 2013, Tate et al., 2008) . Rotaviruses are members of the genus Rotavirus of the family Reoviridae and possess a genome comprising 11 doublestranded RNA segments that encodes 6 structural (VP1-VP4, VP6, VP7) and 6 non-structural (NSP1-NSP6) proteins. Of the documented species (rotavirus A-H) and the tentative species (rotavirus I) (Mihalov-Kovács et al., 2015) , rotavirus A (RVA) is responsible for the majority of seasonal endemic diarrheal disease in young children. Based on sequences of the VP7 and VP4 genes RVAs have been further classified into G and P genotypes, respectively. To date, at least 31 G and 44 P genotypes have been identified worldwide (Matthijnssens et al., 2011 , Trojnar et al., 2013 Rotavirus Classification Working Group (RCWG), 2015) . Recently, a uniform system for rotavirus nomenclature was established for RVA based on nucleotide identities of the 11 rotavirus genome segments and phylogenetic analysis . Globally, the most common combinations in humans are G1, G2, G3, G4 or G9 and the P[4] or P[8] genotypes, especially G1P[8] which accounts for approximately 31% of human strains globally (Bányai et al., 2012 , Dóró et al., 2014 . Notably, genotype distributions vary among geographic regions, and several unusual strain combinations have recently been identified (Dóró et al., 2014 , Nordgren et al., 2012 . As well as humans, rotaviruses infect a wide range of vertebrates including domestic and wild mammals and birds (Estes and Greenberg, 2013) . Rodentia (rodents) is the largest order of (Wilson and Reeder, 2005) . It is therefore no surprise that rodents are also the largest zoonotic source of human infectious diseases (Luis et al., 2013 , Meerburg et al., 2009 . In addition to known pathogens, rodents are also an obvious source for the discovery of novel viruses, and new rodent arenaviruses, coronaviruses, and hantaviruses have been described recently (Gonzalez et al., 2007 . To date, however, relatively little attention has been directed toward the rotaviruses that might circulate in rodent populations (Everest et al., 2011 , Firth et al., 2014 , Greenwood & Sanchez, 2002 , Linhares et al., 1986 , Sachsenröder et al., 2014 . However, as rodents (especially rats) often live in close proximity to humans and domestic animals, and at high densities, they may play an important role in the cross-species transmission of rotaviruses, including to human populations. Zoonotic rotavirus infections in humans are not uncommon, and there are a growing number of reports describing the interspecies transmission of rotavirus among animals and from animals to humans (Ben Hadj Fredj et al., 2013 ., Gautam et al., 2015 , Martella et al., 2010 Medici et al., 2015 . Like influenza viruses, zoonotic rotaviruses can become increasingly "humanized" by reassortment with co-infecting human viruses (Cowley et al., 2013 , Jeong et al., 2014 , Matthijnssens et al., 2009 , Theuns et al., 2015 . Clearly, to better understand the evolution and emergence of rotaviruses it is important to determine the diversity, evolution and origins of rotaviruses in those animals that live in close proximity to humans, including rodents. Wenzhou is located on the southwestern coast of Zhejiang province, China, and is a prefecture-level city ( Figure S1 ). Wenzhou incorporates both urban and rural areas with a total population of 9.12 million, of which 1.69 million can be classed as pediatric. Wenzhou experiences a high level of rotavirus infections, totaling 465,000 cases each year. Longquan, also in Zhejiang province, is located 250 km west of Wenzhou, and is a countylevel city with a population of approximately 280,000. More than 90% of the Longquan's total area is mountainous. Compared with Wenzhou, Longquan experiences a relatively low level of rotavirus infection. In this study we investigated the diversity of rotaviruses in humans and small mammals (rodents and shrews) in Wenzhou and Longquan, as well as the evolutionary relationships between these viruses and those circulating in the local human population. Between October 2013-December 2014 a total of 1099 stool specimens were collected from diarrheal patients visiting the Diarrheal Department of the Second Affiliated Hospital and Yuying Children's Hospital of Wenzhou Medical University, Wenzhou city, Zhejiang province, China. Patient demographic data and clinical symptoms, including any complications, were collected and evaluated. Signed individual written informed consent for fecal sample collection was obtained from the guardian of each patient. Between March 2013-October 2014, 1865 small mammals (rodents and insectivores) were captured in Longwan, Lucheng, Ruian, and Wencheng counties/districts of Wenzhou city, as well as from Longquan city, as described previously (Fig. S1 ). All animals were treated according to the "Rules for Implementation of Laboratory Animal Medicine" from the Ministry of Health, China. Trapped animals were identified by morphological examination and further verified by sequence analysis of the cytochrome b (Cyt-b) gene . All animals were anesthetized before they were sacrificed, and every effort was made to minimize suffering. Stool samples were collected from each animal and used for the detection of rotaviruses. Total RNA was extracted from stool samples using the Bioteke fecal RNA isolation kit (Bioteke, Beijing, China) according to the manufacturer's instructions. All samples were screened for the presence of rotaviruses using nested RT-PCR. Two pairs of primers based on conserved regions of sequences of the VP7 segment of known group A rotaviruses were used (Table S1 ). RT-PCR was performed under the following conditions: incubation at 50°C for 30 min and 94°C for 3 min, 40 cycles of denaturation at 94°C for 40 s, annealing at 46°C for 40 s, and extending at 72°C for 60 s. The second amplification was performed as follows: incubation at 94°C for 3 min, 40 cycles of denaturation at 94°C for 40 s, annealing at 46°C for 40 s, and extending at 72°C for 50 s. Both the VP7 and VP4 genes were recovered from each of the rotavirus positive samples. Additionally, at least one complete genome sequence of each genotype was recovered by RT-PCR using primers based on conserved regions of known viruses. All primer sequences are described in Table S2 . PCR products purified using the QIAquick Gel Extraction kit (Qiagen) were sent to Shanghai Sangon Biological Engineering Technology and Services Co., Ltd. (Shanghai, China) for sequencing. They were sequenced with the PCR primers in both directions on an ABI Prism 3130 automated capillary sequencer (Applied Biosystems, Foster City, CA). Additionally, primer walking was performed to obtain the complete sequence of large genomic segments. Sequencing results were compared with the non-redundant nucleotide database at GenBank using Blastn to confirm they were of rotavirus origin. All rotavirus sequences obtained here from humans and animals have been deposited in GenBank and assigned accession numbers KU243375-KU243694 (Table S3 ). In addition to the sequences recovered here, reference sequences that cover the phylogenetic diversity of rotaviruses were compiled for evolutionary analyses. The following rotavirus data set sizes were used: VP7 ¼84 sequences; VP4 ¼70 sequences; VP1 ¼72 sequences; VP2 ¼73 sequences; VP3 ¼72 sequences; VP6 ¼76 sequences; NSP1 ¼77 sequences; NSP2 ¼76 sequences; NSP3 ¼75 sequences; NSP4 ¼73 sequences; NSP5 ¼74 sequences. The 5′ and 3′ terminal sequences of each segment complementary to the primers were excluded. All viral genome sequences obtained in this study were manually edited using the Seqman program implemented in the DNAStar v7.1 software package (DNASTAR, Inc., USA). Complete coding regions were aligned using the ClustalW method implemented in MEGA v5.05 (Tamura et al., 2011) . DNAStar was also used to calculate nucleotide and amino acid identities. Phylogenetic trees of these data were inferred using the maximum likelihood (ML) method implemented in the PhyML v3.0 package (Guindon et al., 2009 ). The best-fit nucleotide substitution model was determined using jModeltest (Posada 2008) . Accordingly, the VP1, VP2, VP3, NSP1 and NSP4 segments were analyzed using the HKY þIþ Γ model; the VP4, VP6, NSP2 and NSP3 segments were analyzed using the GTRþ Γ model; VP7 was analyzed using the HKYþ Γ model; while NSP5 was analyzed using the GTRþ Iþ Γ model. In addition, 1000 bootstrap replicates were used to assess the support for individual nodes on each tree. During 2013-2014 we studied a total of 1099 diarrheal patients who visited the Diarrheal Department of the Second Affiliated Hospital and Yuying Children's Hospital, Wenzhou city. Of these, 863 (78.53%) were confirmed as experiencing rotavirus infection based on RT-PCR targeting the viral VP7 gene. None of these patients were hospitalized. Of the patients with confirmed rotavirus infection, 99.07% were natives of Wenzhou, 0.58% were from surrounding areas, while 0.35% were from the "floating" (i.e. transient) population from other Chinese provinces. Their ages ranged from 2 months to 12 years and 66.5% were male. To better understand the genetic diversity of rotaviruses in human patients, nucleotide sequences of both the VP7 and VP4 gene segments were recovered from each of rotavirus positive samples (Table 1) . Based on phylogenetic analysis of the VP7 segment sequences (Fig. 1) , we determined that these infections were caused by RVAs at the following frequencies: G1 (8.92%), G2 (4.63%), G3 (50.06%), and G9 (36.16%). This analysis also revealed that G3 viruses (55.45%) were the most common in 2013, but decreased substantially in frequency (9.8%) in 2014, at which point G9 viruses became the predominant genotype (67.64%). Among the P genotypes, P To understand the role played by local small mammals in the evolution and cross-species transmission of rotaviruses, we performed a molecular epidemiological surveillance of rotaviruses in rodents and shrews in the same geographic region as the human cases described above. Accordingly, a total of 1865 small animals including 1416 rodents from 10 species and 449 insectivores of a single species (Suncus murinus, the Asian house shrew) were captured from Longquan and Wenzhou and screened by RT-PCR (Table S4 ). This resulted in the detection of rotaviruses in 7 species of rodents, comprising Rattus norvegicus (7/385, 1.82%), R. tanezumi (1/171, 0.59%), R. losea (1/200, 0.5%), Niviventer confucianus (1/109, 0.92%), Mus musculus (1/35, 2.86%), M. fortis (3/168, 1.79%), and Apodemus agrarius (1/332, 0.3%), with a total prevalence of 1.06% across all species. Similarly, 12 Asian house shrews tested positive for rotavirus RNA, with an overall prevalence of 2.67%. The nucleotide sequences of the VP7 and VP4 genes were also recovered from each of the rotavirus positive samples (Table S5) . Phylogenetic analysis of the VP7 genes revealed the presence of two rotavirus genotypes in rodents and insectivores; G3 (51.72%) and a novel genotype (48.28%) which the RCWG designated as genotype G32 (Fig. 1) Phylogenetic analysis revealed that all VP7 gene sequences from human rotaviruses fell into the G1, G2, G3, and G9 genotypes (Fig. 1) . Within the G1 genotype, those rotaviruses sampled from Wenzhou formed two separated lineages, showing a close evolutionary relationship to viruses sampled outside of China. In the G2 genotype, the viruses identified in Wenzhou were closely related to each other and formed a single distinct lineage, most closely related to human viruses from the USA. As noted previously (Matthijnssens et al., 2009) , the genotype G3 viruses did not form a monophyletic group. In the VP7 gene tree, the G3 viruses identified here and previously were classified into two groups, separated by other genotypes, including the newly identified genotype G32 and the G14. The first group included most human G3 viruses identified here. Notably, within this group, viruses could be further subdivided into two lineages. One lineage comprised the virus RVA/Human/CHN/WZ606/2013/G3P[9] that clustered with the human virus RVA/Human-tc/CHN/L621/2006/G3P[9] from China, as well as a number of animal (feline and porcine) rotaviruses, and therefore indicative of cross-species transmission. The second lineage consisted of only human viruses identified here and previously in and outside of China. The second group, which was more closely related to genotype G14, mainly comprised animal viruses from diverse species and small number of human viruses, again indicative of inter-host transmission. Of particular importance was that the virus RVA/Human-wt/CHN/WZ101/2013/ G3P[3] from Wenzhou was closely related to those viruses identified in mice and rats from neighboring Longquan, compatible with the idea that the human virus originated in the local animal population. All VP7 gene sequences from rodents and shrews were classified into four groups in the VP7 gene tree (Fig. 1 ). The first group included the viruses identified in shrews and mice in Wenzhou. Notably, these viruses comprised a novel genotype -G32that was clearly distinct from known G genotypes ( o82% nucleotide similarity). Within the second group, the virus RVA/Rat-wt/CHN/ RA108/2013/G3P[3] identified in a rat from Ruian was the most distinct, falling in a basal position. The remaining rotaviruses sampled from mice and rats in Wenzhou and Longquan formed two lineages. The viruses RVA/Rat-wt/CHN/WC179/2013/G3P [45] and RVA/Rat-wt/CHN/RA116/2013/G3P[45] identified in rats from Ruian and Wencheng were closely related to each other and clustered with viruses previously identified in horses, rats, cats, and humans. As described above, the rotaviruses sampled from mice and rats in Longquan were closely related to the human virus RVA/Human-wt/CHN/WZ101/2013/G3P [3] . Further analysis of the antigenic regions of VP7 revealed that there were no amino acid differences between the human and rodent viruses (Table S6 ). In the VP4 gene tree (Fig. 1) genotype, and showed a close relationship with a rat virus previously sampled in Germany (Sachsenröder et al., 2014) . To better characterize those rotaviruses identified in humans and small mammals in Zhejiang province, the remaining segments of representative of each virus were sequenced, from which we were able to define their full genotype constellations (Table 2) . Three genotype constellations were identified within the human viruses: (i) Wa-like, (ii) Wa-like x DS-1-like reassortant, (iii) AU-1like (Table 2 ). All the Wa-like viruses identified here possessed the same genomic backbone (I1-R1-C1-M1-A1-N1-T1-E1-H1) (genotypes identified in the trees shown in Figs. 2-4) . To date, only a small number of human Wa-like x DS-1-like reassortants have been described . Interestingly, although we did not identify any DS-1-like viruses, two genetic constellations of Wa-like x DS-1-like reassortants were observed, with the H1 genotype identified in three viruses and the C1, N1, T1, and H1 genotypes identified in another. Notably, one AU-1-like strain possessed the A9 and H6 genotypes. In the phylogenetic trees estimated for the VP1-VP3 and VP6 genes (Fig. 2 ) and the NSP1-NSP5 genes (Figs. 3 and 4) , the human Wa-like viruses (I1-R1-C1-M1-A1-N1-T1-E1-H1) from Wenzhou were closely related to those previously identified in humans from a variety of geographic locations. Notably, they only clustered together in the NSP5 gene tree (Fig. 4) and were separated by viruses sampled from other locations in the remaining 8 gene trees (Figs. 2 and 3) . Hence, these data reveal that there have been multiple independent introductions of these viruses into Wenzhou. Both types of Wa-like x DS-1-like reassortant viruses (Table 1) clustered together in the NSP4 gene tree and with viruses with the genotype constellation (G9-P[4]-I2-R2-C2-M2-A2-N2-T2-E6-H2) in the NSP2 gene tree. However, they exhibited distinct topologies were most closely related to those identified in local rodents that showed the same reassortment pattern (Table 2) . Finally, five segments of the AU-like virus RVA/Human-wt/CHN/ WZ606/2013/G3P[9] were closely related to those viruses sampled from animals, while the remaining four segments were closely related to human viruses from China (Wang et al., 2013a,b). Hence, these very close phylogenetic relationships clearly indicate the interspecies transmission of RVA from animals into humans. A variety of genotype constellations were characterized in the animal rotaviruses sampled here. AU-like viruses were found in both the mouse and rat viruses sampled from Longquan and Wenzhou; they possessed at least two different genotype constellations (I3-R3-C3-M3-A9-N3-T3-E3-H6) and (I3-R3-C3-M10-A22-N3-T3-E3-H13) that were associated with G3P[3] and G3P [45] combinations, respectively. Importantly, rodent viruses with the genotype constellation I3-R3-C3-M3-A9-N3-T3-E3-H6 were closely related to the human virus RVA/Human-wt/CHN/WZ101/ 2013/G3P [3] in all 11 segments , with 94.7% to 99.4% sequence similarity. Also of note was that a novel genotype constellation (G32-P[46]-I24-R18-C17-M17-A28-N17-T19-E24-H19) was identified in shrews and mice in Wenzhou, and was distinct from the known rotaviruses in all 11 gene trees (Figs. 2-4) , with genetic differences ranging from 74% to 89%. Rotavirus-associated diarrheal disease is a major public health problem in China, with several million cases reported in children each year (Liu et al., 2014 , Nan et al., 2014 . The co-circulation of different G and P genotypes of human RVA has been described in multiple regions globally (Bányai et al., 2012 , Dóró et al., 2014 , Liu et al., 2014 , Nan et al., 2014 , Santos & Hoshino, 2005 . We observed high levels of genetic diversity in rotaviruses sampled in Zhejiang province. Indeed, analysis of the 11 genome segments revealed the co-circulation of Wa-like, Wa-like x DS-1 reassortant, and AU-1-like viruses in patients from a single hospital in Wenzhou. Globally, genotypes G1, G2, G3, G4, G9, and G12 combined with P[4], P[6], and P[8] are commonly detected in humans, with G1P [8] the most prevalent (31%) in the pre-rotavirus vaccine era (Bányai et al., 2012) . However, the dominant G and P genotypes, as well as their combinations, exhibit both temporal and geographic variation, even on a restricted scale as in this study . In China, the most commonly detected G and P types are G1 (4 30%) and G3 (4 30%), and P[8] (450%), with the combination G3P[8] (432%) the most prevalent (Liu et al., 2014; Nan et al., 2014; . Similar to previous studies, G3P [8] was most commonly sampled in humans during 2013-2014 in Wenzhou. However, we observed an increasing incidence of G9 (to 67.64% in 2014). One of the most striking observations of our study was the diversity of rotaviruses found in rodents and insectivores in Zhejiang, including a novel genotype constellation identified in rodent and shrews that was distinct in all 11 gene trees (Figs. 1-4) . According to the criteria for genotype demarcation in the genus Rotavirus proposed by the RCWG (Matthijnssens et al., 2011) , this virus is sufficiently genetically distinct that it should be recognized as a distinct genotype constellation. Furthermore, one novel G (G32) and two novel P genotypes (P[45] and P[46]) were identified in rodents captured in Wenzhou, with 420% nucleotide differences from any recognized VP7 and VP4 genes. Finally, at least two different AU-1-like viruses were circulating in rodents and insectivores in this geographic region ( Table 2) . The AU-1-like genotype constellation (I3-R3-C3-M3-A3/A12-N3-T3-E3-H3/H6) is believed to be of feline/canine RVA origin , and AU-1-like reassortants have also been reported in cats (Matthijnssens et al., 2011) . Overall, compared with Wa-like and DS-1-like viruses, AU-1-like and AU-1-like reassortant viruses are at low prevalence in humans . Recently, the genotype constellation G3-P[3]-I8-R3-C3-M3-A9-N3-T3-E3-H6 was identified in bats sampled from China (He et al., 2013) . In this study, three types of AU-1-like genotype constellations were identified in rodents ( Table 2) , indicative of a diverse host range of AU-1-lke genotype constellations with multiple reassortment events involving bat, feline species, and rodents. The potential for novel animal rotaviruses to cross species boundaries and infect humans is clearly of public health importance, and comparable to the process that commonly occurs with hantaviruses and influenza viruses Wang et al., 2013a,b; Yoon et al., 2014) . Strikingly, our phylogenetic analyses identified one human virus (RVA/Human-wt/CHN/WZ101/ 2013/G3P[3]) that was very closely related to those rotaviruses sampled from rodents in the same geographic region, strongly suggestive of animal to human cross-species transmission. Despite the relatively low prevalence of RVAs in rodents and insectivores, such that the likelihood of zoonotic transmission from these animals may be small, additional work to determine whether this rodent virus has infected other humans in this part of China will be important. Because of the close relationship between RVA/Humanwt/CHN/WZ101/2013/G3P[3] and those viruses found in bats, dogs, cats and rodents, this may represent a true animal rotavirus genotype constellation. In addition, although it clustered with human AU-1-like viruses in nine of the segments (G3, P[3], I3, R3, C3, M3, N3, T3, E3) ( Table 2) , its NSP1 and NSP5 clearly fell into the A9 and A6 groups, respectively, suggesting that it may be a reassortant. Previous studies have revealed the presence of multiple RVAs in rodents (Everest et al., 2011; Firth et al., 2014; Greenwood and Sanchez, 2002; Linhares et al., 1986; Sachsenröder et al., 2014) . Our study is of note because rat, mouse and shrew populations were found to harbor (i) known genotypes of RVAs, (ii) novel G and P genotypes, and (iii) a novel genotype constellation. In sum, these data suggest that rodents and insectivores may constitute an important reservoir for rotaviruses, including their potential crossspecies transmission to humans. As rodents and shrews often live in close contact with humans, as well as with farm and companion animals, they clearly provide an important nexus between wildlife communities and human populations. Systematic review of regional and temporal trends in global rotavirus strain diversity in the prerotavirus vaccine era: insights for understanding the impact of rotavirus vaccination programs Feline origin of rotavirus strain Complete genome analysis of a rabbit rotavirus causing gastroenteritis in a human infant Review of global rotavirus strain prevalence data from six years post vaccine licensure surveillance: is there evidence of strain selection from vaccine pressure? Rotaviruses, 6th ed Rotavirus in a wild English red squirrel (Sciurus vulgaris) identified by electron microscopy Detection of zoonotic pathogens and characterization of novel viruses carried by commensal Rattus norvegicus Full genomic characterization and phylogenetic analysis of a zoonotic human G8P[14] rotavirus strain detected in a sample from Guatemala Serological evidence of murine pathogens in wild grey squirrels (Sciurus carolinensis) in North Wales Estimating maximum likelihood phylogenies with PhyML Phylogeny and origins of hantaviruses harbored by bats, insectivores, and rodents. Plos. Pathog. 9 Characterization of a novel G3P[3] rotavirus isolated from a lesser horseshoe bat: a distant relative of feline/canine rotaviruses The evolution and emergence of hantaviruses Whole-genome analysis of a rare human Korean G3P[9] rotavirus strain suggests a complex evolutionary origin potentially involving reassortment events between feline and bovine rotaviruses Isolation and characterization of a novel arenavirus harbored by Rodents and Shrews in Zhejiang province Diversity of rotavirus strains causing diarrhea in o 5 years old Chinese children: a systematic review Rotavirus infection in wild marsupials (Didelphis marsupialis) of the Amazon region Update on the disease burden and circulating strains of rotavirus in China: a systematic review and meta-analysis A comparison of bats and rodents as reservoirs of zoonotic viruses: are bats special? Zoonotic aspects of rotaviruses Full genome-based classification of rotaviruses reveals a common origin between human Wa-like and porcine rotavirus strains and human DS-1-like and bovine rotavirus strains Uniformity of rotavirus strain nomenclature proposed by the Rotavirus Classification Working Group (RCWG) VP6-sequence-based cutoff values as a criterion for rotavirus species demarcation Reassortment of human rotavirus gene segments into G11 rotavirus strains Are human P[14] rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla? Genotype constellation and evolution of group A rotaviruses infecting humans Genetic diversity in three bovine-like human G8P[14] and G10P[14] rotaviruses suggests independent interspecies transmission events Rodent-borne diseases and their risks for public health Candidate new rotavirus species in sheltered dogs Epidemiological and clinical studies of rotavirus-induced diarrhea in China from 1994-2013 jModelTest: phylogenetic model averaging Rotavirus Classification Working Group (RCWG). 2015. Newly assigned genotypes Metagenomic identification of novel enteric viruses in urban wild rats and genome characterization of a group A rotavirus Global distribution of rotavirus serotypes/genotypes and its implication for the development and implementation of an effective rotavirus vaccine MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance and maximum parsimony methods WHO-coordinated Global Rotavirus Surveillance Network., 2012 estimate of worldwide rotavirus-associated mortality in children younger than 5 years before the introduction of universal rotavirus vaccination programmes: a systematic review and meta-analysis Complete genome characterization of recent and ancient Belgian pig group A rotaviruses and assessment of their evolutionary relationship with human rotaviruses Identification of an avian group A rotavirus containing a novel VP4 gene with a close relationship to those of mammalian rotaviruses Discovery, diversity and evolution of novel coronaviruses sampled from rodents in China Ongoing spillover of Hantaan and Gou hantaviruses from rodents is associated with hemorrhagic fever with renal syndrome (HFRS) in China Complex evolutionary patterns of two rare human G3P[9] rotavirus strains possessing a feline/canine-like H6 genotype on an AU-1-like genotype constellation Mammal Species of the World. A Taxonomic and Geographic Reference Evolution and ecology of influenza A viruses Rotavirus incidence and genotype distribution before and after national rotavirus vaccine introduction in Belgium Rotavirus-Specific and Overall Diarrhea Mortality in Chinese Children Younger than 5 Years Evidence for interspecies transmission and reassortment events This work was supported by the National Natural Science Supplementary data associated with this article can be found in the online version at http://dx.doi.org/10.1016/j.virol.2016.04.017.