Rotavirus gastroenteritis in Pakistan, 2018: updated disease burden
BMC Infectious Diseases volume 21, Article number: 426 (2021)
Rotavirus A (RVA) is a significant cause of severe diarrheal illness and one of the common causes of death in children under the age of five. This study was aimed at detecting the prevalence of RVA in Pakistan after rotavirus vaccines were introduced. Fecal samples were obtained from 813 children from different hospitals in Rawalpindi and Islamabad, Pakistan, from January 2018 to December 2018. To obtain additional information from the parents / guardians of the children, a standard questionnaire was used.
Using an enzyme-linked immunosorbent assay kit (ELISA), rotavirus antigen was detected and ELISA positive samples were subjected to reverse transcription PCR (RT-PCR). The findings showed 22% prevalence of RVA in children with acute gastroenteritis (AGE) via ELISA and 21% prevalence via RT-PCR in children with AGE. There was no statistically significant difference between gender, age and RVA infections. The winter, spring and fall/autumn seasons were statistically significant for RVA prevalence.
The present study will provide post vaccine prevalence data for the health policy makers. The implementation of rotavirus vaccines, along with adequate nutrition for babies, clean water supply and maternal hygienic activities during infant feeding, is recommended. Furthermore, continuous surveillance is mandatory in the whole country to calculate the disease burden caused by RVA.
Rotavirus A (RVA) is one of the leading agents among all the AGE-causing pathogens, account for approximately 200,000 deaths each year among children under 5 years of age in underdeveloped nations . It has also been estimated that about half (49%) of these RVA fatalities occurred in four developing countries, one of which is Pakistan with rotavirus infants mortality rate 67.6 per 100,000 children . Four live-attenuated RVA vaccines (Rotarix™, RotaTeq™, Rotavac™ and RotaSiil™) have received prequalification from the World Health Organization (WHO) and are currently available worldwide, and over 100 countries have implemented one of these vaccines in their routine vaccination programs . Rotaviruses are members of the double stranded RNA family Reoviridae and are categorized into 9 confirmed (A-I) and one tentative species (J) on the basis of distinctive antigenic properties . The 70–75 nm diameter of the non-enveloped triple-layered viral particle has 11 segmented double-stranded RNA (dsRNA) genes, coding six structural (VP1-VP4, VP6, VP7) and five or six non-structural proteins (NSP1-NSP5 or NSP6) based on the strain .
In 2009, the World Health Organization (WHO) proposed the implementation of rotavirus vaccines in routine vaccination systems worldwide and, in particular, in countries with a high incidence of diarrhoea-related mortality below 5 years of age . The post-vaccination monitoring studies showed a substantial reduction in RVA prevalence in the countries where the Rotarix™ vaccine was widely used .
The government of Pakistan introduced RVA vaccine (Rotarix™) in the EPI (Expanded program on Immunization) plan of each of the four provinces in 2018 with the support of GAVI, the Vaccine Alliance. The effectiveness of RVA vaccines has yet to be determined in hospital based studies in Pakistan . According to a survey, the RVA vaccine could prevent 3061 deaths per year in Pakistan with the existing routine immunization program, with the highest coverage in the capital city of Islamabad, preventing an additional 1648 deaths per year . In Pakistan, only a few studies have been published on the prevalence of rotavirus infection, so the current study is intended to evaluate the prevalence of rotavirus infection among diarrheal patients admitted to the different health care hospitals in Rawalpindi district.
Material and methods
The said 1 year study 1st January 2018 to 31 December 2018 was conducted from various hospital of capital territory of Pakistan, Molecular Virology Laboratory COMSATS University and Comwave Institute Islamabad. This study was approved by the SUIT (Sarhad University of Science and Information Technology) and PIMS (Pakistan institute of medical sciences) hospital Islamabad, Pakistan “Ethical review committee”. Prior to sample collection an informed consent was obtained from parents /guardians of the children. All experiments were performed in accordance with relevant guidelines and regulations.
A total of 813 fecal samples were collected from various hospitals in Pakistan’s capital territories. Samples were collected from patients under 10 years of age and presented with signs and symptoms of watery diarrhea for 2 to 3 days in accordance with WHO standard case guideline. Patients over the age of ten that had a diarrheal complication that lasted more than 7 days, as well as those that had bloody diarrhea, were removed from the study . Table 1 shows the list of hospitals. A clean, sterile and leak proof plastic container about 1–2 ml or 1–2 g were used for collection of the samples and the sample were transferred to Molecular Virology Laboratory, COMSATS University, Islamabad for further processing.
RVA detection via ELISA
A 10% fecal suspension was prepared in 1 ml of PBS (Phosphate buffer Saline) and tested for rotavirus presence using a commercially available enzyme immunoassay (EIA) kit RIDASCREEN manufactured by R-Biopharm (Germany) as instructed by the manufacturer.
RVA RNA extraction and RT-PCR
A QIAamp Viral RNA mini kit by Qiagen (Hilden, Germany) was used to extract the viral RNA. Samples found to be positive for the ELISA test were confirmed for RVA outer capsid genes (VP7 and VP4) by reverse transcription-PCR (RT-PCR). The primers used were (Beg9 and End9) for VP7 and (VP4–1-17F and Con2Deg) for VP4 . The RT-PCR was carried out under the following conditions: initial reverse transcription step (30 min at 45 °C), polymerase activation (95 °C for 5 min), 40 cycles of amplification (denaturation: 30s at 94 °C; annealing: 45 s at 45 °C for VP4 and 45 s at 50 °C for VP7; product extension for 1 min at 72 °C) and final extension of 10 min at 72 °C.
The PCR product was validated using 2% agarose gel electrophoresis, stained with ethidium bromide and visualized under ultraviolet illumination.
Statistical analyses were carried out using SPSS (Statistical Package for Social Sciences) software (v21.0) . A chi-square test was used to examine the possible relationship between gender difference, age difference, residence, and seasonal distribution and RVA positive and negative status. Prevalence of RVA was determined by dividing the number of cases by the total number of children examined during the study period. The statistical level of significance was set at<0.05.
A total of 175/813 (22%) were found positive via ELISA. Of the 175 RVA antigen-positive samples, 170/175 (97%) were positive by RT-PCR and the remaining 3% were negative by RT-PCR (Table 2).
Demographic features of children with RVA diarrhea
Demographic characteristics of 175 RVA positive fecal samples of paediatric diarrheal patients were examined. A total 566/813 (70%) patients were male while 247/813 (30%) were female. Of the 170 samples positive by RT-PCR, 57 are females (33.5%) and 113 are males (66.5%). However the statistical analysis showed that gender was non-significantly associated with rotavirus infection as shown in the Table 3. All the patients were distributed in four age groups as age group 1 (Baby) (patients of age 1 day to 12 months), age group 2 (Toddler) (age above 1 year to 3 years), age group 3 (Pre-school) (above 3 years to 5 years of age) and age group 4 (Grade School) (age above 5 to 10 years). The non-significant statistical results were obtained in case of rotavirus among different age groups as shown in the Table 3 (p > 0.05). A total of 228/813 (28%) were from rural population while 585/813 (72%) were from urban population. Among the rural residential population 19/228 (8%) were infected while 151/585 (26%) were found infective among urban population residential. The statistical analysis shows that residence of the patients is significant with reference to rotavirus infection as shown in the Table 3 (p<0.05).
Seasonality of RVA infection
A total of 54/813 patients were studied during Winter (December 1st to February 28th), 342/813 during Spring (March 1st to May 31th), 380/813 during Summer (June 1st to August 31th) and 37/813 during Autumn (September 1st to November 30th). The total infected patients during each season and there statistical analysis shows that winter, spring and fall/autumn seasons were statistically significant for RVA prevalence (p<0.05). However summer season was not significantly correlated with the presence of rotavirus infection (p > 0.05) (Table 3).
Rotavirus infection associated diarrhea among children is one of the leading cause of death in developing countries of the world including Pakistan . Besides availing the support of Global Alliance for Vaccine and Immunization (GAVI) and World Health Organization (WHO) in controlling the communicable disease burden like diarrhea, the health system of Pakistan is still not developed. Keeping the above situation in view the current study was designed to evaluate the prevalence of RVA infections among the diarrheal patients of Capital territory of Pakistan during 2018.
In the current study the prevalence of RVA was found 21% via RT-PCR which is comparable with a community based study from Karachi, Pakistan reported 17% of prevalence of RVA among diarrheal patients . This rate of prevalence is lesser than 27% reported from Capital territory of Pakistan by  and previous studies from Pakistan reported (29–34%) prevalence rate during 2008–2014 [14,15,16,17]. Some of the previous studies reported a higher prevalence rate of 63 and 57% from Faisalabad and Karachi respectively [18, 19].
In the present study the RT-PCR was failed to detect RVA among 5 fecal samples positive through ELISA. There may be several reason behind this such as ELISA false positive results, shipment and storage of the samples, degradation in RNA or mutation in gene that could not be detected through general primers and inclusion of older children as most are negative for rotavirus.
In the current study, the prevalence among rural residential population is 19/228 (8%) while 151/585 (26%) were found to be infected among urban residential populations. There is no difference in vaccine coverage between rural and urban areas. People in urban areas have better access to medical care than those in rural areas. As a result, population density and increased transmission may be to blame for the disparity in prevalence rates between the two regions. A pre-vaccination study in Maputo (Southern Mozambique) showed similar results for higher prevalence in urban areas and lower prevalence in rural areas .
Weather patterns have indirect but significant impact on the epidemiology of human rotavirus infections. Low indoor relative humidity and indoor associated conditions crowds can be important factors in the epidemiology of rotavirus disease. Hospitalizations for rotavirus gastroenteritis tend to be more common after cold, dry and warm, rainy months of the year . The said study observed RVA based diarrheal infection throughout the year with significant correlated during fall, winter and autumn. Similar supporting results of circulating RVA throughout the year with peak season during dry months and winter season have already been observed in other study from Pakistan [8, 15,16,17, 19]. The similar pattern of RVA seasonality is observed in other countries of the world [21,22,23].
Non-significant gender-wise distribution of RVA infection was found in the said study which is comparable with other studies from various areas of the same country Pakistan [8, 13, 17]. Toddlers were the most infected population with 39.4% in the current study supported by several studies which shows that the children below 4 months of age are protected from severe RVA associated diarrhea by maternal antibodies in breast feeding which might be reason for higher incidence rate of RVA in children having age above 7 months . It may also be suggested that the toddler group was not vaccinated so that the prevalence of RVA was high in this group. In contrast to our study 71.4% of patients having age below 1 years were reported by several studies from Pakistan by [8, 13, 16,17,18]. The rate of incidence of RVA was low among the children of 3 to 5 years of age which is comparable with sadiq et al., 2019. A similar supporting trends were observed in Bangladesh, India, Vietnam, kenya [25,26,27,28,29]. Low incidence has been seen in older age groups because higher age groups may develop immunity after each infection .
This study had many drawbacks, like several studies that were carried out in the past. Pakistan has implemented rotavirus vaccine in its vaccination program in January 2018. In this post vaccination study, RVA genotyping and other causes of diarrhoea among the respondents were not examined and the burden of rotavirus disease from this area is not representative of the burden throughout the country. In future studies, the inclusion of several surveillance sites representing the entire country and the monitoring of RVA genotypes would include the estimation of the RVA-related disease burden and the impact of the implementation of vaccines on the diversity of RVA genotypes.
To conclude, rotavirus is an important aetiology of acute gastroenteritis in Pakistan among children under 5 years of age. Our study highlights that rotavirus in Pakistan during 2018 causes a high burden of disease in children under 5 years of age. In order to avoid rotavirus infection and to lessen the burden of disease, effective rotavirus vaccination may be considered in Pakistan. It also underscores the importance of hygiene and sanitation in the country. Studies such as this provide evidence that will allow healthcare professionals to establish a more comprehensive diagnostic program that will take into account the potential for viral diseases. Further research to discover the predominant RV genotypes in Pakistan is also required for the future production of vaccines.
Availability of data and materials
The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.
Tate JE, Burton AH, Boschi-Pinto C, Parashar UD. World health Organization–coordinated global rotavirus surveillance network. global, regional, and national estimates of rotavirus mortality in children <5 years of age, 2000–2013. Clin Infect Dis. 2016;62(suppl 2):S96–105. https://doi.org/10.1093/cid/civ1013.
Sadiq A, Bostan N, Bokhari H, Matthijnssens J, Yinda KC, Raza S, et al. Molecular characterization of human group a rotavirus genotypes circulating in Rawalpindi, Islamabad, Pakistan during 2015-2016. PLoS One. 2019;14(7):e0220387. https://doi.org/10.1371/journal.pone.0220387.
Sadiq A, Bostan N. Comparative analysis of g1p [ 8 ] rotaviruses identi fi ed prior to vaccine implementation in pakistan with Rotarix ™ and RotaTeq ™ Vaccine Strains. Front Immunol. 2020;11(October):2–11.
Lestari FB, Vongpunsawad S, Wanlapakorn N, Poovorawan Y. Rotavirus infection in children in Southeast Asia 2008-2018: disease burden, genotype distribution, seasonality, and vaccination. J Biomed Sci. 2020;27(1):1–19. https://doi.org/10.1186/s12929-020-00649-8.
Matthijnssens J, Ciarlet M, Heiman E, Arijs I, Delbeke T, Mcdonald SM, et al. Full genome-based classification of rotaviruses reveals a common origin between human Wa-like and porcine rotavirus strains and human DS-1-like and bovine rotavirus strains. J Virol. 2008;82(7):3204–19. https://doi.org/10.1128/JVI.02257-07.
João ED, Munlela B, Chissaque A, Chilaúle J, Langa J, Augusto O, et al. Molecular epidemiology of rotavirus a strains pre-and post-vaccine (Rotarix®) introduction in Mozambique, 2012–2019: emergence of genotypes G3P and G3P. Pathogens. 2020;9:1–16.
Mukhopadhya I, Murdoch H, Berry S, Hunt A, Iturriza-Gomara M, Smith-Palmer A, et al. Changing molecular epidemiology of rotavirus infection after introduction of monovalent rotavirus vaccination in Scotland. 2017.
Sadiq A, Id NB, Bokhari H, Matthijnssens J. Molecular characterization of human group A rotavirus genotypes circulating in Rawalpindi , Islamabad , Pakistan during 2015–2016. 2019;:1–28.
Rheingans R, Anderson JD, Bagamian KH, Laytner LA, Pecenka CJ, Gilani SSA, et al. Effects of geographic and economic heterogeneity on the burden of rotavirus diarrhea and the impact and cost-effectiveness of vaccination in Pakistan. Vaccine. 2018;36(51):7780–9. https://doi.org/10.1016/j.vaccine.2018.02.008.
WHO. Global networks for surveillance of rotavirus gastroenteritis, 2001–2008. Wkly Epidemiol Rec. 2009;83:421–5 http://www.ncbi.nlm.nih.gov/pubmed/19024780. Accessed 16 Jun 2019.
IBM Corp. IBM SPSS Statistics for Windows, Version 21.0. Armonk: IBM Corp; 2012.
Hussain A. Rotavirus vaccine demand in Pakistan. Iran J Public Health. 2016;45:1658–9.
Kazi AM, Warraich GJ, Qureshi S, Qureshi H, Khan MMA, Zaidi AKM, et al. Sentinel hospital-based surveillance for assessment of burden of rotavirus gastroenteritis in children in Pakistan. PLoS One. 2014;9(10):e108221. https://doi.org/10.1371/journal.pone.0108221.
Afzal A, Akhtar Tariq P, Choudhry S. Rota virus gastroenteritis in children Upto five years of age. 2010. http://www.journalrmc.com/volumes/1394531488.pdf. Accessed 24 Feb 2019.
Alam MM, Khurshid A, Shaukat S, Suleman RM, Sharif S, Angez M, et al. Epidemiology and genetic diversity of rotavirus strains in children with acute gastroenteritis in Lahore, Pakistan. PLoS One. 2013;8:2–9.
Tamim S, Hasan F, Matthijnssens J, Sharif S, Shaukat S, Alam MM, et al. Epidemiology and phylogenetic analysis of VP7 and VP4 genes of rotaviruses circulating in Rawalpindi, Pakistan during 2010. Infect Genet Evol. 2013;14:161–8. https://doi.org/10.1016/j.meegid.2012.10.009.
Umair M, Abbasi BH, Nisar N, Alam MM, Sharif S, Shaukat S, et al. Molecular analysis of group a rotaviruses detected in hospitalized children from Rawalpindi, Pakistan during 2014. Infect Genet Evol. 2017;53:160–6. https://doi.org/10.1016/j.meegid.2017.05.009.
Iftikhar T, Butt A, Nawaz K, Sarwar Y, Ali A, Mustafa T, et al. Genotyping of rotaviruses detected in children admitted to hospital from Faisalabad region, Pakistna. J Med Virol. 2012;84(12):2003–7. https://doi.org/10.1002/jmv.23402.
Habib MI, Kazi SG, Ahmed Khan KM, Zia N. Rota virus diarrhea in hospitalized children. J Coll Physicians Surg Pak. 2014;24(2):114–7. doi:02.2014/JCPSP.114117.
de Deus N, João E, Cuamba A, Cassocera M, Luís L, Acácio S, et al. Epidemiology of rotavirus infection in children from a rural and urban area, in Maputo, southern Mozambique, before vaccine introduction. J Trop Pediatr. 2018;64(2):141–5. https://doi.org/10.1093/tropej/fmx032.
Agbla JMM, Capo-Chichi A, Agbankpé AJ, Dougnon TV, Yadouleton AWM, Houngbégnon O, et al. Epidemiological survey of rotaviruses responsible for infantile diarrhea by the Immunomolecular technique in Cotonou (Benin, West Africa). Int J Microbiol. 2018;2018(7):3602967. https://doi.org/10.1155/2018/3602967.
Bonifacio J, Lupisan S, Roque V, Ducusin MJ, Grabovac V, Batmunkh N, et al. Molecular characterization of rotavirus diarrhea among children aged under five years in the Philippines, 2013–2015. Vaccine. 2018;36(51):7888–93. https://doi.org/10.1016/j.vaccine.2018.08.046.
Myat TW, Thu HM, Kyaw YM, Aye KM, Win MM, Lin H, et al. Sentinel surveillance for rotavirus in children <5 years of age admitted for Diarrheal illness to Yangon Children’s Hospital, Myanmar, 2009–2014. Vaccine. 2018;36:7832–5. https://doi.org/10.1016/j.vaccine.2017.11.002.
Dhital S, Sherchand JB, Pokhrel BM, Parajuli K, Shah N, Kumar Mishra S, et al. Molecular epidemiology of rotavirus causing diarrhea among children less than five years of age visiting national level children hospitals, Nepal. BMC Pediatr. 2017;17(1):101. https://doi.org/10.1186/s12887-017-0858-0.
Haque W, Haque J, Barai D, Rahman S, Moni S, Hossain ME, et al. Distribution of rotavirus genotypes in Dhaka, Bangladesh, 2012–2016: re-emergence of G3P after over a decade of interval. Vaccine. 2018;36(43):6393–400. https://doi.org/10.1016/j.vaccine.2018.08.081.
Giri S, Kumar CPG, Khakha SA, Chawla-Sarkar M, Gopalkrishna V, Chitambar SD, et al. Diversity of rotavirus genotypes circulating in children < 5 years of age hospitalized for acute gastroenteritis in India from 2005 to 2016: analysis of temporal and regional genotype variation. BMC Infect Dis. 2020;20(1):740. https://doi.org/10.1186/s12879-020-05448-y.
Huyen DTT, Hong DT, Trung NT, Hoa TTN, Oanh NK, Thang HV, et al. Epidemiology of acute diarrhea caused by rotavirus in sentinel surveillance sites of Vietnam, 2012–2015. Vaccine. 2018;36(51):7894–900. https://doi.org/10.1016/j.vaccine.2018.05.008.
Fidhow AM, Samwel A, Ng’ang’a Z, Oundo J, Nyangao J, Wences A. Molecular epidemiology and associated risk factors of rotavirus infection among children < 5 yrs hospitalized for acute gastroenteritis in North Eastern, Kenya, 2012. Pan Afr Med J. 2017;28(Suppl 1):1–3.
Tian Y, Chughtai AA, Gao Z, Yan H, Chen Y, Liu B, et al. Prevalence and genotypes of group a rotavirus among outpatient children under five years old with diarrhea in Beijing, China, 2011-2016. BMC Infect Dis. 2018;18(1):497. https://doi.org/10.1186/s12879-018-3411-3.
We acknowledge the Higher Education Commission of Pakistan for providing funding for this work.
This work was supported by the Higher Education Commission of Pakistan.
Ethics approval and consent to participate
This study was approved by the SUIT (Sarhad University of Science and Information Technology) and PIMS (Pakistan institute of medical sciences) hospital Islamabad, Pakistan “Ethical review committee”. Prior to sample collection an informed consent was obtained from parents /guardians of the children. All experiments were performed in accordance with relevant guidelines and regulations.
Consent for publication
All authors do not mention any potential conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Basharat, N., Sadiq, A., Dawood, M. et al. Rotavirus gastroenteritis in Pakistan, 2018: updated disease burden. BMC Infect Dis 21, 426 (2021). https://doi.org/10.1186/s12879-021-06123-6
- Disease burden