Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, June 2004, p. 2206-2213, Vol. 48, No. 6
0066-4804/04/$08.00+0 DOI: 10.1128/AAC.48.6.2206-2213.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
INSERM U444, Epidémiologie et Sciences de l'Information, 75571 Paris Cedex 12,1 Institut Pasteur,2 Assistance Publique, Hôpitaux de Paris, Université Pierre et Marie Curie, Paris, France3
Received 24 September 2003/ Returned for modification 15 December 2003/ Accepted 10 February 2004
|
|
|---|
|
|
|---|
A new conjugate pneumococcal vaccine that protects against both invasive disease and pharyngeal colonization has been used in the United States since early 2000 and is now recommended for young children in other countries. Recent observations in the United States have suggested that this vaccine is highly effective in preventing invasive disease in infants and young children (29). Clinical trials of conjugate vaccines have demonstrated their ability to reduce colonization by vaccine-related serotypes among vaccinated individuals (6). Because asymptomatic carriers of S. pneumoniae, especially children, are the main reservoir for transmission, a decrease in colonization by vaccine-related serotypes in unvaccinated individuals is expected as well, through herd immunity (22). Furthermore, as the vaccine formulation covers most pneumococcal serotypes currently associated with antibiotic resistance, it has also been suggested that conjugate vaccination could play a role in controlling resistance (22), and preliminary studies have shown a reduction in the carriage of resistant strains after vaccination (5, 16).
However, the serotypes of the most frequent and the most resistant strains of S. pneumoniae circulating in some European countries only partially match those used in the vaccine formulation (27). Many are concerned that the differences between the American and European serotype distributions may impede vaccine efficacy in Europe (25, 28).
Given that pneumococcal conjugate vaccines have been used in the general populations of only a few countries and for a short period, it is too early to assess their long-term effects on carriage and on resistance selection from field data. Mathematical modeling provides a useful tool for exploring the impact of vaccination. Major consequences may be investigated in this way. For example, a model suggested that because the vaccine does not protect against all existing bacterial serotypes, an increase in carriage of nonvaccine serotypes was likely (18). This is consistent with recent data suggesting that colonization with nonvaccine serotypes of S. pneumoniae has been increasing in the United States since the introduction of the vaccine and that nonsusceptible nonvaccine serogroups, such as serogroup 35, have emerged (S. Pelton, C. Marchant, D. Christiansen, and A. Loughlin, Abstr. 43rd Intersci. Conf. Antimicrob. Agents Chemother., abstr. 890, 2003).
Here, we focus on the impact of vaccination on antibiotic resistance. We have developed a mathematical model of selection for S. pneumoniae resistance to penicillin G in an age-structured population. In this paper, we study the effect of a vaccination program carried out on children <2 years old with the currently available seven-valent conjugate vaccine and with a hypothetical "optimized" vaccine designed to incorporate most resistant serotypes. More precisely, we investigated two questions: (i) what are the epidemiological characteristics of S. pneumoniae colonization in a vaccinated population, and (ii) how does the distribution of resistance levels among children and adult carriers change after the introduction of vaccination?
(This work was presented, in part, at the 13th European Conference on Clinical Microbiology and Infectious Diseases, Glasgow, 2003, abstract O72.)
|
|
|---|
2,000 strains are typed and evaluated for susceptibility to various antibiotics. We used these data to compute global characteristics for vaccine-type and non-vaccine-type S. pneumoniae. More precisely, we computed global MIC distributions and colonization rates for vaccine-type (or non-vaccine-type) pneumococci in the French population according to age by averaging the corresponding serotype-specific data over all serotypes included (or not included) in the vaccine weighted by the frequency of each serotype.
Vaccines. The currently available seven-valent conjugate vaccine includes serotypes 4, 6B, 9V, 14, 18C, 19F, and 23F. We also imagined an optimized vaccine that would include in its formulation all serotypes currently exhibiting reduced susceptibility to penicillin G (>15% of MICs over 0.1 µg/ml). This amounted to adding four serotypes (serotypes 6A, 15, 19A, and 24) to the seven-valent vaccine (10, 27).
For both these vaccines, we called S. pneumoniae strains with one of the serotypes included in the vaccine type V pneumococci and S. pneumoniae strains with a serotype not included in the vaccine type NV pneumococci.
Model. This work builds on a previously published model of selection for penicillin-resistant S. pneumoniae in the community (26).
The model structure is shown in Fig. 1. Hosts enter the population at birth as noncarriers at a constant rate, µN.
![]() View larger version (41K): [in a new window] |
FIG. 1. Structure of the model. Xnp, number of untreated noncolonized individuals; Ynp, number of treated noncolonized individuals; XV,1, XV,2, etc. (or YV,1, YV,2, etc.), numbers of untreated (or treated) individuals colonized by vaccine-type pneumococci with one of the 10 resistance levels (see text); XNV,1, XNV,2, etc., and YNV,1, YNV,2, etc., numbers of individuals colonized by non-vaccine-type pneumococci. The model is structured into three age classes. A portion of children <2 years old are vaccinated, in which case they can be untreated (V) or treated (W); vaccinated children can be colonized only with non-vaccine-type pneumococci.
|
2 years old), older children (2 to 15 years old), and adults (
15 years old). To incorporate vaccination, the model assumes that a fraction, v, of children <2 years old are vaccinated each year, in which case they are protected against carriage for an average time, dv.
Susceptible hosts may be colonized with either type V or type NV pneumococci. Vaccinated hosts can be colonized only with type NV pneumococci, and dual colonization is excluded. Colonization occurs following contact with hosts carrying type V or NV pneumococci; infectious contacts are more frequent between children than between adults but do not depend on the serotype of the strains. The rates of infectious contacts between individuals belonging to various age classes were calibrated in order to reproduce currently observed proportions of pneumococcal carriers in France in these age classes (26). In the absence of antibiotic treatment, natural decolonization occurs after a time, 1/
, independently of the pneumococcal serotype.
The resistance levels of all bacteria colonizing one host are represented by a single MIC, which takes 1 of 10 possible values: 0.06, 0.125, 0.25, 0.5, 1, 2, 4, 8, 16, or 32.
Independently of their carrier status, hosts are exposed to antibiotics with a certain frequency:
1 for young children,
2 for older children, and
3 for adults. In vaccinated hosts, this rate is reduced to a value of
4. During antibiotic treatment, contacts with carriers are more liable to lead to colonization if the involved strains are resistant, and bacterial colonization is cleared with a probability of
.
In hosts in whom colonization is not eliminated, bacteria with a mutation toward a higher resistance level may replace the original strains (26). We take into account the possibility of the occurrence of both point mutations, leading to a small increase in the MIC for pneumococci, and the acquisition of genetic material from other bacterial species, leading to higher levels of resistance, with the hypothesis that small increments are more frequent than large increments. To capture this feature, we randomly selected an increase in the MIC according to the following seminormal law when antibiotic exposure failed to clear pneumococcal colonization:
![]() |
|
View this table: [in a new window] |
TABLE 1. Probabilities of increment of the penicillin MIC for pneumococcal strains from one level to another following a genetic event when the colonized host is exposed to antibioticsa
|
days on average. Finally, the mortality rate is µ1 for young children, µ2 for older children, and µ3 for adults.
Parameter values consistent with the French situation were estimated from the literature and are detailed in Table 2, together with the relevant references.
|
View this table: [in a new window] |
TABLE 2. List of model parameters and their values
|
Sensitivity analysis. In order to evaluate the robustness of our results, we performed an uncertainty and sensitivity analysis of the model, using the Latin hypercube sampling-partial rank correlation coefficients technique (3). We conducted this analysis for two outcomes: the overall pneumococcal colonization in the population 10 years after the introduction of vaccination and the proportion of resistant strains 20 years after the introduction of vaccination.
|
|
|---|
![]() View larger version (18K): [in a new window] |
FIG. 2. Proportion of carriers of S. pneumoniae in the general population before (5 years) and after (20 years) vaccination of 30% (A) and 80% (B) of children <2 years old. The thick solid line represents the overall carriage rate, while the thin solid line represents carriage of vaccine serotypes, and the dotted line represents carriage of nonvaccine serotypes.
|
15 years for a 30% vaccination coverage level and after only
10 years for an 80% vaccination coverage level. However, almost simultaneously, carriage of type NV pneumococci increased, although this increase never compensated for the decrease in type V strains. The speed of this serotype replacement was also greater when the vaccination coverage was wide. These two opposite trends led to a sharp decrease in overall pneumococcal carriage shortly after the introduction of vaccination, followed by a renewed increase toward a value less than the prevaccination level. Finally, when the vaccination coverage was wide enough, a slight decrease in pneumococcal colonization could be noted after the original bacterial population had been entirely replaced by type NV bacteria.
Vaccination and resistance. Currently observed and predicted distributions of pneumococcal resistance levels after 20 years of vaccination with the seven-valent conjugate vaccine (with a 30% coverage level) are presented in Fig. 3.
![]() View larger version (16K): [in a new window] |
FIG. 3. Currently observed distributions of penicillin MICs in France (A and B) and predicted distributions after 20 years of vaccination with a seven-valent conjugate vaccine (C and D). (A and C) Solid bars represent the distributions of MICs for vaccine serotypes of S. pneumoniae, and open bars represent those for nonvaccine serotypes. (B and D) Overall distributions of MICs for S. pneumoniae. In panel D, the hatched bars represent predictions with vaccination, whereas the stippled bars show predictions without vaccination.
|
With the present high levels of antibiotic exposure in France, resistant strains tend to be selected over time for all serotypes, irrespective of vaccination. Initially, as vaccine serotypes exhibit more decreased susceptibility than others, high levels of resistance are selected faster for type V pneumococci. However, after a while, this progression is slowed down, causing the distribution of resistance levels in type NV pneumococci to catch up. After 20 years, both distributions are very similar, as mostly highly resistant strains are left (Fig. 3C).
Serotype replacement occurs concurrently with the selection of resistance. The portion of type NV pneumococci increases along with their mean resistance level, and no faster, so that the overall distribution is not affected by vaccination.
After 20 years, the overall distribution presents only high resistance levels, irrespective of vaccination (Fig. 3D).
Vaccine "completion." The currently observed and predicted distributions of pneumococcal resistance levels after 20 years of vaccination with the optimized 11-valent conjugate vaccine (with a 30% coverage level) are presented in Fig. 4.
![]() View larger version (16K): [in a new window] |
FIG. 4. Currently observed distributions of penicillin MICs in France (A and B) and predicted distributions after 20 years of vaccination with an optimized (11-valent) conjugate vaccine (C and D). (A and C) Solid bars represent the distributions of MICs for vaccine serotypes of S. pneumoniae; open bars represent those for nonvaccine serotypes. (B and D) Overall distributions of MICs for S. pneumoniae. In panel D, the hatched bars show predictions with vaccination, whereas the stippled bars show predictions without vaccination.
|
0.06 µg/ml) and a peak at high MICs, thus differing from the distribution centered on high MICs obtained without vaccination. Sensitivity analysis. The sensitivity analysis showed that for predicting the overall pneumococcal colonization rate in a vaccinated community, the most critical parameter was the duration of pneumococcal colonization (positively linked) (Table 3), followed by the rates of infectious contacts, especially between adults and between children and adults (positively linked).
|
View this table: [in a new window] |
TABLE 3. Sensitivity analysis of the modela
|
|
View this table: [in a new window] |
TABLE 4. Sensitivity analysis of the modela
|
|
|
|---|
One might think that the reason this seven-valent conjugate vaccine is unable to slow resistance selection is that some of the strains that are not included in its formulation already exhibit some decreased susceptibility to penicillin G. To investigate this issue, we imagined an optimized 11-valent vaccine including all strains that currently exhibit decreased susceptibility to penicillin. With this optimized vaccine, resistance selection was slowed down, and the distribution of MICs in the midterm was affected by the use of the vaccine. However, selection for resistance was not prevented, and the effect of vaccination on resistance decreased with time (data not shown). This can be explained by an increased risk of emergence and transmission of resistant strains among nonvaccine serotypes due to their increased carriage and exposure to antibiotics.
In a Finnish study of the efficacy of a seven-valent pneumococcal vaccine against otitis media (8), episodes attributed to serotypes that were cross-reactive with those included in the vaccine were shown to decrease significantly following vaccination. Therefore, cross-immunity induced by the vaccine could reduce carriage of serotypes 6A, 9N, 18B, 19A, and 23A, as well as that of the seven serotypes included in its formulation, and put the effect of the vaccine midway between our simulations of the 7-valent vaccine and the optimized 11-valent vaccine.
Serotype replacement following vaccination is a consistent prediction from mathematical models (17) and has been observed in several trials (8, 20). However, in some studies, the decrease in carriage was noted in all serotypes of S. pneumoniae rather than in vaccine serotypes only. This finding is consistent with the simulations presented here, in which we observed a transient decrease in the overall carriage shortly after the introduction of vaccination (Fig. 2).
Finally, recent data on pneumococcal colonization and resistance in the United States since the introduction of conjugate vaccination lend support to our predictions. Indeed, these data show an increase in the proportion of nonvaccine S. pneumoniae over the period 2000 to 2003, as well as an increase in antibiotic resistance and the emergence of nonsusceptible nonvaccine serogroups (35 and NT), in total agreement with our predictions (S. Pelton, C. Marchant, D. Christiansen, and A. Loughlin, Abstr. 43rd Intersci. Conf. Antimicrob. Agents Chemother., abstr. 890, 2003).
Hypotheses.
The French Public Hygiene Committee recently recommended the use of a seven-valent pneumococcal conjugate vaccine for all children <2 years old in day care, which represents
30% of all children in this age group. That is why we studied vaccination coverage levels between 30 and 80%, which is consistent with Advisory Committee on Immunization Practices recommendations in the United States of vaccinating all children aged <23 months (1), as well as with the number of doses sold by the manufacturer (15.5 million in 2001) (29).
We supposed that colonized hosts carried either type V or type NV pneumococci, excluding dual colonization. However, simultaneous colonization with up to six different serotypes of S. pneumoniae has been observed in previous studies, with the consequence that serotype replacement under vaccination was favored and vaccine efficacy was decreased (11, 12, 15).
In this study, we did not consider that nonvaccine serotypes had reduced fitness compared to vaccine serotypes. Indeed, even though there may be true differences among the aptitudes for colonization of different serotypes (for instance, differing durations of carriage), there is no evidence that vaccine serotypes as a whole are more adapted for colonization than others (27). We numerically investigated the impact of a difference between carriage durations in vaccine serotypes and other serotypes; we found that this difference would not change our results on resistance selection. However, the speed of serotype replacement was slightly decreased if nonvaccine serotypes had a shorter carriage duration than vaccine serotypes, implying that our results may underestimate the short-term effect of vaccination on the overall pneumococcal colonization rate. Indeed, with a duration of carriage of nonvaccine serotypes 5% (or 10%) shorter than the duration of carriage of vaccine serotypes, serotype replacement would be fully completed 5 years (or 15 years) later, but the proportion of resistant strains after 20 years would remain the same.
Another hypothesis consistent with currently observed epidemiological differences between vaccine and nonvaccine serotypes is that nonvaccine serotypes exhibit a lower likelihood of resistance acquisition. We investigated this possibility by decreasing probabilities of MIC increments, pij (for j > i), for nonvaccine serotypes. We found that this had no impact on the predictions concerning pneumococcal colonization and little impact on those regarding resistance selection. Indeed, even if we suppose that the probabilities of genetic events for nonvaccine serotypes are half the corresponding probabilities for vaccine serotypes, the proportion of resistant strains is only 5% lower than in our predictions after 10 years of vaccination, and this difference is reduced to 1.7% after 20 years.
Sensitivity analysis. Several of the parameters chosen were derived from direct measurments in the community or in vitro, but others required indirect evaluation. We performed a sensitivity analysis using the Latin hypercube sampling technique on the model to evaluate the impact of these uncertainties on the predictions (3).
All in all, uncertainties in estimating the values of four groups of parameters (the carriage duration,
; contact rates, ßij; and treatment characteristics,
i and
) have the most critical effects on the predictions. The analysis also confirmed that a reasonable range of values for the duration of vaccine immunity, as well as the choice of the constants and of the exponent of m in the probability of nondecolonization following treatment,
(m), have little effect on model outcomes.
Developments. Considering the problem of resistance selection in serotypes not included in the vaccine formulation, a dynamic updating of the vaccine could be studied: nonvaccine serotypes would be tested for high levels of resistance at regular time intervals and, if need be, added to the vaccine formulation.
The model as we built it does not allow an evaluation of this solution, as it features only two types of S. pneumoniae and therefore does not allow the detection of newly resistant serotypes among nonvaccine serotypes. However, our study of the increased efficacy of an optimized vaccine is a first step in demonstrating the possibilities of such policies.
|
|
|---|
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
We also introduced the following terms:
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
Model parameters are as shown in Table 1 and Table 2, and the rate of vaccination is v.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»