- Research article
- Open Access
- Open Peer Review
Dynamic models of pneumococcal carriage and the impact of the Heptavalent Pneumococcal Conjugate Vaccine on invasive pneumococcal disease
BMC Infectious Diseasesvolume 10, Article number: 90 (2010)
The 7-valent pneumococcal conjugate vaccine has been introduced in national immunisation programmes of most industrialised countries and recently in two African GAVI eligible countries (Rwanda and The Gambia). However the long term effects of PCV are still unclear, as beneficial direct and herd immunity effects might be countered by serotype replacement.
A dynamic, age-structured, compartmental model of Streptococcus pneumoniae transmission was developed to predict the potential impact of PCV7 on the incidence of invasive disease accounting for both herd immunity and serotype replacement effects. The model was parameterised using epidemiological data from England and Wales and pre and post-vaccination surveillance data from the US.
Model projections showed that serotype replacement plays a crucial role in determining the overall effect of a PCV7 vaccination programme and could reduce, negate or outweigh its beneficial impact. However, using the estimate of the competition parameter derived from the US post-vaccination experience, an infant vaccination programme would prevent 39,000 IPD cases in the 20 years after PCV7 introduction in the UK. Adding a catch-up campaign for under 2 or under 5 year olds would provide a further reduction of 1,200 or 3,300 IPD cases respectively, mostly in the first few years of the programme.
This analysis suggests that a PCV vaccination programme would eradicate vaccine serotypes from circulation. However, the increase in carriage of non-vaccine serotypes, and the consequent increase in invasive disease, could reduce, negate or outweigh the benefit. These results are sensitive to changes in the protective effect of the vaccine, and, most importantly, to the level of competition between vaccine and non-vaccine types. The techniques developed here can be used to assess the introduction of vaccination programmes in developing countries and provide the basis for cost-effectiveness analyses.
Despite the considerable burden of disease caused by Streptococcus pneumoniae infection in developing countries, the available 7-valent pneumococcal conjugate vaccine (PCV7) introduced in 2000 in the US and included in a number of European countries routine immunisation schedules , has yet to be launched in the majority of low-income settings where it is most needed . If on one hand the high price of the vaccine has been a limiting factor, discussions are also ongoing in regards to the feasibility of the programme as well as its long term effectiveness. Though in fact the vaccine has been shown to be effective against invasive and non invasive disease [3–8], the current licensed product contains seven serotypes only (14, 6B, 19F, 23F, 4, 9V, 18C) and it cannot be expected to provide any protection against carriage or disease caused by most other serotypes. On the contrary, evidence from several countries suggests that serotype replacement in carriage generates as a consequence of vaccination [9–12]. An increase in carriage prevalence of serotypes not contained in the vaccine as well as in disease has now been reported suggesting that replacement in carriage may also lead to replacement in disease [13–16].
While surveillance systems can and should closely monitor country-specific trends in the post vaccination era, there is the need to investigate the potential scenarios that the longer term effects of these interventions might produce and to evaluate the economic acceptability of such interventions considering medium to longer term consequences such as herd immunity and serotype replacement.
Transmission dynamic models have now been extensively used to describe the dynamics of infectious diseases  and to predict the effects of vaccination against many childhood infections [18–24]. These models are needed to evaluate the amount of herd immunity that can be generated by these programs . Also, by describing the natural history of the infection, they can provide insights on the magnitude of other indirect effects, such as serotype replacement, and on the aspects that need to be considered when evaluating the overall impact of these mechanisms. Here we use population-based data to develop, parameterise and apply an age-structured transmission dynamic model  and to explore the impact of different PCV7 vaccination strategies on carriage of vaccine- and non-vaccine pneumococcal serotypes (VT and NVT) and on invasive pneumococcal disease (IPD). This model differentiates from previous more theoretical work [27, 28] by the inclusion of an age structure, heterogeneous mixing and the calculation of IPD rates. An age-structured model was considered essential to assess the public health and economic impact of alternative vaccination policies, since the risks of pneumococcal carriage and disease are highly correlated with age, and vaccination programs are targeted at specific age groups (routine vaccination, booster doses, and catch-up campaigns). Our work is also different from other recent and more applied models which looked at the effects of infant vaccination on the epidemiology of S. pneumoniae infection, by using pre and post invasive disease surveillance data from, respectively, Australia and the US [29, 30]. Neither of these two studies, in fact, looked at the impact of vaccination on the ecology of the bacterium and, in particular, on the possibility that non-vaccine serotypes will replace serotypes eliminated by the vaccine. This phenomenon is currently being observed in England and Wales post vaccine introduction http://www.hpa.org.uk and it is having major implications for the overall impact of the programme.
The model presented here generates projections for England and Wales and evaluates the policies that were under consideration at the Department of Health for the UK at the time of introduction of PCV7. However, future changes to the vaccine formulations and/or schedule can be considered with the current model and it can also be adapted to look at other countries programs. Mathematical models can provide insights on the overall impacts produced by alternative public health interventions and should be used by policy makers as an additional tool to inform their decision.
The following three distinctive datasets were used to parameterise the model:
pneumococcal carriage prevalence data from England and Wales  to estimate the pre-vaccination steady state force of infection
the national enhanced surveillance of IPD in England and Wales  to derive the case:carrier ratio (i.e. the proportion of infection that progress into a disease state)
the Active Bacterial Core Surveillance data (ABC)  which, by providing unique at the time of the analysis pre and post vaccination data allowed the estimation of key model parameters such as degree and duration of vaccine protection and competition between vaccine and non-vaccine serotypes.
In the following a brief description of the three datasets is given.
England &Wales longitudinal carriage data
A longitudinal study was carried out by collecting nasopharyngeal swabs every month from 489 individuals including 121 pre-school children (<3 yrs.) for a period of 10 months, from October 2001 to July 2002 in England and Wales [31, 34]. Swabs were cultured for pneumococci according to the WHO guidelines . The number of positive pneumococcal swabs was 932 among the 3,753 swabs taken, (25%), and 34 pneumococcal serotypes were identified (Figure 1). Among these serotypes, 6A, 6B, 14, 19F and 23F were the major serotypes identified from positive swabs (75% in children and 52% in adults). The serotypes included in PCV7 and 6A were identified from 80% of positive swabs from children aged under 5 and 67% of those from subjects aged 5 and older.
England &Wales national surveillance of invasive pneumococcal disease
The annual number of IPD cases was estimated as the average of the overall number of cases that were reported to the reconciled enhanced surveillance system of invasive pneumococcal disease for England & Wales [36, 37] and over the 2003/04 and 2004/2005 epidemiological years. This surveillance system is currently the largest national data source on invasive disease caused by S. pneumoniae infection and on the isolated serotypes (Figure 2). The average annual number of IPD cases during the pre-vaccination period (between 2003/04 and 2004/05) in England and Wales was 6,184.
Active Bacterial Core surveillance of invasive pneumococcal disease
The Active Bacterial Core Surveillance (ABCs) is an active laboratory and population-based surveillance system for invasive bacterial pathogens of public health importance. From the ABCs system, IPD cases were identified from 8 counties in the US (California, Connecticut, Georgia, Maryland, Minnesota, New York, Oregon and Tennessee). The total population is estimated at 18,041,585 (2004 post-censual estimates) and the total number of IPD cases was 22,365 for the period 1998-2004 . As the PCV7 was introduced in mid 2000, the average number of IPD cases during 1998 and 1999 was used as the pre-vaccination figure and the effects of PCV7 are presented in Figure 3.
The model consists of 100 cohorts of individuals (0, 1, 2, 3,..., 99) each corresponding to one year of age and each of equal size (nc = 0.01), with a total stable population (i.e. births equal deaths) which adds up to 1 . Individuals are born into age cohort 0 at the start of the year and live to the age of 99 years, at which point they die. Ageing is thus modelled at the end of each year, with individuals in age group i moving up to become age group i+1. Each age cohort goes through the transmission cycles according to the infection process outlined in the following paragraph.
A compartmental transmission dynamic model is developed to examine the effects of PCV7 on VT and NVT carriage prevalence and, ultimately, on the overall incidence of IPD. Serotypes 4, 6B, 9V, 14, 18C, 19F, 23F as well as serotype 6A are all included in the VT group. Though 6A is not explicitly included in PCV7, carriage of this serotype is assimilated to the VT one as PCV7 was shown to elicit strong opsonic capacity against serotype 6A  and a significant reduction of 6A carriage was observed following vaccination [10, 39, 40]. In the model, carriage of any of the remaining pneumococcal serotypes is considered a NVT carriage episode.
The model consists of a set of ordinary differential equations, which describes the transmission and clearance of carriage with and without vaccination. The model structure is shown in Figure 4 and details of the equations are provided in Appendix 1. A Susceptible-Infected-Susceptible (SIS) type model is chosen (i.e. no immunity after infection) and individuals can get co-infected with one serotype from the VT group and one from the NVT group. Thus, a model that did not consider vaccination would require four compartments: S, non-carriers; V, carriers of VT pneumococci; N, carriers of NVT pneumococci; B, carriers of both VT and NVT pneumococci. To enable consideration of vaccines that provide partial protection against carriage, the model has four additional compartments for individuals that have a degree (i.e. vaccine efficacy against VT carriage acquisition) of vaccine-induced protection (Sv, Vv, Nv, Bv).
Uninfected, unprotected individuals (S) become infected according to the age-specific force of infection of vaccine type (λ Vi ) and non-vaccine type (λ Ni ) pneumococci. Individuals carrying VT become co-colonised at a rate c N λ Ni where c N is the relative susceptibility of an individual colonised with VT to acquiring NVT compared with an uncolonised individual. Similarly NVT carriers become co-infected at a rate c V λ Vi . Once infected, individuals can recover from carriage and go back to a susceptible state (V→S or N→S) or, if they were co-colonised, to carriage of a unique type (B→V or B→N). This happens at an age-dependent rate r i which is the same whether or not the individual was in the colonised state (V or N) or in the co-colonised state (B).
Following vaccination, all vaccinated individuals move into the equivalent vaccine-protected compartment according to their carriage status, and gain partial protection against VT acquisition. The force of infection of VT pneumococci, among vaccine-protected individuals, is reduced by a factor (1-γ i ) where γ i is the degree of protection. The force of infection for NVT among vaccine-protected individuals is the same as among unprotected individuals. Vaccine-protected individuals may lose their protection and go back to the equivalent unprotected compartment at a waning rate w (= 1/δ, where δ is the average duration of protection) that is constant throughout all ages and across all compartments. Vaccination is assumed not to affect infectiousness: vaccinated VT carriers are assumed to be as infectious as unvaccinated VT carriers.
Vaccination programs in the model are implemented as discrete events, rather than as a continuous process at some constant rate: in a routine vaccination program all infants in a cohort are assumed to be vaccinated at the same time according to the vaccination schedule; and in a catch-up campaign individuals in the targeted age group are assumed to be vaccinated simultaneously.
We assume that the risk of disease occurs at the time of carriage acquisition. The proportions of those infected who develop IPD, Case: Carrier ratios, depend on their age and the infecting serotype (VT or NVT).
The average duration of carriage was estimated using a previously developed modelling framework [41, 42] and considering the following age classes: 0-1, 2-4, 5-17 and 18+. Carriage episode was found to be significantly longer in young children (72 days in the 0-1 year of age), and then steadily decreasing with 28 days in the 2-4 years of age, 18 in the 5-17 years and 17 days in the 18+. Age-dependent recovery rates (1/average duration of carriage) were assumed to be the same for VT and NVT (rVi = rNi) [41, 43], and also for vaccinated and unvaccinated individuals (Table 1).
Forces of infection
The forces of infection for VT and NVT were age- and time-dependent, varying as a function of the number of vaccinated and unvaccinated carriers in each age class j and the rate of effective contacts between individuals in the age class i and those in age class j (β ij ), as follows:
The derivation of the transmission rates βij and the underlying assumptions on the mixing patterns adopted in the model are described in Appendix 2 and follow the techniques originally developed for sexually-transmitted infections . In brief, the contact patterns within and between age groups ranges from fully assortative to proportionate, governed by a single parameter, ε (0≤ε≤1) (ε = 0 proportionate; ε = 1 assortative), which is estimated by fitting model projections to the US post vaccination surveillance data as described below.
Mixing and Competition parameters
The introduction of PCV7 for routine use in children aged <5 years in the United States yields valuable information about the impact of vaccination at a population level. We used national surveillance data on the number of cases of invasive disease caused by VT and NVT before and in the four years after the vaccine introduction  to estimate mixing pattern and competition parameters. In particular, the impact on VT carriage among the unvaccinated age groups provided information on the mixing patterns in the population (ε). The speed and extent of the increase in NVT disease provided the estimate of cN (how much VT carriage protects against acquisition of NVT). As there is little information in this data set on cV, we used the value, cV = 0.5, as a baseline value, and analysed the sensitivity using the values of 0 and 1. We optimised ε and cN by fitting to the incidence of IPD in each age group from the US data by minimizing the Poisson deviance. The reported annual coverage level for 1+ doses in children between 19 and 35 months in 2003 of 86% (Dr Cynthia Whitney, personal communication) was assumed as the baseline annual vaccine coverage of routine vaccination from 2001 onward. Vaccination coverage in the first year of the program was set at 43% for 0 (routine) and 1 (catch-up) year olds as the vaccine was introduced around the middle of the year (July 2000). As the vaccine coverage levels varied each year, the coverage assumed in the model was varied as well, from 80% to 90% with an increment of 2% for the sensitivity analysis.
The risk of disease given pneumococcal infection
The age-specific proportions of VT and NVT pneumococcal infections that lead to IPD (case: carrier ratio) were estimated by fitting a model to the age-specific incidence of these diseases attributable to pneumococci as reported from the national enhanced surveillance of pneumococcal disease in the pre-vaccination years . A total of eighteen case: carrier ratios for VT and NVT were estimated for nine age groups, <2 months, 2-11 months, 1-4, 5-14, 15-24, 25-44, 45-64, 65-74, 75+ years.
Vaccine associated parameters
We assume that all individuals who received PCV7 move to the vaccine protected group. The average duration and degree of vaccine-induced protection were estimated by fitting the model to the ABCs data from the pre- and post- vaccination US experience. The sensitivity of the model results to these parameter estimates was assessed.
Modelling PCV7 vaccination
The alternative vaccination strategies investigated here are the ones that were relevant for the UK routine immunisation schedule:
◦ Strategy 1: routine vaccination at 2 and 4 months of age, plus a booster dose at 12 m;
◦ Strategy 2: Strategy 1 plus catch-up campaign targeting all those under the age of 24 m;
◦ Strategy 3: Strategy 1 plus catch-up campaign targeting all those under the age of 60 m.
For each vaccination strategy the model predicted the number of cases of IPD caused by VT and NVT in each age group over time. The coverage of routine and catch-up vaccinations for PCV7 in the UK was assumed to be the same as for the meningococcal serogroup C conjugate vaccine pre-school catch-up (89% for 0 year olds, 84% for 1 year olds, and 76% for the rest) . Catch-up campaigns were implemented at the same time as the routine vaccination was introduced.
There were three stages to the analysis. First, a steady state pre-vaccination model estimated age-stratified values for the forces of infection and case:carrier ratios for VT and NVT. Second, a dynamic post-vaccination model estimated key vaccine parameters, interaction between VT and NVT and the level of assortativeness of mixing pattern from US surveillance data. Third, the dynamic model was used to assess the impact of alternative vaccination strategies in England and Wales.
The pre-vaccination model was programmed in Excel to estimate forces of infection for VT and NVT and for the following age groups: 0-1, 2-4, 5-9, 10-19, 20-39 and 40-99 using the carriage prevalence data available for England and Wales and a given value for the competition parameter (cN). The fully assortative and fully proportionate mixing matrices between the six age-groups were also generated and later used as two extremes in the transmission model. Age-specific case:carrier ratios were then derived fitting the model to the age distribution of IPD cases caused by, respectively, vaccine and non-vaccine serotypes in the pre-vaccination era. The procedure, which required the minimisation of a Poisson deviance using the SOLVER in Excel, was run for both England and Wales and the US to calculate the risk of developing disease when colonised.
The transmission dynamic model was programmed in Berkeley Madonna (R. I. Macey & G. F. Oster, UC Berkeley, CA, USA) and fitted to the pre- and post- vaccination IPD data from the US to estimate degree and duration of protection of the vaccine against invasive disease, and the mixing parameter (ε) for different values of the competition parameter (cN). The forces of infection, case:carrier ratios and extreme mixing matrices were updated for each value of cN. The estimate of cN was derived by finding the value that minimised the deviance.
Once the parameters were generated, the epidemiological model was used to assess the impact of alternative vaccination strategies in England and Wales. The system was solved using the Euler method to integrate ordinary differential equations with fixed time steps of 0.001 years. The model simulated 50 years (five years pre-vaccination and 45 years after vaccination).
Graphical comparison between the observed and predicted VT and NVT IPD cases by age group from 1998 to 2004 in the 8 counties covered by the Active Bacterial Core Surveillance is presented in Figure 5. The baseline value of the PCV7 coverage with 1+ doses was assumed to be 86%. The sensitivity of estimates of the degree of vaccine protection and duration of protection to the assumed PCV7 coverage was investigated. Increasing the vaccine coverage assumed may account for some protection derived by individuals who were only partially vaccinated; reducing the vaccine coverage assumed is equivalent to reducing the take of the vaccine to less than 100%. Parameter estimates are presented in Table 2 according to six vaccine coverage values from 80% to 90% and show that increases in the coverage level produced some reduction in the degree and duration of protection as expected. However, for coverage levels within this range, the estimates of the competition and mixing parameters were not sensitive to the assumed vaccination coverage. These results suggest that individuals who were carrying VT were partly protected against NVT acquisition (cN = 85%, 15% protection), and that population mixing was closer to assortative than proportionate (ε = 0.87) in line with a recent European contact study . The contour plot presented in Figure 6 describes the 95% confidence area of the degree and duration of the vaccine protection while the mixing and competition parameter were fixed (respectively 0.87 and 0.85).
Model projections with vaccination
Seven pairs of the degree and duration of PCV7 vaccine protection were chosen within the 95% confidence area (including the best fitted pair with the baseline US PCV7 coverage) and showed that the UK transmission model predictions of IPD changes are insensitive to the assumptions. The fitting results indicated that the higher the degree is the shorter the duration of protection. The baseline parameters used for the following simulations of England and Wales predictions are the ones estimated when fixing the coverage level to 86% (Table 1). Age-specific case:carrier ratios that are used for the UK predictions are shown in Table 3.
In the long term, the model predicted that all three strategies had a similar impact and eliminated transmission of VT. Strategy 3, which has the largest catch-up, has the most rapid impact on carriage and IPD, while Strategy 1 with no catch-up was the slowest (Figure 7). The model results presented in Figure 7 also showed an increase in the incidence of NVT carriage and IPD after the introduction of vaccination. The extent of this replacement of VT by NVT is governed by the reduction in VT.
In the long term, the model predicts a 63% reduction of all IPD cases in children aged less than 5 years old and 35% reduction in the rest of the population, which amounts to the prevention of 2,300 cases of IPD annually after 20 years of the PCV7 introduction with all three vaccine programme strategies. Following the results for the impact on carriage, Strategy 1 is the least effective in terms of disease prevented. The routine vaccination and booster for one year old (Strategy 1) reduces 6000 IPD cases within the first five years of the programme, and an additional catch-up for children up to 23 months of age (Strategy 2) prevents extra 1,200 IPD cases. Extending the catch-up to age 59 month (Strategy 3) prevents further 3,300 IPD cases. The impact of vaccination on the long term incidence of VT and NVT disease by age group and vaccine programme strategies is illustrated in Figure 8.
If the vaccine was less protective than assumed here, (i.e. lower duration or degree of protection), the program achieved lower coverage, or the population was less assortative, then transmission of vaccine types may not be eliminated. Extensive sensitivity analysis was conducted and the results are presented in Table 4. However, three main conclusions were drawn. Firstly, a reduction in the duration and/or degree of protection would prevent elimination of VT. However, increasing the values of these parameters (90% degree or 20 years) would not provide much further reduction in disease. Secondly, an increase in the value of the competition parameter cN (lower protection from VT carriage) would result in lower replacement of VT by NVT after PCV7 introduction and thus in a small number of IPD cases while a lower competition parameter (i.e. higher protection from VT carriage) would result in more IPD cases than in the pre-vaccination period. Thirdly, lower values of ε (less assortative mixing pattern) would facilitate and hasten elimination of VT transmission with more IPD case reductions. On the contrary, higher values of ε would result in VT continuing to be carried due to less herd immunity effects of the vaccine.
The transmission model of S. pneumoniae infection presented in this paper enables investigation of the potential effects of infant PCV7 vaccination. It builds on previous S. pneumoniae models [27, 28] and, by including age-structure, allows investigation of the impact of vaccination on the age-specific incidence of invasive disease.
Simulation results indicated that all of the vaccination strategies considered in this paper were sufficient to eliminate VT transmission in England and Wales with base case values for vaccination coverage and other parameters. Although the model predicts some replacement of VT with NVT, the annual total of IPD cases was reduced from 6,184 to 3,765 (38% reduction) in the long-term (20 years after the PCV7 introduction) and reduction of 39,000 IPD cases accumulated over the period. Catch-up campaigns for children aged under 2 and 5 prevent additional 1,600 and 4,100 IPD cases respectively compared to routine vaccination alone during the first twenty years of the PCV7 program (1,200 and 3,300 cases reduction in the first five years). However, this conclusion is sensitive to vaccine and other model parameters. Small reductions in vaccine coverage, in the duration or degree of protection, from their base case values are sufficient to prevent elimination of VT being achieved. We cannot therefore exclude the possibility that PCV7 vaccination fails to eliminate VT transmission in the UK. On the other hand, NVT dynamics post vaccination are extremely sensitive to the level of competition between VT and NVT, and the model predicted an overall increase of the number of IPD cases pre to post vaccination if the level of protection of VT from NVT acquisition were almost complete.
The competition parameter, cV (the relative susceptibility of NVT carriers, compared with uncolonised individuals, to acquiring VT carriage), was assumed to be 50% since there is little information from post-vaccination surveillance on the protective effect of NVT carriage on VT acquisition. Sensitivity analysis suggested that the reduction in annual IPD cases after 10 years of PCV7 introduction is not sensitive to this parameter, with a maximum difference of 1% between scenarios (cV = 0 or = 100%). On the other hand, the degree of serotype replacement is sensitive to cN (the relative susceptibility of VT carriers, compared with uncolonised individuals, to acquiring NVT carriage), which was estimated by fitting to the US post vaccination surveillance data. The sensitivity analysis indicated that the lower this value, the more replacement of VT by NVT will occur after the introduction of PCV7 in the community. Based on the UK surveillance after the PCV7 introduction, the replacement appears greater than the replacement in the US http://www.hpa.org.uk and the cause of this difference is currently being investigated. Potential contributing factors include differences in secular trends in prevalence and antimicrobial sensitivity of specific serotypes at the time of introduction , differences in sensitivity of the respective surveillance systems, or differences in clinical practice post PCV introduction.
The mechanism of competition incorporated in the model may influence the projected outcomes . Here, the mechanism assumed was of a reduction of the acquisition rate of other serotypes when already colonised. This mechanism is supported by a recent Danish longitudinal study of pneumococcal infection , which preferred it to a mechanism of competition through an increase in the clearance rate for co-colonised individuals. Our exploration of alternative hypotheses for competition is limited by a lack of co-colonisation data. Clearly, more work is needed in this area and laboratory techniques that enable detection of carriage of multiple serotypes will facilitate studies of the interactions between different pneumococcal serotypes.
The mixing pattern determines herd immunity effects of PCV7 vaccination among the unvaccinated age groups. From the impact observed in the US experience, we estimated that the mixing between individuals is more assortative than proportionate. Thus, similarly to other infections , contacts with others in the same age group play an important role in the transmission of pneumococcal carriage. Population mixing patterns are influential parameters in transmission models. True mixing patterns are unknown, and model parameters are often estimated from the pre-vaccination epidemiology after making arbitrary assumptions about the structure of contacts . Here, we use post vaccination herd immunity effects to inform the level of assortativeness of mixing. Alternative approaches to measure age-specific contact patterns directly through diary-based methods are currently being explored using recently collected contact pattern data [53–56].
The relationship between acquisition of carriage and onset of invasive and non-invasive disease is not well understood. Here we assume that the risk of pneumococcal disease occurs at the time of acquisition of carriage, with age- and type-specific proportions developing IPD. An alternative assumption is that this risk is spread over the entire duration of carriage. We do not envisage that adopting this latter mechanism would markedly affect our results, especially as we have assumed that the duration of carriage is the same for both vaccine and non-vaccine types.
One critical aspect of the current work is that the model categorises pneumococcal serotypes as either vaccine-type (including 6A) or non-vaccine type and assumes that the characteristics are homogeneous within each of these two categories. Firstly, the inclusion of serotype 6A in the vaccine group was implemented due to existing evidence that the 7-valent pneumococcal conjugate vaccine induces cross-protection against 6A [14, 57]. However, the newly identified serotypes 6C and 6D [58, 59], which do not benefit from cross protection [58, 60–62] should be excluded from the vaccine group in future modelling work. Secondly, there is considerable heterogeneity between individual serotypes in, for example, their transmissibility, duration of infection, ability to co-colonise, ability to prevent co-colonization with other serotypes, and potential to cause disease [63–65]. All these factors may influence the response of individual serotypes and the pneumococcal population as a whole to the introduction of infant PCV7 vaccination. Studying these effects would require an individual-based model in order to incorporate many circulating serotypes. Within such a framework it would also be possible to begin to investigate the impact of acquired immunity on pneumococcal transmission by incorporating a mechanism for generating type-specific and/or type-independent immunity [66–69].
The impact of pneumococcal vaccination on non-invasive diseases such as Acute Otitis Media (AOM) and Community Acquired Pneumonia (CAP) is not considered in this paper although some level of protection has been reported against these conditions in both developed [4, 70] and developing countries [5, 71]. Although an individual case of non-invasive pneumococcal disease (NIPD) is generally less severe than IPD, the economic cost of NIPD is higher than that of IPD due to the much greater number of cases in the community. Also, CAP is one of the leading causes of childhood mortality in developing countries, and hence, understanding the impact of pneumococcal vaccine on NIPDs is essential to estimate the overall cost-effectiveness of the program in these settings . The impact of PCV7 vaccination programs on the incidence of pneumococcal AOM and CAP or even on the prevalence of carriage in the general population is currently not well documented. Whilst the capacity of different serotypes to cause AOM is expected to be similar, there are large differences in their risk of causing IPD: some NVT serotypes are very unlikely to cause IPD. The change in prevalence of NVT serotypes that commonly cause IPD can be inferred from the IPD data, but the change in prevalence of NVT serotypes that rarely cause IPD is not observed. The impact on NIPD may not be in line with the effect on IPD. Parameterising models to investigate the effect of the vaccine on AOM or CAP in the community requires more knowledge of the impact of vaccination programs on carriage and NIPD.
The outputs of any model are limited by the validity of the model assumptions and accuracy of the parameter estimates. In this study, considerable attention was paid to estimating appropriate parameter values and to assessing the sensitivity of the results to these values. Data from the longitudinal study in UK families were used to derive the initial prevalence of carriage, recovery rates and the forces of infection. Surveillance data on IPD cases pre- and post- PCV7 introduction in the US were used to estimate vaccine parameters, mixing patterns and to infer the degree of competition between VT and NVT. Moreover, the model developed here may also be useful for other countries considering the introduction of the pneumococcal conjugate vaccine and to estimate effects of alternative vaccine strategies. In future, these models will be improved by refining parameter estimates using the observed effects of the UK programme. The UK experience with the Hib and MenC vaccines [73–75] showed that the post-vaccine era may be far from straightforward and that continued monitoring the duration of vaccine protection and the effects of vaccination on carriage may be critical to the success of the pneumococcal conjugate vaccine program.
The effect of PCV7 introduction on the burden of IPD is assessed using an age structured dynamic model which was parameterised with UK and US data. Model projections show the importance of considering the existence of competition between serotypes contained and not contained in the vaccine formulation in order to being able to realistically predict possible scenarios. The current model will be extended to incorporate non invasive pneumococcal disease and will form the basis for cost-effectiveness analysis for alternative vaccination programmes, especially in developing countries where the burden of pneumococcal disease is extremely high and the financial resources are limited.
A.1. Dynamic model structure in equations
Where i is the age group, νi represents the vaccine coverage in age group i, π i(t) is a function, describing the vaccination program in cohort i to occur at certain time as scheduled, rVi and rNi are, respectively, the proportion that recover from VT and NVT infection, ω is the rate of waning of vaccine-induced protection, c N (c V ) represent the relative susceptibility of VT (NVT) carriers, compared with uncolonised individuals, to acquiring NVT (VT) carriage (i.e. competition parameters), γ is the vaccine efficacy against VT carriage, λ Vi (t) and λ Ni (t) represent the age-specific force of infection for, respectively, VT and NVT pneumococci at time t.
A.2. Mixing pattern
Given values of the pneumococcal transmission parameters (in particular the competition parameters cV and cN), the effective transmission coefficients (βVij and βNij) were calculated by fitting a model to the number of VT and NVT carriers in each age group at steady state (Figure 9). Matrices for VT and NVT were calculated assuming fully assortative (βa) and proportionate (βp) mixing between the six age groups. A parameter ε, the assortativeness of mixing, is introduced to define the weight that is given to each of these matrices :
Hence, if ε = 1, then full weight is given to the fully assortative matrices and mixing patterns is totally assortative whereas if ε = 0 full weight is given to the proportionate matrix. A value of ε is estimated when fitting to the pre and post vaccination disease incidence data for the US.
Observed and fitted prevalence of carriage of pneumococcal vaccine and non-vaccine serotypes by age group in England and Wales are presented in Figure 9. Pneumococcal carriage data were collected as part of a longitudinal pneumococcal family study in the UK that ran from October 2001 to July 2002 .
De Carvalho GH, Muscat M, Monnet DL, Giesecke J, Lopalco PL: Use of seven-valent pneumococcal conjugate vaccine (PCV7) in Europe, 2001-2007. Euro Surveill. 2009, 14: pii-19159
CDC: Progress in introduction of pneumococcal conjugate vaccine--worldwide, 2000-2008. MMWR Morb Mortal Wkly Rep. 2008, 57: 1148-1151.
Black SB, Shinefield HR, Fireman B, Lewis E, Ray P, Hansen JR, et al: Efficacy, safety and immunogenicity of heptavalent pneumococcal conjugate vaccine in children. Pediatr Infect Dis J. 2000, 19: 187-195. 10.1097/00006454-200003000-00003.
Black SB, Shinefield HR, Ling S, Hansen J, Fireman B, Spring D, et al: Effectiveness of heptavalent pneumococcal conjugate vaccine in children younger than five years of age for prevention of pneumonia. Pediatr Infect Dis J. 2002, 21: 810-815. 10.1097/00006454-200209000-00005.
Cutts FT, Zaman SM, Enwere G, Jaffar S, Levine OS, Okoko JB, et al: Efficacy of nine-valent pneumococcal conjugate vaccine against pneumonia and invasive pneumococcal disease in The Gambia: randomised, double-blind, placebo-controlled trial. Lancet. 2005, 365: 1139-1146. 10.1016/S0140-6736(05)71876-6.
Black S, Shinefield H, Baxter R, Austrian R, Bracken L, Hansen J, et al: Postlicensure surveillance for pneumococcal invasive disease after use of heptavalent pneumococcal conjugate vaccine in Northern California Kaiser Permanente. Pediatr Infect Dis J. 2004, 23: 485-489. 10.1097/01.inf.0000129685.04847.94.
O'Brien KL, Moulton LH, Reid R, Weatherholtz R, Oski J, Brown L, et al: Efficacy and safety of seven-valent conjugate pneumococcal vaccine in American Indian children: group randomised trial. Lancet. 2003, 362: 355-361. 10.1016/S0140-6736(03)14022-6.
Reingold A, Hadler J, Farley M, Harrison L, Lynfield R, Thomas A, et al: Direct and Indirect Effects of Routine Vaccination of Children with 7-Valent Pneumococcal Conjugate Vaccine on Incidence of Invasive Pneumococcal Disease --- United States, 1998-2003. MMWR Morb Mortal Wkly Rep. 2005, 54: 893-897.
Huang SS, Hinrichsen VL, Stevenson AE, Rifas-Shiman SL, Kleinman K, Pelton SI, et al: Continued impact of pneumococcal conjugate vaccine on carriage in young children. Pediatrics. 2009, 124: e1-11. 10.1542/peds.2008-3099.
Dagan R, Givon-Lavi N, Zamir O, Sikuler-Cohen M, Guy L, Janco J, et al: Reduction of nasopharyngeal carriage of Streptococcus pneumoniae after administration of a 9-valent pneumococcal conjugate vaccine to toddlers attending day care centers. J Infect Dis. 2002, 185: 927-936. 10.1086/339525.
Huang SS, Platt R, Rifas-Shiman SL, Pelton SI, Goldmann D, Finkelstein JA: Post-PCV7 changes in colonizing pneumococcal serotypes in 16 Massachusetts communities, 2001 and 2004. Pediatrics. 2005, 116: e408-e413. 10.1542/peds.2004-2338.
Sa-Leao R, Nunes S, Brito-Avo A, Frazao N, Simoes AS, Crisostomo MI, et al: Changes in pneumococcal serotypes and antibiotypes carried by vaccinated and unvaccinated day-care centre attendees in Portugal, a country with widespread use of the seven-valent pneumococcal conjugate vaccine. Clin Microbiol Infect. 2009, 15: 1002-1007. 10.1111/j.1469-0691.2009.02775.x.
Singleton RJ, Hennessy TW, Bulkow LR, Hammitt LL, Zulz T, Hurlburt DA, et al: Invasive pneumococcal disease caused by nonvaccine serotypes among alaska native children with high levels of 7-valent pneumococcal conjugate vaccine coverage. JAMA. 2007, 297: 1784-1792. 10.1001/jama.297.16.1784.
Hicks LA, Harrison LH, Flannery B, Hadler JL, Schaffner W, Craig AS, et al: Incidence of pneumococcal disease due to non-pneumococcal conjugate vaccine (PCV7) serotypes in the United States during the era of widespread PCV7 vaccination, 1998-2004. J Infect Dis. 2007, 196: 1346-1354. 10.1086/521626.
CDC: Invasive pneumococcal disease in children 5 years after conjugate vaccine introduction--eight states, 1998-2005. MMWR Morb Mortal Wkly Rep. 2008, 57: 144-148.
Hennessy TW, Singleton RJ, Bulkow LR, Bruden DL, Hurlburt DA, Parks D, et al: Impact of heptavalent pneumococcal conjugate vaccine on invasive disease, antimicrobial resistance and colonization in Alaska Natives: progress towards elimination of a health disparity. Vaccine. 2005, 23: 5464-5473. 10.1016/j.vaccine.2005.08.100.
Anderson R, May R: Infectious diseases of humans. 1991, Oxford, Oxford University Press, 1991
Anderson RM, May RM: Vaccination against rubella and measles: quantitative investigations of different policies. J Hyg (Lond). 1983, 90: 259-325.
Edmunds WJ, Heijden van de OG, Eerola M, Gay NJ: Modelling rubella in Europe. Epidemiol Infect. 2000, 125: 617-634. 10.1017/S0950268800004660.
Trotter CL, Gay NJ, Edmunds WJ: Dynamic models of meningococcal carriage, disease, and the impact of serogroup C conjugate vaccination. Am J Epidemiol. 2005, 162: 89-100. 10.1093/aje/kwi160.
Babad HR, Nokes DJ, Gay NJ, Miller E, Morgan-Capner P, Anderson RM: Predicting the impact of measles vaccination in England and Wales: model validation and analysis of policy options. Epidemiol Infect. 1995, 114: 319-344. 10.1017/S0950268800057976.
Gay NJ: The theory of measles elimination: implications for the design of elimination strategies. J Infect Dis. 2004, 189 (Suppl 1): S27-S35. 10.1086/381592.
Halloran ME, Cochi SL, Lieu TA, Wharton M, Fehrs L: Theoretical epidemiologic and morbidity effects of routine varicella immunization of preschool children in the United States. Am J Epidemiol. 1994, 140: 81-104.
Brisson M, Edmunds WJ, Gay NJ, Law B, De Serres G: Modelling the impact of immunization on the epidemiology of varicella zoster virus. Epidemiol Infect. 2000, 125: 651-669. 10.1017/S0950268800004714.
Edmunds WJ, Medley GF, Nokes DJ: Evaluating the cost-effectiveness of vaccination programmes: a dynamic perspective. Stat Med. 1999, 18: 3263-3282. 10.1002/(SICI)1097-0258(19991215)18:23<3263::AID-SIM315>3.0.CO;2-3.
Schenzle D: An age-structured model of pre- and post-vaccination measles transmission. IMA J Math Appl Med Biol. 1984, 1: 169-191. 10.1093/imammb/1.2.169.
Lipsitch M: Vaccination against colonizing bacteria with multiple serotypes. Proc Natl Acad Sci USA. 1997, 94: 6571-6576. 10.1073/pnas.94.12.6571.
Lipsitch M: Bacterial vaccines and serotype replacement: lessons from Haemophilus influenzae and prospects for Streptococcus pneumoniae. Emerg Infect Dis. 1999, 5: 336-345. 10.3201/eid0503.990304.
Sutton KL, Banks HT, Castillo-Chavez C: Estimation of invasive pneumococcal disease dynamics parameters and the impact of conjugate vaccination in Australia. Math Biosci Eng. 2008, 5: 175-204.
Snedecor SJ, Strutton DR, Ciuryla V, Schwartz EJ, Botteman MF: Transmission-dynamic model to capture the indirect effects of infant vaccination with Prevnar (7-valent pneumococcal conjugate vaccine (PCV7)) in older populations. Vaccine. 2009, 27: 4694-4703. 10.1016/j.vaccine.2009.05.062.
Hussain M, Melegaro A, Pebody RG, George R, Edmunds WJ, Talukdar R, et al: A longitudinal household study of Streptococcus pneumoniae nasopharyngeal carriage in a UK setting. Epidemiol Infect. 2005, 133: 891-898. 10.1017/S0950268805004012.
Melegaro A, Edmunds WJ, Pebody R, Miller E, George R: The current burden of pneumococcal disease in England and Wales. J Infect. 2006, 52: 37-48. 10.1016/j.jinf.2005.02.008.
Poehling KA, Talbot TR, Griffin MR, Craig AS, Whitney CG, Zell E, et al: Invasive pneumococcal disease among infants before and after introduction of pneumococcal conjugate vaccine. JAMA. 2006, 295: 1668-1674. 10.1001/jama.295.14.1668.
Melegaro A, Edmunds WJ: Cost effectiveness analysis of pneumococcal conjugate vaccination in England and Wales. Vaccine. 2004, 22: 4203-4214. 10.1016/j.vaccine.2004.05.003.
O'Brien KL, Nohynek H: Report from a WHO working group: standard method for detecting upper respiratory carriage of Streptococcus pneumoniae. Pediatr Infect Dis J. 2003, 22: 133-140.
George R, Melegaro A: Invasive pneumococcal infection England and Wales, 1999. CDR Weekly. 2001, 11:
Miller E, Waight P, Efstratiou A, Brisson M, Johnson A, George R: Epidemiology of invasive and other pneumococcal disease in children in England and Wales 1996-1998. Acta Paediatr Suppl. 2000, 435: 11-16. 10.1080/080352500300051477.
Lee H, Nahm MH, Burton R, Kim KH: Immune response in infants to the heptavalent pneumococcal conjugate vaccine against vaccine-related serotypes 6A and 19A. Clin Vaccine Immunol. 2009, 16: 376-381. 10.1128/CVI.00344-08.
Dagan R, Givon-Lavi N, Zamir O, Fraser D: Effect of a nonavalent conjugate vaccine on carriage of antibiotic-resistant Streptococcus pneumoniae in day-care centers. Pediatr Infect Dis J. 2003, 22: 532-540. 10.1097/00006454-200306000-00009.
Cheung YB, Zaman SM, Nsekpong ED, Van Beneden CA, Adegbola RA, Greenwood B, et al: Nasopharyngeal Carriage of Streptococcus pneumoniae in Gambian Children who Participated in a 9-valent Pneumococcal Conjugate Vaccine Trial and in Their Younger Siblings. Pediatr Infect Dis J. 2009, 28: 990-995. 10.1097/INF.0b013e3181a78185.
Melegaro A, Choi Y, Pebody R, Gay N: Pneumococcal carriage in United Kingdom families: estimating serotype-specific transmission parameters from longitudinal data. Am J Epidemiol. 2007, 166: 228-235. 10.1093/aje/kwm076.
Melegaro A, Gay NJ, Medley GF: Estimating the transmission parameters of pneumococcal carriage in household. Epidemiol Infect. 2004, 132: 433-441. 10.1017/S0950268804001980.
Cauchemez S, Temime L, Valleron AJ, Varon E, Thomas G, Guillemot D, et al: S. pneumoniae transmission according to inclusion in conjugate vaccines: Bayesian analysis of a longitudinal follow-up in schools. BMC Infect Dis. 2006, 6: 14-10.1186/1471-2334-6-14.
Garnett GP, Anderson RM: Factors controlling the spread of HIV in heterosexual communities in developing countries: patterns of mixing between different age and sexual activity classes. Philos Trans R Soc Lond B Biol Sci. 1993, 342: 137-159. 10.1098/rstb.1993.0143.
Whitney CG, Pilishvili T, Farley MM, Schaffner W, Craig AS, Lynfield R, et al: Effectiveness of seven-valent pneumococcal conjugate vaccine against invasive pneumococcal disease: a matched case-control study. Lancet. 2006, 368: 1495-1502. 10.1016/S0140-6736(06)69637-2.
Health Protection Agency. 2009, [http://www.hpa.org.uk]
Trotter CL, Edmunds WJ: Modelling cost effectiveness of meningococcal serogroup C conjugate vaccination campaign in England and Wales. BMJ. 2002, 324: 809-10.1136/bmj.324.7341.809.
Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, et al: Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 2008, 5: e74-10.1371/journal.pmed.0050074.
Black S: The Volatile Nature of Pneumococcal Serotype Epidemiology: Potential for Misinterpretation. Pediatr Infect Dis J. 2009.
Zhang Y, Auranen K, Eichner M: The influence of competition and vaccination on the coexistence of two pneumococcal serotypes. Epidemiol Infect. 2004, 132: 1073-1081. 10.1017/S0950268804002894.
Auranen K, Mehtala J, Tanskanen A, Kaltoft S: Between-strain competition in acquisition and clearance of pneumococcal carriage--epidemiologic evidence from a longitudinal study of day-care children. Am J Epidemiol. 2010, 171: 169-176. 10.1093/aje/kwp351.
Farrington CP, Kanaan MN, Gay NJ: Estimation of the basic reproduction number for infectious diseases from age-stratified serological survey data. Appl Statist. 2001, 50: 251-292.
Edmunds WJ, O'Callaghan CJ, Nokes DJ: Who mixes with whom? A method to determine the contact patterns of adults that may lead to the spread of airborne infections. Proc R Soc Lond B Biol Sci. 1997, 264: 949-957. 10.1098/rspb.1997.0131.
Hens N, Goeyvaerts N, Aerts M, Shkedy Z, van Damme P, Beutels P: Mining social mixing patterns for infectious disease models based on a two-day population survey in Belgium. BMC Infect Dis. 2009, 9: 5-10.1186/1471-2334-9-5.
Ogunjimi B, Hens N, Goeyvaerts N, Aerts M, van Damme P, Beutels P: Using empirical social contact data to model person to person infectious disease transmission: An illustration for varicella. Math Biosci. 2009, 218: 80-87. 10.1016/j.mbs.2008.12.009.
Wallinga J, Teunis P, Kretzschmar M: Using data on social contacts to estimate age-specific transmission parameters for respiratory-spread infectious agents. Am J Epidemiol. 2006, 164: 936-944. 10.1093/aje/kwj317.
Hsu HE, Shutt KA, Moore MR, Beall BW, Bennett NM, Craig AS, et al: Effect of pneumococcal conjugate vaccine on pneumococcal meningitis. N Engl J Med. 2009, 360: 244-256. 10.1056/NEJMoa0800836.
Park IH, Pritchard DG, Cartee R, Brandao A, Brandileone MC, Nahm MH: Discovery of a new capsular serotype (6C) within serogroup 6 of Streptococcus pneumoniae. J Clin Microbiol. 2007, 45: 1225-1233. 10.1128/JCM.02199-06.
Bratcher PE, Kim KH, Kang JH, Hong JY, Nahm MH: Identification of natural pneumococcal isolates expressing serotype 6D by genetic, biochemical and serological characterization. Microbiology. 2010, 156: 555-560. 10.1099/mic.0.034116-0.
Tocheva AS, Jefferies JM, Christodoulides M, Faust SN, Clarke SC: Increase in serotype 6C pneumococcal carriage, United Kingdom. Emerg Infect Dis. 2010, 16: 154-155. 10.3201/eid1601.090650.
Nahm MH, Lin J, Finkelstein JA, Pelton SI: Increase in the prevalence of the newly discovered pneumococcal serotype 6C in the nasopharynx after introduction of pneumococcal conjugate vaccine. J Infect Dis. 2009, 199: 320-325. 10.1086/596064.
Park IH, Moore MR, Treanor JJ, Pelton SI, Pilishvili T, Beall B, et al: Differential effects of pneumococcal vaccines against serotypes 6A and 6C. J Infect Dis. 2008, 198: 1818-1822. 10.1086/593339.
Brueggemann AB, Griffiths DT, Meats E, Peto T, Crook DW, Spratt BG: Clonal relationships between invasive and carriage Streptococcus pneumoniae and serotype- and clone-specific differences in invasive disease potential. J Infect Dis. 2003, 187: 1424-1432. 10.1086/374624.
Hill PC, Cheung YB, Akisanya A, Sankareh K, Lahai G, Greenwood BM, et al: Nasopharyngeal carriage of Streptococcus pneumoniae in Gambian infants: a longitudinal study. Clin Infect Dis. 2008, 46: 807-814. 10.1086/528688.
Sleeman KL, Griffiths D, Shackley F, Diggle L, Gupta S, Maiden MC, et al: Capsular serotype-specific attack rates and duration of carriage of Streptococcus pneumoniae in a population of children. J Infect Dis. 2006, 194: 682-688. 10.1086/505710.
Goldblatt D, Hussain M, Andrews N, Ashton L, Virta C, Melegaro A, et al: Antibody Responses to Nasopharyngeal Carriage of Streptococcus pneumoniae in Adults: A Longitudinal Household Study. J Infect Dis. 2005, 192: 387-393. 10.1086/431524.
Weinberger DM, Dagan R, Givon-Lavi N, Regev-Yochay G, Malley R, Lipsitch M: Epidemiologic evidence for serotype-specific acquired immunity to pneumococcal carriage. J Infect Dis. 2008, 197: 1511-1518. 10.1086/587941.
Granat SM, Ollgren J, Herva E, Mia Z, Auranen K, Makela PH: Epidemiological evidence for serotype-independent acquired immunity to pneumococcal carriage. J Infect Dis. 2009, 200: 99-106. 10.1086/599364.
Malley R: Antibody and cell-mediated immunity to Streptococcus pneumoniae: implications for vaccine development. J Mol Med. 2010, 88: 135-142. 10.1007/s00109-009-0579-4.
Eskola J, Kilpi T, Palmu A, Jokinen J, Haapakoski J, Herva E, et al: Efficacy of a pneumococcal conjugate vaccine against acute otitis media. N Engl J Med. 2001, 344: 403-409. 10.1056/NEJM200102083440602.
Lucero MG, Nohynek H, Williams G, Tallo V, Simoes EA, Lupisan S, et al: Efficacy of an 11-valent pneumococcal conjugate vaccine against radiologically confirmed pneumonia among children less than 2 years of age in the Philippines: a randomized, double-blind, placebo-controlled trial. Pediatr Infect Dis J. 2009, 28: 455-462. 10.1097/INF.0b013e31819637af.
Weycker D, Richardson E, Oster G: Childhood vaccination against pneumococcal otitis media and pneumonia: an analysis of benefits and costs. Am J Manag Care. 2000, 6: S526-S535.
Trotter CL, Ramsay ME, Slack MP: Rising incidence of Haemophilus influenzae type b disease in England and Wales indicates a need for a second catch-up vaccination campaign. Commun Dis Public Health. 2003, 6: 55-58.
McVernon J, Trotter CL, Slack MP, Ramsay ME: Trends in Haemophilus influenzae type b infections in adults in England and Wales: surveillance study. BMJ. 2004, 329: 655-658. 10.1136/bmj.329.7467.655.
Trotter CL, Andrews NJ, Kaczmarski EB, Miller E, Ramsay ME: Effectiveness of meningococcal serogroup C conjugate vaccine 4 years after introduction. Lancet. 2004, 364: 365-367. 10.1016/S0140-6736(04)16725-1.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/10/90/prepub
We thank Matthew Moore for his insightful comments on this manuscript, Cynthia Whitney and Tamar Pilishvili for supplying the IPD surveillance data in USA, CDC, and Steve Black for valuable discussions. We also thank Pauline Kaye and Usha Gungabissoon for providing the IPD surveillance data in England and Wales. Alessia Melegaro was funded with the POLYMOD grant. Yoon Hong Choi was funded by the Department of Health, UK.
The authors declare that they have no competing interests.
AM conceived the work, developed the programme and produced the initial draft of the manuscript. NJG contributed to the implementation of the model. YHC extended the programme, contributed to the drafting of the manuscript and performed final simulations. WJE and EM contributed to the interpretation of the results. RG provided the carriage data and contributed to the interpretation of the results. All the authors have read and agreed on the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.