Abstract
Background Rotavirus is the major cause of severe gastroenteritis in children aged <5 years. Introduction of Rotarix® rotavirus vaccine (RV1) in Malawi in 2012 has reduced rotavirus-associated hospitalisations and diarrhoeal mortality. However, RV1 impact on gastroenteritis severity remains unknown. We conducted a hospital-based surveillance study to assess RV1 impact on gastroenteritis severity in children aged <5 years, in Malawi.
Methods Stool samples were collected from children hospitalised with acute gastroenteritis from December 2011 – October 2019. Gastroenteritis severity was determined using Ruuska and Vesikari scores. Rotavirus was detected in stool using Enzyme Immunoassay. Rotavirus genotypes were determined using nested RT-PCR. Associations between RV1 vaccination and gastroenteritis severity were investigated using adjusted linear regression.
Results In total, 3,159 children were recruited. After adjusting for Mid-Upper Arm Circumference, age, weight, gender and receipt of other vaccines, all-cause gastroenteritis severity scores were 2.21 units lower (95% CI 1.85, 2.56; p<0.001) among RV1-vaccinated (n=2,224) compared to RV1-unvaccinated children (n=935); the decrease was comparable between rotavirus-positive and rotavirus-negative cases in all age groups. The reduction in severity score was observed against every rotavirus genotype, although the magnitude was smaller among those infected with G12P[6] compared to the remaining genotypes (p=0.011). Other than RV1 vaccination, age was the only variable associated with gastroenteritis severity. Each one-year increment in age was associated with a decrease of 0.43 severity score (95% CI 0.26, 0.60; p<0.001).
Conclusion Our findings provide additional evidence of RV1 impact in a high disease burden, low-income country, lending further support to Malawi’s rotavirus vaccine programme.
Summary In a long-term hospital-based surveillance study in Malawi, we found evidence of the reduction in gastroenteritis severity among hospitalised RV1-vaccinated children infected with both homotypic and heterotypic rotavirus strains and off-target RV1 vaccine effects against non-rotavirus diarrhoeal severity.
1. Introduction
Rotavirus is a leading cause of acute gastroenteritis among children worldwide. Despite the introduction of rotavirus vaccines in many countries, rotavirus is still associated with an estimated 128,500 deaths annually. Over 90% of these cases occur in low and middle-income countries (LMICs) in sub-Saharan Africa and South-east Asia [1].
Rotavirus has a segmented double-stranded ribonucleic acid (dsRNA) genome, surrounded by a triple layered capsid. Most human infections are associated with group A rotaviruses [2], which can further be classified using a dual classification system into G and P types, according to their neutralizing antibody response or nucleotide differences of the genes encoding their outer glycoprotein VP7 and protease-sensitive VP4, respectively [3]. At least 36 G types and 51 P types have been reported (https://rega.kuleuven.be/cev/viralmetagenomics/virus-classification). Of these, genotypes G1P[8], G2P[4], G3P[8], G4P[8], G9P[8] and G12P[8] are the most frequent cause of rotavirus disease in humans worldwide [4, 5].
Malawi introduced the G1P[] Rotarix® rotavirus vaccine (RV1) into its national Expanded Programme on Immunization (EPI) schedule on 28th October 2012, with doses administered at 6 and 10 weeks of age, which was followed by a decline in rotavirus-associated hospitalisations [6] and a reduction in gastroenteritis-related mortality [7]. However, the impact of RV1 on the severity of gastroenteritis has not yet been assessed. We conducted an analysis of the severity of gastroenteritis by comparing the Ruuska and Vesikari disease severity scores [8] in children presenting with rotavirus and non-rotavirus laboratory confirmed gastroenteritis before and after RV1 introduction, and between vaccinated and non-vaccinated children post RV1 introduction.
2. Materials and methods
2.2 Study population
Children under the age of five years who presented with acute gastroenteritis (defined as the passage of at least three looser-than-normal stools in a 24-hour period for less than seven days duration) were enrolled at the Queen Elizabeth Central hospital (QECH), Blantyre, which is the main referral hospital for the southern region of Malawi. To assess the impact of RV1 introduction on the severity of gastroenteritis, data were examined from children hospitalised with gastroenteritis before (December 2011 to October 2012) and after (November 2012 to October 2019) RV1 introduction. These populations were used to (i) assess the impact of RV1 on the severity of all-cause gastroenteritis; and (ii) determine whether gastroenteritis severity differed by rotavirus genotype in RV1-vaccinated and RV1-unvaccinated children.
2.3 Clinical and demographic variables
Gastroenteritis severity was determined using the Ruuska and Vesikari scoring system [8]. The assessment was based on the following parameters: duration and maximum number of episodes of diarrhoea as well as vomiting, fever, and dehydration. Scores 0 – 5 was considered as mild, 6 – 10 as moderate, 11 – 15 as severe and ≥16 as very severe. To assess the impact of RV1 by age, infants were categorised into four age groups (<6 months, 6 – 11 months, 12 – 23 months and 24 – 59 months).
2.4 Rotavirus detection and genotyping
A 10-20% stool suspension in diluent buffer was prepared for each specimen and used to screen for the presence of group A rotavirus using a commercially available Enzyme immunoassay (Rotaclone®, Meridian Bioscience, Ohio, USA). Rotavirus dsRNA was extracted from all rotavirus-positive stool samples using the Viral RNA Mini-Kit (Qiagen, Hilden, Germany). The dsRNA was reverse transcribed to complementary DNA (cDNA) using random primers (Invitrogen, Paisley, UK) and reverse transcriptase enzyme (Superscript III MMLV-RT, Invitrogen, Paisley, UK) [9]. The cDNA was used to assign G genotype (G1, G2, G3, G4, G8, G9, G10, G11, G12) and P genotype (P[4], P[6], P[8], P[9], P[10], P[11], P[14]) using a multiplex heminested RT-PCR as described previously [5].
2.5 Statistical analysis
All statistical analyses were performed in the R environment for statistical computing, version 4.0.2 [10] and GraphPad Prism version 8. Vesikari score distributions were compared between pre-RV1 unvaccinated children, post RV1-unvaccinated children and post RV1-vaccinated children as well as between genotypes using non-parametric Kruskal-Wallis tests. Wilcoxon rank-sum tests were used to compare Vesikari scores between two groups. A linear regression model (Vesikari score = β0 + β1 * [post RV1-unvaccinated] + β2 * [post RV2-vaccinated] + ϵ) was used to estimate the change in Vesikari score that could be associated with receipt of RV1 and RV1 vaccine period. This model was adjusted for the Mid-Upper Arm Circumference (MUAC), age, weight, gender and receipt of the Bacillus Calmette–Guérin (BCG) and Pentavalent (diphtheria, pertussis, tetanus, and hepatitis B and Haemophilus influenzae type b) EPI vaccines. Both unadjusted and adjusted linear regression analysis were used to estimate Vesikari scores by considering G genotypes separately, P genotypes separately and combined G and P genotypes. For instance, the following unadjusted model was used to estimate effects of combined G & P genotypes: Vesikari = β0 + β1 * [post RV1-unvaccinated] + β2 * [post RV2-vaccinated] + β3 [genotype G2P[4]] +β4 [genotype G2P[6]] + β5 [genotype G12P[6]] +β6 [genotype G12P[8]] + β7 [post RV1-unvaccinated * genotype G2P[4]] + β8 [post RV1-unvaccinated * genotype G2P[6]] + β9 [post RV1-unvaccinated * genotype G12P[6]] + β10 [post RV1-unvaccinated * genotype G12P[8]] + β11 [post RV1-vaccinated * genotype G2P[4]] + β12 [post RV1-vaccinated * genotype G2P[6]] + β13 [post RV1-vaccinated * genotype G12P[6]] + β14 [post RV1-vaccinated * genotype G12P[8]] + ϵ. Residuals were computed for the purpose of model diagnostics: homoscedasticity and normality of residuals, linearity of the relationship between the independent and dependent variables.
2.6 Ethics
Ethical approval was obtained from the Research Ethics Committee of the University of Liverpool, Liverpool, UK (000490) and the National Health Sciences Research Committee, Lilongwe, Malawi (#867).
3. Results
2.1 Reduction in gastroenteritis severity following vaccination with RV1
The characteristics of the study participants included in this analysis are summarised in Table 1. In total, 3,159 children were recruited, of which 401 (12.7%) were recruited before RV1 introduction, whereas 2,758 (87.3%) were recruited after RV1 introduction. A total of 80.6% (2,224/2,758) of children recruited in the post-RV1 period were vaccinated. Thus, across the entire study period, 70.4% (2,224/3,159) of the children were vaccinated with RV1.
RV1-vaccinated children presented with less severe gastroenteritis compared with RV1-unvaccinated children during the post-vaccine period (unadjusted Kruskal-Wallis test, p <0.001). There was no difference in the severity of gastroenteritis between RV1-unvaccinated children before and after RV1 introduction (unadjusted Wilcoxon rank-sum test, p = 0.260) (Figure 1a). When RV1-vaccinated children were stratified into rotavirus-positive and rotavirus-negative cases, a decrease in severity score was observed in RV1-vaccinated children for both groups (Figure 1b and 1c). Reductions in gastroenteritis severity three years or later following RV1 introduction were observed in all age groups (Figure 2 and Table S1).
Unadjusted regression analysis confirmed a linear relationship between the reduction in gastroenteritis severity and vaccination with RV1 when pre-RV1 unvaccinated hospitalised gastroenteritis cases were used as a reference group. There was an average estimated reduction of 2.35 (95% confidence interval (CI) 2.03, 2.67; p<0.001) in severity scores among RV1-vaccinated children and no reduction among RV1-unvaccinated children during the post-vaccine period (0.22; 95% CI -0.16, 0.61; p=0.260). Adjusting for the MUAC, age, weight, gender, and EPI vaccination status in the linear regression did not substantially change the estimated reduction in severity scores among RV1-vaccinated (2.21; 95% CI: 1.85, 2.56; p<0.001) and RV1-unvaccinated children (0.05; 95% CI -0.46, 0.36; p=0.820] (Table S2). Unlike the other covariates, there was some evidence of a linear association between age and gastroenteritis severity when children recruited before RV1 introduction and without a RV1 vaccination history were used as a reference group: every increment of 1 year in the age was associated with a decrease of 0.43 (95% CI 0.26, 0.60; p<0.001) in Ruuska and Vesikari score (Table S2). Except for < 6-month-old rotavirus-positive children where severity started to decline from the 2013 – 2014 calendar year, at least a year post-RV1 introduction, substantial decline in severity was observed in all age groups between 2014 and 2015, at least three years after RV1 introduction. Severity scores increased in all cases from 2017 to 2019 (Figure 2).
3.2 Relationship between rotavirus genotype, gastroenteritis severity and vaccination with RV1
The most frequently detected rotavirus genotypes were G1P[8], G2P[4], G2P[6], G12P[6] and G12P[8] [5], which comprised 66.57% (636/1,050) of all genotypes (Table S3). In RV1 unvaccinated children, regardless of the genotype, most gastroenteritis episodes were classified as severe (66.67%, 424/636) and severity scores did not differ by genotype (unadjusted Kruskal-Wallis test, P=0.544). In RV1-vaccinated children, a decrease in severity was observed for infections with all rotavirus genotypes compared to RV1-unvaccinated children, with most vaccinated children having moderate disease. This decrease in severity was less pronounced in cases infected with G12P[6] and G12P[8] genotypes; and RV1-vaccinated children infected with G12P[6] strains had a 2.58 unit (95% CI 0.60, 4.56; p=0.011) higher severity score when compared with the remaining strains after adjusting for age (Figure 3 and Table S4). Stratifying severity scores by G and P genotypes separately suggested that this effect might be associated with the G12 genotype (Figure S2). There were no differences in severity scores between pre-RV1 and post-RV1 unvaccinated groups in G and P genotypes, but a clear difference within the post-RV1 vaccinated group (Figure S2). The decrease in severity among post-RV1 vaccinated children was more pronounced for genotypes G1 and G2 compared to G12 (unadjusted Kruskal-Wallis test, p <0.001). Adjusting for age did not affect the regression outcomes (data not shown).
5. Discussion
Introduction of RV1 into Malawi’s childhood immunisation schedule was associated with a significant reduction in the severity of all-cause gastroenteritis in children under the age of five years hospitalised at QECH in Blantyre. Irrespective of the rotavirus genotype, RV1-vaccinated children presented with less severe rotavirus disease compared to RV1-unvaccinated children during the pre- and post-RV1 introduction periods.
RV1 was developed to prevent children from developing severe gastroenteritis following infection with rotaviruses post-vaccination and not necessarily to prevent them from getting infected with rotaviruses [11, 12]. Thus, our study demonstrates expected, but to date unmeasured, direct impact of RV1 in reducing the severity of disease in children with acute rotavirus gastroenteritis in a low-income setting. The reduction in the severity of gastroenteritis in older children observed in the current study in both RV1-vaccinated and RV1-unvaccinated children could be explained by the potential acquired rotavirus-specific immunity either through rotavirus vaccination [6, 7, 13] or constant exposure to natural rotavirus infection due to the prevailing high force of rotavirus infection in Malawi [14], as the children grow older.
We also identified unexpected off-target vaccine benefit as RV1 recipients presented with less severe gastroenteritis regardless of the presence or absence of rotavirus. This is consistent with previous findings in Malawi and other settings where introduction of RV1 was associated with a reduction in all-cause diarrhoea mortality [7, 15, 16]. In contrast, a study that assessed the performance of rotavirus and oral polio vaccines in Bangladeshi infants did not find any evidence of nonspecific effect of RV1 on diarrheal burden or aetiology for a short period of 18 – 39 weeks post vaccination [17]. In a separate study in the same Bangladeshi population, delayed dosing of RV1 was associated with reduction in the duration of all-cause diarrhoea; however, rotavirus gastroenteritis was the only pathogen for which disease presentation was significantly reduced after a one-year follow-up [18].
Nonspecific vaccine effects were previously observed in infants vaccinated with BCG [19, 20] followed by herpes simplex [21, 22], measles [23] and other live attenuated vaccines such as the trivalent oral polio vaccine [17, 24]. T cell mediated cross-reactivity and trained innate immunity are among the mechanisms that may explain the off-target vaccine benefit [25, 26]. Vaccination with replication-deficient mutant herpes simplex type 1 virus provides protection from bacterial pathogens through induction of more robust CD8+ T cell cytotoxicity, prolonged production of interferon-gamma (IFN-γ), and systemic activation of macrophages [21, 22], whereas BCG elevates innate immune markers, such as IFN-γ, tumor necrosis factor α, interleukin 1βeta (IL-1β) and IL-6 cytokines [27, 28]. Oral live attenuated rotavirus vaccines could potentially employ similar mechanisms as they have been shown to induce innate immune responses [29] and effectively replicate in the gut of vaccinated children, [30-32] which could trigger cross-reactive CD8+ T cell responses. In addition, it is possible that averting rotavirus diarrhoea or experiencing less severe gastroenteritis in RV1-vaccinated children, [13, 33], promotes a healthier gut compared to their RV1-unvaccinated counterparts. The impaired gut integrity and physiology that follows severe rotavirus infection [34], may render them more susceptible to infection and/or more severe disease caused by other enteropathogens. This may potentially contribute, in part, to the indirect effects of rotavirus vaccine previously reported as herd protection [15, 16, 35-37].
RV1-vaccination was associated with a reduction in rotavirus gastroenteritis severity regardless of the infecting rotavirus genotype. Although RV1 effectiveness was previously documented to be lower against genotype G2-associated gastroenteritis compared with disease associated with genotype G1 [5, 38, 39], the reduction in severity in cases associated with G1 and G2 genotypes was similar among RV1-vaccinated children in the current study. However, the decline in gastroenteritis severity was less pronounced among cases associated with the G12 genotype. While the reasons for this observation are unknown, a relatively lower RV1 effectiveness has been reported against some heterotypic rotavirus strains, including G12s in this population [38]. Intriguingly, form 2017, rotavirus and non-rotavirus gastroenteritis severity scores increased in all age groups, although they remained lower than among those in unvaccinated children (Figure S1). Such a trend may be explained by potential changes in the criteria for hospitalisation of children presenting with gastroenteritis at QECH or changes in health seeking behaviour and or access of the communities around Blantyre post 2017. Further investigations are warranted to understand the significance and causes of this trend.
We could not test the trend in reduction of severity between vaccinated and unvaccinated infants because the number of unvaccinated infants decreased during each consecutive year due to increase in RV1 coverage. In addition, our analysis was limited by the availability of gastroenteritis cases from only one surveillance year prior to RV1 introduction, A further limitation is that all children were hospitalised, and hence had at least moderately severe gastroenteritis. Future studies should also examine the benefits of rotavirus vaccination in children with less severe gastroenteritis treated as outpatients or at surrounding health care facilities. Finally, the observed non-specific reduction in all-cause gastroenteritis could not be attributed to individual diarrhoea pathogens as the stool specimens were not routinely screened for pathogens beyond rotavirus.
In conclusion, our study provides evidence of a reduction in gastroenteritis severity among hospitalised children who had been vaccinated with RV1 seven years following its introduction in Malawi’s immunization programme. Rotavirus vaccination reduced the severity of rotavirus gastroenteritis caused by both homotypic and heterotypic rotavirus strains. Furthermore, rotavirus vaccination decreased non-rotavirus gastroenteritis severity, suggesting important off-target vaccine effects. Overall, these data demonstrate previously unmeasured direct and indirect benefits of rotavirus vaccines in Malawian children, providing further support for their continued programmatic use.
Data Availability
The data referred to in this manuscript is freely available upon request.
Supplementary Materials
Supplementary Data S1, Table S1: Gastroenteritis severity based on Vesikari scores among children presenting with gastroenteritis at the Queen Elizabeth Central Hospital in Blantyre, Malawi between the pre-(December 2011 – October 2012) and post-vaccine (November 2012 – October 2019) introduction periods.
Supplementary Data S2, Table S2. Unadjusted and adjusted linear regression model for the estimated reduction in severity scores among hospitalised children at Queen Elizabeth Central Hospital in Blantyre, Malawi after RV1 introduction.
Supplementary Data S3, Table S3. Frequently detected rotavirus strains and their gastroenteritis severity scores in children presenting with gastroenteritis at Queen Elizabeth Central Hospital, Blantyre, Malawi during the pre- and post-vaccine period (December 2011 to October 2019).
Supplementary Data S4, Table S4. Linear regression model for estimated reduction in severity scores and genotype (G and P), adjusted for age, MUAC and EPI vaccination status. Participants from the pre-RV1 period, with genotype G1P[8] and no vaccinations are taken as the reference group.
Supplementary Data 5, Figure S1. Trends in gastroenteritis severity among all children who presented with gastroenteritis at Queen Elizabeth Central Hospital in Blantyre, Malawi during the pre-(December 2011 to October 2012) and post-RV1 (November 2012 to October 2019) periods amongst different age groups. (A) All-cause gastroenteritis. (B) Rotavirus-positive gastroenteritis. (C) Rotavirus-negative gastroenteritis.
Supplementary Data 6, Figure S2. Severity scores in hospitalised children with rotavirus confirmed gastroenteritis at Queen Elizabeth Central Hospital, Blantyre, Malawi (A) Vesikari scores in RV1-vaccinated and RV1-unvaccinated children when only rotavirus P genotypes were considered. (B) Vesikari scores in RV1-vaccinated and RV1-unvaccinated children when only rotavirus P genotypes were considered.
Funding
This work was supported by research grants from the Wellcome Trust (Programme grant number 091909/Z/10/Z, Bill and Melinda Gates Foundation (Grant number: OPP1180423), and CDC funds through WHO (grant number: 2018/815189-0). MIG is partly supported by the NIHR HPRU in Gastrointestinal Infections. K.C.J. is a Wellcome International Training Fellow (201945/Z/16/Z). DH is funded by a National Institute for Health Research (NIHR) Post-doctoral Fellowship (PDF-2018-11-ST2-006).
Disclaimer
The funders had no role in the study design, data collection and interpretation, or the decision to submit the work for publication. The authors did not receive any financial support or other form of reward related to the development of the manuscript. Therefore, findings and conclusions in this report are those of the authors and do not necessarily represent the formal position of the funders. D.H, M.IG., NAC., and K.C.J are affiliated to the National Institute for Health Research (NIHR) Health Protection Research Unit in Gastrointestinal Infections at University of Liverpool, a partnership with Public Health England, in collaboration with University of Warwick. The views expressed are those of the author(s) and not necessarily those of the NIHR, the Department of Health and Social Care or Public Health England.
Author Contributions
M.I.G, N.A.C. and K.C.J conceived and designed the study. K.C.J. collected clinical data and stool samples. C.M., J.J.M., E.C., O.K., K.C.J. performed the laboratory work. J.J.M. M.Y.R.H., D.H., M.I.G., and K.C.J carried out the statistical analysis. J.J.M. and K.C.J. drafted the manuscript. J.J.M., M.Y.R.H., C.M., R.W., A.W.K., C.M.B., E.C, O. K., I.T.S., D.H., M.I.G., N.A.C. and K.C.J contributed to interpretation of the data and writing the manuscript. All authors have read and approved the final manuscript.
Potential conflicts of interest
K. C. J. has received investigator-initiated research grant support from the GSK group of companies. D. H and M. I.-G. has received investigator-initiated research grant support from the GSK group of companies and Sanofi Pasteur Merck Sharpe & Dohme and Merck. NAC has received investigator-initiated grant support for rotavirus research and honoraria for participation in DSMB rotavirus vaccine meetings from the GSK group of companies.
Acknowledgements
We acknowledge the support of the laboratory staff at the Malawi-Liverpool-Wellcome Trust Clinical Research Programme, clinical research team and the study participants. We thank Dr Naor Bar-Zeev, Dr Aisleen Bennett and Dr Louisa Pollock for their role in the collection of clinical data and faecal samples.