- Research article
- Open Access
- Open Peer Review
Changing distribution of age, clinical severity, and genotypes of rotavirus gastroenteritis in hospitalized children after the introduction of vaccination: a single center study in Seoul between 2011 and 2014
BMC Infectious Diseasesvolume 16, Article number: 287 (2016)
This study aimed to explore changes in clinical epidemiology and genotype distribution and their association among hospitalized children with rotavirus gastroenteritis after the introduction of vaccines.
Between November 2010 and October 2014, hospitalized children with acute gastroenteritis were enrolled. Rotavirus genotypes were confirmed through reverse transcription-polymerase chain reaction (RT-PCR), semi-nested PCR, and sequencing. Clinical information including vaccination status and the modified Vesikari scores were collected.
Among 179 children with rotavirus infection, nineteen (10.6 %) were completely vaccinated. During the study period, the number of children between three and 23 months of age decreased significantly compared to the number of children older than 24 months of age (P = 0.010), who showed lower diarrhea severity (duration, P = 0.042; frequency, P = 0.021) but higher vomiting severity (P = 0.007, 0.036) compared to the former. Vaccination status was also significantly associated with lower vomiting severity after adjustment for age (frequency only, P = 0.018). The predominant genotypes were G2P (18.4 %), G1P (14.5 %), and G1PP (12.8 %), and the prevalence of genotypes with uncommon and mixed combinations was more than 50 %. Children infected with G2P strains tended to be older (P = 0.005) and had more severe vomiting (P = 0.018, 0.006) than those with G1P.
Increase in age of infected, hospitalized children was accompanied by change in clinical severity during 2011–2014 after the introduction of vaccines in Seoul. Clinical severity was also associated with vaccination status and genotype. Long-term large scale studies are needed to document the significance of the increase in genotypes of uncommon and mixed combinations.
Currently, rotavirus vaccination plays a key role in protecting children against rotavirus-related hospitalizations in countries implementing rotavirus vaccination [1–3]. In Korea, two rotavirus vaccines, RotaTeq (Merck and Company, Whitehouse Station, NJ, USA) and Rotarix (GlaxoSmithKline Biologicals, Rixensart, Belgium) were introduced in 2007 and 2008, respectively [1, 4]. Although rotavirus vaccines are recommended but not compulsory, the rate of individuals who had received two or three doses of a rotavirus vaccine was reported to be 65.6 % in Seoul and 52.4 % nationwide (2013 Korean National Immunization Survey, http://www.cdc.go.kr). As a result, the proportion of rotavirus infection among hospitalized children with gastroenteritis has been decreasing remarkably [2, 5]. According to a recent Korean Center for Disease Control report (http://www.cdc.go.kr), the incidence of rotavirus gastroenteritis was lower than that of norovirus in 2015. The effectiveness of vaccination against rotavirus-related hospitalization was suggested in six Asian countries including Korea . However, rotavirus gastroenteritis is still a major cause of illness requiring hospitalization even among rotavirus-vaccinated children. In addition, there is still debate on whether diversity of rotavirus strains contributes to low vaccine efficacy [7, 8]. Thus, in addition to efforts to increase vaccination rates, ongoing surveillance of rotavirus genotypes is necessary, particularly in children with vaccine failure .
Regarding the other notable changes in clinical epidemiology after the introduction of rotavirus vaccines, a few studies reported patterns of seasonal change or age change in addition to decreased overall prevalence of rotavirus infection [2, 5, 10]. A delay of the rotavirus season by 2–3 months was observed in Seoul in a multicenter study during 2007–2010 . A few single-center studies in the Korean literature reported that the age of children with rotavirus infection in Seoul or its suburbs increased significantly between 2008 and 2013 compared to the pre-vaccine era [5, 10]. It is unclear whether these changes have been consistent and ongoing. Before the introduction of rotavirus vaccines, rotavirus gastroenteritis usually occurred in children less than 24 months of age and presented with the most severe clinical manifestations including diarrhea among all types of viral gastroenteritis . Because severe diarrhea is usually believed to be prevalent in younger children , the shift of infected age to rotavirus infection might be accompanied by changes in clinical manifestation such as diarrheal severity of rotavirus gastroenteritis. In addition, although debatable, clinical factors such as age and severity of gastroenteritis have also been associated with certain genotypes according to several studies from the pre-vaccine era and, in countries without rotavirus vaccination programs in recent years [12–15]; thus, changes in clinical epidemiology might be accompanied by changes in molecular epidemiology. In previous studies, this association was often observed when the distribution of major circulating genotypes underwent changes .
Determination of circulating rotavirus genotypes is based on characterization of the capsid-protein-coding genes VP7 (G genotype) and VP4 (P genotype). To date, there are 27 G and 37 P genotypes have been identified (http://rotac.regatools.be/classificationinfo.html). Among them, five G genotypes (G1-G4, and G9) and three P genotypes (P, P, and P) are most frequently detected. Surveillance by the World Health Organization (WHO) reported that five combined G/P genotypes (G1P, G2P, G3P, G4P, and G9P) were considered globally prevalent during 2009-2012 . However, there is regional and temporal variation in the most prevalent genotypes. Studies have shown that 70–90 % of all typeable genotypes in Korea belong to one of the five most common combined G/P genotypes or G4P [16–20]. Between the years 1989 and 2009, the most prevalent genotype combination was G1P, though there were temporal changes in genotype predominance. For example, G1P was predominant during 1987–2000  and again between 2004 and 2009 [20, 22, 23]. G2P and G4P, which occurred mainly in neonates, were reported to be the most common strains during the intervening period, 2000–2003 . G9 was first detected in 2002 and was identified in consecutive seasons mostly as G9P [22, 23, 25], becoming the fifth most common genotype in Korea [17, 26]. After the introduction of rotavirus vaccines, G1P has remained predominant, however, the prevalence of G2P seemed to have increased around 2013 in Korea [27, 28]. especially in Seoul, compared to the early vaccination period, although only a few studies on the surveillance of genotypes are available in the literature [27–29].
Therefore, this study explored changes in the clinical epidemiology and genotype distribution among rotavirus-infected children from a single hospital in the southwestern area of Seoul between 2011 and 2014, a period corresponding to 4–7 years after the introduction of rotavirus vaccination. We also examined the correlations between genotypes and clinical factors including age, clinical severity and vaccination status.
Patients and clinical information
This study was conducted prospectively in the pediatric ward of Seoul Metropolitan Government-Seoul National University Boramae Medical Center in Korea with the approval of the Seoul Metropolitan Government-Seoul National University Boramae Medical Center Institutional Review Board (IRB No.02-2011-11, 26-2014-142). Between November 2010 to October 2014, stool samples were collected from children hospitalized with a clinical suspicion of acute gastroenteritis (AGE) without known underlying chronic illnesses with immunosuppression within three days of admission. The patients were admitted via either the emergency department or outpatient clinics. Informed written consent was obtained from the parents or guardians of each child at enrollment. Patients were excluded based on the following criteria: failure to obtain permission, insufficient fecal samples, children with newly diagnosed chronic illnesses including malnutrition, those with a final co-diagnosis as other acute febrile illnesses including pneumonia, or suspected cases of secondary infections (nosocomial gastroenteritis) during their hospital stay.
Clinical information including age of onset, date of hospitalization, fever, duration and maximum number of diarrhea episodes, and duration and maximum number of vomiting episodes before and during hospitalization and at previous visits before admission (outpatient or emergency department) were collected. From these records, modified Vesikari scores  were assessed. Laboratory tests for electrolytes were performed for all children. For most children (>95 %), a stool test for white blood cells and Salmonella species (spp.), Shigella spp., and Cholera spp. were also performed by the hospital laboratory. For children who tested positive for rotavirus stool antigen, a vaccination history was obtained from the individual vaccination cards and the Korean Center for Disease Control website (http://is.cdc.go.kr), which records histories as part of the national vaccination program.
Additionally, to compare the results of clinical data from historical data before the introduction of vaccines in our hospital, electronic medical records from November 2003 to October 2015 were retrospectively reviewed and analyzed for the following information: age, date and year, diagnosis, and results of a rotavirus antigen test.
Stool rotavirus test and genotyping
All stool samples were tested for rotavirus antigen by the hospital laboratory using the enzyme immunoassay kit Bioline rotavirus® (SD standard diagnostics, Youngin, Korea). Stool samples were diluted to 10 % with phosphate-buffered saline, centrifuged, and stored at −70 °C until further processing. For fecal specimens with positive stool antigen tests, viral dsRNA was extracted from stored samples using the QIAamp Viral RNA Mini Kit (Qiagen, Hilden, Germany). For genotyping, the extracted RNA was denatured and reverse transcription polymerase chain reaction (RT-PCR) was performed by using the Qiagen OneStep RT-PCR kit (Qiagene, Hilden, Germany). The VP7 gene was amplified by using primer sets Beg9 and End9 under conditions described previously [30, 31]. Second-round genotyping was analyzed using a pool of multiple primer sets aBT1, aCT2, aET3, aDT4, aAT8, aFT9 and G12, specific to G genotypes 1, 2, 3, 4, 8, 9 and 12, respectively [30–35]. To prevent known mistyping, additional semi-nested PCR for G3 and G9 typing was performed separately with different primers, G3 and G9 . Second-round PCR for G5 and G10 was performed with the reverse primers FT5 and G10 . The VP4 gene was amplified by using the primer sets Con3 and Con2 under conditions described previously . P genotypes were analyzed by semi-nested PCR using a pool of multiple primer sets 2T1, 3T1, 1T1, 4T1, 5T1, and P11, specific to P genotypes 4, 6, 8, 9, 10, and 11, respectively [32, 37]. For P, P and P genotyping, additional semi-nested PCR was performed separately with a mentioned primer sets. For strains that were not P-typeable by genotyping PCR as mentioned above, the primer sets VP4F and VP4R were used to amplify for the VP4 gene under previously described conditions . Amplification products were examined by electrophoresis in a 2 % agarose gel and documented with the Bio-Rad Gel Doc 1000 Documentation System (BioRad, Hercules, CA, USA).
For strains that were not G- or P-typeable by genotyping PCR, semi-nested PCR was performed with VP7 primers or VP4 primers after first round PCR [32, 38]. The resulting second-round amplicons were purified using the QIAquick Gel Extraction Kit (Qiagen) and directly sequenced using ABI Prism BioDye terminator cycle version 3.1 and an automatic DNA sequencer ABI 3730 (Applied Biosystems, Foster City, CA, USA). In addition, to validate the genotyping, one fecal sample corresponding to each of the G genotypes (G1,G2,G3,G4, and G9) and P genotypes (P, P and P) were sequenced using first round PCR products or second-round amplicons. The second-round amplicons of four randomly chosen strains with multiple P genotyping were also sequenced to confirm the results. The resulting nucleotide sequences were analyzed with the Bioedit software package. Genotypes were determined using the NCBI BLAST nucleotide search program in the GenBank database. The sequences of P untypeable strains were submitted to GeneBank under the accession numbers KR611086-KR611099.
We examined the data for normality by using the Kolmogorov-Smirnov test for normality. For age-related analyses, both the monthly age and age groups (four groups: <3 months, 3–23 months, 24–59 months, and ≥60 months) were used. For comparison among years, four seasons were defined like the following as the rotavirus epidemic starts around November: 2011 season (November 2010 to October 2011), 2012 season (November 2011 to October 2012), 2013 season (November 2012 to October 2013), and 2014 season (November 2013 to October 2014). For comparison among periods from retrospective data, three periods were defined like the following: pre-vaccine period (November 2003 to October 2006), early vaccine period (November 2007 to October 2010), and the post-vaccine period (November 2011 to October 2015). Comparisons of categorical data were evaluated using the Chi-square test. For age groups and clinical severity-related categorical data, P value for Linear-by-linear test were used. Continuous variables were summarized using the median and interquartile range (IQR). Comparisons of continuous data were evaluated using the Mann–Whitney test or Kruskal-Wallis test. For the adjustment of genotype-associated factors, data were analyzed using multiple logistic regression with genotypes as the dependent variables and season, monthly age, and vaccination status as explanatory variables. The analysis using ordered logit models was also performed to determine the cumulative odds ratio for clinical severity-related factors as dependent variables and monthly age, vaccination status, and genotype as explanatory variables. Ordinal logistic regression analysis was performed using SAS Version 9.2 (SAS Institute, Cary, NC, USA). All other analyses were performed using SPSS Version 20.0 (IBM, Chicago, IL, USA).
Overall prevalence, seasonal variation, and distribution of age
From 820 samples meeting the inclusion criteria during the study period, 179 samples (21.8 %) positive for rotavirus antigen were genotyped. The results of tests for other pathogen performed in the hospital laboratory were all negative for 179 samples.
From the retrospective data, the prevalence of rotavirus infection (352/2185 = 15.8 %, from November 2011 to October 2015) was decreased significantly compared to that of the pre-vaccine period (300/1152 = 25.9 %, from November 2003 to October 2006) and early-vaccine period (261/1367 = 18.9 %, from November 2007 to October 2010) (P = 0.000). Rotavirus infection was prevalent between January and May, peaking in March without significant variation among the four study seasons (Fig. 1). However, a delay in the peak season during the 2011–2014 seasons compared to the pre-vaccine period was observed (Additional file 1: Figure S1). The median age of the 179 children was 30.0 months [IQR 17.0–50.0], which was similar to the median age of the 312 children (28.6 months [IQR 14.0–48.0]) from the retrospective data during the 2011–2014 seasons. Of 179 children, 144 (80.4 %) were under five years of age. The rotavirus infection was more prevalent in children older than 24 months of age (n = 108) than those 3–23 months of age (n = 59) (Table 1).
The age of children infected with rotavirus was not significantly different between seasons during the study periods (P = 0.073). However, the number of children between three and 23 months of age decreased significantly over time compared to the number of children older than 24 months of age (P = 0.010) (Fig. 2). In addition, the median age of children with rotavirus infection significantly increased compared to the pre-vaccine or early vaccine periods (Additional file 2: Figure S2). The prevalence of the rotavirus infection compared to that of overall acute gastroenteritis was decreased mostly in children between three and 23 months of age, but was rather increased in children under 3 months of age and those over 60 months of age (Additional file 3: Figure S3).
Clinical manifestations by age
Clinical severity based on a modified Vesikari score is summarized in Table 1. According to age, the neonate group had a lower modified Vesikari score than older age groups (P = 0.001). They also had a shorter duration of vomiting (P < 0.001), fewer vomit episodes (P < 0.001), and shorter duration of diarrhea (P = 0.046) than older age groups. After excluding the neonate group (n = 167), children more than two years of age had a significantly longer duration of vomiting (P = 0.007) and a more vomiting episodes (P = 0.036) than children younger than two years of age.
For diarrhea-related severity, children less than two years of age, excluding neonates, had a significantly longer duration of diarrhea (P = 0.042) and more diarrheal episodes (P = 0.021) than children older than two years. There were no differences in other clinical severity-related factors between the groups.
Clinical manifestations by vaccination status
Excluding six children with incomplete vaccination records, 173 samples were analyzed for vaccination status. Of these children, 19 had completed rotavirus vaccination, 14 with RotaTeq and 5 with Rotarix (Table 1). The proportion of vaccinated children did not change over time and was significantly higher in children under five years of age (17.1 %, 19/111) than in children over 5 years (0 %). The median age of completely vaccinated children was 23.0 months [IQR 18.0–29.0], and the median age of unvaccinated children was 32.0 months [IQR 18.0–55.5] (P = 0.025).
As neonatal age significantly affected clinical severity-related factors, we analyzed data after excluding children under three months of age (n = 165). Vaccinated children had shorter durations of vomiting (P = 0.045) and fewer vomiting episodes (P = 0.012) (Table 2). There was no difference in diarrhea-related severity or other clinical severity including the overall Vesikari scores with vaccination status.
Because vomiting-related severity was also affected by age, the analysis using ordered logit models was performed. Vomiting duration was significantly associated with monthly age (cumulative odds ratio, c-OR, 1.017; 95 % CI, 1.005–1.029; P = 0.006). However, there was a significant association between the maximum number of vomiting episodes per 24 h and vaccination status (c-OR, 0.339; 95 % CI, 0.138-0.832; P = 0.018).
Overall prevalence and seasonal variation of rotavirus genotypes
The G genotype was identified in 97.2 % (n = 174) and P genotype in 97.8 % (n = 175) of samples. G1 was the most common G genotype (48.0 %, n = 86), followed by G2 (24.6 %, n = 44), G3 (20.1 %, n = 36), G9 (10.6 %, n = 19), and G4 (3.4 %, n = 6). P was the most prevalent P genotype (68.2 %, n = 122). P was found in 45.3 % (n = 81), and P in 3.9 % (n = 7). Co-infections with P and P genotypes were identified in 19.6 % of cases (n = 35). In G/P combinations, G2P was the most prevalent, found in 18.4 % (n = 33) of cases, followed by G1P (14.5 %, n = 26) and G1PP (13.4 %, n = 24), although the most common genotypes were different among seasons (Table 3).
Overall, the five most common WHO combinations were identified in 41.3 % of cases (n = 74). Conversely, genotypes with unusual P combinations including G1P, G3P and G9P were prevalent in 24.6 % of cases (n = 44). Interestingly, G3P and G9P were more prevalent than G3P and G9P. Co-infection with mixed G/P genotypes was detected in 28.5 % of cases (n = 51). The proportion of uncommon genotypes including mixed type was more than 40-50 %, which was consistent among the four seasons (Table 3). There were no significant differences in the frequencies of other single or multiple G, P, or G/P combinations between seasons.
Rotavirus genotypes by age
There was no age difference between common, uncommon and/or mixed genotypes. However some G, P or G/P genotypes showed an age preference (Table 4).
Children infected with G1 (median, 22.5 months; [IQR 12.5–34.5]) were younger than children with the G2 (median, 41.0; [IQR 25.0–61.9]; P = 0.001), and G9 (median, 39.0 months; [IQR 22.5–56.3]; P = 0.010) genotypes after exclusion of G and P type co-infection cases. Children infected with P (median, 31.5 months [IQR 21.5–53.8]) were older than children with P (median, 26.0 months [IQR 15.3–42.1; P = 0.049) after exclusion of G and P type co-infection cases. However, the age preference of the aforementioned genotypes was not significant after adjustment for seasons in logistic regression analysis.
Among G/P combinations, the most predominant genotypes were G1P and G1G4P in neonates (both 27.3 %, n = 3), G1P at 3–23 months (20.0 %, n = 12), G2P at 24–59 months (21.9 %, n = 16), and G2P and G1PP at ≥60 months (both 28.6 %, n = 10). The prevalence of G2P, G1PP, and untypeable type was significantly higher in older age groups (P = 0.002, 0.010, and 0.025 respectively). G1G4P was identified only in neonates (P < 0.001) (Table 2). The age of children infected with G2P (median, 41.5 months [IQR 25.0–62.75]) or G1PP (median, 44.5 months [IQR 28.5–78.5]) was significantly older than G1P-infected children (median, 22.25 months [IQR 13.88–30.63]; P = 0.001 and P = 0.005, respectively). Compared to G1P, this age preference for older children was still significant after adjusting for seasons and vaccination status in the G2P (adjusted odds ratio, a-OR of monthly age, 1.042; 95 % confidence interval, CI, 1.010–1.075; P = 0.011) and G1PP (a-OR, 1.043; 95 % CI, 1.009–1.077; P = 0.012) genotypes.
Rotavirus genotype by vaccination status
There were no significant differences between the two groups according to the vaccination status in the frequencies of common, uncommon, mixed genotypes, or combinations of G/P genotypes according to vaccination status (Additional file 4: Table S1).
However, among single G/P genotypes, G9 genotype was detected more frequently in the vaccinated group than in the unvaccinated group (n = 5/19; 27.8 % vs n = 14/154; 9.3 %, respectively; P = 0.020) (Additional file 5: Table S2). The association with G9 and vaccination status was still significant (a-OR of vaccination, 5.261; 95 % CI, 1.477–18.734; P = 0.010) after adjustment for seasons and age (a-OR of monthly age, 1.014; 95 % CI, 0.9996–1.028; P = 0.036) in multiple logistic regression.
Rotavirus genotypes and clinical severity
As neonatal age significantly affected clinical severity-related factors, we analyzed data after excluding children under three months of age (n = 167). There was no significant difference in clinical severity among common, uncommon and mixed genotypes. However, some clinical severity was significantly different among some G, P or G/P genotypes (Table 5).
Children (n = 167) infected with the G2 genotype presented with longer durations of vomiting (P = 0.030) and more vomiting episodes (P = 0.007) than those infected with G1, after exclusion of G/P co-infection cases. Children infected with the P genotype had longer durations of diarrhea (P = 0.025), but not more diarrheal episodes (P = 0.070) than those infected with the P genotype, after exclusion of G/P co-infection cases. There were no differences in clinical severity among other G or P genotypes.
For G/P combinations, we compared G1P and other major G/P combinations, including G2P. Children with G2P were older (P = 0.005) and had longer durations of vomiting (P = 0.018), and more vomiting episodes (P = 0.006) than children with G1P. Modified Vesikari scores did not show any significant differences between two genotypes (P = 0.077). Children with G1PP were older than those with G1P (P = 0.003). However, there were no differences in the clinical manifestations of G1PP and G1P (Table 5). Neither were there any significant differences in clinical manifestations between G1P and other G/P combinations, including G1P, G3P, and G9P.
As older age, unvaccinated status, and infection with the G2P genotype compared to G1P were all significantly associated with greater vomiting-related clinical severity, the analysis using ordered logit models was performed. Genotype was found to be a significant factor for both vomiting duration (c-OR of G2P4 genotype after adjustment for age, 4.89; 95 % CI, 1.242–19.258; P = 0.023) and the maximum frequency of vomiting (c-OR of G2P4 genotype after adjustment for vaccination status, 4.68; 95 % CI, 1.658–13.232; P = 0.004); monthly age (c-OR for vomiting duration, 1.015; 95 % CI, 1.001–1.028; P = 0.031) and vaccination status (c-OR for maximum frequency of vomiting, 0.184; 95 % CI, 0.066–0.514; P-0.001) also remained as significant factors for vomiting-related severity after adjustment for genotypes.
This study spanned the fourth to seventh years of surveillance of rotaviruses in a single hospital in southwestern area of Seoul after the introduction of rotavirus vaccines in Korea. Our report contains two notable findings; firstly, it demonstrated a significant increase in the age of hospitalized children with the rotavirus infection after the introduction of the rotavirus vaccination program in Seoul, Korea. We further showed that this shift in infected age could be accompanied by a change in vomiting and diarrheal severity. Secondly, a distinctive change in the distribution of rotavirus genotypes, including a prominent increase in genotypes of uncommon and mixed combinations, was observed during the study periods. A relative increase of the G2P genotype was also observed, which was significantly associated with age and clinical severity in this population. Although these are the results from a single hospital with a modest number of fecal samples, a considerable number of findings in clinical and molecular epidemiology share similar trends observed in other single center studies from Seoul and some multi-center Korean studies conducted in the post-vaccine period [2, 5, 10, 27–29]. To our knowledge, this is one of the first attempts to explore the link between age, clinical severity and genotypes in a region that adopted the rotavirus vaccination, which may contribute to an integrated understanding of the changing patterns in both clinical and molecular epidemiology for rotavirus gastroenteritis after the introduction of vaccines.
The finding of the shift of infected age tended to be more prominent the further in time it was separated from the year in which vaccines were introduced, which was consistent with the results of a few previous Korean studies [2, 10]. In a multicenter study, the prevalence of rotavirus infection compared to that of overall acute gastroenteritis was examined during 2007–2010. The ratio of children between less than 24 months old to more than 24 months old decreased from 1.2 (54.46 %/46.18 %) in 2007 to 0.7 (6.35 %/8.68 %) in 2010 . In another, a single hospital study performed in the suburbs of Seoul, the median age of children with rotavirus infection increased to 3.6 years in 2012 from 1.7 years in 2007 . Although less prominent, a similar shift in age from 14.9 months in the pre-vaccine period to 19 months in the post-vaccine period was also observed in a Greek study , in which the rotavirus vaccination rate was 25–30 %. The shift in infected age may partly reflect the protective effect of vaccines against severe rotavirus infection in vaccinated young children, who were the main focus of the newly introduced rotavirus vaccines in contrast to older children, who could not be vaccinated due to age restrictions or lack of information regarding the vaccines. However, the prevalence of the rotavirus infection, compared to that of overall acute gastroenteritis, actually increased significantly in children, especially those older than 60 months of age, compared to that of the pre-vaccine era, although it decreased significantly in the whole population (Additional file 3: Figure S3). We hypothesize that the changing distribution of the rotavirus genotypes, such as an increase of G2P and G1PP that show older age preference in our study, may also be associated with this possibly transient increase in the age group of infected children, in addition to the age group deviation of vaccine recipients.
Interestingly, we found a significant association between clinical severity and age partly due to shift in age, which could produce a sufficient number of relatively old age children (including approximately 20 % of children older than 60 months) comparable to that of younger age children. Vomiting-related severity increased and diarrhea-associated severity decreased with age in rotavirus gastroenteritis. Although it is controversial if and how clinical severity is associated with age [40, 41], this result is compatible with our speculation, as pediatricians in clinical practice, that vomiting has become the main symptom of rotavirus gastroenteritis, with diarrheal severity becoming less prominent in the last few years. This may be associated with an increase in patient age after the introduction of vaccines. In our study, clinical severity was also related to vaccination status, although the difference was not very impressive, as well as different genotypes in addition to age. On the basis of the distribution of these factors in the study subjects and the clinical severity-related factors examined, we speculate that either older (vomiting-related severity) or younger (diarrhea-related severity) age could be associated with greater clinical severity, or, there could be no association between age and overall severity (as measured by the modified Vesikari score) [40, 41].
Regarding genotype distribution, two notable changes were observed during the study period. First, G2P was the most prevalent genotype overall, though there was seasonal variation. The predominance of G2P, has been reported recently in many countries with or without vaccination including Korea [28, 42]. This change might be related to natural variation rather than vaccine pressure. Second, the prevalence of P combinations with G1, G3 and G9 (G1P, G3P and G9P), along with their mixed P genotype combinations G1PP, G3PP, and G9PP increased to more than 40 % among typeable strains, although the G1, G3, and G9 genotypes were mostly found in combination with P in the past . As a result, the prevalence of the five most common WHO combinations was less than 50 %, which was relatively consistent during all four seasons. Until recently, most studies in Korea identified the five most common WHO combinations plus G4P (mainly in neonates) in more than 70–90 % of typeable strains [16–20]. The genotypes of mixed G/P combinations were not prevalent in most studies. In a pre-vaccine era (2000–2007) study in Korea, 85.3 % were 5 common genotypes plus G4P, and uncommon genotypes and mixed genotypes were 11.8 % (5.82 % and 5.97 %, respectively) . In other meta-analysis during 1989–2009 in Korea, 81.8 % were common genotypes . However, a recent single hospital study in the northeastern area of Seoul reported a high prevalence of genotypes with non-WHO combinations and mixed G/P combinations similar to ours, in which the most prevalent G/P genotype was G9P . It has been known that the occurrence of mixed infection by more than two strains with different genotype combinations could encourage genotype reassortment and the appearance of strains with uncommon genotype combinations [45, 46]. Although increases in naturally occurring reassortment have often been reported , the unusual increase in uncommon and mixed genotypes observed in our study might reflect one of the responses to the introduction of rotavirus vaccines. Increases in previously uncommon G/P combinations, especially those not included in the current vaccines, have also been observed in other countries [9, 47–49]: In a recent study in Brazil, an increase of G3P was reported following the introduction of the G1P attenuated vaccine . In the United States, an increase of G12 and a relatively low vaccine efficacy against this genotype were reported . Although it has not been determined whether vaccine-associated selection pressure may play a role in these observations , studies have reported that increases in uncommon G/P combinations tend to coincide with the implementation of a rotavirus vaccination program [47–49]. Future large scale studies for prolonged periods, both covering clinical and molecular epidemiology, are needed to document the significance of the increase of genotypes of uncommon and mixed combinations in our population.
In our study, the distinctive difference among the major distribution of genotypes was not observed between unvaccinated children or those with vaccine failure: the prevalence of genotypes of uncommon or mixed type combinations or major G/P combinations such as G2P was not significantly different between the two groups. Although the proportion of G9 genotype was significantly higher in the vaccinated group, the significance of this observation could not be determined because of a relatively low significant P value (P = 0.010) in a small number of vaccinated children (n = 19). In previous two Korean studies, although the number of vaccinated children was also small, there was no significant difference in the genotype distribution between vaccinated children and those unvaccinated [27, 28].
We found a significant relationship between certain genotypes and clinical factors, including age and vomiting- or diarrhea-related severity. In addition to the well-known neonatal preference of G4P , we found that G2P, G1PP, and G9 strains preferentially infected older children, after adjustment for season and vaccination status. Among them, G2P was significantly associated with greater vomiting-related severity compared to other genotypes, including G1P, beyond neonatal age. It is not clear why children infected with strains with G2P showed increased vomiting severity compared to those infected with G1P. However, the association between G2P and age or clinical severity itself is mostly consistent with the results of several studies from the pre-vaccine era [12, 13, 51] and in countries without rotavirus vaccination [14, 15]. The literature suggests that the preference of G2P or G2 to infect older children and cause greater clinical severity might be a consequence of a G2 epidemic after a prolonged period of nonprevalence . Cross-protective immunity may be less effective against G2 than other common G types or G2 may possibly be more virulent [12, 52]. These speculations are supported by previous studies on the predominance of G2 in children with second episodes of rotavirus gastroenteritis  or adult outbreaks of rotavirus gastroenteritis .
In Korea, G2P was prevalent before 2002, with its incidence decreasing until recently, and was again prevalent in our study period [24, 28]. Children at ages that are vulnerable to rotavirus gastroenteritis (except rotavirus-naïve children) may be more susceptible to G2P stains than G1P, which was prevalent until recently. Due to rotavirus vaccination, which began in 2007, many younger children have protective immunity against most common and possibly uncommon genotypes. However, many older children, who usually did not receive rotavirus vaccination in the early vaccine period due to an age restrictions or a lack of information regarding the vaccine, may have insufficient protective immunity against G2P strains, probably because of a lack of type-specific immunity and low cross-protective immunity against G2P strains compared to other common genotypes, including G1P. Therefore, infection with G2P strains occur more often in older children than with G1P and other common genotypes, influenced by both the G2P epidemic and vaccines. This association between genotype and age could also be one of the explanations for the significantly older age of unvaccinated children compared to younger children with vaccine failure in addition to age group deviation of the vaccine recipients.
There were some limitations in our study. First, this is not a large-scale multicenter study. Only hospitalized children at a single hospital were included. With regard to genotype-related clinical severity, community-based comparisons may yield more prominent and accurate differences between genotypes . Second, the number of vaccinated children was small, so we could not conclude the genotypic comparison between vaccinated and unvaccinated groups. On the other side, the low proportion of vaccinated children in this study might imply the protective effect of rotavirus vaccine. Third, tests for co-infections with microbial agents other than Salmonella spp., Shigella spp., and Cholera spp. were not performed in our patients. However, a recent study reported that clinical severity was not affected by multiple co-infections with other enteric viruses , with rotavirus showing the most severe clinical manifestations . Finally, other possible risk factors for severity such as previous breastmilk feeding or day care center attendance were not considered. However, children with major potential risk factors such as malnutrition or underlying chronic illnesses affecting host immunity were not included in our study.
This prospective hospital-based study revealed that rotavirus vaccine showed fare efficacy in southwestern area of Seoul between 2011 and 2014, corresponding to the fourth to seventh years of rotavirus vaccination, although rotavirus vaccines were not included in the national vaccination program. Age of children with rotavirus gastroenteritis increased significantly after the introduction of vaccines among hospitalized children in Seoul, partly suggesting a protective effect of vaccines in recipients and a possible association between host age and virus genotypes. Clinical severity of rotavirus infection was associated with age, vaccination status, and certain genotypes in this population. Future large scale studies with long term surveillance of genotypes are needed to document the significance of the increase of genotypes of uncommon and mixed combinations.
a-OR, adjusted odds ratio; c-OR, cumulative odds ratio; IQR, interquartile range; RT-PCR, reverse transcription-polymerase chain reaction; spp, species
World Health Organization. Rotavirus vaccines. WHO position paper - January 2013. Wkly Epidemiol Rec. 2013;88:49–64.
Choi UY, Lee SY, Ma SH, Jang YT, Kim JY, Kim HM, et al. Epidemiological changes in rotavirus gastroenteritis in children under 5 years of age after the introduction of rotavirus vaccines in Korea. Eur J Pediatr. 2013;172:947–52.
Gurgel RQ, Alvarez Ade J, Rodrigues A, Ribeiro RR, Dolabella SS, Da Mota NL, et al. Incidence of rotavirus and circulating genotypes in Northeast Brazil during 7 years of national rotavirus vaccination. PLoS One. 2014;9:e110217.
Matthijnssens J, Joelsson DB, Warakomski DJ, Zhou T, Mathis PK, van Maanen MH, et al. Molecular and biological characterization of the 5 human-bovine rotavirus (WC3)-based reassortant strains of the pentavalent rotavirus vaccine, RotaTeq. Virology. 2010;403:111–27.
Sohn T, Lee C, Kim Y, Kang M, Kim S, Lee S, et al. Clinical and epidemiological study of 1,165 Hospitalized cases of rotaviral gastroenteritis before and after the introduction of rotavirus vaccine, 2006–2013. Korean J Pediatr Infect Dis. 2014;21:174–80.
El Khoury AC, Mast TC, Ciarlet M, Markson LE, Goveia MG. Projecting the effectiveness of RotaTeq(R) against rotavirus-related hospitalizations and deaths in six Asian countries. Hum Vaccin. 2011;7:506–10.
Pera FF, Mutepfa DL, Khan AM, Els JH, Mbewana S, van Dijk AA, et al. Engineering and expression of a human rotavirus candidate vaccine in Nicotiana benthamiana. Virol J. 2015;12:205.
Madhi SA, Cunliffe NA, Steele D, Witte D, Kirsten M, Louw C, et al. Effect of human rotavirus vaccine on severe diarrhea in African infants. N Engl J Med. 2010;362:289–98.
Jain S, Vashistt J, Changotra H. Rotaviruses: is their surveillance needed? Vaccine. 2014;32:3367–78.
Lee H, DY K, Kim J, Choi S. The epidemiological trend of rotavirus gastroenteritis in children in a single center from 2004 to 2012: a retrospective study. Korean J Pediatr Infect Dis. 2014;21:181–90.
Guarino A, Ashkenazi S, Gendrel D, Lo Vecchio A, Shamir R, Szajewska H. European Society for Pediatric Gastroenterology, Hepatology, and Nutrition/European Society for Pediatric Infectious Diseases evidence-based guidelines for the management of acute gastroenteritis in children in Europe: update 2014. J Pediatr Gastroenterol Nutr. 2014;59:132–52.
Bern C, Unicomb L, Gentsch JR, Banul N, Yunus M, Sack RB, et al. Rotavirus diarrhea in Bangladeshi children: correlation of disease severity with serotypes. J Clin Microbiol. 1992;30:3234–8.
Cascio A, Vizzi E, Alaimo C, Arista S. Rotavirus gastroenteritis in Italian children: can severity of symptoms be related to the infecting virus? Clin Infect Dis. 2001;32:1126–32.
Ahmed SF, Mansour AM, Klena JD, Husain TS, Hassan KA, Mohamed F, et al. Rotavirus genotypes associated with acute diarrhea in Egyptian infants. Pediatr Infect Dis J. 2014;33 Suppl 1:S62–8.
Sudarmo SM, Shigemura K, Athiyyah AF, Osawa K, Wardana OP, Darma A, et al. Genotyping and clinical factors in pediatric diarrhea caused by rotaviruses: one-year surveillance in Surabaya. Indonesia Gut Pathog. 2015;7:3.
Min BS, Noh YJ, Shin JH, Baek SY, Kim JO, Min KI, et al. Surveillance study (2000 to 2001) of G- and P-type human rotaviruses circulating in South Korea. J Clin Microbiol. 2004;42:4297–9.
Kang JO, Kilgore P, Kim JS, Nyambat B, Kim J, Suh HS, et al. Molecular epidemiological profile of rotavirus in South Korea, July 2002 through June 2003: emergence of G4P and G9P strains. J Infect Dis. 2005;192 Suppl 1:S57–63.
Huh JW, Kim WH, Yoon MH, Lim YH. Genotypic distribution of rotavirus strains causing severe gastroenteritis in Gyeonggi province, South Korea, from 2003 to 2005. Arch Virol. 2009;154:167–70.
Han TH, Kim CH, Chung JY, Park SH, Hwang ES. Genetic characterization of rotavirus in children in South Korea from 2007 to 2009. Arch Virol. 2010;155:1663–73.
Shim JO, Thai Than V, Ryoo E, Lim I, Yoon Y, Kim K, et al. Distribution of rotavirus G and P genotypes approximately two years following the introduction of rotavirus vaccines in South Korea. J Med Virol. 2013;85:1307–12.
Seo JK, Sim JG. Overview of rotavirus infections in Korea. Pediatr Int. 2000;42:406–10.
Le VP, Kim JY, Cho SL, Nam SW, Lim I, Lee HJ, et al. Detection of unusual rotavirus genotypes G8P and G12P in South Korea. J Med Virol. 2008;80:175–82.
Shim JO, Baek IH, Le VP, Ko EM, Seok WS, Uh Y, et al. Molecular characterization of rotavirus diarrhea among children in South Korea: detection of an unusual G11 strain. Arch Virol. 2011;156:887–92.
Song MO, Kim KJ, Chung SI, Lim I, Kang SY, An CN, et al. Distribution of human group a rotavirus VP7 and VP4 types circulating in Seoul, Korea between 1998 and 2000. J Med Virol. 2003;70:324–8.
Kang JO, Kim CR, Kilgore PE, Choi TY. G and P genotyping of human rotavirus isolated in a university hospital in Korea: implications for nosocomial infections. J Korean Med Sci. 2006;21:983–8.
Than VT, Jeong S, Kim W. A systematic review of genetic diversity of human rotavirus circulating in South Korea. Infect Genet Evol. 2014;28:462–9.
Chung JY, Kim MS, Jung TW, Kim SJ, Kang JH, Han SB, et al. Detection of Rotavirus Genotypes in Korea 5 Years after the Introduction of Rotavirus Vaccines. J Korean Med Sci. 2015;30:1471–5.
Kim JS, Kim HS, Hyun J, Kim HS, Song W, Lee KM, et al. Analysis of rotavirus genotypes in Korea during 2013: an increase in the G2P genotype after the introduction of rotavirus vaccines. Vaccine. 2014;32:6396–402.
Han T, Park S, Chung J, Hwang E. Distribution of human rotavirus genotypes in a tertiary hospital, Seoul, Korea during 2009–2013. Pediatric Infectioon Vaccine. 2015;22:81–90.
Gouvea V, Glass RI, Woods P, Taniguchi K, Clark HF, Forrester B, et al. Polymerase chain reaction amplification and typing of rotavirus nucleic acid from stool specimens. J Clin Microbiol. 1990;28:276–82.
Gomara MI, Cubitt D, Desselberger U, Gray J. Amino acid substitution within the VP7 protein of G2 rotavirus strains associated with failure to serotype. J Clin Microbiol. 2001;39:3796–8.
Iturriza-Gomara M, Kang G, Gray J. Rotavirus genotyping: keeping up with an evolving population of human rotaviruses. J Clin Virol. 2004;31:259–65.
Gouvea V, Santos N, Timenetsky MC. Identification of bovine and porcine rotavirus G types by PCR. J Clin Microbiol. 1994;32:1338–40.
Aladin F, Nawaz S, Iturriza-Gomara M, Gray J. Identification of G8 rotavirus strains determined as G12 by rotavirus genotyping PCR: updating the current genotyping methods. J Clin Virol. 2010;47:340–4.
Banerjee I, Ramani S, Primrose B, Iturriza-Gomara M, Gray JJ, Brown DW, et al. Modification of rotavirus multiplex RT-PCR for the detection of G12 strains based on characterization of emerging G12 rotavirus strains from South India. J Med Virol. 2007;79:1413–21.
Mitui MT, Chandrasena TN, Chan PK, Rajindrajith S, Nelson EA, Leung TF, et al. Inaccurate identification of rotavirus genotype G9 as genotype G3 strains due to primer mismatch. Virol J. 2012;9:144.
Gentsch JR, Glass RI, Woods P, Gouvea V, Gorziglia M, Flores J, et al. Identification of group A rotavirus gene 4 types by polymerase chain reaction. J Clin Microbiol. 1992;30:1365–73.
Simmonds MK, Armah G, Asmah R, Banerjee I, Damanka S, Esona M, et al. New oligonucleotide primers for P-typing of rotavirus strains: Strategies for typing previously untypeable strains. J Clin Virol. 2008;42:368–73.
Trimis G, Koutsoumbari I, Kottaridi C, Palaiologou N, Assimakopoulou E, Spathis A, et al. Hospital-based surveillance of rotavirus gastroenteritis in the era of limited vaccine uptake through the private sector. Vaccine. 2011;29:7292–5.
Albano F, Bruzzese E, Bella A, Cascio A, Titone L, Arista S, et al. Rotavirus and not age determines gastroenteritis severity in children: a hospital-based study. Eur J Pediatr. 2007;166:241–7.
Mathew A, Rao PS, Sowmyanarayanan TV, Kang G. Severity of rotavirus gastroenteritis in an Indian population: report from a 3 year surveillance study. Vaccine. 2014;32 Suppl 1:A45–8.
Matthijnssens J, Zeller M, Heylen E, De Coster S, Vercauteren J, Braeckman T, et al. Higher proportion of G2P rotaviruses in vaccinated hospitalized cases compared with unvaccinated hospitalized cases, despite high vaccine effectiveness against heterotypic G2P rotaviruses. Clin Microbiol Infect. 2014;20:O702–10.
Than VT, Kim W. Prevalence of rotavirus genotypes in South Korea in 1989–2009: implications for a nationwide rotavirus vaccine program. Korean J Pediatr. 2013;56:465–73.
Jeong HS, Lee KB, Jeong AY, Jo MY, Jung SY, Ahn JH, et al. Genotypes of the circulating rotavirus strains in the seven prevaccine seasons from September 2000 to August 2007 in South Korea. Clin Microbiol Infect. 2011;17:232–5.
Fischer TK, Page NA, Griffin DD, Eugen-Olsen J, Pedersen AG, Valentiner-Branth P, et al. Characterization of incompletely typed rotavirus strains from Guinea-Bissau: identification of G8 and G9 types and a high frequency of mixed infections. Virology. 2003;311:125–33.
Freitas ER, Soares CM, Fiaccadori FS, Souza M, Parente JA, Costa PS, et al. Occurrence of group A rotavirus mixed P genotypes infections in children living in Goiania-Goias, Brazil. Eur J Clin Microbiol Infect Dis. 2008;27:1065–9.
Neves MA, Pinheiro HH, Silva RS, Linhares AC, Silva LD, Gabbay YB, et al. High prevalence of G12P rotavirus strains in Rio Branco, Acre, Western Amazon, in the post-rotavirus vaccine introduction period. J Med Virol. 2015;88:782–9.
Payne DC, Selvarangan R, Azimi PH, Boom JA, Englund JA, Staat MA, et al. Long-term Consistency in Rotavirus Vaccine Protection: RV5 and RV1 Vaccine Effectiveness in US Children, 2012–2013. Clin Infect Dis. 2015;61:1792–9.
Soares LD, Guerra SDD, de Oliveira ADL, dos Santos FD, de Menezes EMDC, Mascarenhas JDP, et al. Diversity of Rotavirus Strains Circulating in Northern Brazil After Introduction of a Rotavirus Vaccine: High Prevalence of G3P Genotype. J Med Virol. 2014;86:1065–72.
Shim SY, Jung YC, Le VP, Son DW, Ryoo E, Shim JO, et al. Genetic variation of G4P rotaviruses: evidence for novel strains circulating between the hospital and community. J Med Virol. 2010;82:700–6.
Iturriza-Gomara M, Green J, Brown DW, Ramsay M, Desselberger U, Gray JJ. Molecular epidemiology of human group A rotavirus infections in the United Kingdom between 1995 and 1998. J Clin Microbiol. 2000;38:4394–401.
Griffin DD, Fletcher M, Levy ME, Ching-Lee M, Nogami R, Edwards L, et al. Outbreaks of adult gastroenteritis traced to a single genotype of rotavirus. J Infect Dis. 2002;185:1502–5.
Kang G. Rotavirus genotypes and severity of diarrheal disease. Clin Infect Dis. 2006;43:315–6.
We appreciate the statistical consultation provided by the Medical Research Collaborating Center at the Seoul National University College of Medicine, the Seoul National University Hospital, and Professor So Hee Oh from the Department of Education and Research at the Boramae Medical Center.
Source of funding
Chang JY received a grant from a Basic Science Research Program through the National Research Foundation of Korea (NRF), funded by the Ministry of Science, ICT & Future Planning (NRF-2012R1A1A3010058).
Availability of data and materials
The data cannot be shared as local Institutional Review Board has no policy to share the data without prior permission.
SJO and CJY designed the study, analyzed data and wrote the article. CJY, MJS, and KJS performed the acquisition of data. CJY and SS performed stool examinations, genotyping and sequencing. All authors drafted the article, revised it critically for important intellectual content, and approved the final article.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Seoul Metropolitan Government-Seoul National University Boramae Medical Center Institutional Review Board (IRB No.02-2011-11, 26-2014-142) approved this study. The written consent to participate was obtained from the parents or guardians of each child.
Monthly distribution of rotavirus infections in children during pre-vaccine (2004–2006), early-vaccine (2008–2010), and post-vaccine (2011–2015) periods. (JPG 272 kb)
Age distribution of rotavirus infections in children during pre-vaccine (2004–2006), early-vaccine (2008–2010), and post-vaccine (2011–2015) periods. (JPG 668 kb)
Changes of the rotavirus infection rates by ages during pre-vaccine (2004–2006), early-vaccine (2008–2010), and post-vaccine (2011–2015) periods. (JPG 1221 kb)
Distribution of group A rotavirus G and P genotypes according to vaccination status, n (%) (DOCX 45 kb)
Distribution of group A rotavirus G and P genotypes according to the vaccination status (exclude incomplete vaccination cases, n = 173) (DOCX 34 kb)