Skip to main content
Advertisement
Browse Subject Areas
?

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

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Whole Genomic Analysis of Human G12P[6] and G12P[8] Rotavirus Strains that Have Emerged in Myanmar

  • Tomihiko Ide,

    Affiliation Department of Virology and Parasitology, Fujita Health University School of Medicine, Toyoake, Aichi, Japan

  • Satoshi Komoto ,

    satoshik@fujita-hu.ac.jp

    Affiliation Department of Virology and Parasitology, Fujita Health University School of Medicine, Toyoake, Aichi, Japan

  • Kyoko Higo-Moriguchi,

    Affiliation Department of Virology and Parasitology, Fujita Health University School of Medicine, Toyoake, Aichi, Japan

  • Khaing Win Htun,

    Affiliation Nay Pyi Taw General Hospital (Central Myanmar), Nay Pyi Taw, Myanmar

  • Yi Yi Myint,

    Affiliation Department of Medical Research (Upper Myanmar), Pyin Oo Lwin, Myanmar

  • Theingi Win Myat,

    Affiliation Department of Medical Research (Lower Myanmar), Yangon, Myanmar

  • Kyaw Zin Thant,

    Affiliation Department of Medical Research (Lower Myanmar), Yangon, Myanmar

  • Hlaing Myat Thu,

    Affiliation Department of Medical Research (Lower Myanmar), Yangon, Myanmar

  • Mo Mo Win,

    Affiliation Department of Medical Research (Lower Myanmar), Yangon, Myanmar

  • Htun Naing Oo,

    Affiliation Department of Traditional Medicine (Central Myanmar), Nay Pyi Taw, Myanmar

  • Than Htut,

    Affiliation Ministry of Health (Central Myanmar), Nay Pyi Taw, Myanmar

  • Mitsutaka Wakuda,

    Affiliation Department of Virology and Parasitology, Fujita Health University School of Medicine, Toyoake, Aichi, Japan

  • Francis Ekow Dennis,

    Affiliations Department of Virology II, National Institute of Infectious Diseases, Musashi-Murayama, Tokyo, Japan, Department of Electron Microscopy and Histopathology, Noguchi Memorial Institute for Medical Research, College of Health Sciences, University of Ghana, Legon, Ghana

  • Kei Haga,

    Affiliation Department of Virology II, National Institute of Infectious Diseases, Musashi-Murayama, Tokyo, Japan

  • Yoshiki Fujii,

    Affiliation Department of Virology II, National Institute of Infectious Diseases, Musashi-Murayama, Tokyo, Japan

  • Kazuhiko Katayama,

    Affiliation Department of Virology II, National Institute of Infectious Diseases, Musashi-Murayama, Tokyo, Japan

  • Shofiqur Rahman,

    Affiliation Immunology Research Institute in Gifu, EW Nutrition Japan, Gifu, Japan

  • Sa Van Nguyen,

    Affiliation Immunology Research Institute in Gifu, EW Nutrition Japan, Gifu, Japan

  • Kouji Umeda,

    Affiliation Immunology Research Institute in Gifu, EW Nutrition Japan, Gifu, Japan

  • Keiji Oguma,

    Affiliation Department of Bacteriology, Okayama University Graduate School of Medicine, Dentistry, and Pharmaceutical Sciences, Okayama, Japan

  • Takao Tsuji,

    Affiliation Department of Microbiology, Fujita Health University School of Medicine, Toyoake, Aichi, Japan

  •  [ ... ],
  • Koki Taniguchi

    Affiliation Department of Virology and Parasitology, Fujita Health University School of Medicine, Toyoake, Aichi, Japan

  • [ view all ]
  • [ view less ]

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.

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 [911].

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 [911, 1417]. In Asia, G12 strains are being increasingly identified in diarrheic children in several countries; Bangladesh [10, 1822], Bhutan [23], India [2433], Indonesia [34], Japan [16], Kazakhstan [35], Korea [36], Kyrgyzstan [35], Nepal [3740], Pakistan [41, 42], Sri Lanka [43], Thailand [14, 15, 4446], 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.

thumbnail
Fig 1. Genomic dsRNA profiles of strains A14, A23, A25, P02, P39, and P43.

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).

thumbnail
Fig 2. Genotype natures of the 11 gene segments of six Myanmarese G12 strains, A14, A23, A25, P02, P39, and P43, with those of selected human and animal strains.

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).

thumbnail
Fig 3. Phylogenetic tree constructed from the nucleotide sequences of the VP7 genes of strains A14, A23, A25, P02, P39, and P43, and representative RVA strains.

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).

thumbnail
Fig 4. Phylogenetic tree constructed from the nucleotide sequences of the VP4 genes of strains A14, A23, A25, P02, P39, and P43, and representative RVA strains.

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).

thumbnail
Fig 5. Phylogenetic tree constructed from the nucleotide sequences of the VP6 genes of strains A14, A23, A25, P02, P39, and P43, and representative RVA strains.

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).

thumbnail
Fig 6. Phylogenetic tree constructed from the nucleotide sequences of the VP1 genes of strains A14, A23, A25, P02, P39, and P43, and representative RVA strains.

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).

thumbnail
Fig 7. Phylogenetic tree constructed from the nucleotide sequences of the VP2 genes of strains A14, A23, A25, P02, P39, and P43, and representative RVA strains.

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).

thumbnail
Fig 8. Phylogenetic tree constructed from the nucleotide sequences of the VP3 genes of strains A14, A23, A25, P02, P39, and P43, and representative RVA strains.

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).

thumbnail
Fig 9. Phylogenetic tree constructed from the nucleotide sequences of the NSP1 genes of strains A14, A23, A25, P02, P39, and P43, and representative RVA strains.

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).

thumbnail
Fig 10. Phylogenetic tree constructed from the nucleotide sequences of the NSP2 genes of strains A14, A23, A25, P02, P39, and P43, and representative RVA strains.

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).

thumbnail
Fig 11. Phylogenetic tree constructed from the nucleotide sequences of the NSP3 genes of strains A14, A23, A25, P02, P39, and P43, and representative RVA strains.

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).

thumbnail
Fig 12. Phylogenetic tree constructed from the nucleotide sequences of the NSP4 genes of strains A14, A23, A25, P02, P39, and P43, and representative RVA strains.

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).

thumbnail
Fig 13. Phylogenetic tree constructed from the nucleotide sequences of the NSP5 genes of strains A14, A23, A25, P02, P39, and P43, and representative RVA strains.

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.

Acknowledgments

We wish to thank Mayuko Tomita for her technical assistance.

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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 58. Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985;39: 783–791.
  59. 59. Pietsch C, Liebert UG. Human infection with G12 rotaviruses, Germany. Emerg Infect Dis. 2009;15: 1512–1515. pmid:19788829
  60. 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. 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. 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. 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. 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. 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