Modeling the effects of strain diversity and mechanisms of strain competition on the potential performance of new tuberculosis vaccines

October 21, 2008
105 (42) 16302-16307

Abstract

While bacillus Calmette–Guérin vaccination plays an important role in reducing the morbidity of tuberculosis (TB) infection during childhood, new tuberculosis vaccines are necessary to disrupt the transmission of disease and improve global control of this pathogen. Growing evidence of the presence of meaningful Mycobacterium tuberculosis strain diversity, coupled with the possibility that new vaccines may differentially protect against infection or disease with circulating M. tuberculosis strains, suggest that these vaccines may have complicated effects on disease dynamics. We use a mathematical model to explore the potential effects of strain diversity on the performance of vaccines and find that vaccines offer great promise for improving tuberculosis control, but the expected benefits of mass vaccination will be eroded if strain replacement with M. tuberculosis variants that are not effectively targeted by vaccines occurs. Determining the likelihood of strain replacement will require additional knowledge of the strain specificities of current vaccine candidates, and an improved understanding of the mechanisms of strain interaction, which are responsible for maintaining the diversity of M. tuberculosis within communities.
Despite the availability of effective antibiotics and the widespread use of the Bacille Calmette-Guérin vaccine (bacillus Calmette–Guérin), tuberculosis (TB) remains a major cause of morbidity and mortality globally, leading to nearly nine million new cases per year and two million deaths. Although the bacillus Calmette–Guérin vaccine has been shown to protect against childhood forms of TB, it is far less effective in preventing the adult pulmonary forms of disease that lead to respiratory transmission of the infection. Consequently, a major priority of current global research efforts is the development of new, more effective anti-TB vaccines, capable of having a substantial impact on TB control. To date, this effort has yielded >200 new vaccine candidates, many of which are now undergoing animal testing and early clinical trials; these include recombinant bacillus Calmette–Guérin and other attenuated live mycobacterial vaccines, subunit vaccines, and DNA vaccines (12).
Vaccine discovery efforts have largely focused on the identification of M. tuberculosis (Mtb) antigens that are capable of eliciting effective human T cell responses. Because previous studies had suggested that there was little genetic diversity within Mtb in general and within possible protein targets of host immune surveillance in particular (3, 4), the vaccine development effort has not prioritized issues associated with potential strain variation in immune targets. However, tools for the large-scale genetic typing of mycobacterial strains have generated substantial evidence of important genetic differences among clinical isolates (5). Subsequent studies have revealed a global population structure of M. tuberculosis with several distinct lineages (6, 7) and have provided evidence of heterogeneity in virulence and transmissibility (810), host-specificity (11, 12), and potential for immune system activation (1318). Taken together, these data raise the possibility that antigens identified from a limited set of Mtb strains may be differentially present in clinical isolates, and thus vaccines directed at these antigens may protect less well against some strains of mycobacteria than others.
The hypothesis that a vaccine may be less effective against specific strains of TB was first offered as a possible explanation for the observation that the Beijing family of strains predominated in several countries in East Asia where bacillus Calmette–Guérin coverage was widespread (19). While an initial epidemiological investigation did not find an association between bacillus Calmette–Guérin vaccination status and risk of TB with a Beijing-type strain (20), other work suggests that bacillus Calmette–Guérin vaccination may be associated with disease because of a subgroup of typical Beijing strains (21). The hypothesis that the emergence of the Beijing family of strains may be related to bacillus Calmette–Guérin has continued to circulate (22) and is supported by experiments that have found that bacillus Calmette–Guérin protects relatively poorly against disease because of Beijing-type strains, at least in animal models (9, 23, 24).
Recent evidence also suggests that current subunit vaccine candidates may target antigens that are not uniformly present among Mtb. Hebert et al. investigated the genomic diversity of the pepA and PPE18 genes that code for the protein products composing the Mtb72F subunit vaccine currently under clinical investigation (25). Using Mtb isolates from two different geographical areas, these investigators identified single nucleotide polymorphisms and insertions/deletions within regions reported to be T cell epitopes in a modest proportion of isolates for both sites and consequently suggested that this subunit vaccine may not induce protective immunity to a fraction of clinical Mtb strains.
Strain variation in immune targets could have wide-reaching implication for TB vaccine design and efficacy, because the selective pressure imposed by a widely used but strain-specific TB vaccine might be tremendous (12). Such selection has been demonstrated for several other bacterial pathogens including Streptococcus pneumoniae, Haemophilus influenzae, Neisseria meningitides, Bordetella pertussis [reviewed in ref. 26]. In these cases, vaccination programs have resulted in increases in the prevalence of nonvaccine-type strains after the introduction of vaccination (i.e., strain replacement). Theoretical work also suggests that strain replacement can occur not only because of the differential effectiveness of a vaccine, but also through the use of a pan-effective vaccine when its presence causes additional strain interactions that favor the less-prevalent strains (27). The accumulating evidence of strain variability of Mtb, and the existence of multistrain infection (i.e., coinfection) (28) and superinfection (29), suggest that the population dynamic effects of a vaccine program may be complex.
Several examples of such strain variation in immunogenicity have been detected, and multiple investigators have noted differences in immunomodulatory cytokine production in cell and animal models of infection after exposure to different strains of Mtb. Most notably, infection with the CDC1551 strain appears to trigger a greater protective Th1 immune response than infection with Beijing HN878 (3032), a response that may account for the apparent differences in virulence between the two strains. In a study among TB patients in West Africa, T cell responses to the mycobacterial antigen ESAT-6 were attenuated in sera of those infected with M. africanum compared to those infected with other M. tuberculosis strains (33). Other researchers have detected mycobacterial lipids that are differentially present among Mtb strains and that appear to have a profound effect on the efficacy of the host innate immune response (14, 15). These examples suggest that the immune response provoked by some Mtb strains may be specific to antigenic determinants that are heterogeneously produced within the species.
Neither the specific mechanisms by which new TB vaccines will provide immunity nor the details of how distinct strains of tuberculosis interact and compete are currently well understood. At this time, evidence suggests that distinct strains of Mtb exist and may invoke different immune responses and that some vaccine candidates (especially subunit vaccines) may target only a fraction of clinical strains. To examine the potential effects of strain diversity and vaccine specificity—a combination that has resulted in dramatic examples of strain replacement for other infectious pathogens—we developed simple mathematical models of competing strains of Mtb to study the potential impact of new TB vaccines on (i) the total prevalence of TB and (ii) the relative abundance of strain types. We model vaccines administered before (preexposure) and after infection (postexposure); in both cases we model vaccines that preferentially target the most prevalent strain before the introduction of the vaccine. These models allow us (i) to describe a range of possible qualitative behaviors of the system given our uncertainty about how Mtb strains compete for hosts and how new vaccines will function, and (ii) to identify important unanswered questions about strain interaction and vaccine specificity.

TB Infection Model.

We develop a simple compartmental model of tuberculosis in which two strains are in competition; we arbitrarily designate these as strain 1 and strain 2, with strain 1 being more prevalent in the absence of vaccination. Individuals are born susceptible to infection (X). Upon infection, individuals either become latently infected with probability 1-pi (Li; where i identifies the infecting strain) or progress immediately to infectious TB (Ii) with probability pi. Those in the latent state progress to infectious TB at a rate η. Infectious individuals are removed from the infectious class at rate μT and either die from TB or control their infection. A fraction (r) of infectious individuals who contain but do not sterilize their infection return to a state of slowly progressing latency from which they have a small annual risk of relapse. Individuals are removed from other classes at a much lower rate μ; all removed individuals are replaced by entry of individuals into the susceptible classes representing an assumption of constant population size.
Individuals with a latent infection remain partially susceptible to reinfection. We compare the effect of vaccines in two closely related models: the first model specifies that reinfection by a new strain displaces the old strain (superinfection model, Fig. 1A) while the second model allows that reinfection can result in simultaneous infection with more than one strain type (coinfection model, Fig. 1B). To explore a broad range of potential mechanisms by which distinct strains may activate or subvert host immune responses and to represent differential susceptibility of strains to immunity acquired through infection, we model three mechanisms of within-host strain competition on the population-level impact of vaccination programs:
Fig. 1.
Model structures (A) Superinfection model: Individuals are assigned to unvaccinated (S) or vaccinated (Sv) groups upon entry into population and may be infected with either strain 1 or strain 2. Individuals with latent infection (L) and active disease (I) are indexed by strain type (1 or 2) and vaccination status (vaccinated individuals designated by v). In this version of the model, reinfection of individuals with latent infections can result in strain replacement (dotted arrows). State transitions are as described in the text and by the system of equations presented in the SI Appendix. (B) Coinfection model: In this variant of the model, reinfection of those in the latent state may result in a state of mixed strain coinfection (L12 and L12v).
(1) Individuals who are latently infected with either strain 1 or strain 2 are protected equally well against reinfection with either strain 1 or strain 2 (i.e., self-immunity equals cross-immunity and is symmetric).
(2) Individuals who are latently infected with a particular strain are better protected against reinfection with that same strain than against reinfection with the other strain (i.e., self-immunity is greater than cross-immunity and is symmetric).
(3) Latent infection with strain 1 protects equally well against reinfection with strain 1 or 2, but infection with strain 2 does not protect against reinfection with either strain 1 or 2. That is, strain 1 is immunogenic and primes the immune system effectively to reexposure with any strain, but primary infection with strain 2 does not trigger an effective immune response.
We adopt this agnostic approach toward modeling the means by which strains may either stimulate or be vulnerable to immune responses and these three general mechanisms represent a wide range of strain interactions that are biologically plausible given the available evidence. Mechanism 1 is most likely if observed strain diversity does not translate into functional differences in either the triggering of or the susceptibility to immune responses; mechanism 2 is most likely if each distinct strain triggers an immune response that is most effective against future assault by a similar strain; and mechanism 3 is most likely if there are some strains that do not elicit and may also subvert acquired immune responses—in reality, some combinations of these mechanisms may be at work. Changes to the immunity-related parameter values for these mechanisms of competition are presented in supporting information (SI) Table S1, along with the values and ranges for the other parameters used in simulations. The formula for calculating vaccine impact is provided in the Methods, the descriptions for the model compartments are listed in the legend for Fig. 1, and the differential equations for these models are available in SI Appendix.

Results

We examine the effects of introducing a vaccine into populations with a total median equilibrium TB prevalence of ≈220/100,000; the global TB prevalence was estimated to be 219/100,000 in 2006 (34). At this time in the simulations, strain 2 is the minority strain in the population. The proportion of the total TB prevalence that is strain 2 before vaccination varies with the mechanism of strain interaction assumed (Fig. 2 D–F, dark bars). Under interstrain competition mechanisms 1 (equal and symmetric self- and cross-immunity) and 3 (strain 1 protects against all, strain 2 provides negligible immunity), there is essentially no stable coexistence of strains at equilibrium; hence, we introduce a very small amount of strain 2 into the population at the same time as vaccination, to simulate the sporadic appearance of a vaccine-resistant variant through migration. Under the assumption that individuals in the mixed strain infection states from the coinfection model (L1,2 and L1,2v) can return to the singly infected states upon reinfection (L1, L2, L1v, and L2v), we find that the effects of pre- and postexposure vaccines with equivalent effects on reducing the reproductive number of strain 1 have similar, although not identical, impact on the total equilibrium TB prevalence and on the relative frequency of strain types in both the superinfection and coinfection models. If, on the other hand, the state of mixed infection is relatively protected and thus reinfection is not likely to result in a return to homogenous infection, the coinfection model predicts substantially increased levels of strain diversity. For simplicity, we present only the preexposure vaccine effects in the superinfection model in the main text and figures; the results for the coinfection model (including those in which reinfection of those in the mixed infection states does not result in return to single strain infection states) are available as SI Appendix.
Fig. 2.
Effects of preexposure vaccines in superinfection model. The equilibrium effects of vaccination on total TB prevalence under the three different mechanisms of strain interaction (A, mechanism 1; B, mechanism 2; and C, mechanism 3). These figures show the results of 10,000 simulations with parameters drawn randomly from uniformly distributed ranges listed in Table S1. Red points indicate simulations in which only strain 2 persists after vaccination, blue points indicate simulations in which only strain 1 persists after vaccination, and green points indicate simulations in which both strain 1 and strain 2 persist after vaccination. Vaccine impact is a function of coverage and efficacy and is calculated according to the equation in the text. R2/R1 indicates the relative reproductive number of strain 2 compared to strain 1. Panels D, E, and F show the relative abundance of strain 2 before and after vaccination under the same three mechanisms of strain interaction.
We report the results of 10,000 simulations on parameter sets randomly selected from the uniformly distributed parameter ranges listed in Table S1. We introduce vaccination as an ongoing policy as the TB prevalence approaches equilibrium, and we report the values of TB prevalence and relative strain abundance after the system reaches equilibrium in the presence of vaccination.

Interstrain Competition Mechanism 1.

If latent infection with strain 1 or strain 2 provides equal protection against reinfection with either strain, we find unsurprisingly that, on average, higher total vaccine usage and efficacy is associated with a lower total TB prevalence (Fig. 2A). In scenarios where vaccine coverage and the relative reproductive number of strain 2 are relatively high, we observe an increased probability that strain 2 replaces strain 1. We note that despite this strain replacement, there is still a lower total level of disease after vaccination is implemented, because we enforce that strain 2 always has a lower reproductive number than strain 1 in these simulations. The histogram (Fig. 2D) shows the relative abundance of strain types before and after vaccination and reveals that this mechanism of strain interaction does not support the stable coexistence of the two strain types.

Interstrain Competition Mechanism 2.

If there is strain specificity to the immune response such that latent infection with either strain protects more effectively against reinfection by the same strain than the other strain, we again observe lower mean total equilibrium TB prevalence when vaccination coverage and efficacy are at relatively high levels (Fig. 2B). However, in this scenario where intrastrain competition is stronger than interstrain competition, strain 1 and strain 2 can stably coexist and it is possible that the total TB prevalence can be higher after the introduction of the vaccine (Fig. 2B). Because infection with either strain opposes its own takeover, this mechanism of interstrain competition strongly facilitates strain coexistence (Fig. 2E) and allows that at levels of vaccine impact and relative reproductive number of strain 2 that permit coexistence, use of a vaccine may actually lead to increased disease levels. As vaccine impact is increased to high enough levels that strain 1 is successfully suppressed, this perverse vaccine effect is not observed.

Interstrain Competition Mechanism 3.

If latent infection with strain 1 protects equally well against reinfection with strain 1 or strain 2, but latent infection with strain 2 does not protect against reinfection by either strain, vaccines which specifically target strain 1 may not improve control of TB. This result is not unexpected as under these assumptions vaccination may result in the replacement of an immunogenic strain by a nonimmunogenic one. In particular, when the relative reproductive number of strain 2 is high, implementation of more effective vaccines at higher levels can actually cause a higher equilibrium prevalence of tuberculosis (Fig. 2C), although this perverse effect occurs in only a modest proportion of the simulations in which strain 2 replaces strain 1 (yellow bars, Fig. 2F). When the relative reproductive number of strain 2 is high and there is high vaccine coverage level with a vaccine that can successfully suppress strain 1, the total postvaccination prevalence of TB can exceed the prevaccination levels (Fig. S3 displays the combinations of parameter values that result in this pernicious effect). As with interstrain competition mechanism 1, we observe that stable coexistence of strains is not a generic feature of this mechanism of strain interaction (Fig. 2F).

Effects of Preexposure and Postexposure Vaccination.

Our primary focus is to explore the effects of strain diversity and competition on the effect of vaccines, accordingly, our simulations include only pre- and postexposure vaccines with equivalent ability to reduce the reproductive number of strain 1. This allows us to examine the effect of vaccine mechanism rather than vaccine strength and is not intended to imply that pre- and postexposure vaccines under development are likely to have equal impact on the basic reproductive number of tuberculosis. Other researchers have considered in greater detail the likely benefits of each vaccine approach (3537). While preexposure and postexposure vaccines with equivalent effects on R01 have similar effects on the projected quantity and distribution of tuberculosis strains at equilibrium, they differ in the timing of their impacts on TB epidemics. We find that the postexposure vaccines have a slower impact on disease dynamics than preexposure vaccines with differences persisting for decades (results not shown); this observation is consistent with previous models of TB vaccines.

Discussion

In this study, we use mathematical models to examine a range of possible long-term effects of large-scale vaccination programs with new vaccines against M. tuberculosis. For simplicity, we model the diversity of Mtb using only two strains, and we allow these strains to compete against each other and to be differentially suppressed by vaccination. We examine the sensitivity of the models to combinations of parameter values and explore a range of scenarios in which the mechanisms of competition between strains were allowed to vary. In general, these models suggest that the public health impact of vaccination programs with new TB vaccines will depend crucially on (i) the diversity of circulating Mtb strains, (ii) the mechanisms of competition between strains, and (iii) the strain specificity of these vaccines.
As information about the immunogenicity of strains, the specificity and magnitude of immune responses to previous infection, and the strain specificity of vaccine candidates is limited, we present simple conceptual models that cover a broad range of biologically plausible scenarios. While this approach accurately represents the lack of relevant data and necessarily limits our ability to make quantitative predictions about the effect of new TB vaccines, it does allow us to describe those behaviors that do not appear to depend on these areas of ignorance and identify important uncertainties that are expected to greatly affect the performance of new vaccines. Encouragingly, we find that increasing vaccine coverage and efficacy results in substantial reductions in the predicted TB burden under most assumptions about the strain interactions and vaccine specificity. However, we also find that the benefits of vaccination are reduced if the preferential removal of one competitor allows a previously outcompeted nonvaccine-type strain to emerge. These observations are consistent with results of previous models for other pathogens and with the observed strain replacement that has occurred after strain-specific vaccines have been introduced. We also demonstrate that if a vaccine preferentially targets a strain that provides cross-immunity against a nonvaccine strain, which itself does not provide substantial immunity, the vaccine may cause an absolute increase in the TB burden. Similar perverse effects of vaccines were previously discussed by both McLean (39) and Lipsitch (40), using models for different diseases.
Although vaccination programs reduced TB prevalence in most of our simulations, they also resulted in a shift in strain composition that can potentially undermine some of the benefits of this effect. For example, several studies have reported a positive association between the Beijing genotype and antituberculosis drug resistance within particular locations (4143). A vaccine that reduces the total number of TB cases but results in an increased frequency of Beijing-type strains may lead to the proliferation of strains that are either drug resistant or have an increased potential to become drug resistant.
The mechanisms by which vaccination can shift the competitive balance between competing pathogen strains have been the focus of other mathematical models. In our two-strain models, we consider the effect of a vaccine that selectively suppresses one strain and thus may potentiate the ascent of the other. This approach is similar to other models that have explored the effects of vaccines on the ecology of pathogen strains (44, 45). While this approach has the advantage of simplicity, there are several limitations to the models we present here. Here we assume that the two strains differ in their ability to infect and cause rapid progression to disease, but specify that the duration of infectivity and the rate of slow progression from latency to disease are similar for both strains. We choose to allow the relative infectivity and fraction of infections that lead rapidly to disease to be the traits by which strains differ because, on the basis of other TB models (44, 45), these are characteristics that may have strong impact on disease dynamics. Differences in the expected duration of illness or the rate of endogenous reactivation would likely have resulted in different quantitative results, but similar qualitative conclusions as those presented here. While it is not possible to evaluate the sensitivity of the model to all possible structural assumptions, we have attempted to provide several simple alternatives here and in the SI. While the likelihood of strain coexistence and replacement does show some sensitivity to the structural assumptions used, the central message that M. tuberculosis strain diversity can erode TB vaccine performance is unchanged.
An additional limitation of our approach is that it does not simulate strain evolution, which will occur over the long time scales in which these dynamics unfold. Some modelers have focused on the stochastic effects of mutation on strain diversity and the resultant selection of vaccine-resistant variants under the pressure of vaccination programs (46, 47). These approaches require the specification of parameters governing strain mutation and selection for which we have very little relevant data for Mtb. Thus, the qualitative insights from our models (e.g., Under which combinations of vaccine performance and mechanisms of strain interaction might strain replacement occur?) are more informative than the specific quantitative results (e.g., What is the expected prevalence of disease after vaccine is introduced? How many years after vaccination starts would we expect to see the emergence of nonvaccine-type strains?). We note that any process by which strains mutate and are selected under the pressure of vaccination is likely to erode the benefits of vaccination projected by our simple models.
Previous mathematical models have been used to estimate the impact of new vaccines on the control of TB epidemics and to compare the relative effectiveness of preexposure and postexposure vaccination on the burden of disease (3538, 48). Our models are structured to explore the effects of vaccination on both the total levels of disease and the relative abundance of strain types under various assumptions about between- and within-host competition between strains. While our results largely support previous findings that new vaccines promise to improve disease control, we also find that the shifting population dynamics of Mtb in response to vaccination can limit these benefits.
These simple models identify several important loci of ignorance that may be addressed in future studies. First, larger studies of geographically diverse clinical strains of tuberculosis, as suggested by Hebert et al. (25), would help to clarify to what extent candidate vaccine antigens may be differentially present in clinical strains. If conserved antigens (or combinations of antigens), which are essential for mycobacterial viability and which are capable of eliciting strong T cell responses, can be identified and used to develop new vaccines, the risk of strain replacement would be minimized. Although subunit vaccines targeted at selected antigens are most likely to exert selective pressures, it is also possible that recombinant bacillus Calmette–Guérin or other live attenuated mycobacterial vaccines can have this effect. Limited epidemiological and laboratory evidence consistent with the hypothesis that bacillus Calmette–Guérin vaccination may not offer equal protection against Beijing-type strains suggests that further study of the potential selectivity of these types of vaccines is also warranted. Studies in which vaccinated animals are challenged with diverse Mtb strains (at minimum, a virulent Beijing-type strain in addition to the usual H37Rv or Erdman strain) can shed light on the potential selectivity of vaccine candidates.
Furthermore, these models reveal that the mechanisms by which Mtb strains compete also will affect the performance of vaccines; while molecular strain typing methods allow the detection of reinfection and multiple strain infection events, very little is yet known about how strains differ in their ability to stimulate effective immune responses and how specific these immune responses are. Determining the likelihood of strain replacement by nonvaccine-type strains requires a better understanding of the degree to which there is strain-specific induction and vulnerability to host immune responses and a clearer picture of which other mechanisms are responsible for maintaining the diversity of Mtb within communities.

Methods

Strain Differences.

We assume that strain 2 suffers reproductive deficits compared with strain 1; we specify that strain 2 is less transmissible than strain 1 (β2 < β1) and that a smaller proportion of those infected with strain 2 progress immediately to disease (p2 < p1). For simplicity, we assume that the rate of endogenous reactivation from latency and the duration of infectivity associated with the two strains does not differ for the two strains. The differences between the individual reproductive capacities of the two strains can be calculated by comparing the basic reproduction number for each strain (R0), which is defined as the expected number of secondary cases of infectious disease produced by a single infectious individual entering a completely susceptible population. While models that allow for superinfection and coinfection can allow for a strain with a lower reproductive number to outcompete a strain with a higher reproductive number (49), we do not observe this effect in the ranges of the parameters specified in our models and before the introduction of vaccination, the prevalence of strain 1 is highest. The reproductive number for each strain in this model is defined as:
where i indexes strain type. The first term in this expression is the contribution from those who proceed immediately to disease upon infection and the second is the contribution from those who progress slowly from latency.
When choosing randomized parameters for our simulations, we constrain R01 to be equal to two but allow the parameters p1 and β1 to vary. For each simulation we also randomly select a parameter (z) for strain 2 between 0.8 and 1.0; z is the scalar by which we multiply the parameters p1 and β1 to set the values of p2 and β2. In these models, for values of z below this range, we do not observe emergence of strain 2. We report the relative reproductive capacity of strain 2 compared with strain 1 for each simulation by calculating the ratio R02/R01.

Vaccination.

Here we consider vaccines that protect individuals from infection and/or disease with the most abundant strain (strain 1) and do not have any direct effect on the minority strain (strain 2). The preexposure vaccine reduces the risk of infection, reduces the probability of immediate progression to disease upon infection, and reduces the rate of progression from latency to disease, each by an equal fraction. The postexposure vaccine operates solely by reducing the rate of progression from latency to disease. To simulate the effect of a postexposure vaccine of comparable strength to the preexposure vaccine, we calibrate the reduction of the progression rate necessary to achieve an equivalent impact on the basic reproductive rate of strain 1 in the presence of vaccine. The population level total vaccine impact (VI) can be quantified as
where c is the fraction of the population receiving vaccine, where R01v is the reproductive number for strain 1 in a population where everyone is vaccinated and R01 is the reproductive number for strain 1 in a population where nobody is vaccinated. We assume that vaccination status has no impact on the infectiousness of individuals with active disease and that there is no waning of the vaccine effect.

Acknowledgments.

We thank Gary Schoolnik and Barun Mathema for comments on an early draft and Hana Akselrod for assistance with manuscript preparation. This work was supported by National Institutes of Health Grant K08 AI055985–06 (T.C.).

Supporting Information

Supporting Appendix (PDF)
Supporting Information
Supporting Information (PDF)
Supporting Information

References

1
UD Gupta, VM Katoch, DN McMurray, Current status of TB vaccines. Vaccine 25, 3742–3751 (2007).
2
SH Kaufmann, Envisioning future strategies for vaccination against tuberculosis. Nat Rev Immunol 6, 699–704 (2006).
3
JM Musser, A Amin, S Ramaswamy, Negligible genetic diversity of mycobacterium tuberculosis host immune system protein targets: Evidence of limited selective pressure. Genetics 155, 7–16 (2000).
4
S Sreevatsan, et al., Restricted structural gene polymorphism in the Mycobacterium tuberculosis complex indicates evolutionarily recent global dissemination. Proc Natl Acad Sci USA 94, 9869–9874 (1997).
5
RD Fleischmann, et al., Whole-genome comparison of Mycobacterium tuberculosis clinical and laboratory strains. J Bacteriol 19, 5479–5490 (2002).
6
MM Gutacker, et al., Single-nucleotide polymorphism-based population genetic analysis of Mycobacterium tuberculosis strains from 4 geographic sites. J Infect Dis 193, 121–128 (2006).
7
I Filliol, et al., Global phylogeny of Mycobacterium tuberculosis based on single nucleotide polymorphism (SNP) analysis: Insights into tuberculosis evolution, phylogenetic accuracy of other DNA fingerprinting systems, and recommendations for a minimal standard SNP set. J Bacteriol 188, 759–772 (2006).
8
G Thwaites, et al., Relationship between Mycobacterium tuberculosis genotype and the clinical phenotype of pulmonary and meningeal tuberculosis. J Clin Microbiol 46, 1363–1368 (2008).
9
B López, et al., A marked difference in pathogenesis and immune response induced by different Mycobacterium tuberculosis genotypes. Clin Exp Immunol 133, 30–37 (2003).
10
SE Valway, et al., An outbreak involving extensive transmission of a virulent strain of Mycobacterium tuberculosis. N Engl J Med 338, 633–639 (1998).
11
M Hanekom, et al., Evidence that the spread of Mycobacterium tuberculosis strains with the Beijing genotype is human population dependent. J Clin Microbiol 45, 2263–2266 (2007).
12
S Gagneux, PM Small, Global phylogeography of Mycobacterium tuberculosis and implications for tuberculosis product development. Lancet Infect Dis 7, 328–337 (2007).
13
R van Crevel, et al., Mycobacterium tuberculosis Beijing genotype strains associated with febrile response to treatment. Emerg Infect Dis 7, 880–883 (2001).
14
M Mendelson, S Walters, I Smith, G Kaplan, Strain-specific mycobacterial lipids and the stimulation of protective immunity to tuberculosis. Tuberculosis 85, 407–413 (2005).
15
MB Reed, et al., A glycolipid of hypervirulent tuberculosis strains that inhibits the innate immune response. Nature 431, 84–87 (2004).
16
C Pheiffer, JC Betts, HR Flynn, PT Lukey, P van Helden, Protein expression by a Beijing strain differs from that of another clinical isolate and Mycobacterium tuberculosis H37Rv. Microbiology 151, 1139–1150 (2005).
17
J Dormans, et al., Correlation of virulence, lung pathology, bacterial load and delayed type hypersensitivity responses after infection with different Mycobacterium tuberculosis genotypes in a BALB/c mouse model. Clin Exp Immunol 137, 460–468 (2004).
18
L Rindi, et al., Variation of the expression of Mycobacterium tuberculosis ppe44 gene among clinical isolates. FEMS Immunol Med Microbiol 51, 381–387 (2007).
19
D van Soolingen, et al., Predominance of a single genotype of Mycobacterium tuberculosis in countries of east Asia. J Clin Microbiol 33, 3234–3238 (1995).
20
DD Anh, et al., Mycobacterium tuberculosis Beijing genotype emerging in Vietnam. Emerg Infect Dis 6, 302–305 (2000).
21
K Kremer (Univer of Amsterdam), pp. 234–247, PhD thesis. (2005).
22
F Abebe, G Bjune, The emergence of Beijing family genotypes of Mycobacterium tuberculosis and low-level protection by bacille Calmette-Guérin (BCG) vaccines: Is there a link? Clin Exp Immunol 145, 389–397 (2006).
23
L Tsenova, et al., BCG vaccination confers poor protection against M. tuberculosis HN878-induced central nervous system disease. Vaccine 25, 5126–5132 (2007).
24
L Grode, et al., Increased vaccine efficacy against tuberculosis of recombinant Mycobacterium bovis bacille Calmette-Guérin mutants that secrete listeriolysin. J Clin Invest 115, 2472–2479 (2005).
25
AM Hebert, et al., DNA polymorphisms in the pepA and PPE18 genes among clinical strains of Mycobacterium tuberculosis: Implications for vaccine efficacy. Infect Immun 75, 5798–5805 (2007).
26
M Martcheva, BM Bolker, RD Holt, Vaccine-induced pathogen strain replacement: What are the mechanisms? J R Soc Interface 5, 3–13 (2008).
27
M Iannelli, M Martcheva, X-Z Li, Strain replacement in an epidemic model with super-infection and perfect vaccination. Math Biosci 195, 23–46 (2005).
28
RM Warren, et al., Patients with active tuberculosis often have different strains in the same sputum specimen. Am J Respir Crit Care Med 169, 610–614 (2004).
29
C-Y Chiang, LW Riley, Exogenous reinfection in tuberculosis. Lancet Infect Dis 5, 629–636 (2005).
30
C Manca, et al., Mycobacterium tuberculosis CDC1551 induces a more vigorous host response in vivo and in vitro, but is not more virulent than other clinical isolates. J Immunol 162, 6740–6746 (1999).
31
C Manca, et al., Virulence of a Mycobacterium tuberculosis clinical isolate in mice is determined by failure to induce Th1 type immunity and is associated with induction of IFN-alpha/beta. Proc Natl Acad Sci USA 98, 5752–5757 (2001).
32
C Manca, et al., Differential monocyte activation underlies strain-specific Mycobacterium tuberculosis pathogenesis. Infect Immun 72, 5511–5514 (2004).
33
BC de Jong, et al., Mycobacterium africanum elicits an attenuated T cell response to early secreted antigenic target, 6 kDa, in patients with tuberculosis and their household contacts. J Infect Dis 193, 1279–1286 (2006).
34
; World Health Organization Global Tuberculosis Control:Surveillance, Planning, Financing (World Health Organization, Geneva, WHO/HTM/TB/2008.393. (2008).
35
T Lietman, SM Blower, Potential impact of tuberculosis vaccines as epidemic control agents. Clin Infect Dis 30, S316–S322 (2000).
36
E Ziv, CL Daley, S Blower, Potential public health impact of new tuberculosis vaccines. Emerg Infect Dis 10, 1529–1535 (2004).
37
D Young, C Dye, The development and impact of tuberculosis vaccines. Cell 124, 683–687 (2006).
38
CP Bhunu, W Garira, Z Mukandavire, G Magombedzea, Modelling the effects of pre-exposure and post-exposure vaccines in tuberculosis control. J Theor Biol 254, 633–649 (2008).
39
AR McLean, Vaccination, evolution and changes in the efficacy of vaccines: A theoretical framework. Proc Biol Sci 261, 389–393 (1995).
40
M Lipsitch, Vaccination against colonizing bacteria with multiple serotypes. Proc Natl Acad Sci USA 94, 6571–6576 (1997).
41
F Drobniewski, et al., Rifampin- and multidrug-resistant tuberculosis in Russian civilians and prison inmates: Dominance of the Beijing strain family. Emerg Infect Dis 8, 1320–1326 (2002).
42
HS Cox, et al., The Beijing genotype and drug resistant tuberculosis in the Aral Sea region of Central Asia. Respir Res 6, 134 (2005).
43
, Beijing/W genotype Mycobacterium tuberculosis and drug resistance. Emerg Infect Dis 12, 736–743 (2006).
44
MA Sanchez, SM Blower, Uncertainty and sensitivity analysis of the basic reproductive rate. Tuberculosis as an example. Am J Epidemiol 145, 1127–1137 (1997).
45
C Dye, MA Espinal, Will tuberculosis become resistant to all antibiotics? Proc Biol Sci 268, 45–52 (2001).
46
S Gandon, MJ Mackinnon, S Nee, AF Read, Imperfect vaccination: Some epidemiological and evolutionary consequences. Nature 414, 751–756 (2001).
47
VV Ganusov, R Antia, Imperfect vaccines and the evolution of pathogens causing acute infections in vertebrates. Evolution Int J Org Evolution 60, 957–969 (2006).
48
C Dye, BG Williams, Eliminating human tuberculosis in the twenty-first century. J R Soc Interface 23, 653–662 (2007).
49
MA Nowak, RM May, Superinfection and the evolution of parasite virulence. Proc Roy Soc B 255, 81–89 (1994).

Information & Authors

Information

Published in

The cover image for PNAS Vol.105; No.42
Proceedings of the National Academy of Sciences
Vol. 105 | No. 42
October 21, 2008
PubMed: 18849476

Classifications

Submission history

Received: June 26, 2008
Published online: October 21, 2008
Published in issue: October 21, 2008

Keywords

  1. bacillus Calmette–Guérin
  2. mathematical model
  3. strain replacement

Acknowledgments

We thank Gary Schoolnik and Barun Mathema for comments on an early draft and Hana Akselrod for assistance with manuscript preparation. This work was supported by National Institutes of Health Grant K08 AI055985–06 (T.C.).

Notes

This article contains supporting information online at www.pnas.org/cgi/content/full/0808746105/DCSupplemental.

Authors

Affiliations

Division of Global Health Equity, Brigham and Women's Hospital, Boston MA 02115;
Department of Epidemiology, Harvard School of Public Health, Boston MA 02115;
Caroline Colijn
Department of Epidemiology, Harvard School of Public Health, Boston MA 02115;
Department of Engineering Mathematics, University of Bristol, Bristol BS8 1TR, United Kingdom; and
Megan Murray
Division of Global Health Equity, Brigham and Women's Hospital, Boston MA 02115;
Department of Epidemiology, Harvard School of Public Health, Boston MA 02115;
Division of Infectious Diseases, Massachusetts General Hospital, Boston MA 02114

Notes

To whom correspondence should be addressed. E-mail: [email protected]
Communicated by Barry R. Bloom, Harvard School of Public Health, Boston, MA, September 3, 2008
Author contributions: T.C. and M.M. designed research; T.C. and C.C. performed research; T.C. and C.C. analyzed data; and T.C., C.C., and M.M. wrote the paper.

Competing Interests

The authors declare no conflict of interest.

Metrics & Citations

Metrics

Note: The article usage is presented with a three- to four-day delay and will update daily once available. Due to ths delay, usage data will not appear immediately following publication. Citation information is sourced from Crossref Cited-by service.


Altmetrics




Citations

Export the article citation data by selecting a format from the list below and clicking Export.

Cited by

    Loading...

    View Options

    View options

    PDF format

    Download this article as a PDF file

    DOWNLOAD PDF

    Login options

    Check if you have access through your login credentials or your institution to get full access on this article.

    Personal login Institutional Login

    Recommend to a librarian

    Recommend PNAS to a Librarian

    Purchase options

    Purchase this article to access the full text.

    Single Article Purchase

    Modeling the effects of strain diversity and mechanisms of strain competition on the potential performance of new tuberculosis vaccines
    Proceedings of the National Academy of Sciences
    • Vol. 105
    • No. 42
    • pp. 16059-16407

    Figures

    Tables

    Media

    Share

    Share

    Share article link

    Share on social media