Despite the success of rotavirus vaccines over the last decade, rotavirus remains a leading cause of severe diarrheal disease among young children. Further progress in reducing the burden of disease is inhibited, in part, by vaccine underperformance in certain settings. Early trials suggested that oral poliovirus vaccine (OPV), when administered concomitantly with rotavirus vaccine, reduces rotavirus seroconversion rates after the first rotavirus dose with modest or nonsignificant interference after completion of the full rotavirus vaccine course. Our study aimed to identify a range of individual-level characteristics, including concomitant receipt of OPV, that affect rotavirus vaccine immunogenicity in high- and low-child-mortality settings, controlling for individual- and country-level factors. Our central hypothesis was that OPV administered concomitantly with rotavirus vaccine reduced rotavirus vaccine immunogenicity.
Methods and findings
Pooled, individual-level data from GlaxoSmithKline’s Phase II and III clinical trials of the monovalent rotavirus vaccine (RV1), Rotarix, were analyzed, including 7,280 vaccinated infants (5–17 weeks of age at first vaccine dose) from 22 trials and 33 countries/territories (5 countries/territories with high, 13 with moderately low, and 15 with very low child mortality). Two standard markers for immune response were examined including antirotavirus immunoglobulin A (IgA) seroconversion (defined as the appearance of serum antirotavirus IgA antibodies in subjects initially seronegative) and serum antirotavirus IgA titer, both collected approximately 4–12 weeks after administration of the last rotavirus vaccine dose. Mixed-effect logistic regression and mixed-effect linear regression of log-transformed data were used to identify individual- and country-level predictors of seroconversion (dichotomous) and antibody titer (continuous), respectively. Infants in high-child-mortality settings had lower odds of seroconverting compared with infants in low-child-mortality settings (odds ratio [OR] = 0.48, 95% confidence interval [CI] 0.43–0.53, p < 0.001). Similarly, among those who seroconverted, infants in high-child-mortality settings had lower IgA titers compared with infants in low-child-mortality settings (mean difference [β] = 0.83, 95% CI 0.77–0.90, p < 0.001). Infants who received OPV concomitantly with both their first and their second doses of rotavirus vaccine had 0.63 times the odds of seroconverting (OR = 0.63, 95% CI 0.47–0.84, p = 0.002) compared with infants who received OPV but not concomitantly with either dose. In contrast, among infants who seroconverted, OPV concomitantly administered with both the first and second rotavirus vaccine doses was found to be positively associated with antirotavirus IgA titer (β = 1.28, 95% CI 1.07–1.53, p = 0.009). Our findings may have some limitations in terms of generalizability to routine use of rotavirus vaccine because the analysis was limited to healthy infants receiving RV1 in clinical trial settings.
Our findings suggest that OPV given concomitantly with RV1 was a substantial contributor to reduced antirotavirus IgA seroconversion, and this interference was apparent after the second vaccine dose of RV1, as with the original clinical trials that our reanalysis is based on. However, our findings do suggest that the forthcoming withdrawal of OPV from the infant immunization schedule globally has the potential to improve RV1 performance.
Why was this study done?
- Rotavirus vaccination is the primary method for preventing rotavirus infections and reducing the severity of rotavirus gastroenteritis among young children globally.
- Rotavirus vaccine underperformance (lower vaccine efficacy/effectiveness and immunogenicity) among children in high-child-mortality settings where the incidence of severe disease and death is highest hinders further progress in reducing the rotavirus burden.
- Infant characteristics, such as malnutrition, age at first rotavirus dose, and concomitant receipt of oral poliovirus vaccine (OPV), may impact rotavirus vaccine immunogenicity.
- Identifying individual-level factors that impact vaccine immunogenicity is critical for determining modifiable vaccine strategies or interventions for enhancing vaccine performance.
What did the researchers do and find?
- We analyzed pooled, individual-level data on 7,280 vaccinated infants from 22 GlaxoSmithKline Phase II and III clinical trials of the monovalent rotavirus vaccine (RV1) across 33 high- and low-child-mortality countries/territories.
- We examined the relationship between individual- and country-level predictors and immune response to RV1 vaccination measured via serum antirotavirus immunoglobulin A (IgA) seroconversion and serum antirotavirus IgA titer.
- Our findings suggest that OPV given concomitantly with the first and second doses of rotavirus vaccine reduces antirotavirus IgA seroconversion.
- We did not find the same modifiable characteristics associated with seroconversion to be associated with postvaccine antirotavirus IgA titers among infants who seroconverted.
What do these findings mean?
- Factors that may predict whether an infant responds to RV1 (seroconversion) may not also predict the intensity of the response (titer) given seroconversion.
- Our findings highlight the importance of concomitant OPV and RV1 administration and provide encouraging evidence to suggest that the eventual withdrawal of OPV from national immunization programs could improve rotavirus vaccine performance.
Citation: Baker JM, Tate JE, Leon J, Haber MJ, Lopman BA (2019) Antirotavirus IgA seroconversion rates in children who receive concomitant oral poliovirus vaccine: A secondary, pooled analysis of Phase II and III trial data from 33 countries. PLoS Med 16(12):
Academic Editor: Julie Elissa Bines, The University of Melbourne School of Historical and Philosophical Studies, AUSTRALIA
Received: July 26, 2019; Accepted: November 22, 2019; Published: December 30, 2019
This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Data Availability: The data underlying the results presented in the study were provided by GlaxoSmithKline. Data are available upon request and approval from GlaxoSmithKline at https://www.clinicalstudydatarequest.com.
Funding: JMB received financial support to conduct this research from the Laney Graduate School of Emory University (http://www.graduateschool.emory.edu). BAL’s work was supported by NIH grants R01AI112970 and R01GM124280 (https://grants.nih.gov/grants/oer.htm) and the Modeling the Impact of Rotavirus Vaccination in GAVI Countries project supported by Imperial College London (https://www.imperial.ac.uk). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: I have read the journal’s policy and the authors of this manuscript have the following competing interests: BAL has received personal fees from Takeda Pharmaceuticals, personal fees from CDC Foundation, and personal fees from Hall Booth Smith, P.C., outside the submitted work. All other authors declare no competing interests.
Akaike information criterion; CCID50,
50% cell culture infective dose; CDC,
Centers for Disease Control and Prevention; CI,
confidence interval; ELISA,
enzyme-linked immunosorbent assay; GDP,
gross domestic product; GSK,
immunoglobulin A; IPV,
inactivated poliovirus vaccine; IQR,
interquartile range; LAZ,
length-for-age z-score; OPV,
oral poliovirus vaccine; OR,
odds ratio; ref,
reference group; RV1,
monovalent rotavirus vaccine; SD,
standard deviation; STROBE,
Strengthening the Reporting of Observational studies in Epidemiology; VIF,
variance inflation factor; WHO,
World Health Organization
Globally, rotavirus is the leading cause of severe diarrheal disease among infants and children under 5 years of age, estimated to cause 128,500–215,000 deaths in this age group each year during 2013–2016, when global vaccination coverage was <25% [1,2]. The virus is highly infectious, and improvements in water, sanitation, and hygiene conditions have limited impact in reducing its spread . Rotavirus is a ubiquitous infection among young children , the majority of whom will experience at least one rotavirus infection in the first 2 years of life [4–6]. As such, vaccination is an essential public health measure for preventing infections and reducing the severity of rotavirus gastroenteritis .
Currently, four live attenuated, oral rotavirus vaccines administered during infancy have received World Health Organization (WHO) prequalification: GlaxoSmithKline’s (GSK) monovalent rotavirus vaccine (Rotarix, “RV1”), Merck’s pentavalent vaccine (RotaTeq), the Serum Institute of India’s pentavalent vaccine (Rotasiil), and Bharat Biotech’s monovalent vaccine (Rotavac) . Since the pivotal clinical trial results were released in 2006 [9–11], rotavirus vaccines have been integrated into the national immunization programs of approximately 96 countries . The introduction of these vaccines has led to dramatic reductions in rotavirus disease in many settings [9,10,13–15]. Despite this success, rotavirus remains a leading cause of severe diarrheal disease among young children  and continues to be the predominant cause of hospitalization for severe diarrheal disease, even in some countries where the vaccine is in widespread use [16,17].
A primary obstacle preventing further reductions in the rotavirus burden is vaccine underperformance among children in settings where the incidence of severe disease and death is highest [2,18]. An estimated 85% of rotavirus-related deaths occur among children in Africa and Asia , and in these high-child-mortality settings, vaccine effectiveness against severe gastroenteritis is approximately 50%–60% . In contrast, vaccine efficacy is greater than 90% in low-child-mortality settings [11,13,20]. Rotavirus vaccines are substantially less immunogenic in high-child-mortality settings when compared with low-child-mortality settings [13,20]. Rotavirus vaccine immunogenicity, represented by serum antirotavirus immunoglobulin A (IgA) antibodies, has been shown to vary by country and is inversely associated with mortality rate in children under 5 years of age at the country level .
Reasons for this disparate immunogenicity and efficacy/effectiveness are poorly understood [22,23]. Individual-level host characteristics may play a central role in vaccine immunogenicity by reducing the effective vaccine virus titer (plaque-forming units) delivered to the intestines or reducing the immune response to vaccination [22,24]. For example, immunoglobulin G antibodies received from breast milk might decrease the effective vaccine virus titer delivered to the gut, whereas malnutrition may impair an infant’s immune response to the vaccine . Additional factors that may compromise vaccine performance relate to genetic susceptibility to rotavirus [25–28], environmental enteropathy (chronic inflammation) , and differences in force of infection [13,22,30].
One potentially modifiable factor relates to the interaction with live oral poliovirus vaccine (OPV). Poliovirus vaccines are generally administered on the same infant immunization schedule as rotavirus vaccines. In low-child-mortality settings, inactivated poliovirus vaccine (IPV) is primarily administered, whereas in high-child-mortality settings, live OPV is used . OPV administered concomitantly with rotavirus vaccine has been shown to reduce rotavirus seroconversion [32–36]. Early clinical trials suggested that this inhibitory effect may be strongest with the first rotavirus dose, with more modest, statistically nonsignificant differences observed after completion of the full rotavirus vaccine course . However, more recent evidence suggests that interference from OPV—be it monovalent, bivalent, or trivalent—may persist after two doses of rotavirus vaccine [32,35]. If OPV does substantially interfere with RV1, the global strategy to withdraw trivalent OPV, as implemented in 2016, and then OPV entirely  has the potential to improve rotavirus vaccine performance.
Developing strategies to improve vaccine performance first requires identification of individual-level factors associated with immune response. Isolating such factors across settings may highlight potentially modifiable vaccine strategies or interventions for enhancing vaccine performance and further reducing the burden of rotavirus disease. Rotavirus vaccine clinical trials were powered to assess vaccine efficacy; however, they were not specifically designed to identify individual-level factors associated with vaccine response. In this study, we pooled data from 22 clinical trials in 33 countries (5 countries/territories with high, 13 with moderately low, and 15 with very low child mortality) to identify a range of individual-level characteristics that contribute to RV1 immunogenicity measured via serum antirotavirus IgA in high- and low-child-mortality settings controlling for individual- and country-level factors. A central hypothesis was that OPV administered concomitantly with rotavirus vaccine doses reduced rotavirus vaccine immunogenicity.
GSK clinical trial data
We pooled individual-level data from infants (5–17 weeks of age) enrolled in GSK’s Phase II and III clinical trials of the RV1 vaccine initiated during 2000–2010 (Table 1). RV1 is an oral two-dose vaccine based on a live, attenuated human rotavirus strain (G1P). GSK recommends that the first dose be administered beginning at 6 weeks of age and the second dose be given after an interval of 4 or more weeks and by 24 weeks of age .
The 22 trials included were randomized, double-blind, placebo-controlled trials conducted in a total of 33 countries/territories, including 5 with high, 13 with moderately low, and 15 with very low child mortality (Fig 1). The child mortality strata categories were based on 2004 under-five mortality rates as previously described (S1 Table) [40,41]. Settings with moderately low and very low child mortality were combined and categorized as “low” child mortality for the primary analyses. Clinical research protocols across trials were similar in terms of data collection techniques, vaccine administration, and immunogenicity outcome measures (Table 2). Since the primary aim of the analysis was to examine factors associated with rotavirus vaccine immunogenicity, data were limited to trial participants who received RV1 (n = 53,292, placebo groups excluded). Data were further restricted to infants whose trial participation was completed according to protocol (classified by GSK) and who participated in the rotavirus immunogenicity substudies of the trials (n = 8,309). Lastly, infants who had serum sample collection approximately 4–12 weeks from receipt of his/her last rotavirus vaccine dose were included (n = 7,298).
Fig 1. Clinical trial sites included in the analysis.
Color shading represents the number of participants from each country. Country outline indicates a country’s child mortality classification (see S1 Table for details). For the primary analysis, countries categorized as “very low” and “moderately low” child mortality were combined and categorized as “low” child mortality. Map created using the “rworldmap” package in R.
Explanatory variables, covariates, and end points of interest
All available data for trial participants were provided by GSK, and explanatory variables/covariates were selected for inclusion in the analysis based on existing literature [22,24]. GSK data were supplemented with country-level data on gross domestic product (GDP) per capita in 2004 USD [42,43] to represent a country’s level of development and 2004 under-five mortality rates [44–46]. These country-level variables were considered in an effort to capture potential confounders that remained unmeasured despite the individual-level factors available.
Two standard markers for immune response [47,48] were analyzed as outcomes. First, in accordance with prespecified trial definitions [49,50], seroconversion was defined as the appearance of serum antirotavirus IgA antibodies (i.e., concentrations ≥20 U/mL) in subjects initially (prior to the first rotavirus dose) seronegative. The second end point of interest was serum antirotavirus IgA titer among infants who seroconverted. In all trials, postvaccine antirotavirus titer IgA data were collected approximately 4–12 weeks after the last administration of rotavirus vaccine and were measured using enzyme-linked immunosorbent assay (ELISA) techniques .
Regression models were fit to estimate the effects of individual- and country-level factors on vaccine immunogenicity outcomes while controlling for potential confounders. Mixed-effect logistic regression and mixed-effect linear regression models were used to analyze the antirotavirus IgA seroconversion (dichotomous) and antirotavirus IgA antibody titer (continuous, log transformed) outcomes, respectively. Basic formulas for each model and our detailed analysis plan are shown in S1 Text. Models were fit using the “lme4” package in R software.
The data used for modeling antirotavirus IgA seroconversion included all infants who were either confirmed to be seronegative prior to vaccination via serology sample or who did not have prevaccine serology data available. Prior confirmed rotavirus gastroenteritis was an exclusion criterion in a majority of the included trials, so for the primary analysis in this study, infants without prevaccine serology samples were assumed to be seronegative. The data for modeling antirotavirus IgA titer were restricted to infants who were seropositive after vaccination.
The modeling strategy for both outcomes began with variable specification, incorporating all individual-level characteristics as explanatory variables and controlling for potential country-level covariates selected based on variable distributions and existing literature. Individual-level variables in the initial model included time from last rotavirus dose to serology sample/outcome measure, number of rotavirus vaccine doses, age at first vaccine dose (weeks), vaccine concentration, sex, length-for-age z-score (LAZ) to represent nutritional status, and concomitant OPV. Country-level variables in the initial model included GDP, under-five mortality rate, and child mortality stratum (dichotomous). Because our analysis was restricted to infants who participated in the trials according to protocol, the only instances of missing data were related to LAZ; height was not included in the protocols of a few studies, and infants without these measurements were excluded in the models in which LAZ was used as a predictor (seroconversion models: n = 691 [9%]; IgA titer models: n = 525 [10%]). Random effects were employed using a random intercept for each trial to account for potentially unmeasured differences between trial protocols or environments. The initial model included all explanatory variables of interest and incorporated interaction terms between each main effect and child mortality stratum. The model was applied to combined data from all child mortality strata, and backward elimination was then conducted using an α = 0.10 cutoff for variable inclusion, maintaining a hierarchically well-formulated model throughout.
Next, models were refined by investigating alternative measures of the included variables based on possible relationships identified in bivariate analyses. For example, higher-order terms for continuous variables (such as age and age-squared) were considered when appropriate, and more-detailed stratification of categorical variables was tested. Backward elimination was subsequently performed after each variable was modified. When initial models were prohibitively large and issues with model convergence were encountered, the most informative end models resulting from previous backward elimination procedures were used as the starting point for investigation of refined measures. After all variables were explored, the most parsimonious model with all relevant variables and covariates was selected as the final model using Akaike information criterion (AIC) criteria (lower AIC indicating a better model) and assessed for multicollinearity using variance inflation factors (VIFs). Population attributable fraction and attributable fraction among the exposed were then calculated for variables of interest.
We examined whether the relationships identified using the full, combined dataset remained consistent within each mortality stratum; this would indicate that use of the combined dataset captured differences beyond those driven by child mortality stratum alone. To this end, a sensitivity analysis was conducted by fitting the final seroconversion model separately to data stratified by high- and low-child-mortality settings, as well as settings with moderately low and very low child mortality. In addition, the final seroconversion model was tested using a subset of the data made up of infants who were confirmed to be seronegative prior to vaccination. This subanalysis was conducted to confirm that results produced from analysis of the full seroconversion dataset were consistent with results using infants with both pre- and postvaccine serology. If consistent, this would suggest that inclusion of infants without prevaccine serology data was not causing substantial bias in our estimates. Lastly, a small number of infants received OPV with only their first or second rotavirus vaccine dose but not both. The final seroconversion model was rerun excluding these infants to simplify the OPV variable into three levels: (1) OPV not received concomitantly with the first and second rotavirus doses, (2) OPV received concomitantly with both rotavirus doses, and (3) no OPV received.
All identifiers for trial participants were deidentified and all dates obfuscated by GSK prior to data sharing. This study was determined not to meet the definition of research with human subjects by the Emory University Institutional Review Board. The Centers for Disease Control and Prevention (CDC) was nonengaged.
Data were available for analysis on 7,298 infants whose postvaccine rotavirus serology data were collected 4–12 weeks after receipt of the final vaccine dose (87.8% of the rotavirus immunogenicity cohort). Of these infants, 39.0% (n = 2,849) were from high-child-mortality settings. A total of 18 infants (3 from high-child-mortality settings) were excluded from analysis due to prevaccine serology data indicating these children had prior rotavirus infection (i.e., antirotavirus IgA titer of ≥20 U/mL prior to the first rotavirus vaccine dose). All other infants had either confirmed seronegative status (4,988, 68.5%) or did not have serology data available (2,292, 31.5%) and were assumed to be seronegative based on trial protocol. Data on 7,280 and 5,161 infants were included in the antirotavirus IgA seroconversion and antirotavirus IgA titer modeling, respectively. A basic comparison of the analysis cohort for this study to the entire vaccinated cohort from all 22 trials suggests the populations were demographically similar in terms of sex (analysis cohort versus entire vaccinated cohort: 49.7% versus 49.1% female), height at first dose (56.6 versus 57.4 cm), weight at first dose (5.2 versus 6.5 kilograms), and age at first vaccine dose (9.4 versus 8.8 weeks).
All infants received either two or three doses of rotavirus vaccine (Table 3); nearly all infants in low-child-mortality settings received two doses (99.1%), whereas over half of infants in high-child-mortality settings received a third dose because of study protocols (55.4%). In low-child-mortality settings, one-quarter of infants received a reduced concentration of the vaccine (viral suspension <106.0 50% cell culture infective dose [CCID50]), whereas all infants from high-child-mortality settings received a “standard” concentration (viral suspension ≥106.0 CCID50).
No association with seroconversion was identified for vaccine concentration, sex, age at first rotavirus vaccine dose, time from last rotavirus vaccine dose to postvaccine serology sample, and age at postvaccine serology sample in unadjusted analyses (Tables 4 and 5). In the unadjusted analysis, infants who received three rotavirus vaccine doses were less likely to seroconvert compared with those who received two doses (unadjusted odds ratio [OR] = 0.57, 95% confidence interval [CI] 0.50–0.64, p < 0.001), as nearly three-quarters (73.6%) of children who received three doses were from high-child-mortality settings where seroconversion was less common. Those who received OPV concomitantly with both their first and second doses of rotavirus vaccine had reduced odds of seroconversion compared with those who did not receive OPV concomitantly (unadjusted OR = 0.58, 95% CI 0.51–0.66, p < 0.001). A similar relationship was observed for infants who received OPV concomitantly with only their first or second rotavirus dose, though the sample size for these two categories was very small. Child mortality status (OR comparing high-child-mortality settings to low-child-mortality settings = 0.48, 95% CI 0.43–0.53, p < 0.001) and under-five mortality rate (OR for each one-unit increase in under-five mortality rate [log scale] = 0.69, 95% CI 0.66–0.72, p < 0.001) were negatively associated with seroconversion, whereas GDP was positively associated with seroconversion (OR for each one-unit increase in GDP [log scale] = 1.24, 95% CI 1.20–1.28, p < 0.001), as expected. The unadjusted associations between each variable and vaccine titer were generally similar to those observed for seroconversion.
Individual-level characteristics frequently differed by setting (Table 3). Over 20% of infants in high-child-mortality settings were stunted or severely stunted, whereas fewer than 10% were in low-child-mortality settings (p < 0.001). All infants in high-child-mortality settings received OPV, and over 90% of infants who received OPV received it concomitantly with both their first and second doses of rotavirus vaccine. In contrast, nearly half of infants in low-child-mortality settings did not receive any OPV (because IPV was more commonly administered in these settings). The percentage of children who seroconverted among the OPV categories ranged from 50% among those who received OPV concomitantly with only their first rotavirus vaccine dose to 83% among infants who did not receive OPV (Table 4). The median age at receipt of the first rotavirus dose was 9 weeks (interquartile range [IQR] = 7, 11) and infants were a median of 22 weeks of age (IQR = 21, 26) when their postvaccine serology sample was collected. The time from receipt of the last rotavirus dose to serology sample was slightly longer for infants in low-child-mortality settings (Table 3, median of 8 weeks, IQR = 5, 9) compared with infants in high-child-mortality settings (5 weeks, IQR = 5, 5, p < 0.001).
A majority of infants (70.9%) were seropositive after vaccination (Table 3), with a higher proportion of infants in low-child-mortality settings seroconverting (76.9%) compared with infants in high-child-mortality settings (61.5%, p < 0.001). Antirotavirus IgA seroconversion ranged from 58% in India to over 90% in Hong Kong, Italy, and Chile (Fig 2A). A similar pattern was found with postvaccine antirotavirus IgA titer, in which infants in low-child-mortality settings had geometric mean titers (240, SD = 4) higher than that of infants in high-child-mortality settings (199, SD = 4, p < 0.001). Antirotavirus IgA titer ranged from a median of 34 U/mL in India to 443 U/mL in Japan (Fig 2B).
Fig 2. Percent of infants who seroconverted by country and postvaccine antirotavirus IgA titer (U/mL) by country.
(A) Colored bars represent the percentage of infants in a given country that seroconverted, defined as a serum antirotavirus IgA titer ≥20 U/mL. Error bars indicate the 95% confidence interval. (B) Colored bars represent median and IQR for postvaccine antirotavirus IgA titer. Lines represent 25th percentile − 1.5*IQR and 75th percentile + 1.5*IQR, with dots representing outliers below or above these lines (Tukey box plot method). IgA, immunoglobulin A; IQR, interquartile range.
Antirotavirus IgA seroconversion modeling results
Results of backward elimination and model refinement for antirotavirus IgA seroconversion indicated several main effects and interaction terms contribute to the odds of seroconversion (Table 6). Being in a country with higher GDP and older age at first rotavirus dose were both associated with increased likelihood of antirotavirus IgA seroconversion. In contrast, increased time from last rotavirus dose to serology, low vaccine concentration, and concomitant receipt of OPV with the first and second rotavirus doses were each negatively associated with antirotavirus IgA seroconversion. Infants who received OPV concomitantly with both their first and second doses of rotavirus vaccine had 0.63 times the odds of seroconverting (OR = 0.63, 95% CI 0.47–0.84, p = 0.002) compared with infants who received OPV but not concomitantly with either dose. The final model included interaction terms with child mortality strata for both age at first vaccine dose and LAZ score. Sex, the number of rotavirus vaccine doses received, and under-five mortality rate were dropped from the model during backward elimination.
Among those who received OPV concomitantly with both the first and second rotavirus vaccine doses, 58.2% (95% CI 18.7%–100.0%) of vaccine nonresponse (i.e., those who did not seroconvert) was due to concomitant OPV administration (attributable fraction among the exposed). At the population level, the attributable fraction depends on the proportion of the population exposed; for example, if 95% of an infant population was administered OPV concomitantly with the first two doses of rotavirus vaccine (an estimate based on WHO’s recommended administration of both vaccines at 6 and 10 weeks of age), over half (53.7%, 95% CI 17.6%–99.7%) of vaccine nonresponses would be attributable to concomitant OPV administration with the first two rotavirus vaccine doses.
When the final antirotavirus IgA seroconversion model was applied to each mortality stratum separately (high, low, moderately low, and very low child mortality), similar results were observed (S2 Table). The negative relationship between OPV concomitant with both rotavirus doses 1 and 2 remained within each mortality setting, though with less precision of the parameter estimate. Similarly, when including only infants confirmed to be seronegative prior to rotavirus vaccination, results were consistent with the original analysis (S3 Table). Infants who received OPV concomitantly with rotavirus dose 1 and dose 2 were less likely to seroconvert than others (OR = 0.62, 95% CI 0.46–0.81, p = 0.001). The relationships between other variables and antirotavirus IgA seroconversion displayed only modest changes. Lastly, when the final model was fitted, excluding infants who received OPV concomitantly with only their first or second rotavirus vaccine doses, similar results were observed (S4 Table).
Antirotavirus IgA titer modeling results
Few factors were found to be significantly associated with antirotavirus IgA titers among those children who seroconverted (Table 7). Even when models were refined to explore alternative measures of various factors and when OPV was forced to remain in the model (based on its importance in seroconversion models), the results varied minimally. Increased time from last rotavirus dose to serology sample was negatively associated with antirotavirus IgA titer. Among infants who seroconverted, concomitant OPV with rotavirus doses 1 and 2 was positively associated with antirotavirus IgA titer (β = 1.28, 95% CI 1.07–1.53, p = 0.009), an opposite relationship to that observed among all infants in the seroconversion analysis.
We had the unique ability to pool individual-level data from 22 clinical trials and 33 countries/territories to create a dataset that enabled multilevel assessment of RV1 immunogenicity across settings. The resulting dataset was substantially larger than that of any individual trial or other related study . Our findings suggest that OPV given concomitantly with RV1 reduces antirotavirus IgA seroconversion even after two rotavirus vaccine doses. We did not find the same modifiable characteristics to be associated with postvaccine antirotavirus IgA titers among infants who seroconverted, suggesting that such factors may predict whether an infant responds to RV1 but not the intensity of the response given seroconversion.
These data provide robust evidence that infants who received OPV with both the first and second doses of RV1 were substantially less likely to seroconvert when compared with those not receiving OPV concomitantly. This analysis bolsters existing evidence [32–36] that OPV may interfere with rotavirus seroconversion and suggests that OPV may interfere with seroconversion when given with both the first and second rotavirus vaccine doses. Early evidence indicated that the OPV–rotavirus interaction may be strongest for the first rotavirus dose, with more recent evaluations suggesting lower seroconversion rates after completion of the full rotavirus vaccine course [32,33,35,37]. With this much larger and global dataset, our findings provide support for more recent data [32,35] that suggest a second RV1 dose may not compensate for the reduced initial response. Applying our final seroconversion model to each mortality stratum individually and conducting a separate sensitivity analysis restricting data to only infants with confirmed prevaccine seronegative status provided support, as the direction of the relationship remained within each stratum and despite dropping a substantial portion of the study data. This finding is consistent with the hypothesis the OPV interferes with rotavirus vaccine response, and this effect is not explained by socioeconomic differences among the trial sites. Although our study does not allow for assessment of possible mechanisms behind this relationship, it is consistent with hypotheses related to reduced efficiency of viral entry into mucosal cells, poliovirus causing down-regulation of components of the immune system response to rotavirus [32,52], or stimulation of innate immunity resulting in reduced rotavirus-specific immune response.
Overall, these results raise important programmatic considerations for rotavirus vaccination and evolving OPV eradication strategies. Although it may not be pragmatic to intentionally stagger OPV and rotavirus immunization schedules , the eventual shift from OPV to IPV may result in sizeable increases in RV1 performance in high-child-mortality settings. This notion is further supported by a prior study from Bangladesh that found that each formulation of OPV (monovalent, bivalent, and trivalent) reduced rotavirus seroconversion when administered concomitantly with rotavirus vaccine. Data on the specific formulation of OPV (trivalent versus bivalent) administered to infants in the trials in our study were not available to us, though it is likely that most, if not all, individuals who received OPV received the trivalent form based on the date of administration and national immunization schedules in these settings at the time of the trials.
GDP is known to be strongly associated with rotavirus seroconversion ; however, the effect of GDP is clearly not causal. Rather, GDP serves as a proxy for individual, family, and community factors that are more directly influential in determining vaccine immunogenicity, a portion of which we aimed to identify in this study. We found that GDP remained associated with the probability of seroconversion even after adjusting for child mortality and concomitant OPV vaccination. Comparing the adjusted OR for GDP and seroconversion developed via model selection (Table 6, OR = 1.11, 95% CI 1.04–1.19, p = 0.001) to their crude association (Table 5; OR = 1.24, 95% CI 1.20–1.28, p < 0.001) suggests that a substantial, though incomplete, portion of the crude association was accounted for with our individual-level factors.
We were unable to identify many modifiable factors associated with postvaccine antirotavirus IgA titers among infants who seroconverted, and the relationships between IgA titer and concomitant OPV were opposite to those expected; this suggests that intensity of a seropositive individual’s immune response to RV1 may be complex. Our unadjusted models demonstrated a negative relationship between OPV concomitantly administered with the first two doses of RV1 and antirotavirus IgA titer. However, given seroconversion and after adjusting for child mortality status, this relationship was no longer observed. These results are similar to those of Ramani and colleagues, who found that OPV coadministered with rotavirus vaccine reduced rotavirus seroconversion but did not reduce antirotavirus IgA titer among those who seroconverted . We observed a wide range in antirotavirus IgA titer, which may be due to a complex combination of individual-level factors not included in this analysis. As such, additional studies are needed to determine whether antirotavirus IgA titers among those who seroconverted may serve as a valuable correlate of protection for identifying the mechanisms resulting in differential vaccine immunogenicity and whether they are directly associated with vaccine effectiveness (i.e., the clinical relevance of IgA titers once seroconversion has been reached).
There are important challenges and limitations with this approach. First, although trial protocols were remarkably similar, trials themselves still differed by location, year, and population. We attempted to account for this variability by including a random effect for trial in our models. The required adherence to study protocol and stringent monitoring necessary for a clinical trial means that the results produced from these analyses may not perfectly reflect the findings that would have occurred under more routine, real-world conditions and should, therefore, be interpreted accordingly. Relatedly, the infants included in the trials were all healthy children, potentially limiting the generalizability of the results. Second, this is a secondary analysis of data previously collected for other primary purposes. We lacked data to control for genetic, maternal, socioeconomic, and environmental factors that likely influence individual-level immune response to vaccination. To mitigate residual confounding from factors such as socioeconomic status or environment, a proxy measure (GDP) was included in the models. Third, of the 7,280 infants included in the analysis, 2,292 (31.5%) did not have prevaccine serology data. We addressed this limitation by conducting a sensitivity analysis in which the antirotavirus seroconversion model was applied to only infants confirmed to be seronegative prior to vaccination, and we found the effect of OPV to be similar. Lastly, our analysis was limited to infants who received RV1, and it is possible that the findings and implications of this research may not be generalizable to the other three rotavirus vaccines prequalified by WHO or others currently in development.
Improving rotavirus vaccine performance requires identification of the factors that contribute to vaccine immunogenicity on the individual level. Although we explored a number of potential factors, our findings highlight the importance of concomitant OPV administration and provide encouraging evidence to suggest OPV withdrawal could improve RV1 performance. The ongoing effort by the Global Polio Eradication Initiative to end OPV use creates an ideal natural experiment to confirm our results in the real world. Vaccine immunogenicity data from infants in settings where OPV is currently in use could be compared with immunogenicity among infants after OPV withdrawal. More important still are evaluations of rotavirus vaccine effectiveness administered before and after OPV withdrawal against the clinical end point of rotavirus gastroenteritis. This research provides important programmatic considerations for improving rotavirus vaccine immunogenicity, particularly reduction in concomitant rotavirus vaccine and OPV administration.
S1 Text. Model equations.
Basic model formulas for mixed-effect logistic regression and mixed-effect linear regression of log-transformed data for antirotavirus IgA seroconversion (dichotomous) and antirotavirus IgA titer (continuous) outcomes, respectively. IgA, immunoglobulin A.
The authors would like to acknowledge www.clinicalstudydatarequest.com and GlaxoSmithKline for enabling access to the data.
The findings and conclusions of this report are those of the authors and do not necessarily represent the official position of the US Centers for Disease Control and Prevention.
Troeger C, Khalil IA, Rao PC, Cao S, Blacker BF, Ahmed T, et al. Rotavirus Vaccination and the Global Burden of Rotavirus Diarrhea Among Children Younger Than 5 Years. JAMA Pediatr. 2018;172: 958. pmid:30105384
Tate JE, Burton AH, Boschi-Pinto C, Parashar UD. Global, Regional, and National Estimates of Rotavirus Mortality in Children <5 Years of Age, 2000–2013. Clin Infect Dis. 2016;62: S96–S105. pmid:27059362
Centers for Disease Control and Prevention [Internet]. Atlanta, GA: Centers for Diesease Control and Prevention; [cited 2017 Aug 20]. Rotavirus and the Vaccine (Drops) to Prevent It. Available from: https://www.cdc.gov/vaccines/parents/diseases/child/rotavirus.html.
Velázquez FR, Matson DO, Calva JJ, Guerrero ML, Morrow AL, Carter-Campbell S, et al. Rotavirus Infection in Infants as Protection against Subsequent Infections. N Engl J Med. 1996;335: 1022–1028. pmid:8793926
Gurwith M, Wenman W, Hinde D, Feltham S, Greenberg H. A prospective study of rotavirus infection in infants and young children. J Infect Dis. 1981;144: 218–224. pmid:6268713
Gladstone BP, Ramani S, Mukhopadhya I, Muliyil J, Sarkar R, Rehman AM, et al. Protective Effect of Natural Rotavirus Infection in an Indian Birth Cohort. N Engl J Med. 2011;365: 337–346. pmid:21793745
Santosham M, Steele D. Rotavirus Vaccines—A New Hope. N Engl J Med. 2017;376: 1170–1172. pmid:28328339
World Health Organization [Internet]. Geneva, Switzerland: World Health Organization; 2016 Jun 3 [cited 2019 Feb 3]. Immunization, Vaccines and Biologicals: Rotavirus [Internet]. Available from: http://www.who.int/immunization/diseases/rotavirus/en/.
Vesikari T, Karvonen A, Prymula R, Schuster V, Tejedor J, Cohen R, et al. Efficacy of human rotavirus vaccine against rotavirus gastroenteritis during the first 2 years of life in European infants: randomised, double-blind controlled study. The Lancet. 2007;370: 1757–1763.
Vesikari T, Matson DO, Dennehy P, Van Damme P, Santosham M, Rodriguez Z, et al. Safety and Efficacy of a Pentavalent Human–Bovine (WC3) Reassortant Rotavirus Vaccine. N Engl J Med. 2006;354: 23–33. pmid:16394299
Nelson EAS, Glass RI. Rotavirus: realising the potential of a promising vaccine. The Lancet. 2010;376: 568–570.
Rota Council: Rotavirus Orgnaization of Technical Allies. Rotavirus Deaths & Rotavirus Vaccine Introduction Maps [Internet]. [cited 16 Oct 2016]. Available from: http://rotacouncil.org/toolkit/rotavirus-burden-vaccine-introduction-map/.
Jiang V, Jiang B, Tate J, Parashar UD, Patel MM. Performance of rotavirus vaccines in developed and developing countries. Hum Vaccin. 2010;6: 532–542. pmid:20622508
Linhares AC, Velázquez FR, Pérez-Schael I, Sáez-Llorens X, Abate H, Espinoza F, et al. Efficacy and safety of an oral live attenuated human rotavirus vaccine against rotavirus gastroenteritis during the first 2 years of life in Latin American infants: a randomised, double-blind, placebo-controlled phase III study. The Lancet. 2008;371: 1181–1189.
Ruiz-Palacios GM, Pérez-Schael I, Velázquez FR, Abate H, Breuer T, Clemens SC, et al. Safety and Efficacy of an Attenuated Vaccine against Severe Rotavirus Gastroenteritis. N Engl J Med. 2006;354.
Platts-Mills JA, Amour C, Gratz J, Nshama R, Walongo T, Mujaga B, et al. Impact of Rotavirus Vaccine Introduction and Postintroduction Etiology of Diarrhea Requiring Hospital Admission in Haydom, Tanzania, a Rural African Setting. Clin Infect Dis. 2017;65: 1144–1151. pmid:28575304
Operario DJ, Platts-Mills JA, Nadan S, Page N, Seheri M, Mphahlele J, et al. Etiology of Severe Acute Watery Diarrhea in Children in the Global Rotavirus Surveillance Network Using Quantitative Polymerase Chain Reaction. J Infect Dis. 2017;216: 220–227. pmid:28838152
World Health Organization. Rotavirus vaccines. WHO position paper–January 2013. Releve Epidemiol Hebd. 2013;88: 49–64.
Jonesteller CL, Burnett E, Yen C, Tate JE, Parashar UD. Effectiveness of Rotavirus Vaccination: A Systematic Review of the First Decade of Global Postlicensure Data, 2006–2016. Clin Infect Dis. 2017;65: 840–850. pmid:28444323
Soares-Weiser K, MacLehose H, Bergman H, Ben-Aharon I, Nagpal S, Goldberg E, et al. Vaccines for preventing rotavirus diarrhoea: vaccines in use. In: The Cochrane Collaboration, editor. Cochrane Database of Systematic Reviews. Chichester, UK: John Wiley & Sons; 2012. Available from: http://doi.wiley.com/10.1002/14651858.CD008521.pub3.
Patel M, Glass RI, Jiang B, Santosham M, Lopman B, Parashar U. A Systematic Review of Anti-Rotavirus Serum IgA Antibody Titer as a Potential Correlate of Rotavirus Vaccine Efficacy. J Infect Dis. 2013;208: 284–294. pmid:23596320
Patel M, Shane AL, Parashar UD, Jiang B, Gentsch JR, Glass RI. Oral Rotavirus Vaccines: How Well Will They Work Where They Are Needed Most? J Infect Dis. 2009;200: S39–S48. pmid:19817613
Parker EP, Ramani S, Lopman BA, Church JA, Iturriza-Gómara M, Prendergast AJ, et al. Causes of impaired oral vaccine efficacy in developing countries. Future Microbiol. 2018;13: 97–118. pmid:29218997
Velasquez DE, Parashar U, Jiang B. Decreased performance of live attenuated, oral rotavirus vaccines in low-income settings: causes and contributing factors. Expert Rev Vaccines. 2017; 1–17. pmid:29252042
Clarke E, Desselberger U. Correlates of protection against human rotavirus disease and the factors influencing protection in low-income settings. Mucosal Immunol. 2015;8: 1–17. pmid:25465100
Lee B, Dickson DM, deCamp AC, Ross Colgate E, Diehl SA, Uddin MI, et al. Histo–Blood Group Antigen Phenotype Determines Susceptibility to Genotype-Specific Rotavirus Infections and Impacts Measures of Rotavirus Vaccine Efficacy. J Infect Dis. 2018;217: 1399–1407. pmid:29390150
Payne DC, Currier RL, Staat MA, Sahni LC, Selvarangan R, Halasa NB, et al. Epidemiologic Association Between FUT2 Secretor Status and Severe Rotavirus Gastroenteritis in Children in the United States. JAMA Pediatr. 2015;169: 1040. pmid:26389824
Kambhampati A, Payne DC, Costantini V, Lopman BA. Host Genetic Susceptibility to Enteric Viruses: A Systematic Review and Metaanalysis. Clin Infect Dis. 2016;62: 11–18. pmid:26508510
Naylor C, Lu M, Haque R, Mondal D, Buonomo E, Nayak U, et al. Environmental Enteropathy, Oral Vaccine Failure and Growth Faltering in Infants in Bangladesh. EBioMedicine. 2015;2: 1759–1766. pmid:26870801
Bresee JS, Glass RI, Ivanoff B, Gentsch JR. Current status and future priorities for rotavirus vaccine development, evaluation and implementation in developing countries. Vaccine. 1999;17: 2207–2222. pmid:10403588
World Health Organization [Internet]. Geneva, Switzerland: World Health Organization; 2018 [cited 2019 Jan 30]. Summary of WHO Position Papers- Recommended Routine Immunizations for Children. Available from: https://www.who.int/immunization/policy/Immunization_routine_table2.pdf?ua=1.
Emperador DM, Velasquez DE, Estivariz CF, Lopman B, Jiang B, Parashar U, et al. Interference of Monovalent, Bivalent, and Trivalent Oral Poliovirus Vaccines on Monovalent Rotavirus Vaccine Immunogenicity in Rural Bangladesh. Clin Infect Dis. 2016;62: 150–156. pmid:26349548
Steele AD, De Vos B, Tumbo J, Reynders J, Scholtz F, Bos P, et al. Co-administration study in South African infants of a live-attenuated oral human rotavirus vaccine (RIX4414) and poliovirus vaccines. Vaccine. 2010;28: 6542–6548. pmid:18786585
Li R, Huang T, Li Y, Wang L-H, Tao J, Fu B, et al. Immunogenicity and reactogenicity of the human rotavirus vaccine, RIX4414 oral suspension, when co-administered with routine childhood vaccines in Chinese infants. Hum Vaccines Immunother. 2016;12: 785–793. pmid:27149266
Ramani S, Mamani N, Villena R, Bandyopadhyay AS, Gast C, Sato A, et al. Rotavirus Serum IgA Immune Response in Children Receiving Rotarix Coadministered With bOPV or IPV: Pediatr Infect Dis J. 2016;35: 1137–1139. pmid:27254033
Zaman K, Sack DA, Yunus M, Arifeen SE, Podder G, Azim T, et al. Successful co-administration of a human rotavirus and oral poliovirus vaccines in Bangladeshi infants in a 2-dose schedule at 12 and 16 weeks of age. Vaccine. 2009;27: 1333–1339. pmid:19162114
Patel M, Steele AD, Parashar UD. Influence of oral polio vaccines on performance of the monovalent and pentavalent rotavirus vaccines. Vaccine. 2012;30: A30–A35. pmid:22520134
World Health Organization. Polio vaccines: WHO position paper, March 2016–recommendations. Vaccine. 2017;35: 1197–1199. pmid:27894720
GlaxoSmithKline. Highlights of prescribing information [Internet]. Brentford, UK: GloxoSmithKline; 2016 [cited 2017 Mar 30]. Available from: https://www.gsksource.com/pharma/content/dam/GlaxoSmithKline/US/en/Prescribing_Information/Rotarix/pdf/ROTARIX-PI-PIL.PDF.
World Health Organization. List of Member States by WHO region and mortality stratum [Internet]. World Health Organization; [cited 2017 Aug 12]. Available from: https://www.who.int/whr/2004/annex/topic/en/annex_member_en.pdf.
UNICEF [Internet]. The UN Inter-agency Group for Child Mortality Estimation. New York: UNICEF; 2015 Sep 9 [cited 2019 Jan 18]. Available from: http://www.childmortality.org/.
World Bank Open Data [Internet]. Washington, DC: The World Bank; 2018 [cited 2018 Sep 28]. Available from: https://data.worldbank.org/.
Taiwan (Taiwan Province of China) GDP per Capita (Current Prices, US Dollars) Statistics [Internet]. Economy Watch; 2018 Nov 28 [cited 2018 Nov28]. Available from: http://www.economywatch.com/economic-statistics/Taiwan/GDP_Per_Capita_Current_Prices_US_Dollars/#otheryears.
Mortality rate, under-5 (per 1,000 live births) [Internet]. Washington, DC: The World Bank; 2018 [cited 2018 Sep 28]. Available from: https://data.worldbank.org/indicator/sh.dyn.mort.
Hong Kong- Under-five mortality rate [Internet]. Washington, DC: Knoema; [cited 28 Nov 2018]. Available from: https://knoema.com/atlas/Hong-Kong/topics/Demographics/Mortality/Under-5-mortality-rate.
Wu JC-L, Chiang T-L. Comparing Child Mortality in Taiwan and Selected Industrialized Countries. J Formos Med Assoc. 2007;106: 177–180. pmid:17339165
Zaccaro DJ, Wagener DK, Whisnant CC, Staats HF. Evaluation of vaccine-induced antibody responses: impact of new technologies. Vaccine. 2013;31: 2756–2761. pmid:23583812
Beyer WEP, Palache AM, Lüchters G, Nauta J, Osterhaus ADME. Seroprotection rate, mean fold increase, seroconversion rate: which parameter adequately expresses seroresponse to influenza vaccination? Virus Res. 2004;103: 125–132. pmid:15163500
GlaxoSmithKline. Safety Study Report Rota-023 (444563/023) [Internet]. GlaxoSmithKline; 2004 [cited 2019 May 17]. Available from: https://www.gsk-studyregister.com/study?uniqueStudyId=444563/023.
GlaxoSmithKline. Clinical Study Report for Study 102248 (Rota-037) [Internet]. GlaxoSmithKline; 2009 [cited 2019 May 17]. Available from: https://www.gsk-studyregister.com/study?uniqueStudyId=102248.
ThermoFisher Scientific [Internet]. Waltham, MA: Thermo Fisher Scientific; [cited 2019 Feb 14]. Overview of ELISA. Available from: https://www.thermofisher.com/us/en/home/life-science/protein-biology/protein-biology-learning-center/protein-biology-resource-library/pierce-protein-methods/overview-elisa.html.
Wang H, Moon S, Wang Y, Jiang B. Multiple virus infection alters rotavirus replication and expression of cytokines and Toll-like receptors in intestinal epithelial cells. Virus Res. 2012;167: 48–55. pmid:22497974