Tuberculosis (TB) is a global health emergency that demands concerted management efforts. Recent World Health Organization reports reveal that MDRTB is a substantial problem in every region evaluated4
. However, the impact of the emergence of multidrug-resistant (MDR) strains on global TB control remains unclear. Throughout much of the world, control strategies include case detection through sputum microscopy and standardized treatment regimens administered without drug sensitivity testing. Although this approach is effective in managing TB that is drug-sensitive or resistant to one antibiotic, it neither detects nor cures the majority of cases of MDRTB.
This exclusive focus on drug-sensitive and monoresistant tuberculosis arises from two considerations. The first is the cost and complexity of detecting and treating drug-resistant tuberculosis, and the second is the perception that drug-resistant pathogens are less likely to cause epidemics than drug-sensitive ones. This perception is based on the observation that mutations that confer resistance often alter gene products important for the pathogen's survival and thus exact a ‘fitness cost’5
Mathematical models have been developed to describe the transmission dynamics and the impact of interventions for many infectious diseases6,7
, including drug-sensitive8-11
and drug-resistant tuberculosis1-3,12
. Previous work on the population dynamics of drug resistance has shown that the fitness of drug-resistant strains is a key determinant of the future burden of drug-resistant pathogens in general13,14
and specifically of drug-resistant TB1
. When the average fitness of MDR strains is assumed to be substantially less than drug-sensitive strains, these models predict that the emergence of multidrug resistance will not jeopardize the success of control efforts that focus on the treatment of drug-sensitive disease.
The reproductive fitness of an individual strain of M. tuberculosis
is a complex characteristic determined by the bacterium's ability to infect a susceptible host, persist and proliferate, and be transmitted to a secondary host. Laboratory studies comparing the relative fitness of drug-resistant and drug-sensitive organisms through competition assays have found that although some mutations compromise growth and the ability to survive oxidative stress15
, others have minimal in vitro
effects. Furthermore, when fitness costs were detected, they were often short-lived because compensatory mutations rapidly accumulated to restore biological function16
. The results of empirical studies that have compared the reproductive potential of drug-resistant and drug-sensitive strains in human populations have been similarly heterogeneous; some reports indicate that drug-resistant strains of M. tuberculosis
spread less readily than drug-sensitive strains17
, whereas others show no difference in disease transmission18
Despite the heterogeneity in the fitness of MDR strains suggested by the studies cited, previous models of the transmission dynamics of drug-resistant TB have assumed a single fixed value of average fitness for MDR strains. Given that drug-resistant strains of M. tuberculosis
show a range of fitness19
, we constructed a model that incorporates the emergence of multidrug resistance during therapy, heterogeneity of fitness of MDR strains, and competition during an epidemic (). Following previous TB models, we assumed that susceptible individuals are infected at a rate that is proportional to the prevalence of people with active disease; once infected, individuals progress to disease through a fast route (primary disease) or a slow route (reactivation), and individuals who are infected but do not have active disease may become reinfected (exogenous reinfection). To incorporate competition between MDR and drug-sensitive organisms, we modeled three distinct strains: one that is drug-sensitive and two MDR strains that differ in relative fitness compared to the drug-sensitive strain. The ‘unfit’ MDR strain (designated by subscript ‘U’) has a low fitness (0.3) relative to the drug-sensitive strain, whereas the ‘fit’ MDR strain (subscript ‘F’) has a relative fitness that ranges from 0.8 and 1.2 (). As evidence suggests that most mutations conferring a MDR phenotype will incur an initial fitness cost, we assumed that mutation leads to the unfit MDR phenotype (aU
) 100 times as often as the fit MDR phenotype (aF
Figure 1 Multistrain tuberculosis model structure. Boxes represent state variables; the flow between states is described by the differential equations in Methods. States of infection with MDR strains are identified by subscript ‘U’ (‘unfit’ (more ...)
Definition and values of model parameters
We initially modeled the course of a drug-sensitive tuberculosis epidemic prior to the introduction of combination chemotherapy. This simulation yielded measures of disease consistent with those from selected high-burden regions (). We then modeled a 30-year period before the introduction of Directly Observed Treatment, Short-course (DOTS) treatment, during which drugs were available, but detection and cure rates reflected the suboptimal program performance reported before the implementation of the DOTS strategy. We estimated that approximately 3% of those with initially drug-sensitive disease who fail pre-DOTS treatment acquired multidrug resistance (). The results of this simulation reflect the current epidemiology observed in high-burden countries (). These results also show that the early course of the MDRTB epidemic is insensitive to the relative fitness of the MDR phenotype, with negligible differences in the estimated burden of drug-sensitive or MDRTB over a wide range of fitness costs.
Model fit and simulation results
To explore the long-term dynamics of MDRTB, we then modeled the implementation of the DOTS program beginning 30 years after the initiation of unstructured treatment. We assumed case detection and cure rates for drug-sensitive disease consistent with World Health Organization goals21
, improved treatment efficacy for MDR disease and a reduced proportion of those acquiring the multidrug-resistant phenotype through failed treatment (). As was observed in many countries that saw significant reductions in TB incidence after antibiotics became available, drug-sensitive disease is rapidly brought under control through the implementation of DOTS. In contrast to the short-term results discussed above, the simulated long-term course of tuberculosis epidemics depends critically on the relative fitness of MDR strains ().
Figure 2 Sensitivity of long-term projections of MDRTB epidemics to the relative fitness of MDR strains. Colored traces represent infection and disease with the MDR strains; black traces represent infection and disease with the drug-sensitive strain. Colors represents (more ...)
In this simulation, almost all of the MDR cases generated early in the course of the epidemic harbor the unfit strain; thus, the average fitness is low. When the prevalence of drug-sensitive tuberculosis is reduced through the implementation of DOTS, the impact of acquired drug-resistance on the emerging MDR epidemic becomes less important and the predominant driving force of incident MDR infections shifts to the transmission of MDR strains. Consequently, the fitness difference between the circulating strains becomes the primary determinant of new infections and interstrain competition leads to the ascent of the most transmissible strains ().
The model demonstrates that even when the most fit MDR strain is assumed to be less fit than the drug-sensitive strain, the MDR strain will eventually outcompete the drug-susceptible strain (). However, there exists a fitness threshold for the ‘fit’ MDR strain below which multidrug-resistant disease will not continue to spread as long as DOTS treatment goals are met. In these simulations, this threshold is exceeded when the most fit MDR strains are greater than 70% as fit as the drug-sensitive strains.
This model makes several simplifying assumptions. First, because the relationship between contact patterns and the spread of airborne infections is not well understood, we modeled the population as a homogenously mixing unit. Although this is oversimplified, we based this decision on molecular epidemiological data that suggest substantial transmission of TB occurs following casual contact with infectious TB cases26
and contact tracing fails to identify most secondary cases27
. Second, we modeled the acquisition of multidrug resistance as a single-step process despite the fact that MDR is acquired through sequential accumulation of point mutations. Because monoresistant strains are thus included in the drug-sensitive (non-MDR) population, we may have underestimated the relative fitness of MDR strains in comparison to the non-MDR strains in the model. Finally, we modeled tuberculosis dynamics in a limited set of scenarios consistent with epidemiological data from developing countries. Recent work demonstrates that tuberculosis dynamics are highly sensitive to the rate of infection and suggests that it may not be possible to generalize our results to the setting of a developed country in which the incidence and prevalence of tuberculosis are declining28
Our model suggests that despite a short-term decline in the tuberculosis burden following the implementation of DOTS, the exclusive treatment of drug-sensitive and monoresistant tuberculosis may contribute to the emergence of MDR disease in at least two ways: (i) treatment directly contributes to the pool of MDR cases through acquired resistance, albeit at much lower rates under DOTS than pre-DOTS programs, and (ii) the removal of infectious individuals with predominantly drug-sensitive disease decreases the force of infection and subsequently replenishes the pool of individuals who are fully susceptible to infection with circulating MDR strains.
Previous studies suggest that the fitness of MDR strains is heterogeneous. We demonstrate that future burden of MDRTB is more dependent on the distribution of fitness among circulating strains than on the initial average reproductive fitness. We show that even if highly transmissible MDR strains are only rarely generated through poor treatment, MDRTB may eventually become a major public health threat. Our simulations indicate that DOTS policies should be coupled with strategies to limit the spread of MDRTB (for example, DOTS-plus29
) in order to mitigate the long-term threat of MDR disease. Although an optimal approach would include the concurrent administration of high-quality DOTS programs and the introduction of second-line drug therapy for those with MDRTB, limited resources have made this difficult to achieve. How this combined approach can best be implemented in resource-constrained settings is a topic for urgent operational research.