Rotavirus gastroenteritis in Indian children < 5 years hospitalized for diarrhoea, 2012 to 2016

Background In 2016, the Government of India introduced the oral rotavirus vaccine (ROTAVAC, Bharat Biotech, India) in 4 states of India as part of the Universal Immunization Programme, and expanded to 5 more states in 2017. We report four years of data on rotavirus gastroenteritis in hospitalized children < 5 years of age prior to vaccine introduction. Methods Children from 7 sites in southern and northern India hospitalized for diarrhoea were recruited between July 2012 and June 2016. Stool samples were screened for rotavirus using enzyme immunoassay (EIA). The EIA positive samples were genotyped by reverse-transcription polymerase chain reaction. Results Of the 5834 samples from the 7 sites, 2069 (35.5%) were positive for rotavirus by EIA. Genotyping was performed for 2010 (97.1%) samples. G1P[8](56.3%), G2P[4](9.1%), G9P[4](7.6%), G9P[8](4.2%), and G12P[6](3.7%) were the common genotypes in southern India and G1P[8](36%), G9P[4](11.4%), G2P[4](11.2%), G12P[6](8.4%), and G3P[8](5.9%) in northern India. Conclusions The study highlights the high prevalence of rotavirus gastroenteritis in India and the diversity of rotavirus genotypes across different geographical regions. Pre- vaccine surveillance data is necessary to evaluate the potential change in admission rates for gastroenteritis and circulating rotavirus genotypes after vaccine introduction, thus assessing impact. Electronic supplementary material The online version of this article (10.1186/s12889-019-6406-0) contains supplementary material, which is available to authorized users.


Background
Rotavirus is a genus of the Reoviridae family, and is classified into 7 groups (A-G) [1]. Of the seven groups, group A rotavirus is the most important cause of severe acute diarrhoea in infants and young children globally [1]. Rotavirus consists of 11 segments of double stranded RNA, surrounded by 3 layers of proteins. The outer capsid consists of VP7 (G type) and VP4 (P type) proteins, which are essential for the binary system of rotavirus classification [1]. The four common G types (G1-G4), along with G9 and G12, in association with P [6] and P [8], comprise the common circulating rotavirus genotypes globally [1].
Rotavirus causes an estimated 11.37 million episodes of acute gastroenteritis (AGE) in children < 5 years annually in India, requiring 3.27 million outpatient visits and 872,000 inpatient admissions resulting in total direct costs of Indian Rupee (INR) 10.37 billion per year [2]. In 2011, it is estimated that there were 78,000 deaths due to rotavirus associated AGE in India, with the majority (75.6%, 59,000) in the first 2 years in life [2].
The Indian Rotavirus Strain Surveillance Network, comprising of 10 representative hospitals from 7 Indian cities, estimated rotavirus to cause 40% of hospitalizations due to AGE in children < 5 years during 2005-2009 [3]. The study demonstrated regional trends in rotavirus strain diversity along with the emergence of new strains [3]. A few other prospective hospital based surveillance studies in < 5 years children, conducted between 2009 to 2012, found rotavirus to be associated with approximately 26-39% of AGE cases [4,5]. The climate in north India varies from that in south India. The northern sites have more temperate climate, with a substantial fall in temperature during the winter months. In contrast, south India has a tropical climate. A difference in the seasonal trend in rotavirus associated diarrhoea across northern and southern sites in India has been reported in a few surveillance studies from India [4,6].
Pre-vaccine surveillance data is crucial to study the potential change in admission rates for gastroenteritis and circulating rotavirus genotypes after vaccine introduction. We report the findings of a multicentre hospital-based surveillance from July 2012 to June 2016 from 7 Indian sites on the clinical, epidemiological, and virological features of severe rotavirus disease among Indian children < 5 years of age, with the use of standardized protocols for enrolment and diagnostic evaluation.

Study sites
The study was conducted in hospitals in 7 cities of India, with 5 sites in southern India (Vellore, Trichy, Kolenchery, Hyderabad, Tirupati) and 2 in the north (Delhi, Ludhiana Stephen's Hospital (New Delhi) in the north. These hospitals are tertiary-level health care centres providing both general and specialized medical care. All the testing was performed at CMC, Vellore, using a modification of the World Health Organization (WHO) generic protocol for rotavirus surveillance [9]. Administrative and technical coordination was performed by the Indian Council of Medical Research (ICMR) and the National Institute of Epidemiology (NIE).

Enrolment criteria
Children ≤59 months of age hospitalized for at least 6 h for watery diarrhoea, and given supervised oral, or intravenous, rehydration, were eligible for enrolment. All the eligible children were enrolled after obtaining informed consent from the parents/guardian. An episode of AGE was defined as ≥3 loose stools over a 24 h period. Eligible patients were recruited if they provided a stool sample within 48 h of admission to rule out nosocomial infection. The exclusion criteria for the study were: children > 59 months of age, duration of diarrhoea > 5 days, blood in stool (dysentery), or < 3 episodes in a 24 h period.

Clinical assessment
Detailed clinical information on onset of diarrhoea, vomiting, fever, and dehydration were collected. The information included duration and number of diarrhoeal and/or vomiting episodes, degree of fever, severity of dehydration, and treatment. The Vesikari scoring system was used to assess the severity of diarrhoeal episodes (0-5 for mild, 6-10 for moderate, 11-15 for severe, and 16-20 for very severe) [10].

Sample collection and transport
One diarrhoeal stool sample was collected from each enrolled child. All stool samples were either transported to the testing laboratory within two hours of collection or stored at the sites at 4°C till transport. The stool samples from the sites other than Vellore were transported once every month in boxes with ice packs. On reaching the laboratory, the samples were tested for rotavirus antigen on the day of receipt, and aliquots of samples were stored at -70°C till further testing.

Laboratory procedures
A commercial enzyme immunoassay (EIA) was used to screen all the stool samples for rotavirus antigen (Rotaclone; Meridian Bioscience, Inc.). All EIA positive samples were further tested for VP7 (G type) and VP4 (P type) genes using reverse-transcription polymerase chain reaction (RT-PCR) according to the guidelines of the World Health Organization (WHO) for rotavirus detection and characterization [11]. In brief, viral RNA extraction was performed from 20% stool suspension (W/V) using the QIAamp Viral RNA Mini Kit (Qiagen). This was followed by reverse transcription for forming complementary DNA (cDNA) using Moloney murine reverse transriptase enzyme (Superscript II MMLV-RT, Invitrogen) and random primers (Invitrogen). The cDNA was used as the template for genotyping in a hemi-nested multiplex PCR for VP7 and VP4 genes, using published oligonucleotide primers to identify VP7 genotypes G1, G2, G3, G4, G8, G9, G10, G12, and VP4 genotypes P [4], P [6], P [8], P [9], P [10], and P[11] [3,6,12,13]. For samples negative by genotyping PCR, a VP6 gene specific PCR was performed to confirm rotavirus positivity [3,6]. Sanger sequencing was used to evaluate strains that could not be typed [3,6].

Statistical analysis
All sites submitted summary data on all diarrhoea-related hospital admissions. The forms were scrutinized for completeness, and the data were entered into Excel 2003 (Microsoft). The data were analyzed to evaluate the prevalence of rotavirus associated diarrhoea, genotype diversity, as well as temporal and regional variation in rotavirus strain prevalence from 2012 to 2016. To evaluate the prevalence of rotavirus associated diarrhoea, the proportion of rotavirus positive stool samples out of the total samples collected was calculated. SPSS software (version 21, IBM) was used to calculate pvalue using Fisher's exact test. A p value < 0.05 was considered to be statistically significant. For evaluating the regional variation in rotavirus positivity, the month-wise proportion of rotavirus positivity by EIA in the 2 north Indian sites and 5 south Indian sites were compared over the study period of 4 years.

Ethics
The study was approved by the institutional review board of CMC, Vellore, India (IRB Minute. No. 7311 dated 20.10.2010). The study was also approved by the institutional review boards/ethics committees of all the sites, which included St. Stephen's Hospital (New Delhi), SV Medical College (Tirupati), Christian Medical College (Ludhiana), Pragna hospital (Hyderabad), Child Jesus Hospital (Trichy), and Malankara Orthodox Syrian Church Medical College (Kolenchery). All the study participants were enrolled after obtaining written informed consent from the parents/guardian.

Site wise hospitalization due to rotavirus associated diarrhoea
During the four years of surveillance (July 2012-June 2016), 6576 eligible children hospitalized for AGE were enrolled at the 7 sites (Fig. 1, Table 1). A total of 5834  Table  S1). Stool samples of 35.5% (2069/5834) patients were positive for rotavirus by enzyme immunoassay (EIA). The overall rotavirus positivity was 41.7% (663/1590) in the north compared to 33.1% (1406/4244) in the south ( Table 1). The detection rate of rotavirus ranged from 23.5% in Tirupati to 49.4% in Trichy (Table 1).
Rotavirus-associated diarrhoea was seen throughout the year at all the sites. However, there was distinct seasonality for rotavirus associated diarrhoea, with highest prevalence seen during December-February (Fig. 2).

Characteristics of rotavirus-infected and uninfected children hospitalized with diarrhoea
We evaluated the demographic characteristics of children hospitalized with diarrhoea ( Table 2). There was no significant difference in the mean age (± standard deviation, SD) of children hospitalized with rotavirus-associated diarrhoea (15.0 ± 11.2 months), compared to rotavirus negative cases (14.9 ± 12.5 months). Children hospitalized with rotavirus associated gastroenteritis had more severe disease than children with rotavirus-negative diarrhoea ( Table 2). Children with rotavirus associated diarrhoea had significantly higher vomiting (76.1% versus 66.2%) than rotavirus negative cases ( Table 2). Infants (0-11 months) contributed the major proportion (47.9%) of rotavirus associated diarrhoea, followed by children aged 12 to 23 months (33.5%) and children ≥2 years (18.6%) ( Table 2).

G and P typing at 7 surveillance sites from 2012 to 2016
Of the 2069 samples that were positive for rotavirus by EIA, genotyping was performed for 2010 strains (97.1%); of which 1960 (97.5%) were G typed, and 1955 were P typed. Both G and P types were determined for 1947    (Table 3). We did not find any significant association between the common rotavirus genotypes and severity of diarrhoea (Additional file 2: Table S2). G12 strains were seen in combination with P [4], P [6], P [8] and P[11], and comprised 8.1% (162) of the total strains genotyped. The highest proportion of G12 strains were identified in New Delhi (12.8%) in the north and Kolenchery (10.7%) in the south (Table 3). G3P [8] strains were identified from 2015 onwards and were significantly more common in the north (6.2% in New Delhi, 5.4% in Ludhiana) compared to southern sites (ranged from 0.5% in Kolenchery to 3.3% in Trichy) (p < 0.001). Mixed infections were approximately 8% of all the genotyped strains, while partially typed and untyped samples comprised 1 and 2.1% of the genotyped samples respectively ( Table 3).
Mixed genotypes were more common in the north (10.6%) than the south (6.6%) (p = 0.003). The proportions of partially typed and untyped samples were significantly more common in the southern sites (1.3 and 2.5% in the south, compared to 0.5 and 1.2% in the north respectively) (p = 0.01).

Discussion
This 4 year study emphasizes the role of rotavirus as a major cause for acute gastroenteritis in children < 5 years in India before the introduction of the oral rotavirus vaccine into the national immunization programme. In more than 5800 stool samples, rotavirus was detected in 35.5%, comparable to the 2005-2009 period when rotavirus was detected in stools of 40% of children hospitalized for diarrhoea [3] An earlier review of studies on rotavirus disease in India found approximately 34% (inter study variation: 19-50%) of all AGE related hospitalizations to be due to rotavirus infections [7]. The study also demonstrates the early prevalence of rotavirus gastroenteritis in India, with most disease in the first 2 years of life. Although the percentage of rotavirus-associated diarrhoea remained similar from 2012 to 2016, significant changes in rotavirus strain diversity were observed across all sites. G1P [8] (49.9%) and G2P [4] (9.8%) were the two most common genotypes in this study, which is similar to the previous multi-centre surveillance studies conducted in India between 2005 to 2012 [3,4] Other studies from India have also found these two genotypes as common rotavirus strains causing diarrhoea in children (Table 4). Similar to the surveillance study from 2005 to 2009, an increase in G12 strains, particularly G12P [6], was observed in the southern (Kolenchery, 10.7%) and northern (New Delhi, 12.8%) regions [3] G12 strains were detected first in Philippines in 1987, following which they have been reported from other Asian countries [1,14,15] In India, G12 was reported first in 2003 from the eastern region, and since then it has been reported to cause AGE in children from other regions [3,13,16,17] G9, in association with P [4], P [6], and P [8], was detected in 13.2% of the samples, which is higher than reported in the Indian Rotavirus Network Surveillance study from 2005 to 2009 (7.5%) [3]. G9 was first detected in Philadelphia in 1983, and is now the fifth most common G type worldwide [1,15]. G9 is the most successful among the currently circulating reassortant human rotavirus strains. Similar to G9 rotaviruses,   G12 rotavirus strains also exhibit reassortment activity [15]. The emergence of G12 rotaviruses might be similar to the pattern of evolutionary events that led to the emergence of G9 as a dominant human rotavirus genotype. This view is strengthened by the fact that G12 rotavirus strains are increasingly being detected in Argentina, Nepal, Bangladesh, Japan, and other countries globally [15].
The emergence of G3P [8] in our study across southern and northern sites from 2015 onwards emphasizes the  rapid temporal changes in circulating rotavirus strains. The circulation of uncommon genotypes, such as G1P [4], G1P [6], G2P [6], G2P [8], G3P [4], G3P [6], G4P [6], G9P [6], G12P [4], and G12P [11], and the high proportion of mixed infections (8%) indicate that the children probably acquire infections from a variety of sources, and could serve as sources of new global strains. We did not find any association between the severity of diarrhoea and genotypes of the rotavirus, unlike previous studies which reported that children infected with G1 strains had a greater risk of developing severe diarrhoea than children infected with other rotavirus strains [5,18].
The strengths of this study include the use of a standardized protocol across the sites and laboratory confirmation of rotavirus diarrhoea in a single reference laboratory. The difficulty in extrapolating the data from 7 sites to the whole Indian population can be a potential limitation of the study, but previous studies have shown that data tend to be similar across the country [3].