A user-friendly, open-source tool to project impact and cost of diagnostic tests for tuberculosis
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.001eLife 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.002Introduction
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:
‘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.)
‘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).
‘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.
‘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.
‘Xpert for All’: Xpert MTB/RIF for all patients.
‘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.
‘MODS/TLA’: Sputum smear, plus microcolony-based TB culture (e.g., MODS or thin-layer agar, TLA) for all patients.
‘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).
‘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]).
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:
‘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]);
‘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;
‘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
‘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).
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.
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.
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 protocolUsing 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 protocolAfter 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 protocolThe 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 protocolThe 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 protocolAs 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 protocolWe 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 protocolA 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 protocolEconomic 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 protocolUnder 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.
Outcomes and sensitivity analysis
Request a detailed protocolOur 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 protocolThus, 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 protocolThus, 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 protocolwhere μ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 protocolwhere λ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 protocolwhere λ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 protocolwhere δ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 protocolThis 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.
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 protocolThis 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).
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 protocolThe 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:
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
-
Epidemiological benefits of more-effective tuberculosis vaccines, drugs, and diagnosticsProceedings of the National Academy of Sciences of the United States of America 106:13980–13985.https://doi.org/10.1073/pnas.0901720106
-
Risk of progression to active tuberculosis following reinfection with Mycobacterium tuberculosisClinical Infectious Diseases 54:784–791.https://doi.org/10.1093/cid/cir951
-
Risk factors for tuberculosis in HIV-infected persons. A prospective cohort study. The Gruppo Italiano di Studio Tubercolosi e AIDS (GISTA)The Journal of the American Medical Association 274:143–148.https://doi.org/10.1001/jama.1995.03530020061033
-
The global diabetes model: user friendly version 3.0Diabetes Research and Clinical Practice 50:S15–S46.https://doi.org/10.1016/S0168-8227(00)00215-1
-
Rapid molecular detection of tuberculosis and rifampin resistanceThe New England Journal of Medicine 363:1005–1015.https://doi.org/10.1056/NEJMoa0907847
-
Modelling HIV epidemics in the antiretroviral era: the UNAIDS Estimation and Projection package 2009Sexually Transmitted Infections 86:ii3–ii10.https://doi.org/10.1136/sti.2010.044784
-
Performance of Xpert MTB/RIF for diagnosis of pulmonary and extrapulmonary tuberculosis - a systematic review and meta-anaylsisAmerican Journal of Respiratory and Critical Care Medicine 185:A4715.
-
Which new diagnostics for tuberculosis, and when?The Journal of Infectious Diseases 205:S191–S198.https://doi.org/10.1093/infdis/jis188
-
Human immunodeficiency virus and the prevalence of undiagnosed tuberculosis in African gold minersAmerican Journal of Respiratory and Critical Care Medicine 170:673–679.https://doi.org/10.1164/rccm.200405-590OC
-
The growing burden of tuberculosis: global trends and interactions with the HIV epidemicArchives of Internal Medicine 163:1009–1021.https://doi.org/10.1001/archinte.163.9.1009
-
Meta-analysis of BACTEC MGIT 960 and BACTEC 460 TB, with or without solid media, for detection of mycobacteriaJournal of Clinical Microbiology 42:2321–2325.https://doi.org/10.1128/JCM.42.5.2321-2325.2004
-
Is passive diagnosis enough? the impact of subclinical disease on diagnostic strategies for tuberculosisAmerican Journal of Respiratory and Critical Care Medicine 187:543–551.https://doi.org/10.1164/rccm.201207-1217OC
-
Standard short-course chemotherapy for drug-resistant tuberculosis: treatment outcomes in 6 countriesThe Journal of the American Medical Association 283:2537–2545.https://doi.org/10.1001/jama.283.19.2537
-
Lifetime costs and quality-adjusted life years saved from HIV prevention in the test and treat eraJournal of Acquired Immune Deficiency Syndromes 64:e15–e18.https://doi.org/10.1097/QAI.0b013e3182a5c8d4
-
Revisiting rates of reactivation tuberculosis: a population-based approachAmerican Journal of Respiratory and Critical Care Medicine 182:420–425.https://doi.org/10.1164/rccm.200909-1355OC
-
Location of Xpert(R) MTB/RIF in centralised laboratories in South Africa undermines potential impactThe International Journal of Tuberculosis and Lung Disease 16:701.https://doi.org/10.5588/ijtld.12.0131
-
Microcolony culture techniques for tuberculosis diagnosis: a systematic reviewThe International Journal of Tuberculosis and Lung Disease 16:16–23.https://doi.org/10.5588/ijtld.10.0065
-
Initial drug resistance and tuberculosis treatment outcomes: systematic review and meta-analysisAnnals of Internal Medicine 149:123–134.https://doi.org/10.7326/0003-4819-149-2-200807150-00008
-
Pre-treatment loss to follow-up in low- and lower-middle-income countries and high-burden countries: a systematic review and meta-analysisBulletin of the World Health Organization 92:126–138.https://doi.org/10.2471/BLT.13.124800
-
Comparison of three methods for selecting values of input variables in the analysis of output from a computer codeTechnometrics 21:239–245.
-
Microscopic-observation drug-susceptibility assay for the diagnosis of TBThe New England Journal of Medicine 355:1539–1550.https://doi.org/10.1056/NEJMoa055524
-
Cost analysis of rapid methods for diagnosis of multidrug resistant tuberculosis in different epidemiologic groups in PeruRevista Peruana de Medicina Experimental y Salud Publica 28:426–431.https://doi.org/10.1590/S1726-46342011000300004
-
Fluorescence versus conventional sputum smear microscopy for tuberculosis: a systematic reviewThe Lancet Infectious Diseases 6:570–581.https://doi.org/10.1016/S1473-3099(06)70578-3
-
Sputum processing methods to improve the sensitivity of smear microscopy for tuberculosis: a systematic reviewThe Lancet Infectious Diseases 6:664–674.https://doi.org/10.1016/S1473-3099(06)70602-8
-
Projecting the demographic consequences of adult HIV prevalence trends: the Spectrum Projection PackageSexually Transmitted Infections 80:i14–i18.https://doi.org/10.1136/sti.2004.010157
-
Spectrum: a model platform for linking maternal and child survival interventions with AIDS, family planning and demographic projectionsInternational Journal of Epidemiology 39:i7–i10.https://doi.org/10.1093/ije/dyq016
-
Overview of the Lives Saved Tool (LiST)BMC Public Health 13:S1.https://doi.org/10.1186/1471-2458-13-S3-S1
-
Trial-of-antibiotic algorithm for the diagnosis of tuberculosis in a district hospital in a developing country with high HIV prevalenceThe International Journal of Tuberculosis and Lung Disease 4:513–518.
-
Antiretroviral therapy for tuberculosis control in nine African countriesProceedings of the National Academy of Sciences of the United States of America 107:19485–19489.https://doi.org/10.1073/pnas.1005660107
-
Undiagnosed tuberculosis in a community with high HIV prevalence: implications for tuberculosis controlAmerican Journal of Respiratory and Critical Care Medicine 175:87–93.https://doi.org/10.1164/rccm.200606-759OC
-
Life expectancyhttp://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.
-
BookMultidrug and extensively drug-resistant TB (M/XDR-TB): 2010 global report on surveillance and responseGeneva: WHO.
Article and author information
Author details
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
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)
Further reading
-
- Epidemiology and Global Health
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.
-
- Epidemiology and Global Health
- Microbiology and Infectious Disease
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.