New Research In
Physical Sciences
Social Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
 Agricultural Sciences
 Anthropology
 Applied Biological Sciences
 Biochemistry
 Biophysics and Computational Biology
 Cell Biology
 Developmental Biology
 Ecology
 Environmental Sciences
 Evolution
 Genetics
 Immunology and Inflammation
 Medical Sciences
 Microbiology
 Neuroscience
 Pharmacology
 Physiology
 Plant Biology
 Population Biology
 Psychological and Cognitive Sciences
 Sustainability Science
 Systems Biology
The epidemiological fitness cost of drug resistance in Mycobacterium tuberculosis

Edited by Carl Bergstrom, University of Washington, Seattle, WA, and accepted by the Editorial Board June 22, 2009 (received for review March 5, 2009)
Abstract
The emergence of antibiotic resistance in Mycobacterium tuberculosis has raised the concern that pathogen strains that are virtually untreatable may become widespread. The acquisition of resistance to antibiotics results in a longer duration of infection in a host, but this resistance may come at a cost through a decreased transmission rate. This raises the question of whether the overall fitness of drugresistant strains is higher than that of sensitive strains—essential information for predicting the spread of the disease. Here, we directly estimate the transmission cost of drug resistance, the rate at which resistance evolves, and the relative fitness of resistant strains. These estimates are made by using explicit models of the transmission and evolution of sensitive and resistant strains of M. tuberculosis, using approximate Bayesian computation, and molecular epidemiology data from Cuba, Estonia, and Venezuela. We find that the transmission cost of drug resistance relative to sensitivity can be as low as 10%, that resistance evolves at rates of ≈0.0025–0.02 per case per year, and that the overall fitness of resistant strains is comparable with that of sensitive strains. Furthermore, the contribution of transmission to the spread of drug resistance is very high compared with acquired resistance due to treatment failure (up to 99%). Estimating such parameters directly from in vivo data will be critical to understanding and responding to antibiotic resistance. For instance, projections using our estimates suggest that the prevalence of tuberculosis may decline with successful treatment, but the proportion of cases associated with resistance is likely to increase.
 antibiotic resistance
 approximate Bayesian computation
 bacterial evolution
 molecular epidemiology
 stochastic model
Despite the development of a number of effective treatments over the past half century, tuberculosis remains one of the most destructive bacterial infections in humans. A third of the world's population carries the bacterium, and even though it remains latent in most of these cases, it is responsible for ≈1.6 million deaths each year (1). In addition to this persistent death rate, an increasing number of strains that are resistant to one or several of the available treatments have been observed. Indeed, recent years have seen the emergence of extensively drugresistant strains (XDR TB), that are resistant to almost all available antibiotic treatments (2). In this article we address the question of whether the overall reproductive fitness of drugresistant strains is higher than that of sensitive strains −essential information for predicting the spread of the disease (3).
Potentially mitigating the spread of antibiotic resistance is a possible biological cost to fitness associated with drug resistance. Evidence for this cost comes from in vitro competition studies of drugresistant bacteria (4, 5). For instance, in the case of tuberculosis, the relative growth advantage of drugresistant strains from clinical isolates has been shown to be between 0.71 and 1.16 in an in vitro competition assay (6). This observation suggests that a typical resistant strain in the epidemiological context may have at most a mild transmission cost. Furthermore, mutations that compensate for the fitness cost may occur in the wild (6). Therefore, the fitness of drugresistant strains in the epidemiological context could be as high as that of the sensitive strains or even higher (7, 8). Estimates of relative fitness have also been obtained by using odds ratios derived from molecular epidemiological data (9). These estimates indicate a wide range of fitness costs associated with resistance, including the possibility that some resistant strains bear no fitness cost.
With little fitness cost to the resistance phenotype, one expects transmission of resistant strains to be a substantial contributor to the prevalence of drug resistance. This idea is supported by studies showing resistant isolates to be genetically homogeneous in some geographic areas (10, 11). Other studies argue for the importance of acquired resistance given the genetic heterogeneity of resistant isolates in other areas (12), and that serial samples from patients reveal the relatively frequent occurrence of drug resistance mutations (13). The significance of acquired resistance is also supported by the observation that prior treatment is a major risk factor for multidrugresistant tuberculosis in Europe (14). Because it is difficult to determine the source of resistance for any given case, in practice, cases are categorized as “primary” if patients were never treated or treated for <1 month, and as “acquired” if patients were treated for >1 month (15). This is an imperfect classification for understanding the relative contributions of transmission and the de novo evolution of resistance because the acquired category also includes the possibility of transmitted resistant strains. There is a need to quantify these contributions more carefully.
The emergence of drugresistant Mycobacterium tuberculosis has been studied through mathematical modeling. Through these models parameters such as detection rates and treatment success rates have been identified as important epidemiological factors (7, 8, 16–19). Mathematical models describing drug resistance dynamics rely on assumptions about the relative fitness of drugresistant strains, among other factors. Accurate estimates of the underlying parameters are of critical importance to predicting the spread of drug resistance (20).
One approach to measuring the relative fitness of resistant strains is to use molecular epidemiological data that allow classification of isolates into genotypic classes (clusters). The relative fitness of resistant strains compared with that of sensitive strains has thereby been quantified from a comparison of their genetic clustering (9). Such methods are indirect in that the dynamics of tuberculosis transmission, evolution of resistance, and mutation of molecular markers are not explicitly taken into account. By mathematically modeling these stochastic processes simultaneously and applying modern computational Bayesian methods of inference (21), we improve these estimates of the relative fitness. Additionally we estimate the cost to transmission incurred by resistance, the rate of acquisition of drug resistance due to treatment failure, and we quantify the relative contributions of resistance evolution (acquired) versus transmission of resistant strains (primary).
In this study, we estimate these quantities using detailed molecular data from M. tuberculosis strains isolated from Cuba (22), Estonia (10), and Venezuela (23). Inference in a Bayesian framework involves constructing a likelihood function for the observed data by using a stochastic model. The level of realism attained by using the stochastic model we have developed (see Fig. 1) makes it impractical to work with likelihood functions directly. We therefore use a recently proposed form of approximate Bayesian computation (21). This methodology approximates posterior distributions of the parameters by simulating datasets under the model and comparing these with empirical data through a weighting function involving summary statistics.
Results
The posterior point estimates obtained from the approximate Bayesian computation procedure can be summarized as follows. The transmission cost of resistance is ≈0.1 for the data from Cuba and from Estonia and ≈0.7 for the data from Venezuela. Estimates of the rate of evolution of resistance vary from 0.0025 (Cuba) to 0.020 (Venezuela). Table 1 provides point estimates with credibility intervals of these parameters, and Fig. 2 shows marginal posterior distributions.
Because strain fitness depends on both the transmission rate and the duration of infection, the transmission cost of resistance alone does not capture the overall fitness cost. The fitness of resistant strains relative to sensitive strains gives a measure of the overall advantage (or disadvantage) that resistance confers and can be expressed in terms of the model parameters (Eq. 1). Our results, using stochastic models with epidemiological data, provide statistical evidence that the populationlevel fitness of drugresistant tuberculosis could be as high as that of drugsensitive strains. Our estimates of the relative fitness of drugresistant strains vary from ≈0.3 in Venezuela to ≈1.0 in Cuba and Estonia (Table 1 and Fig. 2).
We examined the relative contributions of treatment failure and the direct transmission of resistant strains to the ascendance of drug resistance. Our results suggest that direct transmission, rather than evolution after treatment failure, is the major factor responsible for the prevalence of resistant strains, with >90% of resistant cases in Estonia and Cuba attributable to transmission (Fig. 2D).
Projections were computed by using our model, simulating forwards in time by using the posterior parameter distribution. In these projections both moderate cure rates (based on 50% treatment success) and optimistic cure rates [based on the WHO target of 85% treatment success (1)] lead to a decline in the median total number of tuberculosis cases (Fig. 3 A and B). The proportion of these cases that are associated with drugresistant strains increases (Fig. 3 C and D). The decline in the total cases is sharper under the assumption of higher cure rates, but in each scenario, all cases eventually become associated with resistant strains. Fig. 3 C and D also gives the proportion of simulations that resulted in extinction of tuberculosis. Increased cure rates lead to local eradication with high probability after 60 y.
The methodology used here has heavy computatational demands that constrain the level of realism of the model. However, a crossvalidation of our results is possible because the data from Cuba consist of the same isolates typed with two different markers. The shared epidemiological history of these two datasets should produce estimates that differ only in mutation rate. This is broadly what we observe. Indeed, we find that the mutation rate of IS6110 is 5–10 times higher than that of spoligotypes [see supporting information (SI) Appendix], agreeing with previous results (24).
Discussion
This study reveals relatively high fitnesses of resistant strains in Cuba and Estonia. These high values may reflect an inconsistency in treatment programs in these countries from the 1950s to the 1990s. There may have been long periods of ineffective treatment—indeed, Estonia now has one of the highest rates of multidrug resistance in the world (1). The intermittent presence of drugs and the resulting transmission of resistant strains would have given drugresistant strains collectively more time within untreated hosts. This may have allowed mutations to evolve that compensate for the fitness cost, resulting in higher relative fitness for resistant strains. In contrast, the lower estimate of the fitness of drugresistant strains in Venezuela is consistent with the early introduction of an effective national tuberculosiscontrol program (1936) (1) and a steady decrease in tuberculosis prevalence since 1939 (23). These results support previous evidence that circulating strains have lower fitness cost (6), suggesting that resistant strains in regions with ineffective control strategies have higher fitness.
We found a high proportion of drugresistant cases due to transmission compared with the evolution of resistance due to treatment failure in Estonia and Cuba. For these countries, we confirm the interpretation of (10) and others that transmission of resistant strains is the major contributor to the drug resistance problem. This observation suggests that it may be effective to allocate resources to containing the transmission of resistant strains. The relative contribution of transmission in Venezuela is likely to be lower, although it is still high enough to consider focusing on blocking transmission.
Factors other than drug resistance, including host and bacterial factors, are likely to confer differential fitnesses among strains. Previous work provides a method for detecting differences in transmission rates among different genotypes (25). That study showed, however, that these differences are generally not large and do not occur frequently. Nevertheless a more complex model in conjunction with more detailed data may provide opportunities to further develop our methods to better characterize the dynamics of resistance. For example, stochastic models with multiple states of resistance such as that of (26) can be used to further our understanding of the emergence of multidrugresistant or extensively drugresistant tuberculosis.
Accurate estimates of parameters such as the cost of resistance are important for the design of public health measures and predicting epidemiological outcomes. Future prevalence of drug resistance has been intensely debated (3, 7, 18), and projections depend on knowledge of strain fitnesses. Under the modeling assumptions of Cohen and Murray (8), if the cost to fitness is between 0 and 0.1 (our estimates for Estonia and Cuba), then a peak incidence of ≈400–600 infectious tuberculosis cases per 100,000 will occur ≈170–220 years after the introduction of antibiotics (figure 2B of ref. 8). Projections from our stochastic model (Fig. 3), add to this forecast by suggesting that with an increased case detection and treatment probability (1) of 70%, the prevalence of tuberculosis may decline and the disease may even be locally eradicated. However, the proportion of cases associated with resistance will increase (see Fig. 3 and SI Appendix).
We remark that these predictions are likely to be optimistic for two reasons. First, they require consistent and sustained success in detection and treatment over long periods of time. Second, if compensatory mutations occur, the transmission cost would continue to decline, leading to fitter resistant strains and serious epidemiological consequences (27). The design of effective treatment programmes will need to take account of both the magnitude of the fitness cost, as well as its future evolution via compensatory mutation.
Materials and Methods
We briefly present the model and methods here and provide the details in the SI Appendix; see also Table S1, Table S2, Table S3, Table S4, Table S5, and Figs. S1–S3 for additional information. To estimate epidemiological parameters from molecular epidemiological data, we (i) develop a stochastic model of the dynamics of tuberculosis transmission, the evolution of drug resistance, and mutation of genetic markers, and (ii) apply approximate Bayesian computation (ABC) to generate posterior distributions of parameters from empirical data. We analyzed data from population level studies of M. tuberculosis isolates from patients in Cuba (22), Estonia (10), and Venezuela (23). These datasets were identified through a comprehensive examination of a large pool of candidate datasets described in previous work (24). The key criterion for selection of a dataset was that the full distribution of genotypes in the sample along with the drugresistance status of every isolate was provided. These datasets used standard methods for molecular typing and drugresistance testing, they were collected over relatively short time periods, and were from regions where information about the time of introduction of antibiotics is available.
The stochastic model extends a previously described birth–deathmutation model (28) by including antibiotic resistance and sensitivity. A population of infectious individuals is divided into several subpopulations (see Fig. 1). This model tracks only infectious cases. Individuals are classified according to whether the strain they are infected with is resistant to at least one antibiotic or not and whether they are undergoing treatment or not. Treated individuals carrying drugresistant strains are further subdivided according to whether resistance was acquired through treatment failure or through transmission. Individuals move from one subpopulation to another through the acquisition of resistance or the commencement of treatment. They join the system through transmission and leave it through death or recovery. Underlying this model structure is a further classification of individuals according to bacterial genotype.
There are six types of event in this model. Transmission of drugsensitive strains occurs at a rate of α per case per year. The transmission rate of resistant strains is given a cost c, so that resistant strains transmit at α (1 − c) per case per year. Note that when the parameter c is positive, it describes a transmission cost of resistance, and when it is negative, there is a transmission benefit to resistance. We assume that the cost to transmission c is independent of markerbased genotypes. Untreated cases are detected and commence treatment at a rate of τ per case per year. The rate at which drugsensitive cases acquire resistance due to treatment failure is ρ per case per year. Treated cases recover due to treatment at a rate of ε_{S} per case per year for sensitive strains and ε_{R} for resistant strains, with ε_{S} > ε_{R}. Recovery in the absence of treatment along with death occur at a combined rate of δ per case per year. Mutation at the marker of interest occurs at the rate μ per case per year. The model assumes that each infectious individual carries a single genotype of M. tuberculosis and that the genetic marker is not under selection. The evolution of the molecular marker thus does not influence the transmission dynamics with respect to drug resistance but does allow genetic variation in M. tuberculosis to be tracked. We simulate this process for the purpose of comparison with molecular data. Of the eight parameters in the model, δ, τ, ε_{S}, and ε_{R} are fixed by using estimates obtained from the literature and the others, α, c, ρ, and μ, are estimated (as shown in Table S1 and discussed in the SI Appendix).
Each simulation of this process begins with a population of 4,000 drugsensitive infectious cases and allows the population to evolve for 40 y. This duration models the period between the introduction of antibiotics in these countries in the 1950s and the collection of the samples in the 1990s. In each simulation, the distribution of genotypes in the initial infectious population is randomly generated according to the infinitealleles model (see the SI Appendix for more details). This has the effect of integrating out uncertainty in the initial distribution. The infinitealleles model has been shown to adequately describe the distribution of M. tuberculosis genotypes into clusters (24).
The relative fitness of drugresistant strains based on this model can be expressed as as derived in the SI Appendix.
The procedure for Bayesian inference uses the ABCPRC algorithm proposed in ref. 21. The algorithm generates a weighted sample from an approximation to the true, but intractable, posterior distribution. The approximation is based on the ability of the model to reproduce informative summary quantities of the data by embedding the intractable target posterior within an augmented model from which sampling is viable. The ability to reproduce summary quantities is evaluated through weighted Euclidean distance between simulated and observed summary statistics. We used 11 summary statistics and weighted them according to their sensitivity to variation in model parameter values. These statistics include the number of distinct genotypes, the proportion of isolates that are drug resistant, and a range of other statistics based on the relationship between resistance phenotype and marker genotype. Table S3 describes these statistics and SI Appendix Section 2 provides further details of the estimation procedure.
The future of drug resistance in M. tuberculosis was explored by using the posterior distribution of parameters obtained through the Bayesian methods and simulating forwards in time. For these projections, we assumed that 70% of cases are detected and treated, which is the target set by the WHO (29), corresponding to τ = 1.21. We consider two scenarios: (i) moderate cure rates based on 50% treatment success for sensitive cases (see Table S1), and (ii) an optimistic cure proportion for sensitive cases of 85%, which is again the WHO target (29). These success proportions correspond to cure rates of (i) ε_{S} = 0.52 and ε_{R} = 0.202 (moderate) and (ii) ε_{S} = 2.95 and ε_{R} = 0.47 (optimistic).
Acknowledgments
We thank C. Bergstrom and L. Simonsen for discussions; S. Alizon, K. Koelle, and T. Stadler for comments on the manuscript; R. Z. Aandahl for technical computational assistance; and C. Sola, S. Samper, and L. Aristimuño for providing the full dataset from Venezuela. This work was supported by the Australian Research Council under the Discovery Project Scheme (DP0556732 and DP0664970).
Footnotes
 ^{2}To whom correspondence should be addressed. Email: m.tanaka{at}unsw.edu.au

Author contributions: S.A.S., A.R.F., and M.M.T. designed research; F.L. performed research; F.L., S.A.S., H.J., A.R.F., and M.M.T. contributed new reagents/analytic tools; F.L., S.A.S., A.R.F., and M.M.T. analyzed data; and F.L., S.A.S., A.R.F., and M.M.T. wrote the paper;.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission. C.B. is a guest editor invited by the Editorial Board.

This article contains supporting information online at www.pnas.org/cgi/content/full/0902437106/DCSupplemental.
References
 ↵
 World Health Organization
 ↵
 ↵
 ↵
 ↵
 Sander P,
 et al.
 ↵
 Gagneux S,
 et al.
 ↵
 ↵
 ↵
 Dye C,
 Williams BG,
 Espinal MA,
 Raviglione MC
 ↵
 Kruuner A,
 et al.
 ↵
 ↵
 Wilson RW,
 et al.
 ↵
 Post FA,
 et al.
 ↵
 Faustini A,
 Hall AJ,
 Perucci CA
 ↵
 Mathema B,
 Kurepina NE,
 Bifani PJ,
 Kreiswirth BN
 ↵
 ↵
 ↵
 Dye C,
 Williams BG
 ↵
 Dye C,
 Espinal MA
 ↵
 ↵
 Sisson SA,
 Fan Y,
 Tanaka MM
 ↵
 ↵
 ↵
 ↵
 Tanaka MM,
 Francis AR
 ↵
 Basu S,
 et al.
 ↵
 ↵
 Tanaka MM,
 Francis AR,
 Luciani F,
 Sisson SA
 ↵
 World Health Organisation
Citation Manager Formats
More Articles of This Classification
Biological Sciences
Population Biology
Related Content
 No related articles found.
Cited by...
 Reversing resistance: different routes and common themes across pathogens
 Selectable Markers for Use in Genetic Manipulation of Extensively DrugResistant (XDR) Acinetobacter baumannii HUMC1
 Tradeoffs between and within scales: environmental persistence and withinhost fitness of avian influenza viruses
 Exact vs. Approximate Computation: Reconciling Different Estimates of Mycobacterium tuberculosis Epidemiological Parameters
 Increasing Viral Dose Causes a Reversal in CD8+ T Cell Immunodominance during Primary Influenza Infection due to Differences in Antigen Presentation, T Cell Avidity, and Precursor Numbers
 Modeling the Dynamic Relationship Between HIV and the Risk of DrugResistant Tuberculosis
 Inferring Epidemiological Parameters on the Basis of Allele Frequencies
 The evolution of drug resistance and the curious orthodoxy of aggressive chemotherapy
 Mycobacterium bovis BCGRussia Clinical Isolate with Noncanonical Spoligotyping Profile