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
Modeling the impact of global tuberculosis control strategies
Related Article
 Tuberculosis: Old problems and new approaches Nov 10, 1998
Abstract
An epidemiological model of tuberculosis has been developed and applied to five regions of the world. Globally, 6.7 million new cases of tuberculosis and 2.4 million deaths from tuberculosis are estimated for 1998. Based on current trends in uptake of the World Health Organization’s strategy of directly observed treatment, shortcourse, we expect a total of 225 million new cases and 79 million deaths from tuberculosis between 1998 and 2030. Active casefinding by using mass miniature radiography could save 23 million lives over this period. A single contact treatment for tuberculosis could avert 24 million cases and 11 million deaths; combined with active screening, it could reduce mortality by nearly 40%. A new vaccine with 50% efficacy could lower incidence by 36 million cases and mortality by 9 million deaths. Support for major extensions to global tuberculosis control strategies will occur only if the size of the problem and the potential for action are recognized more widely.
Each year, there are more than 6.5 million new cases of tuberculosis and more than 2 million deaths from tuberculosis worldwide (1). Tuberculosis is the seventh most important cause of global premature mortality and disability (2) and is projected to remain among the 10 leading causes of disease burden even in the year 2020 (3). Because of a powerful interaction between tuberculosis and HIV, tuberculosis incidence is rising in subSaharan Africa and may rise in Asia. Increasing drugresistance, which could lead to worse treatment outcomes, also has been reported (4).
By using directly observed treatment, shortcourse (DOTS), cure rates of 80 to 90% have been achieved for passively diagnosed cases of smearpositive pulmonary tuberculosis (5). Analyses of national programs in Malawi, Tanzania, and Mozambique and in 10 provinces of China have shown that this strategy is both effective and costeffective (6–8). Based on the success of these programs, the World Health Organization (WHO) adopted DOTS as its strategy for global tuberculosis control, but uptake by national programs has been slow; WHO reports that only 11% of new smearpositive pulmonary tuberculosis cases are enrolled in DOTS programs worldwide (9). Given the pace of uptake of DOTS, the heavy emphasis of WHO on smearpositive pulmonary cases only, and the magnitude of the remaining problem, it is worth examining whether considerable gains against tuberculosis can be achieved through promotion of DOTS alone. In this paper, we evaluate a range of extensions to global control strategies in terms of their potential effects on tuberculosis incidence and mortality by region from 1998 to 2030.
Since the 1960s, simple mathematical models have been used to understand tuberculosis transmission dynamics and to predict the effects of different interventions (10). Recently, there has been renewed interest in using mathematical models to study tuberculosis epidemics (11–14). Blower et al. (11, 12) have used a simple analytical framework to calculate threshold levels of preventive therapy and treatment necessary for tuberculosis eradication. We have elaborated the Blower model to capture some other aspects of the epidemiology of tuberculosis relevant to evaluating control strategies. This model has been applied to real populations, which allows us to study the impact of specific control measures in different regions of the world.
METHODS
The basic structure of the model appears in Fig. 1. The states and parameters in the model are described in the Appendix, with complete technical details available elsewhere†. In the model, individuals enter the uninfected category (U) at birth. An individual may be infected at rate λ and move into one of two categories: fast breakdown (I_{F}), which includes those who progress quickly to clinical tuberculosis over the next 1–5 years; or slow breakdown (I_{S}), which includes those who progress slowly to clinical disease. Only a small fraction of the individuals in the slow breakdown category actually will develop the disease. By modeling the infected state as two categories distinguished by the rate of disease progression, it is possible to capture a composite breakdown rate that declines exponentially from the time of infection‡. Infected individuals who receive isoniazid (INH) preventive therapy and are subject to a reduced rate of breakdown move from I_{F} and I_{S} into the INH slow breakdown category (H_{S}).
Fig. 1 also shows an alternative path from the I_{S} and H_{S} categories to clinical disease, through the superinfected fast breakdown category (S_{F}). Considerable evidence suggests that superinfection occurs and is epidemiologically significant§. The rate of transfer into S_{F} is the risk of infection, λ, times the proportion developing fast breakdown, (1 − p), times (1 − v), where v is the protection from subsequent fast breakdown that a primary infection affords. Some interventions, such as DOTS, have a much larger epidemiological impact if superinfection is common.
Individuals who progress to clinical disease from the infected, superinfected, or preventive therapy categories enter one of the six states for prediagnosed and otherwise untreated cases, represented collectively as C_{U}^{i,j}. The index i identifies three clinical categories: smearpositive pulmonary, smearnegative pulmonary, and extrapulmonary tuberculosis. The index j distinguishes individuals subject to a fast rate of diagnosis from those subject to a slow diagnosis rate. Thus, the casedetection rate and diagnostic delay can be controlled separately for each type of clinical disease, which is important to capture the impact of different interventions. Untreated cases may spontaneously remit at rate ɛ_{U} or may die at a rate that combines the general mortality risk and the tuberculosis fatality rate (25–27).
Once diagnosed, cases may receive no treatment and remain in the same C_{U} category, or they may begin treatment and move into one of the six states labeled C_{T}^{i,k}. Here, i defines clinical type as above while k distinguishes cases receiving good treatment from those receiving bad treatment. The good and bad treatment categories have different cure rates and death rates (25). The bad treatment category includes a range of circumstances in which an individual fails to complete an adequate course of chemotherapy. Treated cases may be cured at a categoryspecific rate, ɛ_{T}^{k}. Recovered cases, either from spontaneous remission or from treatment cure, may relapse (28–30) or become superinfected. Two categories of recovered cases are included (R_{F} and R_{S}), distinguishing those who are subject to a fast relapse rate from those subject to a slow relapse rate.
To capture the profound effects of HIV on the development of tuberculosis, the model used here actually is composed of two submodels—one for the HIV seronegative population and one for the HIV seropositive population. Each submodel follows the structure depicted in Fig. 1. In the HIVpositive submodel, many of the transfer rates are different because of the effects of HIV on the immune system. Individuals move from each category in the HIVnegative submodel to the corresponding category in the HIVpositive submodel at the HIV infection rate, which varies over time. Because the effects of HIV on immune function are not marked with respect to tuberculosis until the CD4 count drops below 500, individuals actually move from the HIVnegative to the HIVpositive submodel after they have been infected with HIV for 3 years. The two submodels also are linked through the infection rate λ because HIVnegative tuberculosis cases can infect HIVpositive individuals, and vice versa. (See Appendix for details on heterogeneous mixing).
The models are solved by using the finite difference method with a time step of 36.5 days. Sensitivity analysis shows that the inaccuracy introduced by this time step is small. Initial conditions have been developed for each state in each region as described below.
Regional Models.
Based on patterns of tuberculosis epidemiology, the world has been divided into five regions, adapted from the World Bank regions used in the Global Burden of Disease Study (8): (i) EME, Established Market Economies; (ii) FSE, Formerly Socialist Economies of Europe; (iii) Asia, composed of India, China, and Other Asia and Islands; (iv) SSA, SubSaharan Africa; and (v) LAC/MEC, a combination of Latin America and the Caribbean with the Middle Eastern Crescent.
In order to elaborate regional models, we first defined plausible ranges for each of the model parameters, based on a review of the published and unpublished literature. Details are included in the Appendix. For each of the five regions, we began by modeling a virgin epidemic; that is, we simulated the introduction of a single infectious case into a population of susceptible individuals. Each epidemic was allowed to run to equilibrium. In LAC/MEC, Asia, and SSA, treatment was introduced into the equilibrium models in order to reproduce the declining epidemics of the last halfcentury. In EME and FSE, where tuberculosis rates have been declining for longer, we have reproduced these declines through a series of incremental changes reflecting reductions in transmission and breakdown, with the addition of treatment starting around the middle of this century.¶
An important determinant of the regional models starting in the 1980s is the HIV epidemic, for which we incorporated estimates developed by the former Global Programme on AIDS (3). In each model, the treatment variables have been adjusted during the 1980s and 1990s to reflect changes in case detection and treatment. The resulting models match closely the trends in the annual risk of infection and levels of incidence and mortality reported by WHO (31–33) and the Global Burden of Disease Study (2). For each region, we have selected only one model; clearly, other combinations of parameter values and shifts in these values over time might also provide adequate fits to known trends in tuberculosis. Uncertainty analysis for a similar model for the United States, however, indicates that evaluations of the impact of different interventions are largely unaffected by the choice of parameter values that fit past trends.‖
Global Control Strategies.
The likely impact of pursuing only the WHO DOTS strategy was evaluated, with a range of assumptions about achievable levels of uptake. Two sets of extensions to the DOTS strategy also were investigated: (i) aggressive applications of existing technologies and (ii) development and application of new technologies. For each scenario, we have estimated incidence and mortality in the HIVnegative and HIVpositive populations, from 1998 to 2030. In order to examine benefits in the longer term, we also have obtained results through the year 2050.
WHO DOTS Strategy.
Three variants of the WHO DOTS strategy were developed based on the extent to which DOTS is adopted, adequately funded, and given sufficient managerial attention. In DOTSM, recent trends in DOTS uptake are expected to continue, while DOTSH and DOTSL refer to a high and a low variant, respectively (Table 1). In all three scenarios, it is assumed that DOTS uptake will continue to increase over the next two decades and will level off by 2020.
In examining the potential impact of alternative strategies, two important sources of future uncertainty were identified: the course of the HIV epidemic and the possible effects of high coverage of Bacille CalmetteGuerín (BCG) vaccination. Three different Global Programme on AIDS projection scenarios for HIV incidence were used (3). We also considered the possibility of sustained BCG efficacy (beyond 15 years after vaccination). A recent metaanalysis suggested a 50% efficacy, but this conclusion is controversial, as the set of trails shows heterogeneity of efficacy (35). Even in populations where BCG appears to be efficacious, there are few data available on its effects after 15 years (36). Without sustained efficacy, BCG is unlikely to lower tuberculosis infection or incidence significantly because it is given primarily at birth, and children are rarely smearpositive (37). With dramatically increasing BCG coverage since the late 1970s, however, the first birth cohorts with high BCG coverage now are beginning to reach age 15 (38). If BCG has sustained efficacy, it might greatly reduce incidence over the next few decades.
In order to investigate the effects of these two sources of uncertainty, we have defined a baseline DOTS scenario (DOTSM with baseline HIV projections and no sustained BCG efficacy), a bestcase DOTS scenario, which assumes DOTSH in the most favorable context (optimistic HIV projections and sustained BCG efficacy of 50%), and a worstcase DOTS scenario, which assumes DOTSL in the least favorable context (pessimistic HIV projections and no sustained BCG efficacy). Each of the alternative strategies described below has been evaluated incrementally to the baseline, bestcase, and worstcase DOTS scenarios.
Existing Technologies.
Expanding implementation of DOTS for smearnegative cases (DOTSNeg). It is assumed that, gradually over the next 5 years, an expanded DOTS strategy could achieve the same cure rates in smearnegative patients as in smearpositive patients.
Active casefinding, symptomatic (ACFSx). The entire population would be screened for respiratory symptoms every 7 years (≈14% of the population each year). It is assumed that the capacity to implement this strategy would develop gradually over the next 10 years. The 3–4% reporting symptoms (39) would be subject to sputum examination. Using this algorithm, we expect to detect twothirds of prevalent smearpositive cases (40). Of detected cases, 90% are assumed to start treatment, with cure rates the same as in DOTSM.
Active case finding by using mass miniature radiography (ACFMMR). The entire population would be screened every 7 years by using mass miniature radiography (MMR), followed by sputum examination for those with suspicious lesions. It is assumed that capacity to implement this strategy will be built over the next 10 years. This screening approach would detect 95% of prevalent smearpositive cases (40), with 90% of detected cases starting treatment. Smearnegative xray suggestives meeting strict clinical criteria also would be treated empirically. We assume that onethird of true smearnegatives would be detected and that 80% of these cases would be enrolled in treatment.
Singlecycle active casefinding (ACFSx1 and ACFMMR1). For both the symptomatic and MMR screening algorithms, we also have examined the potential impact of just one complete cycle of screening. Because capacity will develop gradually, this one cycle would require 12 years to achieve full coverage of the population.
Mass preventive therapy for the entire adult population (MassPT). This intervention would include one cycle of active casefinding by using MMR (with assumptions as above), combined with a Mantoux test. Individuals with a positive skin test but without xray suggestive findings and clinical symptoms would be offered a 6month course of INH preventive therapy. It is assumed that twothirds of the individuals given the Mantoux test would return for the results, 90% of those found positive would start treatment, the efficacy of 6month INH is 67% (41), and onethird of those starting treatment would complete therapy. This single cycle of MassPT would require 12 years for full implementation.
Preventive therapy for HIVpositive patients (PTHIV5 and PTHIV10). Given the current and likely future access to voluntary HIV testing centers, we have studied two scenarios: in the first, 5% of HIVpositive individuals would be screened and offered PT each year whereas the second scenario would extend to 10% of the HIVpositive population (42–44), with 33% treatment completion assumed in both scenarios.
New Technologies.
Improved sputum examination (Smear75). It is assumed that, within five years, a method could be available that would raise the sensitivity of sputum examination from 50 to 75% (45, 46). In each region, the additional patients diagnosed by using this method over conventional smear are assumed to be smearnegatives in terms of transmission potential; extra cases diagnosed would have the same cure rate as smearpositives.
Nucleic acid amplification test (Smear95). Assumptions are as above, except that the test is assumed to be 95% sensitive (47).
Ultra short course chemotherapy (UltraSCC). It is assumed that, within 10 years, a single contact drug regimen could raise cure rates among treated cases to 95%.
Breakdown vaccines (VacBr20, VacBr50, and VacBr80). These vaccines would require 15 years for development and would prevent breakdown from infection to disease. Vaccine efficacies of 20, 50, and 80% have been studied, with assumed coverages of 80%. It is assumed that the vaccine would be administered to the population gradually over a 10year period and then subsequently to every child at birth.
Infection vaccines (VacInf20, VacInf50, and VacInf80). Three analogous scenarios have been elaborated for vaccines that would be given to uninfected persons to prevent infection. It is assumed, however, that infected individuals would be identified by using a Mantoux test and that only twothirds of the individuals tested would return to have their test results read.
RESULTS
Based on the regional models, Table 2 summarizes the regional and global epidemiology of tuberculosis in 1998. Globally, there are 6.7 million new cases and 2.4 million deaths expected in 1998.
Fig. 2 shows the progression of the global tuberculosis epidemic for the high, medium, and low variants of DOTS uptake, assuming baseline HIV projections and no sustained BCG efficacy. In all three scenarios, rising case numbers and deaths caused by interaction with the HIV epidemic will continue until the middle of the next decade. In DOTSM, we expect 225 million new cases of tuberculosis and 79 million deaths between 1998 and 2030. With higher DOTS uptake, cases might be reduced by 25.7 million and deaths by 11.4 million whereas lower uptake would lead to 16.4 million extra cases and 8.0 million extra deaths.
Depending on the course of the HIV epidemic and the possibility of sustained BCG efficacy, there is a range of uncertainty around the potential impact of the DOTS strategy, with the most uncertainty in Asia and SSA. Globally, we expect 171 million new cases and 60 million deaths in the bestcase DOTS scenario and 249 million new cases and 90 million deaths in the worstcase DOTS scenario.
Table 3 describes the benefits that might be derived from extensions to the DOTS strategy. Results are reported for each strategy evaluated incrementally to the baseline DOTS scenario. Expanding the coverage of DOTS to include diagnosed smearnegatives would reduce total tuberculosis mortality by 2% between 1998 and 2030. A single cycle of active casefinding with MMR could prevent 11% of the deaths expected in this period in the baseline scenario. Continuous MMR screening could reduce incidence by 22% and mortality by 29% (Fig. 3). The impact of active casefinding by using a symptomatic screen rather than MMR would be ∼40–50% smaller.
Preventive therapy for HIV patients diagnosed at voluntary testing centers could save nearly 1 million lives if 10% of the HIVpositive population were screened each year. Of the averted deaths, 36% would be in the HIVnegative population because of reduced transmission. The benefits from a single round of mass preventive therapy would be considerable; even with only 33% treatment completion, this strategy would avert 42.6 million cases and 16.7 million deaths (Fig. 3).
Table 3 also provides estimates of cases and deaths averted if new technologies were developed and were applied widely. Improved sputum smear methods or nucleic acid amplification tests could prevent 1.9 to 3.6 million cases and 2.6 to 4.7 million deaths during the period 1998 to 2030. A single contact treatment for active tuberculosis would prevent 23.8 million cases and 10.7 million deaths (Fig. 3). In combination with active casefinding, singledose chemotherapy could reduce tuberculosis mortality by 41% from 1998 to 2030 and by 51% over the period 1998 to 2050.
Even with a 15year time lag for development and clinical trials, vaccines that prevent breakdown or infection would make a major impact between now and 2030. With 50% efficacy, a vaccine that prevents breakdown from infection to disease could prevent 16% of cases and 12% of deaths (Fig. 3). All of the vaccines evaluated prevent a much larger percentage of cases and deaths when the analysis extends farther into the future. Over the period 1998 to 2050, a breakdown vaccine with 50% efficacy would reduce incidence by 29% and mortality by 25%.
We also have evaluated the impact of each strategy incremental to the bestcase and worstcase DOTS scenarios. This uncertainty analysis shows that, for many interventions, the expected percentage reduction in tuberculosis mortality is not affected by the range of DOTS uptake, the course of the HIV epidemic, or the efficacy of BCG (although clearly, the absolute numbers of cases and deaths averted are strongly affected). Notable exceptions include new vaccines, which have much smaller benefits if BCG provides sustained protection and a much greater impact in the worstcase scenario. Cases and deaths averted through the introduction of an ultrashortcourse regimen also depend strongly on DOTS uptake.
DISCUSSION
Even with high uptake of DOTS in the most favorable context, incidence and mortality from tuberculosis in the next three decades will be enormous. It is therefore imperative to consider means of extending the DOTS strategy. Many of the strategies evaluated in this paper would require substantial resources and an extraordinary mobilization of international attention. Initiating and sustaining such a largescale effort would be challenging, but clearly the first step must be a realistic assessment of the potential benefits of these ambitious initiatives.
In this study, we have applied epidemic models to real populations to draw conclusions about one of the world’s leading health challenges. We have not attempted to cost different global control strategies. To translate the results in this paper into policy recommendations, it will be important to proceed with a realistic evaluation of the costs of pursuing active casefinding, mass preventive therapy, and the development of new technologies.
WHO, for many years, has recommended passive casefinding as the primary strategy for detecting tuberculosis (48). The results of this analysis, however, indicate that active casefinding could reduce tuberculosis mortality by onequarter to onethird. Perhaps the time has come for a critical reexamination of the costs, risks, and benefits of active screening.
The concept of INH preventive therapy has been recognized for decades, but in only one case has mass preventive therapy been used to maximum effect (49). In lowincidence communities, the low risk of fatal isoniazid toxicity—0.014% in the U.S. Public Health Service trial (50)—has been a source of great concern. In high incidence regions, however, this small toxicity risk would be far outweighed by the magnitude of the benefits. Even with modest completion rates, a mass preventive therapy campaign could alter significantly the course of the tuberculosis epidemic.
Although aggressive use of existing technologies could produce considerable benefits, tuberculosis likely will kill more than 50 million people over the next three decades, even if these aggressive strategies are pursued. Investing now in new tools could make a substantial difference (51). A single contact treatment for tuberculosis would have an enormous effect and, of all the new technologies, would be the easiest to field. Recent breakthroughs in the treatment of paucibacillary leprosy (52) make consideration of a single contact treatment for tuberculosis worthwhile. By reducing drug and labor costs, such a development probably would prove to be costsaving. Over the next halfcentury, the combination of a single contact treatment and active casefinding could reduce mortality from tuberculosis by >50%. Investments in new vaccines also could have a profound impact over this period.
Tuberculosis program managers and control specialists must develop more ambitious objectives for tuberculosis control over the next halfcentury. Tuberculosis is likely to remain one of the major causes of death in the 21^{st} century, but it is also one of the few for which there are a number of existing and potential new interventions that can make a dramatic difference. A major rethinking of global tuberculosis strategy can occur only if the size of the problem and the potential for action are recognized more widely.
Appendix
The model is determined by the following system of ordinary differential equations: where T is the number of new births; λ is the infection rate; μ is the nontuberculosis death rate; p is the proportion of new infections entering slow breakdown; β_{F} and β_{S} are the fast and slow breakdown rates, respectively; χ is the rate at which infected individuals receive INH preventive therapy; v is the protection from superinfection conferred by primary infection; w is shortterm INH effectiveness; h is longterm INH effectiveness; s^{i} is the proportion of new cases in clinical category i; d^{i,j} is the proportion of prediagnosed cases of clinical type i entering diagnosis category j; δ^{j} is the diagnosis rate for category j; σ is the rate of conversion from smearnegative to smearpositive; ɛ_{U} is the spontaneous cure rate for untreated cases; μ_{U}^{i} is the tuberculosis death rate for untreated cases in clinical category i; g^{i,k} is the proportion of treated cases in clinical category i and treatment category k; ɛ_{T}^{k} is the cure rate for treated cases in treatment category k; μ_{T}^{i,k} is the tuberculosis death rate for treated cases in category i,k; r_{U} is the proportion of spontaneously recovered cases entering the slow relapse category; r_{T}^{k} is the proportion of recovered cases from treatment category k entering the slow relapse category; and ρ_{F} and ρ_{S} are the fast and slow relapse rates, respectively. In the equations for C_{U}^{i,j} and C_{T}^{i,k}, the smear rate σ is multiplied by the number of individuals in the respective category i∗, where i∗ = 2 (smearnegative) for i = 1 (smearpositive) and vice versa, and i∗ = ⊘ for i = 3 (extrapulmonary). The term including σ is added for i = 1, subtracted for i = 2, and equal to zero for i = 3. The result of this formulation is that smearnegative patients convert to smearpositive at a rate of σ.
In this model, λ is a function of the prevalence of infectious sources in the population. That is, λ = KL (Y/N), where λ is the instantaneous rate of infection for a susceptible; K is the number of respiratory contacts per person per year; L is the probability that a respiratory contact with an infectious source will lead to tuberculosis infection; Y is the number of infectious cases in the population (in terms of the model states, Y includes all smearpositive cases plus a small proportion of smearnegatives); and N is the population. Thus, it is assumed that the number of secondary infections caused by a source case in a year (KL) does not depend on population size. This formulation of λ often is used in models of sexually transmitted disease (53).
Because HIVpositive individuals may be more likely to come in contact with other HIVpositive individuals, we have elaborated the expression for λ to allow for heterogeneous mixing: where i refers to the population subgroup (HIV seropositive ≡ 1; HIV seronegative ≡ 2), j is the complementary subgroup, K, L, Y, and N are as described above, ω_{ii} is the proportion of all effective respiratory contacts that are with members of the same subgroup, and ω_{ij} is the proportion of all effective contacts that are with members of the complementary subgroup. We note three identities implied by these definitions: (i) ω_{11} = 1 − ω_{12}; (ii) (KL)_{2}N_{2}ω_{21} = (KL)_{1}N_{1}ω_{12}; and (iii) ω_{22} = 1 − ω_{21}. We have formulated ω_{11} to be proportionate to HIV seroprevalence; that is, ω_{11} = c∗ N_{1}/(N_{1} + N_{2}) for some constant c, which may vary by region. A large value for c implies that HIV prevalence is much higher among the contacts of an HIVpositive individual than among the general population. Thus, proportionate mixing is implied where c = 1 whereas contacts occur exclusively within subgroups when c is the inverse of HIV seroprevalence; that is, c = (N_{1} + N_{2})/N_{1} ⇒ ω_{11} = 1.
Following are the ranges considered for the major parameters of interest in the model: transmission parameter (KL), 5–15 infections per infectious case per year; proportion of new infections breaking down to clinical disease within the first 5 years (1 − p), 5–15% in HIVnegatives and 50–70% in HIVpositives; endogenous reactivation rate (β_{S}), 50–150 persons per 100,000 persons per year in HIVnegatives and 5–10% per year in HIVpositives; proportion of new cases of tuberculosis that are smearpositive pulmonary (s^{1}): 40–50%; protection against breakdown from subsequent infection afforded by primary infection (v), 50–100%; death rate for untreated tuberculosis (μ_{U}), 0.075–0.20 per person per year in HIVnegatives and 0.3–0.6 per person in HIVpositives; cure rate for untreated tuberculosis (ɛ_{U}), 0.05–0.15 per person per year; proportion of recovered cases in fast relapse category (1 − r), 4–15%; slow relapse rate (ρ_{S}), 50–150 persons per 100,000 persons per year (11, 12, 16–18, 25–30). Complete details on the parameter values used in each regional model over time are available in ref. 15.
Footnotes

↵* To whom reprint requests should be addressed. email: jsalomon{at}hsph.harvard.edu.

A Commentary on this article begins on page 13352.

↵† The structure and dynamics of the model, as well as technical details on the modeling of alternative control strategies, are described fully in ref. 15. This report is available at http://www.hsph.harvard.edu/organizations/bdu.

↵‡ Sutherland (16) reported on the time to breakdown in individuals who developed tuberculosis in the Medical Research Council vaccine trial. Nearly 60% developed tuberculosis within 1 year of tuberculin skin test conversion, and the proportion declined exponentially as a function of time since conversion. Krishna Murthy et al. (17) found similar results.

↵§ The significance of superinfection has been a highly controversial issue for many years (18). Evidence in favor of superinfection includes anatomical evidence (19), restriction fragment length polymorphism analyses (20, 21), and epidemiological studies (22–24).

↵¶ See ref. 15 for complete details on parameter values in the regional models.

↵‖ The model for the United States is part of an ongoing project sponsored by the Centers for Disease Control and Prevention. The methodology for undertaking multivariate uncertainty analysis, while at the same time validating a model by its ability to fit past trends in incidence and mortality, is described in ref. 34. In brief, plausible ranges are defined for each parameter in the model, multiple sets of parameter combinations are generated by randomly sampling each range, and the model is run by using each of these parameter combinations. The subset of combinations that optimize a particular goodnessoffit criterion comprise the set of valid models, which then are used to evaluate the future impact of alternative policy scenarios.
ABBREVIATIONS
 WHO,
 World Health Organization;
 DOTS,
 directly observed treatment, shortcourse;
 BCG,
 Bacille CalmetteGuérin;
 MMR,
 mass miniature radiography;
 INH,
 isoniazid
 Copyright © 1998, The National Academy of Sciences
References
 ↵
 World Health Organization Global Tuberculosis Programme
 ↵
 Murray C J L,
 Lopez A D
 ↵
 ↵
 The WHO/IUATLD Global Project on AntiTuberculosis Drug Resistance Surveillance
 ↵
 Porter J D H,
 McAdam K P W J
 Broekmans J F
 ↵
 ↵
 World Bank
 ↵
 ↵
 Waaler H,
 Geser A,
 Andersen S
 ↵
 ↵

 Brewer T F,
 Heymann J,
 Colditz G,
 Wilson M E,
 Auerbach K,
 Kane D,
 Fineberg H V
 ↵
 Vynnycky E,
 Fine P E
 ↵
 Murray C J L,
 Salomon J A
 ↵
 Sutherland I
 ↵
 Krishna Murthy V V,
 Nair S S,
 Gothi G D,
 Chakraborty A K
 ↵
 Styblo K
 ↵
 Canetti G,
 Sutherland I,
 Svandova E
 ↵
 Turett G S,
 Fazal B A,
 Justman J E,
 Alland D,
 Duncalf R M,
 Telzak E E
 ↵
 ↵

 Krishna Murthy V V,
 Chaudhuri K
 ↵
 ↵

 Drolet G J
 ↵
 Berg G
 ↵
 ↵
 Chakraborty A K,
 Gothi G D
 ↵
 Raviglione M C,
 Snider D E,
 Kochi A

 World Health Organization. Global Tuberculosis Programme
 ↵
 Cauthen G M,
 Pio A,
 ten Dam H G
 ↵
 Murray C,
 Sawert H
 ↵
 Colditz G A,
 Brewer T F,
 Berkey C S,
 Wilson M E,
 Burdick E,
 Fineberg H V
 ↵
 ↵
 ↵
 United Nations Children’s Fund
 ↵
 ↵
 ↵
 ↵
 ↵
 ↵
 ↵
 Saceanu C A,
 Pfeiffer N C,
 McLean T
 ↵
 Forbes B A
 ↵
 World Health Organization
 ↵
 Chakraborty A K,
 Channabasavaiah R,
 Krishnamurthy M S,
 Shashidhara A N,
 Motiram G
 ↵
 ↵
 Ad Hoc Committee on Health Research Relating to Future Intervention Options
 ↵
 ↵