Rising challenge of multidrug-resistant tuberculosis in China: a predictive study using Markov modeling
Infectious Diseases of Poverty volume 9, Article number: 65 (2020)
Multidrug-resistant tuberculosis (MDR-TB) is on the rise in China. This study used a dynamic Markov model to predict the longitudinal trends of MDR-TB in China by 2050 and to assess the effects of alternative control measures.
Eight states of tuberculosis transmission were set up in the Markov model using a hypothetical cohort of 100 000 people. The prevalence of MDR-TB and bacteriologically confirmed drug-susceptible tuberculosis (DS-TB+) were simulated and MDR-TB was stratified into whether the disease was treated with the recommended regimen or not.
Without any intervention changes to current conditions, the prevalence of DS-TB+ was projected to decline 67.7% by 2050, decreasing to 20 per 100 000 people, whereas that of MDR-TB was expected to triple to 58/100 000. Furthermore, 86.2% of the MDR-TB cases would be left untreated by the year of 2050. In the case where MDR-TB detection rate reaches 50% or 70% at 5% per year, the decline in prevalence of MDR-TB would be 25.9 and 36.2% respectively. In the case where treatment coverage was improved to 70% or 100% at 5% per year, MDR-TB prevalence in 2050 would decrease by 13.8 and 24.1%, respectively. If both detection rate and treatment coverage reach 70%, the prevalence of MDR-TB by 2050 would be reduced to 28/100 000 by a 51.7% reduction.
MDR-TB, especially untreated MDR-TB, would rise rapidly under China’s current MDR-TB control strategies. Interventions designed to promote effective detection and treatment of MDR-TB are imperative in the fights against MDR-TB epidemics.
Tuberculosis (TB), once called the “white plague”, remains one of the leading causes of deaths from infectious diseases. In 2017 alone, TB caused 10.0 million incidents and resulted in 1.3 million deaths worldwide . Over the past two decades, however, both the incidence and mortality of TB have decreased remarkably thanks to the global efforts against TB. Regardless, drug-resistant tuberculosis, especially multidrug-resistant tuberculosis (MDR-TB), remains a major obstacle in the task of putting an end to the global TB epidemic [2, 3]. The most recent World Health Organization (WHO) report suggests that MDR-TB and rifampicin-resistant tuberculosis (RR-TB) together accounted for 558 000 incidents and 230 000 deaths in 2017 . In the meantime, numerous countries across the world have witnessed a growth in the prevalence of MDR-TB over time, which poses a greater threat to the global prospects of TB control [4, 5].
Given the high priority to TB control, several mathematic models have been developed to estimate the burden of TB [6,7,8,9]. Although a few models have projected the burden of TB in different countries, few gave MDR-TB a specific estimation [10, 11]. Among the proposed forecast models, Markov process is a widely used mathematical stochastic process with the capability to simulate the progress of disease according to certain probabilistic manner between different states [12,13,14].
China has the second highest case burden of MDR-TB worldwide. In 2017, China had an estimated 58 000 MDR-TB/RR-TB incidents, accounting for approximately 10% of the global burden. However, only 22.5% of the above-mentioned estimated MDR-TB/RR-TB cases were laboratory-confirmed and only 10.2% started on treatment, presenting a much lower MDR-TB detection rate and treatment coverage compared to the global average . In the past decades, although China has made remarkable successes in TB control, these achievements might be compromised by the rising prevalence of drug resistance. Illustrating the trajectory of MDR-TB burden will provide a full view of China’s TB epidemics and give insights on policy and strategy development in line with the End TB Strategy, which aims to reduce global TB deaths by 95% and incidence by 90% by the year of 2035 . As the main source of infection transmission, bacteriologically confirmed TB plays a key role in tuberculosis control. Herein, a dynamic Markov model is constructed based on the prevalence of bacteriologically confirmed TB in the national survey to predict the prevalence of MDR-TB and bacteriologically confirmed drug-susceptible TB (DS-TB+) longitudinally under the current China TB/MDR-TB epidemics; and to estimate possible gains from the enhancement of MDR-TB case finding and treatment coverage by 2050.
We constructed a population-level, dynamic and compartmental Markov model in which individuals are grouped into mutually exclusive compartments based on the epidemiological characteristics of TB infections and early model studies [7, 8]. The model, subdivided into MDR-TB and DS-TB+, explicated the occurrence of TB transmission, the development from latent tuberculosis infection (LTBI) to active TB, anti-TB treatment, and treatment outcomes. It simulated the transition among eight states: (1) uninfected; (2) latent infection of drug-susceptible tuberculosis (DS LTBI); (3) latent infection of MDR-TB (MDR LTBI); (4) DS-TB+ untreated; (5) MDR-TB untreated; (6) DS-TB+ treated; (7) MDR-TB treated; (8) deceased (defined in Fig. 1). These eight states were described in Table 1 while the eighteen parameters (from a to r) for transmission of states were presented in Table 2. If patients didn’t receive the recommended treatment, TB was considered untreated. HIV was not explicitly considered because of the low percentage of TB/HIV co-infection (roughly 1%) in China .
The possibility of moving from DS-TB+ untreated to MDR-TB untreated was not considered because the treatment coverage of DS-TB reached nearly 90% in China .
Individuals with active TB disease cannot be re-infected before recovering from the current disease episode .
Model initialization and parameter setting
Simulations of the Markov model were initiated in 2011 for the use of the Fifth National Tuberculosis Epidemiological Survey in 2010, the most recent survey with a sample size of 252 940 by random sampling . At the base year of modeling, a hypothetical cohort of 100 000 Chinese was established and divided into eight TB states. Based on the major surveys in China, the prevalence of LTBI and DS-TB+ were 20% and 66/100 000 respectively, and the proportion of MDR-TB accounted for 6.8% of the total TB cases [16,17,18]. The description of each state, the initial data set and data sources were listed in Table 1.
Transition between states depends on the parameters which determine the pathway of movement within the model in each yearly cycle. Eighteen parameters were defined and extracted from national annual statistics or TB registration, and the estimation of parameters with uncertainty was based on the WHO TB database, the WHO annual report on TB, and published literature [1, 17, 19,20,21,22,23,24,25]. Most of these parameters were held fixed throughout our analysis, except for b and c, which varied with TB prevalence. Table 2 shows details on these probabilities.
Using the model, we projected the impact of three hypothetical scenarios on the epidemics of MDR-TB, including 1) improving detection rate; 2) expanding treatment coverage; 3) increasing detection rate and treatment coverage. Detection rate was the number of diagnosed MDR-TB cases divided by the estimated TB cases. Treatment coverage was expressed as the percentage of notified MDR-TB patients who commenced the WHO recommended MDR-TB treatment. In scenario 1, detection rate was assumed to increase to 50% or 70% at 5% per year from 2020 (baseline is 22.5%). In scenario 2, treatment coverage of MDR-TB was anticipated to reach 70% or 100% at 5% per year from 2020 (baseline is 45.5%). In scenario 3, both detection rate and treatment coverage were assumed to be 70% with a 5% annual increase from the baseline.
We predicted TB prevalence trends from 2019 to 2050, under the current TB control conditions and three intervention schemes. The prevalence of treated MDR-TB was defined as the number of cases enrolled in MDR-TB treatment per 100 000 people. Untreated MDR-TB prevalence was the number of MDR-TB cases left untreated per 100 000 people. The effect of the above-mentioned three scenarios was measured via reduced MDR-TB prevalence compared with the baseline.
Tree Age Pro 2011 software (Tree Age Software, Inc., Williamstown, MA, USA) and Microsoft Excel 2016 (Microsoft Corporation, Redmond, WA, USA) were used to conduct the Markov chain models. GraphPad Prism 6.0 (GraphPad Software, Inc., La Jolla, CA, USA) was used to draw the figures.
TB epidemic trends providing no intervention changes to current TB control conditions
Without any intervention changes to the current conditions, the model projected a considerable decrease in DS-TB+ prevalence and a substantial increase in MDR-TB prevalence from 2019 to 2050. In this case, treated and untreated DS-TB+ were not separated considering the high treatment coverage of DS-TB+ in China. It was estimated that, by 2050, DS-TB+ prevalence would decline to 20 per 100 000, demonstrating a 67.7% reduction from 2019. Conversely, MDR-TB prevalence would increase by three times approximately from 19/100 000 in 2019 to 58/100 000 in 2050. It was further estimated that the prevalence of untreated MDR-TB would increase from 16/10 000 in 2019 to 50/100 000 in 2050, whereas a limited increase would be seen in treated MDR-TB, resulting in a prevalence of 8/100 000 in 2050. By the year of 2050, 74.4% of TB cases would be MDR-TB, 86.2% of which would remain untreated. These untreated MDR-TB cases would account for most TB cases and remain the main reason behind the aggravation of the ongoing TB epidemic. Figure 2 discloses in details the simulation results of the model under baseline condition.
MDR-TB epidemic trends with increased detection rate or treatment coverage
Scenario 1 presented detection rate of MDR-TB increased to 50% or 70% at 5% per year, leading to a projected MDR-TB prevalence with an upward trend much slower than that of the baseline. Under the same scenario, MDR-TB prevalence was estimated to be 43/100 000 or 37/100 000 by 2050, with a detection rate of 50% or 70%, respectively (Fig. 3). Compared to the baseline, MDR-TB prevalence displayed a predicted reduction of 25.9 and 36.2% respectively, whereas 64.9% of MDR-TB cases would remain untreated with the recommended MDR-TB regimens even with a 70% detection rate.
With MDR-TB treatment coverage improved to 70% or 100% at 5% per year as indicated in Scenario 2, MDR-TB prevalence was estimated to be 50 or 44 per 100 000 people by 2050, with a treatment coverage of 70% or 100% (reduced by 13.8 and 24.1%), respectively (Fig. 3). whereas 72.7% of MDR-TB cases would remain untreated due to being undetected even if with a treatment coverage of 100%.
TB epidemic trends with simultaneously increased detection and treatment coverage
In scenario 3 where both detection rate and treatment coverage were set to be 70% with an annual increase of 5% from baseline, the prevalence of MDR-TB would be reduced remarkably and was estimated to reach 28/100 000 or lower, demonstrating a 51.7% decrease compared to the baseline. And the proportion of untreated cases would be 53.6% by 2050. Figure 4 shows the trends of MDR-TB in Scenario 3 compared with the baseline.
A Markov model was constructed in this study to forecast trends of DS-TB+ and MDR-TB in China—a country with heavy burden of TB/MDR-TB. It was estimated that, if the present TB control measures remained unchanged, DS-TB+ prevalence would see a significant decrease, but the MDR-TB epidemic would grow at an alarming rate, indicating a gradual shift in the TB epidemics in China from drug-susceptible tuberculosis (DS-TB) to MDR-TB. In addition, most MDR-TB cases would remain untreated, and thus would continue to transmit the disease. The results suggested that increasing detection rate of MDR-TB could have an obvious effect on MDR-TB prevalence but could not interrupt the emerging trend of MDR-TB prevalence. Meanwhile, with the present MDR-TB detection rate in China being significantly low, MDR-TB prevalence would continue to increase even if treatment coverage improved to 100%. Therefore, detection rate must be enhanced alongside treatment coverage to achieve a significantly reduction in MDR-TB prevalence.
A Markov model was built in the study that separated TB prevalence into DS-TB and MDR-TB sub-categories to achieve in-depth observations and detailed predictions for the TB epidemic in China. The projected DS-TB+ and MDR-TB prevalence was consistent with those projected in previous studies [28, 29], indicating that MDR-TB prevalence would continue to rise. In addition, simulations suggested that most MDR-TB cases would remain untreated, illustrating the urgency for improving MDR-TB detection and treatment. In the meantime, this study constituted as first step to better understanding the trends of both untreated and treated MDR-TB in China and the outcome is also applicable for other high-burden, low-to-middle income countries.
Increasing MDR-TB detection rate from 22 to 70% is challenging but achievable. The low MDR-TB detection rate resulted from a combination of issues, including a lack of laboratory capacity for drug susceptibility test (DST), insufficient funding and inequity in health financing, high cost associated with rapid diagnostics, and low clinical alert for TB . Under China’s current TB control program, health financing, health insurance and health care for MDR-TB are localized and determine the elementary medical services and expenses. The lack of DST resulted in low case detection in some TB-designated health facilities on the county and district level, who in return had to deliver sputum samples to upper-level facilities for DST . In addition, restricted by resources and the capacity in MDR-TB health care, drug-resistance testing is only recommended to patients at high risk of MDR-TB under China’s current MDR-TB control program, whereas these patients merely account for less than 50% of all MDR-TB cases . Even when DST is available, the results can take months to come back, meanwhile, patients would have already been treated with DS-TB regimen . Unfortunately, the use of rapid molecular diagnostics is limited due to their costs and the fact they are beyond health insurance coverage in many provinces. In conclusion, the lack of diagnostic capacity, especially the ability to provide rapid diagnosis, is a major obstacle preventing proper, effective MDR-TB detection in China. In June 2019, the Chinese government announced Action Plan to Stop Tuberculosis (2019–2022). The plan aims to enhance laboratory capacity for MDR-TB diagnosis, requiring a DST coverage higher than 90% among bacteriologically confirmed cases by 2022. The strategy includes developing rapid molecular diagnostics for TB drug susceptibility and training procedures for associated laboratory staff, potentially allowing a significant acceleration in MDR-TB diagnosis.
In addition to case detection, treatment coverage is also essential to the effective control of MDR-TB. Multiple reasons lead to the difficulty of treatment enrollment, including economic hardship, loss to follow-up or death prior to treatment initiation, domestic migration caused by work and the misbelief that MDR-TB could be cured with DS-TB treatment. Whereas, economic hardship is the most common reason due to the high cost for medicines. While the expenses can be partially covered by health insurance, MDR-TB patients are still responsible for a considerable amount of treatment-related expenses that tend to be unaffordable for low-income patients. Findings from previous studies proved that 38% of MDR-TB patients experienced catastrophic healthcare expenses (more than 20% of annual household income) . Consequently, it would be unrealistic to expect patients to bear the high-cost of MDR-TB care by themselves without risking catastrophic expenditures. Whereas, the decline in international donor funding now leaves it to the government to increase the amount of internal funding for the war against the MDR-TB epidemic.
The present study does show restrictions. For example, the general model was built on current estimates without taking new actions and interventions into consideration, such as policy changes, the introduction of new diagnostics, novel medicines, etc. In addition, some parameter values could not be identified accurately and thus, estimates based on several individual studies conducted in different years were used, which may affect the model’s accuracy. Therefore, alternative models projecting the MDR-TB epidemic should be considered and the corresponding result should be compared with the findings in the present study to gain a comprehensive understanding of the trends of TB/MDR-TB in China. Third, the constructed simple Markov model was unable to consider every detail of the transition process. In the present analysis, certain less influential parameters were simplified to perform the estimation, including the progression from DS-TB+ untreated to MDR-TB untreated and TB recurrence.
The present study illustrated the challenges of MDR-TB confronted by China’s TB control program and provided unique insights into MDR-TB intervention. While the prevalence of DS-TB+ displayed a significant decrease, the MDR-TB epidemic would continue to grow rapidly. Improving MDR-TB detection and treatment coverage independently could reduce MDR-TB prevalence altogether, but was unlikely to halt the rise of MDR-TB. The results above indicated that enhanced interventions on promoting MDR-TB detection and treatment simultaneously would be imperative to stop the growing MDR-TB epidemic in China. Such enhanced interventions should involve increased emphasis on DST accessibility for all bacteriologically confirmed TB cases, rapid diagnosis, expanded health insurance coverage, and out-of-pocket payments exemption.
Availability of data and materials
All data generated or analyzed during the present study was included in this published article.
- DS-TB+ :
Bacteriologically confirmed drug-susceptible tuberculosis
World Health Organization
Latent tuberculosis infection
- DS LTBI:
Latent infection of drug-susceptible tuberculosis
- MDR LTBI:
Latent infection of multidrug-resistant tuberculosis
Drug susceptibility test
World Health Organization. Global Tuberculosis Report 2018. Geneva: World Health Organization; 2018. https://apps.who.int/iris/bitstream/handle/10665/274453/9789241565646-eng.pdf. Accessed 30 Mar 2019.
Nathanson E, Nunn P, Uplekar M, Floyd K, Jaramillo E, Lönnroth K, et al. MDR tuberculosis–critical steps for prevention and control. N Engl J Med. 2010;363:1050–8.
Gandhi NR, Nunn P, Dheda K, Schaaf HS, Zignol M, Van Soolingen D, et al. Multidrug-resistant and extensively drug-resistant tuberculosis: a threat to global control of tuberculosis. Lancet. 2010;375:1830–43.
Uplekar M, Weil D, Lonnroth K, Jaramillo E, Lienhardt C, Dias HM, et al. WHO's new end TB strategy. Lancet. 2015;385:1799–801.
Millard J, Ugarte-Gil C, Moore DA. Multidrug resistant tuberculosis. BMJ. 2015;350:882.
Lin HH, Wang L, Zhang H, Ruan Y, Chin DP, Dye C. Tuberculosis control in China: use of modelling to develop targets and policies. Bull World Health Organ. 2015;93:790–8.
Sharma A, Hill A, Kurbatova E, van der Walt M, Kvasnovsky C, Tupasi TE, et al. Estimating the future burden of multidrug-resistant and extensively drug-resistant tuberculosis in India, the Philippines, Russia, and South Africa: a mathematical modelling study. Lancet Infect Dis. 2017;17:707–15.
Menzies NA, Cohen T, Hill AN, Yaesoubi R, Galer K, Wolf E, et al. Prospects for tuberculosis elimination in the United States: results of a transmission dynamic model. Am J Epidemiol. 2018;187:2011–20.
Law S, Piatek AS, Vincent C, Oxlade O, Menzies D. Emergence of drug resistance in patients with tuberculosis cared for by the Indian health-care system: a dynamic modelling study. Lancet Public Health. 2017;2:47–55.
McBryde ES, Meehan MT, Doan TN, Ragonnet R, Marais BJ, Guernier V, et al. The risk of global epidemic replacement with drug-resistant Mycobacterium tuberculosis strains. Int J Infect Dis. 2017;56:14–20.
Yu Y, Shi Y, Yao W. Dynamic model of tuberculosis considering multi-drug resistance and their applications. Infect Dis Model. 2018;3:362–72.
Wu J, Zhou Y, Fu X, Deng M, Zheng Y, Tian G, et al. The burden of chronic hepatitis c in China from 2004 to 2050: an individual-based modeling study. Hepatology. 2019;69:1442–52.
Estes C, Anstee QM, Arias-Loste MT, Bantel H, Bellentani S, Caballeria J, et al. Modeling NAFLD disease burden in China, France, Germany, Italy, Japan, Spain, United Kingdom, and United States for the period 2016–2030. J Hepatol. 2018;69:896–904.
Sun L, Legood R, Sadique Z, dos -Santos-Silva I, Yang L. Cost–effectiveness of risk-based breast cancer screening programme, China. Bull World Health Organ. 2018;96:568.
World Health Organization. The End TB Strategy. Geneva, Switzerland: World Health Organization. 2013. https://www.who.int/tb/End_TB_brochure.pdf?ua=1. Accessed 30 Mar 2019.
Wang LX, Cheng SM, Chen MT, Zhao YL, Zhang H, Jiang SW, et al. The fifth national tuberculosis epidemiological survey in 2010. Chin J Antituberc. 2012;34:485–508.
Gao L, Li X, Liu J, Wang X, Lu W, Bai L, et al. Incidence of active tuberculosis in individuals with latent tuberculosis infection in rural China: follow-up results of a population-based, multicentre, prospective cohort study. Lancet Infect Dis. 2017;17:1053–61.
Li X, Yang Q, Feng B, Xin H, Zhang M, Deng Q, et al. Tuberculosis infection in rural labor migrants in Shenzhen, China: Emerging challenge to tuberculosis control during urbanization. Sci Rep. 2017;7:4457.
National Bureau of Statistics of China. China Statistical Yearbook 2017. China Statistical Press 2017. Available: http://www.stats.gov.cn/tjsj/ndsj/2017/indexch.htm. Accessed 15 Mar 2019.
Huang LQ, Wang WB, Li HD, Xu B, Yang HX. A descriptive study on diagnosis delays and factors impacting accessibility to diagnosis among TB patients in floating population in Shanghai. Chin J Antituberc. 2007;29:127–9.
Sun CF, Hua LQ, Xiong LP. Negative conversion of bacillus positive sputum in patients with pulmonary tuberculosis: an analysis of related factors. Acad J Sec Mil Med Univ. 2004;25:1272–3.
Xu KJ. The cost-effectiveness analysis of active chest X ray screening for tuberculosis in mass population and the study of distinguish NTM and MDR-TB from tuberculosis in China. Zhejiang University; 2014.
Mandal P, Craxton R, Chalmers JD, Gilhooley S, Laurenson IF, McSparron C, et al. Contact tracing in pulmonary and non-pulmonary tuberculosis. QJM. 2012;105:741–7.
Xu J, Hu Y, Jiang W, Hong J, Tang L, Yang M, et al. Prevalence and risk factors of latent tuberculosis infection in close contacts of tuberculosis patients among non-resident populations in Shanghai, China. Chin J Tuberc Respir Dis. 2016;39:25–9.
World Health Organization. WHO Tuberculosis (TB) database downloads. http://www.who.int/tb/country/data/download/en/. Accessed 5 Mar 2019.
Wu P, Lau EH, Cowling BJ, Leung CC, Tam CM, Leung GM. The transmission dynamics of tuberculosis in a recently developed Chinese city. PLoS One. 2010;5(5):10468.
van Rie A, Warren R, Richardson M, Victor TC, Gie RP, Enarson DA, et al. Exogenous reinfection as a cause of recurrent tuberculosis after curative treatment. N Engl J Med. 1999;341(16):1174–9.
Houben RM, Menzies NA, Sumner T, Huynh GH, Arinaminpathy N, Goldhaber-Fiebert JD, et al. Feasibility of achieving the 2025 WHO global tuberculosis targets in South Africa, China, and India: a combined analysis of 11 mathematical models. Lancet Glob Health. 2016;4:806–15.
Xu K, Ding C, Mangan CJ, Li Y, Ren J, Yang S, et al. Tuberculosis in China: a longitudinal predictive model of the general population and recommendations for achieving WHO goals. Respirology. 2017;22:1423–9.
Yang Z, Zhou C, Ning Z, Lu W, Zhao Q, Hu Y, et al. High-risk screening and detection of multidrug-resistant tuberculosis in two prefectures of China: a drug susceptibility surveillance-based secondary data analysis. Glob Health Action. 2018;11:1500763.
Qian L, Yan Q, Lucas H. Drug-resistant tuberculosis control in China: progress and challenges. Infect Dis Poverty. 2016;5(1):1–7.
Chen Y, Yuan Z, Shen X, Wu J, Wu Z, Xu B. Time to multidrug-resistant tuberculosis treatment initiation in association with treatment outcomes in Shanghai, China. Antimicrob Agents Chemother. 2018;62:2259–17.
Li R, Ruan Y, Sun Q, Wang X, Chen M, Zhang H, et al. Effect of a comprehensive programme to provide universal access to care for sputum-smear-positive multidrug-resistant tuberculosis in China: a before-and-after study. Lancet Global Health. 2015;3:217–28.
The study was supported by the Swedish Research Council (No. 540–2013-8797 to SH) and the National Natural Science Foundation of China (No. 81361138019 to BX) joint project (VR-NSFC), and the National Science and Technology Major Project of China (No. 2018ZX10715012–4 to BX).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
About this article
Cite this article
Li, BY., Shi, WP., Zhou, CM. et al. Rising challenge of multidrug-resistant tuberculosis in China: a predictive study using Markov modeling. Infect Dis Poverty 9, 65 (2020). https://doi.org/10.1186/s40249-020-00682-7