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.
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?
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???
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:
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
BCG vaccination proportion on
new births at time t
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
Progression rate from recent
latent infection to remote latent
Progression rate from recent
latent infection to infectious or
Progression rate from remote
latent infection to infectious or
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
(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
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
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
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
Reciprocal for 52, 24, and 12 weeks of isoniazid
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
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
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
? l[j]2 = ?Tl[j]1 + h jmrLl2 ? (? j?l j + ?t2)T l[j]2
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
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
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
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.
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
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%
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.
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
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
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.
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 .
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.
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
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. 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 ).