Introduction

Tuberculosis (TB) remains one of the leading causes of death worldwide from a single infectious agent, and is caused by the Mycobacterium tuberculosis (MTB) bacillus.1 The World Health Organization (WHO) aims to reduce mortality and incidence rates of TB by 95% and 90%, respectively, by 2035, compared with rates in 2015; however, the current TB rate of incidence remains concerningly higher than the rate of TB elimination.2 China was observed to have the third highest TB burden in the world (6.8%) in 2022, behind Indonesia (10%) and India (27%).2 Additionally, TB is also the second leading cause of morbidity and mortality from class A and B infectious diseases in China as per The Infectious Disease Prevention and Control Acts of China.3 Despite recent advances in the management of TB, drug-resistant TB continues to be an onerous public health threat in China.2 Moreover, the prevalence and molecular characteristics of multidrug resistant TB (MDR-TB) are known to vary with geographical environment and bacterial lineage.4 Therefore, a thorough understanding of the genetic basis of drug resistance is necessary for the successful management of MTB-related disease.

Isoniazid (INH) is one of the key drugs for TB treatment due to its early bactericidal effect, and is used to treat both active tuberculosis (TB) and latent TB infection (LTBI) for several decades.5 Both multidrug resistance and resistance to INH without concomitant rifampin resistance are associated with poor response to first-line treatment.6 Resistance to INH is mainly driven by mutations in the katG gene.7 The katG gene encodes for a catalase-peroxidase enzyme in MTB necessary for INH activation into its active form to exert its lethal effect.8 It is well-testified fact that katG gene is most frequently getting mutated and hence promotes INH resistance by diminishing the normal activity of catalase-peroxidase enzyme.9

The ahpC gene encodes alkyl-hydroperoxide reducatse also involved in cellular regulation of oxidative stress.10,11 The loss of katG function in resistant strains is thought to be compensated for by mutations in the oxyR-ahpC intergenic region, which contains the putative promoter of ahpC.12 Mutations in the oxyR-ahpC gene region may lead to the upregulation of ahpC gene expression, enhancing the cell tolerance to isoniazid-induced oxidative stress, and thereby promoting the development of drug resistance.13,14 However, the role of ahpC in INH resistance is controversial. The mutation prevalence in the intergenic region of oxyR-ahpC has shown regional and phenotypic differences.15–17 Previous studies have shown that mutations in the oxyR-ahpC intergenic region may compensate for the loss of catalase-peroxidase activity caused by mutations in katG, and thus may protect MTB bacilli against the toxic effects of organic peroxides.15,17 However, one subsequent study observed that mutations in oxyR-ahpC were not frequently seen in INH-resistant isolates with katG mutations.18 Baker et al reported that mutations in the oxyR-ahpC region did not lead to INH resistance, as these mutations can be found in 23.7% of INH-resistant clinical isolates and in 7.5% of INH-sensitive isolates.16

Chongqing is a high-incidence area for MDR-TB in China. To further verify the role of mutation in the oxyR-ahpC intergenic region in INH resistance, 490 clinical isolates were collected from Chongqing area, the prevalence and distribution of mutations in katG, ahpC and the oxyR-ahpC intergenic region was characterized.

Materials and MethodsClinical MTB Isolates

A total of 490 clinical sputum samples were collected from patients, principally from the Chongqing and surrounding areas, attending the Chongqing Public Health Medical Center between June 2020 and March 2022. Samples were decontaminated and subjected to bacterial culture for 6 weeks using the BACTEC MGIT 960 (Becton Dickinson, Sparks, MD, United States) culture system at Chongqing Public Health Medical Center (Chongqing, China).19 According to the latest WHO guidelines,2 MDR-TB was defined as resistance to at least both isoniazid and rifampicin. Pre-extensively drug-resistant TB (pre-XDR-TB) was defined as MDR-TB with additional resistance to any fluoroquinolone, moxifloxacin, levofloxacin. Furthermore, extensively drug-resistant TB (XDR-TB) is defined as MDR-TB plus resistance to one of the fluoroquinolones and to one of the second-line injectable anti-TB drugs: amikacin, kanamycin. The positive culture samples were subjected to further drug susceptibility testing using proportion method on Lowenstein-Jensen (L-J) medium slants (Encode, Zhuhai, China). The DST panel included first-line drugs (isoniazid, rifampicin) and second-line drugs (moxifloxacin, levofloxacin, amikacin, capreomycin, prothionamide, para-aminosalicylic acid, etc). The DST was conducted using the proportional method on L-J medium for 4—8 weeks.20 Isolates were then frozen in 25% glycerol at −80°C until use.

Mycobacterial Genomic DNA

Single colonies of isolates on L-J medium slants were inoculated into 15 mL centrifuge tubes and suspended in 3 mL of Middlebrook 7H9 liquid medium with 10% (v/v) OADC, and was cultured at 37°C for 3–4 weeks. Cultures were grown to early stationary phase and subsequently used for genomic DNA preparation. Briefly, the mycobacterial cell pellet was washed twice with deionized water, and re-suspended with deionized water containing 0.45% (v/v) Triton X-100 and 0.45% (v/v) Tween-80. The bacterial cells were then incubated at 95°C for 30 min. Subsequently, 0.5 mL of the suspension was centrifuged to collect the supernatant. The extracted supernatant containing MTB genomic DNA was stored at − 20°C until use.

Gene Amplification and Sequencing Analysis

PCR was performed to amplify the fragments of the katG and oxyR-ahpC genes, as previously described.21,22 The katG fragment was amplified with the following primers: katG-F 5′-GGTCGACATTCGCGAGACGTT-3′ and katG-R 5′-TTGTTCCTGCGACGCATCGTG-3′. The genomic DNA was used as a template for amplification. The oxyR-ahpC fragment was amplified with the following primers: oxyR-ahpC-F 5′-GCCTGGGTGTTCGTCACTGGT-3′ and oxyR-ahpC-R 5′- CGCAACGTCGACTGGCTCATA-3′. PCR amplifications were performed using AceTaq Master Mix (Vazyme, P411, Nanjing, China). The reaction mixtures were subjected to 5 min at 95°C, followed by 30 cycles of 10s at 95°C, 30s at 56°C, 45s at 72°C, and were terminated by an additional 5 min at 72°C. Successful gene amplifications were confirmed by ultraviolet transillumination following electrophoresis on a 1% agarose gel stained with ethidium bromide in a 1☓TAE buffer. Sequencing of PCR products was conducted at Tsingke Biotech (Chongqing, China). Gene polymorphisms were aligned with the katG and oxyR-ahpC genes of the H37Rv reference strain (GenBank ID: NC_000962.3) using SnapGene software (SnapGene software, www.snapgene.com, version 4.1).

ResultsDemographic Characteristics of the Patients from Whom the MTB Clinical Isolates Were Collected

A total of 490 clinical isolates from different patients with TB were collected, which included 73 INH-susceptible isolates, 26 INH-mono-resistant isolates, 199 MDR isolates, 190 pre-XDR isolates, and two XDR isolates. Three hundred and forty-one isolates (69.6%) were sampled from male patients and 149 isolates (30.4%) were taken from female patients, resulting in a sex ratio of 2.29. The median patient age was 47 years (range 13–83) (Table 1).

Table 1 Classification of 490 MTB Clinical Isolates Used in This Study

Drug Sensitivity Tests (DSTs)

DSTs were conducted on 490 MTB clinical isolates, and we categorized 417 of the 490 isolates (85.10%) to have resistance against INH, including 26 of 417 (6.24%) mono-INH resistant, 199 of 417 (47.72%) MDR, 190 of 417 (45.56%) pre-XDR, and two of 417 (0.48%) XDR-MTB isolates. The remaining 73 (14.90%) of the 490 isolates were observed to consistently be susceptible to all the fourteen anti-tuberculosis drugs used in the test panel (Figure 1A).

Figure 1 Drug resistance profile of 490 Mycobacterium tuberculosis clinical isolates. (A) Distribution of drug susceptibility profiles among the collected isolates. The pie chart shows the proportion of pan-susceptible, Mono-INH, MDR, pre-XDR, and XDR isolates. (B) Drug resistance patterns among MDR and pre-XDR/XDR isolates. The bar graph depicts the number of MDR (n=199) and pre-XDR/XDR (n=192) isolates that are resistant to each of the fourteen anti-tuberculosis drugs tested.

Resistance to rifapentine, streptomycin, para-aminosalicylic acid, rifabutin, amikacin, ethambutol, capreomycin, prothionamide, clofazimine and clarithromycin was observed in 184 (92.5%), 164 (82.4%), 62 (31.2%),24 (12.1%), 20 (10.1%), 10 (5.0%),3 (1.5%), 1 (0.5%), 2 (1.0%), and 2 (1.0%) isolates, respectively, among the 199 MDR isolates (Figure 1B). Additionally, resistance to rifapentine, streptomycin, para-aminosalicylic acid, amikacin, ethambutol, rifabutin, prothionamide, clarithromycin, clofazimine, capreomycin, and linezolid was observed in 181 (94.3%), 161 (83.9%), 85 (44.3%), 38 (19.8%), 21 (10.9%), 16 (8.3%), 9 (4.7%), 4 (2.1%), 3 (1.6%), and 2 (1.0%) isolates, respectively, among the 190 pre-XDR and 2 XDR isolates (Figure 1B).

Mutations in katG, oxyR and the oxyR-ahpC Intergenic Region

Among the 73 pan-susceptible clinical isolates, 69 of 73 isolates (94.5%) had no mutations, but four (5.5%) harbored two synonymous mutations (Pro375Pro, Phe272Phe) and two non-synonymous mutations (Ile364Val, Asn508Ser) in the katG gene (Table 2). Among the 417 INH-resistant clinical isolates, 382 (91.6%) had a canonical mutation in katG at codon 315, of which S315T (G944C) substitution was the most common (379 of 382, 99.2%).

Table 2 Mutations in katG, oxyR and oxyR-ahpC Intergenic Region of the Clinical Isolates

Among the 26 INH-monoresistant isolates, 18 of 26 (69.23%) had a canonical mutation in katG (S315T) (Table 2). Three other katG mutations were identified in three INH-monoresistant isolates, ie, G→T at 859 (Glu287STOP), G insertion at 1002, and G→T at 1390 (Arg465Leu).

Among the 199 MDR isolates, 190 (95.48%) possessed mutations in katG and oxyR-ahpC (Table 2). The most frequent change observed in katG was S315T (187 of 199, 93.97%). Four isolates were observed to have synonymous mutations in either katG or oxyR, ie, Pro501Pro (one isolate) and Gly428Gly (three isolates), and seven isolates contained a mutation (Ala18Gly) in oxyR. We also observed that four isolates contained double mutations, ie, Ser315Asn with Ile317Val in katG (one isolate), and S315T in katG with −58 G→A in oxyR (three isolates).

With respect to the 192 pre-XDR and XDR isolates, 179 of 192 (93.22%) possessed mutations in katG and oxyR-ahpC (Table 2). Among these 179 isolates, 176 contained the katG S315T mutation (176 of 192, 91.67%). Only one isolate contained a synonymous mutation (Gly428Gly) in katG, and 13 isolates had a mutation (Ala18Gly) in oxyR.

Discussion

In this study, a total of 490 clinical samples were collected from different patients afflicted with TB, with a male–female ratio of 2.29, which is a similar relative proportion of TB prevalence to that observed by the WHO globally, and also to that seen in the national survey of China.2,13

From our sample of 417 INHR clinical isolates, 382 (91.6%) could be explained by mutations in the 315 codon of katG, of which 379 (99.2%) contained the S315T transition. It has been previously shown that MTB strains possessing the S315T mutation are fully virulent in mice,23 and this may be related to the increased capacity for transmission of the S315T strain.24 Consequently, the S315T mutation is more frequently observed in MDR-TB patients than in INH mono-resistant clinical isolates,25 and is considered as a harbinger mutation that often precedes MDR.26 Our results further support this hypothesis. A clear trend was observed in our study, where the proportion of katG 315 mutations in MDR (93.97%) and Pre-XDR/XDR (91.67%) isolates was significantly higher than that in INH-monoresistant (69.23%) isolates. Cohen et al, also showed that the overwhelming majority of MDR-TB and XDR-TB organisms evolved resistance to isoniazid prior to resistance to rifampicin.27 Additionally, analysis of phenotypic drug susceptibility testing by Izu et al, has revealed that isoniazid resistance is acquired before rifampicin resistance.28 Collectively, the preceding observations suggest that some MDR-TB bacilli may have evolved from isoniazid-resistant strains. De novo emergence and reinfection by MDR-TB strains contribute equally to MDR development.27

In the present study, a total of 6 mutation sites (in 10 strains, and other than at the 315 site) in the katG gene were observed in 417 INH-resistant isolates, five of which were also accompanied by mutation at the 315 site. We have thus shown that the katG 315 gene mutation remains a major contributor to the earlier evolution of isoniazid resistance.

The low frequency (1.5%, 3/199 and 1.04%, 2/192) of oxyR-ahpC intergenic mutations observed in our study, even among MDR and pre-XDR/XDR isolates, provides important epidemiological context. While these compensatory mutations are theoretically advantageous for strains bearing the fitness cost of katG mutations,29 their low prevalence in our high MDR/XDR burden setting implies that they may not be the prerequisite for the evolution and success of complex drug-resistant strains. This also indicates that other, as yet uncharacterized, genetic or metabolic pathways might be more critical or common in compensating the katG defect in this population, allowing MDR/XDR strains to maintain transmissibility and virulence.

While our data firmly establish katG S315T (91.6%, 382/417) as the dominant mechanism of INH resistance in this population, the role of other genetic regions remains critical to a comprehensive understanding. Mutations in the katG gene usually confer a high level of INH resistance, which causes decreased INH substrate affinity at the expense of a reduction in overall catalase activity.7 Since the MTB katG enzyme (catalase peroxidase) is one of the major participants in the detoxification of reactive oxygen species (ROS) in MTB, katG mutations that confers INH resistance have been shown to influence bacterial pathogenesis and responses to antibiotic treatment.30–32 This implies that INH resistance-conferring katG mutations have fitness costs for MTB. Thus, it is reasonable to assume that INH-resistant MTB clinical isolates with katG mutations may have other compensatory mutations. Previous studies have shown that mutations in the oxyR-ahpC intergenic region may compensate for the loss of katG or catalase activity caused by katG mutations, which may result in over-expression of ahpC.29 The ahpC gene encodes another important component of the ROS scavenger system in MTB,33 and over-expression of the gene results in low-level INH resistance.11 Additionally, mutations in the oxyR-ahpC intergenic region have been observed to contribute to 10—15% of INH resistance.34–36 However, the frequency of mutation in the oxyR-ahpC area was observed to be rather low (1.5% in MDR and 1.04% in pre-XDR/XDR) in the present study, and all the preceding mutations were identified in MDR isolates with the katG mutation. This finding is consistent with the controversial and likely secondary role of oxyR-ahpC region in INH resistance, supporting the view that these mutations may act as one potential compensatory mechanism for the loss of catalase activity induced by katG mutations in some strains. The complete absence of these mutations in INH-susceptible and INH mono-resistant isolates in our study further supports the model that they do not significantly contribute to isoniazid resistance by themselves, but may play a minor, compensatory role in strains already harboring primary resistance mutations such as katG S315T. However, over-expression of ahpC via mutations in oxyR-ahpC might not be the primary route that compensates for the deficiency caused by the katG mutation, suggesting that the bacterium may compensate for the deficiency caused by the katG mutation via other mechanisms, and this merits further investigation. Beyond confirming the established role of katG and the limited contribution of the oxyR-ahpC intergenic region, our study brings to light a previously underreported genetic change. A considerable proportion of INH-resistant isolates (3.52% in MDR isolates and 6.77% in pre-XDR isolates) with the katG mutation were observed to have a novel mutation in the oxyR gene (Ala18Gly of oxyR), which, to our knowledge, has not been previously observed in the contemporary literature, it is possible that this specific mutation may also compensate for the loss of catalase activity caused by katG mutations. Although the functional consequence of this specific amino acid change is unknown, its location in the oxyR gene, which encodes a redox-sensitive transcription regulator, suggests it could potentially alter the expression of oxidative stress response genes including ahpC. This hypothesis positions the oxyR Ala18Gly mutation as a candidate compensatory mutation, warranting dedicated functional validation in future studies.

A key limitation of this study is the lack of paired phenotypic data (such as catalase activity assays or in vitro growth kinetics) so that we could perform direct functional comparisons between strains carrying only katG mutations and strains carrying both katG and oxyR/oxyR-ahpC mutations. These experiments will be critical to determine the compensatory effects of these mutations. Furthermore, a retrospective analysis of DST data showed that no specific resistance patterns were found among isolates carrying both katG and oxyR-ahpC mutations compared with isolates carrying only katG mutations (Supplementary Tables S1 and S2). This observation suggests that potential compensatory effects of the oxyR-ahpC mutation may not translate directly into broad changes in the phenotypic resistance spectrum detected by standard DST, but may instead fine-tune bacterial fitness. Future studies are needed to combine these phenotypic assays to elucidate the precise biological costs and compensatory mechanisms involved.

In our study, although the incidence of katG mutations was higher, the oxyR-ahpC intergenic mutation rate was significantly lower (1.5% for MDR and 1.04% for pre-XDR/XDR), indicating that ahpC upregulation through this promoter region is not a major compensatory pathway in the clinical MTB strains prevalent in Chongqing. This observation raises the question of how these adapted drug-resistant strains overcome the vulnerability to oxidative stress caused by katG deficiency. It is possible that other, more potent genetic or metabolic adaptations are at work. These may include mutations in other regulators of the oxidative stress response, enhanced activity of additional detoxifying enzymes such as superoxide dismutase.17,37 Identifying these alternative mechanisms represents an important goal for future studies of the evolution of drug-resistant tuberculosis.

Conclusion

The present study supports the utility of the katG 315 mutation as a predictive molecular marker for INH resistance. Our findings support the primacy of katG S315T in molecular diagnosis in Chongqing and suggest that routine screening of the oxyR-ahpC intergenic region may have limited added value due to its low prevalence. Mutations in the oxyR-ahpC intergenic region exist selectively and in low frequency in INH-resistant isolates with katG mutations, thus supporting their putative role in compensating for the catalase deficiency caused by katG mutations. However, mutation in the oxyR-ahpC intergenic region seems unlikely to be the principal route which compensates for the catalase deficiency caused by katG mutations in MTB clinical isolates. Furthermore, one novel mutation oxyR Ala18Gly identified in this study warrants further investigation to elucidate its potential compensatory function. Future research directions include: 1) verifying the compensatory function by testing the in vitro antioxidant capacity of the oxyR Ala18Gly mutant; 2) conducting transcriptomic studies in strains carrying mutations in the oxyR-ahpC intergenic region to explore the expression of a wider range of oxidative stress regulators; 3) comprehensive genome sequencing to identify other unknown compensatory mutations. The outcomes of the present study contribute substantially to our knowledge with respect to the spectrum of gene mutations that may participate in INH resistance, and provide tantalizing clues for the elucidation of the drug-resistance mechanisms that are employed by MTB.

Ethics Approval and Consent to Participate

This study used archived strains isolated from clinical samples, no human subjects and no identifiable human data were included. The authors obtained approval for this study and exemption from the requirement for written informed patient consent from the Ethics Committee of Chongqing Public Heath Medical Center (Permit Number: 2022-093-01-KY). During the study, patient data confidentiality and compliance with the Declaration of Helsinki were followed.

Funding

This work was supported by the Chongqing medical scientific research project (Joint project of Chongqing Health Commission and Science and Technology Bureau) (2024MSXM120, 2024QNXM020, 2025QNXM052).

Disclosure

The authors declare that this research was conducted in the absence of any commercial or financial relationships that could be construed as potential conflicts of interest.

References

1. Li J, Liu XQ, Jiang SW, et al. Improving tuberculosis case detection in underdeveloped multi-ethnic regions with high disease burden: a case study of integrated control program in China. Infect Dis Poverty. 2017;6(1):151. doi:10.1186/s40249-017-0365-4

2. World Health Organization. Global Tuberculosis Report 2024. Geneva: World Health Organization; 2024. Licence: CC BY-NC-SA 3.0 IGO.

3. He WC, Ju K, Gao YM, et al. Spatial inequality, characteristics of internal migration, and pulmonary tuberculosis in China, 2011-2017: a spatial analysis. Infect Dis Poverty. 2020;9(1):159. doi:10.1186/s40249-020-00778-0

4. Yang C, Luo T, Sun G, et al. Mycobacterium tuberculosis Beijing strains favor transmission but not drug resistance in China. Clin Infect Dis. 2012;55(9):1179–1187. doi:10.1093/cid/cis670

5. Murray JF, Schraufnagel DE, Hopewell PC. Treatment of tuberculosis. A historical perspective. Ann Am Thorac Soc. 2015;12(12):1749–1759. doi:10.1513/AnnalsATS.201509-632PS

6. Menzies D, Benedetti A, Paydar A, et al. Standardized treatment of active tuberculosis in patients with previous treatment and/or with mono-resistance to isoniazid: a systematic review and meta-analysis. PLoS Med. 2009;6(9):e1000150. doi:10.1371/journal.pmed.1000150

7. Lyu J, Wu W, Cheng P, et al. A chip for detecting tuberculosis drug resistance based on Polymerase Chain Reaction (PCR)-magnetic bead molecule platform. Front Microbiol. 2018;9:2106. doi:10.3389/fmicb.2018.02106

8. Zhang Y, Heym B, Allen B, Young D, Cole S. The catalase-peroxidase gene and isoniazid resistance of Mycobacterium tuberculosis. Nature. 1992;358(6387):591–593. doi:10.1038/358591a0

9. Ando H, Kondo Y, Suetake T, et al. Identification of katG mutations associated with high-level isoniazid resistance in Mycobacterium tuberculosis. Antimicrob Agents Chemother. 2010;54(5):1793–1799. doi:10.1128/AAC.01691-09

10. Dookie N, Rambaran S, Padayatchi N, Mahomed S, Naidoo K. Evolution of drug resistance in Mycobacterium tuberculosis: a review on the molecular determinants of resistance and implications for personalized care. J Antimicrob Chemother. 2018;73(5):1138–1151. doi:10.1093/jac/dkx506

11. Heym B, Stavropoulos E, Honoré N, et al. Effects of overexpression of the alkyl hydroperoxide reductase AhpC on the virulence and isoniazid resistance of Mycobacterium tuberculosis. Infect Immun. 1997;65(4):1395–1401. doi:10.1128/iai.65.4.1395-1401.1997

12. Kelley CL, Rouse DA, Morris SL. Analysis of ahpC gene mutations in isoniazid-resistant clinical isolates of Mycobacterium tuberculosis. Antimicrob Agents Chemother. 1997;41(9):2057–2058. doi:10.1128/AAC.41.9.2057

13. Dalla Costa ER, Ribeiro MO, Silva MSN, et al. Correlations of mutations in katG, oxyR-ahpC and inhA genes and in vitro susceptibility in Mycobacterium tuberculosis clinical strains segregated by spoligotype families from tuberculosis prevalent countries in South America. BMC Microbiol. 2009;9:39. doi:10.1186/1471-2180-9-39

14. Norouzi F, Moghim S, Farzaneh S, Fazeli H, Salehi M, Nasr Esfahani B. Significance of the coexistence of non-codon 315 katG, inhA, and oxyR-ahpC intergenic gene mutations among isoniazid-resistant and multidrug-resistant isolates of Mycobacterium tuberculosis: a report of novel mutations. Pathog Glob Health. 2022;116(1):22–29. doi:10.1080/20477724.2021.1928870

15. Rinder H, Thomschke A, Rüsch-Gerdes S, et al. Significance of ahpC promoter mutations for the prediction of isoniazid resistance in Mycobacterium tuberculosis. Eur J Clin Microbiol Infect Dis. 1998;17(7):508–511. doi:10.1007/BF01691135

16. Baker LV, Brown TJ, Maxwell O, et al. Molecular analysis of isoniazid-resistant Mycobacterium tuberculosis isolates from England and Wales reveals the phylogenetic significance of the ahpC −46A polymorphism. Antimicrob Agents Chemother. 2005;49(4):1455–1464. doi:10.1128/AAC.49.4.1455-1464.2005

17. Sherman DR, Mdluli K, Hickey MJ, et al. Compensatory ahpC gene expression in isoniazid-resistant Mycobacterium tuberculosis. Science. 1996;272(5268):1641–1643. doi:10.1126/science.272.5268.1641

18. Sreevatsan S, Pan X, Zhang Y, Deretic V, Musser JM. Analysis of the oxyR-ahpC region in isoniazid-resistant and -susceptible Mycobacterium tuberculosis complex organisms recovered from diseased humans and animals in diverse localities. Antimicrob Agents Chemother. 1997;41(3):600–606. doi:10.1128/AAC.41.3.600

19. Huang TS, Lee SSJ, Tu HZ, et al. Use of MGIT 960 for rapid quantitative measurement of the susceptibility of Mycobacterium tuberculosis complex to ciprofloxacin and ethionamide. J Antimicrob Chemother. 2004;53(4):600–603. doi:10.1093/jac/dkh120

20. Woods GL, Weinstein MP, Woods GL. Susceptibility testing for mycobacteria. Clin Infect Dis. 2000;31(5):1209–1215. doi:10.1086/317441

21. Choi JH, Lee KW, Kang HR, et al. Clinical efficacy of direct DNA sequencing analysis on sputum specimens for early detection of drug-resistant Mycobacterium tuberculosis in a clinical setting. Chest. 2010;137(2):393–400. doi:10.1378/chest.09-0150

22. Silva MSN, Senna SG, Ribeiro MO, et al. Mutations in katG, inhA, and ahpC genes of Brazilian isoniazid-resistant isolates of Mycobacterium tuberculosis. J Clin Microbiol. 2003;41(9):4471–4474. doi:10.1128/JCM.41.9.4471-4474.2003

23. Pym AS, Saint-Joanis B, Cole ST. Effect of katG mutations on the virulence of Mycobacterium tuberculosis and the implication for transmission in humans. Infect Immun. 2002;70(9):4955–4960. doi:10.1128/IAI.70.9.4955-4960.2002

24. Gagneux S, Burgos MV, DeRiemer K, et al. Impact of bacterial genetics on the transmission of isoniazid-resistant Mycobacterium tuberculosis. PLoS Pathog. 2006;2(6):e61. doi:10.1371/journal.ppat.0020061

25. van Soolingen D, de Haas PE, van Doorn HR, Kuijper E, Rinder H, Borgdorff MW. Mutations at amino acid position 315 of the katG gene are associated with high-level resistance to isoniazid, other drug resistance, and successful transmission of Mycobacterium tuberculosis in the Netherlands. J Infect Dis. 2000;182(6):1788–1790. doi:10.1086/317598

26. Manson AL, Cohen KA, Abeel T, et al. Genomic analysis of globally diverse Mycobacterium tuberculosis strains provides insights into the emergence and spread of multidrug resistance. Nat Genet. 2017;49(3):395–402. doi:10.1038/ng.3767

27. Srinivasan V, Ha VTN, Vinh DN, et al. Sources of multidrug resistance in patients with previous isoniazid-resistant tuberculosis identified using whole genome sequencing: a longitudinal cohort study. Clin Infect Dis. 2020;71(10):e532–e539. doi:10.1093/cid/ciaa254

28. Izu A, Cohen T, Degruttola V. Bayesian estimation of mixture models with prespecified elements to compare drug resistance in treatment-naïve and experienced tuberculosis cases. PLoS Comput Biol. 2013;9(3):e1002973. doi:10.1371/journal.pcbi.1002973

29. Wan L, Guo Q, Wei JH, et al. Accuracy of a reverse dot blot hybridization assay for simultaneous detection of the resistance of four anti-tuberculosis drugs in Mycobacterium tuberculosis isolated from China. Infect Dis Poverty. 2020;9(1):38. doi:10.1186/s40249-020-00652-z

30. Ofori-Anyinam B, Hamblin M, Coldren ML, et al. Catalase activity deficiency sensitizes multidrug-resistant Mycobacterium tuberculosis to the ATP synthase inhibitor bedaquiline. Nat Commun. 2024;15(1):9792. doi:10.1038/s41467-024-53933-8

31. Nieto R LM, Mehaffy C, Creissen E, et al. Virulence of Mycobacterium tuberculosis after acquisition of isoniazid resistance: individual nature of katG mutants and the possible role of AhpC. PLoS One. 2016;11(11):e0166807. doi:10.1371/journal.pone.0166807

32. Vijay S, Bao NLH, Vinh DN, et al. Rifampicin tolerance and growth fitness among isoniazid-resistant clinical Mycobacterium tuberculosis isolates from a longitudinal study. Elife. 2024;13:RP93243. doi:10.7554/eLife.93243

33. Wilson TM, Collins DM. ahpC, a gene involved in isoniazid resistance of the Mycobacterium tuberculosis complex. Mol Microbiol. 1996;19(5):1025–1034. doi:10.1046/j.1365-2958.1996.449980.x

34. Farhat MR, Sultana R, Iartchouk O, et al. Genetic determinants of drug resistance in Mycobacterium tuberculosis and their diagnostic value. Am J Respir Crit Care Med. 2016;194(5):621–630. doi:10.1164/rccm.201510-2091OC

35. Miotto P, Tessema B, Tagliani E, et al. A standardised method for interpreting the association between mutations and phenotypic drug resistance in Mycobacterium tuberculosis. Eur Respir J. 2017;50(6):1701354. doi:10.1183/13993003.01354-2017

36. Zhang H, Li D, Zhao L, et al. Genome sequencing of 161 Mycobacterium tuberculosis isolates from China identifies genes and intergenic regions associated with drug resistance. Nat Genet. 2013;45(10):1255–1260. doi:10.1038/ng.2735

37. Sherman DR, Mdluli K, Hickey MJ, Barry CE, Stover CK. AhpC, oxidative stress and drug resistance in Mycobacterium tuberculosis. Biofactors. 1999;10(2–3):211–217. doi:10.1002/biof.5520100219