Mathematical modelling of the impact of treating latent tuberculosis infection in the elderly in a city with intermediate tuberculosis burden

Scientific Reports, Mar 2019

Hong Kong is a high-income city with intermediate tuberculosis (TB) burden primarily driven by endogenous reactivations. A high proportion of remote latently infected people, particularly elderly, hinders the effectiveness of current strategies focusing on passive TB detection. In this study, we developed a mathematical model to evaluate the impact of treating latent TB infection (LTBI) in the elderly in addition to current TB control strategies. The model was calibrated using the annual age-stratified TB notifications from 1965–2013 in Hong Kong. Our results showed that at present, approximately 75% of annual new notifications were from reactivations. Given the present treatment completion rate, even if only a low to moderate proportion (approximately 20% to 40%) of elderly people were screened and treated for LTBI, the overall TB incidence could be reduced by almost 50%, to reach the 2025 milestone of the global End TB Strategy. Nevertheless, due to a high risk of hepatotoxicity in elderly population, benefit-risk ratios were mostly below unity; thus, intervention programs should be carefully formulated, including prioritising LTBI treatment for high-risk elderly groups who are closely monitored for possible adverse side effects.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://www.nature.com/articles/s41598-019-41256-4.pdf

Mathematical modelling of the impact of treating latent tuberculosis infection in the elderly in a city with intermediate tuberculosis burden

Mathematical modelling of the impact of treating latent tuberculosis infection in the elderly in a city with intermediate tuberculosis burden Ka Chun Chong Greta Chun Huen tam 0 Maggie Haitian Wang 0 1 Katherine Min Jia Chi Chiu Leung 0 1 2 Wing WaiYew 2 Benny ChungYing Zee Pui Hong Chung 0 steven Yuk Fai Lau 0 Xiaoran Han 0 eng Kiong Yeoh 0 0 Jc School of Public Health and Primary care, the chinese University of Hong Kong , Hong Kong, china 1 clinical t rials and Biostatistics Laboratory, Shenzhen Research institute, the chinese University of Hong Kong , Shenzhen, china 2 Stanley Ho c entre for emerging infectious Diseases, t he c hinese University of Hong Kong, Hong Kong , c hina OPEN Published: xx xx xxxx Hong Kong is a high-income city with intermediate tuberculosis (TB) burden primarily driven by endogenous reactivations. A high proportion of remote latently infected people, particularly elderly, hinders the effectiveness of current strategies focusing on passive TB detection. In this study, we developed a mathematical model to evaluate the impact of treating latent TB infection (LTBI) in the elderly in addition to current TB control strategies. The model was calibrated using the annual age-stratified TB notifications from 1965-2013 in Hong Kong. Our results showed that at present, approximately 75% of annual new notifications were from reactivations. Given the present treatment completion rate, even if only a low to moderate proportion (approximately 20% to 40%) of elderly people were screened and treated for LTBI, the overall TB incidence could be reduced by almost 50%, to reach the 2025 milestone of the global End TB Strategy. Nevertheless, due to a high risk of hepatotoxicity in elderly population, benefit-risk ratios were mostly below unity; thus, intervention programs should be carefully formulated, including prioritising LTBI treatment for high-risk elderly groups who are closely monitored for possible adverse side effects. - 1,2, Latent tuberculosis infection (LTBI) due to Mycobacterium tuberculosis may persist for many years. Certain risk factors may predispose to ?endogenous reactivation? of tuberculosis (TB); alternatively, ?exogenous reinfection? can occur1. Hong Kong is a high-income region with intermediate TB burden that has had decreasing TB notifications since the 1950s; however, this decrease has slowed since the 1970s2?4. Although cases resulting from primary infection and reinfection have decreased in Hong Kong, current measures using passive screening and Directed Observed Treatment, Short-course (DOTS) are insufficient as it only helps reducing recent transmissions but is unable to prevent active TB development from reactivation from remote infection, which has become the main source of TB incidence as with many intermediate burden areas1?3. The epidemiology of TB in Hong Kong is similar with some countries, with persistently higher incidence and risk of reactivation2,4,5. An early local study showed that the estimated proportion of recently transmitted TB was only 15% to 20%6. A study of 2243 Hong Kong residents in elderly care homes noted a high prevalence of reactivated TB7. Reactivation occurs predominantly in the elderly population1,2; thus, older-age groups account for the current TB burden, with 42.5% of notifications from people aged ?6 5 years in 2016 . Under the World Health Organization (WHO) End TB Strategy, countries should target 50% and 95% reductions in the 2015 TB incidence by 2025 and 2035, respectively. To accelerate the decrease in TB incidence, targeted LTBI diagnosis and treatment have become the focus of future interventions in low-to-moderate incidence regions dominated by reactivation cases8. In Singapore, the TB incidence rate increased after 2008 and remained at 47/100,000 in 2017, in which an estimate of LTBI prevalence of 10.7% was found9,10. In Japan, the notification rates dropped slowly in recent years and an increasing number of LTBI notification from 2007 to 2016 was reported11. Not just in Asia, a recent investigation indicated around 1.7 billion individuals were infected with LTBI in 2014 which closed to a quarter of the world population12. Although United States is a low TB burden country, the elimination has not yet been met. One of the primary reasons is that the prevalence of LTBI decreased slowly between 2000 and 2012 suggesting interventions for LTBI control is needed13. Although LTBI diagnosis and treatment are recommended, the LTBI control varied between countries. While there is no practical programs for LTBI screening and treatment in Hong Kong, active screening of elderly in community is recommended in Japan14. In the United States, screening is recommended to older people being newly admitted to long-term care facilities15. To inform future LTBI management strategies, we investigated the public health impact of an expanded LTBI screening and treatment program targeting the elderly, as the high number of reactivated cases is a barrier for TB prevention. LTBI treatment is available particularly for elderly with co-morbidities that imply a higher risk of progression to active diseases but acceptance and treatment completion rates vary, possibly due to adverse reactions associated with treatment, particularly hepatotoxicity induced by isoniazid (INH) therapy16?18. To our knowledge, the effectiveness of LTBI treatment for elderly people in Hong Kong remains unclear. The following issues need to be addressed: while considering an active screening strategy in the elderly, would LTBI treatment alleviate a reactivation-dominated TB burden? To assess the balance between the treatment benefit and the risk from hepatitis, what are the benefit-risk ratios? To address these questions, we developed a mathematical model to evaluate the impact of LTBI treatment interventions in the elderly. Materials and Methods Data collection. Annual age-stratified TB notifications (all forms) from 1965?2013, collected from the Department of Health, were used for model fitting. Population data were collected from Census and Statistics Department reports19,20. Model parameter assumptions were obtained from scientific literature (Table?1). The population projection was based on the Census and Statistics Department statistical report21. The ethics approval was obtained from The Joint Chinese University of Hong Kong ? New Territories East Cluster Clinical Research Ethics Committee. Data were retrospectively collected in aggregate form and publicly available so no informed consent was required. All methods were performed in accordance with the relevant guidelines and regulations. Disease transmission model. We extended the age-structured compartmental model to study TB dynamics and the impact of interventions2,3,22,23. Compartments were built based on transmission characteristics. Essential features of the model are shown in Fig.?1. The differential equations describe the transition process of number of individuals between classes given different rates of changes. For model equations, the rate of change in a state is denoted by an upper dot (?) which is the first derivative with respect to time. We stratified the individuals into two age groups: below 60 years old (youngers, a = 1) and greater than or equal to 60 years old (elders, a = 2). For example, S?2 represents the rate of change in the number of susceptibles in the elderly. The age groups are discretized given a younger-to-elder transition rate ?. Denoting the current number of individuals among a population, who are susceptibles (Sa), recent latent TB infections (Lsa), remote latent TB infections (Lla), infectious (Iia), non-infectious (Ina), successfully treated (Tca), treatment defaulted (Tda), receive LTBI treatments adhering with duration j (Tl[j]a), and recovered (Ra) in a-th age group, the model equations are described as follows: Susceptible. According to population dynamics, susceptible individuals were defined as those without TB infection who were not vaccinated or had vaccine failure. The number of susceptibles increased with the number of new birth (?t) by time t. As BCG vaccination was given to newborns and school children since early 1950s, a proportion of them is immunized at a rate ?t at time t given a vaccine efficacy ?. They are infected with an age-specified force of infection ?ta. We assumed only sputum smear-positive individuals can transmit TB in the transmission equations. Net rate of natural mortality and migration in a-th age group (?ta) is included in all compartments. The discretized differential equations read: S?1 = (1 ? ??t)?t ? (?1t + ? + ?1t)S1 S?2 = ?S1 ? (?2t + ?2t)S2 ?1t = ?2?1k NIkit , ?2t = ?2?2k NIkit , k=1 k k=1 k ? M = ????11 ?12??? ? ???21 ?22?? and where ?ak is the age-specified transmission rate from age group k-th to a-th and Ntk is the population size of age group k-th at time t. Because of different contact patterns among youngers and elders, the ?ak forms the ?who acquires infection from whom? (WAIFW) mixing matrix24: Nta ?t ?ta ?t ? ?ta ?ak ? ? ?s ?i ?a ? q p ?c ?d ? ? ? m r hj ?j ?l j Parameter Definition Proportion developing recent latent infection stage from remote 0.0133 latent and recovered stages Population size of a-th age group at time t Number of new births at time t Net rate of natural mortality and migration of a-th age group at time t BCG vaccination proportion on new births at time t BCG efficacy Force of infection for age group a at time t Age-specified transmission rate from age group k-th to a-th Transition rate from youngers to elders Progression rate from recent latent infection to remote latent Progression rate from recent latent infection to infectious or non-infectious stage Progression rate from remote latent infection to infectious or non-infectious stage Proportion of infectious TB from either fast or slow latent TB Proportion of infectious and non-infectious individuals being detected and treated Efficacy for treatment completion Efficacy for treatment defaulted Progression rate from treatment to recovered stage Natural annual recovery rate Relapse rate Screened proportion (acceptability) of LTBI in elders Detection rate of LTBI LTBI treatment adherence level for duration j Progression rate from LTBI treatment with duration j to recovered stage 0.2 per year 0.03 per year22,51,52 ? ? 0.0005 per year for younger adults. A multiplier K to this value was assumed for elders and estimated by the MCMC approach2,22,51,52 Reciprocal for 6 months Test for different proportions 77%, an average sensitivity of the tuberculin skin test16. Enhancing sensitivity of 90% (T-spot TB test) was also tested. Duration j can be 52, 24, and 12 weeks; the current adherence levels were respectively assumed to be 65% (complete treatment), 15%, and 5%16,29,30. We varied the treatment completion rate from 0?80% to assess the changes (24- and 12-weeks completion rate was unchanged). Reciprocal for 52, 24, and 12 weeks of isoniazid therapy16,27,28 Latently infected (recent). Infected individuals move into the recent latent class. Recovered and remote latent TB subjects can also be re-infected at a rate ? and move to this class. The recent latent subjects can develop active disease at a rate ?i or progress to remote latent class at a rate ?s. L? s1 = ?t1S1 + ??t1(R1 + Ll1) ? (?s + ?i + ? + ?t1)Ls1 L? s2 = ?Ls1 + ?t2S2 + ??t2(R2 + Ll2) ? (?s + ?i + ?t2)Ls2 Assumptions Latently infected (remote). Individuals who are not active TB ? 5 years of the recent latent period move to the remote latent stage2,23. We defined TB reactivation as a TB infection that developed ? 5 years after initial infection. A constant reactivation rate was assumed, and we estimated a fractional increase (K) in progression rate from elderly-young individuals. Those progressed to remote latent can develop active disease in either infectious or non-infectious at a rate ?a. As for intervention testing, a proportion m of the latent subjects in elderly are screened at a sensitivity r and treated with LTBI treatment with adherence level of duration js i.e. 52 (whole course), 24, and 12 weeks. A number of treated LTBIs returns to this compartment given a proportion of LTBI treatment failure (1- ?l j) due to limited efficacy levels of different durations. L? l1 = ?sLs1 ? (??t1 + ?1 + ? + ?t1)Ll1 L? l2 = ?Ll1 + ?sLs2 + ?j? j(1 ? ?l j)T l[j]2 ? (??t2 + ?2 + ?jh jmr + ?t2)Ll2 Infectious and non-infectious. Given that infectious TB cases are predominantly sputum smear-positive (microscopic examination), these subjects are defined as being infectious with a proportion ?. Those smear negative or extrapulmonary TB cases are defined as being non-infectious. Primary progression to active TB can occur ? 5 years of the first infection or result from reinfection from recovered and remote latent individuals i.e. the infectious and non-infectious individuals are progressed from recent and remote latently infected. They can be detected and treated with a proportion q, and the remaining individuals (1- q) will be naturally cured at a rate ?. Given limited efficacy for completed treatment (1- ?c) and defaulted treatment (1- ?d), treated individuals can progress to these compartments at a rate ?. Recovered individuals can also relapse and move into these classes. I?i1 = ?[?R1 + ?iLs1 + ?1Ll1 + (1 ? ?c)?T c1 + (1 ? ?d)?T d1] ? [q + ? + (1 ? q)? + ?t1]I i1 I?i2 = ?I i1 + ?[?R2 + ?iLs2 + ?2Ll2 + (1 ? ?c)?T c2 + (1 ? ?d)?T d2] ? [q + (1 ? q)? + ?t2]I i2 and and I?n1 = (1 ? ?)[?R1 + ?iLs1 + ?1Ll1 + (1 ? ?c)?T c1 + (1 ? ?d)?T d1] ? [q + ? + (1 ? q)? + ?t1]I n1 I?n2 = ?I n1 + (1 ? ?)[?R2 + ?iLs2 + ?2Ll2 + (1 ? ?c)?T c2 + (1 ? ?d)?T d2] ? [q + (1 ? q)? + ?t2]I n2 Treatment completion and defaulted. After passive screening, infectious and non-infectious individuals are treated with anti-TB drugs and they can complete treatment (at least 6 months) or default treatment respectively with proportions p and (1-p). The individuals can recover or be failure to treatment at a rate ?. T? c1 = pq(I i1 + I n1) ? (? + ? + ?t1)T c1 T? c2 = ?T c1 + pq(I i2 + I n2) ? (? + ?t2)T c2 LTBI treatment. The screened proportion (m) of remote latent subjects are treated with LTBI therapy with an efficacy ?l j and a recovery rate ?j corresponding to different durations. The hjs are the tested proportions of subjects in different adherence levels for duration j. Subjects with different adherence are stratified into j compartments i.e. Tl[j]. ? l[j]1 = h jmrLl1 ? (? j?l j + ? + ?t )T l[j]1 T 1 ? l[j]2 = ?Tl[j]1 + h jmrLl2 ? (? j?l j + ?t2)T l[j]2 T Recovered. Infected subjects can be cured naturally or by treatments and they can move into the recovered class. Recovered individuals can relapse at a rate ? to either infectious or non-infectious. They can also be re-infected to the latent stage at a rate ?. R? 1 = (1 ? q)?(I i1 + I n1) + ?(?cT c1 + ?dT d1) ? (??1 + ? + ? + ?t1)R1 R? 2 = ?R1 + (1 ? q)?(I i2 + I n2) + ?(?cT c2 + ?dT d2) ? (??2 + ? + ?t2)R2 Natural deaths were considered in all compartments. We assumed no imported cases in the model setting, as migrant infections comprise a low proportion of TB in Hong Kong. Model calibration for baseline scenario. In the initial stages of model calibration, we assumed all of the population in 1965 born before 1952 (age >13 years) were susceptible or infected with TB i.e. ~60% of total population as BCG vaccination has been given to newborns since 195219,20. Following Wu et al.3, we assumed 70% prevalence of LTBI with 20% recent latent infections. The initial numbers of active TB cases were assumed at 2.5 times the number of TB notifications. The baseline scenario was calibrated by fitting annual age-stratified TB notification data to the model. The first 5 years of data (1965?1969) were discarded from model fitting to remove the burn-in period of the model simulations. Fixed parameters are summarised in Table?1. The model generated age-stratified number of TB new cases (for infectious and non-infectious) was fitted to the observed annual notifications by using the Markov Chain Monte Carlo (MCMC) method in a Bayesian inferential framework. A likelihood function was formed by assuming Poisson-distributed TB notifications: L(?) = ?t=1 e??t? xt t xt! where xt and ?t are respectively the observed and the expected new cases generated by the transmission model at time t. In the MCMC estimation, flat prior distributions were assumed for the parameters (?) i.e. ?11, ?12, ?21, ?22, and K. A random walk Metropolis algorithm was used to obtain the posterior distributions. Totally 10,000 iterations were used as the burn-in period and the subsequent 100,000 iterations were used to draw the posterior estimates. Step sizes were selected to obtain acceptance proportions between 20% and 40%. The convergence of the MCMC chains was diagnosed with time-series trace plots and tested by Gelman and Rubin convergence measure25. Auto-correlation of the MCMC chains was checked by auto-correlation function plots. Assuming the Bacillus Calmette?Gu?rin (BCG) vaccination was implemented at a constant rate of 99.4% (average over past 5 years) and no LTBI was treated, a forecast of TB incidence (number of cases/total population size per 100,000) to year 2030 was simulated for the baseline scenario. The variation in births, mortality, and net movement was in accordance with projected statistics. To assess the goodness of fit for the model calibration, R-square, mean absolute error (MAE), and mean absolute percentage error (MAPE) were used. In addition, a retrospective validation process using data from 1965?2003 for model fitting and data from 2004?2013 for validation was conducted to examine the consistency in terms of predictive power. The deterministic equations of the mathematical model and calibration were programmed using R software version 3.0.3 (R Foundation for Statistical Computing, Vienna, Austria). Scenarios for active screening and treatment of TB in the elderly. In Hong Kong, diagnosis of LTBI is commonly established using the tuberculin skin test (TST), followed by 6?12 months treatment with INH offered to all patients. The T-Spot. TB test is an alternative method for LTBI diagnosis, but is not currently in the local strategic protocol16. A shorter course of weekly INH and rifapentine regimen is sometimes offered to elderly patients with a recent study indicated it is as effective as INH treatment alone26. Given the scenario that LTBIs in the elderly group were screened annually and those found positive were treated, we examined the impact of varying proportions of LTBI patients who were willing to be screened (acceptability) and who completed treatment (adherence) on the annual TB incidence starting from 2014. Additional compartments (Tl[j]) were added to the baseline model for LTBI patients who adhered to the j-th treatment duration: 52 (whole course), 24, and 12 weeks (Fig.?1). The diagnostic sensitivity for LTBIs in the elderly was fixed at 76% from TST and 93% from T-spot. TB test17. The current practice for the duration of INH therapy in Hong Kong was fixed at 12 months although WHO recommended 6 months duration of INH is adequate treatment for LTBI27,28. Adherence levels affected treatment effectiveness in the elderly population. Present adherence levels were assumed to be 65%, 15%, and 5% for 52, 24, and 12 weeks? treatment, respectively16,29,30, with corresponding efficacies of 75%, 65%, and 21%, respectively16,29. The averted incidence (i.e. reduction of TB incidence when the LTBI intervention was applied, compared with no intervention scenario in a year) was determined according to differing acceptability and treatment completion rates. As hepatitis is a major safety indicator of INH treatment, the benefit-risk ratio between the number of active TB cases prevented and the number of hepatitis events by year was determined16,29. The number of hepatitis events was calculated as the product between the number of cases treated and the elderly hepatitis rates assumed to be 2.3%, 2.0%, and 1.7% for 52, 24, and 12 weeks of INH treatment, respectively16,29,31. sensitivity analysis. A sensitivity analysis was conducted by obtaining the partial rank correlation coefficients (PRCCs) between the parameter values and the averted TB incidence to explore the sensitivity of the selected model parameters32,33. Plausible distributions of fixed parameters (i.e. BCG vaccination rate, BCG efficacy, proportion developing recent latent infection stage, progression rates from recent and remote latent infection to infectious or non-infectious stage, proportion of infectious TB, proportion of infectious and non-infectious individuals being detected and treated, proportion of treated individuals able to complete TB treatments, efficacies for treatment completion or treatment defaulted, natural recovery rate, and relapse rate) were assumed given plausible ranges of values (Table?1). Fitted parameters were assumed to follow their posterior distributions. In addition, multivariate sensitivity analysis was conducted to assess the impact of varying all selected parameters on the uncertainty of the annual TB incidence. The 95% credible intervals (CI) of the simulation outcomes were drawn over 10,000 realisations. Results Baseline scenario. Figure?2 shows the observed and fitted TB incidence by age and estimated number of reactivated cases. All MCMC chains of the estimated parameters were well-converged given the Gelman and Rubin convergence measure approach to unity. Posterior distributions of the model parameters are shown in Supplementary Information File. The R-square, MAE, and MAPE of the model fitness were 96.4%, 9.6/100,000, and 8.8% respectively showing a high level of agreement between observed data and model simulated TB infection trend. In the retrospective validation using 2004?2013 data, we showed that the R-square, MAE, and MAPE of the model simulation were 89.5%, 6.8/100,000, and 8.2% respectively which indicated a high consistency in model prediction. In the baseline scenario, TB incidence in younger individuals gradually decreased from 222/100,000 in 1970 to 38/100,000 in 2010, while incidence in elderly patients remained at ~30?40/100,000 throughout the study period (Fig.?2). Reactivated cases in the elderly steadily increased from 8/100,000 in 1970 to 30/100,000 in 2010. In the elderly, the proportion of reactivations increased from 31% in 1970 to almost 75% in 2010. From the estimates, the reactivation rate was 3 times higher in elderly patients than young patients (K = 3.33). LTBI interventions. Figure?3 shows the impact of LTBI interventions on changes in annual TB incidence by year, according to various acceptability rates treatment in the elderly, considering the present level of treatment adherence (65% completion rate). If no elderly cases were treated, TB cases decreased slowly from 2014?2030 with ~10/100,000 incidence reduction. However, if a small number of elderly LTBI cases were screened, i.e. 10% annually after 2014, the reduction rate could be similar in 2022, with a 15/100,000 averted incidence in 2030. A 20% treatment acceptability could induce a 20/100,000 decrease in incidence from 2014?2021, while a 40% acceptability could reduce the annual incidence by almost 50% in 2025 (~40/100,000). Increasing the acceptability for LTBI screening to 80% resulted in a substantial decrease in incidence (>20/100,000), but did not generate a large difference in later years, i.e. ~5/100,000 reduction in cases from 2025?2030. LTBI treatment did not substantially impact the overall control when acceptability was set at >40% (Fig.?3B). The difference in TB incidence between 2025 and 2030 was only ~6/100,000 with a treatment acceptability between 40?80%. If the LTBI detection rate was increased to 93% (e.g. using T-spot. TB), the additional effect on incidence reduction was minimal (Fig.?3B). Figure?4 shows the averted incidence in 2020, 2025, and 2030 with varying acceptability and LTBI treatment completion rates. In general, maintaining the treatment completion rate at >60% and increasing acceptability resulted in observable effects on the incidence averted. Compared with the no intervention scenario, if the treatment completion rate was enhanced to 80%, a 40% acceptability obtained ~20 and 30/100,000 averted incidence in 2020 and 2030, respectively. Nevertheless, despite maintaining a high treatment completion rate, a higher acceptability (>40%) could not draw observable long-term impact. An 80% acceptability could yield ~30/100,000 cases averted in 2025, but only yield 4/100,000 additional cases averted in 2030. Meanwhile, the LTBI intervention became less effective for incidence reduction when the treatment completion rate was low. When the treatment completion rate was as low as 20%, increasing the acceptability to 40% in the elderly yielded only ~12/100,000 averted incidence in 2030. Benefit-risk ratio. Given a relatively high hepatitis rate, benefit-risk ratios were kept at a low level when the present level of treatment adherence was maintained (Fig.?3C). Even if the LTBI intervention was given annually with high acceptability, benefit-risk ratios were below unity before 2025, i.e. the number of active TB cases prevented were fewer than the number of hepatitis events. However, preventing a TB case reduced further transmission chains, and the intervention could obtain a benefit-risk ratio greater than unity in 2030, given >30% acceptability. Benefit-risk ratios remained at low levels and did not vary substantially (range, 0.2?0.7) in the early years with varying acceptability and treatment completion rates (Fig.?5). Meanwhile, if current treatment completion rate of the annual LTBI intervention was maintained or increased to >65%, a remarkable gain in benefit-risk ratios was observed in later years, even with moderate acceptability (Fig.?5B,C). sensitivity analysis. A sensitivity analysis was conducted to explore the sensitivity of the model parameters and the tornado diagram showed that the simulated averted incidences were not sensitive when the parameters varied (Fig.?6A). Among the parameters, the progression rate from recent latent infection and proportion of infectious TB had the most influence on the simulated averted incidence (PRCC = 0.107 and 0.077, respectively) as the effect of the TB treatment largely relied on the generated number of infectious cases based on these rate parameters. The absolute numbers of the PRCCs from all the other parameters were all <0.05. Figure?6B,C shows the simulated scenarios for changes in TB incidence when model parameters varied multivariately. In general, the overall trend of the impact from LTBI interventions was consistent with the abovementioned results. The error intervals of the 95% CIs were approximately ?15, ?12, and ?10/100,000 for the TB incidence in 2020, 2025, and 2030, respectively. Discussion The high number of LTBI cases in elderly individuals and the associated risk of reactivation poses a challenge to TB control, with management in this population likely to improve control. Particularly in some regions with low to moderate TB burden, slow decline in TB has been attributed to the cases that develop through reactivation. With a high reactivation rate occurring in the elderly population, a timely evaluation of the effectiveness of LTBI treatment for elderly people is able to offer a valuable advice for TB elimination. In this study, we assessed the impact of LTBI screening and treatment varying the acceptability and adherence from elderly people. Our results showed that at present, reactivations comprise approximately 75% of all TB cases. Given the present treatment completion rate, if ~20?40% of elderly people were willing to be screened and treated for LTBI, the overall TB incidence could be reduced by almost 50% long-term. If the treatment completion rate can be maintained at >65%, increasing acceptability could substantially impact the averted incidence. Nevertheless, this effect tended to taper off in later years. The number of active TB cases prevented were fewer than the number of hepatitis events that occurred in the early years when the LTBI intervention was provided annually, due to hepatotoxicity. Despite high acceptability and treatment adherence, the intervention did not prevent >2 active cases per hepatitis event long-term. As few diagnostic and treatment guidelines targeting TB in the elderly have been established2,34, our study profiles the effectiveness of LTBI interventions in the elderly, and our findings can be used when developing strategies optimised for LTBI screening and treatment. Similar to other studies2,3, our results showed that TB incidence, comprised primarily of elderly reactivation cases, cannot be adequately controlled with current strategies using passive case finding and DOTS. Corroborating an earlier study in Hong Kong3, we showed that approximately 75% of new TB cases were from reactivations, consistent with findings from China35. The TB burden in Japan is similar to Hong Kong, in which 71.8% of new cases in 2016 were attributable to people aged >60 years, who likely had reactivated infection36. In contrast to Hong Kong, the Japanese government has targeted groups at high risk of developing TB in the LTBI treatment program, provided that side-effects and other drug effects were carefully monitored in the elderly receiving the chemoprophylaxis37. To investigate the impact of LTBI measures, Huynh et al.35 developed an individual-based model to investigate whether TB interventions could help achieve the WHO?s End-TB target in China. It was reported that preventative therapy in elderly patients who were screened and treated appropriately for LTBI could help reduce the TB incidence to 60% in 2025. Although the TB epidemiology in China is different from that in Hong Kong, their findings are comparable with those obtained in our study. Although we showed that TB incidence can be substantially reduced long-term, the current strategy for detection and treatment of LTBI will not achieve the End-TB target of 90% reduction or maintenance of ~10/100,000 in Hong Kong by 20358, the current level for low-burden countries38. In general, only a low proportion of patients with LTBIs develop disease. In this study, we found that benefit-risk ratios remained at low levels when the elderly population were screened and treated for LTBI. This indicates that offering LTBI treatments to a general elderly population should be carefully considered. However, LTBI screening and treatments should be prioritized to specific high-risk groups3. WHO defines high-risk groups as individuals with particular pathological conditions, those who are in contact with active TB patients, and those in high-risk living or work settings39. In Hong Kong, a study showed 68.6% of the elderly home residents were positive with the TST screening results7. The estimated LTBI prevalence was comparatively higher than an earlier local study and that in other countries40. Despite these results, the prevalence of hepatotoxicity could be greater among local elderly patients owing to the high prevalence of chronic hepatitis B carriers in Hong Kong7. The optimal approach for delivering interventions, likely to involve stratifying patients according to hepatotoxicity risk and vigilant monitoring during treatment, would require prudent evaluations. Apart from drug toxicity, we showed that poor adherence to LTBI treatment can adversely impact the reduction of TB incidence. We showed that the present LTBI treatment completion rate was able to achieve a desirable averted incidence, and this appears comparable to the 75% rate obtained in a similar setting that offered self-administered INH over 6 months in Australia41. Investigations have shown that the elderly were even more likely to adhere to INH treatment than younger patients42,43. In addition, the estimates in the study offer a valuable reference for the number of LTBIs needed to be approach in achieving a greater benefit with required level of acceptance and tolerance as they can vary by regions and populations according to the Public Health England44. Apart from that, the study findings may give additional values for implying a revision of related policy. For example, a recent published clinical trial showed that 4-month of rifampin was not inferior to the typical 9-month of INH for the prevention of active TB and it was associated with a higher rate of treatment completion (i.e. ~15%) and reduction of hepatotoxic events45. Based on our estimates, such increase in treatment completion would yield >30/100,000 cases averted and a benefit-risk ratio greater than unity in 2025 even with a low acceptability in the elderly population. The current study has several major limitations. In the mathematical model, we did not capture demographic and epidemiological characteristics such as sex, TB-associated mortality, and drug-resistant TB. While we did not modulate TB-associated mortality as an outcome in the simulations, local studies reported that approximately 4% mortality is currently associated with a new TB notification34. Using this figure, >100 deaths could be prevented in 2030 if 20% of elderly patients with LTBI were detected and treated. Moreover, we acknowledged simply using hepatitis events to determine a benefit-risk ratio may not be comprehensive enough when comparing with typical benefit-risk analysis as there are other potential side effects from the treatment (e.g. rash, nausea, malaise, and fever), even though it has been employed in different investigations16,29. For example, Sadatsafavi et al. (2013) conducted a throughout cost-benefit analysis and the risks including all adverse events, TB reactivation, and deaths were quantified in terms of quality-adjusted life-years46. However, we tend to use a simple benefit-risk ratio to study a trade-off between effectiveness and safety of an intervention program within the scope of our study. Also, our modelling results could be overestimated, as a proportion of drug-resistant TB cases might be included in the analyses. In Hong Kong, the prevalence of INH- and/or rifampicin-resistant TB is low; the rates for mono-drug-, multi-drug-, and extensively drug-resistant TB were 7.97% (325 cases), 0.834% (34 cases), and 0.0491% (2 cases) respectively47. No evidence has suggested that prophylactic treatment for LTBI causes drug-resistant bacilli strains48, and drug resistance was more likely due to transmission of resistant strains. In Hong Kong, household transmissions of MDR-TB were found, but were rare49. In conclusion, this study showed that for a city with intermediate TB burden primarily driven by reactivations, if LTBI interventions can be provided to the elderly and maintained at the current treatment completion rate of 65%, treating a low-moderate proportion of individuals with LTBIs annually will reduce the TB incidence by 50%, to reach the 2025 milestone of the global End-TB Strategy. Nevertheless, due to a high risk of hepatotoxicity in the elderly population, benefit-risk ratios were comparatively low; thus, intervention programs for optimising screening and treatment strategies for LTBI should be carefully considered. Our findings underscore the need for new anti-TB drugs/drug regimens with greater efficacy and safety for treating LTBI. 1 0 Acknowledgements The authors thank the external reviewer and editors for their expertise advice on constructing the manuscript. This project forms part of a series of studies commissioned by the Food and Health Bureau of the Hong Kong SAR Government and funded by the Health and Medical Research Fund [Phase 3, 2015-2020]. The opinions and assertions contained in the article are private views of the authors and do not necessarily reflect those of the Centre for Health Protection, Department of Health, Hong Kong. The work is supported by Health and Medical Research Fund [Number: CU-16-C11] and National Natural Science Foundation of China [31871340]. Author Contributions C.C.L. collected the data. K.M.J. manipulated the data. K.C.C., C.C.L., G.C.H.T., M.H.W. and E.K.Y. analyzed, interpreted the data and wrote the manuscript. S.Y.F.L. and H.X. analyzed and verified the results. B.C.Y.Z., W.W.Y., P.H.C. and E.K.Y. revised the manuscript. All authors read and approved the final manuscript. Additional Information Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-019-41256-4. Competing Interests: The authors declare no competing interests. Publisher?s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article?s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article?s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. 1 1 1. Department of Health, Hong Kong. Tuberculosis Manual . ( 2006 ). 2. Vynnycky , E. , Borgdorff , M. W. , Leung , C. C. , Tam , C. M. & Fine , P. E. M. Limited impact of tuberculosis control in Hong Kong: attributable to high risks of reactivation disease . Epidemiol. Infect . 136 , 943 - 952 ( 2008 ). 3. Wu , P. et al. The transmission dynamics of tuberculosis in a recently developed Chinese city . PLOS ONE 5 , e10468 ( 2010 ). 4. Centre for Health Protection, Hong Kong. Statistics on Tuberculosis. Statistics on communicable diseases . 2017 . Available at, http:// www.chp.gov.hk/en/notifiable1/10/26/43.html. ( 2017 ). 5. Public Health England. Tuberculosis in England 2016 report. ( 2016 ). 6. Chan-Yeung , M. et al. Molecular and conventional epidemiology of tuberculosis in Hong Kong: a population-based prospective study . J. Clin. Microbiol . 41 , 2706 - 2708 ( 2003 ). 7. Chan-Yeung , M. et al. Prevalence of tuberculous infection and active tuberculosis in old age homes in Hong Kong . J. Am. Geriatr. Soc . 54 , 1334 - 1340 ( 2006 ). 8. World Health Organization. Implementing the end TB strategy: the essentials . 2015 . Available at, http://www.who.int/entity/tb/ publications/2015/end_tb_essential. pdf?ua=1 . ( 2017 ). 9. Yap , P. et al. Prevalence of and risk factors associated with latent tuberculosis in Singapore: A cross-sectional survey . Int. J. Infect. Dis . 72 , 55 - 62 ( 2018 ). 10. World Health Organization. Tuberculosis country profiles . 2018 . Available at, https://www.who.int/tb/country/data/profiles/en. ( 2018 ). 11. Tuberculosis Surveillance Center. Tuberculosis in Japan - annual report 2017 . ( 2017 ). 12. Houben , R. M. & Dodd , P. J. The global burden of latent tuberculosis infection: a re-estimation using mathematical modelling . PLOS Med . 13 , e1002152 ( 2016 ). 13. Mancuso , J. D. , Diffenderfer , J. M. , Ghassemieh , B. J. , Horne , D. J. & Kao , T. C. The prevalence of latent tuberculosis infection in the United States . Am. J. Respir. Crit. Care Med . 194 , 501 - 509 ( 2016 ). 14. Ghotbi , N. , Nishimura , S. & Takatsuka , N. Japan's national tuberculosis control strategies with economic considerations . Environ. Health Prev. Med . 10 , 213 ( 2005 ). 15. Centers for Disease Control and Prevention. Prevention and control of tuberculosis in facilities providing long-term care to the elderly. Recommendations of the Advisory Committee for Elimination of Tuberculosis . MMWR. 39 , 7 - 20 ( 1990 ). 16. Leung , C. C. et al. Treatment of latent infection with M. tuberculosis: update 2010 . Eur. Respir. J. erj00793-2010 ( 2010 ). 17. Horsburgh , C. R. Jr . et al. Latent TB infection treatment acceptance and completion in the United States and Canada . Chest 137 , 401 - 409 ( 2010 ). 18. Hirsch-Moverman , Y. , Daftary , A. , Franks , J. & Colson , P. W. Adherence to treatment for latent tuberculosis infection: systematic review of studies in the US and Canada . Int. J. Tuberc. Lung. Dis . 12 , 1235 - 1254 ( 2008 ). 19. Census and Statistical Department Hong Kong. Demographic Trends in Hong Kong 1981 - 2011 . 2012 . Available at, http://www. censtatd.gov.hk/hkstat/sub/sp20.jsp?productCode= B1120017 . ( 2017 ). 20. Census and Statistical Department, Hong Kong. Hong Kong Life Tables 2009 - 2064 . 2015 . Available at, http://www.censtatd.gov.hk/ hkstat/sub/sp190.jsp?productCode= B1120016 . ( 2017 ). 21. Census and Statistical Department, Hong Kong. Hong Kong Population Projections 2015 - 2064 . 2015 . Available at, https://www. censtatd.gov.hk/hkstat/sub/sp190.jsp?productCode= B1120015 . ( 2017 ). 22. Dowdy , D. W. , Dye , C. & Cohen , T. Data needs for evidence-based decisions: a tuberculosis modeler's 'wish list' . Int. J. Tuberc. Lung. Dis . 17 , 866 - 877 ( 2013 ). 23. Vynnycky , E. & Fine , P. E. M. The natural history of tuberculosis: the implications of age-dependent risks of disease and the role of reinfection . Epidemiol. Infect . 119 , 183 - 201 ( 1997 ). 24. Dietz , K. & Schenzle , D. Proportionate mixing models for age-dependent infection transmission . J. Math. Biol . 22 , 117 - 120 ( 1985 ). 25. Gelman , A. & Rubin , D. B. Inference from iterative simulation using multiple sequences . Stat. Sci. 7 , 457 - 472 ( 1992 ). 26. Colditz , G. A. et al. Efficacy of BCG vaccine in the prevention of tuberculosis: meta-analysis of the published literature . JAMA . 271 , 698 - 702 ( 1994 ). 27. Cohn , D. L. et al. Targeted tuberculin testing and treatment of latent tuberculosis infection . Am. J. Respir. Crit. Care Med . 161 , 221 - 247 ( 2000 ). 28. Ziakas , P. D. & Mylonakis , E. 4 months of rifampin compared with 9 months of isoniazid for the management of latent tuberculosis infection: a meta-analysis and cost-effectiveness study that focuses on compliance and liver toxicity . Clin. Infect. Dis . 49 , 1883 - 1889 ( 2009 ). 29. International Union Against Tuberculosis Committee on Prophylaxis. Efficacy of various durations of isoniazid preventive therapy for tuberculosis: five years of follow-up in the IUAT trial . Bull. World Health Organ . 60 , 555 ( 1982 ). 30. Sprinson , J. E. et al. Evaluation of tuberculosis contact investigations in California . Int. J. Tuberc. Lung. Dis . 7 , S363 - S368 ( 2003 ). 31. Kopanoff , D. E. , Snider , D. E. Jr. & Caras , G. J. Isoniazid-related hepatitis: a US Public Health Service cooperative surveillance study . Am. Rev. Respir. Dis . 117 , 991 - 1001 ( 1978 ). 32. Marino , S. , Hogue , I. B. , Ray , C. J. & Kirschner , D. E. A methodology for performing global uncertainty and sensitivity analysis in systems biology . J. Theor. Biol . 254 , 178 - 196 ( 2008 ). 33. Shrestha , S. et al. Impact of Targeted Tuberculosis Vaccination Among a Mining Population in South Africa: A Model-Based Study . Am. J. Epidemiol . 186 , 1362 - 1369 ( 2017 ). 34. Li , J. et al. The strategic framework of tuberculosis control and prevention in the elderly: a scoping review towards End TB targets . Infect. Dis. Poverty 6 , 70 ( 2017 ). 35. Huynh , G. H. et al. Tuberculosis control strategies to reach the 2035 global targets in China: the role of changing demographics and reactivation disease . BMC Med . 13 , 88 ( 2015 ). 36. Seto , J. et al. Mycobacterium tuberculosis transmission among elderly persons , Yamagata Prefecture, Japan , 2009 - 2015 . Emerg. Infect. Dis . 23 , 448 ( 2017 ). 37. Prevention Committee of the Japanese Society for Tuberculosis, & Treatment Committee of the Japanese Society for Tuberculosis. Treatment guidelines for latent tuberculosis infection . Kekkaku:[Tuberculosis] 89 , 21 ( 2014 ). 38. World Health Organization. Global tuberculosis report 2017 . ( 2017 ). 39. World Health Organization. Guidelines on the management of latent tuberculosis infection . 2015 . Available at, http://apps.who.int/ iris/bitstream/10665/136471/1/9789241548908_eng. pdf?ua=1&ua=1 . ( 2017 ) 40. Vega , R. A. , Conde , J. G. & Diaz , M. Prevalence of tuberculin reactivity and prevalence of risk factors for the development of active tuberculosis in a nursing home in Puerto Rico . P. R. Health Sci. J . 15 , 27 - 31 ( 1996 ). 41. Dobler , C. C. & Marks , G. B. Completion of treatment for latent tuberculosis infection with monthly drug dispensation directly through the tuberculosis clinic . PLOS ONE 7 , e48900 ( 2012 ). 42. Colson , P. W. et al. Acceptance of treatment for latent tuberculosis infection: prospective cohort study in the United States and Canada . Int. J. Tuberc. Lung. Dis . 17 , 473 - 479 ( 2013 ). 43. Bock , N. N. , Metzger , B. S. , Tapia , J. R. & Blumberg , H. M. A tuberculin screening and isoniazid preventive therapy program in an inner-city population . Am. J. Respir. Crit. Care Med . 159 , 295 - 300 ( 1999 ). 44. Public Health England. Tuberculosis in England 2017 report (presenting data to end of 2016 ); 2017 . Available at, https://assets.publishing. service.gov.uk/government/uploads/system/uploads/attachment_data/file/686185/TB_Annual _Report_2017_v1.1 .pdf. ( 2018 ). 45. Menzies , D. et al. Four months of rifampin or nine months of isoniazid for latent tuberculosis in adults . N. Engl. J. Med . 379 , 440 - 453 ( 2018 ). 46. Sadatsafavi , M. et al. A quantitative benefit-risk analysis of isoniazid for treatment of latent tuberculosis infection using incremental benefit framework . Value Health 16 , 66 - 75 ( 2013 ). 47. Centre for Health Protection Hong Kong . Antimicrobial resistance profiles of Mycobacterium tuberculosis . Statistics on laboratory surveillance, Antimicrobial susceptibility testing . 2017 , Available at, http://www.chp.gov.hk/en/data/1/10/641/695/4844.html. ( 2017 ) 48. Hong Kong Chest Service/Tuberculosis Research Centre, Madras/British Medical Research Council. A double-blind placebocontrolled clinical trial of three antituberculosis chemoprophylaxis regimens in patients with silicosis in Hong Kong . Am. Rev. Respir. Dis . 145 , 36 - 41 ( 1992 ). 49. Leung , E. C. C. et al. Transmission of multidrug-resistant and extensively drug-resistant tuberculosis in a metropolitan city . Eur. Respir. J. erj00712-2012 ( 2012 ). 50. Department of Health, Hong Kong. Annual report 1949-2013 tuberculosis & chest service . ( 2013 ). 51. Horsburgh , C. R. Jr . et al. Revisiting rates of reactivation tuberculosis: a population-based approach . Am. J. Respir. Crit. Care Med . 182 , 420 - 425 ( 2010 ). 52. Cohen , T. , Colijn , C. , Finklea , B. & Murray , M. Exogenous re-infection and the dynamics of tuberculosis epidemics: local effects in a network model of transmission . J. R. Soc. Interface 4 , 523 - 531 ( 2006 ). 53. Tam , C. The DOTS Strategy in Hong Kong . Hong Kong Med. Diary 11 , 3 - 4 ( 2006 ). 54. Tam , C. M. , Leung , C. C. , Noertjojo , K. , Chan , S. L. & Chan-Yeung , M. Tuberculosis in Hong Kong-patient characteristics and treatment outcome . Hong Kong Med. J . 9 , 83 - 90 ( 2003 ). 55. Wong , M. Y. , Leung , C. C. , Tam , C. M. & Lee , S. N. Directly observed treatment of tuberculosis in Hong Kong . Int. J. Tuberc. Lung. Dis . 9 , 443 - 449 ( 2005 ).


This is a preview of a remote PDF: https://www.nature.com/articles/s41598-019-41256-4.pdf

Ka Chun Chong, Chi Chiu Leung, Wing Wai Yew, Benny Chung Ying Zee, Greta Chun Huen Tam, Maggie Haitian Wang, Katherine Min Jia, Pui Hong Chung, Steven Yuk Fai Lau, Xiaoran Han, Eng Kiong Yeoh. Mathematical modelling of the impact of treating latent tuberculosis infection in the elderly in a city with intermediate tuberculosis burden, Scientific Reports, 2019, DOI: 10.1038/s41598-019-41256-4