A user-friendly, open-source tool to project impact and cost of diagnostic tests for tuberculosis

  1. David W Dowdy  Is a corresponding author
  2. Jason R Andrews
  3. Peter J Dodd
  4. Robert H Gilman
  1. Johns Hopkins Bloomberg School of Public Health, United States
  2. Johns Hopkins University, United States
  3. Massachusetts General Hospital, United States
  4. University of Sheffield, United Kingdom

Abstract

Most models of infectious diseases, including tuberculosis (TB), do not provide results customized to local conditions. We created a dynamic transmission model to project TB incidence, TB mortality, multidrug-resistant (MDR) TB prevalence, and incremental costs over 5 years after scale-up of nine alternative diagnostic strategies. A corresponding web-based interface allows users to specify local costs and epidemiology. In settings with little capacity for up-front investment, same-day microscopy had the greatest impact on TB incidence and became cost-saving within 5 years if delivered at $10/test. With greater initial investment, population-level scale-up of Xpert MTB/RIF or microcolony-based culture often averted 10 times more TB cases than narrowly-targeted strategies, at minimal incremental long-term cost. Xpert for smear-positive TB had reasonable impact on MDR-TB incidence, but at substantial price and little impact on overall TB incidence and mortality. This user-friendly modeling framework improves decision-makers' ability to evaluate the local impact of TB diagnostic strategies.

https://doi.org/10.7554/eLife.02565.001

eLife digest

Tuberculosis is an infectious bacterial disease caused predominantly by the microorganism Mycobacterium tuberculosis. Although the number of deaths from tuberculosis has been falling in recent years, the disease still kills more than 1 million people every year, mainly in developing countries. Tuberculosis can be treated with antibiotics, but the emergence of bacteria that are resistant to existing drugs is threatening efforts to eradicate the disease.

Preventing the spread of tuberculosis is heavily dependent on accurate diagnosis of individuals with the disease. This is challenging because the initial symptoms are often mild, usually just a cough, which means that someone can spread the disease to many others over a period of several months before the symptoms become worse—fever, night sweats, and weight loss—and they realize that they are sick. Multiple diagnostic strategies are available, from the relatively low-tech—examining sputum samples under a microscope to detect tuberculosis bacteria—to more sophisticated tests that can detect bacterial DNA and determine whether the bacteria are drug-resistant in less than 2 hr.

Choosing which diagnostic strategy to adopt can be challenging because the optimal solution in a region will depend on the specific local conditions. To overcome this problem, Dowdy et al. have developed a computer program that enables decision-makers to input four key parameters that describe the tuberculosis situation in their region, and to obtain 5-year projections of the rate of new infections, mortality, and total costs likely to result from adopting any of nine different diagnostic strategies.

The four parameters are the number of new cases of tuberculosis each year (incidence), the proportion of new cases that are multi-drug resistant, the proportion of the adult population that has HIV, and the local costs of various diagnostic techniques and treatments. Since the entire computer program is written in a freely available open-source programming language (Python), any user can tweak these parameters to provide a more precise fit to their own region. Alternatively, the standard version of the program can be run directly from a website without any need to interact with computer code.

This model is the first to enable local decision-makers to evaluate the impact of different diagnostic strategies for tuberculosis under the conditions specific to their region. The model predicts, for example, that in areas where there is little money available for up-front investment, same-day microscopy analysis of sputum samples and starting patients on treatment is the most cost-effective strategy for reducing the rate of new infections. Given the wide variation in conditions within even small geographical areas, this more flexible approach should lead to the more efficient use of resources and may, ultimately, help to reduce the spread of tuberculosis.

https://doi.org/10.7554/eLife.02565.002

Introduction

Infectious disease transmission models are important tools for translating the best current knowledge of the natural history and epidemiology of infectious diseases into projections of epidemiological impact (e.g., incidence, mortality) and costs under alternative strategies for disease control (Garnett et al., 2011). Currently, most published transmission models are either loosely calibrated to reflect global/regional outcomes or more tightly fit to specific epidemiological settings; in either case, model results may be difficult for local decision-makers in the majority of public health settings to utilize. Simplified models designed for in-country use by decision-makers, most notably the Spectrum suite of models supported by the Futures Institute (Stover et al., 2010), have been used to inform decision-making in the fields of reproductive health and human immunodeficiency virus (HIV) for over a decade (Stover, 2004). Estimates from the Spectrum models are now routinely incorporated into official global and country-level estimates of HIV disease burden (Brown et al., 2010) and intervention impact (Farnham et al., 2013). Other simplified models are readily available for impact projections related to non-infectious diseases, where transmission assumptions are less important (Betz Brown et al., 2000; Walker et al., 2013). However, to date, simple, user-friendly transmission models have not been widely used for decision-making related to many infectious diseases other than HIV. Diagnosis of active tuberculosis (TB) is an example of a public health intervention for which transmission models may provide guidance on both global (Dowdy et al., 2006; Abu-Raddad et al., 2009) and country-specific levels (Menzies et al., 2012). Specifically, an unprecedented number of new diagnostic strategies for active TB are now recommended by the World Health Organization (WHO), including same-day microscopy (World Health Organization, 2011), microcolony-based culture techniques (Leung et al., 2012) such as the microscopic-observation drug-susceptibility (MODS) assay (Moore et al., 2006), line-probe assays for drug susceptibility testing (Bwanga et al., 2009), and Xpert MTB/RIF (‘Xpert’), a molecular assay capable of providing results (including rifampin resistance) in 90 min with minimal human resource requirements (Boehme et al., 2010, 2011). TB program decision-makers must repeatedly determine when to invest in scaling up a novel diagnostic test, which test(s) to promote, and whether the implementation strategy should differ by epidemiological situation (Cobelens et al., 2012). Without transmission models to provide locally relevant estimates of cost and impact under alternative implementation strategies, such decisions will be made without systematically considering the implications of available scientific evidence.

To aid in this decision-making process, we created a flexible, simple modeling tool that allows non-expert users to define their local situation according to three key epidemiological parameters (TB incidence, proportion of new TB cases that are multidrug-resistant [MDR], and adult human immunodeficiency virus [HIV] prevalence) and local unit costs of TB diagnosis and treatment. This tool then incorporates those estimates into a combined decision analysis-transmission framework to generate 5-year projections of TB incidence, mortality, and control costs for nine diagnostic strategies (Figures 1 and 2). These strategies are:

  1. ‘Baseline’: Sputum smear microscopy for each diagnostic attempt, with liquid-media TB culture only to evaluate smear-positive cases with a history of previous TB treatment for drug resistance. (Cultures in all scenarios trigger drug-susceptibility testing if positive.)

  2. ‘TB culture if previously treated’: Sputum smear microscopy used for patients without a history of TB treatment; smear plus liquid-media culture used to diagnose TB in any previously treated individual with symptoms (regardless of smear status).

  3. ‘Xpert if HIV-positive’: Xpert MTB/RIF for HIV-infected patients only, with a positive test for rifampin resistance triggering treatment for MDR-TB. Xpert is assumed to be deployed at the district level, such that results cannot generally be provided during the same clinical encounter (Lawn et al., 2012). This strategy is conceived as a ‘best-case’ scenario for HIV-targeted TB testing: if individuals unaware of their HIV status are not tested with Xpert, this strategy will overestimate effectiveness, and if those unaware of their status are tested, it will underestimate costs.

  4. ‘Xpert if Smear-Positive’: Xpert MTB/RIF for smear-positive patients only (i.e., for rapid DST), with a positive test for rifampin resistance triggering treatment for MDR-TB.

  5. ‘Xpert for All’: Xpert MTB/RIF for all patients.

  6. ‘Xpert with Culture DST Confirmation’: Xpert MTB/RIF for all patients, but treatment for MDR-TB only initiated if rifampin resistance is confirmed by culture.

  7. ‘MODS/TLA’: Sputum smear, plus microcolony-based TB culture (e.g., MODS or thin-layer agar, TLA) for all patients.

  8. ‘Same-Day Microscopy’: Double the per-test cost of sputum smear microscopy, in exchange for the ability to provide results to patients in the same clinical encounter (e.g., with peripheral, unbatched reading of sputum smears).

  9. ‘Same-Day Xpert’: Double the per-test cost of Xpert MTB/RIF, in exchange for the ability to provide results to patients in the same clinical encounter (e.g., peripheral deployment, with greater costs reflecting lower volume per machine [Vassall et al., 2011]).

Overview of user-friendly model.

Users are asked, via open-source computer script or Web interface, to select one of the nine diagnostic strategies and to provide unit costs and three basic epidemiological parameters (TB incidence, MDR-TB prevalence among new cases, and adult HIV prevalence). The selected diagnostic strategy is used to populate a decision tree that calculates (a) the probability of missed diagnosis, unsuccessful treatment, and successful treatment, (b) costs, and (c) diagnostic delays. These outputs depend on patients' TB (yes/no, and drug susceptibility status), HIV, and TB treatment history status. The selected epidemiological parameters are then used to populate a dynamic transmission model, creating a steady-state population that reflects local TB epidemiology. The decision tree—which inputs user-defined unit costs—is then incorporated into the transmission model to project outcomes under the selected diagnostic scenario. Users can sequentially select multiple diagnostic scenarios for comparison, and the computer script (though not the Web interface) allows users to manipulate input parameters at their discretion.

https://doi.org/10.7554/eLife.02565.003
Transmission model of TB diagnosis.

Boxes represent sub-populations in the model, and arrows represent rates of movement between those sub-populations. Parallel structures exist for: (a) HIV-infected vs HIV-uninfected; (b) never-treated vs previously treated (for TB); and (c) among TB-infected individuals, drug-susceptible vs isoniazid-monoresistant vs rifampin-resistant (including MDR). ‘Pre-diagnostic’ TB refers to individuals who are infectious but have not yet begun to seek care. Mortality occurs from all sub-populations (not shown), and at a higher rate among those with HIV and active TB.

https://doi.org/10.7554/eLife.02565.004

For purposes of illustration, we evaluated each of these diagnostic strategies in four emblematic epidemiological settings, defined by TB incidence, MDR-TB prevalence among new cases, and adult HIV prevalence:

  1. ‘Reference/High-Incidence Setting’ (e.g., Southeast Asia): TB incidence 250 per 100,000/year (twice the global incidence [World Health Organization, 2012]), MDR-TB prevalence of 3.7% in new TB cases (the estimated global prevalence [World Health Organization, 2012]), adult HIV prevalence of 0·83% (the estimated global prevalence [UNAIDS, 2012]);

  2. ‘Low-Incidence Setting’ (e.g., United States, Western Europe): TB incidence at entry of 8.9 per 100,000/year and declining, with similar MDR-TB (as a proportion of new cases) and HIV prevalence as above;

  3. ‘High MDR Setting’ (e.g., former Soviet Union): TB incidence 100 per 100,000/year, MDR-TB prevalence of 10% in new TB cases, adult HIV prevalence of 0.83%; and

  4. ‘High HIV Setting’ (e.g., sub-Saharan Africa): Adult HIV prevalence of 20%, TB incidence 500 per 100,000, and MDR-TB prevalence among new cases of 3.7%.

In each setting, we used a uniform set of costs for purposes of comparison (Table 1). Although these four settings form the basis of the results presented here, decision-makers can use the open-source model program (included as Supplementary file 1, written in the open-source programming language Python, Version 2.7, www.python.org) to re-define any parameter in Table 1 according to their best local knowledge; a manual for doing so is also included as Supplementary file 2. Capacity to create non-equilibrium settings (e.g., declining TB incidence, increasing MDR-TB) is included. We also provide a web-based interface (flexdx.modeltb.org) that allows users to input local TB incidence, MDR-TB prevalence, HIV prevalence, and unit costs; this interface provides customized results without the requirement to manipulate programming code. The program corresponding to the web version is also available on a public repository (https://github.com/JJPennington/FlexDx-TB-Web-Django).

Table 1

Model input parameters*

https://doi.org/10.7554/eLife.02565.005
ParameterValueReference(s)/Rationale
TB and HIV Transmission
Transmission rate, per smear-positive/highly infectious person-yearCalibrated to user-defined TB incidence
Proportional reduction in per-case transmission rate, MDR-TBCalibrated to user-defined MDR-TB prevalence
Proportional reduction in fitness, isoniazid-monoresistant TB25% of MDR-TB reductionAssumption
HIV incidence rate, per yearCalibrated to user-defined HIV prevalence
Relative transmission rate from smear-negative/less infectious TB0.22(Behr et al., 1999)
Proportion of pulmonary TB that is smear-positive/highly infectious
 HIV-negative0.63(Steingart et al., 2006a; Steingart et al., 2006b)
 HIV-infected0.50(Getahun et al., 2007)
TB Progression
Endogenous reactivation rate
 HIV-negative0.0005/year(Horsburgh et al., 2010)
 HIV-infected0.05/year(Antonucci et al., 1995)
Proportion of recent infections resulting in rapid progression
 HIV-negative0.14(Vynnycky and Fine, 1997; Dye et al., 1998)
 HIV-infected0.470.75 without ART, (Daley et al., 1992)
75% reduction if on ART, (Williams et al., 2010) 50% ART coverage
Reduction in TB rapid progression probability due to latent TB infection (HIV-negative only)0.79(Andrews et al., 2012)
TB Mortality and Resolution
Life expectancy at age 1555 years(World Bank, 2012)
Annual mortality from HIV0.05/year(UNAIDS, 2012)
Annual mortality from TB
 HIV-negative, smear-positive/highly infectious0.23/year(Tiemersma et al., 2011)
 HIV-negative, smear-negative/less infectious0.07/year(Tiemersma et al., 2011)
 HIV-infected1.0/year(Corbett et al., 2003; Corbett et al., 2007; Wood et al., 2007)
Rate of spontaneous TB resolution (HIV-negative only)
 Smear-positive/highly infectious0.1/year(Tiemersma et al., 2011)
 Smear-negative/less infectious0.27/year(Tiemersma et al., 2011)
TB Treatment Outcomes and Emergence of Drug Resistance
Probability of failure or relapse (within 1 year)
 Drug-susceptible0.04(World Health Organization, 2012)
 INH-monoresistant, first-line therapy0.21(Menzies et al., 2009b)
 INH-monoresistant, retreatment or 2nd-line0.16(Menzies et al., 2009b)
 MDR-TB, first-line or retreatment0.50(Espinal et al., 2000)
 MDR-TB, second-line therapy0.30(World Health Organization, 2010)
Proportion of one-year recurrence due to failure
 Drug-susceptible0.14(Lew et al., 2008)
 INH-monoresistant0.33
 MDR-TB0.56
Probability of acquired drug resistance (per treatment course)
 Susceptible becoming INH-monoresistant0.001(Menzies et al., 2009a; Menzies et al., 2009b)
 Susceptible becoming MDR-TB0.002
 INH-monoresistant becoming MDR-TB0.045
 If treated with 2 effective drugs for >6 mos0.017
Behavioral Parameters
Infectious months before starting to seek care
 HIV-negative9 months(Dowdy et al., 2013)
 HIV-infected1 month(Corbett et al., 2004)
Diagnostic frequency while seeking care5.0/year(Storla et al., 2008; Sreeramareddy et al., 2009)
Probability of treatment in a TB patient whose microbiological test is negative0.25(Wilkinson et al., 2000; Dowdy et al., 2008)
Loss to follow-up between diagnostic presentation and treatment initiation
 Sputum smear or GXP (not same-day)0.15(MacPherson et al., 2014)
 Culture (microcolony or commercial liquid)0.25(Dowdy et al., 2008)
 Same-day diagnosis0Assumption
Diagnostic Accuracy
Sensitivity for smear-negative/less-infectious TB
 Sputum smear microscopy0
 Xpert MTB/RIF0.72(Brownell et al., 2012)
 Culture (microcolony or commercial liquid)0.85(Cruciani et al., 2004; Leung et al., 2012)
Specificity for TB(Steingart et al., 2006; Boehme et al., 2011; Leung et al., 2012)
 Sputum smear microscopy0.98
 Xpert MTB/RIF0.98
 Microcolony culture0.98
Sensitivity for drug resistance (if TB detected)
 Microcolony culture (rifampin and isoniazid)0.98(Minion et al., 2010)
 Xpert MTB/RIF (rifampin only)0.94(Boehme et al., 2011)
Specificity for drug resistance (if TB detected)
 Microcolony culture (isoniazid)0.96(Minion et al., 2010)
 Microcolony culture (rifampin)0.99(Minion et al., 2010)
 Xpert MTB/RIF (rifampin)0.98(Boehme et al., 2011)
Diagnostic Delay and non-TB Care-Seeking
Days from presentation to treatment initiation
 Sputum smear or Xpert MTB/RIF7 daysAssume 1 week
 Microcolony or commercial liquid culture30 days(Boehme et al., 2011)
Months of therapy before a failing regimen will be changed, or before default and recurrence6 monthsAssumption
Annual rate of diagnostic evaluation for TB, among people who do not have active TB0.01/year10% of suspects have TB, high-incidence setting
Cost Parameters (user-defined; values below for comparison purposes only)
Per-patient cost of TB therapy
 First-lineUS$500User-defined
 RetreatmentUS$1000(Vassall et al., 2011)
 Second-line/MDRUS$5000(Vassall et al., 2011)
 Outpatient visit (diagnosis or follow-up)US$10(Vassall et al., 2011)
Per-test cost:
 Sputum smearUS$2(Vassall et al., 2011)
 Same-day sputum smearUS$10Assumption
 Xpert MTB/RIFUS$15(Vassall et al., 2011)
 Same-day Xpert MTB/RIFUS$30Assumption
 Microcolony culture (with DST)US$5(Solari et al., 2011)
 Commercial liquid-media cultureUS$20(Vassall et al., 2011)
 Commercial liquid-media culture + DSTUS$40(Vassall et al., 2011)
  1. *

    In the actual model program (Supplementary file 1), users can change any parameter based on local values.

  2. For reference, the transmission rate (in infections per person-year during diagnosis-seeking active TB) is 36.9 in the reference scenario, 14.0 in the low-incidence scenario, 25.4 in the high MDR scenario, and 12.9 in the high HIV scenario. Corresponding proportional reductions in MDR-TB transmission rate are 0.23, 0.23, 0.21, and 0.19; and HIV incidence estimates (per 1000 adult person-years) are 0.7, 0.6, 0.6, and 18.9.

Results

Model validation

To validate the model, we compared selected model outcomes to published global estimates. In the high-incidence setting, our model estimated TB mortality at 14% of incidence (95% uncertainty range: 7–20%, WHO global estimate 14% [World Health Organization, 2012]), HIV-associated TB at 13% of all incident TB (95% uncertainty range: 4–14%, global estimate 13% [World Health Organization, 2012]), previously treated cases at 13% of all incident cases (95% uncertainty range: 9–34%, global estimate 14% [World Health Organization, 2012]), and duration of TB disease at 1.2 years (95% uncertainty range 0.8–2.1, global estimate 1.4 years [World Health Organization, 2012]). Our model estimated that MDR-TB prevalence in previously treated cases was 15.4% (WHO estimate 20% [World Health Organization, 2012]), but unlike our model, WHO notifications often count failure and recurrence after default (in the same person) as two separate cases. At steady-state in the model, 80% of incident TB was due to recent infection rather than reactivation.

Comparison of diagnostic strategies in the high incidence scenario

Figure 3 shows the projected incremental 5-year cost and impact of each of nine selected diagnostic strategies (described in greater detail in the ‘Materials and methods’ section), in the high-incidence scenario. In general, both the cost and impact of targeted strategies (culture for retreatment, Xpert for HIV-positive, Xpert for smear-positive) on incidence were small relative to broader diagnostic strategies (Xpert for all, MODS/TLA, same-day Xpert). The incremental cost-effectiveness, in terms of cost per case averted (i.e., slope of the line from origin to each point in Figure 3), was similar across all strategies with the exception of same-day smear, which was cost-saving and had greater effectiveness than the baseline. Same-day microscopy remained cost-saving by year 5 as long as same-day results could be provided at less than five times the per-smear cost of routine results (i.e., <$10/test). Among the targeted strategies, Xpert for smear-positives was the most expensive but had the greatest impact on MDR-TB cases averted, whereas Xpert for HIV-infected individuals and culture for retreatment cases offered smaller gains at lower cost. Among the broad strategies, culture confirmation of rifampin-resistant tests on Xpert saved costs with little reduction in effectiveness. There was little difference between culture-confirmed Xpert and MODS/TLA, and same-day Xpert was the most expensive and most effective strategy.

Incremental 5-year cost and impact of TB diagnostic strategies, high-incidence setting.

Shown are cumulative projected 5-year costs and impact (averted TB cases [panel A] or MDR-TB cases [panel B]) of each diagnostic strategy described in the Introduction, incremental to the baseline strategy, per 100,000 population. Strategies with greater impact appear to the right on the x-axis; more costly strategies appear higher on the y-axis. The same-day smear strategy is cost-saving but shown at an incremental cost of $0 for simplicity.

https://doi.org/10.7554/eLife.02565.006

Comparison of diagnostic strategies across settings

The projected impacts and costs of the nine diagnostic strategies relative to the baseline strategy, in each of four selected epidemiological settings, are shown in Figure 4. In all four settings, the ranking of diagnostic strategies remained similar for all outcomes, although the cost of broader diagnostic strategies relative to targeted strategies fell substantially over 5 years in higher-incidence settings as the broader strategies generated declines in TB incidence. In the low-incidence setting, where a higher proportion of TB treatments are false-positive and more incident TB is also due to reactivation (60% of all new cases), the relative cost of improved TB diagnosis was the highest, while the relative impact was the least. In the high HIV setting, the impact of diagnostic interventions on TB incidence was diminished relative to the high-incidence setting, though the impact on TB mortality was similar. Additionally in this setting, the Xpert for HIV-positive strategy was substantially more costly, but also more effective, than in other settings.

Relative impact of diagnostic strategies in emblematic settings.

Shown are projected changes in TB incidence, MDR-TB incidence, TB mortality, and costs (in Year 1 and Year 5 after immediate implementation), relative to baseline (Strategy 1) after implementing each of the diagnostic strategies described in the text. Epidemiological outcomes are measured at the end of Year 5. Panel A (high incidence) shows a setting with TB incidence of 250 per 100,000/year, stable MDR-TB prevalence of 3.7% among new cases, adult HIV prevalence of 0.83%, and cost of $500 to treat one case of TB with first-line therapy. In panel B (low incidence), the TB incidence is reduced to 8.3 per 100,000/year (implemented by gradual decline in incidence over 50 years). In panel C (high MDR), MDR-TB prevalence among new cases is set at 3.7% in the beginning of year 1, increasing to 10.7% by the end of year 5. In panel D (high HIV), adult HIV prevalence is set to 20% and TB incidence is set to 500 per 100,000/year.

https://doi.org/10.7554/eLife.02565.007

Sensitivity analyses

The impact of the ‘Xpert for all’ strategy on TB incidence (selected a priori as the primary outcome for sensitivity analysis) was most sensitive to three parameters, both in terms of absolute effects on impact estimates and partial rank correlation coefficients. These three parameters were: (1) the proportion of TB patients who would be empirically treated even if microbiologic testing yielded a negative result (‘empiric treatment proportion’), (2) duration of infectiousness before seeking care (‘pre-diagnostic delay’), and (3) rate of reactivation from latent infection to active disease (‘reactivation rate’). For this latter parameter, in-depth investigation revealed that the key determinant was not the reactivation rate per se, but rather the proportion of active TB representing recent vs remote infection. If the empiric treatment proportion was increased from 25% to 37.5%, the projected reduction in TB incidence fell from 20% to 13%. By contrast, when only 12.5% of false-negative patients were started on therapy, ‘Xpert for all’ achieved a 31% reduction in incidence. Corresponding reductions in incidence with ‘Xpert for all’ (20% at baseline) included: 27% if pre-diagnostic delay was shortened from 9 months to 4.5 months, 14% if pre-diagnostic delay was lengthened to 13.5 months, 15% if the reactivation rate was doubled from 0.05% to 0.1% per year, and 23% if reactivation was halved to 0.025% per year. The projected 5-year reduction in incidence for this strategy did not fall outside the range of 12–25% under variation of any other model parameter, up to 50% of that parameter's baseline value (Table 1).

Both costs and incremental costs were more sensitive to the cost of TB treatment than the cost of diagnostics. For example, doubling the cost of first-line therapy (from $500 to $1000 per person) augmented the incremental year 1 costs for ‘Xpert for all’ from a 57% increase over baseline to a 77% increase, and doubling the cost of MDR therapy (from $5000 to $10,000 per person) generated an 86% increase in incremental year 1 costs, whereas doubling the unit cost of Xpert (from $15 to $30) only resulted in a 72% increase. Unit costs other than those for first-line treatment, MDR treatment, and the diagnostic modality under study in each scenario were not important determinants of incremental costs.

Discussion

To date, most transmission models of infectious disease control interventions present results that are not directly usable by decision makers because they are not customizable to local conditions. We present a flexible, user-friendly model of TB diagnosis and transmission that allows users without modeling expertise to define an epidemiologic setting (according to TB incidence, MDR-TB prevalence, and HIV prevalence) and unit costs, evaluating various diagnostic strategies in that setting in terms of population-level costs and impact. While this model cannot precisely replicate the epidemiological situation in any given location, it applies a standardized methodology across a wide range of settings, thereby illustrating important interactions between epidemiological parameters and projected impact. We provide both a web interface for rapid calculations and full model code whereby users can change any model parameter for ‘personalized’ sensitivity analysis. Our model results suggest that the rank-ordering of diagnostic strategies may be relatively stable across epidemiological settings, but that the actual population-level costs and impact differ dramatically. This model also serves as an example of how epidemiologists can provide decision-makers across a wide variety of local settings with rapid access to customizable ‘first-pass’ projections of cost and impact from transmission models without the need to construct tightly fitted models to represent all epidemiological settings.

Our results provide important guidance to TB decision-makers. Specifically, in settings where little additional up-front investment is possible (<25% increase in TB control budget), same-day microscopy has the greatest impact on TB incidence with no increase in overall cost at 5 years, provided that same-day microscopy can be feasibly delivered for less than $10 per patient. For settings in which containing MDR-TB is the most important consideration, Xpert for smear-positives has the greatest effectiveness for this outcome, but at substantial price and very little impact on overall TB incidence and mortality. Xpert for HIV-positives and culture for retreatment cases both offer meaningful, albeit small, gains; the cost and impact of these strategies are the highest in HIV-endemic settings. In settings where more initial investment is possible (about 50% increase in TB control budget), broader scale-up of either Xpert or MODS/TLA for all TB diagnosis can offer substantially greater benefits, both in terms of reduced incidence (often 10 times more TB cases averted than with the more narrowly-targeted strategies above) and long-term costs (in that the incremental cost of these strategies declines greatly by Year 5). In general, culture confirmation of positive Xpert results before committing to a course of second-line therapy is preferred. Finally, where the greatest impact is sought, combination of Xpert for all plus infrastructure for same-day diagnosis achieves this aim in all settings, but at the highest cost.

Although this model represents a highly simplified framework, it compares well to other global estimates to which it was not fit (e.g., WHO estimates of TB mortality, previously treated TB, HIV-associated TB, and MDR-TB prevalence among previously treated cases). Its results are also similar to those of other mathematical models of TB diagnosis that are fit to specific locations. For example, Menzies et al. (2012) modeled scale-up of Xpert for people living with HIV in southern Africa, estimating a 6% reduction in incidence, 21% reduction in TB mortality, 25% reduction in MDR-TB incidence that declines to <10% over a longer time frame, and 40% increase in costs (not including costs of antiretroviral therapy). Our corresponding estimates for the Xpert for HIV-positive strategy in the high-HIV setting are 5% reduction in incidence, 25% reduction in TB mortality, 7% reduction in MDR-TB incidence, and increase in costs from 41% (year 1) to 28% (year 5). For purposes of deciding between alternative strategies of scaling up TB diagnostics, the estimates from these two models are likely to provide similar guidance.

For any such simplified transmission model to have an impact on decision-making, it is important to consider who the eventual users of such a model might be, and to develop the model in consultation with those groups. In this case, we identified a number of potential end-user organizations including technical assistance agencies (e.g., KNCV Tuberculosis Foundation), funding bodies (e.g., Global Fund for AIDS, Tuberculosis, and Malaria), non-governmental organizations (e.g., Medecins Sans Frontiers), National Tuberculosis Programs, and regional offices of the World Health Organization. We then invited representatives of these organizations to attend a 1-day workshop (April 2014, The Hague, The Netherlands) describing methods and challenges related to modeling of TB diagnostics in general, including this transmission model in particular. We developed ‘hands-on’ exercises for participants to better learn the model and solicited specific feedback, which is being incorporated into the model structure and web interface in ongoing fashion. We next intend to develop a series of informal ‘case studies’ whereby the use of this model for in-country decision-making can be demonstrated and disseminated.

As with any modeling analysis, this research has important limitations. In order to provide sufficient flexibility and generalizability, we make a number of strong assumptions that include a constant population, homogeneous mixing, no change in parameter values over time, and simplistic incorporation of HIV and drug resistance. Without making such simplifying assumptions, it is impossible to deliver a flexible modeling framework that can generate transparent, customizable, rapid results (i.e., without complex statistical fitting to each individual epidemiological scenario). This model can replicate user-defined TB incidence, MDR-TB prevalence, and HIV incidence but does not describe the breakdown of these values in key subpopulations (e.g., congregate settings or geographical ‘hotspots’), nor does it incorporate operational aspects of the TB diagnostic system in any one setting. Thus, this model does not ameliorate the need for more detailed models in settings where precise estimates are needed; rather it provides access to ‘first-pass’ estimates in settings (i.e., the vast majority of local decision-makers) where such tightly calibrated model projections are not available. This model focuses on transmission dynamics and thus does not include pediatric TB and extrapulmonary TB; these largely non-infectious disease manifestations remain very important components of the TB epidemic that are not captured here. We validate our model against global estimates and other models; ideally, further validation would be performed over time using field data across a wide variety of settings. Our model allows users to define three key epidemiological parameters (as well as other model assumptions within the program), but data to inform even these three parameters—as well as unit costs—are unlikely to be available on sub-national levels. As a result, users wishing to adapt the model to a smaller geographic scale will need to perform additional data-gathering exercises to inform these estimates if they wish to maximize the utility of the model.

Nevertheless, even if high-quality data are not available at the local level, this model allows decision-makers to estimate epidemiological and economic values according to reasonable assumptions (e.g., comparison of TB notification rates or budget line-items to those in the published literature or at the national level), vary those assumptions in real-time, and obtain corresponding projections of comparative impact that incorporate the best available current data on epidemiological or natural history parameters (e.g., TB progression and reactivation). Future efforts might also provide flexibility to specify operational characteristics (e.g., health system capacity) as well.

A final concern is that, by providing users the ability to specify TB incidence, MDR-TB prevalence, and adult HIV prevalence, a number of scenarios can be created (e.g., low TB incidence and very high adult HIV prevalence) that are not epidemiologically realistic. Although there is danger in allowing uninformed users to make projections for such scenarios (and the model will reject or alert users to highly implausible values), we believe that this risk is outweighed by the benefit of providing full flexibility to model epidemiological scenarios (e.g., sub-district level data) that will never be captured by a limited number of closely-calibrated TB transmission models.

In summary, we have created a flexible modeling framework that allows users without modeling expertise to generate simulated populations with locally relevant values for TB incidence, MDR-TB prevalence, adult HIV prevalence, and TB treatment costs. By comparing an array of diagnostic options across emblematic epidemiological scenarios, we provide guidance to decision-makers who seek to ascertain the optimal diagnostic strategy to achieve their selected disease control targets, and to do so using a standardized methodology. Success in the fight against infectious disease generally, and TB specifically, depends on our ability to place global knowledge in the hands of local decision-makers, enabling them to choose those interventions that are likely to have the greatest impact, given existing resources and local epidemiological realities. This flexible modeling framework of diagnostic interventions takes an important step in that direction.

Materials and methods

Transmission model

Request a detailed protocol

Using previously published models of TB diagnostics as a guide (Dye et al., 1998; Abu-Raddad et al., 2009; Menzies et al., 2012; Dowdy et al., 2013), we constructed a transmission model of TB using ordinary differential equations. This model categorizes patients according to HIV status (positive or negative), TB treatment status (never treated or previously treated), TB disease status (as shown in Figure 2), and among those who are infected with TB, drug resistance status (susceptible, isoniazid-monoresistant, and rifampin-resistant including MDR), and level of infectiousness (smear-negative/less-infectious and smear-positive/highly infectious). Individuals enter the model at age 15, and TB with no pulmonary component (i.e., not infectious) is not included. For purposes of transparent communication, we chose a population size of 100,000 (to match standard reporting of TB outcomes) and assumed a constant population with no net population growth or immigration/emigration. After constructing the transmission framework, we used decision analysis to estimate (a) the probability of each diagnostic outcome; (b) the diagnostic delay; and (c) the cost of TB diagnosis and treatment under each of the nine diagnostic strategies, assuming immediate implementation at the beginning of a given year (‘Year 1’). Two separate authors (DWD and PJD) independently coded the model; these models gave comparable results.

Model initiation

Request a detailed protocol

After setting probabilities and costs under each diagnostic strategy, we then created a flexible modeling structure capable of generating epidemiological scenarios as a function of three variables, which the user specifies: HIV prevalence (assuming a global mean level of antiretroviral therapy coverage), TB incidence, and MDR-TB prevalence among new TB cases. We accomplished this by allowing three key parameters to vary across model scenarios: annual HIV incidence, rate of TB transmission per smear-positive/highly infectious person-year, and relative per-case transmission rate of rifampin-resistant TB. We also allow users to specify all relevant unit costs for TB diagnosis and treatment; other model parameters were estimated from existing literature (Table 1). We then created a program that numerically generates a unique steady-state population meeting the user-defined values; this population serves as the baseline strategy (Strategy 1 above) at the beginning of the time frame under evaluation. The program has flexibility to create its steady-state population 50 years in the past, allowing it to replicate the protracted, slow declines in TB incidence as seen in many lower-incidence settings; this is done automatically for any scenario with a target TB incidence less than 50 per 100,000/year.

Model compartments

Request a detailed protocol

The mathematical model consists of the following TB compartments:

  • Uhp, Uninfected

  • Lhdp, Latently infected

  • Ehdp, Early active (infectious status i = 0)

  • Ahdip, Late active

  • Phdip, Active ‘pre-treatment’: diagnosis in progress, will lead to appropriate therapy

  • Ihdip, Active ‘inappropriate treatment’: receiving therapy that ends in default or failure

In these compartments, the subscript h refers to HIV status (h = 0 if HIV-uninfected, 1 if HIV-infected), d refers to drug resistance status (d = 0 if drug-susceptible, 1 if isoniazid [INH]-monoresistant, and 2 if multidrug-resistant [MDR]), i refers to infectious status (i = 0 if smear-negative/less infectious and 1 if smear-positive/highly infectious), and p refers to previous treatment status (p=0 if never treated, 1 if previously treated). Infectious status can be conceptualized as an individual's sputum smear status, if two smears were to be performed in a quality-assured laboratory at any given point in time.

Model structure

Request a detailed protocol

The model assumes an adult population of stable size with no immigration or emigration: the number of individuals entering the uninfected compartment U is set as equal to the number who die (whether from TB or other causes) from all other compartments. Pediatric and purely extrapulmonary TB are not explicitly considered because the diagnostic considerations for these manifestations are different. In the short-term, however, to the extent that these forms of TB are non-infectious and equally fatal as adult pulmonary TB, their effects on TB incidence and mortality may be approximated by dividing the model's projected incidence and mortality by (1 − proportion of TB that is not adult pulmonary), to obtain a new incidence/mortality estimate. Thus, if 20% of all TB in a given location is extrapulmonary or paediatric, the rough projected total TB incidence would be (projected TB incidence)/(0.8).

In this model, we consider latent TB infection to be asymptomatic and non-infectious, with a constant rate of reactivation and ongoing risk of exogenous reinfection leading to active TB; individuals successfully treated for TB are assumed to return to this compartment upon initiation of effective therapy (i.e., therapy that will result in completion, with no relapse for 2 years). Upon developing active TB, individuals enter a ‘pre-diagnostic’ phase that is characterized by a low level of infectiousness and mortality (equivalent to smear-negative TB) and during which individuals do not actively seek diagnosis. The duration of this phase (9 months) was selected a priori based on an existing model (Dowdy et al., 2013) in which the total duration of disease after incorporating this phase reflected the global ratio of prevalence to incidence, as estimated by the World Health Organization. We compared the total duration of disease to this ratio as part of our model validation and assumed that this ‘pre-diagnostic’ phase is much shorter for HIV-infected vs HIV-uninfected individuals. Upon completing this ‘pre-diagnostic’ phase, individuals progress to a diagnosis-seeking phase of active disease, which is characterized by separation into highly infectious (‘smear-positive’) and less infectious (‘smear-negative’) compartments. Among HIV-uninfected individuals, the highly infectious compartment also carries higher mortality risk. Diagnosis-seeking active TB implies active seeking of diagnosis at a defined rate; the probability that any single diagnostic attempt will result in effective therapy is calculated as a function of diagnostic sensitivity, probability of empiric therapy (i.e., without bacteriological confirmation), prior treatment status, and losses to follow-up, as described below. Each diagnostic attempt, if successful, leads either to effective therapy (which is initiated after a defined diagnostic delay) or to ineffective therapy (defined as leading to failure or default). In order to focus on differences between the nine selected strategies above in a tractable framework, we subsumed all other diagnostic tests and procedures (e.g., chest X-ray, antibiotic trials) as a probability of non-microbiologic diagnosis, without attempting to specify the associated cost or diagnostic delay. Once effective therapy is initiated, it is assumed to immediately render the individual non-infectious, with no residual risk of mortality. Upon initiation of ineffective therapy, individuals are assumed to remain infectious (at the ‘smear-negative’/less infectious level) for a defined period before either failing (followed by another round of therapy, which can be either appropriate/curative or ineffective) or default. Reasoning that default will occur, on average, at the midpoint between receipt of fully-ineffective and fully-effective therapy, half of defaulters are presumed to develop recurrent TB (which is assumed to occur immediately), while the other half return to the latent TB compartment (from which reactivation or reinfection remains possible). All individuals who relapse within 1 year are included as failures; thus, no specific parameter for relapse is incorporated. Individuals who are effectively treated, or whose disease is contained without therapy, return to the latent compartment following the convention of other TB models (Dye et al., 1998; Abu-Raddad et al., 2009).

Role of HIV coinfection

Request a detailed protocol

As the goal of this model is to focus on TB-related interventions, HIV is modeled as occurring at a defined annual incidence, calibrated to achieve a given user-defined prevalence at baseline. We do not explicitly model HIV infection in dynamic fashion (i.e., the HIV incidence rate does not depend on the number of HIV-infected individuals in the model). HIV infection is assumed to affect all parameters related to TB disease, including the level of immune protection afforded by latent infection (assumed zero if HIV-infected), mortality rate (increased), rate of reactivation from latent TB (increased), duration of ‘pre-diagnostic’ TB (decreased), and risk of ‘primary’ progression to active disease upon infection (increased). For purposes of maintaining a simple model structure, we do not explicitly model CD4 counts or antiretroviral therapy (ART), but instead assume the global average of ART coverage, as estimated by UNAIDS (2012). We weight all HIV-related parameters according to this estimated probability of ART receipt; this probability can be modified by users.

Drug resistance

Request a detailed protocol

We assume infection with TB strains of three different drug resistance levels: fully susceptible, INH-monoresistant, and MDR. Dual infection with multiple strains is not considered in this model, but superinfection (i.e., reinfection of a latently infected individual with a different strain, resulting in primary progression to disease with the reinfecting strain) is allowed, as is acquisition of resistance (i.e., change of strain from more susceptible to less susceptible) as a result of therapy.

Individuals without TB

Request a detailed protocol

A key consideration with any TB diagnostic strategy is the role of the diagnostic test as applied to individuals who do not have TB (i.e., specificity). We included this element by considering that a small proportion of the population without TB would present with TB-like symptoms each year. This proportion (selected such that 10% of individuals being evaluated in the high-incidence setting actually have underlying TB) remains constant across all scenarios, such that the pre-test probability of TB is higher in settings of high TB incidence, and declines over time as successful TB control strategies are employed. Individuals without TB who are (inappropriately) treated for TB incur costs of TB therapy and are also marked as ‘previously treated’ for purposes of diagnostic evaluation in the future.

Economic evaluation

Request a detailed protocol

Economic parameters are estimated using a unit-costing approach, whereby the unit cost of a TB diagnostic attempt and a TB treatment course (separately for first-line, category two, and second-line therapy) is enumerated under each scenario, and this cost is multiplied by the number of diagnostic attempts and treatment courses performed. For simplicity, we adopt the perspective of a TB control program for our costing; additional costs of HIV care and general health services (e.g., hospitalization) are not included. The decision tree below is used to estimate the cost per diagnostic attempt or treatment course, conditional on an individual's HIV, drug resistance, and previous treatment status. Individuals who default or die on therapy are assumed to incur half the cost of a treatment course. Given the short time horizon (5 years), the focus on costs and outcomes (rather than cost-effectiveness per se), and the desire to compare costs in year 1 and at the end of year 5 in equivalent terms, we did not discount future costs or outcomes for this analysis. All costs are reported in US dollars, assuming the year of costs that is specified by the user.

Decision analysis: probability of diagnostic success

Request a detailed protocol

Under each scenario, we use decision analysis to ascertain the following quantities related to each diagnostic attempt:

  • Probability of receiving successful treatment.

  • Probability of receiving ineffective treatment (resulting in failure or default, including the probability of acquired resistance).

  • Cost of treatment (conditional on whether treatment is successful or ineffective).

  • Cost of diagnosis.

  • Diagnostic delay incurred before treatment initiated.

These quantities are calculated conditional on each patient's infectious (smear) status, drug resistance status, HIV status, and prior treatment status. These probabilities are calculated for each diagnostic attempt, with the result fed back into the transmission model for purposes of appropriately allocating flows between compartments. Inputs into the decision model include the probabilities of failure/recurrence, probability of empiric therapy, loss to follow-up before treatment, diagnostic accuracy, diagnostic delay, and economic parameters as shown in Table 1. Outputs from the decision tree appear as parameters in the model, as described in Table 2 and the following equations.

Table 2

Model parameters and symbolic representations

https://doi.org/10.7554/eLife.02565.008
ParameterRepresentationBaseline value (see Table 1)
Transmission rate (transmission events per highly infectious person-year)βCalibrated to TB incidence
Proportional reduction in per-case transmission rate
 Drug-susceptible TBφ01.0
 Isoniazid-monoresistant TBφ125%* of φ2
 MDR-TBφ2Calibrated
HIV incidence rate, per yearθCalibrated to HIV prevalence
Relative transmission rate from smear-negative/less infectious TBζ0.22
Proportion of pulmonary TB that is smear-positive/highly infectious
 HIV-negativeψ00.63
 HIV-infectedψ10.50
Endogenous reactivation rate, per year
 HIV-negativeε00.005
 HIV-infectedε10.05
Proportion of recent infections resulting in rapid progression
 HIV-negativeπ00.14
 HIV-infectedπ10.47
Reduction in TB rapid progression probability due to latent TB infection
 HIV-negativeι0.79
 HIV-infectedNot included0
Baseline mortality rate, per yearμbl1/55 = 0.018
Additional HIV-related mortality rate, per yearμh0.05
Additional untreated TB-related mortality rate, per year
 HIV-negative, smear-positive/highly infectiousμt10.23
 HIV-negative, smear-negative/less infectiousμt00.07
 HIV-infectedμth1.0
Rate of spontaneous TB resolution, per year
 Smear-positive/highly infectiousν10.1
 Smear-negative/less infectiousν00.27
 HIV-infectedNot included0
Rate of starting diagnosis-seeking in active TB, per year
 HIV-negativeδe01.33 (9 months)
 HIV-infectedδe112 (1/month)
Rate of progression: ineffective therapy to repeat therapy (failure) or active TB (relapse), per yearδf6/12 = 0.5
Rate of diagnostic evaluation for TB, per year
 Late active TBInput into decision tree5.0
 No active TBτ00.01
Decision tree outputs (in addition to unit costs):Vary by intervention
 Successful diagnosis rate of late active TB, per yearσhdip
 Rate of movement from successful diagnosis to treatment (1/diagnostic delay), per yearρhdip
 Ineffective diagnosis rate of late active TB, per yearκhdip
 Rate of diagnosis and treatment leading to new resistance, per year
 susceptible to INH-monoresistantαsihip
 susceptible to MDRαsmhip
 INH-monoresistant to MDRαimhip
  1. *

    Calculated such that (1−φ1) = 0.25*(1−φ2).

Outcomes and sensitivity analysis

Request a detailed protocol

Our primary outcomes under each scenario were TB incidence, TB mortality, MDR-TB incidence, and incremental TB diagnostic and treatment costs (during year 1 and at the end of the 5-year period) relative to the baseline strategy. By giving flexibility to change model inputs, we provide users the ability to conduct any sensitivity analysis desired. However, for illustrative purposes, we also conducted a series of one-way sensitivity analyses in which each model parameter in Table 1 was varied by ±50% of its listed value (for proportions, 50% of the difference between the value and either zero or one). Our primary outcome for sensitivity analysis was the change in TB incidence, comparing the ‘Xpert for all’ strategy to baseline in the high incidence setting.

We also conducted multivariable uncertainty analyses by calculating partial rank correlation coefficients (Kendall, 1942) between each natural history parameter and the outcomes of TB incidence, TB mortality, and 5-year costs. In addition, we constructed 95% uncertainty intervals around our estimates of outcomes in each individual country by simultaneously varying each model parameter by ±10% over a uniform distribution and each target value (i.e., TB incidence, HIV prevalence, and MDR-TB prevalence) over its reported uncertainty range. In this fashion, we constructed 10,000 simulations using Latin Hypercube Sampling (McKay et al., 1979) and took 95% uncertainty ranges as the 2.5th and 97.5th percentiles of outcomes in these simulations; these ranges are provided in the web-based version of the model for each country.

Model equations

Rates of flow between compartments are governed by the system of ordinary differential equations listed in Equations 1–6.

We first define the forces of infection (according to resistance status) and total mortality for simplicity.

Forces of infection d)

Request a detailed protocol
λd(t)=[βN(t)]*φd*{ζ*h,p[Ehdp(t)+Ahd0p(t)+Phd0p(t)+Ihd0p(t)]+h,p[Ahd1p(t)+Phd1p(t)]}

Thus, TB infection is modelled as a density-dependent process, a function of the transmission rate (β), total number of individuals in the population N(t), number of individuals with ‘less infectious’ TB (E, early active; Ahd0p, late active smear-negative; Phd0p, ‘pre-treatment’ smear-negative; and I, active on ineffective treatment) weighted by the relative transmission rate ζ, and the number of individuals with fully infectious TB (Ahd1p, late active smear-positive; Phd1p, ‘pre-treatment’ smear-positive). Three separate forces of infection are calculated for the three separate strains of drug resistance, with relative transmission weights of φd. We use λtot to denote the sum of these three forces.

Total mortality (μtot)

Request a detailed protocol
μtot(t)=μbl*N(t)+μh*p{U1p(t)+d[L1dp(t)+E1dp(t)+i(A1dip(t)+P1dip(t)+I1dip(t))]}+μt0*p,d{E0dp(t)+A0d0p(t)+P0d0p(t)+i(I0dip(t))}+μt1*p,h,d{A0d1p(t)+P0d1p(t)}+μth*p,d{E1dp(t)+i(A1dip(t)+P1dip(t)+I1dip(t))}

Thus, total mortality is the sum of:

  • baseline mortality μbl (experienced by all individuals N),

  • HIV-associated mortality μh (experienced by all individuals with HIV, h = 1),

  • ‘less infectious’ (i.e., lower) TB-associated mortality μt0 (experienced by all HIV-uninfected individuals with early active TB, E, smear-negative active and ‘pre-treatment’ TB, A0d0p and P0d0p, and ineffectively treated TB, I0dip),

  • ‘highly infectious’ (i.e., higher) TB-associated mortality μt1 (experienced by all HIV-uninfected individuals with smear-positive active and ‘pre-treatment’ TB, A0d0p and P0d0p), and

  • TB/HIV–associated mortality μth (experienced by all HIV-infected individuals with any form of TB).

Uninfected compartment (U)

Request a detailed protocol
(1) dUhp(t)dt=Ih=0,p=0*μtot(t)[λ0(t)+λ1(t)+λ2(t)+μbl(t)]*Uhp(t)Ih=0*[θ*U0p(t)]+Ih=1*[θ*U0p(t)μh*U1p(t)]Ip=0*[τ0*(1sh)*Uh0(t)]+Ip=1*[τ0*(1sh)*Uh0(t)]

where μtot is the sum of all mortality (to maintain a constant population), λd is the force of infection for a given drug resistance strain, μbl is the baseline mortality rate, μh is the HIV-specific mortality rate, Ieq is an indicator function (= 1 if the condition eq is met, 0 otherwise), θ is the HIV incidence rate, τ0 is the rate of seeking diagnosis among people without TB, and sh is the specificity of the diagnostic test. Thus, uninfected individuals exit through infection and death, acquire HIV according to the HIV incidence rate, and become previously treated for TB (inappropriately) according to the specificity of the test. The HIV-uninfected, not previously treated compartment is replenished at a rate that matches total mortality and thereby maintains a constant population.

Latently infected compartment (L)

Request a detailed protocol
(2) dLhdp(t)dt=λd(t)*(1πh)*{Uhp(t)+d[Lhdp(t)*(1Ih=0*ι)]}+Ip=1*i,p[ρhdip*Phdip(t)]+Ih=0*{ν0*E0dp+i[νi*(A0dip+P0dip+I0dip)]}{[(λ0(t)+λ1(t)+λ2(t))*(1Ih=0*ι)+εh+μbl+Ih=1*μh]*Lhdp(t)}Ih=0*[θ*L0dp(t)]+Ih=1*[θ*L0dp(t)]

where λd is the strain-specific force of infection, πh is the proportion of recent infections that progress rapidly to active TB, ι is the relative reduction in rapid progression after infection among people with latent TB, Ieq is an indicator function (= 1 if the condition eq is met, 0 otherwise), ρhdip is the rate of treatment after successful diagnosis is initiated, νi is the spontaneous recovery rate, εh is the endogenous reactivation rate, μbl is the baseline mortality rate, μh is the HIV-specific mortality rate, and θ is the HIV incidence rate. Thus, individuals enter the latently infected compartment through initial TB infection (without rapid progression), reinfection (without rapid progression, and accounting for immune protection), initiation of successful treatment, or spontaneous resolution. Individuals completing treatment only enter the previously treated compartment (p=1). Individuals exit this compartment through TB reinfection, endogenous reactivation, and death, and they acquire HIV infection at a constant rate.

Early active compartment (E)

Request a detailed protocol
(3) dEhdp(t)dt=λd(t)*πh*{Uhp(t)+d[Lhdp(t)*(1Ih=0*ι)]}+Lhdp(t)*εh{[μbl+Ih=0*(μt0+ν0+δe0)+Ih=1*(μh+μth+δe1)]*Ehdp(t)}Ih=0*[θ*E0dp(t)]+Ih=1*[θ*E0dp(t)]

where λd is the force of infection, πh is the proportion of recent infections that progress rapidly to active TB, ι is the relative reduction in rapid progression after infection among people with latent TB, Ih=0 is an indicator function of HIV status (= 1 if h = 0, 0 otherwise), εh is the endogenous reactivation rate, μbl is the baseline mortality rate, μt0 is the TB-specific mortality rate for less-infectious TB, ν0 is the spontaneous recovery rate for less-infectious TB, δeh is the HIV-specific rate of progression to late active TB, μh is the HIV-specific mortality rate, μth is the TB-specific mortality rate for people living with HIV, and θ is the HIV incidence rate. Thus, individuals enter this compartment through rapid progression of recent infection or endogenous reactivation of latent infection. They exit through progression to late active disease, spontaneous cure (if HIV-uninfected), or death, and they acquire HIV infection at a constant rate.

Late active compartment (A)

Request a detailed protocol
(4) dAhdip(t)dt=δeh*[Ii=1*ψh+Ii=0*(1ψh)]*Ehdp(t)+δf*Ihdip(t)Ahdip(t)*[μbl+Ih=0*(μti+νi+σ0dip+κ0dip+Id=0*αsi0ip+Id=0*αsm0ip+Id=1*αim0ip)+Ih=1*(μh+μth+σ1dip+κ1dip+Id=0*αsi1ip+Id=0*αsm1ip+Id=1*αim1ip)]Ih=0*[θ*A0dip(t)]+Ih=1*[θ*A0dip(t)]

where δeh is the rate of progression from early active TB, Ieq is an indicator function (= 1 if the condition eq is met, 0 otherwise), ψh is the proportion of active TB that is highly infectious (smear-positive), δf is the rate (1/duration) of ineffective therapy, μbl is the baseline mortality rate, μti is the TB-specific mortality rate for non-HIV-associated TB, μh is the HIV-specific mortality rate, μth is the TB-specific mortality rate for people living with HIV, νi is the spontaneous recovery rate, σhdip is the rate of diagnosis ultimately leading to successful treatment, κhdip is the rate of placing individuals on ineffective treatment that does not generate acquired resistance, αsihip is the rate of placing individuals on ineffective treatment that generates INH monoresistance, αsmhip and αimhip are the rates of placing individuals on ineffective treatment that generates MDR-TB, and θ is the HIV incidence rate. Thus, individuals enter this compartment through progression from early active disease or relapse/failure after ineffective treatment. They exit through spontaneous recovery, diagnosis leading to successful treatment, initiation of ineffective treatment (which can, in turn, generate acquired drug resistance), or death, and they acquire HIV infection at a constant rate.

Active ‘pre-treatment’ compartment (P)

Request a detailed protocol

This compartment consists of individuals who have initiated a diagnostic attempt that will lead to successful treatment, yet remain infectious while awaiting initiation of treatment. Inclusion of this compartment is designed to capture the effects of diagnostic delays.

(5) dPhdip(t)dt=σhdip*Ahdip(t)Phdip(t)*[μbl+Ih=0*(μti+νi+ρ0dip)+Ih=1*(μh+μth+ρ1dip)]Ih=0*[θ*P0dip(t)]+Ih=1*[θ*P0dip(t)]

where σhdip is the rate of initiating successful diagnostic attempts, Ih=0 is an indicator function of HIV status (= 1 if h = 0, 0 otherwise), μbl is the non-TB mortality rate, μti is the TB-specific mortality rate for non-HIV-associated TB, μh is the HIV-specific mortality rate, μth is the TB-specific mortality rate for people living with HIV, νi is the spontaneous recovery rate, ρhdip is the rate of starting therapy after initiating a successful diagnostic attempt (1/diagnostic delay), and θ is the HIV incidence rate. Thus, individuals enter this compartment by initiation of successful diagnostic attempts and exit through initiation of treatment, spontaneous cure, or death. They acquire HIV infection at a constant rate.

Active ‘inappropriately treated’ compartment (I)

Request a detailed protocol

This compartment contains all individuals being treated whose treatment course ends in default or failure. Unlike the previous (successful treatment) compartment, diagnostic delay is not explicitly incorporated into this compartment, as to do so would prolong the duration of time until individuals who default re-enter the ‘late active’ compartment. Inclusion of a diagnostic delay before this compartment does not materially affect results. Individuals exit this compartment either in default after partial therapy—which has a defined probability of achieving cure despite not being completed—or failure. Failure leads immediately to another course of treatment, which can either be successful (i.e., transition to the latent compartment) or unsuccessful (i.e., remain in the inappropriate treatment compartment—which results in transition to the ‘previously treated’ compartment, p=1, not shown below).

(6) dIhdip(t)dt=κhdip*Ahdip(t)+Id=1*αsihdip*Ah0ip(t)+Id=2*[αsmhip*Ah0ip(t)+αimhip*Ah1ip(t)]Ihdip(t)*[μbl+Ih=0*(μt0+νi+δf)+Ih=1*(μh+μth+δf)]Ih=0*[θ*I0dip(t)]+Ih=1*[θ*I0dip(t)]

where κhdip is the rate of placing individuals on ineffective treatment that does not generate acquired resistance, αsihip is the rate of placing individuals on ineffective treatment that generates INH monoresistance from drug-susceptible TB, αsmhip is the rate of placing individuals on ineffective treatment that generates MDR-TB from drug-susceptible TB, αimhip is the rate of placing individuals on ineffective treatment that generates MDR-TB from INH-monoresistant TB, μbl is the non-TB mortality rate, μti is the TB-specific mortality rate for non-HIV-associated TB, μh is the HIV-specific mortality rate, μth is the TB-specific mortality rate for people living with HIV, νi is the spontaneous recovery rate, δf is the rate (1/duration) of ineffective therapy, and θ is the HIV incidence rate. Thus, individuals enter this compartment through ineffective treatment from the late active compartment (conditional on whether that treatment also generates new drug resistance) and exit through completion of a course of ineffective therapy, spontaneous resolution, or death. They acquire HIV infection at a constant rate.

Model fitting

Request a detailed protocol

The equations above comprise a series of 100 ordinary differential equations. In order to generate an equilibrium population according to user specifications of TB incidence, MDR-TB prevalence, and HIV prevalence, it is necessary to solve for the roots of a system of these 100 equations, plus three additional equations to account for the user inputs. We accomplish this using the ‘fsolve’ routine in SciPy (www.scipy.org). To solve for these three additional equations, we first match each user input to a single parameter: TB incidence to the transmission rate β, MDR-TB prevalence to the relative reduction in transmission for MDR-TB φ2, and HIV prevalence to the HIV incidence rate θ. We then constrain the system of equations such that the total population remains constant at 100,000, and we add the following three equations to the system:

dβ/dt=(calculated TB incidence)(user-defined TB incidence)
dφ2/dt=(calculated MDR-TB prevalence)(user-defined MDR-TB prevalence)
dθ/dt=(calculated HIV prevalence)(user-defined HIV prevalence)

By solving for the roots at which this system of equations equals zero, we generate an equilibrium (steady-state) population that also defines β, φ2, and θ such that the user-defined targets are also met.

Additional complexity is added to account for non-equilibrium in TB incidence and MDR-TB prevalence. We accomplish this by fitting an equilibrium defined by the user-specified target of HIV prevalence, and the user-specified ‘initial’ values of TB incidence and MDR-TB prevalence. This equilibrium is set to be 50 years in the past; at this time, the system of equations is solved as above. We then allow for the parameters β and φ2 to be altered from their original values (corresponding to the equilibrium condition) such that a second set of targets are attained after 50 years (start of the analysis). If these calculated values at the end of 50 years do not match the user-defined targets, the parameter values are changed accordingly, and the model is re-run from equilibrium until the appropriate parameter value is identified that generates the user-defined TB incidence and MDR-TB prevalence, within a relative tolerance of 0.05 in each variable. These parameters may then be further modified such that they change over the analysis frame of five years (e.g., to describe an epidemic with increasing MDR-TB prevalence over time).

In the primary version of the model, this is only automated for low-incidence scenarios in which the user inputs a TB incidence of less than 50 per 100,000/year. In such cases, the model generates an equilibrium population 50 years in the past with a TB incidence of 50 per 100,000/year and reduces the transmission parameter β until the user-defined TB incidence is achieved 50 years later (i.e., the start of the analytic period). This accounts for the fact that most low-incidence settings have substantially more latent TB infection than would be estimated by an equilibrium population with a very low TB incidence—thereby more appropriately reflecting a higher proportion of TB due to reactivation rather than recent infection. However, we include code (described in the user manual below) that allows users to define high incidence scenarios that are likewise not at equilibrium, as well as ‘emerging MDR’ scenarios in which the prevalence of MDR-TB is increasing rather than stable.

References

    1. Brownell R
    2. Metcalfe J
    3. Millman AJ
    4. Miller C
    5. Cattamanchi AP
    (2012)
    Performance of Xpert MTB/RIF for diagnosis of pulmonary and extrapulmonary tuberculosis - a systematic review and meta-anaylsis
    American Journal of Respiratory and Critical Care Medicine 185:A4715.
    1. McKay MD
    2. Beckman RJ
    3. Conover WJ
    (1979)
    Comparison of three methods for selecting values of input variables in the analysis of output from a computer code
    Technometrics 21:239–245.
  1. Book
    1. UNAIDS
    (2012)
    UNAIDS report on the global AIDS epidemic 2012
    Geneva: UNAIDS.
    1. Wilkinson D
    2. Newman W
    3. Reid A
    4. Squire SB
    5. Sturm AW
    6. Gilks CF
    (2000)
    Trial-of-antibiotic algorithm for the diagnosis of tuberculosis in a district hospital in a developing country with high HIV prevalence
    The International Journal of Tuberculosis and Lung Disease 4:513–518.
    1. World Bank
    (2012) Life expectancy
    http://www.google.com/publicdata/explore?ds=d5bncppjof8f9_&ctype=l&strail=false&nselm=h&hl=en&dl=en#!ctype=l&strail=false&bcs=d&nselm=h&met_y=sp_dyn_le00_in&scale_y=lin&ind_y=false&rdim=region&ifdim=region&tdim=true&hl=en_US&dl=en&ind=false, Accessed July 11.
  2. Book
    1. World Health Organization
    (2010)
    Multidrug and extensively drug-resistant TB (M/XDR-TB): 2010 global report on surveillance and response
    Geneva: WHO.
  3. Book
    1. World Health Organization
    (2011)
    Same-day diagnosis of tuberculosis by microscopy: policy statement
    Geneva: WHO.
  4. Book
    1. World Health Organization
    (2012)
    Global tuberculosis report 2012
    Geneva: WHO.

Article and author information

Author details

  1. David W Dowdy

    1. Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Baltimore, United States
    2. Center for Tuberculosis Research, Johns Hopkins University, Baltimore, United States
    Contribution
    DWD, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    For correspondence
    ddowdy1@jhmi.edu
    Competing interests
    The authors declare that no competing interests exist.
  2. Jason R Andrews

    Division of Infectious Diseases, Department of Medicine, Massachusetts General Hospital, Boston, United States
    Present address
    Division of Infectious Diseases, Department of Medicine, Stanford University School of Medicine, Palo Alto, United States
    Contribution
    JRA, Conception and design, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  3. Peter J Dodd

    School of Health and Related Research, University of Sheffield, Sheffield, United Kingdom
    Contribution
    PJD, Conception and design, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  4. Robert H Gilman

    Department of International Health, Johns Hopkins Bloomberg School of Public Health, Baltimore, United States
    Contribution
    RHG, Conception and design, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.

Funding

National Institutes of Health (NIH) (R21AI101152)

  • David W Dowdy

Canadian Institutes of Health Research (MOP 271997)

  • David W Dowdy

The funder had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We are grateful to Jeff Pennington for development of the web interface, and to Adithya Cattamanchi, Madhukar Pai, Jonathan Golub, and Richard Chaisson for helpful comments.

Copyright

© 2014, Dowdy et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,897
    views
  • 160
    downloads
  • 13
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. David W Dowdy
  2. Jason R Andrews
  3. Peter J Dodd
  4. Robert H Gilman
(2014)
A user-friendly, open-source tool to project impact and cost of diagnostic tests for tuberculosis
eLife 3:e02565.
https://doi.org/10.7554/eLife.02565

Share this article

https://doi.org/10.7554/eLife.02565

Further reading

    1. Epidemiology and Global Health
    Xiaoning Wang, Jinxiang Zhao ... Dong Liu
    Research Article

    Artificially sweetened beverages containing noncaloric monosaccharides were suggested as healthier alternatives to sugar-sweetened beverages. Nevertheless, the potential detrimental effects of these noncaloric monosaccharides on blood vessel function remain inadequately understood. We have established a zebrafish model that exhibits significant excessive angiogenesis induced by high glucose, resembling the hyperangiogenic characteristics observed in proliferative diabetic retinopathy (PDR). Utilizing this model, we observed that glucose and noncaloric monosaccharides could induce excessive formation of blood vessels, especially intersegmental vessels (ISVs). The excessively branched vessels were observed to be formed by ectopic activation of quiescent endothelial cells (ECs) into tip cells. Single-cell transcriptomic sequencing analysis of the ECs in the embryos exposed to high glucose revealed an augmented ratio of capillary ECs, proliferating ECs, and a series of upregulated proangiogenic genes. Further analysis and experiments validated that reduced foxo1a mediated the excessive angiogenesis induced by monosaccharides via upregulating the expression of marcksl1a. This study has provided new evidence showing the negative effects of noncaloric monosaccharides on the vascular system and the underlying mechanisms.

    1. Epidemiology and Global Health
    2. Microbiology and Infectious Disease
    Amanda C Perofsky, John Huddleston ... Cécile Viboud
    Research Article

    Influenza viruses continually evolve new antigenic variants, through mutations in epitopes of their major surface proteins, hemagglutinin (HA) and neuraminidase (NA). Antigenic drift potentiates the reinfection of previously infected individuals, but the contribution of this process to variability in annual epidemics is not well understood. Here, we link influenza A(H3N2) virus evolution to regional epidemic dynamics in the United States during 1997—2019. We integrate phenotypic measures of HA antigenic drift and sequence-based measures of HA and NA fitness to infer antigenic and genetic distances between viruses circulating in successive seasons. We estimate the magnitude, severity, timing, transmission rate, age-specific patterns, and subtype dominance of each regional outbreak and find that genetic distance based on broad sets of epitope sites is the strongest evolutionary predictor of A(H3N2) virus epidemiology. Increased HA and NA epitope distance between seasons correlates with larger, more intense epidemics, higher transmission, greater A(H3N2) subtype dominance, and a greater proportion of cases in adults relative to children, consistent with increased population susceptibility. Based on random forest models, A(H1N1) incidence impacts A(H3N2) epidemics to a greater extent than viral evolution, suggesting that subtype interference is a major driver of influenza A virus infection ynamics, presumably via heterosubtypic cross-immunity.