Skip to main content

Regional distribution of Mycobacterium tuberculosis infection and resistance to rifampicin and isoniazid as determined by high-resolution melt analysis



Identifying the transmission mode and resistance mechanism of Mycobacterium tuberculosis (MTB) is key to prevent disease transmission. However, there is a lack of regional data. Therefore, the aim of this study was to identify risk factors associated with the transmission of MTB and regional patterns of resistance to isoniazid (INH) and rifampicin (RFP), as well as the prevalence of multidrug-resistant tuberculosis (MDR-TB).


High-resolution melt (HRM) analysis was conducted using sputum, alveolar lavage fluid, and pleural fluid samples collected from 17,515 patients with suspected or confirmed MTB infection in the downtown area and nine counties of Luoyang City from 2019 to 2021.


Of the 17,515 patients, 82.6% resided in rural areas, and 96.0% appeared for an initial screening. The HRM positivity rate was 16.8%, with a higher rate in males than females (18.0% vs. 14.1%, p < 0.001). As expected, a positive sputum smear was correlated with a positive result for HRM analysis. By age, the highest rates of MTB infection occurred in males (22.9%) aged 26–30 years and females (28.1%) aged 21–25. The rates of resistance to RFP and INH and the incidence of MDR were higher in males than females (20.5% vs. 16.1%, p < 0.001, 15.9% vs. 12.0%, p < 0.001 and 12.9% vs. 10.2%, p < 0.001, respectively). The HRM positivity rate was much higher in previously treated patients than those newly diagnosed for MTB infection. Notably, males at the initial screening had significantly higher rates of HRM positive, INH resistance, RFP resistance, and MDR-TB than females (all, p < 0.05), but not those previously treated for MTB infection. The HRM positivity and drug resistance rates were much higher in the urban vs. rural population. By multivariate analyses, previous treatment, age < 51 years, residing in an urban area, and male sex were significantly and positively associated with drug resistance after adjusting for smear results and year of testing.


Males were at higher risks for MTB infection and drug resistance, while a younger age was associated with MTB infection, resistance to INH and RFP, and MDR-TB. Further comprehensive monitoring of resistance patterns is needed to control the spread of MTB infection and manage drug resistance locally.

Peer Review reports


The incidence of Mycobacterium tuberculosis (MTB) infection is reportedly decreasing, although the actual number of cases could be far greater due to under-reporting and misdiagnosis [1]. Despite the long history of drug development, MTB infection continues as a public health concern. Notably, multidrug-resistant tuberculosis (MDR-TB) and rifampicin (RFP)-resistant TB are particularly tenacious due to the toxicity of anti-TB drugs, high mortality rate, and substantial economic burden [2]. The World Health Organization estimated that there were approximately 465,000 cases of drug-resistant tuberculosis (DR-TB) worldwide in 2019 and only about one-third of MDR-TB cases received appropriate treatment [1, 3]. Inappropriate treatment leads to the development of drug resistance and poor outcomes [4, 5]. Therefore, clarification of regional distributions and risk factors associated with the development of drug resistance could improve utilization of limited healthcare resources to effectively control the spread of MTB infection. Although phenotypic antibiotic susceptibility testing remains the current gold standard for diagnosis of DR-TB [6], this method is time-consuming due to the relatively long detection cycle. Therefore, molecular diagnostic techniques should be more widely applied in clinical practice for detection and resistance analysis for rapid and effective control of MTB infection. High-resolution melt (HRM) analysis [7] can be used to simultaneously analyze sensitivity to first- and second-line drugs as supplemental epidemiological data of regional MTB infection and drug resistance.

Drug resistance, pathogenic toxicity, inappropriate treatment, and rapid urbanization are all associated with the continued spread and treatment outcome of MTB infection [8,9,10]. Since transmission and resistance patterns vary geographically [11], incomplete surveillance data can hinder timely treatment of high-risk groups and assessment of drug resistance [12]. Although treatment history is considered an important risk factor, increasing numbers of patients are initially diagnosed with DR-TB, suggesting the need to increase monitoring of drug resistance [13, 14]. The incidence of MTB infection is regularly reported in many regions, complementing global surveillance and monitoring databases. However, data regarding transmission characteristics, molecular resistance patterns, and risk factors for DR-TB in in this region remain limited. Therefore, the aim of this study was to identify risk factors related to MTB transmission and drug resistance to provide theoretical support for local management of MTB infection in China.

Materials and methods

From January 2019 to December 2021, HRM analysis was performed using sputum, pleural fluid, and bronchoalveolar lavage fluid samples collected from all patients with clinically confirmed or suspected MTB infection. The patients included those with clinical symptoms or latent infection, immunopathological features of MTB infection, imaging signs of infection, close contacts with excreting patients, strong positive tuberculin test reactions and history of MTB infection residing in Luoyang City and county districts (Yanshi, Luanchuan, Songxian, Ruyang, Luoning, Mengjin, Xin’an, Yiyang, and Yichuan).

Prior to analysis, duplicate data were excluded by name, sex, age, and region. Individuals aged < 15 and > 86 years were combined into a single age category, while all others were classified into 5-year age groups. Variables included the date of testing, HRM-positive results, INH resistance (INH-R), RFP resistance (RFP-R), multidrug resistant (MDR), sex, age, region (county-districts are uniformly classified as township areas), sputum smear results, and systematic recorded or self-reported treatment history. MDR-TB was defined as resistance to at least INH and RIF. The patients were classified according to the history of treatment. Newly diagnosed patients, those who have never received anti-MTB treatment, those undergoing standard anti-MTB treatment regimens and have not completed the course of treatment, those receiving irregular anti-MTB therapy for less than 1 month. Retreatment patients, those who received anti-MTB drugs in the past for one month or more, those failed or relapsed in the initial treatment [15, 16]. Initial resistance refers to a newly diagnosed case of DR-MTB. Acquired resistance refers to resistance in previously treated patients who received first-line drugs for at least 4 weeks. Patients newly diagnosed with MDR, as confirmed by HRM analysis, received drug-resistant or MDR-TB management regimen and were not included in the previously treated group. In this study, resistance to RFP was confirmed by the presence of the 81-bp core region (rifampicin-resistance-determining region) of the rpoB gene (codons 426–452 of MTB correspond to codons 507–533 of Escherichia coli) [17] and resistance to INH was confirmed by presence of the codon katG315, inhA promoter region (locus − 17 ~ -8 locus), AhpC promoter region (locus − 44 ~ -30 and − 15 ~ 3), and codon inhA94.

Statistical analysis

The HRM analysis results, HRM-positive results, and rates of RFP-R, INH-R, and MDR-TB were stratified by sex for the newly diagnosed and previously treated groups, urban and rural groups, and sputum-positive and -negative groups. The HRM-positive groups, as well as the RFP-R, INH-R, and MDR-TB groups were then stratified by age to confirm the time pattern of DR-MTB infection. Meanwhile, the prevalence of specific mutations related to drug resistance was also assessed. Multivariate logistic regression models were created with all variables that may be associated with RFP-R, INH-R, and MDR-TB to assess potential correlations with other covariates. The odds ratio (OR) and 95% confidence interval (CI) were calculated. Pearson’s chi-square analysis was performed to identify significant differences in the detection and resistance rates of MTB under different conditions. A probability (p) value of < 0.05 was considered statistically significant. All analyses were performed using STATA/SE 15.1 software (StataCorp LLC, College Station, TX, USA).


The frequency of HRM analysis increased each year during the study period, from 3662 cases (2645 males, 1017 females) in 2019 to 7893 cases (5249 males, 2599 females) in 2021. From 2019 to 2021, HRM analysis was conducted with samples from 17,515 patients, which included 12,054 (68.8%) males and 5461 (31.2%) females. More than 95.0% were newly diagnosed patients, 82.6% resided in rural area, and only 12.3% of the sputum smears were positive for acid-fast bacilli. Of the 17,515 patients included for initial HRM analysis, 16.8% had positive results, with significantly more males than females (18.0% vs. 14.1%, p < 0.001), of which 19.3% were resistant to INH, 14.9% were resistant to RFP, and the incidence of MDR-TB was 12.2%. The HRM positivity rate was significantly higher in previously treated patients than newly diagnosed patients (51.5% vs. 15.3%, p < 0.001). Among the newly diagnosed patients, the HRM positivity rate was significantly higher in males than females (16.9% vs. 11.9%, p < 0.001), whereas among the previously treated patients, the HRM positivity rate was significantly higher in females than males (63.8% vs. 45.6%, p < 0.001). Moreover, the HRM positivity rate was significantly higher in patients residing in urban rather than rural areas (35.3% vs. 12.9%, p < 0.001), with significantly higher rates in males than females in both urban and rural areas (36.9% vs. 31.5%, p = 0.004; 13.9% vs. 10.6%, p < 0.001). The HRM positivity rate was significantly higher in patients with positive vs. negative sputum smear results (76.0% vs. 7.9%, p < 0.001), although there was no significant difference in the detection rate between sexes (Table 1).

Both males and females aged 51–70 years were more likely to undergo HRM analysis. HRM-positive males mainly appeared in 51–65 years old, accounting for 32.5%. In contrast, HRM-positive females mainly concentrated in 21–35 years old, with a proportion of 31.2% (Fig. 1 A and 1B, respectively). Among males, the HRM positivity rate peaked at 26–30 years to 22.9%, then basically stabilized, but remained above 15% in all age categories. In contrast, among females, the HRM positivity rate peaked at 21–25 years to 28.1%, then decreased rapidly to 13.4% at 36–40 years. The median (interquartile range, IQR) age of males who underwent HRM testing and had positive results were 55 (39–67) and 54 (37–65) years, respectively, while in females were 54 (IQR 33, 66) and 48 years (IQR 29, 67), respectively (data not shown).

Fig. 1
figure 1

Number of HRM tests and positive results stratified by age category for males (A) and females (B). The annual count of HRM test results for each age classification is expressed in bars.The proportion of participating HRM tests and HRM-positive results , as well as the proportion of smears positive for each age category are indicated by lines (second y-axis)

Fig. 2
figure 2

MTB positive results and the proportion of RFP-R(A), INH-R(B) and MDR(C) by age and gender category. The total number of MTB positive by age are expressed in bars, and the proportion for RFP-R, INH-R and MDR are expressed in lines (second y-axis). RFP-R rifampicin-resistant, INH-R isoniazid-resistant

During study period, the number of HRM-positive cases increased from 2019 to 2020 and to 2021, although the HRM positivity rate decreased in both males (24.6%, 16.5%, and 15.9, respectively) and females (21.1%, 12.8%, and 12.2%, respectively). Similarly, although the number of INH-R, RFP-R, and MDR-TB cases had increased from 2019 to 2021, the overall incidence of drug resistance tended to decrease in both males and females.

Table 1 Population characteristics of the study stratified by sex

When stratified by age, the highest rates of RFP-R, INH-R, and MDR occurred in males aged 36–40 years and females aged 46–50 years. Among males, 25.9% were resistant to RFP, 25.8% were resistant to INH, and 21.3% were diagnosed with MDR-TB. Among females, 19.6% were resistant to RFP, 26.1% were resistant to INH, and 18.2% were diagnosed with MDR-TB. As shown in Fig. 2 A, the RFP-R rate was lower in males than females at the ages of 41–45 and 71–75 years, but higher than in females in all other age groups. To identify factors responsible for this difference, analysis of the characteristics of HRM-positive populations in different age groups found that the proportions of urban-dwelling males aged 41–45 and 71–75 years (34.4% and 26.6%, respectively) were lower than the mean proportion of the urban-dwelling male population in each age range (36.9%), while the proportions of rural-dwelling males (65.6% and 73.4%, respectively) were higher than the mean proportion of each age range (63.1%). Meanwhile, the proportions of previously treated females aged 41–45 and 71–75 years (18.2% and 24.6%, respectively) were higher than the mean proportions of previously treated females at each age range (16.7%), while the proportions of newly diagnosed females (81.8% and 75.4%, respectively) were lower than the mean proportions of newly diagnosed females at each age range (83.3%). The detection rates of MDR-TB and INH-R among all age groups showed a trend of drug resistance consistent with RFP-R (Fig. 2B C). The data for drug resistance in patients aged < 15 and > 85 years in this study were limited and, thus, will be included in future studies

The rates of RFP-R, INH-R, and MDR-TB were 18.9%, 14.0% and 11.4%, respectively, among newly diagnosed patients, 22.7%, 20.8%, and 17.5% among previously treated patients, 23.2%, 18.9%, and 15.5% among patients residing in urban areas, 17.1%, 12.5%, and 9.4% among patients residing in rural areas, 20.4%, 16.3%, and 13.0% among patients with positive sputum results, and 18.1%, 12.4%, and 10.6% among patients with negative sputum results. RFP-R, INH-R, and MDR-TB were related to sex, age, previous treatment of TB, and region (Table 2)

Table 2 Resistance of RFP, INH and MDR under different factors by sex

Scatter plots were created with age plotted along the y-axis and different variables along the x-axis (Fig. 3). As indicated by the scatter plots presented in Fig. 3 A, there were significant differences between males and females in the incidences of INH-R (20.5% vs. 16.1%, χ2 = 6.827, p = 0.009), RFP-R (15.9% vs. 12.0%, χ2 = 6.814, p = 0.009), and MDR (12.9% vs. 10.2%, χ2 = 3.953, p = 0.047). As shown in Fig. 3B, there were significant differences in the detection rates of RFP-R and MDR-TB between the previously treated and newly diagnosed groups (all, p < 0.001), but not INH-R (p = 0.080). Hence, Pearson’s chi-square analysis was also performed to assess difference in the rates of RFP-R, INH-R, and MDR-TB between males and females in consideration of previous treatment or not for MTB infection. The results revealed statistically significant differences in the rates of INH-R, RFP-R, and MDR-TB between newly diagnosed males and females, but not those previously treated for MTB infection, as shown in Fig. 3 C. Also, the differences between the detection rates of INH-R, RFP-R and MDR-TB in the retreatment females and those in newly diagnosed were statistically significant (p = 0.007, p < 0.001 and p < 0.001, respectively), whereas only RFP-R was significantly different between males previously treated and newly diagnosed (p = 0.029). The results of Pearson’s chi-square analysis revealed significant differences in the rates of RFP-R, INH-R, and MDR-TB between patients residing in urban vs. rural areas (all, p < 0.001). There were also significant differences in resistance rates between males and females residing in urban (p = 0.039, 0.036, and 0.036, respectively), but not rural areas (p = 0.092, 0.095, and 0.418, respectively) (Fig. 3D and E, respectively). As indicated in Fig. 3 F, only RFP-R was statistically significant between positive and negative sputum smears (p = 0.004). Also, there were significant differences in the rates of INH-R, RFP-R, and MDR-TB between males and females with positive sputum results (p = 0.010, 0.005, and 0.011, respectively). As shown in Fig. 3G, there was a significant difference in the rate of RFP-R between males with positive sputum results and females with negative sputum results (p = 0.016). Analysis of drug resistance rates among different age groups showed that the incidences of RFP-R and MDR-TB, but not INH-R, were significantly higher in patients aged < 51 vs. >50 years (p < 0.001, p < 0.001, and p = 0.125, respectively), while the difference in INH resistance rate can be reflected between males aged < 51 years and females aged > 50 year (p = 0.003), and between males and females aged > 50 years (p = 0.042), (Fig. 3 H and 3I, respectively).

Fig. 3
figure 3

Represented correlations among patient age, treatment history, region, smear, and gender by scatter plots. The dots indicate the age of patients of DR-TB and the median age is indicated below each variable. A B D F H indicate differences in resistance to INH, RFP and MDR-TB in the context of different gender, treatment history, regional distribution, sputum smear results, age groups greater than 50 years or less than 51 years, respectively. C E G I indicate differences in detection to RFP-R, INH-R and MDR-TB by sex in the context of different TB treatment history, region, sputum smear results, less than 51 or older than 50 years

Multivariate logistic regression analysis was used to determine whether sex, age, previous TB treatment, region, sputum smear results, and year of testing were associated with the incidence of INH-R, RFP-R, and MDR-TB. The results obtained with an adjusted model revealed that male, age < 51 years, urban area, and previous treatment for TB were positively associated with an increased risk of MTB resistance (Fig. 4). After taking into account the sputum results and year of detection, the risks of RFP-R, INH-R, and MDR-TB were increased in males by 1.47-fold (95% CI = 1.15–1.91), 1.37-fold (95% CI = 1.10–1.73), and 1.56-fold (95% CI = 1.24–1.96), respectively, as compared to females, and were increased by 1.79-fold (95% CI = 1.33–2.37), 1.35-fold (95% CI = 1.03–1.77), and 1.79-fold (95% CI = 1.31–2.42) in patients with a history of MTB infection as compared to those newly diagnosed. In addition, risks of drug resistance were increased by 1.52-fold (95% CI = 1.22–1.89), 1.47-fold (95% CI = 1.21–1.78), and 1.52-fold (95% CI = 1.20–1.93) in patients residing in urban vs. rural areas, and by 1.49-fold (95% CI = 1.20–1.83), 1.14-fold (95% CI = 0.95–1.38), and 1.56-fold (95% CI = 1.24–1.96) in patients aged < 51 vs. >50 years (Table 3). See End of Text.

Fig. 4
figure 4

Graphic expression of the final model of multivariate logistic regression for factors with RFP-R, INH-R, MDR in cases of HRM-positive tuberculosis. (2941 cases of HRM positive from 2019 to 2021)

Table 3 Final model of multivariate logistic regression for RFP-R, INH-R and MDR-TB associated with sex, age, history of previous treatment, region, smear and year of HRM assay. Total HRM positive n = 2,941

Among the gene mutations associated with RFP-R, 91.5% of cases had one mutation, 7.6% had two, and 1.1% had three. The sites with the highest mutation rates of the rpoB gene were codons 529–533 and 521–528, which accounted for 56.7% and 19.7% of all mutations, respectively. The most common double mutation occurred at codons 521–528 and 513–520 followed by codons 521–528 and 529–533, which accounted for 36.4% and 24.2%, respectively. Meanwhile, the most common triple mutations involved codons 507–512, 521–528, and 513–520 (three cases), followed by codons 507–512, 521–528, and 529–533 (two cases). Among the mutations associated with INH-R, single and double mutations were identified in 95.4% and 4.6% of patients, respectively. The most common point mutations involved katG315 and the inhA promoter region (-17 ~ -8 site), which accounted for 88.5% and 16.1% of all mutation types, respectively, of which katG315 mutations accounted for 64.9% and katG315 deletions for 8.1%. The AhpC promoter region (-44 ~ -30 and − 15 ~ 3 sites) or the inhA promoter region (-17 ~ -8 site) combined with katG315 accounted for the highest double mutation rates of 36.8% and 36.8%, respectively. The molecular patterns associated with MDR-TB were mutations to katG315 combined with codon 529–533 or 521–528 and the inhA promoter region (-17 ~ -8 site) combined with codon 529–533, which accounted for 44.8%, 14.2%, and 13.7% of all mutation patterns, respectively (Fig. 5).

Fig. 5
figure 5

Prevalence of resistance-conferring mutations in genes of MTB strains. The gene coverage region and mutation points is INH A AhpC promoter region (-44 ~ -30 and-15 ~ 3 sites), B inhA94 codon, C inhA promoter region (-17 ~ -8 site), D KatG315 codon. RFP(covering the rpoB gene) A’ 507–512 codon, B’ 521-528codon, C’ 513-520codon, D’ 529-533codon


Prevention and control of MTB infection should be prioritized among high-risk groups [18]. Actually, before HRM testing, clinical screening was performed to rule out low-risk groups of TB infection. The vast majority (95%) of these populations were newly diagnosed patients and from rural areas (82.6%), indicating that individuals residing in rural areas are at higher risk for MTB infection.

HRM analysis is a highly reliable method for identification of MTB infection and highly accurate for detection of drug resistance, particularly to first-line anti-TB drugs [19]. The accuracy of HRM analysis is reportedly greater than 80.0% for resistance to first-line anti-TB drugs and more than 89.0% for second-line anti-TB drugs [20]. Therefore, this method is considered a viable option for rapid diagnosis of DR-TB. In the present study, the HRM positivity rate was higher in males than females (18.0% vs. 14.1%, respectively), indicating that males are at greater risk for MTB infection, which is consistent with previous studies conducted in other regions of China [21]. A possible explanation for this difference in susceptibility between sexes may be related to poorer lifestyle habits of males [22]. The positivity rate of HRM was much higher in groups with positive sputum smears and a history of TB treatment, indicating that positive results of HRM analysis may be closely related to the load of MTB. The HRM positivity rate was much lower in the rural vs. urban group (12.9% vs. 35.3%, respectively), possibly due to differences in healthcare availability in rural areas. Overcrowding in urban areas results in the rapid spread of MTB infection and the continued development of reporting networks and the addition of healthcare centers make the reporting of TB positive cases in urban areas more timely and comprehensive may explain why the detection rate is higher in urban vs. rural areas. In addition, the proportions of previously treated cases and positive sputum smears were lower in rural than urban areas (3.8% vs. 5.8% and 9.5% vs. 25.8%, respectively) (results not shown).

In this study, there were notable differences in the incidence of MTB infections among age groups. Although both males and females participated in HRM testing, most were aged 51–70 years, the highest rate of MTB infection occurred in males aged 26–30 years and females aged 21–25 years, suggesting that young adults might be more susceptible to MTB infection than previously reported, which may be related to hormone-mediated immune responses, increased risk of lung disease, and non-standardized treatments [23, 24]. In addition, young adults are more likely to gather in groups, thereby facilitating the spread of MTB infection [25]. Although young adults are generally considered to be at low risk for MTB infection, additional testing of this group would help to clarify the risk factors associated with the transmission of MTB [26]. Notably, about 2.5 times more males aged 51–65 years undergo HRM testing as compared to females was a main reason why the males MTB positive groups were main at this ages. The rate of MTB infection in males remained greater than 15% in all other age groups, shows, compared with females, the tuberculosis in males throughout the life course is all the focus of management.

Overall, the ratio of MTB-positive males to females was 2.8, while the ratios of males to females with INH-R, RFP-R, and MDR-TB were 3.6, 3.8, and 3.6, respectively, indicating that males were at greater risks for DR-TB not only due to the high ratio of MTB infection in males than females, consistent with a previous study conducted in China [27]. Between 2019 and 2021, the HRM positivity rate as well as the rates of INH-R, RFP-R, MDR tended to decrease, demonstrating that strict control strategies and management of confirmed cases of MTB infection helped to prevent the spread of DR-TB, as least to a certain extent. The decline in drug resistance also reflects the effectiveness of the control strategy recommended by the World Health Organization (i.e., “Directly Observed Treatment, Short-Course”).

Drug resistance was more common in patients previously treated for MTB infection than newly diagnosed cases. Initial drug resistance reflects the ability of local DR-TB to spread and highlights the need for prevention and control strategies. Among newly diagnosed patients, drug resistance tended to decrease over the study period, indicating successful management of primary DR-TB in Luoyang City from 2019 to 2021. However, the MDR and RFP-R rates were higher locally than nationally, indicating the need for continued management. The rate of drug resistance among previously treated patients reflects the rationality for treatment strategies at the local level to control the spread of DR-TB strains. In Luoyang City, the resistance rate was significantly higher in patients previously treated for MTB infection than newly diagnosed cases and the detection rates of INH-R, RFP-R, and MDR-TB tended to increase over the 3-year study period (data not shown), suggesting a continued need to strengthen compliance with medication and management strategies. Phenotypic drug susceptibility methods are time-consuming, leading to a dependence on clinical empiric medication, which also can promote acquired resistance. Therefore, faster and more accurate methods are needed to screen for drug susceptibility, especially in developing countries.

The initial and acquired drug resistance rates were both higher in urban than rural areas, indicating that residents of cities are at greater risk for DR-TB due to the relatively higher population density and mobility [28,29,30]. In addition, the greater availability of healthcare could contribute to the high drug resistance rates in cities. Although the rate of drug resistance was higher in patients with positive vs. negative sputum results, the results of this study showed that the difference in INH-R and MDR detection rates between patients with positive vs. negative sputum results were all not significant. Also, according to the multivariate regression model, the OR value was relatively close to 1 for the association of positive sputum results with an increased risk of drug resistance. Therefore, further studies with larger cohorts are needed to determine whether a positive sputum result is actually a risk factor for drug resistance.

In the multivariate model, history of treatment and younger age were associated with an increased risk for resistance to RFP and INH [31, 32]. Previous studies have found that females are at an increased risk for RFP-R [33]. In contrast, the results of the present study found that males were at an increased risk for RFP-R and INH-R, as well as those aged < 51 years, consistent with the findings of previous studies [33], but inconsistent with the age group of drug resistance in studies conducted in Anhui and Taiwan, China [34, 35]. The higher drug resistance rate in females than males aged 41–45 and 71–75 years may be due to fewer cases of MTB infections among urban males and the relatively high rate among previously treated females in these age groups.

In Luoyang City, the most common genetic mutation associated with RFP-R occurred at codons 529–533 and 521–528, which may be associated with mutations at codon 531 (TCG) and codon 526 (CAC), respectively [36,37,38]. Consistent with previous reports [39, 40], two mutations, katG315 and the inhA promoter region, were associated with INH-R. These findings suggest that both mutations are potential molecular biomarkers for screening of DR-TB. RFP and INH are key anti-TB drugs, thus resistance is likely to negatively affect patient care and the use of hospital resources. RFP-R is an indicator of the development of MDR-TB. The molecular basis for the spread of MDR-TB is mainly associated with katG315, codon 529–533, codon 521–528, and the inhA promoter region. So, blocking the mechanism of transmission via katG315 and codon 529–533 could hinder the emergence and spread of DR-TB.

The aim of this study was to assess HRM detection and TB resistance of confirmed and suspected cases of TB, which only found that the spread of MTB and drug resistance patterns were associated with high-risk populations, at least to a certain extent. However, data regarding TB transmission and resistance analysis in the region were limited, thereby warranting further studies. Besides, gene sequencing was not conducted to confirm the accuracy of the single nucleotide polymorphisms detected by HRM analysis or multiple mutations associated with drug resistance.


The majority of resistant strains were isolated from males, who are at a highest risk for MTB infection in adulthood. Therefore, testing of this population is particularly important. The highest rate of MTB infection occurred in males aged 26–30 years and females aged 21–25 years. Prevention and control of MTB infection in rural populations cannot be taken lightly, although urban populations are also at a high risk of MTB infection and high resistance rates. Previous treatment is a risk factor for drug resistance. In addition, mutations to katG315 and codon 529–533 were associated with DR-TB and, thus, should be considered when selecting a treatment strategy. Overall, populations at high risk for MTB infection include females aged < 51 years, males at all ages, individuals previously treated for MTB infection, and those residing in urban areas.

Availability of data and materials

The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.


  1. Chakaya J, Khan M, Ntoumi F, Aklillu E, Fatima R, Mwaba P, Kapata N, Mfinanga S, Hasnain SE, Katoto PDMC, et al: Global tuberculosis report 2020 - reflections on the global TB burden, treatment and prevention efforts. Int J Infect Dis 2021, 113 Suppl 1.

  2. Lange C, Abubakar I, Alffenaar J-WC, Bothamley G, Caminero JA, Carvalho ACC, Chang K-C, Codecasa L, Correia A, Crudu V, et al. Management of patients with multidrug-resistant/extensively drug-resistant tuberculosis in Europe: a TBNET consensus statement. Eur Respir J. 2014;44(1):23–63.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Ullah I, Shah AA, Basit A, Ali M, Khan A, Ullah U, Ihtesham M, Mehreen S, Mughal A, Javaid A. Rifampicin resistance mutations in the 81 bp RRDR of rpoB gene in mycobacterium tuberculosis clinical isolates using Xpert MTB/RIF in Khyber Pakhtunkhwa, Pakistan: a retrospective study. BMC Infect Dis. 2016;16:413.

    Article  PubMed  PubMed Central  Google Scholar 

  4. 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):321–38.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Yang C, Luo T, Shen X, Wu J, Gan M, Xu P, Wu Z, Lin S, Tian J, Liu Q, et al. Transmission of multidrug-resistant mycobacterium tuberculosis in Shanghai, China: a retrospective observational study using whole-genome sequencing and epidemiological investigation. Lancet Infect Dis. 2017;17(3):275–84.

    Article  PubMed  CAS  Google Scholar 

  6. Kohli M, Schiller I, Dendukuri N, Yao M, Dheda K, Denkinger CM, Schumacher SG, Steingart KR. Xpert MTB/RIF Ultra and Xpert MTB/RIF assays for extrapulmonary tuberculosis and rifampicin resistance in adults. Cochrane Database Syst Rev. 2021;1:CD012768.

    PubMed  Google Scholar 

  7. Ramirez MV, Cowart KC, Campbell PJ, Morlock GP, Sikes D, Winchell JM, Posey JE. Rapid detection of multidrug-resistant mycobacterium tuberculosis by use of real-time PCR and high-resolution melt analysis. J Clin Microbiol. 2010;48(11):4003–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Global, regional, and national burden of tuberculosis, 1990–2016: results from the global burden of diseases, injuries, and risk factors 2016 study. Lancet Infect Dis 2018, 18(12):1329–1349.

  9. Prasad BM, Satyanarayana S, Chadha SS. Lessons learnt from active tuberculosis case finding in an urban slum setting of Agra city, India. Indian J Tuberc. 2016;63(3):199–202.

    Article  PubMed  Google Scholar 

  10. Limenih YA, Workie DL. Survival analysis of time to cure on multi-drug resistance tuberculosis patients in Amhara region, Ethiopia. BMC Public Health. 2019;19(1):165.

    Article  PubMed  PubMed Central  Google Scholar 

  11. McLaren ZM, Brouwer E, Ederer D, Fischer K, Branson N. Gender patterns of tuberculosis testing and disease in South Africa. Int J Tuberc Lung Dis. 2015;19(1):104–10.

    Article  PubMed  CAS  Google Scholar 

  12. Snow KJ, Sismanidis C, Denholm J, Sawyer SM, Graham SM. The incidence of tuberculosis among adolescents and young adults: a global estimate. Eur Respir J 2018, 51(2).

  13. Asgedom SW, Teweldemedhin M, Gebreyesus H. Prevalence of multidrug-resistant tuberculosis and associated factors in ethiopia: a systematic review. J Pathog. 2018;2018:7104921.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Chen MP, Miramontes R, Kammerer JS. Multidrug-resistant tuberculosis in the United States, 2011–2016: patient characteristics and risk factors. Int J Tuberc Lung Dis. 2020;24(1):92–9.

    Article  PubMed  CAS  Google Scholar 

  15. WHO. revised definitions and reporting framework for tuberculosis. Euro Surveill. 2013;18(16):20455.

    Google Scholar 

  16. Long Q, Qu Y, Lucas H. Drug-resistant tuberculosis control in China: progress and challenges. Infect Dis Poverty. 2016;5:9.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Andre E, Goeminne L, Cabibbe A, Beckert P, Kabamba Mukadi B, Mathys V, Gagneux S, Niemann S, Van Ingen J, Cambau E. Consensus numbering system for the rifampicin resistance-associated rpoB gene mutations in pathogenic mycobacteria. Clin Microbiol Infect. 2017;23(3):167–72.

    Article  PubMed  CAS  Google Scholar 

  18. Yuen CM, Kammerer JS, Marks K, Navin TR, France AM. Recent transmission of tuberculosis - United States, 2011–2014. PLoS ONE. 2016;11(4):e0153728.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Keikha M, Karbalaei M. High resolution melting assay as a reliable method for diagnosing drug-resistant TB cases: a systematic review and meta-analysis. BMC Infect Dis. 2021;21(1):989.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Liu Z, Gao Y, Yang H, Bao H, Qin L, Zhu C, Chen Y, Hu Z. Impact of hypoxia on drug resistance and growth characteristics of mycobacterium tuberculosis clinical isolates. PLoS ONE. 2016;11(11):e0166052.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Minero G, Bagnasco M, Fock J, Tian B, Garbarino F, Hansen MF. Automated on-chip analysis of tuberculosis drug-resistance mutation with integrated DNA ligation and amplification. Anal Bioanal Chem. 2020;412(12):-.

  22. Okethwangu D, Birungi D, Biribawa C, Kwesiga B, Turyahabwe S, Ario AR, Zhu B-P. Multidrug-resistant tuberculosis outbreak associated with poor treatment adherence and delayed treatment: Arua District, Uganda, 2013–2017. BMC Infect Dis. 2019;19(1):387.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Zhao Y, Ying H, Demei J, Xie J. Tuberculosis and sexual inequality: the role of sex hormones in immunity. Crit Rev Eukaryot Gene Expr. 2012;22(3):233–41.

    Article  PubMed  CAS  Google Scholar 

  24. Borgdorff MW, Nagelkerke NJ, Dye C, Nunn P. Gender and tuberculosis: a comparison of prevalence surveys with notification data to explore sex differences in case detection. Int J Tuberc Lung Dis. 2000;4(2):123–32.

    PubMed  CAS  Google Scholar 

  25. Tao N-N, He X-C, Zhang X-X, Liu Y, Yu C-B, Li H-C. Trends and characteristics of drug-resistant tuberculosis in rural Shandong, China. Int J Infect Dis 2017, 65.

  26. Snow KJ, Nelson LJ, Sismanidis C, Sawyer SM, Graham SM. Incidence and prevalence of bacteriologically confirmed pulmonary tuberculosis among adolescents and young adults: a systematic review. Epidemiol Infect. 2018;146(8):946–53.

    Article  PubMed  Google Scholar 

  27. An Q, Song W, Liu J, Tao N, Liu Y, Zhang Q, Xu T, Li S, Liu S, Li Y, et al. Primary drug-resistance pattern and trend in elderly tuberculosis patients in Shandong, China, from 2004 to 2019. Infect Drug Resist. 2020;13:4133–45.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Marx FM, Fiebig L, Hauer B, Brodhun B, Glaser-Paschke G, Magdorf K, Haas W. Higher rate of tuberculosis in second generation migrants compared to native residents in a metropolitan setting in Western Europe. PLoS ONE. 2015;10(6):e0119693.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Mizukoshi F, Kobayashi N, Kirikae F, Ohta K, Tsuyuguchi K, Yamada N, Inoue Y, Horiba M, Kawata N, Ichinose A, et al: Molecular epidemiology of drug-resistant mycobacterium tuberculosis in Japan. mSphere 2021, 6(4):e0097820.

  30. Chen H, He L, Huang H, Shi C, Ni X, Dai G, Ma L, Li W. Mycobacterium tuberculosis lineage distribution in Xinjiang and Gansu Provinces, China. Sci Rep. 2017;7(1):1068.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Pradipta IS, Forsman LD, Bruchfeld J, Hak E, Alffenaar J-W. Risk factors of multidrug-resistant tuberculosis: A global systematic review and meta-analysis. J Infect. 2018;77(6):469–78.

    Article  PubMed  Google Scholar 

  32. Nair SA, Raizada N, Sachdeva KS, Denkinger C, Schumacher S, Dewan P, Kulsange S, Boehme C, Paramsivan CN, Arinaminpathy N. Factors associated with tuberculosis and rifampicin-resistant tuberculosis amongst symptomatic patients in India: A Retrospective Analysis. PLoS ONE. 2016;11(2):e0150054.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Seifert M, Aung HT, Besler N, Harris V, Mar TT, Colman RE, Rodwell TC, Aung ST. Age and sex distribution of mycobacterium tuberculosis infection and rifampicin resistance in Myanmar as detected by Xpert MTB/RIF. BMC Infect Dis. 2021;21(1):781.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. Chuang P-H, Wu M-H, Fan S-Y, Lin K-Y, Jou R. Population-based drug resistance surveillance of multidrug-resistant tuberculosis in Taiwan, 2007–2014. PLoS ONE. 2016;11(11):e0165222.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Yao S, Yan J, Li L, Ma D, Liu J, Wang Q, Wang A, Bao F, Zhang Z, Bao X: Determining drug resistance and risk factors for multidrug-resistant tuberculosis in sputum smear-positive tuberculosis outpatients in Anhui Province, China, 2015–2016. Infect Drug Resist 2020, 13:1023–1032.

  36. Miotto P, Tessema B, Tagliani E, Chindelevitch L, Starks AM, Emerson C, Hanna D, Kim PS, Liwski R, Zignol M, et al: A standardised method for interpreting the association between mutations and phenotypic drug resistance in. Eur Respir J 2017, 50(6).

  37. Jagielski T, Bakuła Z, Brzostek A, Minias A, Stachowiak R, Kalita J, Napiórkowska A, Augustynowicz-Kopeć E, Żaczek A, Vasiliauskiene E, et al: characterization of mutations conferring resistance to rifampin in mycobacterium tuberculosis clinical strains. Antimicrob Agents Chemother 2018, 62(10).

  38. Rahman A, Sahrin M, Afrin S, Earley K, Ahmed S, Rahman SMM, Banu S. Comparison of xpert MTB/RIF Assay and GenoType MTBDRplus DNA probes for detection of mutations associated with rifampicin resistance in mycobacterium tuberculosis. PLoS ONE. 2016;11(4):e0152694.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Georghiou SB, Seifert M, Catanzaro D, Garfein RS, Valafar F, Crudu V, Rodrigues C, Victor TC, Catanzaro A, Rodwell TC. Frequency and distribution of tuberculosis resistance-associated mutations between Mumbai, Moldova, and Eastern Cape. Antimicrob Agents Chemother. 2016;60(7):3994–4004.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Fenner L, Egger M, Bodmer T, Altpeter E, Zwahlen M, Jaton K, Pfyffer GE, Borrell S, Dubuis O, Bruderer T, et al. Effect of mutation and genetic background on drug resistance in Mycobacterium tuberculosis. Antimicrob Agents Chemother. 2012;56(6):3047–53.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references


Special thanks to Luoyang CDC for providing part of the data for the study.


This work was supported by the National Nature Science Foundation of China (32072899), Henan Science and Technology Development Plan (212102310741, 212102310887) and Joint construction project of Henan Province (LHGJ20200589).

Author information

Authors and Affiliations



XY and ZZ directed data analysis and article structure conception. ZW carried out paper writing and raw data collation and analysis. TG, TJ, QZ, and YH participating in data collation and analysis. XZ and SK coached and participated in table design and production. LL provided support for relevant technologies.

Corresponding authors

Correspondence to Zhanqin Zhao or Yun Xue.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Review Committee of The First Affiliated Hospital and Clinical Medical College, Henan University of Science and Technology (Approval Registration Number − 2022/03/B070). The need for informed consent was waived by this Ethics Committee considering the retrospective nature of the study and the anonymous clinical data used for analysis. At the same time, all analyses were performed in accordance with guidelines and regulations for handling de-identified data.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence 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 licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, Z., Guo, T., Jiang, T. et al. Regional distribution of Mycobacterium tuberculosis infection and resistance to rifampicin and isoniazid as determined by high-resolution melt analysis. BMC Infect Dis 22, 812 (2022).

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI:


  • Multidrug-resistant tuberculosis
  • Transmission
  • Molecular epidemiology
  • Resistance model