Molecular Epidemiology and Clinical Features of Rotavirus Infection Among Pediatric Patients in East Java, Indonesia During 2015–2018: Dynamic Changes in Rotavirus Genotypes From Equine-Like G3 to Typical Human G1/G3

Group A rotavirus (RVA) is the most important cause of severe gastroenteritis among children worldwide, and effective RVA vaccines have been introduced in many countries. Here we performed a molecular epidemiological analysis of RVA infection among pediatric patients in East Java, Indonesia, during 2015–2018. A total of 432 stool samples were collected from hospitalized pediatric patients with acute gastroenteritis. None of the patients in this cohort had been immunized with an RVA vaccine. The overall prevalence of RVA infection was 31.7% (137/432), and RVA infection was significantly more prevalent in the 6- to 11-month age group than in the other age groups (P < 0.05). Multiplex reverse transcription-PCR (RT-PCR) revealed that the most common G-P combination was equine-like G3P[8] (70.8%), followed by equine-like G3P[6] (12.4%), human G1P[8] (8.8%), G3P[6] (1.5%), and G1P[6] (0.7%). Interestingly, the equine-like strains were exclusively detected until May 2017, but in July 2017 they were completely replaced by a typical human genotype (G1 and G3), suggesting that the dynamic changes in RVA genotypes from equine-like G3 to typical human G1/G3 in Indonesia can occur even in the country with low RVA vaccine coverage rate. The mechanism of the dynamic changes in RVA genotypes needs to be explored. Infants and children with RVA-associated gastroenteritis presented more frequently with some dehydration, vomiting, and watery diarrhea, indicating a greater severity of RVA infection compared to those with non-RVA gastroenteritis. In conclusion, a dynamic change was found in the RVA genotype from equine-like G3 to a typical human genotype. Since severe cases of RVA infection were prevalent, especially in children aged 6 to 11 months or more generally in those less than 2 years old, RVA vaccination should be included in Indonesia’s national immunization program.

Introduction

Acute gastroenteritis is a common and important cause of childhood morbidity and mortality in developing countries (Elliott, 2007). Diarrhea, the most common symptom of acute gastroenteritis, remains the third leading cause of overall morbidity and the leading cause of infant mortality in Indonesia, and 38–67% of children with diarrhea admitted to hospitals in Indonesia were infected with group A rotavirus (RVA) (Soenarto et al., 1981; Sudarmo et al., 2015; Nirwati et al., 2016). RVA, which is known as ‘winter diarrhea,’ exhibits distinct seasonality in temperate areas, whereas the effect of season on RVA infection is not extreme in the tropics (Levy et al., 2009).

RVA belongs to the Reoviridae family and is a non-enveloped virus, consisting of 11 segments of double-stranded RNA and encoding six structural and six non-structural proteins (NSPs). RVA is classified into 35 G and 50 P types on the basis of VP7 (glycoprotein) and VP4 (protease-sensitive) protein (Rojas et al., 2017), respectively. RVA genotypes G1P[8], G2P[4], G3P[8], G4P[8], G9P[8], and G12P[8] are commonly found worldwide (Santos and Hoshino, 2005; Dóró et al., 2014). These genotypes were also commonly detected in Indonesia during 2010–2015 (Sudarmo et al., 2015; Mulyani et al., 2018). There are also three genotype constellations of human RVA, termed Wa-like (G1-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1), DS-1-like (G2-P[4]- I2- R2-C2-M2-A2-N2-T2-E2-H2) and the less common AU-1-like (G3-P[9]- I3- R3-C3-M3-A3-N3-T3-E3-H3) (Matthijnssens et al., 2008). The G1, G3, G4, G9 and G12 are usually grouped into Wa-like constellation, and G2 and G8 are grouped into DS-1-like constellation. We recently reported equine-like G3P[8]/P[6] strains with a DS-1 like backbone (equine-like G3-P[8]/[6]- I2- R2-C2-M2-A2-N2-T2-E2-H2), resulting from a rare human/equine reassortment, circulating among pediatric inpatients with acute gastroenteritis at a private hospital in the suburbs of Surabaya, Indonesia (Utsumi et al., 2018). All the tested samples were determined to be the equine-like G3 genotype, which was initially isolated in Australia and Japan, later emerging globally (Malasao et al., 2015; Arana et al., 2016; Cowley et al., 2016; Doro et al., 2016; Guerra et al., 2016; Komoto et al., 2018; Pietsch and Liebert, 2018). A report from Australia pointed out that equine-like G3P[8] was more common in the areas where Rotarix^® was being used (Roczo-Farkas et al., 2016), and another study found that the equine-like G3P[8] strain was responsible for the most episodes of severe gastroenteritis (Cowley et al., 2016).

Two kinds of RVA vaccines (Rotarix^® and RotaTeq^®) were licensed in 2006, and 84 countries have introduced Rotarix^® or RotaTeq^® into their national immunization programs (Abou-Nader et al., 2018). The introduction of RVA vaccines into human populations may impose a selective pressure on circulating RVA strains especially in the post-vaccine era (Matthijnssens et al., 2010), although the current RVA vaccines are highly effective. Post-vaccination surveillance studies in the countries where the Rotarix^® vaccine was commonly used, showed significant reduction of RVA prevalences and a genotype shift in proportion from G1P[8] to G2P[4] strains (Dóró et al., 2014; Mukhopadhya et al., 2017). Indonesia has not included an RVA vaccine in its national immunization program yet, even though RVA vaccines have been commercially available since 2011. Little is known about the coverage and impact of RVA vaccines in Indonesia.

The purpose of the present study was to evaluate the transition of prevalences and genotypes of RVA and the severity of RVA-associated acute gastroenteritis among children admitted to a referral hospital in East Java, Indonesia.

Materials and Methods

Study Populations and Stool Samples

A total of 432 stool samples were collected from pediatric patients who were admitted to a government-owned hospital with a clinical diagnosis of acute gastroenteritis between September 2015 and March 2018. This hospital is a referral hospital in East Java, located in Surabaya, the second-largest city in Indonesia. The study population ranged in age from 1 month to 15 years and 11 months (median age, 12.5 months). The patients were classified into 7 age groups: <6, 6–11, 12–23, 24–35, 36–47, 48–59, and ≥60 months. The criterion for acute gastroenteritis was defined as the occurrence of ≥3 times watery or looser-than-normal stools per day and lasts for less than 14 days. To assess the shifting of RVA genotypes from 2013 to 2018 in East Java, 4 RVA-positive stool samples collected in 2013 and 3 in 2014 from pediatric patients with acute gastroenteritis at the same hospital were randomly included. None of the patients in this cohort had been immunized with an RVA vaccine. Demographic data (age, sex) and clinical features (type of stool, vomiting, diarrhea, fever, dehydration) were recorded. The severity of dehydration was determined using the WHO criteria (World Health Organization [WHO], 2014). The seasonality regional rainfall data of Surabaya were obtained from the ClimatView organized by Japan Meteorological Agency¹, because rainfall has been considered as the indicator for seasonality in Tropic areas. Written informed consent was obtained from parents or guardians of all the children. The study protocol was reviewed and approved by the ethics committees of the government-owned hospital (No. 188) and Airlangga University (No. 2054) in Indonesia, and of Kobe University (No. 1857) in Japan.

Rotavirus Detection by Immunochromatography Test

All the stool samples were screened for RVA antigen using the Dipstick “Eiken” Rota kit (Eiken Chemical, Co., Tokyo, Japan).

RNA Extraction and RT-PCR Genotyping

Viral RNA was extracted from 140 μl of the supernatant with a QIAamp Viral RNA mini kit (Qiagen, Hilden, Germany). Positive samples determined to be positive by an immunochromatography test were subjected to genotyping in the VP7 (G typing) and VP4 (P typing) genes by multiplex reverse transcription-PCR (RT-PCR) with newly designed primers for equine-like G3 in the VP7 gene (Fujii et al., 2019b), and with the primers for VP4 that have been designed and used in our lab (Table 1). Equine-like G3 in the VP7 was typed with the primer G3e-757F (Table 1).

TABLE 1

Table 1. Primer set for VP7 genotyping.

Nucleotide Sequencing and Phylogenetic Analyses

Twenty-three (4 collected in 2013, 3 collected in 2014, and 16 collected in 2015–2018) RVA-positive specimens detected by RT-PCR contained sufficient RNA for further whole genome characterizations and were analyzed by next-generation sequencing (NGS). The cDNA library preparation and Illumina MiSeq (Illumina, San Diego, CA, United States) sequencing were performed as previously described (Doan et al., 2017). The full-length nucleotide sequence of each genome segment of the RVAs was obtained using CLC Genomics Workbench v7.0.3 software (CLC Bio, Aarhus, Denmark) with the assembled contigs as query sequences and 11 genome segments of a reference RVA as the target sequences. The genotypes of the 11 genome segments of the 23 Indonesian strains were determined using the RotaC v2.0 web-based classification tool². The genome sequences were aligned with reference sequences using CLUSTAL X (version 1.83) software, and the phylogenetic trees were constructed by the neighbor-joining method. To confirm the reliability of the phylogenetic tree analysis, bootstrap resampling and reconstruction were carried out 1000 times. These analyses were carried out using Molecular Evolutionary Genetic Analysis (MEGA4) software (Tamura et al., 2007). The gene sequences described in the present study have been deposited in the GenBank database under accession numbers (LC434517–LC434532, LC434537–LC434541, LC430883–LC430888, LC469324–LC469341 and LC469398–LC469603).

Statistical Analysis

Statistical analysis was performed by the chi-squared test or Fisher’s exact test for categorical variables using SPSS Statistics 17.0 (Advanced Analytics, Inc., Tokyo, Japan). Results were considered statistically significant at P < 0.05.

Results

Prevalence of RVA

RVA was detected from 137 samples by immunochromatography testing, and the results were confirmed by PCR. Thus, 31.7% of the 432 pediatric patients with acute diarrhea were found to be infected with RVA. The prevalences of RVA infection in 2015, 2016, and 2017 were 31.6, 36.3, and 24.2%, respectively. The relationship between seasonal distribution of RVA infection and rainfall is shown in Figure 1. RVA infections were seen throughout the year and peaked in November 2017, followed by June 2016. Approximately 88% (121/137) of the RVA-positive samples were from children less than 2 years old. RVA infection was significantly most prevalent in the 6- to 11-month age group (P < 0.05).

FIGURE 1

Figure 1. Distribution of equine-like G3P[8]/P[6] (black columns) and global G1/G3P[8]/[6] (gray columns) rotavirus strains collected in East Java, Indonesia, between September 2015 and March 2018. Numbers of rotavirus-positive samples are expressed as percentages of total rotavirus-positive samples in each month of the sampling period.

Distribution, Transition of G and P Genotypes of RVA and Complete Genome Constellations of RVA Strains

Of the 137 pediatric patients with RVA-associated gastroenteritis, the G genotype was determined in all the cases, while the P genotype was determined in 130 cases. One hundred thirty-seven samples were categorized into three G genotypes (typical human G1, G3, and equine-like G3), and 130 samples were categorized into two P genotypes (P[6] and P[8]). The equine-like G3 strains were determined as the dominant genotype (87.6%), followed by human G1 (9.5%) and human G3 (2.9%). Among the P genotypes identified here, P[8] was the dominant genotype (80.3%), followed by P[6] (14.6%), and 7 strains (5.1%) were non-typeable, because PCR products were not amplified by RT-PCR using 2^nd PCR primer set for VP4 genotyping (Tables 1, 2). The most common G-P combination was equine-like G3P[8], accounting for 80.0%, followed by equine-like G3P[6] (12.4%), human G1P[8] (8.8%), G3P[6] (1.5%), G3P[8] (1.5%), and G1P[6] (0.7%) (Table 2). The G-P combinations of the 7 additional samples used as controls in this study were classified into equine-like G3P[8] (collected in 2013) and human G1P[8] (collected in 2013 and 2014). Most of the RVA strains detected were clustered with equine-like G3 during 2015–2017, whereas human genotype G1 became predominant during 2017–2018 (Table 3). In this study, equine-like G3 strains were detected until May 2017, but in July 2017 they were completely replaced by typical human genotypes G1 and G3 (Figure 1).

TABLE 2

Table 2. Total genotype distribution of RNA strains detected among diarrheal children in East Java during 2015–2018.

TABLE 3

Table 3. Yearly genotype distribution of RVA strains detected among diarrheal children in East Java during 2015–2016.

The genotype constellation for 11 genome segments of 23 representative strains from 2013 to 2018 was shown in Figure 2. Typical human G1P[8] and equine-like G3P[8] with DS-1-like backbone was exhibited in 2013 and 2014, then epidemic genotype was completely replaced by equine-like G3P[8]/[6] with DS-1-like backbone during 2015–2017 (May). Since July 2017, epidemic genotype has dynamically shifted from equine-like G3P[8]/[6] to typical human G1/G3P[8]/[6] regardless of DS-1-like or Wa-like genotype constellation.

FIGURE 2

Figure 2. Genotype constellation comparison by the study period. The 23 Indonesian RVA strains in this study were analyzed by whole genome sequences. Green and red indicate Wa-like and DS-1 like gene segments, respectively. The P[6] VP4 genotype is colored yellow, brown is used to indicate a gene segment of equine origin, and blue is used to indicate a gene segment of bovine origin.

Nucleotide Sequences and Phylogenetic Analysis of Equine-Like G3 and Other Human RVAs

In the VP7 gene, representative strains were classified into equine-like G3 lineage together with the strains previously reported as equine-like G3 in another location in East Java (Utsumi et al., 2018) and in other countries such as Australia, Thailand, Spain, Japan, and Brazil (Malasao et al., 2015; Arana et al., 2016; Cowley et al., 2016; Doro et al., 2016; Guerra et al., 2016; Komoto et al., 2016) (Figure 3). Two control strains collected in 2013 were also classified into equine-like G3 lineage (Figure 3). The rest of the representative strains were classified into human G1 and G3 genotypes together with the strains collected in 2013 and 2014 (Figure 3). The nucleotide sequence identities of the VP7 gene between the equine-like G3 strains in this study were extremely high (99.0–100%) regardless of [P] genotype. Additionally, its identities between equine-like G3 strains in this study and those in other countries were high (98.7–99.6%), suggesting that relationship of equine-like G3P[8] VP7 gene and equine-like G3P[6] VP7 gene were genetically close. In the VP4 gene, the representative strains in this study were classified into P[8] and P[6] genotypes (Figure 4). In the VP4 gene, equine-like G3 and human G1/G3 strains tended to have their own group in P[8] genotype with high nucleotide identities (≥99.1%), respectively. The nucleotide identities between equine-like G3 and human G1/G3 group were 74.8–75.6%. However, RVA strains showed high nucleotide identity (92.5–99.9%) within P[6] genotype regardless of whether equine-like G3 or typical human G1/G3 type. Genotypes determined by phylogenetic analysis were consistent with those determined by RT-PCR.

FIGURE 3

Figure 3. Phylogenetic analysis of RVA VP7 (G genotype) gene. The tree was constructed using the neighbor-joining method. Bootstrap values (>70) are shown at the branch nodes. The names of the RVA strains detected in this study in 2015, 2016, 2017, 2018, and 2013–2014 are highlighted in blue, red, brown, purple, and green, respectively.

FIGURE 4

Figure 4. Phylogenetic analysis of RVA VP4 (P genotype) nucleotide sequences. The tree was constructed using the neighbor-joining method. Bootstrap values (>70) are shown at the branch nodes. RVA strains detected in this study are indicated by boldface.

The phylogenetic tree based on the NSP4 gene showed that eight strains were classified as bovine-like strain and the remaining 15 were classified as human strain (Supplementary Figure S8), while one (STM387) was not amplified. Those bovine-like strains were identified in 2012–2014 in Australia, in 2013 in Japan and in Indonesia (Malasao et al., 2015; Cowley et al., 2016; Utsumi et al., 2018). The representative RVA strains in the VP6, VP1–3, NSP1–NAP3, and NSP5 genes formed a cluster with those of equine-like DS-1-like strains from other countries (Supplementary Figures S1–S9).

Clinical Data and Severity of RVA Infection

The demographic and major clinical characteristics associated with disease severity are shown in Table 4. The typical symptoms associated with acute viral gastroenteritis are diarrhea, vomiting, fever, and dehydration. The prevalences of vomiting, watery stool, and some dehydration were significantly greater in the RVA-positive children than in the RVA-negative children (Table 4). An increase in the frequency of vomiting (≧5 in 24 h) was highly associated with the RVA-positive children (54.2%), though the difference with the RVA-negative children was not statistically significant (Table 4).

TABLE 4

Table 4. Clinical features observed among diarrheic children positive for rotavirus infection.

We analyzed the demographic and clinical features among children who experienced gastroenteritis with equine-like G3 strains or non-equine-like G3 strains (human G1 and G3). The differences in demographic and clinical characteristics (age, sex, duration, and frequency of diarrhea and vomiting, and dehydration) were not statistically significant between patients with equine-like G3 RVA infection and patients with non-equine-like G3 RVA infection (data not shown). Most of the patients infected with RVA were fully recovered or recovered, regardless of equine-like G3 or typical human G1/G3. There was no statistic differences in prognosis of infants and young children between equine-like and typical human RVA, although the number of samples were too small to compare.

Discussion

We previously reported that equine-like G3P[8]/[6] DS-1-llike RVA strain was predominant (100%) at a private hospital in the suburbs of Surabaya during 2015–2016 (Utsumi et al., 2018). Here we demonstrated that an unusual RVA genotype (equine-like G3P[8]/[6]) was also the predominant strain in this study population during 2015–2018, who were inpatients at the national referral hospital of East Java located in Surabaya, with decreased prevalence (87.6%). Two of the strains (D13 and D63), collected in 2013 in the same location in East Java and included in this study as controls, were classified into equine-like G3. Our results suggest that equine-like G3 strains, first emerged in Indonesia in 2013 at the same time as the first detection in Australia and Japan, have spread in the eastern part of Java Island, Indonesia. Phylogenetic analysis of the VP7 gene demonstrated that the equine-like G3P[8]/P[6] strains are highly homologous to other strains (99.9–97.3%) from Europe, Asia, South America, and Australia, which several research groups reported as equine-like RVA strains.

The equine-like G3P[8] DS-1 like strain was a rare human/equine reassortant and clustered with equine G3 strains at a high genetic distance by the phylogenetic analysis in VP7 gene and was named as equine-like G3P[8] (Cowley et al., 2016). The equine-like G3 genotype has been found in Australia, Asia, Europe, and South America since 2013, and it has been spreading worldwide, including newly reported cases from Germany (Pietsch and Liebert, 2018) and the United States (Perkins et al., 2017). Equine-like G3P[8]/[6] RVA strains in Indonesia were associated with large numbers of gastroenteritis cases as well as with cases in Australia and Spain (Arana et al., 2016; Cowley et al., 2016). The global emergence of equine-like G3 DS-1-like strains has raised a question whether vaccines have induced selective pressure on zoonotic strains (Luchs et al., 2019). It is important to determine if the global emergence of these novel reassortants may have any implications in RVA vaccine effectiveness.

Typical human G1P[8]/[6] and G3P[8]/[6] strains (12.4%), but not equine-like G3, have been found since July 2017 in the present study. It is noteworthy that the equine-like strains were exclusively detected until May 2017, but in July 2017 they were completely replaced by a typical human RVA (G1 or G3) (Figure 1 and Table 3). Not typical Wa-like strain, but unusual G1P[8] DS-1-like and equine-like G3P[8] DS-1-like strains co-circulated during 2013–2014, whereas equine-like G3P[8]/[6] DS-1-like distinctly circulated during 2015–2017 and G1/G3P[8]/[6] Wa-like/DS-1-like strains were found during July 2017- March 2018 (Figure 2). Now the genotypes has shifted to the typical human G1P[8] and G3P[8], which were previously prevalent in Indonesia. Genotype shift from dominant G1P[8] Wa-like to G2P[4] DS-1-like has been observed in the first 2 years after introduction of Rotarix^® (Gurgel et al., 2007; Zeller et al., 2010; Kirkwood et al., 2011; Mukhopadhya et al., 2017). Therefore, it is possible that introduction of RVA vaccine in 2011 in Indonesia resulted in genotype shift in Indonesia. However, there is no plausible explanation for this dynamic transition which occurred only in 2 months. To our knowledge, this is the first report of dynamic changes in rotavirus genotypes from equine-like G3 to typical human G1/G3.

Report from Japan demonstrated G1P[8] DS-1-like was the most prevalent strain during 2012–2014 under vaccination coverage of 39% in 2012–2013 (Fujii et al., 2019a). Another reports from Slovenia demonstrated the most prevalent genotype has shifted from G1P[8] to G4P[8] and G2P[4] under vaccination coverage of 27% in 2007–2013 (Steyer et al., 2014). Our finding also suggests that viral strain change can be influenced by the vaccination even in the country with low vaccine coverage like Indonesia. Further investigation may give us an insight to the direction of the global RVA epidemic in the future.

The overall prevalence of RVA infection in this study was 31.7%, lower than the previously reported data (40%) in East Java in 2013 (Sudarmo et al., 2015). Co-infection with bacteria was suspected in this cohort, since approximately 70% of pediatric patients with rotavirus were found to be co-infected with bacteria in 2013 in this government-owned hospital (Sudarmo et al., 2015). A previous report revealed that the overall prevalence of RVA infection in Java (Central and West), Bali, and West Nusa Tenggara, Indonesia, were 48.7% (2010–2015), 53.7% (2010), 53.3% (2011), 45.1% (2012), 41.7% (2013), 38.7% (2014), and 52.5% (2015) (Mulyani et al., 2018). These results demonstrate a gradual reduction in overall RVA infection by 2014 in those areas. In addition, the prevalence of RVA infection in this study (East Java) was similar to that in Central Java during 2006–2010 (Nirwati et al., 2016). Taken together, our results suggest that vaccination can affect prevalence of RVA infection even in the area with low vaccine coverage. RVA vaccination has not been included in the national immunization program in Indonesia, even though RVA vaccines have been commercially available since 2011 (Mulyani et al., 2018). None of the patients in this study had been immunized with an RVA vaccine. The RVA vaccine coverage in Indonesia is considered rather low since RVA vaccines are too expensive for ordinary citizens (Sudarmo et al., 2015). For this reason, in evaluating RVA vaccine efficacy, careful consideration on the low vaccination coverage should be given.

Introduction of RVA vaccination may be closely associated with the strain selection of RVA genotypes (Roczo-Farkas et al., 2016; Hoque et al., 2018). In addition, in Australia the equine-like RVA strain was found to be more common in areas where Rotarix^® vaccine was being administered (Roczo-Farkas et al., 2016). A new study targeting infants with RVA vaccines in private facilities of Surabaya revealed that approximately 70% of infants were immunized with the Rotarix^® vaccine (Gunawan et al., 2019). These findings raise the possibility that the Rotarix^® vaccine is associated with the emergence of the equine-like strains and that RVA vaccines may induce selective pressure that favors certain genotypes (Roczo-Farkas et al., 2016). The G2P[4] genotype that re-emerged after RVA vaccine introduction in certain areas (Roczo-Farkas et al., 2017; Khandoker et al., 2018) was not found in this study period, but in 2013 it was the predominant genotype in the same location in East Java (Sudarmo et al., 2015).

Infants and children with RVA-associated gastroenteritis exhibited some dehydration, vomiting, and watery diarrhea, which indicate the severity of RVA infection, more frequently than those with non-rotavirus gastroenteritis (Subekti et al., 2002; Mulyani et al., 2018). There were significantly high rates of infection among children younger than 2 years old, especially those between 6 and 11 months old. These findings are well in agreement with those of previous studies (Subekti et al., 2002; Sowmyanarayanan et al., 2012; Gupta et al., 2018). The lower prevalence in older children could be due to acquired immunity through previous exposures (Velazquez, 2009).

In the present study, RVA infection was detected continuously throughout a whole year without any seasonal peak. This was consistent with a study conducted in 2010–2015 in Indonesia (Mulyani et al., 2018) but inconsistent with a study from the same location in 2001–2004, showing peaks in the dry season (June to August) (Wilopo et al., 2009). Other temperate countries typically have RV peaks in winter (dry and cool conditions). These results suggest that considerable attention to local climatic conditions may help our understanding of the prevalence of RVA infection (Levy et al., 2009).

The severity of symptoms in RVA patients with equine-like G3 genotype is similar to that in RVA patients with the typical human RVA G1 and G3 in this study, inconsistent with studies conducted in Australia and Japan, indicating that the equine-like G3P[8] strain was associated with the most episodes of severe gastroenteritis (Cowley et al., 2016; Komoto et al., 2018). Since the number of global strains in this study was limited, it is necessary to accumulate more data on the severity according to each RVA genotype. Additionally, all of the pediatric patients with acute gastroenteritis in this study were hospitalized, suggesting that their symptoms were more severe than those of non-hospitalized cases. Since we failed to collect samples from non-hospitalized children, it is difficult to draw any conclusion about the severity of RVA-associated acute gastroenteritis among children depending on genotype. Further studies involving both hospitalized and non-hospitalized cases of rotavirus gastroenteritis are needed.

Strain diversity occurs over time along with temporal and genetic changes in the viruses (Gupta et al., 2018). The accumulation of knowledge about the features of viral strains may provide insight into what impact the introduction of rotavirus vaccination has on virus epidemiology and circulating genotypes. Thus, continuous monitoring of rotavirus genotypes, especially the emergence of any unusual rotavirus strains, along with data on prevalence and clinical features, will be crucial for assessing the disease burden and for achieving success in future vaccination programs in Indonesia.

Conclusion

A dynamic change in the rotavirus genotype from equine-like G3 to a typical human genotype (G1 or G3) was observed in this study. Since severe cases of RVA infection were prevalent especially in children aged 6 to 11 months and, more generally, those less than 2 years old, RVA vaccination should be included in Indonesia’s national immunization program. Careful monitoring on the emergence of unusual RVA strains is also needed.

Ethics Statement

Written informed consent was obtained from parents or guardians of all the children. The study protocol was reviewed and approved by the ethics committees of the government-owned hospital and Airlangga University in Indonesia, and Kobe University in Japan.

Author Contributions

AA was collected patient’s data. SS, RR, AD, and DR collected samples. TU was responsible for writing the manuscript. RW, ZD, and LY performed all experiments. So, J, and ML collected samples. CM, LD, and TA gave assistance for the research. YD, YF, HS, KK, and IS gave assistance for the research and analysis. IS supervised the study and helped to draft the manuscript. All authors critically revised the manuscript, read and approved the submitted version.

Funding

This research was supported by the Japan Initiative for Global Research Network on Infectious Diseases (J-GRID) from Ministry of Education, Culture, Sport, Science & Technology in Japan, Japan Agency for Medical Research and Development (AMED: 18fm0108004h0004), Shinryokukai.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We sincerely thank Prof. Dato’ Sri Tahir through Tahir’s professorship, and the Ministry of Research, Technology and Higher Education (RISTEKDIKTI) of Indonesia for their advice to carry out this study in Indonesia.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.00940/full#supplementary-material

Footnotes

1. ^ www.data.jma.go.jp/gmd/cpd/monitor/climatview/frame.php
2. ^ http://rotac.regatools.be/

References

Abou-Nader, A. J., Sauer, M. A., Steele, A. D., Tate, J. E., Atherly, D., Parashar, U. D., et al. (2018). Global rotavirus vaccine introductions and coverage: 2006 - 2016. Hum. Vaccin. Immunother. 14, 2281–2296. doi: 10.1080/21645515.2018.1470725

PubMed Abstract | CrossRef Full Text | Google Scholar

Arana, A., Montes, M., Jere, K. C., Alkorta, M., Iturriza-Gomara, M., and Cilla, G. (2016). Emergence and spread of G3P[8] rotaviruses possessing an equine-like VP7 and a DS-1-like genetic backbone in the Basque Country (North of Spain), 2015. Infect. Genet. Evol. 44, 137–144. doi: 10.1016/j.meegid.2016.06.048

PubMed Abstract | CrossRef Full Text | Google Scholar

Cowley, D., Donato, C. M., Roczo-Farkas, S., and Kirkwood, C. D. (2016). Emergence of a novel equine-like G3P[8] inter-genogroup reassortant rotavirus strain associated with gastroenteritis in Australian children. J. Gen. Virol. 97, 403–410. doi: 10.1099/jgv.0.000352

PubMed Abstract | CrossRef Full Text | Google Scholar

Doan, Y. H., Suzuki, Y., Fujii, Y., Haga, K., Fujimoto, A., Takai-Todaka, R., et al. (2017). Complex reassortment events of unusual G9P[4] rotavirus strains in India between 2011 and 2013. Infect. Genet. Evol. 54, 417–428. doi: 10.1016/j.meegid.2017.07.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Dóró, R., Laszlo, B., Martella, V., Leshem, E., Gentsch, J., Parashar, U., et al. (2014). Review of global rotavirus strain prevalence data from six years post vaccine licensure surveillance: is there evidence of strain selection from vaccine pressure? Infect. Genet. Evol. 28, 446–461. doi: 10.1016/j.meegid.2014.08.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Doro, R., Marton, S., Bartokne, A. H., Lengyel, G., Agocs, Z., Jakab, F., et al. (2016). Equine-like G3 rotavirus in Hungary, 2015 - Is it a novel intergenogroup reassortant pandemic strain? Acta Microbiol. Immunol. Hung. 63, 243–255. doi: 10.1556/030.63.2016.2.8

PubMed Abstract | CrossRef Full Text | Google Scholar

Elliott, E. J. (2007). Acute gastroenteritis in children. BMJ 334, 35–40. doi: 10.1136/bmj.39036.406169.80

PubMed Abstract | CrossRef Full Text | Google Scholar

Fujii, Y., Doan, Y. H., Suzuki, Y., Nakagomi, T., Nakagomi, O., and Katayama, K. (2019a). Study of complete genome sequences of rotavirus A epidemics and evolution in Japan in 2012-2014. Front. Microbiol. 10:38. doi: 10.3389/fmicb.2019.00038

PubMed Abstract | CrossRef Full Text | Google Scholar

Fujii, Y., Doan, Y. H., Wahyuni, R. M., Lusida, M. I., Utsumi, T., Shoji, I., et al. (2019b). Improvement of rotavirus genotyping method by using the semi-nested multiplex-PCR with new- primer set. Front. Microbiol. 10:647. doi: 10.3389/fmicb.2019.00647

PubMed Abstract | CrossRef Full Text | Google Scholar

Guerra, S. F., Soares, L. S., Lobo, P. S., Penha, E. T., Sousa, E. C., Bezerra, D. A., et al. (2016). Detection of a novel equine-like G3 rotavirus associated with acute gastroenteritis in Brazil. J. Gen. Virol. 97, 3131–3138. doi: 10.1099/jgv.0.000626

PubMed Abstract | CrossRef Full Text | Google Scholar

Gunawan, E., Utsumi, T., Wahyuni, R. M., Dinana, Z., Sudarmo, S. M., Shoji, I., et al. (2019). Post-vaccinated asymptomatic rotavirus infections: a community profile study of children in Surabaya, Indonesia. J. Infect. Public Health [Epub ahead of print].

PubMed Abstract | Google Scholar

Gupta, S., Krishnan, A., Sharma, S., Kumar, P., Aneja, S., and Ray, P. (2018). Changing pattern of prevalence, genetic diversity, and mixed infections of viruses associated with acute gastroenteritis in pediatric patients in New Delhi, India. J. Med. Virol. 90, 469–476. doi: 10.1002/jmv.24980

PubMed Abstract | CrossRef Full Text | Google Scholar

Gurgel, R. Q., Cuevas, L. E., Vieira, S. C., Barros, V. C., Fontes, P. B., Salustino, E. F., et al. (2007). Predominance of rotavirus P[4]G2 in a vaccinated population, Brazil. Emerg. Infect. Dis. 13, 1571–1573. doi: 10.3201/eid1310.070412

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoque, S. A., Kobayashi, M., Takanashi, S., Anwar, K. S., Watanabe, T., Khamrin, P., et al. (2018). Role of rotavirus vaccination on an emerging G8P[8] rotavirus strain causing an outbreak in central Japan. Vaccine 36, 43–49. doi: 10.1016/j.vaccine.2017

CrossRef Full Text | Google Scholar

Khandoker, N., Thongprachum, A., Takanashi, S., Okitsu, S., Nishimura, S., Kikuta, H., et al. (2018). Molecular epidemiology of rotavirus gastroenteritis in Japan during 2014-2015: characterization of re-emerging G2P[4] after rotavirus vaccine introduction. J. Med. Virol. 90, 1040–1046. doi: 10.1002/jmv.25067

PubMed Abstract | CrossRef Full Text | Google Scholar

Kirkwood, C. D., Boniface, K., Barnes, G. L., and Bishop, R. F. (2011). Distribution of rotavirus genotypes after introduction of rotavirus vaccines, Rotarix^® and RotaTeq^®, into the national immunization program of Australia. Pediatr. Infect. Dis. 30, S48–S53. doi: 10.1097/INF.0b013e3181fefd90

PubMed Abstract | CrossRef Full Text | Google Scholar

Komoto, S., Ide, T., Negoro, M., Tanaka, T., Asada, K., Umemoto, M., et al. (2018). Characterization of unusual DS-1-like G3P[8] rotavirus strains in children with diarrhea in Japan. J. Med. Virol. 90, 890–898. doi: 10.1002/jmv.25016

PubMed Abstract | CrossRef Full Text | Google Scholar

Komoto, S., Tacharoenmuang, R., Guntapong, R., Ide, T., Tsuji, T., Yoshikawa, T., et al. (2016). Reassortment of human and animal rotavirus gene segments in emerging DS-1-like G1P[8] rotavirus strains. PLoS One 11:e0148416. doi: 10.1371/journal.pone.0148416

PubMed Abstract | CrossRef Full Text | Google Scholar

Levy, K., Hubbard, A. E., and Eisenberg, J. N. (2009). Seasonality of rotavirus disease in the tropics: a systematic review and meta-analysis. Int. J. Epidemiol. 38, 1487–1496. doi: 10.1093/ije/dyn260

PubMed Abstract | CrossRef Full Text | Google Scholar

Luchs, A., da Costa, A. C., Cilli, A., Komninakis, S. C. V., Carmona, R. C. C., Boen, L., et al. (2019). Spread of the emerging equine-like G3P[8] DS-1-like genetic backbone rotavirus strain in Brazil and identification of potential genetic variants. J. Gen. Virol. 100, 7–25. doi: 10.1099/jgv.0.001171

PubMed Abstract | CrossRef Full Text | Google Scholar

Malasao, R., Saito, M., Suzuki, A., Imagawa, T., Nukiwa-Soma, N., Tohma, K., et al. (2015). Human G3P[4] rotavirus obtained in Japan, 2013, possibly emerged through a human-equine rotavirus reassortment event. Virus Genes 50, 129–133. doi: 10.1007/s11262-014-1135-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Matthijnssens, J., Ciarlet, M., Heiman, E., Arijs, I., Delbeke, T., McDonald S. M., et al. (2008). 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. 82, 3204–3219. doi: 10.1128/JVI.02257-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Matthijnssens, J., Heylen, E., Zeller, M., Rahman, M., Lemey, P., and Van Ranst, M. (2010). Phylodynamic analyses of rotavirus genotypes G9 and G12 underscore their potential for swift global spread. Mol. Biol. Evol. 27, 2431–2436. doi: 10.1093/molbev/msq137

PubMed Abstract | CrossRef Full Text | Google Scholar

Mukhopadhya, I., Murdoch, H., Berry, S., Hunt, A., Iturriza-Gomara, M., Smith-Palmer, A., et al. (2017). Changing molecular epidemiology of rotavirus infection after introduction of monovalent rotavirus vaccination in Scotland. Vaccine 35, 156–163. doi: 10.1016/j.vaccine.2016.11.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Mulyani, N. S., Prasetyo, D., Karyana, I. P. G., Sukardi, W., Damayanti, W., Anggraini, D., et al. (2018). Diarrhea among hospitalized children under five: a call for inclusion of rotavirus vaccine to the national immunization program in Indonesia. Vaccine 36, 7826–7831. doi: 10.1016/j.vaccine.2018.05.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Nirwati, H., Wibawa, T., Aman, A. T., Wahab, A., and Soenarto, Y. (2016). Detection of group A rotavirus strains circulating among children with acute diarrhea in Indonesia. Springerplus 5:97. doi: 10.1186/s40064-016-1724-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Perkins, C., Mijatovic-Rustempasic, S., Ward, M. L., Cortese, M. M., and Bowen, M. D. (2017). Genomic characterization of the first equine-like G3P[8] rotavirus strain detected in the United States. Genome Announc. 5:e1341-17. doi: 10.1128/genomeA.01341-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Pietsch, C., and Liebert, U. G. (2018). Molecular characterization of different equine-like G3 rotavirus strains from Germany. Infect. Genet. Evol. 57, 46–50. doi: 10.1016/j.meegid.2017

PubMed Abstract | CrossRef Full Text | Google Scholar

Roczo-Farkas, S., Kirkwood, C. D., Bines, J. E., Enteric Virus Group, Murdoch Childrens Research Institute, and Royal Children’s Hospital (2017). Australian rotavirus surveillance program: annual report, 2016. Commun. Dis. Intell. Q. Rep. 41, E455–E471.

Google Scholar

Roczo-Farkas, S., Kirkwood, C. D., Bines, J. E., and the Australian Rotavirus Surveillance Group (2016). Australian rotavirus surveillance program annual report, 2015. Commun. Dis. Intell. Q. Rep. 40, E527–E538.

Google Scholar

Rojas, M. A., Goncalves, J. L. S., Dias, H. G., Manchego, A., and Santos, N. (2017). Identification of two novel rotavirus A genotypes, G35 and P[50], from peruvian alpaca faeces. Infect. Genet. Evol. 55, 71–74. doi: 10.1016/j.meegid.2017.08.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Santos, N., and Hoshino, Y. (2005). Global distribution of rotavirus serotypes/genotypes and its implication for the development and implementation of an effective rotavirus vaccine. Rev. Med. Virol. 15, 29–56. doi: 10.1002/rmv.448

PubMed Abstract | CrossRef Full Text | Google Scholar

Soenarto, Y., Sebodo, T., Ridho, R., Alrasjid, H., Rohde, J. E., Bugg, H. C., et al. (1981). Acute diarrhea and rotavirus infection in newborn babies and children in Yogyakarta, Indonesia, from June 1978 to June 1979. J. Clin. Microbiol. 14, 123–129.

PubMed Abstract | Google Scholar

Sowmyanarayanan, T. V., Ramani, S., Sarkar, R., Arumugam, R., Warier, J. P., Moses, P. D., et al. (2012). Severity of rotavirus gastroenteritis in Indian children requiring hospitalization. Vaccine 30(Suppl. 1), A167–A172. doi: 10.1016/j.vaccine.2011.07.145

PubMed Abstract | CrossRef Full Text | Google Scholar

Steyer, A., Sagadin, M., Kolenc, M., and Poljšak-Prijatelj, M. (2014). Molecular characterization of rotavirus strains from pre- and post-vaccination periods in a country with low vaccination coverage: the case of Slovenia. Infect. Genet. Evol. 28, 413–425. doi: 10.1016/j.meegid.2014.06.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Subekti, D., Lesmana, M., Tjaniadi, P., Safari, N., Frazier, E., Simanjuntak, C., et al. (2002). Incidence of Norwalk-like viruses, rotavirus and adenovirus infection in patients with acute gastroenteritis in Jakarta, Indonesia. FEMS Immunol. Med. Microbiol. 33, 27–33. doi: 10.1111/j.1574-695X.2002.tb00568.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sudarmo, S. M., Shigemura, K., Athiyyah, A. F., Osawa, K., Wardana, O. P., Darma, A., et al. (2015). Genotyping and clinical factors in pediatric diarrhea caused by rotaviruses: one-year surveillance in Surabaya, Indonesia. Gut Pathog. 7:3. doi: 10.1186/s13099-015-0048-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Tamura, K., Dudley, J., Nei, M., and Kumar, S. (2007). MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol. Biol. Evol. 24, 1596–1599. doi: 10.1093/molbev/msm092

PubMed Abstract | CrossRef Full Text | Google Scholar

Utsumi, T., Wahyuni, R. M., Doan, Y. H., Dinana, Z., Soegijanto, S., Fujii, Y., et al. (2018). Equine-like G3 rotavirus strains as predominant strains among children in Indonesia in 2015-2016. Infect. Genet. Evol. 61, 224–228. doi: 10.1016/j.meegid.2018.03.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Velazquez, F. R. (2009). Protective effects of natural rotavirus infection. Pediatr. Infect. Dis. J. 28, S54–S56. doi: 10.1097/INF.0b013e3181967c03

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilopo, S. A., Soenarto, Y., Bresee, J. S., Tholib, A., Aminah, S., Cahyono, A., et al. (2009). Rotavirus surveillance to determine disease burden and epidemiology in Java, Indonesia, August 2001 through April 2004. Vaccine 27(Suppl. 5), F61–F66. doi: 10.1016/j.vaccine.2009.09.004

PubMed Abstract | CrossRef Full Text | Google Scholar

World Health Organization [WHO] (2014). Distance Learning Course Module 4: Diarrhoea. Integrated Management of Childhood Illness.

Zeller, M., Rahman, M., Heylen, E., De Coster, S., De Vos, S., Arijs, I., et al. (2010). Rotavirus incidence and genotype distribution before and after national rotavirus vaccine introduction in Belgium. Vaccine 28, 7507–7513. doi: 10.1016/j.vaccine.2010.09.004

PubMed Abstract | CrossRef Full Text | Google Scholar