Abstract Background The existing estimate of the global burden of latent TB infection (LTBI) as “one-third” of the world population is nearly 20 y old. Given the importance of controlling LTBI as part of the End TB Strategy for eliminating TB by 2050, changes in demography and scientific understanding, and progress in TB control, it is important to re-assess the global burden of LTBI. Methods and Findings We constructed trends in annual risk in infection (ARI) for countries between 1934 and 2014 using a combination of direct estimates of ARI from LTBI surveys (131 surveys from 1950 to 2011) and indirect estimates of ARI calculated from World Health Organisation (WHO) estimates of smear positive TB prevalence from 1990 to 2014. Gaussian process regression was used to generate ARIs for country-years without data and to represent uncertainty. Estimated ARI time-series were applied to the demography in each country to calculate the number and proportions of individuals infected, recently infected (infected within 2 y), and recently infected with isoniazid (INH)-resistant strains. Resulting estimates were aggregated by WHO region. We estimated the contribution of existing infections to TB incidence in 2035 and 2050. In 2014, the global burden of LTBI was 23.0% (95% uncertainty interval [UI]: 20.4%–26.4%), amounting to approximately 1.7 billion people. WHO South-East Asia, Western-Pacific, and Africa regions had the highest prevalence and accounted for around 80% of those with LTBI. Prevalence of recent infection was 0.8% (95% UI: 0.7%–0.9%) of the global population, amounting to 55.5 (95% UI: 48.2–63.8) million individuals currently at high risk of TB disease, of which 10.9% (95% UI:10.2%–11.8%) was isoniazid-resistant. Current LTBI alone, assuming no additional infections from 2015 onwards, would be expected to generate TB incidences in the region of 16.5 per 100,000 per year in 2035 and 8.3 per 100,000 per year in 2050. Limitations included the quantity and methodological heterogeneity of direct ARI data, and limited evidence to inform on potential clearance of LTBI. Conclusions We estimate that approximately 1.7 billion individuals were latently infected with Mycobacterium tuberculosis (M.tb) globally in 2014, just under a quarter of the global population. Investment in new tools to improve diagnosis and treatment of those with LTBI at risk of progressing to disease is urgently needed to address this latent reservoir if the 2050 target of eliminating TB is to be reached.

Author Summary Why Was This Study Done? Addressing the latent TB infection reservoir is critical to achieving TB elimination.

The current estimate that “one-third” of the global population is infected with tuberculosis is widely cited but has not been formally estimated for nearly 20 y.

Changes in demography, the size and distribution of TB burden, as well as new scientific insights and the availability of new data mean a re-estimation is needed. What Did the Researchers Do and Find? We generated an annual risk of infection between 1934 and 2014 and applied this to a country-level demographic model, quantifying uncertainty wherever possible.

We estimated that approximately 1.7 billion individuals were infected with LTBI in 2014; just under a quarter of the global population.

If left unaddressed, the current LTBI burden alone will likely prevent achieving the global TB targets for TB elimination. What Do These Findings Mean? For long-term TB control to be successful, an aggressive approach to LTBI is needed.

Research and development should focus on developing better tools to identify individuals who will benefit from LTBI treatment.

Estimates would be strengthened by additional empirical data from new population-based studies of LTBI prevalence.

Citation: Houben RMGJ, Dodd PJ (2016) The Global Burden of Latent Tuberculosis Infection: A Re-estimation Using Mathematical Modelling. PLoS Med 13(10): e1002152. https://doi.org/10.1371/journal.pmed.1002152 Academic Editor: John Z. Metcalfe, University of California San Francisco, UNITED STATES Received: April 5, 2016; Accepted: September 9, 2016; Published: October 25, 2016 Copyright: © 2016 Houben, Dodd. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability: All relevant data are within the paper and its Supporting Information files. Funding: RMGJH was partly funded by the Bill and Melinda Gates Foundation (TB Modelling and Analysis Consortium: OPP1084276). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing interests: The authors have declared that no competing interests exist. Abbreviations: ARI, annual risk in infection; INH, isoniazid; LTBI, latent TB infection; M.tb, Mycobacterium tuberculosis; PLHIV, people living with HIV; SDG, Sustainable Development Goals; TST, tuberculin skin test; WHO, World Health Organisation

Introduction Infection with Mycobacterium tuberculosis (M.tb) is the precursor to TB disease, which is responsible for 1.5 million deaths each year—more than any other infectious disease [1]. Once infected, the individual is at highest risk of developing TB disease within the first two years, but can remain at risk for their lifetime [2]. The population carrying a latent TB infection (LTBI) is commonly quoted as “one-third” of the global population, a reservoir of approximately 2.3 billion individuals [3–6]. As the global community looks to meet ambitious targets for reduction (90% reduction in TB incidence by 2035) and even elimination of TB (less than 1 incident case per 1,000,000 per year) by 2050 [7], our ability to address the LTBI reservoir will be critical in our chance to succeed. Despite its clear importance to global TB control efforts, the most recent attempt to estimate the global burden of LTBI was in 1998 [3]. Since then, the size and distribution of the global population [8] and TB burden [1] has changed dramatically, as has our understanding of prevalent disease as a driver of infection [9,10]. Global population growth from around 6 billion in 1998 to over 7 billion in 2014 has been mainly driven by areas with the highest TB burden, such as Southeast Asia and sub-Saharan Africa [1,8]. The previous estimation method relied on a fixed relationship between TB burden to the annual risk of acquiring LTBI, the so-called “Styblo rule” [3]. Since then, two groups have shown this long-held rule of thumb substantially overestimates infection risk in modern populations [9,10]. Given these changes and the drive towards eliminating TB [7,11], an updated estimate of the global burden of LTBI that incorporates the available data and applies current scientific insights is urgently needed [4,12]. An updated estimate of the size and distribution of the LTBI reservoir should also address questions about the likely contribution of LTBI to TB disease over the coming decades. Specifically, how many active TB cases would arise from the currently infected individuals alone if all transmission was halted now? Updated LTBI burden estimates also indicate the population in need of interventions and new tools, thus catalyzing new research and potential investment from commercial partners in, for example, vaccines, and tools for the diagnosis and treatment of LTBI [4]. Critical questions include the number of those with LTBI at highest risk of developing disease, i.e., those infected within the past 2 y [12]. This population is a focus of proposed “test and treat” strategies in TB, which would use an RNA expression profile test to identify individuals most likely to develop TB [13], improving on the low predictive value of existing tests for LTBI [14]. As resistance to TB drugs is rising, an estimate of the proportion of LTBI that involves isoniazid (INH)-resistant strains is important, since INH remains the cornerstone of most treatment regimens for LTBI [4]. Finally, as TB becomes rarer, the epidemiology of LTBI will have a renewed and increasing importance for monitoring the progress of control efforts. In this paper, we estimate the global burden of LTBI and its distribution by country, geographical region, and age group. We also estimate the number of recent infections and the number of recent infections with INH-resistant strains. Finally, we predict the TB incidence in 2035 and 2050 solely due to the existing LTBI reservoir.

Methods To estimate the burden of LTBI, we reconstructed country trends in annual risk of M.tb infection (ARI) and combined these historical projections with demographic data to estimate burden of recent and all-time infection by age. ARI trends were modelled for 168 countries (comprising >99.9% of the world population) using a flexible non-parametric regression framework accounting for measurement uncertainty and applied to two sources of data. The first source of data was direct estimates of ARI from tuberculin skin test (TST) surveys. We abstracted data on estimated ARI, survey sample size, and mean age for 100 country-years in 24 countries from Cauthen et al. [15] and undertook a systematic review of nationally representative ARI estimates in the years 1990–2014 (see Methods and Figure A in S1 Text), yielding further data on 31 country-years in 19 countries, to give a total of 37 unique countries with TST survey data. Historically, ARIs have been estimated from TST surveys without presentation of uncertainty. In the Methods section of S1 Text, we show how sample size and mean age can be used to approximately and conservatively quantify measurement precision. We used this method for studies not stating ARI estimate precision. The second source of data was indirect: combining WHO estimates of TB prevalence (5,373 country-years for 218 countries) with an uncertain representation of the revised Styblo ratio that accounts for uncertainty and relates the prevalence of smear positive disease to ARI [1,9,10]. A previously published study of childhood tuberculosis [16] characterized this ratio in the modern era by fitting a log-normal distribution to data from reviews of studies in which both ARI and prevalence estimates were available [9,10]. To estimate the proportion of prevalent TB that is smear positive for each country, we averaged estimates of smear positivity for 0–4, 5–14, and ≥15 y age groups from a recent systematic review [17] against the proportion of cases in these age-groups calculated using the model of Dodd et al. [16]. To calculate the impact of HIV on the proportion of TB in a country that is smear positive, we first calculated the fraction of prevalent TB in people living with HIV (PLHIV) by adjusting the WHO estimates of HIV prevalence in incident TB for each year using estimates of the duration of prevalence in PLHIV and HIV-uninfected individuals [1]. The mean smear positivity of HIV-TB was then reduced by a fraction reported in Corbett et al. [18]. The uncertainty of each ingredient in these ARI calculations was propagated using the delta method. More details are described in the Methods section of S1 Text. Gaussian process regression with a linear trend was applied to the data on ARI (on a log scale), using the measurement precision calculated for each data point. This implicitly combines the data feeding into WHO prevalence estimates (and Styblo ratio) and the TST data feeding into ARI estimates, with the assumption of a normal approximation to the likelihood. A sensitivity analysis re-analysed these data with a constant rather than linear trend assumption. To allow a comparison with the 1998 estimate, we assumed a constant ARI before 1934. For each country, 200 simulated ARI trajectories from 1934 to 2014 were used to compute the cumulative hazard of infection for individuals by age. The cumulative hazard was converted into a probability of infection and combined with UN Population Division estimates of country demography in 2014 to give our estimates of all-time infection. We computed the probability of infection for the first time within 2 y, using the cumulative hazard up to 2012 to calculate the fraction at each age who had escaped infection until then, and the cumulative hazard from 2012 to 2014 to calculate the fraction of these who were then infected. To calculate the fraction of those at a given age infected or re-infected during the last 2 y, we introduced a beta distribution characterizing an uncertain partial protection against re-infection of 79% (70%–86%) from Andrews et al. [19]. See the Methods section of S1 Text for details. As estimates of protection have varied [20,21], we conducted a sensitivity analysis using a protection of 50% with the same variance. Estimates of infection prevalence were summarized by region, age, and medians mapped by country. To calculate the number of infections within the last 2 y with resistance to isoniazid, we combined our estimates of infection within 2 y in each country with a recent analysis of the proportion of new infections in each country that are isoniazid-resistance using data from the Global Project on Anti-tuberculosis Drug Resistance Surveillance at WHO [22]. This proportion was treated as uncertain and sampled from the output of this analysis. A conceptual overview diagram for the methods is presented in Figure R in S1 Text. Finally, we estimated the regional prevalence of latent infection in 2035 and 2050 under the assumption of no M.tb transmission after 2014, using UN Population Division demographic projections, and calculated the likely implications of existing M.tb infections for future TB incidence, assuming a 0.15% per year remote activation rate [21,23,24]. Results are reported as medians together with 95% uncertainty intervals (95% UI), calculated as the 2.5% to 97.5% percentile range. No specific funding was received for this work.

Discussion The global burden of LTBI is just under a quarter of the population—around 1.7 billion individuals—with substantial geographical and age variation. We estimate 56 million people are at high risk of developing TB disease because of a recent (re-)infection, 11% of whom are carrying an isoniazid-resistant strain. With reasonable assumptions for reactivation risks, incident TB disease arising from the 2014 LTBI reservoir alone would prohibit reaching the 2035 and 2050 End TB Strategy goals. Using our method to calculate LTBI prevalence in 1997 yielded 27%, suggesting the difference between our estimate for 2014 and the 1997 estimate [3] is mainly due to changes in methodology, including the revised Styblo rule [9,10]. Our results also matched a recent survey-based estimate of 13 million latent infections for the United States [25]. One limitation of this study was that we were not able to consider the effects due to different tuberculin strains and cut-off choices for TST tests or the potential impact of BCG vaccination on these test outcomes. Recent work has suggested that LTBI may be a dynamic and heterogeneous state [12] and highlighted the difference between M.tb infection and positive tests for infection. Our work does not make this distinction but is in line with literature relating TB prevalence to infection risk and risks of progression to disease. In addition, we assumed lifelong LTBI in common with previous estimates and consistent with observation [2]. While it may be biologically plausible that some individuals clear their LTBI in absence of treatment, there is no published evidence to support a separate model scenario. Our approach to quantifying the typical infectiousness of prevalent TB cases by country entailed some necessary simplifications. Notably, we assumed the same TB case-detection rate applied independent of HIV status and anti-retroviral treatment status. However, these assumptions had relatively little influence on ARI estimates and their uncertainty (see S1 Text, page 8), with the dominant contributions to uncertainty coming from the TB prevalence estimates themselves and uncertainty in the Styblo ratio. We also assumed that the same ARI applied to all individuals in a country and neglected migration between countries. A major strength of this study is its treatment of uncertainty. We were able to characterize and include measurement precision for data on ARI derived from TST surveys as well as indirectly from prevalence estimates. While uncertainty in ARI estimates grew substantially for the earliest years, the number of people alive today who were alive then was small, limiting the impact of this uncertainty on our overall estimate. This also limited the impact of our assumption of linear trends extrapolating backward from data; our very conservative assumption of flat trends resulted in less than a five-percentage-point difference in our overall prevalence estimate. Our estimate of TB incidence in 2035 and 2050 from currently existing LTBI uses a single reactivation rate. This parameter has been estimated at widely varying levels, [23,24,26], and we have not attempted to include potential geographic heterogeneity or anticipate trends due to changes in population health, for example, through achievements of the Sustainable Development Goals (SDG) [27] or due to cohort aging. However, our estimates of TB disease incidence in 2035 and 2050 are proportional to this parameter, and even large reductions do not alter the conclusion that incidence will exceed the 2050 elimination target. It is a limitation that there were not more recent direct national estimates of ARI, which would have strengthened our estimates. More empirical data on the epidemiology of LTBI, including from the use of modern tests that are less prone to biases of interpretation and cross-reaction [28], should improve our understanding of these features and generate more data upon which to base estimates. As active TB becomes rarer, surveys of infection will become an increasingly appealing option for monitoring trends, but this relies on understanding the relationship between infection and other burden metrics. Policies to address the LTBI reservoir have to balance the potential of harm versus the benefit for the individual [29]. In the current landscape of diagnostic tools, WHO recommends LTBI testing and treatment only in high-risk groups, such as people living with HIV, and close contacts of TB cases. While these guidelines are sensible, it is clear that a more aggressive approach is needed to reduce the threat to long-term TB control targets stemming from a LTBI reservoir of approximately 1.7 billion individuals. Future work with this model could inform current policies by estimating the burden of LTBI in specific risk groups, such as people living with HIV or diabetes. A test to more precisely identify those at substantial risk of progressing to disease could enable targeted LTBI treatment beyond known risk groups [12]. Emerging tests based on RNA signatures may come to provide a more practicable method of identifying individuals for LTBI treatment [13]. Among biomedical interventions, a vaccine that prevents progression to disease from LTBI could make a major contribution, depending on global availability [30]. Beyond the biomedical perspective, improvements in social and economic conditions globally have been associated with reductions in TB burden in historic and contemporary contexts [31,32], and could also contribute to reducing the TB burden originating from the LTBI reservoir. Treatment for LTBI still relies heavily on isoniazid, either as monotherapy or as part of a combination regimen [4,11]. We found that just under 11% of all recent M.tb infections are likely to be isoniazid resistant, with much higher rates in some regions, and this proportion is likely to increase. While less common, rifampicin resistance also has the potential to threaten the usefulness of rifampicin-containing prophylactic regimens. New treatments that bypass the rising resistance to isoniazid and rifampicin are needed to fully operationalise interventions to test and treat LTBI.

Conclusion We estimate that approximately 1.7 billion individuals were latently infected with M.tb globally in 2014, just under a quarter of the global population. Investment in new tools to improve diagnosis and treatment of those with LTBI at risk of progressing to disease is urgently needed to address this latent reservoir if the TB community is to reach the 2050 target of eliminating TB.

Acknowledgments We thank Richard White and other members from the London School of Hygiene and Tropical Medicine TB modelling group for their valuable comments on the initial design and results from this work.

Author Contributions Conceptualization: RMGJH PJD. Investigation: RMGJH PJD. Methodology: RMGJH PJD. Writing – original draft: RMGJH PJD. Writing – review & editing: RMGJH PJD.