ABSTRACT
Streptococcus pneumoniae remains a leading cause of morbidity and mortality. Pneumococcal conjugate vaccines (PCVs) are effective but target only a fraction of the more than 90 pneumococcal serotypes. As a result, the introduction of PCVs has been followed by the emergence of non-vaccine serotypes. With higher-valency PCVs currently under development, there is a need to understand and predict patterns of serotype replacement to anticipate future changes. In this study, we evaluated patterns of change in serotype prevalence post-PCV introduction in Israel. We found that the assumption that non-vaccine serotypes increase by the same proportion overestimates changes in serotype prevalence in Jewish and Bedouin children. Furthermore, pre-vaccine prevalence was positively associated with increases in prevalence over the study period. From our analyses, serotypes 12F, 8, 16F, 33F, 9N, 7B, 10A, 22F, 24F, and 17F were estimated to have gained the most cases of invasive pneumococcal disease through serotype replacement in the Jewish population. However, this model also failed to quantify some additional cases gained, suggesting that changes in carriage in children alone may be insufficient to explain serotype replacement in disease. Understanding of serotype replacement is important as higher-valency vaccines are introduced.
Streptococcus pneumoniae (pneumococcus) is a leading cause of morbidity and mortality globally and causes a range of diseases including pneumonia, meningitis, and septicemia [1]. Nasopharyngeal carriage of pneumococcus is a precursor for disease [2]. Carriage is common in young children, who are the main reservoir of transmission to adults and other children [2, 3]. The likelihood of developing disease following exposure to pneumococcus depends on bacterial factors (i.e., serotype and other virulence factors) and host factors (age, immunocompetence, recent viral infection) [4, 5].
Pneumococcal conjugate vaccines (PCVs) protect against carriage and invasive infections due to serotypes targeted by the vaccine [6–9]. However, PCVs target only a fraction of the more than 90 pneumococcal serotypes. As a result, despite an overall decline in pneumococcal disease, the introduction of PCVs is followed by the emergence of non-vaccine serotypes among healthy carriers and, consequently, as causes of disease [6, 10, 11]. This ‘serotype replacement’ has reduced the overall impact of PCVs on disease rates and has prompted the development of PCVs that target larger numbers of serotypes. Several new PCVs are under development that target 15, 20, or more serotypes [12–14]. These vaccines will again likely disrupt the balance of serotypes in the nasopharynx and lead to further serotype replacement. As these next-generation conjugate vaccines move towards licensure, it is important to understand likely patterns of serotype replacement so that the marginal benefit of higher-valency vaccines can be anticipated. This is useful as countries decide which product to recommend in a given setting.
When evaluating the potential impact of new vaccines, there is a need to predict the overall benefits, which are influenced by both declines in disease incidence due to vaccine-targeted serotypes and by increases in disease incidence due to non-vaccine serotypes. At the same time, it is important to predict which serotypes are likely to emerge to help with future adjustments to the vaccine. To produce such projections, it is necessary to make assumptions about which serotypes are likely to increase and by how much. There are several proposed frameworks for accomplishing this goal. Nurhonen and Auranen predicted the impact of alternative PCV formulations with the assumption that all non-vaccine serotypes increase by the same proportion following vaccine introduction and that the invasiveness of the non-vaccine serotypes as a group is constant over time [15]. This approach can successfully capture the overall magnitude of the benefit of a vaccine by lumping vaccine and non-vaccine serotypes into broad groups [16–18]. However, when the goal is to examine emerging serotypes, it is important to consider individual serotypes. Recent work has used genomic data with models of negative frequency-dependent selection to quantify post-vaccine expansion patterns [19, 20]. However, whole genome sequence data are not always widely available, and their interpretation is complex. A complementary approach would be to identify characteristics of serotypes that emerge following vaccine introduction and to use this information to inform projections.
In this study, we characterized the expansion patterns of non-vaccine serotypes among healthy carriers following the introduction of PCVs in Israel. We evaluated correlates of serotype expansion, compared the expansion patterns in Israel to those that would be expected if the commonly-employed assumption that all non-vaccine serotypes expand by the same amount were correct, and assessed the implications for disease projections if this assumption is wrong. Finally, we provided a framework to rank serotypes not currently included in PCV13 based on the likelihood that they may be important causes of disease in the future.
METHODS
Data, setting, & participants
PCV7 was introduced for use in Israel in the national immunization program in July 2009 with a catch-up program in children <2 years old and was replaced by PCV13 in November 2010 without further catch-up. Carriage data were collected as part of prospective surveillance in southern Israel from November 2009-June 2016, described in detail elsewhere [21]. Every weekday during the study period, nasopharyngeal cultures were collected from the first four Jewish and first four Bedouin children under five presenting at a pediatric emergency department in Be’er Sheva, Israel for any reason. Information was recorded regarding ethnicity, result of the pneumococcal culture, pneumococcal serotype (if applicable), clinical diagnosis, number of PCV doses received to date (PCV7 and/or PCV13), and the year and month when the swab was recorded.
Invasive pneumococcal disease (IPD) data were also collected during the same time period as part of nationwide surveillance. These data were collected for individuals of all ages. Variables were collected regarding age group (<5, 5–17, 18–39, 40–64, and 65+ years), serotype, date (July 2009- June 2016), and number of swabs. Both the carriage and IPD datasets were aggregated into counts by serotype and epidemiological year (seven July-June years); the number of swabs negative for pneumococcus in the carriage dataset was also tallied. Year 0 in the carriage dataset is only ¾ of an epidemiological year, due to the fact that data collection started in November (instead of July) that year.
Other data
Additional serotype-specific characteristics were used in regression models to evaluate associations with baseline carriage prevalence and with changes in prevalence. These serotype-specific variables included relative density of growth reached at an early time point during in vitro growth [22], case fatality rate, negative charge carbon, total carbon per polysaccharide repeat, and capsule thickness [23]. In instances of missing data, missing values were imputed in the model assuming the data were missing completely at random (S1 Text).
Carriage model
The goal was to quantify changes in the prevalence of carriage of individual serotypes and to evaluate correlates of baseline prevalence or changes in prevalence. To do this, we fit a hierarchical Bayesian regression model to the data. The data for individual serotypes were sparse, so we used a hierarchical prior structure to provide stabilized estimates of the prevalence ratios for all serotypes [24]. In this model, pneumococcal carriage at time (epidemiological year) t follows a multinomial distribution with the number of successes being the number of detections of each non-vaccine serotype (i.e., each non-vaccine serotype representing a unique category), and the number of trials being the total number of nasopharyngeal swabs at time t. The “reference” category in this model includes all negative swabs and detections of PCV13 serotypes combined. The change in prevalence in the post-vaccine period varies by serotype and is a function of serotype-specific covariates. In this model, time is defined as a three-category variable, with the three categories corresponding to the pre- (November 2009-June 2010), early post-PCV (July 2010-June 2012), and late post-PCV (July 2012-June 2016) vaccination periods. This time structure provided a better fit compared to linear time, and was compared to several other functions of time. The three-category structure of time allows serotype prevalence to level out to its final post-vaccine prevalence, facilitating a better comparison of pre- to post-vaccine prevalence. The overall model is defined as:
n.carrt = pneumococcal carriage count at time period t
N.swabt = total # swabs at time period t
previit = pneumococcal prevalence for non-vaccine serotype i at time period t (m=54 non-vaccine serotypes)
previrt =pneumococcal prevalence from the “reference” group at time period t (“reference” swab is defined as negative for any pneumococcus or positive for vaccine-serotypes)
i = 1, 2,…, m
t= 0, 1,…, 6 (year)
time1= 0, 1, 1, 1, 1, 1, 1
time2= 0, 0, 0, 1, 1, 1, 1
xij= the jth serotype-level covariate for non-vaccine serotype i
All prior distributions were chosen to be weakly informative. The intercepts (relative prevalence at time 0) and slopes (early and late post-vaccine prevalence compared to time 0) were allowed to vary by serotype and were estimated hierarchically by having the serotype-specific parameters centered around global parameters that were common to all non-vaccine serotypes. Conjugate priors were used in all cases except for the standard deviation parameters, where uniform priors were used [25].
Different assumptions about model structure and covariates were evaluated for this model. Variations in random and fixed effects were considered, different structures of time were evaluated, and covariance matrices were applied to model both independence and dependence of the intercepts and slopes. The best model formulation was chosen using the optimal deviance information criterion (DIC) [26], using decreases of at least 10 to be considered an improved model. Using Markov chain Monte Carlo sampling techniques, we obtained 150,000 posterior samples from the joint posterior distribution following a burn-in period of 10,000 iterations. The Gelman-Rubin [27] and Geweke [28] diagnostics were used to investigate convergence of individual parameters and effective sample sizes were calculated to ensure we collected enough posterior samples post-convergence to make accurate inference. Additional details about the model, including structure and model diagnostics, can be found in the supplementary text (S1 Text).
Serotype-specific prevalence, relative risk ratios (RRRs; the multiplicative change for a one-year increase for a specific serotype prevalence relative to the reference prevalence), slopes, and intercepts were estimated from the final model. Where relevant, these estimates were compared to the corresponding values from the constant proportional change model, which assumes that serotype replacement in carriage is complete and that all serotypes increase by the same factor (i.e., prevalence ratio is the same for all serotypes).
Linking carriage with disease
We sought to quantify the impact of serotype replacement on IPD. This is challenging because many serotypes have secular trends or exhibit epidemic patterns even in the absence of the vaccines. Simply comparing the cases of IPD before and after vaccination can sometimes be misleading. Instead, we chose a more stable approach which used the estimates of the serotype-specific prevalence ratios from our carriage model to estimate the additional cases gained as a result of increased carriage in children (cases.gained.explained). The number of cases that would be expected in the absence of changes in carriage in children is calculated by dividing the observed number of cases at a given time point by the prevalence ratio for the same serotype. The cases gained (explained by increased carriage) estimate is calculated by subtracting this counterfactual for the number of cases from the observed number of cases for each serotype, age category and time point:
where the prevalence ratios (prev.ratio) from time 0 to time 6 were estimated from the model. We also calculated the number of cases gained that were not explained by increased carriage in children (e.g., could be due to secular trends; cases.gained.unexplained), by subtracting the observed number of IPD cases at time 0 from the number of cases that would have been expected in the absence of changes in carriage in children (i.e., the number of IPD cases in the first time period):
Using the values of cases gained as a result of increased carriage in children, we sorted the serotypes from highest to lowest. All of these calculations were carried out using results from the Jewish carriage model only, since the IPD dataset predominantly included Jewish individuals.
In some instances, the expected number of IPD cases exceeded the number of observed cases. In such cases, we adjusted the expected number of cases to be the number of observed cases. In other words, we let the maximum number of expected cases be the number of observed cases, resulting in zero cases gained in that instance.
To assess whether changes in carriage in younger children or older children better correlated with changes in IPD in adults, we additionally fit the carriage model separately for younger children (<12 months) and older children (different subsets of 24–60-month-old children), and compared the number of unexplained vs explained cases gained between the two age groups.
All analyses were performed using R version 3.4.0 [29]. The carriage model was fit using JAGS version 4.3.0 [30]. All code, with simulated data, can be found online at https://github.com/mailephillips/post-pcv-expansion.
RESULTS
Characteristics of the data
The Israel carriage dataset contained 10,396 observations, with 57.1% of samples from Bedouin children and the rest from Jewish children. Of all swabs, 47.7% were positive for pneumococcus, with 42.6% positive among Jewish children and 51.5% positive among Bedouin children. In 2009, just after vaccine introduction, 63.5% of colonized children had a vaccine-targeted serotype and 36.5% had a non-vaccine serotype. The most commonly-carried non-PCV13 serotypes in 2009 were 15B/C, 15A, 16F, 10A, 38, 35B, 23B, 23A, 21, and 10B.
Serotype expansion patterns in pneumococcal carriage
Some typical assumptions when projecting serotype replacement are that all non-vaccine serotypes increase by the same amount (i.e., same prevalence ratio), that the rank order of non-vaccine serotypes is maintained following vaccine introduction, and that non-vaccine serotypes completely replace vaccine serotypes in carriage. When we compare the observed pre-vaccine prevalence of serotypes to the post-vaccine prevalence, serotypes increase by different proportions and do not maintain rank order, suggesting that the assumption of constant prevalence ratio does not hold (Figure 1). However, our analysis suggests that instead some serotypes expand more than expected, and others remain stable or even decline (Figure 2).
The observed log-prevalence of pneumococcus for each serotype is shown at Year 0 and Year 6, separately for Jewish and Bedouin children. If a serotype increased from or decreased to 0, the log prevalence is shown as “not detected” at the bottom of the graph.
Log-transformed ratios comparing late post- to pre-vaccine prevalence (last year/first year) are shown for each serotype with 95% credible intervals from the Jewish and Bedouin carriage models. The serotypes are shown from highest to lowest log prevalence ratios.
The ten serotypes with the largest relative increase in carriage prevalence among Jewish children, from highest to lowest, were 7B, 17F, 9N, 10A, 23B, 21, 31, 35B, 20, and 8. Among Bedouin children, the largest relative increases were 9N, 7B, 10F, 35F, 10B, 35C, 2, 29, 13, and 31 (Figure 2, Table S2). The majority of the prevalence ratios were larger than one in both models (Jewish median prevalence ratio: 1.35, Bedouin median: 1.36). In both models, many serotypes had changes that were smaller than expected based on the model assuming proportional increases (constant proportional increase: 1.48 Jewish, 2.03 Bedouin). The serotype-specific intercepts and other serotype-specific parameters estimated from the model can be found in Table S2 and Figures S5–7.
Correlates of serotype-specific changes
While serotype-specific prevalence of the non-PCV13 serotypes mostly increased after vaccination, the degree to which each serotype changed varied. The distribution of the change in prevalence post-vaccination was skewed, but we were able to predict where on the distribution each serotype was located using serotype-specific characteristics. We evaluated several serotype-specific correlates of fitness to try to explain this pattern. The best correlate of the post-vaccine increase was pre-vaccine prevalence. The relationship between pre-vaccine prevalence and the change in serotype prevalence over time improved model fit when included with a covariance matrix compared to all other combinations of predictors. When fit to carriage data from Jewish children, the model with the covariance matrix was the best model, based on DIC. When fit to carriage data from Bedouin children, the model with independent slopes (no predictors) and the model with the covariance matrix both fit the data best based on DIC. There was an overall positive association between pre-vaccine and the change to late post-PCV prevalence in the Jewish dataset (0.04, 95% credible interval: –0.68, 0.68), and a negative association in the Bedouin dataset (−0.21, 95% credible interval: –0.69, 0.34). Combinations of relative density, case fatality rate, negative charge carbon, total carbon per polysaccharide repeat, and capsule thickness were not as strongly associated with the prevalence ratio and resulted in worse model fits in both datasets.
Ranking serotypes based on estimated cases gained
From 2009–2015, there were 3,968 total cases of IPD. We estimate that during this period, there was an excess of 624 (95% credible interval: 397, 850) cases gained across all age groups as a result of increased carriage of non-vaccine serotypes in children (Figure 3, Table 1). From highest to lowest, serotypes 12F, 8, 16F, 33F, 9N, 7B, 10A, 22F, 24F, and 17F were estimated to have gained the most cases of IPD through serotype replacement (each with an estimated 21 to 112 additional cases gained). These 10 serotypes combined had a total of 1,223 reported IPD cases during this period. An estimated 446 of these cases could be attributed to increases in carriage prevalence (Table 1). Fourteen serotypes did not account for any additional cases gained across this study period.
Results are shown for each age group and year of the study. Each full bar represents the total number of observed IPD cases for that stratum. The different patterned portions of the bars denote the fraction of total cases attributable to different groups of serotypes: observed vaccine-targeted (VT) serotype cases (solid gray), expected non-vaccine-targeted (NVT) serotype cases had PCV13 not been introduced (diagonal black lines), estimated additional cases gained as a result of increased carriage in children, explained by this model (horizontal black lines), and estimated additional cases gained that are not explained by this model (gray diagonal gray lines).
The results of the observed carriage at time 0, observed prevalence ratio (time6/time0), model prevalence ratio (time6/time0), and estimated cases gained through serotype replacement by age groups are shown for each of the serotypes used in the carriage model. Serotypes are ordered from most to least cases gained among all ages. Note that the total cases gained in the last row adds up to more than the number of cases above for each serotype, because the serotypes shown are only those used in the carriage model. The values for total cases gained also included calculations for additional serotypes (not shown) that were responsible for IPD, but were not present in carriage in children under 5 during this time period.
Serotype 12F, in particular, was responsible for 112 cases gained as a result of increased carriage in children out of the total 394 IPD cases for 12F. In all age groups, 12F was among the three highest ranked serotypes among cases gained as a result of increased carriage in children. Serotype 8 was similarly highly ranked for cases gained as a result of increased carriage in children in all age groups except for children <5 years. Serotype 16F was highly ranked among cases gained among adults 18 years and older, driving up its overall ranking.
Adults 65+ and children under five years were responsible for the most IPD cases (Table 1, Figure 3). In general, observed vaccine-type cases decreased over the study period, whereas the additional cases gained (both explained and unexplained) increased over the study period. Plots for individual serotype IPD cases over the study period can be found in Figure S8. Serotype 12F in particular was responsible for the most estimated cased gained that are not explained by increased carriage in children.
We had intended to compare results from carriage models with children less than one year and different subsets of 24–60 months of age to assess whether changes in carriage in younger children or older children better correlated with changes in IPD in adults; however, the data were too sparse for a formal comparison.
Cases not explained by carriage
The number of additional cases gained over this time period that were not explained by increased carriage in children was also substantial. Over the study period, there were an estimated 462 (95% credible interval: 312, 641) cases that were not predicted based on changes in carriage in children. The serotypes with the most unexplained cases gained were 12F, 24F, 22F, 10A, and 15BC (Table 2). Serotype 12F had the most unexplained cases, with an estimated 156 (95% credible interval: 50, 268) cases.
The median and 95% credible intervals for the unexplained excess cases during the full time period and across all ages is shown for each serotype.
DISCUSSION
With several higher-valency conjugate vaccines under development, there is a need to understand and predict the patterns of serotype replacement to anticipate future changes. In this study, we created a model that quantifies changes in serotype prevalence post-PCV introduction in Israel based on carriage data. We found that the common assumption that non-vaccine serotypes increase by the same proportion underestimates changes in serotype prevalence in Jewish children and overestimates them in Bedouin children. The change in prevalence post-vaccination was too variable across serotypes to use only one value to represent the proportional change for all serotypes. Instead, we found that the change in the prevalence of serotypes over time is not the same for all serotypes but has a positively skewed distribution. Using pre-vaccine prevalence data, we can estimate where on the distribution of post-vaccine changes that serotypes will fall. We can then combine observed IPD estimates and model-estimated prevalence ratios to get an estimate of serotype-specific replacement due to increased carriage. Using the estimated prevalence ratios, along with observed IPD data, we were able to rank serotypes based on cases gained as a result of increased carriage in children. However, there were also additional cases gained over this study period that were not predicted by increased carriage, suggesting the need to include additional mechanisms in the estimates.
The estimated number of cases that were not explained by the increased carriage in children illuminated some of the limitations of using carriage alone to predict disease. Changes in carriage prevalence explained some of the variations in IPD in adults, but not all. A large number of excess IPD cases occurred during this time period that were not predicted by increased carriage. Serotype 12F in particular had a large number of unexplained cases gained. Surprisingly, serotype 8 did not have many unexplained cases. In fact, we would have expected there to be more IPD cases due to serotype 8 overall, either explained or unexplained.
Changes like serotype-specific epidemics cycles or changes in ascertainment could account for some of this variability. Moreover, carriage patterns in specific age groups (e.g., older children, adults) might better reflect the changes in IPD seen in these other age groups. We had intended to carry out this comparison, but our data were too sparse in subsets of the population. Future work could investigate this difference further.
We found that pre-vaccine prevalence was correlated with where on the distribution of slopes serotypes are located, with serotypes having higher pre-vaccine prevalence associated with larger increases after vaccination. In other words, more prevalent serotypes before introduction of PCV-13 were associated with larger increases over the study period.
Our method of ranking serotypes varied between age groups, but overall the same serotypes had the highest number of cases gained through serotype replacement after vaccination. Unsurprisingly, serotype 8 was among the highest serotypes in this ranking system. It has recently been seen to be one of the most important serotypes to emerge in other populations, particularly in older age groups, which is consistent with our findings. Also not surprising were the high rankings of serotypes 12F and 33F, which have been noted to increase in recent years.
Recent work has shown that models that incorporate genomics or clonal groups can forecast which serotypes will be successful after vaccination [19, 20, 31]. Future work could combine the results of this study with new genomics studies, directly incorporating information on the fitness of a serotype (e.g., as calculated by Azarian) in these models, or this type of model structure could be used to benchmark alternative modeling approaches.
A fundamental limitation of using carriage data to forecast changes in IPD is that changes in serotypes that are rarely carried in children cannot be easily tracked using carriage data. It is not clear if serotypes with these characteristics, like serotypes 1, 5, and 8 might be transmitted directly among older age groups (bypassing children) or if they might be transmitted between symptomatic individuals.
This study had some other limitations. Many serotypes are rare, and as a result we did not have estimates for all 90+ serotypes for all variables. For the carriage data, we had count data for only 67 serotypes. However, this lack of data for the other serotypes is not likely to impact the results greatly because they are so rare in both carriage and disease and would probably not factor into the model.
This analysis was limited to the available data in a single country. The dataset contained only seven years of data, and only included carriage data from children. The model structure additionally had some assumptions. We attempted to address these structural assumptions by comparing several models (variations in random and fixed effects, covariance matrices, predictors, and time). When compared to the observed data, the model appeared to fit, suggesting that these assumptions were valid. Future research could investigate these assumptions further, using other populations and additional time periods. Additionally, the IPD data available were comprised of mostly Jewish individuals. As a result, we were unable to extrapolate further for the Bedouin population past fitting the carriage model. Further analyses could explore higher transmission populations such as the Bedouin population.
Though we identified some of the top serotypes using our method of ranking, it is important to note that looking at individual serotypes can be sometimes misleading because the data are so noisy. Previous work demonstrates that estimating changes for specific serotypes can be subject to a high degree of noise, sometimes leading to inaccurate projections [17]. Future work could incorporate additional predictors of serotype prevalence growth as they become available.
Having a strong understanding of the patterns of serotype replacement could be important as newer higher-valent vaccines are developed and introduced. With higher-valency vaccines, we need to be able to predict both the decline in vaccine-serotypes and the increase in non-vaccineserotypes. We can use this model to quantify changes in serotype replacement in other populations using the same analysis. This model may also help to optimize future serotype compositions.
Data Availability
All code and simulated versions of the data are available online at https://github.com/mailephillips/post-pcv-expansion. Real data were not available for public use, so data were simulated for the purposes of sharing the code online.
SUPPLEMENTAL INFORMATION
S1 Text. Carriage Model
Figure S1. Observed and expected prevalence plot, by year (Jewish children). Each panel shows the observed (black points) and model-fitted estimates of prevalence (different-colored dots and bars represent point estimates and 95% credible intervals), by year of the study period.
Figure S2. Observed and expected prevalence plot, by year (Bedouin children). Each panel shows the observed (black points) and model-fitted estimates of prevalence (different-colored dots and bars represent point estimates and 95% credible intervals), by year of the study period.
Figure S3. Observed versus fitted prevalence ratios (Jewish children). The plot shows the observed (black points) and model-fitted prevalence ratios (different-colored dots and bars represent point estimates and 95% credible intervals). The size of the different-colored dots indicates the number of isolates over the study period from the raw data (smaller dots = less isolates, larger dots = more isolates). Ratios shown are for the first to last years of the study period
. In instances where the observed prevalence in both Year 0 and Year 6 were 0, the observed prevalence ratio is indicated to be one. In cases where the observed prevalence in Year 0 was 0 and but the observed prevalence in Year 6 was nonzero, no observed point is shown (serotypes 20, 27, 29, and 7B).
Figure S4. Observed versus fitted prevalence ratios (Bedouin children). The plot shows the observed (black points) and model-fitted prevalence ratios (different-colored dots and bars represent point estimates and 95% credible intervals). The size of the different-colored dots indicates the number of isolates over the study period from the raw data (smaller dots = less isolates, larger dots = more isolates). Ratios shown are for the first to last years of the study period
.
Figure S5. Estimates for log relative risk ratios for pre-PCV period. Each serotype-specific log relative risk ratio with its 95% credible interval is shown for the pre-PCV period (β0,i), estimated from the model. The serotypes are shown from highest to lowest log relative risk ratios, separately for Jewish (A) and Bedouin (B) children.
Figure S6. Estimates for log relative risk ratios for early post-PCV period. Each serotype-specific log relative risk ratio with its 95% credible interval is shown for the early post-PCV period (β1,i), estimated from the model. The serotypes are shown from highest to lowest log relative risk ratios, separately for Jewish (A) and Bedouin (B) children.
Figure S7. Estimates for log relative risk ratios for late post-PCV period. Each serotype-specific log relative risk ratio with its 95% credible interval is shown for the late post-PCV period (β1,i), estimated from the model. The serotypes are shown from highest to lowest log relative risk ratios, separately for Jewish (A) and Bedouin (B) children.
Figure S8. Estimated additional cases gained as a result of increased carriage in children, by serotype. Each panel shows an individual serotype over the study period for all ages. Each full bar represents the total number of observed IPD cases for that stratum. The different colors denote the fraction of total cases attributable to different groups of serotypes: expected non-vaccine-targeted serotype cases had PCV13 not been introduced (solid red), estimated additional cases gained that are not explained by this model (hatched blue), and the estimated additional cases gained as a result of increased carriage in children (hatched red).
Table S1. Penalized deviance for all model combinations fit. The deviance information criterion (DIC) used to determine the best fit model is shown for all combinations of models for both Jewish and Bedouin datasets. An “X” indicates whether that particular feature was used in building the specified model.
Table S2. Serotype-specific parameters estimated from model: Jewish children under 5 in Israel, November 2009-July 2016. Estimates for the serotype-specific pre-PCV relative risk ratios (RRRs) (eβ0,i), early post-PCV RRRs (eβ1,i), and late post-PCV RRRs (eβ2,i) are shown for the Jewish carriage model, with their 95% credible intervals.
Table S3. Serotype-specific parameters estimated from model: Bedouin children under 5 in Israel, November 2009-July 2016. Estimates for the serotype-specific pre-PCV relative risk ratios (RRRs) (eβ0,i), early post-PCV RRRs (eβ1,i), and late post-PCV RRRs (eβ2,i) are shown for the Bedouin carriage model, with their 95% credible intervals.
Funding
This work was supported by the National Institutes of Health/National Institute of Allergy and Infectious Diseases [grant numbers R01-AI123208, R01-AI137093] and the Bill and Melinda Gates Foundation (OPP1176267). The funding agency was not involved in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication. The corresponding author had full access to all the data in the study and had final responsibility for the decision to submit for publication.
Conflicts of Interest
DMW has received consulting fees from Pfizer, Merck, GSK, and Affinivax, and is Principal Investigator on a grant from Pfizer to Yale University. RD has received consulting fees from Pfizer, MSD and MeMed; research grants from Pfizer and MSD; speaker fees from Pfizer. GRY has received consulting fees and research funding from Pfizer and research support from GSK. All other co-authors declare no potential conflict of interest.
