Whole-Transcriptome and -Genome Analysis of Extensively Drug-Resistant Mycobacterium tuberculosis Clinical Isolates Identifies Downregulation of ethA as a Mechanism of Ethionamide Resistance

A quick, easy and efficient protocol for extracting high-quality RNA from Mycobacterium tuberculosis using a spin column commercial kit

RNA extraction from Mycobacterium tuberculosis has been a historically challenging task for researchers due to the thick lipids associated with the cell wall of this "notorious" pathogen that is responsible for Tuberculosis (TB) outbreaks. Several studies have successfully extracted RNA from M. tuberculosis using a Trizol reagent combined with organic solvents. Recently, our laboratory has successfully extracted high quality total RNA using a commercial kit from clinical strains belonging to F15/LAM4/KZN, Beijing and F11 strain families and H37Rv laboratory strain by exploiting high speed homogenizer for cell lysis and spin columns for RNA purification. The quality and integrity of the extracted RNA was analyzed and confirmed through the Nanodrop, Bioanalyzer and RNA 3-(N-morpholino) propanesulfonic acid (MOPS) gel electrophoresis. Furthermore, to confirm the integrity of small RNA (sRNA) molecules due to their vulnerability to degradation, the RNA samples were converted to cDNA and sRNAs were amplified and confirmed through PCR. This detailed RNA extraction protocol proposes to carve a new path into TB transcriptome research without the use of organic solvent for downstream purification steps while yielding high quality RNA that can be used to understand M. tuberculosis transcriptome regulation.

Whole-genome sequencing and transcriptome-characterized in vitro evolution of aminoglycoside resistance in Mycobacterium tuberculosis

Antibiotic resistance of Mycobacterium tuberculosis (Mtb) is a major public health concern worldwide. Therefore, it is of great significance to characterize the mutational pathways by which susceptible Mtb evolves into drug resistance. In this study, we used laboratory evolution to explore the mutational pathways of aminoglycoside resistance. The level of resistance in amikacin inducing Mtb was also associated with changes in susceptibility to other anti-tuberculosis drugs such as isoniazid, levofloxacin and capreomycin. Whole-genome sequencing (WGS) revealed that the induced resistant Mtb strains had accumulated diverse mutations. We found that rrs A1401G was the predominant mutation in aminoglycoside-resistant clinical Mtb isolates from Guangdong. In addition, this study provided global insight into the characteristics of the transcriptome in four representative induced strains and revealed that rrs mutated and unmutated aminoglycoside-resistant Mtb strains have different transcriptional profiles. WGS analysis and transcriptional profiling of Mtb strains during evolution revealed that Mtb strains harbouring rrs A1401G have an evolutionary advantage over other drug-resistant strains under the pressure of aminoglycosides because of their ultra-high resistance level and low physiological impact on the strain. The results of this study should advance our understanding of aminoglycoside resistance mechanisms.

Host determinants of TB disease as targets for HDTs development: A study protocol

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), remains a major global health challenge. New diagnostic and therapeutic strategies are required to curb TB transmission. Here we outline a protocol to validate inflammatory proteins as potential biomarkers of TB disease and to evaluate the candidate genes as potential targets for host-directed therapy (HDT) development. Blood will be isolated from healthy, latent TB infected (LTBI) individuals and TB patients, and expression profiles of genes of interest will be determined using qPCR. A human monocytic cell line will be utilized to knock down genes of interest and to evaluate their contribution to Mtb infection. Pharmaceutical interception of target genes will be performed in peripheral blood mononuclear cells (PBMCs) infected with Mtb. This work will result in identification of TB associated inflammatory markers that can also be targeted for HDT development.

Omics analysis of Mycobacterium tuberculosis isolates uncovers Rv3094c, an ethionamide metabolism-associated gene

Global control of the tuberculosis epidemic is threatened by increasing prevalence of drug resistant M. tuberculosis isolates. Many genome-wide studies focus on SNP-associated drug resistance mechanisms, but drug resistance in 5-30% of M. tuberculosis isolates (varying with antibiotic) appears unrelated to reported SNPs, and alternative drug resistance mechanisms involving variation in gene/protein expression are not well-studied. Here, using an omics approach, we identify 388 genes with lineage-related differential expression and 68 candidate drug resistance-associated gene pairs/clusters in 11 M. tuberculosis isolates (variable lineage/drug resistance profiles). Structural, mutagenesis, biochemical and bioinformatic studies on Rv3094c from the Rv3093c-Rv3095 gene cluster, a gene cluster selected for further investigation as it contains a putative monooxygenase/repressor pair and is associated with ethionamide resistance, provide insights on its involvement in ethionamide sulfoxidation, the initial step in its activation. Analysis of the structure of Rv3094c and its complex with ethionamide and flavin mononucleotide, to the best of our knowledge the first structures of an enzyme involved in ethionamide activation, identify key residues in the flavin mononucleotide and ethionamide binding pockets of Rv3094c, and F221, a gate between flavin mononucleotide and ethionamide allowing their interaction to complete the sulfoxidation reaction. Our work broadens understanding of both lineage- and drug resistance-associated gene/protein expression perturbations and identifies another player in mycobacterial ethionamide metabolism.

Airway microecology in rifampicin-resistant and rifampicin-sensitive pulmonary tuberculosis patients

BACKGROUND

Pulmonary tuberculosis is a chronic infectious disease of the respiratory system. It is still one of the leading causes of death from a single infectious disease, but it has been stuck in the study of a single pathogen. Recent studies have shown that many diseases are associated with disruption of the native microbiota. In this study we investigated the occurrence of tuberculosis and the correlation between drug resistance and respiratory flora. High-throughput 16 S rRNA gene sequencing was used to characterize the respiratory microbiota composition of 30 tuberculosis (TB) affected patients and compared with 30 healthy (H) controls. According to their Gene Xpert results, 30 pulmonary tuberculosis patients were divided into 12 persons in the drug-sensitive group (DS0) and 18 persons in the drug-resistant group (DR0). The microbial flora of the two were compared with the H group.

RESULTS

The data generated by sequencing showed that Firmicutes, Proteus, Bacteroides, Actinomyces and Fusobacterium were the five main bacterial phyla detected, and they constituted more than 96% of the microbial community. The relative abundances of Fusobacterium, Haemophilus, Porphyromonas, Neisseria, TM7, Spirochetes, SR1, and Tenericutes in the TB group was lower than that of the H group, and Granulicatella was higher than the H group. The PcoA diagrams of the two groups had obvious clustering differences. The Alpha diversity of the TB group was lower than that of the H group, and the Beta diversity was higher than that of the H group (P < 0.05). The relative abundance of Streptococcus in the DS0 group was significantly higher than that in the DR0 group (P < 0.05).

CONCLUSION

Pulmonary tuberculosis can cause disorders of the respiratory tract microbial flora, in which the relative abundance of Streptococcus was significantly different between rifampicin-sensitive and rifampicin-resistant patients.

Variants associated with Bedaquiline (BDQ) resistance identified in Rv0678 and efflux pump genes in Mycobacterium tuberculosis isolates from BDQ naïve TB patients in Pakistan

Background

Mutations in the Rv0678, pepQ and atpE genes of Mycobacterium tuberculosis (MTB) have been reported to be associated with reduced antimycobacterial susceptibility to bedaquiline (BDQ). Resistance conferring mutations in treatment naïve MTB strains is likely to have implications for BDQ based new drug regimen that aim to shorten treatment duration. We therefore investigated the genetic basis of resistance to BDQ in MTB clinical isolates from BDQ naïve TB patients from Pakistan. In addition, mutations in genes associated with efflux pumps were investigated as an alternate mechanism of resistance.

Methods

Based on convenience sampling, we studied 48 MTB clinical isolates from BDQ naïve TB patients. These isolates (from our strain bank) included 38 MDR/pre-XDR/XDR (10 BDQ resistant, 8 BDQ intermediate and 20 BDQ susceptible) and 10 pan drug susceptible MTB isolates. All strains were subjected to whole genome sequencing and genomes were analysed to identify variants in Rv0678, pepQ, atpE, Rv1979c, mmpLS and mmpL5 and drug resistance associated efflux pump genes.

Results

Of the BDQ resistant and intermediate strains 44% (8/18) had variants in Rv0678 including; two reported mutations S63R/G, six previously unreported variants; L40F, R50Q and R107C and three frameshift mutations; G25fs, D64fs and D109fs.

Variants in efflux pumps; Rv1273c (G462K), Rv0507c (R426H) and Rv1634c (E198R) were found to be present in drug resistant isolates including BDQ resistant and intermediate isolates. E198R in efflux pump gene Rv1634c was the most frequently occurring variant in BDQ resistant and intermediate isolates (n = 10).

Conclusion

We found RAVs in Rv0678 to be commonly associated with BDQ resistance. Further confirmation of the role of variants in efflux pump genes in resistance is required so that they may be incorporated in genome-based diagnostics for drug resistant MTB.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-022-02475-4.

Alterations in molecular response of Mycobacterium tuberculosis against anti-tuberculosis drugs

BACKGROUND

Tuberculosis (TB), an infectious disease caused by Mycobacterium tuberculosis, has plagued humans since the early middle-ages. More than one million deaths are recorded annually due to TB, even in present times. These deaths are primarily attributed to the constant appearance of resistant TB strains. Even with the advent of new therapeutics and diagnostics techniques, tuberculosis remains challenging to control due to resistant M. tuberculosis strains. Aided by various molecular changes, these strains adapt to stress created by anti-tuberculosis drugs.

MATERIALS AND METHODS

The review thus is an overview of ongoing research in the genome and transcriptome of antibiotic-resistant TB. It explores omics-based research to identify mutation and utilization of differential gene expression.

CONCLUSIONS

This study shows several mutations distinctive in the first- and second-line drug-resistant M. tuberculosis strains. It also explores the expressional differences of genes involved in the fundamental process of the cells and how they help in drug resistance. With the development of transcriptomics-based studies, a new insight has developed to inquire about gene expression changes in drug resistance. This information on expressional pattern changes can be utilized to design the basic platform of anti-TB treatments and therapeutic approaches. These novel insights can be instrumental in disease diagnosis and global containment of resistant TB.

Molecular Determinants of Ethionamide Resistance in Clinical Isolates of Mycobacterium tuberculosis

Background: Ethionamide and prothionamide are now included in group C of the WHO recommended drugs for the treatment of tuberculosis resistant to rifampicin and multidrug-resistant tuberculosis. The clinical relevance of ethionamide and prothionamide has increased with the wide spread of resistant tuberculosis. Methods: We retrospectively analyzed 349 clinical isolates obtained between 2016 and 2020. The susceptibility to ethionamide was tested using both the Bactec^TM MGIT^TM 960 system and the Sensititre^TM MYCOTB plate. Results: The MIC of ethionamide increases with the total resistance of the isolates in a row from susceptible to XDR strains. A significant part of the isolates have a MIC below the breakpoint: 25%, 36%, and 50% for XDR, pre-XDR, and MDR strains. Sensitivity and specificity of detection of mutations were 96% and 86% using MGIT resistance as a reference. Conclusions: Phenotypic methods for testing ethionamide are imperfectly correlated, and the isolates with MIC of 5 mg/L have the intermediate resistance. A significant proportion of resistant TB cases are susceptible and eligible for ethionamide treatment. Resistance could be explained using only analysis of loci ethA, P_fabG1, and inhA for most isolates in the Moscow region. The promoter mutation P_fabG1 c(-15)t predicts resistance to ethionamide with high specificity but low sensitivity.

Identification of Hub Genes Associated With Tuberculous Pleurisy by Integrated Bioinformatics Analysis

Improving the understanding of the molecular mechanism of tuberculous pleurisy is required to develop diagnosis and new therapy strategies of targeted genes. The purpose of this study is to identify important genes related to tuberculous pleurisy. In this study, the expression profile obtained by sequencing the surgically resected pleural tissue was used to explore the differentially co-expressed genes between tuberculous pleurisy tissue and normal tissue. 29 differentially co-expressed genes were screened by weighted gene co-expression network analysis (WGCNA) and differential gene expression analysis methods. According to the functional annotation analysis of R clusterProfiler software package, these genes are mainly enriched in nucleotide−sugar biosynthetic process (biological process), ficolin−1−rich granule lumen (cell component), and electron transfer activity (molecular function). In addition, in the protein-protein interaction (PPI) network, 20 hub genes of DEGs and WCGNA genes were identified using the CytoHubba plug-in of Cytoscape. In the end, RPL17 was identified as a gene that can be the biomarker of tuberculous pleurisy. At the same time, there are seven genes that may have relationship with the disease (UBA7, NDUFB8, UQCRFS1, JUNB, PSMC4, PHPT1, and MAPK11).

Transcriptomic Profile of Mycobacterium smegmatis in Response to an Imidazo[1,2-b][1,2,4,5]tetrazine Reveals Its Possible Impact on Iron Metabolism

Tuberculosis (TB), caused by the Mycobacterium tuberculosis complex bacteria, is one of the most pressing health problems. The development of new drugs and new therapeutic regimens effective against the pathogen is one of the greatest challenges in the way of tuberculosis control. Imidazo[1,2-b][1,2,4,5]tetrazines have shown promising activity against M. tuberculosis and M. smegmatis strains. Mutations in MSMEG_1380 lead to mmpS5–mmpL5 operon overexpression, which provides M. smegmatis with efflux-mediated resistance to imidazo[1,2-b][1,2,4,5]tetrazines, but the exact mechanism of action of these compounds remains unknown. To assess the mode of action of imidazo[1,2-b][1,2,4,5]tetrazines, we analyzed the transcriptomic response of M. smegmatis to three different concentrations of 3a compound: 1/8×, 1/4×, and 1/2× MIC. Six groups of genes responsible for siderophore synthesis and transport were upregulated in a dose-dependent manner, while virtual docking revealed proteins involved in siderophore synthesis as possible targets for 3a. Thus, we suggest that imidazo[1,2-b][1,2,4,5]tetrazines may affect mycobacterial iron metabolism.

Clustering of genes from microarray data using hierarchical projective adaptive resonance theory: a case study of tuberculosis

We propose the hierarchical Projective Adaptive Resonance Theory (PART) algorithm for classification of gene expression data. This algorithm is realized by combing transposed quasi-supervised PART and unsupervised PART. We develop the corresponding validation statistics for each process and compare it with other clustering algorithms in a case study of tuberculosis (TB). First, we use sample-based transposed quasi-supervised PART to obtain optimal clustering results of samples distinguished by time post-infection and the representative genes for each cluster including up-regulated, down-regulated and stable genes. The up- and down-regulated genes show more than 90% similarity to the result derived from Linear Models for Microarray Data and are verified by weighted k-nearest neighbor model on TB projection. Second, we use gene-based unsupervised PART algorithm to cluster these representative genes where functional enrichment analysis is conducted in each cluster. We further confirm the main immune response of human macrophage-like THP-1 cells against TB within 2 days is type I interferon-mediated innate immunity. This study demonstrates how hierarchical PART algorithm analyzes microarray data. The sample-based quasi-supervised PART extracts representative genes and narrows down the shortlist of disease-relevant genes and gene-based unsupervised PART classifies representative genes that help to interpret immune response against TB.

IRF1 as a potential biomarker in Mycobacterium tuberculosis infection

Pulmonary tuberculosis (PTB) is a major global public health problem. The purpose of this study was to find biomarkers that can be used to diagnose tuberculosis. We used four NCBI GEO data sets to conduct analysis. Among the four data sets, GSE139825 is lung tissue microarray, and GSE83456, GSE19491 and GSE50834 are blood microarray. The differential genes of GSE139825 and GSE83456 were 68 and 226, and intersection genes were 11. Gene ontology (GO) analyses of 11 intersection genes revealed that the changes were mostly enriched in regulation of leucocyte cell-cell adhesion and regulation of T-cell activation. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of DEGs revealed that the host response in TB strongly involves cytokine-cytokine receptor interactions and folate biosynthesis. In order to further narrow the range of biomarkers, we used protein-protein interaction to establish a hub gene network of two data sets and a network of 11 candidate genes. Eventually, IRF1 was selected as a biomarker. As validation, IRF1 levels were shown to be up-regulated in patients with TB relative to healthy controls in data sets GSE19491 and GSE50834. Additionally, IRF1 levels were measured in the new patient samples using ELISA. IRF1 was seen to be significantly up-regulated in patients with TB compared with healthy controls with an AUC of 0.801. These results collectively indicate that IRF1 could serve as a new biomarker for the diagnosis of pulmonary tuberculosis.

Transcriptome analysis and identification of candidate genes involved in glyphosate resistance in the fungus Fusarium verticillioides

Glyphosate is a broad-spectrum herbicide that has been widely used for nonselective weed control in soybean fields. In the present study, RNA-seq of an Fusarium verticillioides isolate exhibiting resistance to 120 mM glyphosate revealed gene expression occurring in the presence of glyphosate and led to the identification and screening of candidate genes. A transcriptome analysis revealed 5,548 and 5,361 differentially expressed genes (DEGs) in the glyphosate resistant (GR) Fusarium verticillioides isolate treated with 45 and 90 mM glyphosate, respectively. The gene ontology (GO) pathways associated with these differentially expressed genes primarily included metabolic process, amine metabolic process, cellular aromatic compound metabolism and stress response. The primary Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathways included biosynthesis of secondary metabolites, carbon metabolism, glycolysis/gluconeogenesis, and nitrogen metabolism. The glyphosate degradation-related gene fv04, which belongs to the 3-isopropylalate dehydratase of the aconitase superfamily, was cloned to generate the prokaryotic expression vector pET-29b-fv04, which could be stably expressed in E. coli and promote the degradation of 52.3% of 500 mg/L glyphosate in 72 h. The results of the present study provide new ideas and insights for the acquisition of glyphosate resistance resources.

Heterologous Expression of ethA and katG in Mycobacterium marinum Enables the Rapid Identification of New Prodrugs Active against Mycobacterium tuberculosis

Screening strategies for antituberculosis compounds using Mycobacterium tuberculosis are time consuming and require biosafety level 3 (BSL3) facilities, which makes the development of high-throughput assays difficult and expensive. Mycobacterium marinum, a close genetic relative of M. tuberculosis, possesses several advantages as a suitable model for tuberculosis drug screening. However, despite the high genetic similarity, there are some obvious differences in susceptibility to some tuberculosis drugs between these two species, especially for the prodrugs ethionamide and isoniazid. In this study, we aimed to improve M. marinum as a model for antituberculosis drug identification by heterologous expression of two common drug activators, EthA and KatG. These two activators were overexpressed in M. marinum, and the strains were tested against ethionamide, isoniazid, and a library of established antimycobacterial compounds from TB Alliance to compare drug susceptibility. Both in vitro and in vivo using zebrafish larvae, these genetically modified M. marinum strains showed significantly higher susceptibility against ethionamide and isoniazid, which require activation by EthA and KatG. More importantly, a strain overexpressing both ethA and katG was potentially more susceptible to approximately 20% of the antituberculosis hit compounds from the TB Alliance library. Most of these compounds were activated by EthA in M. marinum Four of these compounds were selected for further analysis, and three of them showed obvious EthA-dependent activity against M. tuberculosis Overall, our developed M. marinum strains are valuable tools for high-throughput discovery of potential novel antituberculosis prodrugs.

Genetic Variation Putatively Associated with Mycobacterium tuberculosis Resistance to Perchlozone, a New Thiosemicarbazone: Clues from Whole Genome Sequencing and Implications for Treatment of Multidrug-Resistant Tuberculosis

Perchlozone ([PCZ] 4-thioureido-iminomethylpyridinium perchlorate) is a new thiosemicarbazone approved for the treatment of multidrug-resistant tuberculosis (MDR-TB) in Russia and some other countries. The ethA and hadABC mutations may confer PCZ resistance. At the same time, ethA mutations are known to mediate resistance to ethionamide (ETH) and prothionamide (PTH). We aimed to study the genetic variation underlying Mycobacterium tuberculosis resistance to PCZ through whole genome sequencing (WGS) of consecutive isolates recovered during long-term treatment. This prospective study included patients admitted in 2018–2019 to the regional tuberculosis dispensary, Kaliningrad, Russia, whose treatment regimen included PCZ. Multiple M. tuberculosis isolates were recovered during PCZ treatment, and the bacterial DNA was subjected to WGS followed by bioinformatics analysis. We identified mutations in the genes putatively associated with PCZ resistance, ethA, and hadA. The most frequent one was a frameshift ethA 106 GA > G (seven of nine patients) and most of the other mutations were also likely present before PCZ treatment. In one patient, a frameshift mutation ethA 702 CT > C emerged after six months of PCZ treatment. A frequent presence of cross-resistance mutations to PCZ and ETH/PTH should be taken into consideration when PCZ is included in the treatment regimen of MDR-TB patients.

Insertion and deletion evolution reflects antibiotics selection pressure in a Mycobacterium tuberculosis outbreak

In genome evolution, genetic variants are the source of diversity, which natural selection acts upon. Treatment of human tuberculosis (TB) induces a strong selection pressure for the emergence of antibiotic resistance-conferring variants in the infecting Mycobacterium tuberculosis (MTB) strains. MTB evolution in response to treatment has been intensively studied and mainly attributed to point substitutions. However, the frequency and contribution of insertions and deletions (indels) to MTB genome evolution remains poorly understood. Here, we analyzed a multi-drug resistant MTB outbreak for the presence of high-quality indels and substitutions. We find that indels are significantly enriched in genes conferring antibiotic resistance. Furthermore, we show that indels are inherited during the outbreak and follow a molecular clock with an evolutionary rate of 5.37e-9 indels/site/year, which is 23 times lower than the substitution rate. Inherited indels may co-occur with substitutions in genes along related biological pathways; examples are iron storage and resistance to second-line antibiotics. This suggests that epistatic interactions between indels and substitutions affect antibiotic resistance and compensatory evolution in MTB.

Whole Genome and Transcriptome Sequencing of Two Multi-Drug Resistant Mycobacterium tuberculosis Strains to Facilitate Illustrating Their Virulence in vivo

Mycobacterium tuberculosis clinical strains usually possess traits different from the laboratory strains like H37Rv, especially those clinical drug resistant strains. With whole genome and transcriptome sequencing, we depicted the feature of two multi-drug resistant Mtb strains in resistance and virulence. Compared with H37Rv, the differential expressed genes (DEGs) of the MDR strains showed featured enrichment in arginine biosynthesis, fatty acid biosynthesis, and metabolism pathway. In the subset of virulence genes, the overlapping DEGs of the MDR strains exhibited downregulation of the cluster in type VII secretion system. In the mice experiment, the MDR strains showed attenuated but distinct virulence, both in survival rate and pathology. Taken together, the whole genome and transcriptome analysis could help understand the unique feature of the MDR strains both in resistance and virulence.

Differential Expression of Resistant and Efflux Pump Genes in MDR-TB Isolates

Background:

Numerous investigations demonstrate efflux as a worldwide bacterial mode of action which contributes to the resistance of drugs. The activity of antibiotics, which subjects to efflux, can be improved by the combined usage of efflux inhibitors. However, the efflux role to the overall levels of antibiotic resistance of clinical M. tuberculosis isolates is inadequately comprehended and is still disregarded by many.

Method:

Here, we assessed the contribution of resistant genes associated with isoniazid (INH) and rifampin (R) resistance to the levels of drug resistance in the (27) clinical isolates of MDR-TB. Additionally, the role of the resistance for six putative drug efflux pump genes to the antibiotics was investigated. The level of katG expression was down-regulated in 24/27 (88.88%) of MDR-TB isolates. Of the 27 MDR-TB isolates, inhA, oxyR-ahpC, and rpoB showed either overexpression or up-regulation in 8 (29.62%), 4 (14.81 %), and 24 (88.88%), respectively. Moreover, the efflux pump genes drrA, drrB, efpA, Rv2459, Rv1634, and Rv1250 were overexpressed under INH/RIF plus fresh pomegranate juice (FPJ) stress signifying the efflux pumps contribution to the overall levels of the resistance of MDR-TB isolates.

Conclusion:

These results displayed that the levels of drug resistance of MDR-TB clinical isolates are due to combination among drug efflux pump and the presence of mutations in target genes, a truth which is often ignored by the specialists of tuberculosis in favour of the almost undoubted significance of drug target- gene mutations for the resistance in M. tuberculosis.

Bacterial Genome-Wide Association Identifies Novel Factors That Contribute to Ethionamide and Prothionamide Susceptibility in Mycobacterium tuberculosis

In Mycobacterium tuberculosis, recent genome-wide association studies have identified a novel constellation of mutations that are correlated with high-level drug resistances. Interpreting the functional importance of the new resistance-associated mutations has been complicated, however, by a lack of experimental validation and a poor understanding of the epistatic factors influencing these correlations, including strain background and programmatic variation in treatment regimens. Here we perform a genome-wide association analysis in a panel of Mycobacterium tuberculosis strains from China to identify variants correlated with resistance to the second-line prodrug ethionamide (ETH). Mutations in a bacterial monooxygenase, Rv0565c, are significantly associated with ETH resistance. We demonstrate that Rv0565c is a novel activator of ETH, independent of the two known activators, EthA and MymA. Clinically prevalent mutations abrogate Rv0565c function, and deletion of Rv0565c confers a consistent fitness benefit on M. tuberculosis in the presence of partially inhibitory doses of ETH. Interestingly, Rv0565c activity affects susceptibility to prothionamide (PTH), the ETH analog used in China, to a greater degree. Further, clinical isolates vary in their susceptibility to both ETH and PTH, to an extent that correlates with the total expression of ETH/PTH activators (EthA, MymA, and Rv0565c). These results suggest that clinical strains considered susceptible to ETH/PTH are not equally fit during treatment due to both Rv0565c mutations and more global variation in the expression of the prodrug activators.IMPORTANCE Phenotypic antibiotic susceptibility testing in Mycobacterium tuberculosis is slow and cumbersome. Rapid molecular diagnostics promise to help guide therapy, but such assays rely on complete knowledge of the molecular determinants of altered antibiotic susceptibility. Recent genomic studies of antibiotic-resistant M. tuberculosis have identified several candidate loci beyond those already known to contribute to antibiotic resistance; however, efforts to provide experimental validation have lagged. Our study identifies a gene (Rv0565c) that is associated with resistance to the second-line antibiotic ethionamide at a population level. We then use bacterial genetics to show that the variants found in clinical strains of M. tuberculosis improve bacterial survival after ethionamide exposure.

Genome-Wide Transcriptional Responses of Mycobacterium to Antibiotics

Antibiotics can stimulate or depress gene expression in bacteria. The analysis of transcriptional responses of Mycobacterium to antimycobacterial compounds has improved our understanding of the mode of action of various drug classes and the efficacy and effect of such compounds on the global metabolism of Mycobacterium. This approach can provide new insights for known antibiotics, for example those currently used for tuberculosis treatment, as well as help to identify the mode of action and predict the targets of new compounds identified by whole-cell screening assays. In addition, changes in gene expression profiles after antimycobacterial treatment can provide information about the adaptive ability of bacteria to escape the effects of antibiotics and allow monitoring of the physiology of the bacteria during treatment. Genome-wide expression profiling also makes it possible to pinpoint genes differentially expressed between drug sensitive Mycobacterium and multidrug-resistant clinical isolates. Finally, genes involved in adaptive responses and drug tolerance could become new targets for improving the efficacy of existing antibiotics.