- Research article
- Open Access
- Open Peer Review
Tuberculosis control strategies to reach the 2035 global targets in China: the role of changing demographics and reactivation disease
BMC Medicinevolume 13, Article number: 88 (2015)
In the last 20 years, China ramped up a DOTS (directly observed treatment, short-course)-based tuberculosis (TB) control program with 80% population coverage, achieving the 2015 Millennium Development Goal of a 50% reduction in TB prevalence and mortality. Recently, the World Health Organization developed the End TB Strategy, with an overall goal of a 90% reduction in TB incidence and a 95% reduction in TB deaths from 2015–2035. As the TB burden shifts to older individuals and China’s overall population ages, it is unclear if maintaining the current DOTS strategy will be sufficient for China to reach the global targets.
We developed an individual-based computational model of TB transmission, implementing realistic age demographics and fitting to country-level data of age-dependent prevalence over time. We explored the trajectory of TB burden if the DOTS strategy is maintained or if new interventions are introduced using currently available and soon-to-be-available tools. These interventions include increasing population coverage of DOTS, reducing time to treatment, increasing treatment success, and active case finding among elders > 65 years old. We also considered preventative therapy in latently infected elders, a strategy limited by resource constraints and the risk of adverse events.
Maintenance of the DOTS strategy reduces TB incidence and mortality by 42% (95% credible interval, 27-59%) and 41% (5-64%), respectively, between 2015 and 2035. A combination of all feasible interventions nears the 2035 mortality target, reducing TB incidence and mortality by 59% (50-76%) and 83% (73-94%). Addition of preventative therapy for elders would enable China to nearly reach both the incidence and mortality targets, reducing incidence and mortality by 84% (78-93%) and 92% (86-98%).
The current decline in incidence is driven by two factors: maintaining a low level of new infections in young individuals and the aging out of older latently infected individuals who contribute incidence due to reactivation disease. While further reducing the level of new infections has a modest effect on burden, interventions that limit reactivation have a greater impact on TB burden. Tools that make preventative therapy more feasible on a large scale and in elders will help China achieve the global targets.
Significant progress in tuberculosis (TB) control has been achieved worldwide over the last two decades. Global TB mortality has fallen by 45%, and TB incidence is declining . Recently, the World Health Organization (WHO) established an ambitious post-2015 global strategy, the End TB Strategy . This strategy outlines a 2025 milestone of 50% reduction in incidence and 75% reduction in mortality, and an overall 2035 target of 90% reduction in incidence and 95% reduction in mortality. In order to reach these targets, countries will likely need to redouble their TB control efforts and perhaps adopt new TB control strategies .
Between 1992 and 2012, China made impressive progress in TB control. Prior to 1992, most TB patients were treated in private hospitals, where patients typically received low-quality care - improper treatment was widespread, and only approximately 20% of patients had supervised TB treatment. In addition, nearly 50% experienced interrupted or shortened treatment and there was little follow-up of patients who dropped out or relapsed after a treatment episode [4-7]. Starting in 1992, China ramped up a high-quality directly observed treatment, short-course (DOTS)-based strategy in Center for Disease Control (CDC) public health clinics in 13 provinces covering half the population, requiring hospitals to refer suspected TB patients to the CDC system. In the early 2000s, the DOTS program was expanded nationwide and an Internet-based disease reporting system was introduced [8-11], further increasing referrals from the hospital to the CDC system. By 2010 it was estimated that approximately 80% of all TB patients were confirmed and treated within the CDC system [8,9], where the treatment success rate was estimated to be 85% .
Simultaneously with the DOTS ramp-up, using serial nationwide prevalence surveys, China documented a 65% reduction in prevalence of smear-positive TB over the period from 1990 to 2010 [6-9]. These gains enabled China to achieve the 2015 global TB control target of halving TB prevalence five years before the target date . In spite of these gains, China still has nearly a million incident TB cases annually and nearly a quarter of the world’s multi-drug-resistant (MDR) TB cases .
Given the success of the DOTS ramp-up in the CDC public health clinics, it is possible that simply maintaining this level of TB control will enable China to reach the 2035 global targets. However, as the risk of infection has fallen, the burden of disease has shifted to older population groups, which is reinforced by a growing elderly population in China [12-14]. Thus, reactivation disease poses a growing challenge to TB control, as has been observed in Hong Kong . This could limit the continued impact of the DOTS strategy, which does not specifically address reactivation disease in the short term.
To assess the potential impact of reactivation disease and the aging demographics in China, we developed an individual-based transmission model of TB in the Chinese population, explicitly considering the changing demographics of China and the changes in the patient treatment seeking pathway over the last 20 years. This individual-based model builds on the modeling efforts of previous groups [14,16-23], enabling flexibility in the use of available demographics data and changing pathways to care over time.
Using realistic age structure and fitting to age-dependent prevalence data, we independently estimate the contribution of reactivation disease as a driver of TB incidence. In this context, we explore the impact of a basis set of interventions, based on available or soon-to-be-available tools. This basis set represents a set of independent interventions which, in combination, represent all intervention strategies that could possibly be implemented using the available or soon-to-be-available tools. All interventions were modeled optimistically with 100% nationwide coverage starting in 2015. These interventions were parameterized based on what was deemed feasible to implement in the Chinese health care system; they include:
increased access to care, by increasing patient referrals from private hospitals to the CDC such that all TB patients are confirmed and treated in the CDC.
reduced time to treatment by reducing provider and diagnostic delay using new diagnostics and/or streamlining the diagnostic pathway. This has been achieved in smaller pilot studies with varying levels of success and can reduce the time to treatment up to 50% [24-29].
increased treatment success within the CDC, using new drugs/drug regimens which would be effective in patients with both drug-sensitive (DS) and MDR TB. Several new drug regimens are currently in clinical trials or have recently been approved, including bedaquiline, the REMox regimen, and the PaMZ regimen, which have the potential to increase the treatment success rate regardless of MDR status [30-35].
active case finding in elders > 65 years old, by combining active TB screening with the annual health screening done in this group . Individuals found by active case finding would be provided with care in the CDC health system.
preventative therapy in elders > 65 years old, where latent screening of patients is done in combination with the annual health screening done in this age group. Although in reality latent screening would be done in combination with active case finding, we first model preventative therapy alone in order to explore the basis set of interventions. The combination with active case finding is also modeled. We also note that while preventative therapy is included in the list of interventions because it is a currently available tool, it is not deemed feasible in China due to the age of the screened population, and the relative risks of hepatic adverse events in this age group [37,38]. In addition, general population screening is not considered feasible at this time due to the overall scale of the population.
We use the model to estimate the contribution of reactivation disease to overall burden at the current time. We also quantify the impact of maintaining the DOTS strategy or expanding the TB control strategy to include new interventions. From this analysis, we discuss the feasibility of China reaching the 2035 global targets using existing tools for TB control.
The present study utilizes the Disease Transmission Kernel (DTK) model developed by the Institute for Disease Modeling group at Intellectual Ventures. The model and all necessary input files are available by request at the Institute for Disease Modeling website . Additional file 1 details the model structure, assumptions, and a complete list of model inputs. The model schematic is included in Figure 1. The fit of the mean trajectory to the data during the calibration time period is shown in Figure 2.
We use a simulated population of 500,000 individuals in 1990 that grows to 720,000 by 2035. This represents a 0.05% sampling of the true population size in China. The absolute population value was chosen to reflect available computational resources. The change in population is driven by the UN Population database estimates for medium fertility and age-disaggregated non-TB mortality . The overall population growth is shown in Additional file 1: Figure S1.
We use a 100-year burn-in period before the calibration period, which is from 1990–2012. The purpose of the burn-in period is to enable simulated individuals to have age-appropriate TB exposure by the end of the burn-in period and the beginning of the model calibration period. At the beginning of the simulation, the population is initialized with 5% prevalence of latent TB with no treatment available until 40 years before the start of the calibration period, at which time hospital treatment is made available to 90% of patients (the remaining 10% of patients have no access to care). The burn-in period is continued for another 40 years until the start of the calibration period (1990), at which time the prevalence and incidence of TB are near steady state, reflecting the observed relative stability in TB prevalence measured in 1990 and 2000 in provinces which did not have any TB-specific interventions during the 1990s (these provinces received DOTS-based treatment during the 2000s) [4-11]. Thus, the absolute value of the population TB prevalence at the beginning of the burn-in period did not affect the steady state equilibrium which defined the start of the calibration period. We allow the acquisition of MDR beginning five years before the end of the burn-in period, reflecting the introduction of rifampin in 1985 and the known non-steady state growth of MDR in 1990 . In order to model the two-step expansion of DOTS [4-11] and because the rate of interprovincial migration is less than 5% , we model two distinct homogeneous mixing pools (representing the provinces that received DOTS in the 1990s and 2000s, respectively) with the population evenly divided between them. Acknowledging the reality and challenge of parameterizing the age-stratified population mixing and age-stratified migration rates at the country level, we assume transmission is age-independent in this analysis.
TB natural history
We stratify disease progression based on age group (adults > 15 and children < 15 years old) and hold constant all disease parameters throughout the simulation, to clearly delineate the role of the aging population structure in China. We did not further stratify the disease parameters in elders . The summary list of relevant parameters is included in Table 1. Individuals are born uninfected and all individuals are vaccinated with bacille Calmette-Guerin (BCG), which reduces the probability of infection by 50% for 10 years [44-47]. We also include an age-dependent reduced susceptibility to infection of 40% from age 2–10, reflecting the BCG-independent reduced susceptibility to TB during the “safe school years” in this age group [48,49]. Once infected with latent TB, individuals are binned into either a latent fast or latent slow stage, representing the distribution of rates with which individuals may progress from latent to active disease. For individuals with the slow rate of progression, this probability of progression is calibrated. Individuals in the latent stage are not infectious and no additional mortality is associated with this phase . Individuals binned into the latent fast stage progress to active symptomatic with a median time of 4 months [51,52]. We calibrate the rate of reactivation in latently infected individuals binned into the latent slow group. Reflecting the spectrum of disease from latent to active disease [53-59], individuals progress from latent disease to an active presymptomatic phase before progressing to a symptomatic phase. The presymptomatic phase represents a period of reduced infectivity where patients may have objective signs of TB (such as increased cough) but lack subjective symptoms (that is, they do not notice their increased cough). The duration of this stage was estimated to be 1 year [56-59]. No additional mortality is associated with this stage.
Active symptomatic disease is split into three clinically defined states: smear-positive pulmonary, smear-negative pulmonary, and extrapulmonary. Smear-positive disease is more likely in adults (65% of all active symptomatic cases) but represents a smaller fraction of active disease in children [35,48,60,61]. Only pulmonary forms are infectious, with smear-negative disease being only 15% as infectious as smear-positive disease and extrapulmonary disease being non-infectious [62-64]. We calibrated the contact rate, the number of individuals who would be infected by an untreated individual with smear-positive disease. The duration of symptomatic disease has a mean of 5.5 years if no treatment is available  and is not modeled as variable with clinical presentation. At the end of the symptomatic disease duration, a smear-positive infection will either self-resolve or result in TB-related mortality . The relative ratio of natural recovery to TB-related mortality varies with clinical presentation and does not vary with age. In smear-negative patients and extrapulmonary patients, the disease can either self-resolve or result in TB-related mortality .
Acquired and transmitted MDR TB
MDR TB and DS TB are independently tracked in the model. We do not track acquisition of resistance to individual drugs or further resistance on top of MDR (that is, extensively drug-resistant (XDR) TB), as these were not expected to have a significant effect on our analysis given the relatively small contribution of MDR to overall incidence. MDR TB can be acquired during treatment for DS TB, occurring at the rate specified in the section on TB treatment: treatment outcomes. We assume that the MDR strain is 85% as fit as the DS strain [66-70].
TB treatment: pathways to care
Two treatment pathways are modeled: the private hospital system and the CDC system with its public health TB clinics [8,11], as shown in Figure 1B. Parameterization of the time to treatment and treatment outcomes was based on expert opinion from the Chinese CDC.
Once an individual enters the active symptomatic phase, they start seeking treatment with a median duration of 4 months [24-28,71]. We do not model provider visits which result in no treatment or the time for patient referral from hospitals to the public health system, so the time to treatment encompasses the overall time of treatment seeking until the start of treatment, including patient delay, diagnostic delay, and provider delay. Within the CDC, the sensitivity of diagnosis, regardless of smear status, is thought to be over 95%, using a combination of symptom screening, X-ray, and ultimately culture . Because of this high sensitivity, we do not disaggregate the many levels of the diagnostic pathway, as this was not expected to have a significant impact on our analysis of the role of reactivation. Thus, at the time of treatment initiation, all patients who are TB positive are provided with treatment.
We model only the provision of first- and second-line treatment, and do not explicitly model each individual drug or the adherence patterns to each drug (see TB treatment: treatment outcomes). Most patients receive first-line treatment. Reflecting the current level of MDR testing and provision of second-line treatment, only 1.3% of MDR-positive, smear-positive patients within the CDC treatment pathway receive second-line treatment. Although this proportion is likely slightly higher in retreatment patients and there may be additional time delay associated with MDR testing, we do not explicitly model this, as it is not expected to have a large effect on our analysis of the baseline or new interventions. The treatment outcomes of MDR patients receiving first- or second-line treatment are listed in Additional file 1: Table S4.
At the end of the treatment duration, those who are cured return to the susceptible pool, with a 65% reduced susceptibility to reinfection [73-75]. Individuals who are initially thought to be cured but will ultimately relapse are tracked separately from true cure. These individuals progress through a latent non-infectious phase having a median duration of 9 months , before progressing directly to the active symptomatic phase. At the start of the active symptomatic phase, they seek treatment again. Those who fail a treatment reseek care with a rate dependent on which health system they received their most recent treatment - within the CDC the median time to retreatment is 3 months, reflecting better follow-up in the CDC. If they remain in the hospital system, the median time to retreatment is set at 22 months. Individuals who failed treatment in the hospital sector have a 20% probability of shifting to the CDC for retreatment. These values were set based on expert opinion from the Chinese CDC and to achieve a parsimonious fit with the known data on the proportion of new and retreatment patients within the CDC and hospital systems.
TB treatment: treatment outcomes
Treatment within the CDC system is generally of higher quality than that available in private hospitals. We model only DOTS and second-line combination therapy, and do not disaggregate individual drugs. We also do not model individual treatment adherence patterns, instead including those who drop out as treatment failures. All treatment outcomes were based on data available from the Chinese National TB Control Program and expert opinion [4-7].
In the hospital system, treatment outcomes were set as follows: 55% have long-term cure, 26% fail during treatment, 11% were initially cured but then relapse, and 8% die during treatment. Among those who fail or relapse, there is a 10% probability of developing MDR. In the CDC system, treatment outcomes were set as follows: 82% have long-term cure, 9% fail during treatment, 9% were initially cured but then relapse, and 1% die during treatment. Those who fail or relapse have a 2% probability of developing MDR. These treatment outcomes were slightly poorer if patients were treatment experienced, and are detailed in Additional file 1: Table S3. The parameterization for new treatment based on expected treatment outcomes using new drugs is also described in Additional file 1: Table S3.
MDR patients who received DOTS (that is, due to lack of MDR diagnosis) had treatment outcomes set as follows: 35% have long-term cure, 40% fail during treatment, 10% were initially cured but then relapse, and 20% die during treatment. For the small portion of MDR patients who received second-line treatment in the CDC, the treatment outcomes were slightly higher: 60% have long-term cure, 15% fail during treatment, 10% were initially cured but then relapse, and 15% die during treatment.
TB treatment: DOTS ramp -up and the shifting access to care
During the calibration period from 1990–2012, we model the historical ramp-up of DOTS according to historically observed patterns [4-11]. During the 1990s, DOTS was implemented in the CDC system through public health clinics in 13 provinces covering half of China’s population, and subsequently expanded nationwide during the 2000s. Treatment within the CDC system was generally of higher quality than that available in private hospitals. (see TB treatment: treatment outcomes). Each of the DOTS ramp-ups was modeled as a linear expansion occurring over three years. From 1992–1995, in the provinces where DOTS expansion occurred in the 1990s, the proportion of patients who did not receive care was reduced from 10% to 5%, and among patients who received care, 60% of them were shifted from the hospital to the CDC. In 2002–2005, changes in the treatment pathways were expanded to the entire country. Country-wide, the proportion of patients with no access to care was reduced to 5%. Of those who did receive care, 80% of patients were shifted from the hospital to the CDC. These transitions are described in Additional file 1: Table S1 and S2.
The simulation is calibrated to the TB burden (age-dependent prevalence, smear-positive prevalence, and overall prevalence) in China from 1990–2010 as estimated by the Ministry of Health prevalence surveys done in 1990, 2000, and 2010. We also calibrate to the percentage of MDR in new and retreatment patients (survey done in 2007) and the estimated percentage of MDR in all patients (estimated by the Ministry of Health prevalence surveys) [1,7,8,41,77,78].
The model parameters that were calibrated were the contact rate (the infectiousness of a person with smear-positive TB, the average probability of TB transmission per timestep from an infected individual), the fraction of latently infected adults who were classified as fast progressors, and the rate of progression from latent to active in individuals classified as slow progressors. None of these parameters have been directly measured on the country level in China and were deemed most likely to exert the dominant effects on the population-level TB burden and the size of the latent reservoir.
All other parameters were informed by available literature and expert opinion, and were held fixed during the calibration. This reduces the total number of dimensions in the calibration parameter space to a computationally tractable size. This set of fixed parameters includes the total duration of active disease and the relative rate of seeking care in naïve and retreatment patients, which was informed by available information on the time to initial treatment, time to retreatment, and total number of provider visits. While these durations are relevant for analyzing specific interventions which shorten the pathway to care, these values are not expected to have a large effect on our analysis of the relative importance of new transmission and reactivation disease. As our data is restricted to the country-level data, we do not explore individual heterogeneity in the model parameters in this study.
Bayesian calibration using incremental mixture importance sampling [79,80] was employed for parameter space exploration and calibration. Initially, a uniform prior is used for each of the parameters across a plausible parameter range informed by the literature and expert opinion, as detailed in Table 1. We use Latin hypercube sampling to initially explore 1,000 parameter combinations, after which we use the likelihood weighted posterior to select the next round of 60 parameter combinations, for a total of 100 rounds, resulting in 7,000 total parameter combinations sampled from the three-dimensional parameter space.
The likelihood was calculated as the product of the likelihoods for each data point (see Additional file 1 for details on the likelihood calculation and calibration methods). Subsequently, to evaluate the baseline and the impact of new interventions, the likelihood weighted parameter space is resampled 100 times, resulting in a total of 26 unique parameter combinations. These parameter combinations were rerun using 10 random number seeds and averaged together to reduce the stochastic noise. The weighted mean of these parameter combinations thus includes both parameter and stochastic uncertainty. The 26 unique parameter combinations were used to estimate the Bayesian 95% credible interval (see Additional file 1 for details on the calculation).
All scenarios are modeled optimistically with 100% country-wide implementation in 2015 and impact measured by 2035. This represents the upper limit of impact that could be achieved if these intervention strategies were implemented as described.
The first intervention strategy is to increase access to care by shifting the patients currently treated in public hospitals (20% of all patients) to the CDC system, where the quality of treatment is higher. This means 100% of TB patients in China would be treated in a CDC DOTS program.
The second intervention strategy involves reducing the time to treatment by reducing provider and diagnostic delay using new diagnostics and/or streamlining the diagnostic pathway. We reduced the time to treatment by 33% for all patients, regardless of which system they used to receive treatment. Treatment naïve patients’ time to treatment was reduced from a median of 128 days to 88 days, for all patients seen in the CDC and hospital system. The time to retreatment was not changed [24-29].
Another strategy is increased treatment success within the CDC, using new drugs/drug regimens which would be effective in both DS and MDR patients [30-35]. When given to DS patients, the estimated treatment outcomes are 92% long-term cure, 3.5% fail during treatment, 3.5% are initially cured but then relapse, and 1% die during treatment. This was slightly lowered for MDR and treatment experienced patients. See Additional file 1: Table S3 for the full breakdown of treatment outcome.
Active case finding in elders > 65 years old, by combining TB screening with the annual health screening done in this population, is another intervention strategy. Although in reality this would be done throughout the year, this was modeled as a single yearly occurrence for all individuals over age 65, where anyone who was in the active symptomatic phase immediately received treatment from the CDC .
The final strategy is preventative therapy in elders > 65 years old, where screening of patients is done in combination with active case finding. This was also modeled as a single yearly occurrence for all individuals over age 65. The diagnostic test to identify latently infected individuals was not explicitly modeled. Latent treatment was parameterized as a 9-month regimen with a cure rate of 80% . Individuals who were not cured returned to the latent phase and were eligible to be retreated in subsequent iterations.
Results and discussion
Model structure and calibration
We developed a dynamic microsimulation transmission model which tracks individuals from birth to death, including acquisition of latent infection, progression to active disease, treatment seeking behavior, and treatment status. A model schematic is included in Figure 1 illustrating the progression of disease and treatment pathways. Key model inputs are shown in Table 1, and the fit of the mean trajectory to the data during the calibration time period is shown in Figure 2. The decline in incidence over the calibration period is the result of both aging out of the latent reservoir and reduced incidence stemming from infectious individuals seeking retreatment, as high-quality DOTS has lowered the proportion of treated individuals who fail treatment. Additional model outputs are available in Additional file 1.
The model estimates that if the status quo in TB diagnosis and treatment is maintained, TB incidence and mortality will decline slowly by 42% (27-59%) and 41% (5-64%, 95% credible interval), respectively, between 2015 and 2035 (Figure 2, Table 2). The combination of continued aging out of the latent reservoir and a low annual risk of infection is consistent with further aging of the epidemic (Additional file 1: Figure S4). MDR is projected to remain relatively stable at below 10% of overall incidence (Additional file 1: Figure S4).
The baseline projection represents a mean of several points resampled from the calibration parameter space. Notably, there are significant tradeoffs among the calibration parameters that can result in similarly good fits to the data but which estimate different future trends in incidence (Figure 3). Simulations having a higher contact rate and a lower fast progressor fraction (orange curve, based on calibration points within the orange box in Figure 3A) have a larger latent reservoir, and incidence is dominated by slow progressors reactivating from this reservoir. In contrast, simulation which have a lower contact rate and a higher fast progressor fraction (purple curves, based on calibration points within the purple box in Figure 3A) have a smaller latent reservoir; in this parameter space, the total incidence is lower and is projected to decline to a lower level if the DOTS strategy is maintained. The model estimate for the future trend in mortality is consistent with the model estimate in the trend of overall incidence. As variation in the calibration parameters does not affect the individual case fatality rate, simulations that estimate a higher overall incidence also estimate a higher overall mortality.
Future intervention strategies
All of the modeled interventions are parameterized based on feasibility within the existing health care ecosystem of private hospitals and public CDC clinics. The relative impact of these interventions is described in Table 2 and shown in Figure 4.
Increased access to care could be achieved by increasing patient referrals from private hospitals to the CDC. This would increase from 80% to 100% the percentage of TB patients who are confirmed and treated in the CDC clinics, where treatment success rates are higher (see Additional file 1 for detailed tables of the treatment success rates). It would also provide all patients with better follow-up if retreatment is necessary, reducing the infectiousness from treatment experienced individuals. The model estimates that this will result in an estimated reduction in TB incidence and mortality by 47% (31-63%) and 65% (54-79%), respectively, over 20 years. This intervention is the most effective single intervention that is feasible.
Reducing the time to treatment shortens the duration of infectiousness in treatment naïve individuals, and could be achieved by using new diagnostics and/or streamlining the diagnostic pathway [24-29]. A one-third reduction in mean time to treatment from 128 days to 88 days for treatment naïve patients in both the hospital and CDC system would result in a limited impact in both TB incidence and mortality compared to baseline.
Improving treatment success via more effective drug regimens and better treatment monitoring is the second best intervention. This would be effective for both DS and MDR patients, and reduces the need for retreatment in all patients. By doing so, this intervention also reduces the infectiousness stemming from treatment experienced individuals. The model estimates that this will result in an estimated reduction in TB incidence and mortality by 49% (35-64%) and 60 (43-77%), respectively, over 20 years.
Active case finding in elders reduces the time to treatment, primarily for treatment naïve individuals but also affecting those who are treatment experienced. In addition, because it is combined with the annual health screening already done for elders > 65 years old, this intervention would be targeted rather than used in the general population. It is estimated that this would result in a decline in TB incidence and mortality of 48% (34-64%) and 58% (40-72%) over 20 years.
Preventative therapy in elders > 65 years old would be the most effective single intervention if it could be made feasible in this age group by addressing the relative risks of hepatic adverse events [37,38]. Preventative therapy is modeled as a 9-month treatment with an overall treatment cure rate of 80% . This strategy directly reduces the size of the latent reservoir. In the first few years, the bulk of the latent reservoir is treated, resulting in a rapid decline in incidence from 2015–2025, and reaching the 2025 milestone for incidence decline. However, the subsequent decline in incidence from 2025 to 2035 is not steep enough to reach the 2035 global target.
The combination of all feasible interventions (all interventions except for preventative therapy) has a larger impact than any of the individual feasible interventions alone but still a smaller impact than preventative therapy. Between 2015–2035, implementing all feasible interventions is estimated to result in a 59% (50-76%) and 83% (73-94%) decline in incidence and mortality, nearly achieving both the 2025 milestone for mortality and the 2035 mortality target. With regards to incidence, the steepest decline is seen between 2015–2025, nearly reaching the 2025 milestone, but the total decline from 2015–2035 does not reach the 2035 incidence target. Including preventative therapy with all feasible interventions has the greatest impact of all the modeled interventions, and is likely to enable China to nearly reach the 2035 incidence and mortality global targets. Between 2015–2035, the projected change in incidence and mortality is −84% (78-93%) and −92% (86-98%).
Intervention impact is sensitive to calibration parameters
Recognizing that the specific calibration parameters present strong differences in the estimated trend of incidence and mortality (Figure 3), we directly compared the impact of the feasible interventions from different areas in the calibrated parameter space (Figure 5). Comparing model outputs which use the orange and purple boxes in parameter space, the projection of incidence from 2015–2035 is quite divergent, both at baseline and with all feasible interventions. The absolute impact of implementing all feasible interventions is larger if the higher contact rate is used. This is consistent with a higher absolute amount of recent transmission that is driven by the higher contact rate. The estimate for absolute incidence with all feasible interventions given a higher baseline contact rate (green line from orange baseline) is higher than the estimate from baseline at a lower contact rate (purple line). This suggests that in addition to modeling specific interventions, it is necessary to more completely specify the calibrated parameter space.
The combination of an aging demographic in China and the increasing role of reactivation disease represents a growing challenge to TB control as China considers its post-2015 strategy. We have constructed a mathematical model of TB transmission at the country level in China, taking into account aging of the population and estimating the contribution of reactivation to overall incidence. The nationwide roll-out of the DOTS program reduced the annual risk of infection(ARI) [81,82] by improving treatment outcomes and reducing infectiousness from treatment experienced individuals. Given the high population coverage of DOTS in the CDC public health clinics, we estimate that new transmission is not the major driver of overall TB incidence. Rather, reactivation disease, combined with the growing elderly population, will be the major determinant of the decline in TB incidence and mortality over the next two decades.
Our work shows that if the status quo DOTS strategy is maintained, the TB burden in China will decline but will not reach the 2025 milestones, even if an additional decade is provided. However, additional data are necessary to better specify what the baseline incidence trajectory might be, as projections from different points in the calibration parameter space are divergent over the next 20 years. This could include data on the percentage of the population that is latently infected, stratified by age. The current model estimate for the latent fraction is in line with that estimated in rural areas  and in smaller (non-country level) studies, which have primarily focused on city subpopulations or high-risk groups such as health care workers . It is possible that further data collection and analysis in this area can be used to make improved projections of incidence and mortality. Data on the proportion of incidence due to reactivation versus recent transmission could also help to specify the incidence trajectory, though country-specific data in this area is limited and considerable additional data collection would likely be necessary. As shown in this work, this improved specification of the model is equally important to the model estimate of new interventions.
The best single intervention is a system innovation whereby all TB patients would receive their initial treatment in a DOTS program. A combination of all feasible interventions, including expanded access to high-quality care, improved treatment quality, shorter treatment delay, and active case finding in elders will reduce incidence and mortality by 59% (50-76%) and 83% (73-94%), respectively, by 2035. This nearly achieves the 2035 mortality target and suggests that the 2025 incidence milestone can be achieved by 2035 if all feasible interventions are implemented. Further, if preventative therapy can be made feasible in elders, this would be a transformational intervention which is very likely to enable China to reach the 2035 targets.
Overall, our analysis of intervention strategies, selected based on availability of current tools and the structure of the Chinese health care system, suggests that if all feasible interventions are implemented, China may come close to achieving the 2035 mortality target but is unlikely to achieve the 2035 incidence target. Tools aimed at reducing reactivation from the latent reservoir will be critical to quickly reduce incidence in China. This could include a better drug regimen for preventative therapy and/or better monitoring [85,86]. Alternatively, non-TB-specific interventions that might affect the secular trend in TB disease, which were not explicitly modeled in this study, including improved nutrition and better living standards, hold forth the possibility of limiting the rate of reactivation [87-91].
Our model is limited by our assumptions regarding model structure and implementation of new interventions which affect our estimates of the TB burden. First, we assume the disease parameters, including the rate of disease progression and the infectiousness of active disease, remain constant during the entire simulation period. We do not explicitly model how secular trends in transmission related to changing living patterns might have affected the disease parameters [87-91]. Although this would directly influence our results, empiric data in this area are sparse and could either raise or lower our estimate of the role of reactivation disease. A growing population of healthy elders would have a lower likelihood of reactivation, while a growing population of elders living with immune-modulating diseases such as diabetes could result in an overall higher likelihood of reactivation with age . Secular trends in transmission due to changing living conditions, urbanization, and migration could also change the TB transmission rate, directly affecting our estimates of how new interventions might lower TB prevalence. We restrict our analysis in this study to the role of the aging population assuming the disease and transmission parameters remain constant.
Second, we have assumed homogeneous mixing, a simplifying assumption which does not account for spatial differences, non-uniform age-based mixing, and inter-individual heterogeneities in susceptibility and infectivity. Data from the TB prevalence surveys have shown a difference in TB prevalence between rural and urban areas , likely due to differences in socioeconomic status, contact patterns, demographic structure, and access to care. Rural-to-urban migrants may face additional delays in care seeking due to geographic and financial restrictions [92,93]. Further, age-dependent mixing is a notable contributor to the age-dependent incidence observed for a variety of respiratory illnesses [94-96], and also varies in urban and rural areas. Our current analysis is restricted to the country level, and does not explicitly model the population heterogeneity which could help or hinder the attainment of the TB control targets. In addition, we have not modeled heterogeneity among individuals and during the course of the infection.
We are optimistic that ongoing data collection in well-instrumented sites will enable better quantification of these unknowns, improving the quality and utility of TB modeling to inform TB control programs. Third, because we have not explicitly specified how these interventions would be operationalized, we do not address the relative cost of the interventions.
Our work suggests that after ramping up DOTS to a high population coverage, reactivation from the latent reservoir plays a growing role in driving incidence, and interventions using existing tools to further reduce new transmission will have a limited impact. To eliminate TB as a public health problem in the Chinese setting, transformative approaches that can limit or prevent reactivation of latent TB infection will likely be needed.
annual risk of infection
bacilli Calmette-Guerin vaccine
directly observed treatment short-course
disease transmission kernel
World Health Organization
extensively drug resistant
WHO: Global tuberculosis report 2014. http://www.who.int/tb/publications/global_report/en/.
WHO: Global strategy and targets for tuberculosis prevention, care and control after 2015. http://www.who.int/tb/post2015_tbstrategy.pdf?ua=1.
Lonnroth K, Castro KG, Chakaya JM, Chauhan LS, Floyd K, Glaziou P, et al. Tuberculosis control and elimination 2010–50: cure, care, and social development. Lancet. 2010;375:1814–29.
He GX, van der Hof S, van der Werf MJ, Guo H, Hu YL, Fan JH, et al. Inappropriate tuberculosis treatment regimens in Chinese tuberculosis hospitals. Clin Infect Dis. 2011;52:e153–6.
China Tuberculosis Control Collaboration. Results of directly observed short-course chemotherapy in 112 842 Chinese patients with smear-positive tuberculosis. Lancet. 1996;347:358–62.
Chen X, Zhao F, Duanmu H, Wan L, Wang L, Du X, et al. The DOTS strategy in China: results and lessons after 10 years. Bull WHO. 2002;80:430–6.
China Tuberculosis Control Collaboration. The effect of tuberculosis control in China. Lancet. 2004;364:417–22.
Wang L, Liu J, Chin DP. Progress in tuberculosis control and the evolving public-health system in China. Lancet. 2007;369:691–66.
Wang L, Zhang H, Ruan Y, Chin DP, Xia Y, Cheng S, et al. Tuberculosis prevalence in China, 1990–2010; a longitudinal analysis of national survey data. Lancet. 2014;2057–2064:383(9934).
Wan L, Cheng S, Chin DP. A new disease reporting system increases TB case detection in China. Bull WHO. 2007;85:401.
Wang L, Cheng S, Xu M, Huang F, Xu W, Li R, et al. Model collaboration between hospitals and public health system to improve tuberculosis control in China. IJTULD. 2009;13:486–92.
Mori T, Leung CC. Tuberculosis in the global aging population. Infect Dis Clin North Am. 2010;24:751–68.
Chatterji S, Kowal P, Mathers C, Naidoo N, Verdes E, Smith JP, et al. The health of aging populations in China and India. Health Aff. 2008;27:1052–63.
Dye C, Williams BG. The population dynamics and control of tuberculosis. Science. 2010;328:856–61.
Wu P, Lau EHY, Cowling BJ, Leung CC, Tam CM, Leung GM. The transmission dynamics of tuberculosis in a recently developed Chinese city. PLoS One. 2010;5:e10468.
Abu-Raddad LJ, Sabatelli L, Achterberg JT, Sugimoto JD, Longini IM, Dye C, et al. Epidemiological benefits of more-effective tuberculosis vaccines, drugs and diagnostics. Proc Natl Acad Sci U S A. 2009;106:13980–5.
Blower SM, Mclean AR, Porco TC, Small PM, Hopewell PC, Sanchez MA, et al. The intrinsic transmission dynamics of tuberculosis epidemics. Nature Med. 1995;1:815–21.
Dowdy DW, Lotia I, Azman AS, Creswell J, Sahu S, Khan AJ. Population-level impact of active tuberculosis case finding in an Asian megacity. PLoS One. 2013;8:e77517.
Brooks-Pollock E, Cohen T, Murray M. The impact of realistic age structure in simple models of tuberculosis transmission. PLoS One. 2010;5:e8479.
Mehra M, Cossrow N, Kambili C, Underwood R, Makkar R, Potluri R. Assessment of tuberculosis burden in China using a dynamic disease simulation model. IJTLD. 2013;17:1186–94.
Houben RM, Sumner T, Grant AD, White RG. Ability of preventative therapy to cure latent Mycobacterium tuberculosis infection in high burden settings. PNAS. 2014;111:5325–30.
Dye C, Garnett GP, Sleeman K, Williams BG. Prospects for worldwide tuberculosis control under the WHO DOTS strategy. Lancet. 1998;352:1886–91.
Borgdorff MW, Floyd K, Broekmans JF. Interventions to reduce tuberculosis mortality and transmission in low and middle-income countries. Bull WHO. 2002;80:217–27.
Sreeramareddy CT, Panduru KV, Menten J, Van den Ende J. Time delays in diagnosis of pulmonary tuberculosis: a systematic review of literature. BMC Infect Dis. 2009;9:91–101.
Li Y, Ehiri J, Tang S, Li D, Bian Y, Lin H, et al. Factors associated with patient, and diagnostic delays in Chinese TB patients: a systematic review and meta-analysis. BMC Med. 2013;11:156–70.
Storla DG, Yumer S, Bjune GA. A systematic review of delay in the diagnosis and treatment of tuberculosis. BMC Public Health. 2008;8:15.
Wei X, Zou G, Yin J, Walley J, Sun Q. Comparing patient care seeking pathways in three models of hospital and TB programme collaboration in China. BMC Inf Dis. 2013;13:93.
Xu B, Jian QW, Xiu Y, Diwan VK. Diagnostic delays in access to tuberculosis care in counties with or without the National Tuberculosis Control Programme in rural China. IJTLD. 2005;9:784–90.
Boehme CC, Nicol MP, Nabeta P, Michael JS, Gotuzzo E, Tahirli R, et al. Feasibility, diagnostic accuracy and effectiveness of decentralized use of the Xpert MTB/RIF test for diagnosis of tuberculosis and multidrug resistance: a multicenter implementation study. Lancet. 2011;377:1495–505.
Lienhardt C, Vernon A, Raviglione MC. New drugs and new regimens for the treatment of tuberculosis: review of the drug development pipeline and implications for national programmes. Curr Opin Pulm Med. 2010;16:186–93.
Diacon AH, Dawson R, von Groote-Bidlingmaier F, Symons G, Venter A, Donald PR, et al. 14-day bactericidal activity of PA-824, bedaquiline, pyrazinamide, and moxifloxacin combinations: a randomized trial. Lancet. 2012;380:986–93.
Gillespie SH, Crook AM, McHugh TD, Mendel CM, Meredith SK, Murray SR, et al. Four-month moxifloxacin-based regimens for drug-sensitive tuberculosis. NEJM. 2014;371:1577–87.
Diacon AH, Pym A, Grobusch MP, de los Rios JM, Gotuzzo E, Vasilyeva I, et al. Multidrug-resistant tuberculosis and culture conversion with bedaquiline. NEJM. 2014;371:723–32.
Gler MT, Skripconoka V, Sanchez-Garavito E, Xiao H, Cabrera-Rivero JL, Vargas-Vasquez DE, et al. Delamanid for multidrug-resistant pulmonary tuberculosis. NEJM. 2012;366:2151–60.
Jindani A, Harrison TS, Nunn AJ, Phillips PP, Churchyard GJ, Charalambous S, et al. High-dose rifapentine with moxifloxacin for pulmonary tuberculosis. NEJM. 2014;371:1599–608.
Sun X, Chen Y, Tong X, Feng Z, Wei L, Zhou D, et al. The use of annual physical examinations among the elderly in rural China: a cross-sectional study. BMC Health Serv Res. 2014;14:16.
Smith BM, Schwartzman K, Bartlett G, Menzies D. Adverse events associated with treatment of latent tuberculosis in the general population. CMAJ. 2011;183:E173–9.
Smieja M, Marchetti C, Cook D, Smaill FM. Isoniazid for preventing tuberculosis in non-HIV infected persons. Cochrane Database Syst Rev. 2000;2:CD001363.
Institute for Disease Modeling: Software. http://idmod.org/software.
United Nations. Department of Economic and Social Affairs, Population Division. http://www.un.org/popin/data.html.
Zhao Y, Xu S, Wang L, Chin DP, Wang S, Jiang G, et al. National survey of drug-resistant tuberculosis in China. NEJM. 2012;366:2161–70.
Fan C. Interprovincial migration, population redistribution, and regional development in China: 1990 and 2000 census comparisons. Prof Geogr. 2005;57:295–311.
Perez-Guzman C, Vargas MH, Torres-Cruz A, Villarreal-Velarde H. Does aging modify pulmonary tuberculosis? Chest. 1999;116:856–7.
Colditz GA, Brewer TF, Berkey CS, Wilson ME, Burdick E, Fineberg HV, et al. Efficacy of BCG vaccine in the prevention of tuberculosis: meta-analysis of the published literature. JAMA. 1994;271:698–702.
Fine PE. Variation in protection by BCG: implications of and for heterologous immunity. Lancet. 1995;346:1339–45.
Sterne JA, Rodrigues LC, Guedes IN. Does the efficacy of BCG decline with time since vaccination? IJTLD. 1998;2:200–7.
Trunz BB, Fine PEM, Dye C. Effect of BCG vaccination on childhood tuberculous meningitis and military tuberculosis worldwide: a meta-analysis and assessment of cost-effectiveness. Lancet. 2006;367:1173–90.
Cruz AT, Starke JR. Clinical manifestations of TB in children. Pediatr Respir Rev. 2007;8:107–17.
Marais BJ, Gie RP, Schaaf HS, Hesseling AC, Obihara CC, Starke JJ, et al. The natural history of childhood intra-thoracic tuberculosis: a critical review of literature from the pre-chemotherapy era. IJTLD. 2004;8:392–402.
Vynnycky E, Fine PM. Lifetime risks, incubation period, and serial interval of tuberculosis. Am J Epidemiol. 2000;152:247–63.
Comstock GW. Epidemiology of tuberculosis. Am Rev Respir Dis. 1982;125:8–15.
Achkar JM, Jenny-Avital ER. Incipient and subclinical tuberculosis: defining early disease states in the context of host immune response. J Infect Dis. 2011;204:S1179–86.
Barry 3rd CE, Boshoff HI, Dartois V, Dick T, Ehrt S, Flynn J, et al. The spectrum of latent tuberculosis: rethinking the biology and intervention strategies. Nat Rev Microbiol. 2010;7:845–55.
Lin PL, Flynn JL. Understanding latent tuberculosis: a moving target. J Immunol. 2010;185:15–22.
Azman AS, Golub JE, Dowdy DW. How much is tuberculosis screening worth? Estimating the value of active case finding for tuberculosis in South Africa, China, and India. BMC Med. 2014;12:216.
Dowdy DW, Basu S, Andrews JR. Is passive diagnosis enough? The impact of subclinical disease on diagnostic strategies for tuberculosis. Am J Respir Crit Care Med. 2013;187:543–51.
Menzies NA, Cohen T, Lin HH, Murray M, Salomon JA. Population health impact and cost-effectiveness of tuberculosis diagnosis with Xpert MTB/RIF: a dynamic simulation and economic evaluation. PLoS Med. 2012;9:e1001347.
Zwerling A, White RG, Vassall A, Cohen T, Dowdy DW, Houben RMGJ. Modeling of novel diagnostic strategies for active tuberculosis - a systematic review: current practices and recommendations. PLoS One. 2014;9:e110558.
Swaminathan S, Rekha B. Pediatric Tuberculosis: global overview and challenges. CID. 2010;50:S184–94.
Shah NS, Cavanaugh JS, Pratt R, Cain KP, Wells C, Laserson K, et al. Epidemiology of smear-negative pulmonary tuberculosis in the United States, 1993–2008. IJTLD. 2012;16:1234–40.
Tostmann A, Kik SV, Kalisvaart NA, Sebek MM, Verver S, Boeree MJ, et al. Tuberculosis transmission by patients with smear-negative pulmonary tuberculosis in a large cohort in the Netherlands. Clin Infect Dis. 2008;47:1135–42.
Behr MA, Warren SA, Salamon H, Hopewell PC, Ponce de Leon A, Daley CL, et al. Transmission of Mycobacterium tuberculosis from patients smear-negative for acid-fast bacilli. Lancet. 1999;353:444–9.
Hernandez-Garduno E, Cook V, Kunimoto D, Elwood RK, Black WA, Fitzgerald JM. Transmission of tuberculosis from smear negative patients: a molecular epidemiology study. Thorax. 2004;59:286–90.
Tiemersma EW, van der Werf MJ, Borgdorff MW, Williams BG, Nagelkerke NJ. Natural history of tuberculosis: duration and fatality of untreated pulmonary tuberculosis in HIV negative patients: a systematic review. PLoS One. 2011;6:e17601.
Dye C, Williams BG, Espinal MA, Raviglione MC. Erasing the world’s slow stain: strategies to beat multidrug-resistant tuberculosis. Science. 2002;295:2042–6.
Cohen T, Murray M. Modeling epidemics of multidrug-resistant M. tuberculosis of heterogeneous fitness. Nat Med. 2004;10:1117–21.
Espinal MA. The global situation of MDR-TB. Tuberculosis. 2003;83:44–51.
Borrell S, Gagneux S. Infectiousness, reproductive fitness and evolution of drug-resistant Mycobacterium tuberculosis. Int J Tuberc Lung Dis. 2009;13:1456–66.
Cohen T, Sommers B, Murray M. The effect of drug resistance on the fitness of Mycobacterium tuberculosis. Lancet Inf Dis. 2003;3:13–21.
Public Health England: Tuberculosis (TB) in the UK: annual report. https://www.gov.uk/government/publications/tuberculosis-tb-in-the-uk.
World Health Organization. Toman’s tuberculosis: case detection, treatment and monitoring. Geneva: World Health Organization; 2004.
Vynnycky E, Fine EM. The natural history of tuberculosis: the implications of age-dependent risks of disease and the role of reinfection. Epidemiol Infect. 1997;119:183–201.
Sutherland I, Svandova E, Radhakrishna S. Development of clinical tuberculosis following infection with tubercle bacilli: 1. A theoretical model for the development of clinical tuberculosis following infection, linking from data on the risk of tuberculous infection and the incidence of clinical tuberculosis in the Netherlands. Tubercle. 1982;63:255–68.
Andrews JR, Noubary F, Walensky RP, Cerda R, Losina E, Horsburgh CR. Risk of progression to active tuberculosis following reinfection with Mycobacterium tuberculosis. Clin Inf Dis. 2012;54:784–91.
Cao JP, Zhang LY, Zhu Q, Chin DP. Two year follow up of directly observed intermittent regimens for smear-positive pulmonary tuberculosis in China. IJTLD. 1998;2:360–4.
Ministry of Public Health of the People’s Republic of China. Nationwide random survey for the epidemiology of tuberculosis in 1990. Beijing: Ministry of Public Health of the People’s Republic of China; 1990.
Ministry of Public Health of the People’s Republic of China. Report on nationwide random survey for the epidemiology of tuberculosis in 2000. Beijing: Ministry of Public Health of the People’s Republic of China; 2000.
Raftery AE, Bao L. Estimating and projecting trends in HIV/AIDS generalized epidemics using incremental mixture importance sampling. Biometrics. 2010;66:1162–73.
Steele RJ, Raftery AE, Emond MJ. Computing normalizing constants for finite mixture models via incremental mixture importance sampling (IMIS). J Comput Graph Stat. 2006;15:712–34.
Zhang LX, Tu DH, He GX, Ma ZQ, Nagelkerke NJ, Borgdorff MW, et al. Risk of tuberculosis infection and tuberculous meningitis after discontinuation of BCG in Beijing. Am J Respir Crit Care Med. 2000;162:1314–7.
van Leth F, van der Werf MJ, Borgdorff MW. Prevalence of tuberculous infection and incidence of tuberculosis; a reassessment of the Styblo rule. Bull WHO. 2008;86:20–6.
Gao L, Lu W, Bai L, Wang X, Xu J, Catanzaro A, et al. Latent tuberculosis infection in rural China: baseline results of a population-based, multicentre, prospective cohort study. Lancet Inf Dis. 2015;14:71085.
Zhang X, Jia H, Liu F, Pan L, Xing A, Gu S, et al. Prevalence and risk factors for latent tuberculosis infection among health care workers in China: a cross-sectional study. PLoS One. 2013;8:e66412.
Leung CC, Rieder HL, Lange C, Yew WW. Treatment of latent infection with Mycobacterium tuberculosis: update 2010. ESRJ. 2011;37:690–711.
Robertson BD, Altmann D, Barry C, Bishai B, Cole S, Dick T, et al. Detection and treatment of subclinical tuberculosis. Tuberculosis. 2012;92:447–52.
Raviglione MC, Sudre P, Rieder HL, Spinaci S, Kochi A. Secular trends of tuberculosis in Western Europe. Bull WHO. 1993;71:297–306.
Cegielski JP, McMurray DN. The relationship between malnutrition and tuberculosis: evidence from studies in humans and experimental animals. IJTLD. 2004;8:286–98.
Wu P, Cowling BJ, Schooling CM, Wong IOL, Johnston JM, Leung CC, et al. Age-period cohort analysis of tuberculosis notifications in Hong Kong from 1961 to 2005. Thorax. 2008;63:312–6.
Vynnycky E, Fine EM. Interpreting the decline in tuberculosis: the role of secular trends in effective contact. Int J Epi. 1999;28:327–34.
Oxlade O, Schwartzman K, Behr MA, Benedetti A, Pai M, Heyman J, et al. Global tuberculosis trends: a reflection of changes in tuberculosis control or in population health? IJTLD. 2009;13:1238–46.
Shi L. Health care in China: a rural urban comparison after the socioeconomic reforms. Bull WHO. 1993;71:723–36.
Long Q, Li Y, Wang Y, Yue Y, Tang C, Tang S, et al. Barriers to accessing TB diagnosis for rural-to-urban migrants with chronic cough in Chongqing, China: a mixed methods study. BMC Health Serv Research. 2008;8:202.
Read JM, Lessler J, Riley S, Wang S, Tan LJ, Kwok KO, et al. Social mixing patterns in rural and urban areas of southern China. Proc Royal Society B. 2014;281:20140268.
Wallinga J, Teunis P, Kretzschmar M. Using data on social contacts to estimate age-specific transmission parameters for respiratory-spread infectious agents. Am J Epidmiol. 2006;164:936–44.
Eames KTD, Tilston NL, Brooks-Pollock E, Edmunds WJ. Measured dynamic social contact patterns explain the spread of H1N1v influenza. PloS Comp Bio. 2012;8:e1002425.
Murray CJL et al. Global, regional, and national incidence and mortality for HIV, tuberculosis and malaria during 1990–2013: a systemic analysis for the Global Burden of Disease Study 2013. Lancet. 2014;384:1005–70.
Daley CL, Small PM, Schecter GF, Schoolnik GK, McAdam RA, Jacobs Jr WR, et al. An outbreak of tuberculosis with accelerated progression among persons infected with the human immunodeficiency virus - an analysis using RFLP. NEJM. 1992;326:231–5.
Shea KM, Kammerer JS, Winston CA, Navin TR, Horsburgh Jr CR. Estimated rate of reactivation of latent tuberculosis infection in the United States, overall and by population subgroup. Am J Epidemiol. 2014;179:216–25.
Kline SE, Hedemark LL, Davies SF. Outbreak of tuberculosis among regular patrons of a neighborhood bar. NEJM. 1995;333:222–7.
Valway SE, Sanchez MPC, Shinnick TF, Orme I, Agerton T, Hoy D, et al. An outbreak involving extensive transmission of a virulent strain of Mycobacterium tuberculosis. NEJM. 1998;338:633–9.
Kiers A, Drost AP, van Soolingen D, Veen J. Use of DNA fingerprinting in international source case finding during a large outbreak of tuberculosis in The Netherlands. IJTLD. 1997;1:239–47.
Vynnycky E, Fine PE. The annual risk of infection with Mycobacterium tuberculosis in England and Wales since 1901. Int J Tuberc Lung Dis. 1997;1:389–96.
Nelson LJ, Schneider E, Wells CD, Moore M. Epidemiology of childhood tuberculosis in the United States, 1993–2001: the need for continued vigilance. Pediatrics. 2004;114:333–41.
Starke JR. Pediatric tuberculosis: time for a new approach. Tuberculosis. 2003;83:208–12.
Ozvaran MK, Baran R, Tor M, Dilek I, Demiryontar D, Arinc S, et al. Extrapulmonary tuberculosis in non-human immunodeficiency virus-infected adults in an endemic region. Ann Thorac Med. 2007;2:118–21.
Yang Z, Kong Y, Wilson F, Foxman B, Fowler AH, Marrs CF, et al. Identification of risk factors for extrapulmonary tuberculosis. Clin Inf Dis. 2004;38:199–205.
CDC: Tuberculosis infection control. http://www.cdc.gov/tb/education/corecurr/pdf/chapter7.pdf.
Verver, Warren RM, Munch Z, Vynnycky E, van Helden PD, Richardson M, et al. Transmission of tuberculosis in a high incidence urban community in South Africa. Int J Epidemiol. 2004;33:351–7.
The authors would like to thank Bill and Melinda Gates for their active support of this work and their sponsorship through the Global Good Fund. Productive discussions with colleagues at the Institute for Disease Modeling, including Hil Lyons and Stewart Chang, and at the Tuberculosis Program of the Bill & Melinda Gates Foundation are also greatly appreciated.
The authors declare that they have no competing interests.
GHH, DPC, DJK, and LW designed and conceptualized the study. GHH and DJK wrote the code, conducted analyses, and wrote the draft of the manuscript. DPC, BGW, PAE, RL, and LW reviewed and provided feedback in response to the manuscript draft. All authors read and approved the final manuscript.
Supplementary material describing model structure, calibration and calculation of the credible interval.