Clinically prevalent mutations in Mycobacterium tuberculosis alter propionate metabolism and mediate multidrug tolerance

Genomic analysis of lineage-specific transmission of multidrug resistance tuberculosis in China

OBJECTIVES

We investigated the genetic diversities and lineage-specific transmission dynamics of multidrug-resistant tuberculosis (MDR-TB), with the goal of determining the potential factors driving the MDR epidemics in China.

METHODS

We curated a large nationwide Mycobacterium tuberculosis (M. tuberculosis) whole genome sequence data set, including 1313 MDR strains. We reconstructed the phylogeny and mapped the transmission networks of MDR-TB across China using Bayesian inference. To identify drug-resistance variants linked to enhanced transmissibility, we employed ordinary least-squares (OLS) regression analysis.

RESULT

The majority of MDR-TB strains in China belong to lineage 2.2.1. Transmission chain analysis has indicated that the repeated and frequent transmission of L2.2.1 plays a central role in the establishment of MDR epidemic in China, but no occurrence of a large predominant MDR outbreak was detected. Using OLS regression, the most common single nucleotide polymorphisms (SNPs) associated with resistance to isoniazid (katG_p.Ser315Thr and katG_p.Ser315Asn) and rifampicin (rpoB_p.Ser450Leu, rpoB_p.His445Tyr, rpoB_p.His445Arg, rpoB_p.His445Asp, and rpoB_p.His445Asn) were more likely to be found in L2 clustered strains. Several putative compensatory mutations in rpoA, rpoC, and katG were significantly associated with clustering. The eastern, central, and southern regions of China had a high level of connectivity for the migration of L2 MDR strains throughout the country. The skyline plot showed distinct population size expansion dynamics for MDR-TB lineages in China.

CONCLUSION

MDR-TB epidemic in China is predominantly driven by the spread of highly transmissible Beijing strains. A range of drug-resistance mutations of L2 MDR-TB strains displayed minimal fitness costs and may facilitate their transmission.

Propionate prevents loss of the PDIM virulence lipid in Mycobacterium tuberculosis

Phthiocerol dimycocerosate (PDIM) is an essential virulence lipid of Mycobacterium tuberculosis. In vitro culturing rapidly selects for spontaneous PDIM-negative mutants that have attenuated virulence and increased cell wall permeability, thus impacting the relevance of experimental findings. PDIM loss can also reduce the efficacy of the BCG Pasteur vaccine. Here we show that vancomycin susceptibility can rapidly screen for M. tuberculosis PDIM production. We find that metabolic deficiency of methylmalonyl-CoA impedes the growth of PDIM-producing bacilli, selecting for PDIM-negative variants. Supplementation with odd-chain fatty acids, cholesterol or vitamin B12 restores PDIM-positive bacterial growth. Specifically, we show that propionate supplementation enhances PDIM-producing bacterial growth and selects against PDIM-negative mutants, analogous to in vivo conditions. Our study provides a simple approach to screen for and maintain PDIM production, and reveals how discrepancies between the host and in vitro nutrient environments can attenuate bacterial pathogenicity.

Identification of bacterial determinants of tuberculosis infection and treatment outcomes: a phenogenomic analysis of clinical strains

BACKGROUND

Bacterial diversity could contribute to the diversity of tuberculosis infection and treatment outcomes observed clinically, but the biological basis of this association is poorly understood. The aim of this study was to identify associations between phenogenomic variation in Mycobacterium tuberculosis and tuberculosis clinical features.

METHODS

We developed a high-throughput platform to define phenotype-genotype relationships in M tuberculosis clinical isolates, which we tested on a set of 158 drug-sensitive M tuberculosis strains sampled from a large tuberculosis clinical study in Ho Chi Minh City, Viet Nam. We tagged the strains with unique genetic barcodes in multiplicate, allowing us to pool the strains for in-vitro competitive fitness assays across 16 host-relevant antibiotic and metabolic conditions. Relative fitness was quantified by deep sequencing, enumerating output barcode read counts relative to input normalised values. We performed a genome-wide association study to identify phylogenetically linked and monogenic mutations associated with the in-vitro fitness phenotypes. These genetic determinants were further associated with relevant clinical outcomes (cavitary disease and treatment failure) by calculating odds ratios (ORs) with binomial logistic regressions. We also assessed the population-level transmission of strains associated with cavitary disease and treatment failure using terminal branch length analysis of the phylogenetic data.

FINDINGS

M tuberculosis clinical strains had diverse growth characteristics in host-like metabolic and drug conditions. These fitness phenotypes were highly heritable, and we identified monogenic and phylogenetically linked variants associated with the fitness phenotypes. These data enabled us to define two genetic features that were associated with clinical outcomes. First, mutations in Rv1339, a phosphodiesterase, which were associated with slow growth in glycerol, were further associated with treatment failure (OR 5·34, 95% CI 1·21-23·58, p=0·027). Second, we identified a phenotypically distinct slow-growing subclade of lineage 1 strains (L1.1.1.1) that was associated with cavitary disease (OR 2·49, 1·11-5·59, p=0·027) and treatment failure (OR 4·76, 1·53-14·78, p=0·0069), and which had shorter terminal branch lengths on the phylogenetic tree, suggesting increased transmission.

INTERPRETATION

Slow growth under various antibiotic and metabolic conditions served as in-vitro intermediate phenotypes underlying the association between M tuberculosis monogenic and phylogenetically linked mutations and outcomes such as cavitary disease, treatment failure, and transmission potential. These data suggest that M tuberculosis growth regulation is an adaptive advantage for bacterial success in human populations, at least in some circumstances. These data further suggest markers for the underlying bacterial processes that contribute to these clinical outcomes.

FUNDING

National Health and Medical Research Council/A∗STAR, National Institutes of Allergy and Infectious Diseases, National Institute of Child Health and Human Development, and the Wellcome Trust Fellowship in Public Health and Tropical Medicine.

Substrate and functional characterization of the lysine acetyltransferase MsKat and deacetylase MsCobB in Mycobacterium smegmatis

Tuberculosis (TB) is a serious cause of infectious death worldwide. Recent studies have reported that about 30% of the Mtb proteome was modified post-translationally, indicating that their functions are essential for drug resistance, mycobacterial survival, and pathogenicity. Among them, lysine acetylation, reversibly regulated by acetyltransferase and deacetylase, has important roles involved in energy metabolism, cellular adaptation, and protein interactions. However, the substrate and biological functions of these two important regulatory enzymes remain unclear. Herein, we utilized the non-pathogenic M. smegmatis strain as a model and systematically investigated the dynamic proteome changes in response to the overexpressing of MsKat/MsCobB in mycobacteria. A total of 4179 proteins and 1236 acetylated sites were identified in our data. Further analysis of the dynamic changes involved in proteome and acetylome showed that MsKat/MsCobB played a regulatory role in various metabolic pathways and nucleic acid processes. After that, the quantitative mass spectrometric method was utilized and proved that the AMP-dependent synthetase, Citrate synthase, ATP-dependent specificity component of the Clp protease, and ATP-dependent DNA/RNA helicases were identified to be the substrates of MsKat. Overall, our study provided an important resource underlying the substrates and functions of the acetylation regulatory enzymes in mycobacteria. SIGNIFICANCE: In this study, we systematically analyzed the dynamic molecular changes in response to the MsKat/MsCobB overexpression in mycobacteria at proteome and lysine acetylation level by using a TMT-based quantitative proteomic approach. Pathways related with glycolysis, degradation of branched chain amino acids, phosphotransferase system were affected after disturbance of the two regulates enzymes involved in lysine acetylation. We also proved that AMP-dependent synthetase Clp protease, ATP-dependent DNA/RNA helicases and citrate synthase was the substrate of MsKat according to our proteomic data and biological validation. Together, our study underlined the substrates and functions of the acetylation regulatory enzymes in mycobacteria.

Rifampicin tolerance and growth fitness among isoniazid-resistant clinical Mycobacterium tuberculosis isolates: an in-vitro longitudinal study

Antibiotic tolerance in Mycobacterium tuberculosis leads to less effective bacterial killing, poor treatment responses and resistant emergence. Therefore, we investigated the rifampicin tolerance of M. tuberculosis isolates, with or without pre-existing isoniazid-resistance. We determined the in-vitro rifampicin survival fraction by minimum duration of killing assay in isoniazid susceptible (IS, n=119) and resistant (IR, n=84) M. tuberculosis isolates. Then we correlated the rifampicin tolerance with bacterial growth, rifampicin minimum inhibitory concentrations (MICs) and isoniazid-resistant mutations. The longitudinal IR isolates collected from patients were analyzed for changes in rifampicin tolerance and associated emergence of genetic variants. The median duration of rifampicin exposure reducing the M. tuberculosis surviving fraction by 90% (minimum duration of killing-MDK90) increased from 1.23 (95%CI 1.11; 1.37) and 1.31 (95%CI 1.14; 1.48) to 2.55 (95%CI 2.04; 2.97) and 1.98 (95%CI 1.69; 2.56) days, for IS and IR respectively, during 15 to 60 days of incubation. This indicated the presence of fast and slow growing tolerant sub-populations. A range of 6 log 10 -fold survival fraction enabled classification of tolerance as low, medium or high and revealed IR association with increased tolerance with faster growth (OR=2.68 for low vs. medium, OR=4.42 for low vs. high, P -trend=0.0003). The high tolerance in IR isolates was specific to those collected during rifampicin treatment in patients and associated with bacterial genetic microvariants. Furthermore, the high rifampicin tolerant IR isolates have survival potential similar to multi-drug resistant isolates. These findings suggest that IR tuberculosis needs to be evaluated for high rifampicin tolerance to improve treatment regimen and prevent the risk of MDR-TB emergence.

Sub-MIC vancomycin enhances the antibiotic tolerance of vancomycin-intermediate Staphylococcus aureus through downregulation of protein succinylation

Bacteria develop tolerance after transient exposure to antibiotics, and tolerance is a significant driver of resistance. The purpose of this study is to evaluate the mechanisms underlying tolerance formation in vancomycin-intermediate Staphylococcus aureus (VISA) strains. VISA strains were cultured with sub-minimum inhibitory concentrations (sub-MICs) of vancomycin. Enhanced vancomycin tolerance was observed in VISA strains with distinct genetic lineages. Western blot revealed that the VISA protein succinylation (Ksucc) levels decreased with the increase in vancomycin exposure. Importantly, Ksucc modification, vancomycin tolerance, and cell wall synthesis were simultaneously affected after deletion of SacobB, which encodes a desuccinylase in S. aureus. Several Ksucc sites were identified in MurA, and vancomycin MIC levels of murA mutant and Ksucc-simulated (MurA(K69E) and MurA(K191E)) mutants were reduced. The vancomycin MIC levels of K65-MurA(K191E) in particular decreased to 1 mg/L, converting VISA strain K65 to a vancomycin-susceptible S. aureus strain. We further demonstrated that the enzymatic activity of MurA was dependent on Ksucc modification. Our data suggested the influence of vancomycin exposure on bacterial tolerance, and protein Ksucc modification is a novel mechanism in regulating vancomycin tolerance.

Compensatory evolution in NusG improves fitness of drug-resistant M. tuberculosis

Drug-resistant bacteria are emerging as a global threat, despite frequently being less fit than their drug-susceptible ancestors1-8. Here we sought to define the mechanisms that drive or buffer the fitness cost of rifampicin resistance (RifR) in the bacterial pathogen Mycobacterium tuberculosis (Mtb). Rifampicin inhibits RNA polymerase (RNAP) and is a cornerstone of modern short-course tuberculosis therapy9,10. However, RifR Mtb accounts for one-quarter of all deaths due to drug-resistant bacteria11,12. We took a comparative functional genomics approach to define processes that are differentially vulnerable to CRISPR interference (CRISPRi) inhibition in RifR Mtb. Among other hits, we found that the universally conserved transcription factor NusG is crucial for the fitness of RifR Mtb. In contrast to its role in Escherichia coli, Mtb NusG has an essential RNAP pro-pausing function mediated by distinct contacts with RNAP and the DNA13. We find this pro-pausing NusG-RNAP interface to be under positive selection in clinical RifR Mtb isolates. Mutations in the NusG-RNAP interface reduce pro-pausing activity and increase fitness of RifR Mtb. Collectively, these results define excessive RNAP pausing as a molecular mechanism that drives the fitness cost of RifR in Mtb, identify a new mechanism of compensation to overcome this cost, suggest rational approaches to exacerbate the fitness cost, and, more broadly, could inform new therapeutic approaches to develop drug combinations to slow the evolution of RifR in Mtb.

Drug-resistant tuberculosis: a persistent global health concern

Drug-resistant tuberculosis (TB) is estimated to cause 13% of all antimicrobial resistance-attributable deaths worldwide and is driven by both ongoing resistance acquisition and person-to-person transmission. Poor outcomes are exacerbated by late diagnosis and inadequate access to effective treatment. Advances in rapid molecular testing have recently improved the diagnosis of TB and drug resistance. Next-generation sequencing of Mycobacterium tuberculosis has increased our understanding of genetic resistance mechanisms and can now detect mutations associated with resistance phenotypes. All-oral, shorter drug regimens that can achieve high cure rates of drug-resistant TB within 6-9 months are now available and recommended but have yet to be scaled to global clinical use. Promising regimens for the prevention of drug-resistant TB among high-risk contacts are supported by early clinical trial data but final results are pending. A person-centred approach is crucial in managing drug-resistant TB to reduce the risk of poor treatment outcomes, side effects, stigma and mental health burden associated with the diagnosis. In this Review, we describe current surveillance of drug-resistant TB and the causes, risk factors and determinants of drug resistance as well as the stigma and mental health considerations associated with it. We discuss recent advances in diagnostics and drug-susceptibility testing and outline the progress in developing better treatment and preventive therapies.

Iron-related gene mutations driving global Mycobacterium tuberculosis transmission revealed by whole-genome sequencing

BACKGROUND

Iron plays a crucial role in the growth of Mycobacterium tuberculosis (M. tuberculosis). However, the precise regulatory mechanism governing this system requires further elucidation. Additionally, limited studies have examined the impact of gene mutations related to iron on the transmission of M. tuberculosis globally. This research aims to investigate the correlation between mutations in iron-related genes and the worldwide transmission of M. tuberculosis.

RESULTS

A total of 13,532 isolates of M. tuberculosis were included in this study. Among them, 6,104 (45.11%) were identified as genomic clustered isolates, while 8,395 (62.04%) were classified as genomic clade isolates. Our results showed that a total of 12 single nucleotide polymorphisms (SNPs) showed a positive correlation with clustering, such as Rv1469 (ctpD, C758T), Rv3703c (etgB, G1122T), and Rv3743c (ctpJ, G676C). Additionally, seven SNPs, including Rv0104 (T167G, T478G), Rv0211 (pckA, A302C), Rv0283 (eccB3, C423T), Rv1436 (gap, G654T), ctpD C758T, and etgB C578A, demonstrated a positive correlation with transmission clades across different countries. Notably, our findings highlighted the positive association of Rv0104 T167G, pckA A302C, eccB3 C423T, ctpD C758T, and etgB C578A with transmission clades across diverse regions. Furthermore, our analysis identified 78 SNPs that exhibited significant associations with clade size.

CONCLUSIONS

Our study reveals the link between iron-related gene SNPs and M. tuberculosis transmission, offering insights into crucial factors influencing the pathogenicity of the disease. This research holds promise for targeted strategies in prevention and treatment, advancing research and interventions in this field.

Ongoing evolution of the Mycobacterium tuberculosis lactate dehydrogenase reveals the pleiotropic effects of bacterial adaption to host pressure

The bacterial determinants that facilitate Mycobacterium tuberculosis (Mtb) adaptation to the human host environment are poorly characterized. We have sought to decipher the pressures facing the bacterium in vivo by assessing Mtb genes that are under positive selection in clinical isolates. One of the strongest targets of selection in the Mtb genome is lldD2, which encodes a quinone-dependent L-lactate dehydrogenase (LldD2) that catalyzes the oxidation of lactate to pyruvate. Lactate accumulation is a salient feature of the intracellular environment during infection and lldD2 is essential for Mtb growth in macrophages. We determined the extent of lldD2 variation across a set of global clinical isolates and defined how prevalent mutations modulate Mtb fitness. We show the stepwise nature of lldD2 evolution that occurs as a result of ongoing lldD2 selection in the background of ancestral lineage-defining mutations and demonstrate that the genetic evolution of lldD2 additively augments Mtb growth in lactate. Using quinone-dependent antibiotic susceptibility as a functional reporter, we also find that the evolved lldD2 mutations functionally increase the quinone-dependent activity of LldD2. Using 13C-lactate metabolic flux tracing, we find that lldD2 is necessary for robust incorporation of lactate into central carbon metabolism. In the absence of lldD2, label preferentially accumulates in dihydroxyacetone phosphate (DHAP) and glyceraldehyde-3-phosphate (G3P) and is associated with a discernible growth defect, providing experimental evidence for accrued lactate toxicity via the deleterious buildup of sugar phosphates. The evolved lldD2 variants increase lactate incorporation to pyruvate while altering triose phosphate flux, suggesting both an anaplerotic and detoxification benefit to lldD2 evolution. We further show that the mycobacterial cell is transcriptionally sensitive to the changes associated with altered lldD2 activity which affect the expression of genes involved in cell wall lipid metabolism and the ESX- 1 virulence system. Together, these data illustrate a multifunctional role of LldD2 that provides context for the selective advantage of lldD2 mutations in adapting to host stress.

Polyketide synthases mutation in tuberculosis transmission revealed by whole genomic sequence, China, 2011-2019

Introduction: Tuberculosis (TB) is an infectious disease caused by a bacterium called Mycobacterium tuberculosis (Mtb). Previous studies have primarily focused on the transmissibility of multidrug-resistant (MDR) or extensively drug-resistant (XDR) Mtb. However, variations in virulence across Mtb lineages may also account for differences in transmissibility. In Mtb, polyketide synthase (PKS) genes encode large multifunctional proteins which have been shown to be major mycobacterial virulence factors. Therefore, this study aimed to identify the role of PKS mutations in TB transmission and assess its risk and characteristics. Methods: Whole genome sequences (WGSs) data from 3,204 Mtb isolates was collected from 2011 to 2019 in China. Whole genome single nucleotide polymorphism (SNP) profiles were used for phylogenetic tree analysis. Putative transmission clusters (≤10 SNPs) were identified. To identify the role of PKS mutations in TB transmission, we compared SNPs in the PKS gene region between "clustered isolates" and "non-clustered isolates" in different lineages. Results: Cluster-associated mutations in ppsA, pks12, and pks13 were identified among different lineage isolates. They were statistically significant among clustered strains, indicating that they may enhance the transmissibility of Mtb. Conclusion: Overall, this study provides new insights into the function of PKS and its localization in M. tuberculosis. The study found that ppsA, pks12, and pks13 may contribute to disease progression and higher transmission of certain strains. We also discussed the prospective use of mutant ppsA, pks12, and pks13 genes as drug targets.

Retinoic acid signaling in development and differentiation commitment and its regulatory topology

Retinoic acid (RA), the derivative of vitamin A/retinol, is a signaling molecule with important implications in health and disease. It is a well-known developmental morphogen that functions mainly through the transcriptional activity of nuclear RA receptors (RARs) and, uncommonly, through other nuclear receptors, including peroxisome proliferator-activated receptors. Intracellular RA is under spatiotemporally fine-tuned regulation by synthesis and degradation processes catalyzed by retinaldehyde dehydrogenases and P450 family enzymes, respectively. In addition to dictating the transcription architecture, RA also impinges on cell functioning through non-genomic mechanisms independent of RAR transcriptional activity. Although RA-based differentiation therapy has achieved impressive success in the treatment of hematologic malignancies, RA also has pro-tumor activity. Here, we highlight the relevance of RA signaling in cell-fate determination, neurogenesis, visual function, inflammatory responses and gametogenesis commitment. Genetic and post-translational modifications of RAR are also discussed. A better understanding of RA signaling will foster the development of precision medicine to improve the defects caused by deregulated RA signaling.

Association between fatty acid metabolism gene mutations and Mycobacterium tuberculosis transmission revealed by whole genome sequencing

BACKGROUND

Fatty acid metabolism greatly promotes the virulence and pathogenicity of Mycobacterium tuberculosis (M.tb). However, the regulatory mechanism of fatty acid metabolism in M.tb remains to be elucidated, and limited evidence about the effects of gene mutations in fatty acid metabolism on the transmission of M.tb was reported.

RESULTS

Overall, a total of 3193 M.tb isolates were included in the study, of which 1596 (50%) were genomic clustered isolates. Most of the tuberculosis isolates belonged to lineage2(n = 2744,85.93%), followed by lineage4(n = 439,13.75%) and lineage3(n = 10,0.31%).Regression results showed that the mutations of gca (136,605, 317G > C, Arg106Pro; OR, 22.144; 95% CI, 2.591-189.272), ogt(1,477,346, 286G > C ,Gly96Arg; OR, 3.893; 95%CI, 1.432-10.583), and rpsA (1,834,776, 1235 C > T, Ala412Val; OR, 3.674; 95% CI, 1.217-11.091) were significantly associated with clustering; mutations in gca and rpsA were also significantly associated with clustering of lineage2. Mutation in arsA(3,001,498, 885 C > G, Thr295Thr; OR, 6.278; 95% CI, 2.508-15.711) was significantly associated with cross-regional clusters. We also found that 20 mutation sites were positively correlated with cluster size, while 11 fatty acid mutation sites were negatively correlated with cluster size.

CONCLUSION

Our research results suggested that mutations in genes related to fatty acid metabolism were related to the transmission of M.tb. This research could help in the future control of the transmission of M.tb.

Association between two-component systems gene mutation and Mycobacterium tuberculosis transmission revealed by whole genome sequencing

BACKGROUND

Two-component systems (TCSs) assume a pivotal function in Mycobacterium tuberculosis (M.tuberculosis) growth. However, the exact regulatory mechanism of this system needs to be elucidated, and only a few studies have investigated the effect of gene mutations within TCSs on M.tuberculosis transmission. This research explored the relationship between TCSs gene mutation and the global transmission of (M.tuberculosis).

RESULTS

A total of 13531 M.tuberculosis strains were enrolled in the study. Most of the M.tuberculosis strains belonged to lineage4 (n=6497,48.0%), followed by lineage2 (n=5136,38.0%). Our results showed that a total of 36 single nucleotide polymorphisms (SNPs) were positively correlated with clustering of lineage2, such as Rv0758 (phoR, C820G), Rv1747(T1102C), and Rv1057(C1168T). A total of 30 SNPs showed positive correlation with clustering of lineage4, such as phoR(C182A, C1184G, C662T, T758G), Rv3764c (tcrY, G1151T), and Rv1747 C20T. A total of 19 SNPs were positively correlated with cross-country transmission of lineage2, such as phoR A575C, Rv1028c (kdpD, G383T, G1246C), and Rv1057 G817T. A total of 41 SNPs were positively correlated with cross-country transmission of lineage4, such as phoR(T758G, T327G, C284G), kdpD(G1755A, G625C), Rv1057 C980T, and Rv1747 T373G.

CONCLUSIONS

Our study identified that SNPs in genes of two-component systems were related to the transmission of M. tuberculosis. This finding adds another layer of complexity to M. tuberculosis virulence and provides insight into future research that will help to elucidate a novel mechanism of M. tuberculosis pathogenicity.

Genome-Wide Association Studies (GWAS) Approaches for the Detection of Genetic Variants Associated with Antibiotic Resistance: A Systematic Review

Antibiotic resistance is a significant threat to public health worldwide. Genome-wide association studies (GWAS) have emerged as a powerful tool to identify genetic variants associated with this antibiotic resistance. By analyzing large datasets of bacterial genomes, GWAS can provide valuable insights into the resistance mechanisms and facilitate the discovery of new drug targets. The present study aimed to undertake a systematic review of different GWAS approaches used for detecting genetic variants associated with antibiotic resistance. We comprehensively searched the PubMed and Scopus databases to identify relevant studies published from 2013 to February 2023. A total of 40 studies met our inclusion criteria. These studies explored a wide range of bacterial species, antibiotics, and study designs. Notably, most of the studies were centered around human pathogens such as Mycobacterium tuberculosis, Escherichia coli, Neisseria gonorrhoeae, and Staphylococcus aureus. The review seeks to explore the several GWAS approaches utilized to investigate the genetic mechanisms associated with antibiotic resistance. Furthermore, it examines the contributions of GWAS approaches in identifying resistance-associated genetic variants through binary and continuous phenotypes. Overall, GWAS holds great potential to enhance our understanding of bacterial resistance and improve strategies to combat infectious diseases.

Mycobacterial RNase E cleaves with a distinct sequence preference and controls the degradation rates of most Mycolicibacterium smegmatis mRNAs

The mechanisms and regulation of RNA degradation in mycobacteria have been subject to increased interest following the identification of interplay between RNA metabolism and drug resistance. Mycobacteria encode multiple ribonucleases predicted to participate in mRNA degradation and/or processing of stable RNAs. RNase E is hypothesized to play a major role in mRNA degradation because of its essentiality in mycobacteria and its role in mRNA degradation in gram-negative bacteria. Here, we defined the impact of RNase E on mRNA degradation rates transcriptome-wide in the nonpathogenic model Mycolicibacterium smegmatis. RNase E played a rate-limiting role in degradation of the transcripts encoded by at least 89% of protein-coding genes, with leadered transcripts often being more affected by RNase E repression than leaderless transcripts. There was an apparent global slowing of transcription in response to knockdown of RNase E, suggesting that M. smegmatis regulates transcription in responses to changes in mRNA degradation. This compensation was incomplete, as the abundance of most transcripts increased upon RNase E knockdown. We assessed the sequence preferences for cleavage by RNase E transcriptome-wide in M. smegmatis and Mycobacterium tuberculosis and found a consistent bias for cleavage in C-rich regions. Purified RNase E had a clear preference for cleavage immediately upstream of cytidines, distinct from the sequence preferences of RNase E in gram-negative bacteria. We furthermore report a high-resolution map of mRNA cleavage sites in M. tuberculosis, which occur primarily within the RNase E-preferred sequence context, confirming that RNase E has a broad impact on the M. tuberculosis transcriptome.

KatG catalase deficiency confers bedaquiline hyper-susceptibility to isoniazid resistant Mycobacterium tuberculosis

Multidrug-resistant tuberculosis (MDR-TB) is a growing source of global mortality and threatens global control of tuberculosis (TB) disease. The diarylquinoline bedaquiline (BDQ) recently emerged as a highly efficacious drug against MDR-TB, defined as resistance to the first-line drugs isoniazid (INH) and rifampin. INH resistance is primarily caused by loss-of-function mutations in the catalase KatG, but mechanisms underlying BDQ's efficacy against MDR-TB remain unknown. Here we employ a systems biology approach to investigate BDQ hyper-susceptibility in INH-resistant Mycobacterium tuberculosis . We found hyper-susceptibility to BDQ in INH-resistant cells is due to several physiological changes induced by KatG deficiency, including increased susceptibility to reactive oxygen species and DNA damage, remodeling of transcriptional programs, and metabolic repression of folate biosynthesis. We demonstrate BDQ hyper-susceptibility is common in INH-resistant clinical isolates. Collectively, these results highlight how altered bacterial physiology can impact drug efficacy in drug-resistant bacteria.

The PDIM paradox of Mycobacterium tuberculosis: new solutions to a persistent problem

Phthiocerol dimycocerosate (PDIM) is an essential virulence lipid of Mycobacterium tuberculosis. In vitro culturing rapidly selects for spontaneous mutations that cause PDIM loss leading to virulence attenuation and increased cell wall permeability. We discovered that PDIM loss is due to a metabolic deficiency of methylmalonyl-CoA that impedes the growth of PDIM-producing bacilli. This can be remedied by supplementation with odd-chain fatty acids, cholesterol, or vitamin B12. We developed a much-needed facile and scalable routine assay for PDIM production and show that propionate supplementation enhances the growth of PDIM-producing bacilli and selects against PDIM-negative mutants, analogous to in vivo conditions. Our results solve a major issue in tuberculosis research and exemplify how discrepancies between the host and in vitro nutrient environments can attenuate bacterial pathogenicity.

Ongoing evolution of the Mycobacterium tuberculosis lactate dehydrogenase reveals the pleiotropic effects of bacterial adaption to host pressure

The bacterial determinants that facilitate Mycobacterium tuberculosis (Mtb) adaptation to the human host environment are poorly characterized. We have sought to decipher the pressures facing the bacterium in vivo by assessing Mtb genes that are under positive selection in clinical isolates. One of the strongest targets of selection in the Mtb genome is lldD2 , which encodes a quinone-dependent L-lactate dehydrogenase (LldD2) that catalyzes the oxidation of lactate to pyruvate. Lactate accumulation is a salient feature of the intracellular environment during infection and lldD2 is essential for Mtb growth in macrophages. We determined the extent of lldD2 variation across a set of global clinical isolates and defined how prevalent mutations modulates Mtb fitness. We show the stepwise nature of lldD2 evolution that occurs as a result of ongoing lldD2 selection in the background of ancestral lineage defining mutations and demonstrate that the genetic evolution of lldD2 additively augments Mtb growth in lactate. Using quinone-dependent antibiotic susceptibility as a functional reporter, we also find that the evolved lldD2 mutations functionally increase the quinone-dependent activity of LldD2. Using 13 C-lactate metabolic flux tracing, we find that lldD2 is necessary for robust incorporation of lactate into central carbon metabolism. In the absence of lldD2 , label preferentially accumulates in methylglyoxal precursors dihydroxyacetone phosphate (DHAP) and glyceraldehyde-3-phosphate (G3P) and is associated with a discernible growth defect, providing experimental evidence for accumulated lactate toxicity via a methylglyoxal pathway that has been proposed previously. The evolved lldD2 variants increase lactate incorporation to pyruvate but also alter flux in the methylglyoxal pathway, suggesting both an anaplerotic and detoxification benefit to lldD2 evolution. We further show that the mycobacterial cell is transcriptionally sensitive to the changes associated with altered lldD2 activity which affect the expression of genes involved in cell wall lipid metabolism and the ESX-1 virulence system. Together, these data illustrate a multifunctional role of LldD2 that provide context for the selective advantage of lldD2 mutations in adapting to host stress.

A statistical genomics framework to trace bacterial genomic predictors of clinical outcomes in Staphylococcus aureus bacteremia

Outcomes of severe bacterial infections are determined by the interplay between host, pathogen, and treatments. While human genomics has provided insights into host factors impacting Staphylococcus aureus infections, comparatively little is known about S. aureus genotypes and disease severity. Building on the hypothesis that bacterial pathoadaptation is a key outcome driver, we developed a genome-wide association study (GWAS) framework to identify adaptive mutations associated with treatment failure and mortality in S. aureus bacteremia (1,358 episodes). Our research highlights the potential of vancomycin-selected mutations and vancomycin minimum inhibitory concentration (MIC) as key explanatory variables to predict infection severity. The contribution of bacterial variation was much lower for clinical outcomes (heritability <5%); however, GWASs allowed us to identify additional, MIC-independent candidate pathogenesis loci. Using supervised machine learning, we were able to quantify the predictive potential of these adaptive signatures. Our statistical genomics framework provides a powerful means to capture adaptive mutations impacting severe bacterial infections.

The Relative Positioning of Genotyping and Phenotyping for Tuberculosis Resistance Screening in Two EU National Reference Laboratories in 2023

The routine use of whole genome sequencing (WGS) as a reference typing technique for Mycobacterium tuberculosis epidemiology combined with the catalogued and extensive knowledge base of resistance-associated mutations means an initial susceptibility prediction can be derived from all cultured isolates in our laboratories based on WGS data alone. Preliminary work has confirmed, in our low-burden settings, these predictions are for first-line drugs, reproducible, robust with an accuracy similar to phenotypic drug susceptibility testing (pDST) and in many cases able to also predict the level of resistance (MIC). Routine screening for drug resistance by WGS results in approximately 80% of the isolates received being predicted as fully susceptible to the first-line drugs. Parallel testing with both WGS and pDST has demonstrated that routine pDST of genotypically fully susceptible isolates yields minimal additional information. Thus, rather than re-confirming all fully sensitive WGS-based predictions, we suggest that a more efficient use of available mycobacterial culture capacity in our setting is the development of a more extensive and detailed pDST targeted at any mono or multi-drug-resistant isolates identified by WGS screening. Phenotypic susceptibility retains a key role in the determination of an extended susceptibility profile for mono/multi-drugresistant isolates identified by WGS screening. The pDST information collected is also needed to support the development of future catalogues of resistance-associated mutations.

Insight into Population Structure and Drug Resistance of Pediatric Tuberculosis Strains from China and Russia Gained through Whole-Genome Sequencing

This study aimed to determine phenotypic and genotypic drug resistance patterns of Mycobacterium tuberculosis strains from children with tuberculosis (TB) in China and Russia, two high-burden countries for multi/extensively-drug resistant (MDR/XDR) TB. Whole-genome sequencing data of M. tuberculosis isolates from China (n = 137) and Russia (n = 60) were analyzed for phylogenetic markers and drug-resistance mutations, followed by comparison with phenotypic susceptibility data. The Beijing genotype was detected in 126 Chinese and 50 Russian isolates. The Euro-American lineage was detected in 10 Russian and 11 Chinese isolates. In the Russian collection, the Beijing genotype and Beijing B0/W148-cluster were dominated by MDR strains (68% and 94%, respectively). Ninety percent of B0/W148 strains were phenotypically pre-XDR. In the Chinese collection, neither of the Beijing sublineages was associated with MDR/pre-XDR status. MDR was mostly caused by low fitness cost mutations (rpoB S450L, katG S315T, rpsL K43R). Chinese rifampicin-resistant strains demonstrated a higher diversity of resistance mutations than Russian isolates (p = 0.003). The rifampicin and isoniazid resistance compensatory mutations were detected in some MDR strains, but they were not widespread. The molecular mechanisms of M. tuberculosis adaptation to anti-TB treatment are not unique to the pediatric strains, but they reflect the general situation with TB in Russia and China.

Evolution and epidemic success of Mycobacterium tuberculosis in eastern China: evidence from a prospective study

BACKGROUND

Lineage distribution of Mycobacterium tuberculosis (Mtb) isolates is strongly associated with geographically distinct human populations, and its transmission can be further impacted by the bacterial genome. However, the epidemic success of Mtb isolates at an individual level was unknown in eastern China. Knowledge regarding the emergence and transmission of Mtb isolates as well as relevant factors may offer a new solution to curb the spread of the disease. Thus, this study aims to reveal the evolution and epidemic success of Mtb isolates in eastern China.

RESULTS

Of initial 1040 isolates, 997 were retained after removing duplicates and those with insufficient sequencing depth. Of the final samples, 733 (73.52%) were from Zhejiang Province, and 264 (26.48%) were from Shanghai City. Lineage 2 and lineage 4 accounted for 80.44% and 19.56%, with common ancestors dating around 7017 years ago and 6882 years ago, respectively. Sub-lineage L2.2 (80.34%) contributed the majority of total isolates, followed by L4.4 (8.93%) and L4.5 (8.43%). Additionally, 51 (5.12%) isolates were identified to be multidrug-resistant (MDR), of which 21 (29.17%) were pre-extensively drug-resistant (pre-XDR). One clade harboring katG S315T mutation may date back to 65 years ago and subsequently acquired mutations conferring resistance to another five antibiotic drugs. The prevalence of compensatory mutation was the highest in pre-XDR isolates (76.19%), followed by MDR isolates (47.06%) and other drug-resistant isolates (20.60%). Time-scaled haplotypic density analyses suggested comparable success indices between lineage 2 and lineage 4 (P = 0.306), and drug resistance did not significantly promote the transmission of Mtb isolates (P = 0.340). But for pre-XDR isolates, we found a higher success index in those with compensatory mutations (P = 0.025). Mutations under positive selection were found in genes associated with resistance to second-line injectables (whiB6) and drug tolerance (prpR) in both lineage 2 and lineage 4.

CONCLUSIONS

Our study demonstrates the population expansion of lineage 2 and lineage 4 in eastern China, with comparable transmission capacity, while accumulation of resistance mutations does not necessarily facilitate the success of Mtb isolates. Compensatory mutations usually accompany drug resistance and significantly contribute to the epidemiological transmission of pre-XDR strains. Prospective molecular surveillance is required to further monitor the emergence and spread of pre-XDR/XDR strains in eastern China.

Structural insights into RNase J that plays an essential role in Mycobacterium tuberculosis RNA metabolism

Ribonucleases (RNases) are responsible for RNA metabolism. RNase J, the core enzyme of the RNA degradosome, plays an essential role in global mRNA decay. Emerging evidence showed that the RNase J of Mycobacterium tuberculosis (Mtb-RNase J) could be an excellent target for treating Mtb infection. Here, crystal structures of Mtb-RNase J in apo-state and complex with the single-strand RNA reveal the conformational change upon RNA binding and hydrolysis. Mtb-RNase J forms an active homodimer through the interactions between the β-CASP and the β-lactamase domain. Knockout of RNase J slows the growth rate and changes the colony morphologies and cell length in Mycobacterium smegmatis, which is restored by RNase J complementation. Finally, RNA-seq analysis shows that the knockout strain significantly changes the expression levels of 49 genes in metabolic pathways. Thus, our current study explores the structural basis of Mtb-RNase J and might provide a promising candidate in pharmacological treatment for tuberculosis.

High-throughput phenogenotyping of Mycobacteria tuberculosis clinical strains reveals bacterial determinants of treatment outcomes

Background

Combatting the tuberculosis (TB) epidemic caused by Mycobacterium tuberculosis ( Mtb ) necessitates a better understanding of the factors contributing to patient clinical outcomes and transmission. While host and environmental factors have been evaluated, the impact of Mtb genetic background and phenotypic diversity is underexplored. Previous work has made associations between Mtb genetic lineages and some clinical and epidemiological features, but the bacterial traits underlying these connections are largely unknown.

Methods

We developed a high-throughput functional genomics platform for defining genotype-phenotype relationships across a panel of Mtb clinical isolates. These phenotypic fitness profiles function as intermediate traits which can be linked to Mtb genetic variants and associated with clinical and epidemiological outcomes. We applied this approach to a collection of 158 Mtb strains from a study of Mtb transmission in Ho Chi Minh City, Vietnam. Mtb strains were genetically tagged in multiplicate, which allowed us to pool the strains and assess in vitro competitive fitness using deep sequencing across a set of 14 host-relevant antibiotic and metabolic conditions. Phylogenetic and monogenic associations with these intermediate traits were identified and then associated with clinical outcomes.

Findings

Mtb clinical strains have a broad range of growth and drug response dynamics that can be clustered by their phylogenetic relationships. We identified novel monogenic associations with Mtb fitness in various metabolic and antibiotic conditions. Among these, we find that mutations in Rv1339 , a phosphodiesterase, which were identified through their association with slow growth in glycerol, are further associated with treatment failure. We also identify a previously uncharacterized subclade of Lineage 1 strains (L1.1.1.1) that is phenotypically distinguished by slow growth under most antibiotic and metabolic stress conditions in vitro . This clade is associated with cavitary disease, treatment failure, and demonstrates increased transmission potential.

Interpretation

High-throughput phenogenotyping of Mtb clinical strains enabled bacterial intermediate trait identification that can provide a mechanistic link between Mtb genetic variation and patient clinical outcomes. Mtb strains associated with cavitary disease, treatment failure, and transmission potential display intermediate phenotypes distinguished by slow growth under various antibiotic and metabolic conditions. These data suggest that Mtb growth regulation is an adaptive advantage for host bacterial success in human populations, in at least some circumstances. These data further suggest markers for the underlying bacterial processes that govern these clinical outcomes.

Funding

National Institutes of Allergy and Infectious Diseases: P01 AI132130 (SS, SMF); P01 AI143575 (XW, SMF); U19 AI142793 (QL, SMF); 5T32AI132120-03 (SS); 5T32AI132120-04 (SS); 5T32AI049928-17 (SS) Wellcome Trust Fellowship in Public Health and Tropical Medicine: 097124/Z/11/Z (NTTT) National Health and Medical Research Council (NHMRC)/A*STAR joint call: APP1056689 (SJD) The funding sources had no involvement in study methodology, data collection, analysis, and interpretation nor in the writing or submission of the manuscript.

Research in context

Evidence before this study: We used different combinations of the words mycobacterium tuberculosis, tuberculosis, clinical strains, intermediate phenotypes, genetic barcoding, phenogenomics, cavitary disease, treatment failure, and transmission to search the PubMed database for all studies published up until January 20 th , 2022. We only considered English language publications, which biases our search. Previous work linking Mtb determinants to clinical or epidemiological data has made associations between bacterial lineage, or less frequently, genetic polymorphisms to in vitro or in vivo models of pathogenesis, transmission, and clinical outcomes such as cavitary disease, treatment failure, delayed culture conversion, and severity. Many of these studies focus on the global pandemic Lineage 2 and Lineage 4 Mtb strains due in part to a deletion in a polyketide synthase implicated in host-pathogen interactions. There are a number of Mtb GWAS studies that have led to novel genetic determinants of in vitro drug resistance and tolerance. Previous Mtb GWAS analyses with clinical outcomes did not experimentally test any predicted phenotypes of the clinical strains. Published laboratory-based studies of Mtb clinical strains involve relatively small numbers of strains, do not identify the genetic basis of relevant phenotypes, or link findings to the corresponding clinical outcomes. There are two recent studies of other pathogens that describe phenogenomic analyses. One study of 331 M. abscessus clinical strains performed one-by-one phenotyping to identify bacterial features associated with clearance of infection and another details a competition experiment utilizing three barcoded Plasmodium falciparum clinical isolates to assay antimalarial fitness and resistance. Added value of this study: We developed a functional genomics platform to perform high-throughput phenotyping of Mtb clinical strains. We then used these phenotypes as intermediate traits to identify novel bacterial genetic features associated with clinical outcomes. We leveraged this platform with a sample of 158 Mtb clinical strains from a cross sectional study of Mtb transmission in Ho Chi Minh City, Vietnam. To enable high-throughput phenotyping of large numbers of Mtb clinical isolates, we applied a DNA barcoding approach that has not been previously utilized for the high-throughput analysis of Mtb clinical strains. This approach allowed us to perform pooled competitive fitness assays, tracking strain fitness using deep sequencing. We measured the replicative fitness of the clinical strains in multiplicate under 14 metabolic and antibiotic stress condition. To our knowledge, this is the largest phenotypic screen of Mtb clinical isolates to date. We performed bacterial GWAS to delineate the Mtb genetic variants associated with each fitness phenotype, identifying monogenic associations with several conditions. We then defined Mtb phenotypic and genetic features associated with clinical outcomes. We find that a subclade of Mtb strains, defined by variants largely involved in fatty acid metabolic pathways, share a universal slow growth phenotype that is associated with cavitary disease, treatment failure and increased transmission potential in Vietnam. We also find that mutations in Rv1339 , a poorly characterized phosphodiesterase, also associate with slow growth in vitro and with treatment failure in patients. Implications of all the available evidence: Phenogenomic profiling demonstrates that Mtb strains exhibit distinct growth characteristics under metabolic and antibiotic stress conditions. These fitness profiles can serve as intermediate traits for GWAS and association with clinical outcomes. Intermediate phenotyping allows us to examine potential processes by which bacterial strain differences contribute to clinical outcomes. Our study identifies clinical strains with slow growth phenotypes under in vitro models of antibiotic and host-like metabolic conditions that are associated with adverse clinical outcomes. It is possible that the bacterial intermediate phenotypes we identified are directly related to the mechanisms of these outcomes, or they may serve as markers for the causal yet unidentified bacterial determinants. Via the intermediate phenotyping, we also discovered a surprising diversity in Mtb responses to the new anti-mycobacterial drugs that target central metabolic processes, which will be important in considering roll-out of these new agents. Our study and others that have identified Mtb determinants of TB clinical and epidemiological phenotypes should inform efforts to improve diagnostics and drug regimen design.

Spatiotemporal perspectives on tuberculosis chemotherapy

Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis (TB), accounts for over ten million infections and over 1.5 million deaths every year [1]. Upon infection, the seesaw between Mtb and our immune systems creates microenvironments that are compositionally distinctive and changing over time. While the field has begun to better understand the spatial complexity of TB disease, our understanding and experimental dissection of the temporal dynamics of TB and TB drug treatment is much more rudimentary. However, it is the combined spatiotemporal heterogeneity of TB disease that creates niches and time windows within which the pathogen can survive and thrive during treatment. Here, we review the emerging data on the interactions of spatial and temporal dynamics as they relate to TB disease and treatment. A better understanding of the interactions of Mtb, host, and antibiotics through space and time will elucidate treatment failure and potentially identify opportunities for new TB treatment regimens.

The implication of Mycobacterium tuberculosis-mediated metabolism of targeted xenobiotics

Drug metabolism is generally associated with liver enzymes. However, in the case of Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis (TB), Mtb-mediated drug metabolism plays a significant role in treatment outcomes. Mtb is equipped with enzymes that catalyse biotransformation reactions on xenobiotics with consequences either in its favour or as a hindrance by deactivating or activating chemical entities, respectively. Considering the range of chemical reactions involved in the biosynthetic pathways of Mtb, information related to the biotransformation of antitubercular compounds would provide opportunities for the development of new chemical tools to study successful TB infections while also highlighting potential areas for drug discovery, host-directed therapy, dose optimization and elucidation of mechanisms of action. In this Review, we discuss Mtb-mediated biotransformations and propose a holistic approach to address drug metabolism in TB drug discovery and related areas.

Mycobacterium tuberculosis-mediated metabolism of xenobiotics poses an important research question for antitubercular drug discovery. Identification of the metabolic fate of compounds can inform requisite structure–activity relationship strategies early on in a drug discovery programme towards improving the properties of the compound.

Navigating a Path to Rifampicin Resistance in Tuberculosis

For model bacteria, genetic drug resistance usually arises from antibiotic-tolerant subpopulations, but whether this is true for the globally important pathogen Mycobacterium tuberculosis-the cause of tuberculosis-is not known. Here, we discuss a recent article by Sebastian et al. (J. Sebastian, A. Thomas, C. Levine, R. Shrestha, et al., mBio 14:e0279522, 2023, 10.1128/mbio.02795-22) which leverages a robotic transwell microtiter experimental system coupled with deep sequencing of a barcoded library of M. tuberculosis to answer this question for rifampicin resistance. The authors investigate two distinct forms of antibiotic-tolerant subpopulations-classical tolerance, characterized by prolonged minimum duration of killing, and "differentially detectable" (DD) bacilli that are viable but can be recovered only in liquid medium as opposed to plating. They demonstrate that, indeed, resistance arises preferentially from both rifampicin-tolerant subpopulations, though earlier in the DD population. Use of barcoded libraries and parallel culture systems shows promise in investigating phenotypes mediated by minority subpopulations of bacteria such as development of antibiotic resistance.

Mycobacterium tuberculosis Requires the Outer Membrane Lipid Phthiocerol Dimycocerosate for Starvation-Induced Antibiotic Tolerance

Tolerance of Mycobacterium tuberculosis to antibiotics contributes to the long duration of tuberculosis (TB) treatment and the emergence of drug-resistant strains. M. tuberculosis drug tolerance is induced by nutrient restriction, but the genetic determinants that promote antibiotic tolerance triggered by nutrient limitation have not been comprehensively identified. Here, we show that M. tuberculosis requires production of the outer membrane lipid phthiocerol dimycocerosate (PDIM) to tolerate antibiotics under nutrient-limited conditions. We developed an arrayed transposon (Tn) mutant library in M. tuberculosis Erdman and used orthogonal pooling and transposon sequencing (Tn-seq) to map the locations of individual mutants in the library. We screened a subset of the library (~1,000 mutants) by Tn-seq and identified 32 and 102 Tn mutants with altered tolerance to antibiotics under stationary-phase and phosphate-starved conditions, respectively. Two mutants recovered from the arrayed library, ppgK::Tn and clpS::Tn, showed increased susceptibility to two different drug combinations under both nutrient-limited conditions, but their phenotypes were not complemented by the Tn-disrupted gene. Whole-genome sequencing revealed single nucleotide polymorphisms in both the ppgK::Tn and clpS::Tn mutants that prevented PDIM production. Complementation of the clpS::Tn ppsD Q291* mutant with ppsD restored PDIM production and antibiotic tolerance, demonstrating that loss of PDIM sensitized M. tuberculosis to antibiotics. Our data suggest that drugs targeting production of PDIM, a critical M. tuberculosis virulence determinant, have the potential to enhance the efficacy of existing antibiotics, thereby shortening TB treatment and limiting development of drug resistance. IMPORTANCE Mycobacterium tuberculosis causes 10 million cases of active TB disease and over 1 million deaths worldwide each year. TB treatment is complex, requiring at least 6 months of therapy with a combination of antibiotics. One factor that contributes to the length of TB treatment is M. tuberculosis phenotypic antibiotic tolerance, which allows the bacteria to survive prolonged drug exposure even in the absence of genetic mutations causing drug resistance. Here, we report a genetic screen to identify M. tuberculosis genes that promote drug tolerance during nutrient starvation. Our study revealed the outer membrane lipid phthiocerol dimycocerosate (PDIM) as a key determinant of M. tuberculosis antibiotic tolerance triggered by nutrient starvation. Our study implicates PDIM synthesis as a potential target for development of new TB drugs that would sensitize M. tuberculosis to existing antibiotics to shorten TB treatment.

Single Cell Killing Kinetics Differentiate Phenotypic Bacterial Responses to Different Antibacterial Classes

With the spread of multidrug-resistant bacteria, there has been an increasing focus on molecular classes that have not yet yielded an antibiotic. A key capability for assessing and prescribing new antibacterial treatments is to compare the effects antibacterial agents have on bacterial growth at a phenotypic, single-cell level. Here, we combined time-lapse microscopy with microfluidics to investigate the concentration-dependent killing kinetics of stationary-phase Escherichia coli cells. We used antibacterial agents from three different molecular classes, β-lactams and fluoroquinolones, with the known antibiotics ampicillin and ciprofloxacin, respectively, and a new experimental class, protein Ψ-capsids. We found that bacterial cells elongated when treated with ampicillin and ciprofloxacin used at their minimum inhibitory concentration (MIC). This was in contrast to Ψ-capsids, which arrested bacterial elongation within the first two hours of treatment. At concentrations exceeding the MIC, all the antibacterial agents tested arrested bacterial growth within the first 2 h of treatment. Further, our single-cell experiments revealed differences in the modes of action of three different agents. At the MIC, ampicillin and ciprofloxacin caused the lysis of bacterial cells, whereas at higher concentrations, the mode of action shifted toward membrane disruption. The Ψ-capsids killed cells by disrupting their membranes at all concentrations tested. Finally, at increasing concentrations, ampicillin and Ψ-capsids reduced the fraction of the population that survived treatment in a viable but nonculturable state, whereas ciprofloxacin increased this fraction. This study introduces an effective capability to differentiate the killing kinetics of antibacterial agents from different molecular classes and offers a high content analysis of antibacterial mechanisms at the single-cell level. IMPORTANCE Antibiotics act against bacterial pathogens by inhibiting their growth or killing them directly. Different modes of action determine different antibacterial responses, whereas phenotypic differences in bacteria can challenge the efficacy of antibiotics. Therefore, it is important to be able to differentiate the concentration-dependent killing kinetics of antibacterial agents at a single-cell level, in particular for molecular classes which have not yielded an antibiotic before. Here, we measured single-cell responses using microfluidics-enabled imaging, revealing that a novel class of antibacterial agents, protein Ψ-capsids, arrests bacterial elongation at the onset of treatment, whereas elongation continues for cells treated with β-lactam and fluoroquinolone antibiotics. The study advances our current understanding of antibacterial function and offers an effective strategy for the comparative design of new antibacterial therapies, as well as clinical antibiotic susceptibility testing.

GWAS and functional studies suggest a role for altered DNA repair in the evolution of drug resistance in Mycobacterium tuberculosis

The emergence of drug resistance in Mycobacterium tuberculosis (Mtb) is alarming and demands in-depth knowledge for timely diagnosis. We performed genome-wide association analysis using 2237 clinical strains of Mtb to identify novel genetic factors that evoke drug resistance. In addition to the known direct targets, we identified for the first time, a strong association between mutations in DNA repair genes and the multidrug-resistant phenotype. To evaluate the impact of variants identified in the clinical samples in the evolution of drug resistance, we utilized knockouts and complemented strains in Mycobacterium smegmatis and Mtb. Results show that variant mutations compromised the functions of MutY and UvrB. MutY variant showed enhanced survival compared with wild-type (Rv) when the Mtb strains were subjected to multiple rounds of ex vivo antibiotic stress. In an in vivo guinea pig infection model, the MutY variant outcompeted the wild-type strain. We show that novel variant mutations in the DNA repair genes collectively compromise their functions and contribute to better survival under antibiotic/host stress conditions.

A research program-linked, course-based undergraduate research experience that allows undergraduates to participate in current research on mycobacterial gene regulation

Undergraduate instructional biology laboratories are typically taught within two paradigms. Some labs focus on protocols and techniques delivered in "cookbook" format with defined experimental outcomes. There is increasing momentum to alternatively employ student-driven, open-ended, and discovery-based strategies, often via course-based undergraduate research experiences (CUREs) using crowd-sourcing initiatives. A fraction of students also participate in funded research in faculty research labs, where they have opportunities to work on projects designed to expand the frontiers of human knowledge. These experiences are widely recognized as valuable but are not scalable, as most institutions have many more undergraduates than research lab positions. We sought to address this gap through our department's curriculum by creating an opportunity for students to participate in the real-world research process within a laboratory course. We conceived, developed, and delivered an authentic, guided research experience to students in an upper-level molecular biology laboratory course. We refer to this model as a "research program-linked CURE." The research questions come directly from a faculty member's research lab and evolve along with that research program. Students study post-transcriptional regulation in mycobacteria. We use current molecular biology methodologies to test hypotheses like "UTRs affect RNA and protein expression levels," "there is functional redundancy among RNA helicases," and "carbon starvation alters mRNA 5' end chemistries." We conducted standard assessments and developed a customized "Skills and Concepts Inventory" survey to gauge how well the course met our student learning outcomes. We report the results of our assessments and describe challenges addressed during development and execution of the course, including organizing activities to fit within an instructional lab, balancing breadth with depth, and maintaining authenticity while giving students the experience of obtaining interpretable and novel results. Our data suggest student learning was enhanced through this truly authentic research approach. Further, students were able to perceive they were participants and contributors within an active research paradigm. Students reported increases in their self-identification as scientists, and a positive impact on their career trajectories. An additional benefit was reciprocation back to the funded research laboratory, by funneling course alumni, results, materials, and protocols.

Tools to develop antibiotic combinations that target drug tolerance in Mycobacterium tuberculosis

Combination therapy is necessary to treat tuberculosis to decrease the rate of disease relapse and prevent the acquisition of drug resistance, and shorter regimens are urgently needed. The adaptation of Mycobacterium tuberculosis to various lesion microenvironments in infection induces various states of slow replication and non-replication and subsequent antibiotic tolerance. This non-heritable tolerance to treatment necessitates lengthy combination therapy. Therefore, it is critical to develop combination therapies that specifically target the different types of drug-tolerant cells in infection. As new tools to study drug combinations earlier in the drug development pipeline are being actively developed, we must consider how to best model the drug-tolerant cells to use these tools to design the best antibiotic combinations that target those cells and shorten tuberculosis therapy. In this review, we discuss the factors underlying types of drug tolerance, how combination therapy targets these populations of bacteria, and how drug tolerance is currently modeled for the development of tuberculosis multidrug therapy. We highlight areas for future studies to develop new tools that better model drug tolerance in tuberculosis infection specifically for combination therapy testing to bring the best drug regimens forward to the clinic.

Tuberculosis treatment failure associated with evolution of antibiotic resilience

The widespread use of antibiotics has placed bacterial pathogens under intense pressure to evolve new survival mechanisms. Genomic analysis of 51,229 Mycobacterium tuberculosis (Mtb)clinical isolates has identified an essential transcriptional regulator, Rv1830, herein called resR for resilience regulator, as a frequent target of positive (adaptive) selection. resR mutants do not show canonical drug resistance or drug tolerance but instead shorten the post-antibiotic effect, meaning that they enable Mtb to resume growth after drug exposure substantially faster than wild-type strains. We refer to this phenotype as antibiotic resilience. ResR acts in a regulatory cascade with other transcription factors controlling cell growth and division, which are also under positive selection in clinical isolates of Mtb. Mutations of these genes are associated with treatment failure and the acquisition of canonical drug resistance.

Immune cell interactions in tuberculosis

Despite having been identified as the organism that causes tuberculosis in 1882, Mycobacterium tuberculosis has managed to still evade our understanding of the protective immune response against it, defying the development of an effective vaccine. Technology and novel experimental models have revealed much new knowledge, particularly with respect to the heterogeneity of the bacillus and the host response. This review focuses on certain immunological elements that have recently yielded exciting data and highlights the importance of taking a holistic approach to understanding the interaction of M. tuberculosis with the many host cells that contribute to the development of protective immunity.

Anti-tuberculosis treatment strategies and drug development: challenges and priorities

Despite two decades of intensified research to understand and cure tuberculosis disease, biological uncertainties remain and hamper progress. However, owing to collaborative initiatives including academia, the pharmaceutical industry and non-for-profit organizations, the drug candidate pipeline is promising. This exceptional success comes with the inherent challenge of prioritizing multidrug regimens for clinical trials and revamping trial designs to accelerate regimen development and capitalize on drug discovery breakthroughs. Most wanted are markers of progression from latent infection to active pulmonary disease, markers of drug response and predictors of relapse, in vitro tools to uncover synergies that translate clinically and animal models to reliably assess the treatment shortening potential of new regimens. In this Review, we highlight the benefits and challenges of 'one-size-fits-all' regimens and treatment duration versus individualized therapy based on disease severity and host and pathogen characteristics, considering scientific and operational perspectives.

Unraveling the mechanisms of intrinsic drug resistance in Mycobacterium tuberculosis

Tuberculosis (TB) is among the most difficult infections to treat, requiring several months of multidrug therapy to produce a durable cure. The reasons necessitating long treatment times are complex and multifactorial. However, one major difficulty of treating TB is the resistance of the infecting bacterium, Mycobacterium tuberculosis (Mtb), to many distinct classes of antimicrobials. This review will focus on the major gaps in our understanding of intrinsic drug resistance in Mtb and how functional and chemical-genetics can help close those gaps. A better understanding of intrinsic drug resistance will help lay the foundation for strategies to disarm and circumvent these mechanisms to develop more potent antitubercular therapies.

The evolving biology of Mycobacterium tuberculosis drug resistance

Tuberculosis, caused by Mycobacterium tuberculosis (Mtb) is an ancient disease that has remained a leading cause of infectious death. Mtb has evolved drug resistance to every antibiotic regimen ever introduced, greatly complicating treatment, lowering rates of cure and menacing TB control in parts of the world. As technology has advanced, our understanding of antimicrobial resistance has improved, and our models of the phenomenon have evolved. In this review, we focus on recent research progress that supports an updated model for the evolution of drug resistance in Mtb. We highlight the contribution of drug tolerance on the path to resistance, and the influence of heterogeneity on tolerance. Resistance is likely to remain an issue for as long as drugs are needed to treat TB. However, with technology driving new insights and careful management of newly developed resources, antimicrobial resistance need not continue to threaten global progress against TB, as it has done for decades.

Mycobacterium tuberculosis functional genetic diversity, altered drug sensitivity, and precision medicine

In the face of the unrelenting global burden of tuberculosis (TB), antibiotics remain our most effective tools to save lives and control the spread of Mycobacterium tuberculosis (Mtb). However, we confront a dual challenge in our use of antibiotics: simplifying and shortening the TB drug regimen while also limiting the emergence and propagation of antibiotic resistance. This task is now more feasible due to the increasing availability of bacterial genomic data at or near the point of care. These resources create an opportunity to envision how integration of bacterial genetic determinants of antibiotic response into treatment algorithms might transform TB care. Historically, Mtb drug resistance studies focused on mutations in genes encoding antibiotic targets and the resulting increases in the minimal inhibitory concentrations (MICs) above a breakpoint value. But recent progress in elucidating the effects of functional genetic diversity in Mtb has revealed various genetic loci that are associated with drug phenotypes such as low-level MIC increases and tolerance which predict the development of resistance and treatment failure. As a result, we are now poised to advance precision medicine approaches in TB treatment. By incorporating information regarding Mtb genetic characteristics into the development of drug regimens, clinical care which tailors antibiotic treatment to maximize the likelihood of success has come into reach.

Failing upwards: Genetics-based strategies to improve antibiotic discovery and efficacy in Mycobacterium tuberculosis

Therapeutic advances in the 20th century significantly reduced tuberculosis (TB) mortality. Nonetheless, TB still poses a massive global health challenge with significant annual morbidity and mortality that has been amplified during the COVID-19 pandemic. Unlike most common bacterial infectious diseases, successful TB treatment requires months-long regimens, which complicates the ability to treat all cases quickly and effectively. Improving TB chemotherapy by reducing treatment duration and optimizing combinations of drugs is an important step to reducing relapse. In this review, we outline the limitations of current multidrug regimens against TB and have reviewed the genetic tools available to improve the identification of drug targets. The rational design of regimens that sterilize diverse phenotypic subpopulations will maximize bacterial killing while minimizing both treatment duration and infection relapse. Importantly, the TB field currently has all the necessary genetic and analytical tools to screen for and prioritize drug targets in vitro based on the vulnerability of essential and non-essential genes in the Mtb genome and to translate these findings in in vivo models. Combining genetic methods with chemical screens offers a formidable strategy to redefine the preclinical design of TB therapy by identifying powerful new targets altogether, as well as targets that lend new efficacy to existing drugs.

Transcriptional regulation and drug resistance in Mycobacterium tuberculosis

Bacterial drug resistance is one of the major challenges to present and future human health, as the continuous selection of multidrug resistant bacteria poses at serious risk the possibility to treat infectious diseases in the near future. One of the infection at higher risk to become incurable is tuberculosis, due to the few drugs available in the market against Mycobacterium tuberculosis. Drug resistance in this species is usually due to point mutations in the drug target or in proteins required to activate prodrugs. However, another interesting and underexplored aspect of bacterial physiology with important impact on drug susceptibility is represented by the changes in transcriptional regulation following drug exposure. The main regulators involved in this phenomenon in M. tuberculosis are the sigma factors, and regulators belonging to the WhiB, GntR, XRE, Mar and TetR families. Better understanding the impact of these regulators in survival to drug treatment might contribute to identify new drug targets and/or to design new strategies of intervention.

Transcontinental spread and evolution of Mycobacterium tuberculosis W148 European/Russian clade toward extensively drug resistant tuberculosis

Transmission-driven multi-/extensively drug resistant (M/XDR) tuberculosis (TB) is the largest single contributor to human mortality due to antimicrobial resistance. A few major clades of the Mycobacterium tuberculosis complex belonging to lineage 2, responsible for high prevalence of MDR-TB in Eurasia, show outstanding transnational distributions. Here, we determined factors underlying the emergence and epidemic spread of the W148 clade by genome sequencing and Bayesian demogenetic analyses of 720 isolates from 23 countries. We dated a common ancestor around 1963 and identified two successive epidemic expansions in the late 1980s and late 1990s, coinciding with major socio-economic changes in the post-Soviet Era. These population expansions favored accumulation of resistance mutations to up to 11 anti-TB drugs, with MDR evolving toward additional resistances to fluoroquinolones and second-line injectable drugs within 20 years on average. Timescaled haplotypic density analysis revealed that widespread acquisition of compensatory mutations was associated with transmission success of XDR strains. Virtually all W148 strains harbored a hypervirulence-associated ppe38 gene locus, and incipient recurrent emergence of prpR mutation-mediated drug tolerance was detected. The outstanding genetic arsenal of this geographically widespread M/XDR strain clade represents a “perfect storm” that jeopardizes the successful introduction of new anti-M/XDR-TB antibiotic regimens.

An outbreak of the multidrug-resistant Mycobacterium tuberculosis lineage W148 has spread widely across Russia, Central Asia and Europe. Here, the authors use whole genome sequences of ~700 isolates of this lineage collected over ~20 years to analyze its spread, evolution of drug resistance, and impact of compensatory mutations.

High fluoroquinolone resistance proportions among multidrug-resistant tuberculosis driven by dominant L2 Mycobacterium tuberculosis clones in the Mumbai Metropolitan Region

Background

Multidrug-resistant (MDR) Mycobacterium tuberculosis complex (MTBC) strains are a serious health problem in India, also contributing to one-fourth of the global MDR tuberculosis (TB) burden. About 36% of the MDR MTBC strains are reported fluoroquinolone (FQ) resistant leading to high pre-extensively drug-resistant (pre-XDR) and XDR-TB (further resistance against bedaquiline and/or linezolid) rates. Still, factors driving the MDR/pre-XDR epidemic in India are not well defined.

Methods

In a retrospective study, we analyzed 1852 consecutive MTBC strains obtained from patients from a tertiary care hospital laboratory in Mumbai by whole genome sequencing (WGS). Univariate and multivariate statistics was used to investigate factors associated with pre-XDR. Core genome multi locus sequence typing, time scaled haplotypic density (THD) method and homoplasy analysis were used to analyze epidemiological success, and positive selection in different strain groups, respectively.

Results

In total, 1016 MTBC strains were MDR, out of which 703 (69.2%) were pre-XDR and 45 (4.4%) were XDR. Cluster rates were high among MDR (57.8%) and pre-XDR/XDR (79%) strains with three dominant L2 (Beijing) strain clusters (Cl 1–3) representing half of the pre-XDR and 40% of the XDR-TB cases. L2 strains were associated with pre-XDR/XDR-TB (P < 0.001) and, particularly Cl 1–3 strains, had high first-line and FQ resistance rates (81.6–90.6%). Epidemic success analysis using THD showed that L2 strains outperformed L1, L3, and L4 strains in short- and long-term time scales. More importantly, L2 MDR and MDR + strains had higher THD success indices than their not-MDR counterparts. Overall, compensatory mutation rates were highest in L2 strains and positive selection was detected in genes of L2 strains associated with drug tolerance (prpB and ppsA) and virulence (Rv2828c). Compensatory mutations in L2 strains were associated with a threefold increase of THD indices, suggesting improved transmissibility.

Conclusions

Our data indicate a drastic increase of FQ resistance, as well as emerging bedaquiline resistance which endangers the success of newly endorsed MDR-TB treatment regimens. Rapid changes in treatment and control strategies are required to contain transmission of highly successful pre-XDR L2 strains in the Mumbai Metropolitan region but presumably also India-wide.

Understanding the contribution of metabolism to Mycobacterium tuberculosis drug tolerance

Treatment of Mycobacterium tuberculosis (Mtb) infections is particularly arduous. One challenge to effectively treating tuberculosis is that drug efficacy in vivo often fails to match drug efficacy in vitro. This is due to multiple reasons, including inadequate drug concentrations reaching Mtb at the site of infection and physiological changes of Mtb in response to host derived stresses that render the bacteria more tolerant to antibiotics. To more effectively and efficiently treat tuberculosis, it is necessary to better understand the physiologic state of Mtb that promotes drug tolerance in the host. Towards this end, multiple studies have converged on bacterial central carbon metabolism as a critical contributor to Mtb drug tolerance. In this review, we present the evidence that changes in central carbon metabolism can promote drug tolerance, depending on the environment surrounding Mtb. We posit that these metabolic pathways could be potential drug targets to stymie the development of drug tolerance and enhance the efficacy of current antimicrobial therapy.

Central carbon metabolism remodeling as a mechanism to develop drug tolerance and drug resistance in Mycobacterium tuberculosis

Suboptimal efficacy of the current antibiotic regimens and frequent emergence of antibiotic-resistant Mycobacterium tuberculosis (Mtb), an etiological agent of tuberculosis (TB), render TB the world’s deadliest infectious disease before the COVID-19 outbreak. Our outdated TB treatment method is designed to eradicate actively replicating populations of Mtb. Unfortunately, accumulating evidence suggests that a small population of Mtb can survive antimycobacterial pressure of antibiotics by entering a “persister” state (slowly replicating or non-replicating and lacking a stably heritable antibiotic resistance, termed drug tolerance). The formation of drug-tolerant Mtb persisters is associated with TB treatment failure and is thought to be an adaptive strategy for eventual development of permanent genetic mutation-mediated drug resistance. Thus, the molecular mechanisms behind persister formation and drug tolerance acquisition are a source of new antibiotic targets to eradicate both Mtb persisters and drug-resistant Mtb. As Mtb persisters are genetically identical to antibiotic susceptible populations, metabolomics has emerged as a vital biochemical tool to differentiate these populations by determining phenotypic shifts and metabolic reprogramming. Metabolomics, which provides detailed insights into the molecular basis of drug tolerance and resistance in Mtb, has unique advantages over other techniques by its ability to identify specific metabolic differences between the two genetically identical populations. This review summarizes the recent advances in our understanding of the metabolic adaptations used by Mtb persisters to achieve intrinsic drug tolerance and facilitate the emergence of drug resistance. These findings present metabolomics as a powerful tool to identify previously unexplored antibiotic targets and improved combinations of drug regimens against drug-resistant TB infection.

Are all antibiotic persisters created equal?

Antibiotic persisters are a sub-population of bacteria able to survive in the presence of bactericidal antibiotic despite the lack of heritable drug resistance mechanisms. This phenomenon exists across many bacterial species and is observed for many different antibiotics. Though these bacteria are often described as “multidrug persisters” very few experiments have been carried out to determine the homogeneity of a persister population to different drugs. Further, there is much debate in the field as to the origins of a persister cell. Is it formed spontaneously? Does it form in response to stress? These questions are particularly pressing in the field of Mycobacterium tuberculosis, where persisters may play a crucial role in the required length of treatment and the development of multidrug resistant organisms. Here we aim to interpret the known mechanisms of antibiotic persistence and how they may relate to improving treatments for M. tuberculosis, exposing the gaps in knowledge that prevent us from answering the question: Are all antibiotic persisters created equal?

A convolutional neural network highlights mutations relevant to antimicrobial resistance in Mycobacterium tuberculosis

Long diagnostic wait times hinder international efforts to address antibiotic resistance in M. tuberculosis. Pathogen whole genome sequencing, coupled with statistical and machine learning models, offers a promising solution. However, generalizability and clinical adoption have been limited by a lack of interpretability, especially in deep learning methods. Here, we present two deep convolutional neural networks that predict antibiotic resistance phenotypes of M. tuberculosis isolates: a multi-drug CNN (MD-CNN), that predicts resistance to 13 antibiotics based on 18 genomic loci, with AUCs 82.6-99.5% and higher sensitivity than state-of-the-art methods; and a set of 13 single-drug CNNs (SD-CNN) with AUCs 80.1-97.1% and higher specificity than the previous state-of-the-art. Using saliency methods to evaluate the contribution of input sequence features to the SD-CNN predictions, we identify 18 sites in the genome not previously associated with resistance. The CNN models permit functional variant discovery, biologically meaningful interpretation, and clinical applicability.

Loss of RNase J leads to multi-drug tolerance and accumulation of highly structured mRNA fragments in Mycobacterium tuberculosis

Despite the existence of well-characterized, canonical mutations that confer high-level drug resistance to Mycobacterium tuberculosis (Mtb), there is evidence that drug resistance mechanisms are more complex than simple acquisition of such mutations. Recent studies have shown that Mtb can acquire non-canonical resistance-associated mutations that confer survival advantages in the presence of certain drugs, likely acting as stepping-stones for acquisition of high-level resistance. Rv2752c/rnj, encoding RNase J, is disproportionately mutated in drug-resistant clinical Mtb isolates. Here we show that deletion of rnj confers increased tolerance to lethal concentrations of several drugs. RNAseq revealed that RNase J affects expression of a subset of genes enriched for PE/PPE genes and stable RNAs and is key for proper 23S rRNA maturation. Gene expression differences implicated two sRNAs and ppe50-ppe51 as important contributors to the drug tolerance phenotype. In addition, we found that in the absence of RNase J, many short RNA fragments accumulate because they are degraded at slower rates. We show that the accumulated transcript fragments are targets of RNase J and are characterized by strong secondary structure and high G+C content, indicating that RNase J has a rate-limiting role in degradation of highly structured RNAs. Taken together, our results demonstrate that RNase J indirectly affects drug tolerance, as well as reveal the endogenous roles of RNase J in mycobacterial RNA metabolism.