- Open Access
Nomogram to predict multidrug-resistant tuberculosis
Annals of Clinical Microbiology and Antimicrobials volume 19, Article number: 27 (2020)
Multidrug-resistant tuberculosis (MDR-TB) is burgeoning globally, and has been a serious challenge in TB management. Clinically, the ability to identify MDR-TB is still limited, especially in smear-negative TB. The aim of this study was to develop a nomogram for predicting MDR-TB.
Demographics and clinical characteristics of both MDR-TB and drug-susceptible TB patients were utilized to develop a nomogram for predicting MDR-TB. The LASSO regression method was applied to filter variables and select predictors, and multivariate logistic regression was used to construct a nomogram. The discriminatory ability of the model was determined by calculating the area under the curve (AUC). Moreover, calibration analysis and decision curve analysis (DCA) of the model were performed. This study involved a second analysis of a completed prospective cohort study conducted in a country with a high TB burden.
Five variables of TB patients were selected through the LASSO regression method, and a nomogram was built based on these variables. The predictive model yielded an AUC of 0.759 (95% CI, 0.719–0.799), and in the internal validation, the AUC was 0.757 (95% CI, 0.715–0.793). The predictive model was well-calibrated, and DCA showed that if the threshold probability of MDR-TB was between 70 and 90%, using the proposed nomogram to predict MDR-TB would obtain a net benefit.
In this study, a nomogram was constructed that incorporated five demographic and clinical characteristics of TB patients. The nomogram may be of great value for the prediction of MDR-TB in patients with sputum-free or smear-negative TB.
Tuberculosis (TB) continues to be a heavy burden globally, and alarmingly, the epidemic of resistance is burgeoning . Multidrug-resistant TB (MDR-TB) was defined as resistance to at least isoniazid and rifampin. There were approximately 458,000 new prevalent cases of multidrug-resistant TB (MDR-TB) globally in 2017 . The mortality of MDR-TB in endemic countries or regions using traditional regimens reached 40% . It is much costlier to treat MDR-TB than drug-susceptible TB (DS-TB) [1, 3]. There is a difference in regimen and management between MDR-TB and DS-TB, and therefore it is crucial to identify MDR-TB. Although the diagnostic capability of MDR-TB increased by the phenotypic drug susceptibility testing (DST), array-based platforms and line probe assays [4, 5], missed diagnosis or a delayed diagnosis of MDR-TB is still common in clinical practice, especially in sputum-free or smear-negative TB .
In general, MDR-TB is mainly acquired (caused by improper treatment programs or poor patient treatment compliance); however, primary transmission of MDR-TB could be a dominant mode of spread in epidemics [7, 8]. It has been revealed that several clinical, environmental, and socioeconomic characteristics were different between cases with MDR-TB and DS-TB . In addition, it has previously been pointed out that the pathogenicity of DR and DS Mycobacterium tuberculosis differs in both animal models and human patients [9,10,11]. To accurately identify MDR-TB and start an appropriate treatment regimen is essential in MDR-TB control. However, to the best of our knowledge, there is currently no model available for the prediction of MDR-TB.
In the present study, we developed a nomogram for predicting MDR-TB based on demographic and clinical characteristics of TB patients, which were collected from a completed 3-year prospective cohort study.
Materials and methods
Study population and ethics
This was an analytical study of a previously completed prospective cohort study that was conducted in Peru, which is a country with a high TB burden as defined by the WHO . In the present study, a total of 700 confirmed TB patients were enrolled between September 2010 and September 2013 in Peru . MDR-TB referred to resistance to at least rifampicin and isoniazid, whereas DS-TB referred to susceptibility to both rifampicin and isoniazid . Patients who expectorated sputum had their sputum tested by DST for rifampicin and isoniazid using the microscopic observation drug susceptibility assay (MODS) [12, 13], in which MODS testing indicated resistance to rifampicin and isoniazid would be confirmed by a second test. All TB patients underwent a human immunodeficiency virus (HIV) test . Ethical approval of this secondary analysis was obtained from the Jinhua Municipal Central Hospital (Jinhua, China). Informed consent was waived because the data used in this study were publicly available from a public database , and patient information was anonymous.
The variables of TB patients collected for further analysis were as follows: gender, smoking, alcohol use, coexisting diabetes, HIV infection status, previous TB history, socioeconomic status (divided into three levels based on the scoring system used in the Peruvian National Census) , employment status (unemployed, working, student, or unknown), secondary education status (completed or not), history of incarceration, sputum smear grade, mean cough duration, hospitalization history, side effects of treatment (yes/no), and spoligotypes (stratified based on SpolDB4 database) .
Continuous data and categorical data were expressed as the median (interquartile) and the number (proportion), respectively. Group comparisons between MDR-TB and DS-TB were performed using an unpaired t test or Wilcoxon rank sum test, Pearson Chi square test or the Fisher’s exact test, as appropriate. The least absolute shrinkage and selection operator (LASSO) regression technique was used for data dimension and predictor selection. Multivariable logistic regression analysis was used to develop a predictive model and a nomogram of MDR-TB. The discriminatory capacity of the model was determined by calculating the area under the curve (AUC). The bootstrapping method (resampling = 500) was employed for internal validation . The calibration of the model was evaluated by using the Hosmer–Lemeshow test, and the clinical usefulness of the model was assessed by decision curve analysis (DCA) . Statistical analysis was conducted using R software (version 3.5.1; R Foundation for Statistical Computing, Vienna, Austria), and P < 0.05 was considered statistically significant.
Among the study participants, 30.4% (213/700) were MDR-TB patients. Demographic and clinical characteristics of study participants are shown in Table 1. Of the 14 variables collected from patients, 5 variables were selected based on non-zero coefficients calculated by LASSO regression analysis (Fig. 1). These variables included completed secondary education, previous history of TB, any side effect of treatment, history of hospitalization, and spoligotype.
To develop a predictive model for MDR-TB, multivariable logistic regression analysis was performed based on the aforementioned 5 variables selected by the LASSO regression technique. The AUC for the predictive model was 0.759 (95% confidence interval [CI], 0.719–0.799), and the internal validation using the bootstrap method (resampling = 500) was 0.757 (95% CI 0.715–0.793) (Fig. 2). To present the predictive model, a nomogram was constructed, thereby providing a convenient, personalized tool to predict the probability of MDR-TB (Fig. 3).
The proposed model was well-calibrated (Fig. 4). The Hosmer–Lemeshow test yielded a nonsignificant P value of 0.452, thereby suggesting that there was no statistical departure from a perfect fit between the predicted and observed values.
To assess its clinical usefulness, DCA was also performed. The decision curve showed that when based on the nomogram in this study, the threshold probability of MDR-TB in TB patients was of 70–90% (Fig. 5), and application of this nomogram to predict MDR-TB would add significantly more benefit than either the treat-all scheme or the treat-none scheme.
In the present study, a nomogram for predicting MDR-TB among TB patients was built. This nomogram incorporated 5 variables, including completed secondary education, previous history of TB, any side effect of treatment, history of hospitalization, and spoligotype. The nomogram showed good discriminatory ability, calibration, and clinical usefulness.
DR-TB has always been a public health crisis and a health security threat . Alarmingly, in the last few decades, the number of detected cases of MDR-TB dramatically increased in several countries, including India, China, Indonesia, the Philippines, Pakistan, Nigeria, Bangladesh, and South Africa . Therefore, the possibility of unsuspected drug-resistance should always be considered when assessing a TB patient in any country or region in clinical practice. The WHO estimated that there were 558,000 new cases with resistance to rifampicin, of which approximately 82% had MDR-TB in 2017 . Although DR-TB accounts for only 5% of TB, the costs involving the diagnosis and treatment of DR-TB are approximately 1/3 of that of global TB [1, 3]. Noteworthy, in endemic countries that utilizes traditional regimens for TB patients, those with DR-TB were associated with a staggering mortality, 40% in MDR-TB and 60–70% in extensively DR-TB (XDR-TB) . At present, only 55% of MDR-TB patients are successfully treated worldwide .
Although MDR poses a challenge for treatment, in general, most MDR-TB patients can be cured by early detection of drug-resistance and using appropriately designed treatment regimens. MODS is a traditional method used to test rifampicin or isoniazid resistance in DST, and remains a major method used in some countries and regions. The main drawback of this method is that it is time-consuming and easy to miss the diagnosis. In recent years, the use of phenotypic DST has significantly expanded since the rapid test Xpert MTB/RIF was first recommended by the WHO in 2010 . The advantage of this method is that it simultaneously detects TB and drug resistance to rifampicin in a very short time. Other molecular detection tests include array-based platforms and line probe assays [4, 5]. However, the aforementioned techniques for the identification of MDR-TB relies on sputum samples and are still unavailable in resource-limited settings . Therefore, for MDR-TB patients with no sputum or without drug-resistance testing, or with sputum-negative TB, it is extremely easy to delay treatment or receive inappropriate treatment.
Regarding the risk of TB infection, several factors have been reported, such as smoking, poverty status, coexisting diabetes, HIV positive, and other conditions that lead to immunosuppression [1, 9]. In addition, previous studies have developed several models for predicting TB infection based on nosocomial populations, and it was pointed out by the authors that these models could improve the diagnosis of TB [18,19,20]. However, to our knowledge, few studies have described MDR-TB risk prediction models. In the current study, a nomogram for predicting MDR-TB was built based on 5 variables of TB patients. The variables included in the nomogram were filtered by LASSO regression analysis, which is considered superior to select predictors by univariate analysis [21, 22]. Furthermore, we evaluated the clinical significance of these predictors. Regarding “Education level of the patient”, it has been reported that economic status (poverty status) was associated with TB incidence [1, 9]. In general, the economic status is related to completion of education, and the economic level of those who completed education is higher compared to those who did not completed education. Therefore, we speculate that the “Education level” may associate with MDR-TB infection. Regarding the “History of hospital admission” and “side effects to drug”, it has been reported that poor adherence increases MDR-TB incidence , patients with a history of hospital admission may be better at treating adherence than patients who have never been hospitalized. Chaotic treatment (improper treatment programs) for TB has also resulted in an increase in MDR-TB incidence . TB patients who received more chaotic treatment may experience more side effects to the drug. Hence, “History of hospital admission” and “side effects to drug” may be related to a risk of MDR-TB infection. These 5 predictors are easily available clinically. The nomogram showed good discriminatory ability and calibration, and the DCA evaluation showed its clinical usefulness. Since phenotypic DST has not yet been fully covered in TB-endemic areas, and many TB patients are sputum-free or sputum-negative TB, this cost-free nomogram may be helpful in screening MDR-TB in these TB patients.
This predictive model has several limitations. Firstly, this nomogram was built based on a 3-year prospective study conducted in Peru . There is a regional difference in the prevalence of MDR-TB ; therefore, whether the nomogram is applicable to other regions or countries requires further multicenter verification. Secondly, it cannot be determined whether MDR-TB was acquired using this nomogram. Thirdly, some patient’s variables, such as regimens and duration of treatment were not included in our analysis because they were not available in the original dataset. Despite these limitations, this study was the first to develop a nomogram for predicting MDR-TB risk in TB patients.
In this study, a nomogram for MDR-TB risk prediction was built, which incorporated five demographic and clinical characteristics of TB patients. This may be a useful complement to the current identification of MDR-TB, especially in patients with MDR-TB from whom test samples are not available.
Availability of data and materials
The data used in this study can be downloaded from ‘DATADRYAD’ database (http://www.Datadryad.org).
Least absolute shrinkage and selection operator
Area under the curve
Decision curve analysis
Human immunodeficiency virus
World Health Organization. Global tuberculosis report 2018. New York: World Health Organization; 2018; http://www.who.int/tb/en/. Accessed 31 Mar 2019.
O’Donnell MR, Jarand J, Loveday M, Padayatchi N, Zelnick J, Werner L, Naidoo K, Master I, Osburn G, Kvasnovsky C, Shean K, Pai M, Van der Walt M, Horsburgh CR, Dheda K. High incidence of hospital admissions with multidrug-resistant and extensively drug-resistant tuberculosis among South African health care workers. Ann Intern Med. 2010;153(8):516e22.
Pooran A, Pieterson E, Davids M, Theron G, Dheda K. What is the cost of diagnosis and management of drug resistant tuberculosis in South Africa? PLoS ONE. 2013;8(1):e54587.
Dheda K, Barry CE, Maartens G. Tuberculosis. Lancet. 2016;387(10024):1211e26.
Dheda K, Gumbo T, Gandhi NR, Murray M, Theron G, Udwadia Z, Migliori GB, Warren R. Global control of tuberculosis: from extensively drug-resistant to untreatable tuberculosis. Lancet Respir Med. 2014;2(4):321e38.
Dheda K, Chang KC, Guglielmetti L, Furin J, Schaaf HS, Chesov D, Esmail A, Lange C. Clinical management of adults and children with multidrug-resistant and extensively drug-resistant tuberculosis. Clin Microbiol Infect. 2017;23(3):131–40.
Dheda K, Gumbo T, Maartens G, Dooley KE, McNerney R, Murray M, Furin J, Nardell EA, London L, Lessem E, Theron G, van Helden P, Niemann S, Merker M, Dowdy D, Van Rie A, Siu GK, Pasipanodya JG, Rodrigues C, Clark TG, Sirgel FA, Esmail A, Lin HH, Atre SR, Schaaf HS, Chang KC, Lange C, Nahid P, Udwadia ZF, Horsburgh CR Jr, Churchyard GJ, Menzies D, Hesseling AC, Nuermberger E, McIlleron H, Fennelly KP, Goemaere E, Jaramillo E, Low M, Jara CM, Padayatchi N, Warren RM. The epidemiology, pathogenesis, transmission, diagnosis, and management of multidrug-resistant, extensively drug-resistant, and incurable tuberculosis. Lancet Respir Med. 2017;5:291–360.
Kendall EA, Fofana MO, Dowdy DW. Burden of transmitted multidrug resistance in epidemics of tuberculosis: a transmission modelling analysis. Lancet Respir Med. 2015;3(12):963e72.
Grandjean L, Gilman RH, Martin L, Soto E, Castro B, Lopez S, Coronel J, Castillo E, Alarcon V, Lopez V, San Miguel A, Quispe N, Asencios L, Dye C, Moore DA. Transmission of multidrug-resistant and drug-susceptible tuberculosis within households: a prospective cohort study. PLoS Med. 2015;12(6):e1001843.
Middlebrook G, Cohn ML. Some observations on the pathogenicity of isoniazid-resistant variants of tubercle bacilli. Science. 1953;118(3063):297–9.
Gagneux S. Fitness cost of drug resistance in Mycobacterium tuberculosis. Clin Microbiol Infect. 2009;15(Suppl 1):66–8.
Moore DA, Mendoza D, Gilman RH, Evans CA, Hollm Delgado MG, Guerra J, Caviedes L, Vargas D, Ticona E, Ortiz J, Soto G, Serpa J, Tuberculosis Working Group in Peru. Microscopic observation drug susceptibility assay, a rapid, reliable diagnostic test for multidrug-resistant tuberculosis suitable for use in resource-poor settings. J Clin Microbiol. 2004;42(10):4432–7.
Caviedes L, Lee TS, Gilman RH, Sheen P, Spellman E, Lee EH, Berg DE, Montenegro-James S. Rapid, efficient detection and drug susceptibility testing of Mycobacterium tuberculosis in sputum by microscopic observation of broth cultures. The Tuberculosis Working Group in Peru. J Clin Microbiol. 2000;38(3):1203–8.
Grandjean L, Gilman RH, Martin L, Soto E, Castro B, Lopez S, Coronel J, Castillo E, Alarcon V, Lopez V, San Miguel A, Quispe N, Asencios L, Dye C, Moore DA. Data from: transmission of multidrug-resistant and drug-susceptible tuberculosis within households: a prospective cohort study. Dryad Digital Repository. 2015. https://doi.org/10.5061/dryad.br760.
Collins GS, Reitsma JB, Altman DG, Moons KG. Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): the TRIPOD statement. Ann Intern Med. 2015;162(1):55–63.
Van Calster B, Wynants L, Verbeek JFM, Verbakel JY, Christodoulou E, Vickers AJ, Roobol MJ, Steyerberg EW. Reporting and interpreting decision curve analysis: a guide for investigators. Eur Urol. 2018;74(6):796–804.
Seung KJ, Keshavjee S, Rich ML. Multidrug-resistant tuberculosis and extensively drug-resistant tuberculosis. Cold Spring Harb Perspect Med. 2015;5(9):a017863.
Rakoczy KS, Cohen SH, Nguyen HH. Derivation and validation of a clinical prediction score for isolation of inpatients with suspected pulmonary tuberculosis. Infect Control Hosp Epidemiol. 2008;29(10):927–32.
Huang YQ, Liang CH, He L, He L, Tian J, Liang CS, Chen X, Ma ZL, Liu ZY. Development and validation of a radiomics nomogram for preoperative prediction of lymph node metastasis in colorectal cancer. J Clin Oncol. 2016;34(18):2157–64.
Pegues CF, Johnson DC, Pegues DA, Spencer M, Hopkins CC. Implementation and evaluation of an algorithm for isolation of patients with suspected pulmonary tuberculosis. Infect Control Hosp Epidemiol. 1996;17(7):412–41.
Tibshirani R. The lasso method for variable selection in the Cox model. Stat Med. 1997;16(4):385–95.
Balachandran VP, Gonen M, Smith JJ, DeMatteo RP. Nomograms in oncology: more than meets the eye. Lancet Oncol. 2015;16(4):e173–80.
We appreciate all the participants involved in the study and the data providers of the study .
This study was supported by the Medical and Health Science and Technology Plan Project of Zhejiang Province (2018KY861), the Chinese Medicine Science and Technology project of Jinhua City (2017jzk05), and the Science and Technology Project of Jinhua City (2020XG-15, 2020-4-013).
Ethics approval and consent to participate
This study was approved by the ethics committee of Jinhua Municipal Central Hospital. The requirement for informed consent was waived because the data were anonymous.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wang, S., Tu, J. Nomogram to predict multidrug-resistant tuberculosis. Ann Clin Microbiol Antimicrob 19, 27 (2020). https://doi.org/10.1186/s12941-020-00369-9