Incidence Density and Predictors of Multidrug-Resistant Tuberculosis Among Individuals With Previous Tuberculosis History: A 15-Year Retrospective Cohort Study

Background: To date, too little attention has been paid to monitoring and estimating the risk of incident multidrug-resistant tuberculosis (MDR-TB) among individuals with a previous tuberculosis history (PTBH). The purpose of this study was to assess the incidence of and risk factors for MDR-TB in those individuals. Methods: Between 2005 and 2020, a large, retrospective, population-based cohort study was performed in Hangzhou, China. A multivariable Cox regression model was used to evaluate independent predictors of incident MDR-TB among individuals with PTBH. Results: The incidence density of MDR-TB was 22.6 per 1,000 person-years (95% confidence level and an interval of 20.9–24.3) for individuals with PTBH. The incidence of MDR-TB increased significantly in individuals who • were under 60 years old. • were male. • had a history of direct contact. • came from low-income families. • worked in high-risk occupations. • lived in rural areas. • had a retreatment TB history. • had an unfavorable outcome in their previous treatment (P < 0.05). In addition, we found that the following factors were significantly linked to the MDR-TB risk among individuals with PTBH (P < 0.05): • sociodemographic factors such as the 21–30 and 31–40 year age groups, or a history of direct contact. • clinical factors like passive modes of TB case finding (PMTCF), human immunodeficiency virus infection, unfavorable treatment outcomes, retreated TB history, non-standardized treatment regimens of retreatment TB patients, and duration of pulmonary cavities (DPC). • microbiological factors, such as duration of positive sputum culture. We also found that the 21–30 year age group, low family income, and PMTCF were significantly linked to incident MDR-TB only in males with PTBH, whilst the 41–50 year age group, extended treatment course, and DPC were significantly associated with female MDR-TB only. Conclusion: The incidence of MDR-TB was high, with a higher rate among subjects with a history of direct contact and unfavorable treatment outcomes. There was a gender difference in the incidence density and risk factors of MDR-TB among individuals with PTBH. Long-term monitoring and gender-specific risk-factor modifications should be given to individuals with PTBH.

Background: To date, too little attention has been paid to monitoring and estimating the risk of incident multidrug-resistant tuberculosis (MDR-TB) among individuals with a previous tuberculosis history (PTBH). The purpose of this study was to assess the incidence of and risk factors for MDR-TB in those individuals.
Methods: Between 2005 and 2020, a large, retrospective, population-based cohort study was performed in Hangzhou, China. A multivariable Cox regression model was used to evaluate independent predictors of incident MDR-TB among individuals with PTBH.
Results: The incidence density of MDR-TB was 22.6 per 1,000 person-years (95% confidence level and an interval of 20.9-24.3) for individuals with PTBH. The incidence of MDR-TB increased significantly in individuals who • were under 60 years old.
• had a history of direct contact.
• came from low-income families.
• worked in high-risk occupations.
• lived in rural areas.
• had a retreatment TB history.
• had an unfavorable outcome in their previous treatment (P < 0.05). In addition, we found that the following factors were significantly linked to the MDR-TB risk among individuals with PTBH (P < 0.05): • sociodemographic factors such as the 21-30 and 31-40 year age groups, or a history of direct contact.
• clinical factors like passive modes of TB case finding (PMTCF), human immunodeficiency virus infection, unfavorable treatment outcomes, retreated TB history, non-standardized treatment regimens of retreatment TB patients, and duration of pulmonary cavities (DPC).

INTRODUCTION
Multidrug-resistant tuberculosis (MDR-TB) is associated with severe and fatal respiratory disease in humans. To date, MDR-TB continues to be a serious global public health issue (1). Although the prevalence of MDR-TB in China has increased dramatically in the past decades (2), prevention and early detection of MDR-TB are still inadequate (3).
Continuous monitoring indicates that some MDR-TB patients have a previous history of tuberculosis (PTBH) (e.g., a case after completing TB treatment happens MDR-TB) (4,5). The most recent MDR-TB epidemic revealed an increased risk of morbidity in individuals with PTBH (6,7). In our investigation of MDR-TB, we found that the surveillance and management of MDR-TB were less sensitive (i.e., no strategies for early detection) among individuals with PTBH.
Our literature review concluded that a clear understanding of how previous TB cases contributed to future MDR-TB onset was still lacking. Previous reports focused primarily on the prevalence and predictors of MDR-TB in TB patients (5,(8)(9)(10). Few researchers have studied the incidence and risk factors of MDR-TB among individuals with PTBH (10)(11)(12).
To reduce the morbidity and mortality of MDR-TB, the government and research institutions must address the potential disease burden and multifactor pathogenesis of MDR-TB among individuals with PTBH. We conducted an extensive, retrospective, population-based cohort study to determine the association between PTBH and incident MDR-TB in Hangzhou City, China, between October 1, 2005, and September 30, 2020. The purpose of this study was to (1) assess the incidence density of MDR-TB among individuals with PTBH, and (2) define specific risk factors for MDR-TB in this population.

Study Design and Settings
The study was conducted in Hangzhou, China. The study cohort included individuals with PTBH, including newly diagnosed TB history (NDTH) and retreated TB history (RTH) with drugresistance test (DRT) who were recruited in the study between October 1, 2005, and September 30, 2020. For this study, we selected MDR-TB cases from all TB designated hospitals in Hangzhou City. MTB-TB cases were diagnosed by clinicians through Gene Xpert (GX) and traditional drug susceptibility testing (DST) (13).
Subjects were divided into incident MDR-TB group (i.e., the exposure group) and non-incident MDR-TB group (such as the control group) based on DRT results for participants during the follow-up period. Subjects were selected if they • had a history of TB but did not have MDR-TB during previous treatment episodes; • were surviving during the study; • had a history of TB treatment; • had a confirmed treatment outcome; and • were available for follow-up.

Subjects were excluded if
• they had a history of MDR-TB infection before the previous TB treatment episodes; • no DST results were reported; • they were a TB patient being treated (i.e., a case with an anti-TB drug therapy during the course of study); • no treatment outcome could be identified; FIGURE 1 | Flow chart of the study population. TB, tuberculosis; MDR-TB, multidrug-resistant tuberculosis; NDTH, newly diagnosed tuberculosis history; RTH, retreated tuberculosis history; DST, drug susceptibility testing.
• they were not treated (i.e., a subject who refused treatment after TB diagnosis) (Figure 1).

Sample Size Calculation
To calculate the sample size of cohort study, we used the following formula (14): where N = sample size; α = alpha (expected significant level, two-tailed test); β = 1-power (expected power, two-tailed test); Z statistics (Z) -Z statistics for confidence level; Z α = standard normal variate for level of significance; Z β = standard normal variate for power or type 2 error as explained in earlier section; m is the number of control subjects per experimental subject; p 0 is the probability of event in controls (p 0 can be estimated as the population prevalence of the event under investigation); p 1 is the probability of an event in experimental subjects. In this study, the investigators presented their results with a 95% confidence interval (CI), Z 0.05 = 1.96 (α = 0.05), Z 0.10 = 1.64 (β = 0.10), m = 1, p 0 = 0.0021% (15), and p 1 = 8.0% (16). According to the sample size calculation, our study had to take at least 153 samples in each group. However, subjects were recruited in this study using the "all-comers" principle (17), provided they met the inclusion and exclusion criteria. This approach meant the sample was more likely to be representative of the general population in Hangzhou.

Data Collection
Data was collected from 10 TB hospitals that provided TB treatment and management in Hangzhou. Laboratory test data included sputum smear, culture, and DST during the baseline and follow-up visits.
All data was collected from the NTSS and entered into an electronic database. Standard participant reporting included the collection by trained investigators of sociodemographic, clinical, and microbiological information, along with initial and followup visits.

Variables and Definitions
In the present study, case definitions and classifications were consistent with the World Health Organization's (WHO) revised TB definitions and reporting framework (18). Incident MDR-TB or non-incident MDR-TB was defined as the main outcome variable. Table 1 shows the definitions used in the study. Based on WHO and national guidelines, we defined and classified the primary covariate variables, and sputum smear, culture, and DST results (18)(19)(20).

Incidence Density
The incidence density of MDR-TB per 1,000 person-years (PYs) was calculated as the number of cases who developed new onset of MDR-TB during the follow-up period divided by the PYs. The PYs were evaluated by using an approximate algorithm [i.e., PYs = (the annual average number of individuals with PTBH) × (total observation years); the annual average number of individuals with PTBH referred to the meaning of individuals with PTBH at the beginning of 2 years in a row] (21). The incidence density was presented with 95% CI using the Poisson distribution. The densities were compared by calculating incidence density ratios with a 95% CI.

Laboratory Methods
The present study mainly used the methods of traditional detection (e.g., sputum smear and culture) for diagnosing TB cases from 2005 to 2014 (22). Furthermore, the TB diagnosis mainly included molecular biological (i.e., a GX method) and traditional detection methods between 2015 and 2020 (13). Based on the DST result, an MDR-TB case was diagnosed in the designated laboratory. The DST was performed on all culture-positive isolates against first-line (isoniazid, rifampicin, pyrazinamide, ethambutol, and streptomycin) and second-line anti-TB drugs (such as kanamycin and ofloxacin) (20). The methods included conventional DST and GX Mycobacterium TB (MTB)/rifampicin. According to standard procedures, a solid or automated liquid culture media system (BACTEC MGIT 960; Becton Dickinson, Sparks, Maryland, USA) was used for the DST (20). The GX MTB/rifampicin (Cepheid, Sunnyvale, USA) was a semi-nested polymerase chain reaction system with fully automated real time, based on molecular beacon technology (23).
The DST was conducted during follow-up. The detections of conventional DST and GX MTB/rifampicin of collected samples were performed at the TB Program Laboratory of Hangzhou Center for Disease Control and Prevention (a Biosafety Level 3 laboratory with proficiency testing approved by the National Reference Laboratory in China). We excluded MDR-TB cases with laboratory cross-contamination. Borderline TB drugs resistance was considered resistant.

Statistical Analysis
In this study, a binary variable including incident MDR-TB and non-incident MDR-TB categories was defined as the outcome variable. The analysis of subjects' characteristics was conducted by using the descriptive method. Continuous variables were described using the mean with standard deviation (SD) while categorical data were analyzed using percentage (proportion). Markov chain Monte Carlo (MCMC) method of multiple imputation (MI) was used to impute censored missing data in this study (24). The data was missing at random in this study. The missing-data mechanism was called missing at random according to the probability that a value was missing did not depend on the missing value but did depend on observed quantities (values of variables that were measured). The following variables were used in our MI: family income, comorbidities, treatment course, and duration of negative sputum smear. The MI procedure in R software was used to perform 200 imputations of each variable by using a MCMC approach. We evaluated the initial estimates for MCMC through the expectation-maximization algorithm.
Univariable and multivariable Cox proportional hazards regression models were used to analyze the risk of incident MDR-TB among individuals with PTBH. Subjects who could not be evaluated during the study were excluded from the analysis. A univariable Cox regression analysis was conducted to identify the factors associated with incident MDR-TB. Variables were analyzed using hazard ratio (HR) generated using the same Cox regression analysis.
Subsequently, independent predictors associated with incident MDR-TB were evaluated using HR generated by a multivariable Cox regression model. All variables with a P < 0.05 were included the Cox analysis using the backward stepwise method based on the minimum statistics of the Akaike information criterion. Variables with a P < 0.05 were considered statistically significant in the Cox analysis and included in the final Cox proportional hazards regression model. R software (version i 386 4.0.4; www.R-project.org, 2021) was used for all statistical analyses.

MDR-TB case A patient infected with TB resistant to at least H and R.
Incident MDR-TB An MDR-TB patient is diagnosed between the initiating of previous anti-TB treatment and the date of study end.
Non-incident MDR-TB Incident MDR-TB has not occurred between the initiating of previous anti-TB treatment and the date of study end.

Low-income level
The annual household income is below RMB 150,000 Yuan.

Middle level and above income
The annual household income is greater than or equal to RMB 150,000 Yuan.
High-risk occupation Such as factory workers, unemployed persons, rural migrant workers.
Low-risk occupation Such as peasants, teachers, students, service workers, nursing workers, attendants, doctors, fisherfolk, drivers, office workers, and retirees.
Active modes of TB case finding The suspected or confirmed patients are found according to early clinical symptoms, initiative care-seeking, and high-risk population screening.
Passive modes of TB case finding The suspected or confirmed patients are found according to physical health examination, differential diagnosis of other diseases, and the screening of close contacts.
A history of direct contact There is a history of close contact with MDR-TB cases within the past 12 months before the onset of MDR-TB.
Previously treated TB patients It is defined as patients who have a previous history of treatment with anti-TB drugs for 1 month or more. It includes relapse patients, treatment after failure patients, treatment after loss to follow-up patients, and other previously treated patients.

Retreated TB patients
Refers to patients with initial treatment failure (e.g., NDTPs with sputum positive are still sputum-test-positive at the end of the 5th month or after a course of treatment), relapse cases (e.g., TB has a relapse for the cured NDTPs or NDTPs with completing treatment), returned cases (e.g., re-entry after abandoning therapy), chronic cases, and other (i.e., loss to follow-up, discontinued treatment, and unknown or undocumented therapy outcomes) cases.
TB patients were not treated Refers to refusing treatment after diagnosis.

Favorable treatment outcome
The standard definition is as follows: (1) TB cases of previous sputum positive are cured after completing treatment; (2) TB cases of previous sputum negative have completed TB treatment, who have presented a remarkable improvement of clinical manifestation.

Unfavorable treatment outcome
Here it refers to previous TB patients who are not cured or have not a remarkable improvement of clinical manifestation after completing treatment, or have not completed TB therapy (i.e., TB cases with discontinued treatment).
Individuals with PTBH It is defined as follows: (a) subjects have a history of anti-TB treatment (such as a history of newly diagnosed TB treatment or a history of retreated TB treatment) before their inclusion in this study; (b) after a period of anti-TB treatment, a patient is represented as a favorable treatment outcome or an unfavorable treatment outcome.
Different PTBH Including newly diagnosed TB history and retreated TB history.

2HRZE/4HR
NDTPs are treated by using first-line drug treatment (i.e., R, H, E, and Z) during the intensive treatment phase of 2 months, and using R and H during a 4-month consolidation period; the frequency of TB therapy is once a day.
2H3R3Z3/4H3R3 NDTPs are treated by using first-line drug treatment (i.e., R, H, and Z) during the intensive treatment phase of 2 months, and using R and H during a 4-month consolidation period; the frequency of TB therapy is once every 3 days.
2H3R3Z3E3/4H3R3 NDTPs are treated by using first-line drug treatment (i.e., R, H, E, and Z) during the intensive treatment phase of 2 months, and using R and H during a 4-month consolidation period; the frequency of TB therapy is once every 3 days.
2HREZ/4H3R3 NDTPs are treated by using first-line drug treatment (i.e., R, H, E, and Z) during the intensive treatment phase of 2 months, and using R and H during a 4-month consolidation period; the treatment frequency of NDTPs is once a day and once every 3 days during the intensive and consolidation treatment phases, respectively.

2HRZES/6HRE
RTPs are treated by using TRs (i.e., R, H, E, Z, and S) during the intensive treatment phase of 2 months, and using H, R, and E during a 6-month consolidation period; the frequency of TB therapy is once a day.

3HRZE/6HRE
RTPs are treated by using TRs (i.e., R, H, E, and Z) during the intensive treatment phase of 3 months, and using H, R, and E during a 6-month consolidation period; the frequency of TB therapy is once a day.

3HRZES/6HRE
RTPs are treated by using TRs (i.e., R, H, E, Z, and S) during the intensive treatment phase of 3 months, and using H, R, and E during a 6-month consolidation period; the frequency of TB therapy is once a day.
2H3R3Z3E3S3/6H3R3E3 RTPs are treated by using TRs (i.e., R, H, E, Z and S) during the intensive treatment phase of 2 months, and using H, R, and E during a 6-month consolidation period; the frequency of TB therapy is once every 3 days.

Individualized TRs
According to the clinical experience of doctors, it is defined as the TRs of 4-6 anti-TB drugs (such as 3 or 4 first-line drugs treatment, an injectable second-line drug, and/or a fluoroquinolone).
Extended treatment course Refers to greater than a 6-month course of treatment for NDTPs or greater than a 9-month course of treatment for RTPs.
Standardized treatment course Refers to a 6-month course of NDTPs treatment and the 8-month treatment course or the 9-month treatment course for RTPs.  subjects is shown in Figure 1. The basic characteristics of the subjects are listed in Table 2.

MDR-TB
Of the 13,551 subjects with PTBH, including 12,172 subjects who completed the follow-up study, 89 deaths, the loss to followup of 1,004 subjects, and 286 subjects with missing data, were included in the analytic cohort. Of the total, the vast majority was of the Han nationality [13,344 (98.47%)]. Most of the subjects had educational levels of high school and below [9,872 (72.85%)] and low family incomes [8,827 (65.14%)]. The mean age in the study population was 50.96 ± 21.03, and the ratio of males to females was 2.25-1. The main reason for this was that TB occurred mainly in the elderly population or males in China (25).

Incidence Density of MDR-TB Among Individuals With PTBH
The subjects were observed for a period of 15 PYs in total. Among the 13,551 individuals without MDR-TB at the baseline, there were 673 new cases of MDR-TB. As shown in Table 3, the incidence density of MDR-TB was 22.6 per 1,000 PYs (95% CI 20.9-24.3) among individuals with PTBH in Hangzhou.
We then examined the incidence density of MDR-TB based on age, gender, nationality, education levels, registered household, family income, occupational risk, residences, history of direct contact, different PTBH, types of MDR-TB diagnosis, and treatment outcomes. Among them, subjects with a history of direct contact had the highest incidence of MDR-TB [57.5 per 1,000 PYs (95% CI 47.2-69.4)]. In contrast, subjects with favorable treatment outcomes had the lowest incidence [8.6 per 1,000 PYs (95% CI 7.6-9.7)]. The incidence of MDR-TB increased significantly for subjects under 60 years of age [crude relative risk 2.38 (95% CI 1.88-3.00)]. It also increased for those with a history of direct contact, male, low family income, high-risk occupation, rural areas, repeated treatment for TB, and unfavorable treatment outcomes (P < 0.05) ( Table 3).  covariate and the MDR-TB risk in subjects. Twenty-one of the 38 tested covariates were associated with a high risk of incident MDR-TB in the study population (P < 0.05). The significant covariates were:

Predictors of Incident MDR-TB Among Individuals With PTBH
• sociodemographic characteristics, including the 21-30, 31-40, and 41-50 year age groups, male, a history of direct contact, low family income (LFI), high-risk occupation, high school education and below, residence in rural areas, and migrant individuals with PTBH.   The remaining 17 covariates were not associated with incident MDR-TB among individuals with PTBH (P > 0.05): • three age groups: 0-10, 11-20, and 51-60 years.
• a history of direct contact (e.g., unknown).
• TB patients with severe infection.
• family members' management or self-management for TB patients. • non-standardized TRs for newly diagnosed TB patients.
• frequencies of chest X-ray examination.
• duration of abnormal imaging findings.
• duration without chest radiological findings.
• duration of negative sputum culture.
• duration of negative sputum smear.
To further explore the independent predictors of incident MDR-TB in individuals with PTBH, we conducted a multivariable Cox regression analysis. Table 5 summarizes the results of the analysis for this population. The Cox analysis demonstrated a significant MDR-TB risk in subjects with PTBH and the following characteristics: • two age groups: 21-30 and 31-40 years.
• a history of direct contact.
From this model, we could also see that a history of direct contact was the most striking predictor for incident MDR-TB in this population (OR 5.04, 95% CI: 1.53-10.22, P < 0.001).
To explore possible gender differences in the association between incident MDR-TB and predictors, we performed genderstratified, multiple regression analyses. Analysis stratified by gender showed that the 21-30 year age group, LFI, and PMTCF were significantly linked to incident MDR-TB only in males, whilst the 41-50-year age group, ETC, and DPC were significantly associated with female MDR-TB only. Seven variables (i.e., the 31-40 year age group, a history of direct contact, HIV infection, RTH, unfavorable treatment outcome, non-standardized TRs of RTPs, and DPSC) were related to incident MDR-TB in both genders ( Table 5). The models also indicated that of all of the independent predictors, a history of direct contact was the strongest impact factor for both male MDR-TB (OR 5.42, 95% CI: 1.98-9.97, P = 0.002) and female MDR-TB (OR 3.63, 95% CI: 1.67-7.03, P < 0.001) ( Table 5).

DISCUSSION
In this study, we conducted a 15-year retrospective cohort study to explore the incidence and risk factors of MDR-TB in individuals with PTBH. Our findings may provide more reliable evidence in developing prevention and control strategies for MDR-TB among individuals with PTBH. Because the pathogenic mechanism of MDR-TB has not been fully clarified, we can only decrease the risk of incident MDR-TB by modifying potential risk factors. Thus, we anticipate these results will be useful in reducing the disease burden of MDR-TB and improving risk monitoring and management of MDR-TB among individuals with PTBH.

Incidence Density of MDR-TB
The incidence density of MDR-TB reported in this study is high, with 22.6 per 1,000 PYs from individuals with PTBH in Hangzhou, China. This finding broadly supports the work of other studies in this area linking MDR-TB with previously treated patients (26). One implication of this finding is the possibility that the government and TB control officials should immediately take measures to monitor and manage individuals with PTBH. As demonstrated in our study, the incidence of MDR-TB increased significantly in subjects under 60 years old and characterized by a history of direct contact, male, low family income, high-risk occupation, rural area of residence, RTH, and unfavorable treatment outcomes. Up to now, China has experienced significant challenges when facing a high incidence of MDR-TB (1). One possible explanation for this may be that the high-risk MDR-TB population is often not well-monitored and managed. In addition, because the TB prevalence level stays high in China, which can, in part, explain the high incidence of MDR-TB (27). To reduce the burden of MDR-TB, interventions should be initiated for the subjects in those high-risk groups.
In our study, the highest incidence density was found in close contact with the MDR-TB case among individuals with PTBH. A possible explanation for this might be that most subjects did not use personal protective measures when they closely contacted the MDR-TB case (28). That strong association implies that interventions like early isolation and treatment of MDR-TB patients, personal protective measures for susceptible persons, and early detection of close contacts are urgently needed to reduce MDR-TB incidence among individuals with PTBH. It is noteworthy that the incidence of MDR-TB was higher amongst individuals with unfavorable treatment outcomes in their previous treatment episodes. In the current study, comparing NDTH with RTH showed that the incidence of MDR-TB was higher in individuals with RTH. A similar trend was observed in previous studies that evaluated the prevalence of MDR-TB in TB patients (6,10,29). Enhanced attention and long-term monitoring should be given to patients who have been treated previously, especially those with a history of direct contact, unfavorable treatment outcomes or RTH.

Predictors of Incident MDR-TB
To our best knowledge, this is the first study to examine the associations between MDR-TB and a broad range of potential risk factors (i.e., epidemiological, clinical, and bacteriological factors) among a large sample of the Chinese population with PTBH. A notable finding in this study was that 10 independent predictors are associated with the increased risk of MDR-TB among individuals with PTBH. These findings were also reported in many previous studies of TB patients (6,10,29).
Notably, this study also observed that RTPs are treated using non-standardized TRs (i.e., 2H3R3Z3E3S3/6H3R3E3, 3HRZE/6HRE, and individualized TRs), which increased the risk of incident MDR-TB dramatically. Thus, it can be seen that the acquired infection due to poor TRs during TB treatment is one of the main risk factors for MDR-TB. This association might be attributed to the increased chance of drug resistance with longer exposure to anti-TB drugs (30). Furthermore, the treatment of non-standardized TRs will likely lead to an increased MDR-TB risk by using substandard anti-TB drugs (i.e., poor drug dosages) in TB cases (31). To decrease the risk of MDR-TB, standardized TRs must be implemented by RTPs. However, in reviewing the literature, we found that there were two distinct standpoints for non-standardized TRs. On the one hand, from a programmatic perspective, one would like to see standardized treatment regimens (32). On the other hand, from a personalized medicine perspective, if proper DST or sequencing is available, personalized regimens might be indicated for specific patients (33). Thus, the selection of TRs for TB patients should be based on scientific and precise assessment. Surprisingly, the present study found a relationship between the PMTCF (such as physical examination, contact examination, and differential diagnosis of other diseases) and incident MDR-TB. A possible explanation for this association is that the increasing risk of MDR-TB may originate from TB cases who were not diagnosed due to a failure of diagnostics at baseline diagnosis. Further analysis reveals that the main reason for this is that there is inadequate sputum smear and culture, the lack of DST, low sensitivity for detecting TB, and poor patient compliance. This reminds us that adequate medical support (i.e., sputum smear, sputum culture, and DST) should be considered by the government. The detection sensitivity for TB should be strengthened in the medical institution. In addition, a health communication schedule with the knowledge of TB diagnosis and treatment should be provided to increase patient compliance. Another possible reason is that delayed diagnosis and treatment of TB may increase the risk of MDR-TB (34). As far as we know, PMTCF may lead to the delayed finding of TB cases. If the TB case could not be found timely, the diagnosis and treatment of TB would be delayed. Maybe the TB case would develop into serious TB leading to an extended course of treatment, which might become a risk factor associated with MDR-TB (35). This study provides important evidence that active modes of TB case finding are beneficial for preventing or reducing MDR-TB morbidity.
Most importantly, we used a stratified comparison to analyze gender differences in the predictors of MDR-TB that were found in individuals with PTBH. There were many differences in predictors for incident MDR-TB between the sexes. These findings might be attributed to certain gender disparities with sociodemographic, therapeutic and managed factors implicated in the development of MDR-TB (32,36). In fact, the 21-30 year age group, LFI, and PMTCF in males with PTBH are associated with incident MDR-TB. Several possible reasons for this increased risk are as follows. First, it may be due to the high-intensity of work and study leading to decreased immunity. Second, TB cannot obtain timely diagnosis, therapy, and management because the males usually need to go out to work to support their family (i.e., the floating population). We also found that the 41-50 year age group, ETC, and DPC contributed to female MDR-TB only. Although the underlying mechanisms are unclear, a possible explanation for these results may be the lack of adequate nutrition due to reducing body weight in females with PTBH (37). Women might have a low body mass index (BMI) because of poor nutrition and therefore develop MDR-TB. In addition, they might develop MDR-TB due to anxiety for the BMI and poor health (38). If the anxiety would be further developed, the risk of anorexic might be greatly increased. Subsequently, it would promote the occurrence of MDR-TB. Similarly, these factors might lead to ETC and longer DPC for females in their previous treatment episodes (39). Interestingly, the current study found that the risk of MDR-TB had the difference of age groups between genders. It could therefore be assumed that a health promotion program should be performed in males with the 21-40 years age group and females with the 31-50 years age group for the control and prevention of MDR-TB. In summary, gender-specific intervention programs against MDR-TB among individuals with PTBH should be considered to reduce the MDR-TB risk by modifying risk factors.
Finally, it should be reminded that the present study indicated that the risk of MDR-TB among individuals with PTBH was attributed mainly to acquired infection via exposure to a case of MDR in the household or community. The present finding showed that there was an urgent need to contain the epidemic of MDR-TB through potential intervention strategies (such as early detection, early isolation, early diagnosis, and early anti-TB treatment of MDR-TB case, and personal protective measures of the susceptible population) in individuals with PTBH. Moreover, the government would need to carry out a health education program on the knowledge of TB infection control in the household or community. Although these findings were also reported in many previous studies of TB patients (5)(6)(7)(8)(9)(10)(11). For example, unfavorable treatment outcomes remain a primary concern for the control of MDR-TB among individuals with PTBH. To monitor and manage the risk of MDR-TB, the focus must be on high-risk factors in individuals with PTBH.
Our studies had several limitations. First, it was a retrospective study, which might cause recall bias. To reduce recall bias, we collected the first follow-up records [i.e., epidemiological history including clear contact history and unclear contact history, demographic data (such as age, gender, occupation, and areas), and clinical data] after the diagnosis of MDR-TB through the TB information system. Secondly, we failed to include some longestablished impact factors (e.g., the frequency or intensity of exposure). Thirdly, there are some potential confounders like the psychological health of TB case, the stability, malabsorption, and quality of TB drug could not be controlled. Fourthly, our study might have a selection bias. For example, we might not include all MDR-TB cases, due to non-access or access at other, non-hospital facilities. To reduce selection bias, we have retrieved and collected medical records of MDR-TB cases from hospitals outside Hangzhou through the NTSS. However, we might not include MDR-TB cases with non-access or nonhospital facilities. According to our investigation, there was little for patients with non-access or non-hospital facilities in China. Fifth, the present study has the limitations of using electronic data (i.e., non-standardized, especially over a long period). Finally, in order to reduce limitations, we performed additional measures such as the multicenter research, follow-up of the participants, and the inclusion of additional predictors. Moreover, we verified the validity and reliability of electronic data by using the face-to-face interview among 100 participants with random selection. Despite several limitations of this study, there are limited reports on impact factors of MDR-TB among individuals with PTBH. We feel these findings would instill some cognition on the novel strategy of prevention and control of incident MDR-TB.

CONCLUSION
In conclusion, our study showed that, in individuals with PTBH, the incidence of MDR-TB is high, with higher rates among subjects with a history of direct contact and unfavorable treatment outcome. Enhanced attention and long-term monitoring should be given to patients who have been treated previously, especially those with a history of direct contact or previous unfavorable treatment outcomes. This study broadens our knowledge of MDR-TB as a growing public health issue in China and underscores the necessity of healthcare plans for the prevention of MDR-TB.
We identified four primary risk factors for MDR-TB among individuals with PTBH. They included a history of direct contact, unfavorable treatment outcomes, non-standardized TRs for RTPs, and retreated TB history. We also found that there was a gender difference in risk factors of MDR-TB among individuals with PTBH. Therefore, we offer the following recommendations: • provide a personal protection program of close contacts for MDR-TB. • increase the rate of favorable treatment outcomes in TB cases.
• strengthen early detection of TB.
• modify gender-specific risk factors.
Such measures are critical to preventing the spread of infection from MDR-TB in individuals with PTBH.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by the Ethics Committee of the Hangzhou Center for Disease Control and Prevention. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.