Figures
Abstract
G12 rotaviruses are emerging rotavirus strains causing severe diarrhea in infants and young children worldwide. However, the whole genomes of only a few G12 strains have been fully sequenced and analyzed. In this study, we sequenced and characterized the complete genomes of six G12 strains (RVA/Human-tc/MMR/A14/2011/G12P[8], RVA/Human-tc/MMR/A23/2011/G12P[6], RVA/Human-tc/MMR/A25/2011/G12P[8], RVA/Human-tc/MMR/P02/2011/G12P[8], RVA/Human-tc/MMR/P39/2011/G12P[8], and RVA/Human-tc/MMR/P43/2011/G12P[8]) detected in six stool samples from children with acute gastroenteritis in Myanmar. On whole genomic analysis, all six Myanmarese G12 strains were found to have a Wa-like genetic backbone: G12-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1 for strains A14, A25, P02, P39, and P43, and G12-P[6]-I1-R1-C1-M1-A1-N1-T1-E1-H1 for strain A23. Phylogenetic analysis showed that most genes of the six strains examined in this study were genetically related to globally circulating human G1, G3, G9, and G12 strains. Of note is that the NSP4 gene of strain A23 exhibited the closest relationship with the cognate genes of human-like bovine strains as well as human strains, suggesting the occurrence of reassortment between human and bovine strains. Furthermore, strains A14, A25, P02, P39, and P43 were very closely related to one another in all the 11 gene segments, indicating derivation of the five strains from a common origin. On the other hand, strain A23 consistently formed distinct clusters as to all the 11 gene segments, indicating a distinct origin of strain A23 from that of strains A14, A25, P02, P39, and P43. To our knowledge, this is the first report on whole genome-based characterization of G12 strains that have emerged in Myanmar. Our observations will provide important insights into the evolutionary dynamics of spreading G12 rotaviruses in Asia.
Citation: Ide T, Komoto S, Higo-Moriguchi K, Htun KW, Myint YY, Myat TW, et al. (2015) Whole Genomic Analysis of Human G12P[6] and G12P[8] Rotavirus Strains that Have Emerged in Myanmar. PLoS ONE 10(5): e0124965. https://doi.org/10.1371/journal.pone.0124965
Academic Editor: Karol Sestak, Tulane University, UNITED STATES
Received: January 20, 2015; Accepted: March 20, 2015; Published: May 4, 2015
Copyright: © 2015 Ide 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 nucleotide sequence data presented in this paper have been deposited in the DDBJ and EMBL/GenBank data libraries. The accession numbers for the nucleotide sequences of the VP1-4, VP6-7, and NSP1-5 genes of strains A14, A23, A25, P02, P39, and P43 are LC019041-LC019051, LC019052-LC019062, LC019063-LC019073, LC019074-LC019084, LC019085-LC019095, and LC019096-LC019106, respectively.
Funding: This study was supported in part by the MEXT-Supported Program for the Strategic Research Foundation at Private Universities, 2010-2014 (KT), and Grants-in-Aid for Research on Emerging and Re-emerging Infectious Diseases from the Ministry of Health, Labor, and Welfare of Japan (KT). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Group A rotavirus (RVA), a member of the family Reoviridae, is the major cause of severe dehydrating diarrhea in the young of humans and many animal species worldwide. In humans, RVA infections are associated with high morbidity and mortality, being responsible for an estimated annual 453,000 deaths among children <5 years of age [1]. More than half of these deaths were estimated to occur in developing countries in Asia and Africa [1, 2]. In Myanmar, a high disease burden of RVA infection in children (<5 years old) was found in studies conducted in Yangon, RVAs being identified in 53–57% of fecal samples from diarrheal disease cases in the years 2002–2005 [3, 4].
The RVA virion is a triple-layered, non-enveloped icosahedron that encloses an 11-segment genome of double-stranded (ds)RNA [5]. RVA has two outer capsid proteins, VP7 and VP4, which are implicated independently in neutralization, and define the G and P genotypes, respectively. To date, RVAs have been classified into at least 27 G and 37 P genotypes [6, 7]. Among them, 5 G (G1-4 and G9) and 3 P (P[4], P[6], and P[8]) genotypes are commonly associated with human infections [8, 9]. Over the last decade, the global emergence of unusual G12 strains has been a matter of concern, and G12 seems to be the sixth major human G genotype [9–11].
The first G12 strain, L26 (G12P[4]), was detected in children <2 years old with acute diarrhea in the Philippines in 1987 [12, 13]. More than 10 years later, G12 strains began to emerge globally, predominantly in combination with either a P[6] or P[8] genotype and less commonly with a P[4] or P[9] genotype [9–11, 14–17]. In Asia, G12 strains are being increasingly identified in diarrheic children in several countries; Bangladesh [10, 18–22], Bhutan [23], India [24–33], Indonesia [34], Japan [16], Kazakhstan [35], Korea [36], Kyrgyzstan [35], Nepal [37–40], Pakistan [41, 42], Sri Lanka [43], Thailand [14, 15, 44–46], and Vietnam [47], indicating the continued spread of G12 strains in Asia. In 2011, we detected 29 G12 strains from diarrheic children in Central Myanmar enrolling a total of 54 RVA-positive stool specimens by PCR-based G and P genotyping (Ide et al., in preparation).
A whole genome-based genotyping system was recently proposed for RVAs based on assignment to all the 11 gene segments (i.e., G/P and non-G/P genes) [48]. In the new genotyping system, the acronym Gx-P[x]-Ix-Rx-Cx-Mx-Ax-Nx-Tx-Ex-Hx, where x is an integer, defines the genotype of the VP7-VP4-VP6-VP1-VP2-VP3-NSP1-NSP2-NSP3-NSP4-NSP5 genes of a given RVA strain. The Wa-like strains are characterized by non-G/P genotypes (I1-R1-C1-M1-A1-N1-T1-E1-H1), and tend to have G/P genotypes G1P[8], G3P[8], G4P[8], and G9P[8] [49]. In contrast, the DS-1-like strains are characterized by non-G/P genotypes (I2-R2-C2-M2-A2-N2-T2-E2-H2), and tend to have G/P genotype G2P[4]. The third minor AU-1-like strains are characterized by non-G/P genotypes (I3-R3-C3-M3-A3-N3-T3-E3-H3), and tend to have G/P genotype G3P[9]. Whole genome-based analysis is a reliable method for obtaining conclusive data on the origin of an RVA strain, and for tracing its evolutionary pattern [48, 50]. To date, the whole genome sequences of only a few G12 strains have been fully sequenced and characterized, which provided evidence that the recently detected Asian G12 strains with the Wa-like genotype backbone are distantly related to prototype Asian G12 strains L26 (G12P[4]) and T152 (G12P[9]) [14] having the DS-1-like and AU-1-like genotype backbone, respectively [10]. However, as the overall genomic constellation and the exact evolutionary pattern of Asian G12 strains remain to be elucidated, whole genomic analysis of Myanmarese G12 strains might be useful for obtaining a more precise understanding of the evolutionary pattern of emerging G12 strains in Asia. In the present study, we analyzed the whole genomes of six G12 strains that have emerged in Myanmar. Moreover, in this study, deep sequencing with the next generation sequencing (NGS) Illumina MiSeq platform was carried out to obtain the complete nucleotide sequences of the whole genomes of these six G12 strains.
Materials and Methods
Ethics statement
The study was approved by the Ethical Committee on Medical Research Involving Human Subjects of the Department of Medical Research (Central Myanmar). In this study, written informed consent for the testing of stool samples for RVAs and characterization of identified RVA strains was obtained from children’s parents/guardians.
Virus strains
The full-genomic sequences were determined for strains A14, A23, A25, P02, P39, and P43, which were identified as the pathogens causing diarrhea in six stool specimens from children with acute diarrhea during the RVA strain surveillance in the Pediatric Infectious Disease Wards of the Defense Services Obstetrics, Gyneacology and Children’s Hospital, Central Myanmar in 2011, which involved a total of 54 RVA-positive fecal samples (Ide et al., in preparation). Stool samples containing strains A14, A23, A25, P02, P39, and P43 were kept at −30°C until use.
Virus isolation
Stool samples suspended in PBS containing 0.5 mM MgCl2 and 1 mM CaCl2 were inoculated onto monkey kidney cell line MA104 for virus isolation [51], and the cultures were serially passaged two more times in MA104 cells. The viral dsRNAs were extracted from the cell cultures using a QIAamp Viral RNA Mini Kit (Qiagen). The extracted dsRNAs were used for (i) polyacrylamide gel electrophoresis (PAGE) analysis, and (ii) whole genomic analysis. For PAGE analysis, the dsRNAs were electrophoresed in a 10% polyacrylamide gel for 16 h at 20 mA at room temperature, followed by silver staining [52] to determine the genomic dsRNA profiles. For whole genomic analysis, viral dsRNAs were subjected to Illumina MiSeq sequencing as described below.
cDNA library building and Illumina MiSeq sequencing
Preparation of a cDNA library and Illumina MiSeq sequencing were performed as reported previously [49, 53]. Briefly, a 200 bp fragment library ligated with bar-coded adapters was constructed for individual strains using an NEBNext Ultra RNA library Prep Kit for Illumina v1.2 (New England Biolabs) according to the manufacturer’s instructions. Library purification was performed using Agencourt AMPure XP magnetic beads (Beckman Coulter). The quality of the purified cDNA libraries was assessed on a MultiNA MCE-202 bioanalyzer (Shimadzu). Nucleotide sequencing was performed on an Illumina MiSeq sequencer (Illumina) using a MiSeq Reagent Kit v2 (Illumina) to generate 151 paired-end reads. Data analysis was carried out using CLC Genomics Workbench v7.0.3 (CLC Bio). Contigs were assembled from the obtained sequence reads by de novo assembly. Using the assembled contigs as query sequences, the Basic Local Alignment Search Tool (BLAST) non-redundant nucleotide database was searched to obtain the full-length nucleotide sequence of each gene segment of the six strains. To fill in missing sequence gaps for the NSP2 gene of strain A23, and the NSP5 genes of strains A14 and A25, specific primers were used to process viral dsRNAs of the three strains using Sanger sequencing [54]. The primer pairs used for cDNA amplification of the NSP2 and NSP5 genes are (+) 5’-GGCTTTTAAAGCGTCTCAGTC-3’ and (-) 5’-GGTCACATAAGCGCTTTCTATTC-3’, and (+) 5’-GGCTTTTAAAGCGCTACAGTG-3’ and (-) 5’-GGTCACAAAACGGGAGTGGGG-3’, respectively [54]. The final full-length genomes for the six strains were assembled using a combination of Illumina MiSeq and Sanger sequencing data. The nucleotide sequences were translated into amino acid sequences using GENETYX v11 (GENETYX).
SNP detection
Single nucleotide polymorphisms (SNPs) were called using CLC Genomics Workbench v7.0.3 using all available sequencing data with at least 100 sequence reads covering each nucleotide position. A variant frequency threshold of ≥1% per site was used and approximate variant P-values were calculated.
Determination of RVA genotypes
The genotype of each of the 11 gene segments of the six strains was determined using the RotaC v2.0 automated genotyping tool (http://rotac.regatools.be/) [55] according to the guidelines proposed by the Rotavirus Classification Working Group (RCWG).
Phylogenetic analyses
Sequence comparisons were carried out as described previously [54]. Briefly, multiple alignment of each viral gene was performed using CLUSTAL W [56]. Phylogenetic trees were constructed using the maximum likelihood method and the Kimura 2-parameter substitution model using MEGA6.06 [57]. The reliability of the branching order was estimated from 1000 bootstrap replicates [58]. The results of phylogenetic analyses were validated using several other genetic distance models, such as the Jukes-Cantor, Tamura 3-parameter, Hasegawa-Kishino-Yano, and Tamuta-Nei ones (data not shown).
Nucleotide sequence accession numbers
The nucleotide sequence data presented in this paper have been deposited in the DDBJ and EMBL/GenBank data libraries. The accession numbers for the nucleotide sequences of the VP1-4, VP6-7, and NSP1-5 genes of strains A14, A23, A25, P02, P39, and P43 are LC019041-LC019051, LC019052-LC019062, LC019063-LC019073, LC019074-LC019084, LC019085-LC019095, and LC019096-LC019106, respectively.
Results and Discussion
Isolation of strains A14, A23, A25, P02, P39, and P43 in cell culture
To obtain enough genomic material for whole genome sequencing of the emerging G12 strains in Myanmar, we primarily attempted to isolate strains A14, A23, A25, P02, P39, and P43 using the MA104 cell line; all six strains could be cell culture-adapted. Virion dsRNAs were extracted and then analyzed by PAGE. Fig 1 shows the profiles of the viral dsRNAs from human strain KU (G1P[8]), as a reference (lane 1), and strains A14 (lane 2), A23 (lane 3), A25 (lane 4), P02 (lane 5), P39 (lane 6), and P43 (lane 7) from the cell cultures. The individual dsRNA migration pattern in PAGE of cell culture-adapted strains A14, A23, A25, P02, P39, and P43 was identical to that of the original strains A14, A23, A25, P02, P39, and P43 present in stool samples, respectively (data not shown). They all showed a long electropherotype. Cell culture-adapted strains A14, A23, A25, P02, P39, and P43 were named RVA/Human-tc/MMR/A14/2011/G12P[8], RVA/Human-tc/MMR/A23/2011/G12P[6], RVA/Human-tc/MMR/A25/2011/G12P[8], RVA/Human-tc/MMR/P02/2011/G12P[8], RVA/Human-tc/MMR/P39/2011/G12P[8], and RVA/Human-tc/MMR/P43/2011/G12P[8], respectively, according to the guidelines for the uniformity of RVAs proposed by the RCWG. Of note was that strains A14, A25, P02, P39, and P43 showed an almost identical electropherotype, suggesting a close genetic relatedness among the five strains.
Lane 1, dsRNAs of strain KU (G1P[8]) extracted from a cell culture; lanes 2–7, dsRNAs of strains A14 (lane 2), A23 (lane 3), A25 (lane 4), P02 (lane 5), P39 (lane 6), and P43 (lane 7) extracted from cell cultures. The numbers on the left indicate the order of the genomic dsRNA segments of strain KU.
https://doi.org/10.1371/journal.pone.0124965.g001
Nucleotide sequencing and whole-genome-based genotyping of strains A14, A23, A25, P02, P39, and P43
In order to gain an insight into the genetic variability among strains A14, A23, A25, P02, P39, and P43, and the genetic relatedness with other RVA strains worldwide, the full-genome sequences of all the 11 segments of these six strains were determined using the NGS Illumina MiSeq platform. The whole genomes of the six strains were amplified using a sequence-independent primer set and then sequenced successfully. Illumina MiSeq sequencing yielded 2.0 x 106 reads (~143 bp average length), 1.9 x 106 reads (~144 bp average length), 1.7 x 106 reads (~144 bp average length), 1.8 x 106 reads (~145 bp average length), 1.7 x 106 reads (~142 bp average length), and 1.7 x 106 reads (~145 bp average length) for strains A14, A23, A25, P02, P39, and P43, respectively. Missing sequence information in the contigs of the NSP2 gene of strain A23, and the NSP5 genes of strains A14 and A25 was filled in by Sanger sequencing of RT-PCR products. This approach allowed determination of the complete or nearly complete nucleotide sequences of all the 11 segments of the six strains, while deep sequencing revealed the presence of multiple SNPs at <1.8% frequency in the six strains.
The 11 genes of strains A14, A23, A25, P02, P39, and P43 were assigned as G12-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1 (strains A14, A25, P02, P39, and P43), and G12-P[6]-I1-R1-C1-M1-A1-N1-T1-E1-H1 (strain A23) (Fig 2). Strains A14, A23, A25, P02, P39, and P43 were confirmed to have the G12P[8], G12P[6], G12P[8], G12P[8], G12P[8], and G12P[8] genotypes, respectively, as determined by PCR-based genotyping (Ide et al., in preparation). Comparison of the complete genotype constellations of strains A14, A23, A25, P02, P39, and P43 with those of other G12 and non-G12 strains is shown in Fig 2. All the six Myanmarese G12 strains exhibited typical Wa-like genotype constellations, which are commonly found in the G12 strains recently detected worldwide (Rahman et al., 2007). Furthermore, as suggested by the genomic dsRNA profiles observed on PAGE analysis (Fig 1), strains A14, A25, P02, P39, and P43 exhibited extremely high nucleotide sequence identities (99.1–100%) to one another for all the 11 gene segments. On the other hand, the nucleotide sequence similarities of the VP4 and other 10 gene segments (VP7, VP6, VP1-3, and NSP1-5) of strain A23 to those of strains A14, A25, P02, P39, and P43 were comparatively low (75.3 and 93.1–98.7%, respectively).
Strains A14, A23, A25, P02, P39, and P43 are shown in red. Gray shading indicates the 10 gene segments (VP7, VP6, VP1-3, and NSP1-5) with genotypes identical to those of strains A14, A23, A25, P02, P39, and P43. Green shading indicates the VP4 gene segments with a P[8] genotype identical to those of strains A14, A25, P02, P39, and P43. Blue shading indicates the VP4 gene segments with a P[6] genotype identical to that of strain A23. “−” indicates that no sequence data were available in the DDBJ and EMBL/GenBank data libraries. aThe gene segments that are most similar to those of strain A14. bThe gene segments that are most similar to those of strain A23. cThe gene segments that are most similar to those of strain A25. dThe gene segments that are most similar to those of strain P02. eThe gene segments that are most similar to those of strain P39. fThe gene segments that are most similar to those of strain P43.
https://doi.org/10.1371/journal.pone.0124965.g002
Phylogenetic analyses
We next constructed phylogenetic trees using the full-genome sequence for each of the 11 gene segments because phylogenetic analysis of RVA nucleotide sequences provides direct evidence of their relatedness to those of other strains, even within the same genotype [48].
The VP7 genes of strains A14, A25, P02, P39, and P43 exhibited the maximum nucleotide sequence identities (99.7–99.8%) with that of Thai human strain CU331-NR (G12P[6]) [45] (Fig 2). On phylogenetic analysis, strains A14, A25, P02, P39, and P43 were found to form a cluster with strain CU331-NR and several human G12 strains from Thailand in G12 lineage-3, in which the majority of globally circulating G12 strains cluster [9, 10] (Fig 3). On the other hand, the VP7 gene of strain A23 showed the highest nucleotide sequence similarity (99.1%) with South African human strain MRC-DPRU4090 (G12P[6]) (Fig 2). Phylogenetically, strain A23 was found to be closely related with strain MRC-DPRU4090 in a common branch with German human strain GER172-08 (G12P[6]) [59] in G12 lineage-3 (Fig 3).
In the tree, the positions of strains A14, A23, A25, P02, P39, and P43 are shown in red. Bootstrap values of <75% are not shown. Scale bars, 0.05 substitutions per nucleotide.
https://doi.org/10.1371/journal.pone.0124965.g003
The VP4 genes of strains A14, A25, P02, P39, and P43 showed the highest nucleotide sequence identities (99.6%) with the cognate genes of Thai human strain CU769-KK (G1P[8]) [46], Russian human strain Nov09-B28 (G1P[8]), and/or Canadian human strain RT010-09 (G3P[8]) (Fig 2). On phylogenetic analysis, strains A14, A25, P02, P39, and P43 clustered near these, and several human G1 and G3 strains from different countries of the world in P[8] lineage-3 (Fig 4). In contrast, the VP4 gene of strain A23 exhibited the maximum nucleotide sequence identity (99.2%) with South African human strain MRC-DPRU4090 (G12P[6]) (Fig 2). Phylogenetically, strain A23 was found to be closely related with strain MRC-DPRU4090 in a common branch with several human G12 strains from different parts of the world in P[6] lineage-2 (Fig 4).
See legend of Fig 3. Scale bars, 0.05 substitutions per nucleotide.
https://doi.org/10.1371/journal.pone.0124965.g004
The VP6 genes of strains A14, A25, P02, P39, and P43 exhibited the highest nucleotide sequence identities (99.3–99.6%) with the VP6 gene of Australian human strain CK00089 (G1P[8]) (Fig 2). On phylogenetic analysis, strains A14, A25, P02, P39, and P43 were found to be closely related with strain CK00089 near several human G12 strains (Fig 5). On the other hand, the VP6 gene of strain A23 showed the maximum nucleotide similarity (99.6%) with South African human strain MRC-DPRU4090 (G12P[6]) (Fig 2). Phylogenetically, strain A23 was found to be clustered with strain MRC-DPRU4090 and several human G1 strains (Fig 5).
See legend of Fig 3. Scale bars, 0.02 substitutions per nucleotide.
https://doi.org/10.1371/journal.pone.0124965.g005
The VP1 genes of strains A14, A25, P02, P39, and P43 showed the maximum nucleotide sequence identities (99.6–99.7%) with the cognate gene of Australian human strain CK00089 (G1P[8]) (Fig 2). On phylogenetic analysis, strains A14, A25, P02, P39, and P43 were found to be clustered near strain CK00089, and several human G1, G11, and G12 strains (Fig 6). On the other hand, the VP1 gene of strain A23 exhibited the highest nucleotide sequence similarity (99.1%) with South African human strains MRC-DPRU2130-05 (G12P[6]) and MRC-DPRU840 (G1P[8]) (Fig 2). Phylogenetically, strain A23 clustered near these, and several human G1 and G12 strains (Fig 6).
See legend of Fig 3. Scale bars, 0.05 substitutions per nucleotide.
https://doi.org/10.1371/journal.pone.0124965.g006
The VP2 genes of strains A14, A25, P02, P39, and P43 showed the highest nucleotide sequence similarities (99.6–99.7%) with the VP2 genes of Thai human strains CU331-NR (G12P[6]) and/or CU460-KK (G12P[8]) [45] (Fig 2). On phylogenetic analysis, strains A14, A25, P02, P39, and P43 were found to be closely related with these strains in a common branch with South African human strain MRC-DPRU2132 (G1P[8]) (Fig 7). In contrast, the VP2 gene of strain A23 exhibited the maximum nucleotide sequence identity (99.4%) with South African human strain MRC-DPRU4090 (G12P[6]) (Fig 2). On phylogenetic analysis, strain A23 was found to be closely related with strain MRC-DPRU4090 in a common branch with German human strain GER172-08 (G12P[6]) and Bangladeshi human strain MMC71 (G1P[8]) [60] (Fig 7).
See legend of Fig 3. Scale bars, 0.05 substitutions per nucleotide.
https://doi.org/10.1371/journal.pone.0124965.g007
The VP3 genes of strains A14, A25, P02, P39, and P43 showed the maximum nucleotide sequence identities (99.7–99.8%) with those of Thai human strains CU331-NR (G12P[6]) and/or CU460-KK (G12P[8]) (Fig 2). On phylogenetic analysis, strains A14, A25, P02, P39, and P43 were found to be clustered with these strains, and several human G1 and G12 strains (Fig 8). On the other hand, the VP3 gene of strain A23 exhibited the highest nucleotide sequence similarity (99.5%) with Canadian human strain RT103-07 (G1P[8]) (Fig 2). Phylogenetically, strain A23 was found to be closely related with South African human strain MRC-DPRU4090 (G12P[6]) in a common branch with strain RT103-07 and several human G12 strains from Africa (Fig 8).
See legend of Fig 3. Scale bars, 0.05 substitutions per nucleotide.
https://doi.org/10.1371/journal.pone.0124965.g008
The NSP1 genes of strains A14, A25, P02, P39, and P43 exhibited the highest nucleotide sequence identities (99.5–99.6%) with the cognate gene of Australian human strain CK00089 (G1P[8]) (Fig 2). On phylogenetic analysis, strains A14, A25, P02, P39, and P43 were found to be closely related with strain CK00089 in a common branch with Thai human strains CU331-NR (G12P[6]) and CU460-KK (G12P[8]) (Fig 9). On the other hand, the NSP1 gene of strain A23 showed the maximum nucleotide sequence similarity (99.0%) with Argentinian human strains (Arg6627 (G12P[8]) and Arg7500 (G12P[8]) [61]), Paraguayan human strains (1670SR (G1P[8]), 1679SR (G1P[8]), and 1692SR (G1P[8])), and South African human strain MRC-DPRU799 (G1P[8]) (Fig 2). Phylogenetically, strain A23 was found to be clustered near these strains, and several human G1 and G12 strains (Fig 9).
See legend of Fig 3. Scale bars, 0.05 substitutions per nucleotide.
https://doi.org/10.1371/journal.pone.0124965.g009
The NSP2 genes of strains A14, A25, P02, P39, and P43 exhibited the maximum nucleotide sequence identities (99.6–99.8%) with those of Thai human strains CU331-NR (G12P[6]) and CU460-KK (G12P[8]) (Fig 2). On phylogenetic analysis, strains A14, A25, P02, P39, and P43 were found to be closely related with these strains (Fig 10). On the other hand, the NSP2 gene of strain A23 showed the highest nucleotide sequence identity (99.8%) with Ugandan human strain MRC-DPRU4616 (G12P[6]), and comparable identity (99.7%) with South African human strain MRC-DPRU4090 (G12P[6]) (Fig 2). Phylogenetically, strain A23 was found to be closely related with these strains in a common branch with German human strain GER172-08 (G12P[6]) and Russian human strain Nov10-N921 (G9P[6]) [62] (Fig 10).
See legend of Fig 3. Scale bars, 0.02 substitutions per nucleotide.
https://doi.org/10.1371/journal.pone.0124965.g010
The NSP3 genes of strains A14, A25, P02, P39, and P43 showed the highest nucleotide sequence identities (99.3–99.5%) with the NSP3 gene of Thai human strain CU331-NR (G12P[6]) (Fig 2). On phylogenetic analysis, strains A14, A25, P02, P39, and P43 were found to be closely related with strain CU331-NR in a common branch with several human G1, G3, G9, and G12 strains (Fig 11). In contrast, the NSP3 gene of strain A23 exhibited the maximum nucleotide sequence identity (99.2%) with German human strain GER172-08 (G12P[6]), and comparable identity (99.1%) with South African human strains MRC-DPRU4090 (G12P[6]), MRC-DPRU799 (G1P[8]), and MRC-DPRU840 (G1P[8]) (Fig 2). On phylogenetic analysis, strain A23 was found to be closely related with strain MRC-DPRU4090 in a common branch with these, and several human G1 and G12 strains (Fig 11).
See legend of Fig 3. Scale bars, 0.05 substitutions per nucleotide.
https://doi.org/10.1371/journal.pone.0124965.g011
The NSP4 genes of strains A14, A25, P02, P39, and P43 showed the maximum nucleotide sequence similarities (99.4–99.6%) with the cognate genes of Thai human strain CU769-KK (G1P[8]), Chinese human strain LL11131077 (G9P[8]), and/or Russian human strain Nov09-B28 (G1P[8]) (Fig 2). On phylogenetic analysis, strains A14, A25, P02, P39, and P43 were found to be clustered near these, and several human G1, G3, G9, and G12 strains (Fig 12). In contrast, the NSP4 gene of strain A23 exhibited the highest nucleotide sequence identity (99.5%) with Indian human-like bovine strain B85-MP (G3P[3]), South African human strain MRC-DPRU799 (G1P[8]), and American human strain USA2007719739 (G1P[8]) (Fig 2). Phylogenetically, A23 was found to be clustered with these, and several human-like bovine and human strains within the human-like E1 subcluster (Fig 12).
See legend of Fig 3. Scale bars, 0.02 substitutions per nucleotide.
https://doi.org/10.1371/journal.pone.0124965.g012
The NSP5 genes of strains A14, A25, P02, P39, and P43 exhibited the maximum nucleotide sequence identities (99.6–99.8%) with that of Australian human strain CK00089 (G1P[8]) (Fig 2). On phylogenetic analysis, strains A14, A25, P02, P39, and P43 were found to be clustered with strain CK00089, and several human G1 and G12 strains (Fig 13). On the other hand, the NSP5 gene of strain A23 showed the maximum nucleotide sequence similarity (99.6%) with South African human strains MRC-DPRU4090 (G12P[6]) and MRC-DPRU1191 (G12P[8]) (Fig 2). On phylogenetic analysis, strain A23 was found to be clustered near these, and several human G9 and G12 strains (Fig 13).
See legend of Fig 3. Scale bars, 0.02 substitutions per nucleotide.
https://doi.org/10.1371/journal.pone.0124965.g013
In the present study, we analyzed the whole genomes of six Asian G12 strains that have emerged in Myanmar (strains A14, A23, A25, P02, P39, and P43) identified in six stool specimens from children with gastroenteritis. All six Myanmarese G12 strains were found to have a Wa-like backbone: G12-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1 (strains A14, A25, P02, P39, and P43) and G12-P[6]-I1-R1-C1-M1-A1-N1-T1-E1-H1 (strain A23). Whole genomic analyses of currently circulating G12 strains from different parts of the world have revealed the predominant possession of a Wa-like genetic backbone by them [10]. Because a Wa-like genotype backbone is believed to be one of the stable genetic constellations of human RVAs, the Wa-like backbone might have conferred an evolutionary advantage on G12 strains, making it possible for them to spread across the globe in an incredibly short time period [48, 50, 63]. Phylogenetic analysis revealed that in most genes the six strains examined in this study were genetically related to globally circulating human G1, G3, G9, and G12 strains. Furthermore, strains A14, A25, P02, P39, and P43 were very closely related to one another as to all the 11 genes, indicating the derivation of the five strains from a common origin. On the other hand, strain A23 consistently formed distinct clusters for all the 11 gene segments, indicating its distinct origin from that of strains A14, A25, P02, P39, and P43. While eight of the 11 gene segments (VP7, VP4, VP2-3, and NSP1-4) of strains A14, A25, P02, P39, and P43 were very closely related to those of Asian strains, all the 11 gene segments of strain A23 were very closely associated with African strains. Thus, the G12 strains that have emerged in Myanmar appeared to have originated from not a single, but at least two distinct ancestral G12 strains.
Notably, the NSP4 gene of strain A23 having E1 genotype showed the closest relationship with the NSP4 genes of human-like bovine strains from India, as well as human strains (Fig 12). Genotype E1 has been recently detected in bovines in India and Korea, in addition to more common bovine genotype E2 [64, 65]. However, the exact origin of the NSP4 gene of strain A23 could not be ascertained due to a lack of a sufficient number of representative strains as references. In any case, this observation suggests the occurrence of reassortment between human and bovine strains. To our knowledge, this is the first description on full genome-based characterization of the G12 strains that have emerged in Myanmar.
The increasing detection of G12 strains in Asia implies the continued spread of G12 strains in Asia. Because it has not been examined as to whether or not the two available RVA vaccines (Rotarix (GlaxoSmithKline) and RotaTeq (Merck)) are effective for prevention against unusual strains, such as G12P[6] strains that share no G/P genotype specificities with strains included in these vaccines, continuing RVA surveillance of the G12 genotype may be required. Simultaneous monitoring of RVA strains in humans and animals is also essential for a better understanding of RVA ecology. Furthermore, whole genome-based analyses are essential to understand the evolutionary dynamics of G12 strains spreading in Asia.
Author Contributions
Conceived and designed the experiments: TI SK KT. Performed the experiments: TI SK KHM MW KT. Analyzed the data: TI SK KT. Contributed reagents/materials/analysis tools: KWH YYM TWM KZT HMT MMW HNO TH FED KH YF KK SR SVN KU KO TT. Wrote the paper: TI SK KT.
References
- 1. Tate JE, Burton AH, Boschi-Pinto C, Steele AD, Dugue J, Parashar UD, et al. 2008 estimate of worldwide rotavirus-associated mortality in children younger than 5 years before the introduction of universal rotavirus vaccination programmes: a systemic review and meta-analysis. Lancet Infect Dis. 2012;12: 136–141. pmid:22030330
- 2. Kahn G, Fitzwater S, Tate J, Kang G, Ganguly N, Nair G, et al. Epidemiology and prospects for prevention of rotavirus disease in India. Indian Pediatr. 2012;49: 467–474. pmid:22796685
- 3. Moe K, Hummelman EG, Oo WM, Lwin T, Htwe TT. Hospital-based surveillance for rotavirus diarrhea in children in Yangon, Myanmar. J Infect Dis. 2005;192 Suppl 1: S111–S113. pmid:16088793
- 4. Moe K, Thu HM, Oo WM, Aye KM, Shwe TT, Mar W, et al. Genotyping of rotavirus isolates collected from children less than 5 years of age admitted for diarrhoea at the Yangon Children’s Hospital, Myanmar. Vaccine. 2009;27 Suppl 5: F89–F92. pmid:19931728
- 5.
Estes MK, Greenberg HB. Rotaviruses. In: Knipe DM, Howley PM, editors. Fields Virology. 6th ed. Philadelphia: Lippincott Williams & Wilkins; 2013. pp. 1347–1401.
- 6. Matthijnssens J, Ciarlet M, McDonald SM, Attoui H, Bányai K, Brister JR, et al. Uniformity of rotavirus strain nomenclature proposed by the Rotavirus Classification Working Group (RCWG). Arch Virol. 2011;156: 1397–1413. pmid:21597953
- 7. Trojnar E, Sachsenröder J, Twardziok S, Reetz J, Otto PH, Johne R. Identification of an avian group A rotavirus containing a novel VP4 gene with a close relationship to those of mammalian rotaviruses. J Gen Virol. 2013;94: 136–142. pmid:23052396
- 8. Santos N, Hoshino Y. Global distribution of rotavirus serotypes/genotypes and its implication for the development and implementation of an effective rotavirus vaccine. Rev Med Virol. 2005;15: 29–56. pmid:15484186
- 9. Matthijnssens J, Heylen E, Zeller M, Rahman M, Lemey P, Van Ranst M. Phylodynamic analyses of rotavirus genotypes G9 and G12 underscore their potential for swift global spread. Mol Biol Evol. 2010;27: 2431–2436. pmid:20522727
- 10. Rahman M, Matthijnssens J, Yang X, Delbeke T, Arijs I, Taniguchi K, et al. Evolutionary history and global spread of the emerging G12 human rotaviruses. J Virol. 2007;81: 2382–2390. pmid:17166908
- 11. Matthijnssens J, Bilcke J, Ciarlet M, Martella V, Bányai K, Rahman M, et al. Rotavirus disease and vaccination: impact on genotype diversity. Future Microbiol. 2009; 4: 1303–1316. pmid:19995190
- 12. Taniguchi K, Urasawa T, Kobayashi N, Gorziglia M, Urasawa S. Nucleotide sequence of VP4 and VP7 genes of human rotaviruses with subgroup I specificity and long RNA pattern: implication for new G serotype specificity. J Virol. 1990;64: 5640–5644. pmid:2170690
- 13. Urasawa S, Urasawa T, Wakasugi F, Kobayashi N, Taniguchi K, Lintag IC, et al. Presumptive seventh serotype of human rotavirus. Arch Virol. 1990;113: 279–282. pmid:2171461
- 14. Pongsuwanna Y, Guntapong R, Chiwakul M, Tacharoenmuang R, Onvimala N, Wakuda M, et al. Detection of a human rotavirus with G12 and P[9] specificity in Thailand. J Clin Microbiol. 2002;40: 1390–1394. pmid:11923362
- 15. Wakuda M, Nagashima S, Kobayashi N, Pongsuwanna Y, Taniguchi K. Serologic and genomic characterization of a G12 human rotavirus in Thailand. J Clin Microbiol. 2003;41: 5764–5769. pmid:14662978
- 16. Shinozaki K, Okada M, Nagashima S, Kaiho I, Taniguchi K. Characterization of human rotavirus strains with G12 and P[9] detected in Japan. J Med Virol. 2004;73: 612–616. pmid:15221908
- 17. Kawai K, O’Brien MA, Goveia MG, Mast TC, El Khoury AC. Burden of rotavirus gastroenteritis and distribution of rotavirus strains in Asia: a systemic review. Vaccine. 2012;30: 1211–1254.
- 18. Paul SK, Kobayashi N, Nagashima S, Ishino M, Watanabe S, Alam MM, et al. Phylogenetic analysis of rotaviruses with genotypes G1, G2, G9 and G12 in Bangladesh: evidence for a close relationship between rotaviruses from children and adults. Arch Virol. 2008;153: 1999–2012. pmid:18839059
- 19. Paul SK, Hossain MA, Ahmed MU, Alam MM, Musa AK, Shamsuzzaman AK, et al. Prevalence of VP7 and VP4 genotypes of human group A rotavirus in infants and children with acute diarrhea in a northern city of Bangladesh. Mymensingh Med J. 2009;18: 190–197. pmid:19623146
- 20. Paul SK, Ahmed MU, Hossain MA, Mahmud MC, Bhuiyan MR, Saha SK, et al. Molecular characterization of group A human rotavirus among hospitalized children and adults in Bangladesh: finding of emerging G12 strain. Mymensingh Med J. 2010;19: 16–26. pmid:20046166
- 21. Ahmed K, Ahmed S, Mitui MT, Rahman A, Kabir L, Hannan A, et al. Molecular characterization of VP7 gene of human rotaviruses from Bangladesh. Virus Genes. 2010;40: 347–356. pmid:20217207
- 22. Ahmed S, Hussain M, Akhter S, Islam T, Ahmed SU, Kabir ML. Genotypes of rotavirus diarrhoea in a children hospital of Bangladesh. Mymensingh Med J. 2012;21: 497–502. pmid:22828550
- 23. Wangchuk S, Mitui MT, Tshering K, Yahiro T, Bandhari P, Zangmo S, et al. Dominance of emerging G9 and G12 genotypes and polymorphism of VP7 and VP4 of rotaviruses from Bhutanese children with severe diarrhea prior to the introduction of vaccine. PLoS One. 2014;9: e110795. pmid:25330070
- 24. Das S, Varghese V, Chaudhury S, Barman P, Mahapatra S, Kojima K, et al. Emergence of novel human group A rotavirus G12 strains in India. J Clin Microbiol. 2003;41: 2760–2762. pmid:12791925
- 25. Samajdar S, Varghese V, Barman P, Ghosh S, Mitra U, Dutta P, et al. Changing pattern of human group A rotaviruses: emergence of G12 as an important pathogen among children in eastern India. J Clin Virol. 2006;36: 183–188. pmid:16679056
- 26. Ramani S, Banerjee I, Gladstone BP, Sarkar R, Selvapandian D, Le Fevre AM, et al. Geographic information systems and genotyping in identification of rotavirus G12 infections in residents of an urban slum with subsequent detection in hospitalized children: emergence of G12 genotype in South India. J Clin Microbiol. 2007;45: 432–437. pmid:17135437
- 27. Ray P, Sharma S, Agarwal RK, Longmei K, Gentsch JR, Paul VK, et al. First detection of G12 rotaviruses in newborns with neonatal rotavirus infection at all India Institute of Medical Sciences, New Delhi, India. J Clin Microbiol. 2007;45: 3824–3827. pmid:17728476
- 28. Sharma S, Ray P, Gentsch JR, Glass RI, Kalra V, Bhan MK. Emergence of G12 rotavirus strains in Delhi, India, in 2000 to 2007. J Clin Microbiol. 2008;46: 1343–1348. pmid:18272705
- 29. Kang G, Arora R, Chitambar SD, Deshpande J, Gupte MD, Kulkarni M, et al. Multicenter, hospital-based surveillance of rotavirus disease and strains among Indian children aged <5 years. J Infect Dis. 2009;200 Suppl 1: S147–S153. pmid:19817593
- 30. Chakravarti A, Chauhan MS, Sharma A, Verma V. Distribution of human rotavirus G and P genotypes in a hospital setting from Northern India. Southeast Asian J Trop Med Public Health. 2010;41: 1145–1152. pmid:21073035
- 31. Mukherjee A, Chattopadhyay S, Bagchi P, Dutta D, Singh NB, Arora R, et al. Surveillance and molecular characterization of rotavirus strains circulating in Manipur, north-eastern India: increasing prevalence of emerging G12 strains. Infect Genet Evol. 2010;10: 311–320. pmid:20083233
- 32. Kang G, Desai R, Arora R, Chitamabar S, Naik TN, Krishnan T, et al. Diversity of circulating rotavirus strains in children hospitalized with diarrhea in India, 2005–2009. Vaccine. 2013;31: 2879–2883. pmid:23624096
- 33. Mullick S, Mukherjee A, Ghosh S, Pazhani GP, Sur D, Manna B, et al. Community based case-control study of rotavirus gastroenteritis among young children during 2008–2010 reveals vast genetic diversity and increased prevalence of G9 strains in Kolkata. PLoS One. 2014;9: e112970. pmid:25401757
- 34. Wulan WN, Listiyaningsih E, Samsi KM, Agtini MD, Kasper MR, Putnam SD. Identification of a rotavirus G12 strain, Indonesia. Emerg Infect Dis. 2010;16: 159–161. pmid:20031072
- 35. Vainio K, Latipov R, Utegenova E, Kasymbekova K, Juraev R, Asilova M, et al. Rotavirus genotype distribution in Kyrgyzstan and Kazakhstan, 2007–2009. APMIS. 2013;121: 447–455. pmid:23078218
- 36. Le VP, Kim JY, Cho SL, Nam SW, Lim I, Lee HJ, et al. Detection of unusual rotavirus genotypes G8P[8] and G12P[6] in South Korea. J Med Virol. 2008;80: 175–182. pmid:18041003
- 37. Uchida R, Pandey BD, Sherchand JB, Ahmed K, Yokoo M, Nakagomi T, et al. Molecular epidemiology of rotavirus diarrhea among children and adults in Nepal: detection of G12 strains with P[6] or P[8] and a G11P[25] strain. J Clin Microbiol. 2006;44: 3499–3505. pmid:17021073
- 38. Pun SB, Nakagomi T, Sherchand JB, Pandey BD, Cuevas LE, Cunliffe NA, et al. Detection of G12 human rotaviruses in Nepal. Emerg Infect Dis. 2007;13: 482–484. pmid:17552107
- 39. Sherchand JB, Nakagomi O, Dove W, Nakagomi T, Yokoo M, Pandey BD, et al. Molecular epidemiology of rotavirus diarrhea among children aged <5 years in Nepal: predominance of emergent G12 strains during 2 years. J Infect Dis. 2009;200 Suppl 1: S182–S187. pmid:19817599
- 40. Ansari S, Sherchand JB, Rijal BP, Parajuli K, Mishra SK, Dahal RK, et al. Characterization of rotavirus causing acute diarrhoea in children in Kathmandu, Nepal, showing the dominance of serotype G12. J Med Microbiol. 2013;62: 114–120. pmid:23038804
- 41. Alam MM, Malik SA, Shaukat S, Naeem A, Sharif S, Angez M, et al. Genetic characterization of rotavirus subtypes in Pakistan-first report of G12 genotype from Pakistan under WHO-Eastern Mediterranean region. Virus Res. 2009;144: 280–284. pmid:19720243
- 42. Tamim S, Hasan F, Matthijnssens J, Sharif S, Shaukat S, Alam MM, et al. Epidemiology and phylogenetic analysis of VP7 and VP4 genes of rotaviruses circulating in Rawalpindi, Pakistan during 2010. Infect Genet Evol. 2013;14: 161–168. pmid:23187023
- 43. Ahmed K, Batuwanthudawe R, Chandrasena TG, Mitui MT, Rajindrajith S, Galagoda G, et al. Rotavirus infections with multiple emerging genotypes in Sri Lanka. Arch Virol. 2010;155: 71–75. pmid:19921392
- 44. Pongsuwanna Y, Guntapong R, Tacharoenmuang R, Prapanpoj M, Kameoka M, Taniguchi K. A long-term survey on the distribution of the human rotavirus G type in Thailand. J Med Virol. 2010;82: 157–163. pmid:19950228
- 45. Khananurak K, Vutithanachot V, Simakachorn N, Theamboonlers A, Chongsrisawat V, Poovorawan Y. Prevalence and phylogenetic analysis of rotavirus genotypes in Thailand between 2007 and 2009. Infect Genet Evol. 2010;10: 537–545. pmid:20156600
- 46. Theamboonlers A, Maiklang O, Thongmee T, Chieochansin T, Vuthitanachot V, Poovorawan Y. Complete genotype constellation of human rotavirus group A circulating in Thailand. Infect Genet Evol. 2014;21: 295–302. pmid:24296010
- 47. Tra My PV, Rabaa MA, Vinh H, Holms EC, Hoang NV, Vinh NT, et al. The emergence of rotavirus G12 and the prevalence of enteric viruses in hospitalized pediatric diarrheal patients in southern Vietnam. Am J Trop Med Hyg. 2011;85: 768–775. pmid:21976585
- 48. Matthijnssens J, Ciarlet M, Heiman E, Arijs I, Delbeke T, McDonald SM, et al. 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. J Virol. 2008;82: 3204–3219. pmid:18216098
- 49. Dennis FE, Fujii Y, Haga K, Damanka S, Lartey B, Agbemabiese CA, et al. Identification of novel Ghanaian G8P[6] human-bovine reassortant rotavirus strain by Next Generation Sequencing. PLoS One. 2014;9: e100699. pmid:24971993
- 50. Ghosh S, Kobayashi N. Whole-genomic analysis of rotavirus strains: current status and future prospects. Future Microbiol. 2011;6: 1049–1065. pmid:21958144
- 51. Komoto S, Maeno Y, Tomita M, Matsuoka T, Ohfu M, Yodoshi T, et al. Whole genomic analysis of a porcine-like human G5P[6] rotavirus strain isolated from a child with diarrhoea and encephalopathy in Japan. J Gen Virol. 2013;94: 1568–1575. pmid:23515025
- 52. Komoto S, Sasaki J, Taniguchi K. Reverse genetics system for introduction of site-specific mutations into the double-stranded RNA genome of infectious rotavirus. Proc Natl Sci USA. 2006;103: 4646–4651. pmid:16537420
- 53. Komoto S, Pongsuwanna Y, Ide T, Wakuda M, Guntapong R, Dennis FE, et al. Whole genomic analysis of porcine G10P[5] rotavirus strain P343 provides evidence for bovine-to-porcine interspecies transmission. Vet Microbiol. 2014;174: 577–583. pmid:25457370
- 54. Komoto S, Apondi EW, Shah M, Odoyo E, Nyangao J, Tomita M, et al. Whole genomic analysis of human G12P[6] and G12P[8] rotavirus strains that have emerged in Kenya: Identification of porcine-like NSP4 genes. Infect Genet Evol. 2014;27: 277–293. pmid:25111611
- 55. Maes P, Matthijnssens J, Rahman M, Van Ranst M. RotaC: a web-based tool for the complete genome classification of group A rotaviruses. BMC Microbiol. 2009;9: 238. pmid:19930627
- 56. Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22: 4673–4680. pmid:7984417
- 57. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30: 2725–2729. pmid:24132122
- 58. Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985;39: 783–791.
- 59. Pietsch C, Liebert UG. Human infection with G12 rotaviruses, Germany. Emerg Infect Dis. 2009;15: 1512–1515. pmid:19788829
- 60. Ghosh S, Paul SK, Yamamoto D, Nagashima S, Kobayashi N. Full genomic analyses of human rotavirus strains possessing the rare P[8]b VP4 subtype. Infect Genet Evol. 2011;11: 1481–1486. pmid:21640848
- 61. Stupka JA, Degiuseppe JI, Parra GI; Argentinean National Rotavirus Surveillance Network. Increased frequency of rotavirus G3P[8] and G12P[8] in Argentina during 2008–2009: whole-genome characterization of emerging G12P[8] strains. J Clin Virol. 2012;54: 162–167. pmid:22410133
- 62. Zhirakovskaia EV, Aksanova RKh, Gorbunova MG, Tikunov Alu, Kuril’shchikov AM, Sokolov SN, et al. Genetic diversity of group A rotavirus isolates in Western Siberia in 2007–2011. Mol Gen Mikrobiol Virusol. 2012;4: 33–41. pmid:23248851
- 63. Matthijnssens J, Van Ranst M. Genotype constellation and evolution of group A rotaviruses infecting humans. Curr Opin Virol. 2012;2: 426–433. pmid:22683209
- 64. Malik YS, Kumar N, Sharma K, Ghosh S, Bányai K, Balasubramanian G, et al. Molecular analysis of non structural rotavirus group A enterotoxin gene of bovine origin from India. Infect Genet Evol. 2014;25: 20–27. pmid:24747605
- 65. Park JG, Kim DS, Matthijnssens J, Kwon HJ, Zeller M, Alfajaro MM, et al. Comparison of pathogenicities and nucleotide changes between porcine and bovine reassortant rotavirus strains possessing the same genotype constellation in piglets and calves. Vet Microbiol. 2014;172: 51–62. pmid:24861840