Mapping and manipulating the Mycobacterium tuberculosis transcriptome using a transcription factor overexpression-derived regulatory network

Computational drug repositioning identifies niclosamide and tribromsalan as inhibitors of Mycobacterium tuberculosis and Mycobacterium abscessus

Tuberculosis (TB) is still a major global health challenge, killing over 1.5 million people each year, and hence, there is a need to identify and develop novel treatments for Mycobacterium tuberculosis (M. tuberculosis). The prevalence of infections caused by nontuberculous mycobacteria (NTM) is also increasing and has overtaken TB cases in the United States and much of the developed world. Mycobacterium abscessus (M. abscessus) is one of the most frequently encountered NTM and is difficult to treat. We describe the use of drug-disease association using a semantic knowledge graph approach combined with machine learning models that has enabled the identification of several molecules for testing anti-mycobacterial activity. We established that niclosamide (M. tuberculosis IC90 2.95 μM; M. abscessus IC90 59.1 μM) and tribromsalan (M. tuberculosis IC90 76.92 μM; M. abscessus IC90 147.4 μM) inhibit M. tuberculosis and M. abscessus in vitro. To investigate the mode of action, we determined the transcriptional response of M. tuberculosis and M. abscessus to both compounds in axenic log phase, demonstrating a broad effect on gene expression that differed from known M. tuberculosis inhibitors. Both compounds elicited transcriptional responses indicative of respiratory pathway stress and the dysregulation of fatty acid metabolism.

Genetically encoded transcriptional plasticity underlies stress adaptation in Mycobacterium tuberculosis

Transcriptional regulation is a critical adaptive mechanism that allows bacteria to respond to changing environments, yet the concept of transcriptional plasticity (TP) - the variability of gene expression in response to environmental changes - remains largely unexplored. In this study, we investigate the genome-wide TP profiles of Mycobacterium tuberculosis (Mtb) genes by analyzing 894 RNA sequencing samples derived from 73 different environmental conditions. Our data reveal that Mtb genes exhibit significant TP variation that correlates with gene function and gene essentiality. We also find that critical genetic features, such as gene length, GC content, and operon size independently impose constraints on TP, beyond trans-regulation. By extending our analysis to include two other Mycobacterium species -- M. smegmatis and M. abscessus -- we demonstrate a striking conservation of the TP landscape. This study provides a comprehensive understanding of the TP exhibited by mycobacteria genes, shedding light on this significant, yet understudied, genetic feature encoded in bacterial genomes.

Adaptation of the Mycobacterium tuberculosis transcriptome to biofilm growth

Mycobacterium tuberculosis (M. tb), the causative agent of tuberculosis (TB), is a leading global cause of death from infectious disease. Biofilms are increasingly recognized as a relevant growth form during M. tb infection and may impede treatment by enabling bacterial drug and immune tolerance. M. tb has a complicated regulatory network that has been well-characterized for many relevant disease states, including dormancy and hypoxia. However, despite its importance, our knowledge of the genes and pathways involved in biofilm formation is limited. Here we characterize the biofilm transcriptomes of fully virulent clinical isolates and find that the regulatory systems underlying biofilm growth vary widely between strains and are also distinct from regulatory programs associated with other environmental cues. We used experimental evolution to investigate changes to the transcriptome during adaptation to biofilm growth and found that the application of a uniform selection pressure resulted in loss of strain-to-strain variation in gene expression, resulting in a more uniform biofilm transcriptome. The adaptive trajectories of transcriptomes were shaped by the genetic background of the M. tb population leading to convergence on a sub-lineage specific transcriptome. We identified widespread upregulation of non-coding RNA (ncRNA) as a common feature of the biofilm transcriptome and hypothesize that ncRNA function in genome-wide modulation of gene expression, thereby facilitating rapid regulatory responses to new environments. These results reveal a new facet of the M. tb regulatory system and provide valuable insight into how M. tb adapts to new environments.

Evolution and emergence of Mycobacterium tuberculosis

Tuberculosis (TB) remains one of the deadliest infectious diseases in human history, prevailing even in the 21st century. The causative agents of TB are represented by a group of closely related bacteria belonging to the Mycobacterium tuberculosis complex (MTBC), which can be subdivided into several lineages of human- and animal-adapted strains, thought to have shared a last common ancestor emerged by clonal expansion from a pool of recombinogenic Mycobacterium canettii-like tubercle bacilli. A better understanding of how MTBC populations evolved from less virulent mycobacteria may allow for discovering improved TB control strategies and future epidemiologic trends. In this review, we highlight new insights into the evolution of mycobacteria at the genus level, describing different milestones in the evolution of mycobacteria, with a focus on the genomic events that have likely enabled the emergence and the dominance of the MTBC. We also review the recent literature describing the various MTBC lineages and highlight their particularities and differences with a focus on host preferences and geographic distribution. Finally, we discuss on putative mechanisms driving the evolution of tubercle bacilli and mycobacteria in general, by taking the mycobacteria-specific distributive conjugal transfer as an example.

Regulatory Intersection of Two-component System and Ser/Thr Protein Kinase Signaling in Mycobacterium tuberculosis

Phosphosignaling in bacteria is mediated by two distinct systems, the two-component systems (TCSs) and the protein Ser/Thr/Tyr, or O-phosphorylation systems. These two arms of phosphosignaling are currently thought to be largely independent from one another. We mined a deep Mycobacterium tuberculosis (Mtb) phosphoproteome and identified over 170 O-phosphorylation sites on histidine kinases and response regulators of TCSs, suggesting that the two signaling pathways extensively intersect. Several TCSs were phosphorylated on multiple sites, and many by multiple Ser/Thr protein kinases, suggesting convergent and cooperative regulatory interactions. To test in which way these O-phosphorylation sites affect TCS activity, we reconstituted the NarSL phosphorelay in vitro. The Ser/Thr protein kinase PknL phosphorylated the histidine kinase NarS and activated its autophosphorylating activity. A phosphoablative mutation at the PknL phosphorylation site Thr380 resulted in low autophosphorylating activity, whereas a phosphomimetic mutation strongly activated autophosphorylation. The phosphomimetic mutation also resulted in more efficient phosphotransfer from NarS to the response regulator NarL and suppression of gene expression. These data show control of NarSL signaling by STPKs through a phosphoswitch and point to extensive, functional crosstalk between TCSs and O-phosphorylation.

Discovering common pathogenetic processes between COVID-19 and tuberculosis by bioinformatics and system biology approach

Introduction

The coronavirus disease 2019 (COVID-19) pandemic, stemming from the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has persistently threatened the global health system. Meanwhile, tuberculosis (TB) caused by Mycobacterium tuberculosis (M. tuberculosis) still continues to be endemic in various regions of the world. There is a certain degree of similarity between the clinical features of COVID-19 and TB, but the underlying common pathogenetic processes between COVID-19 and TB are not well understood.

Methods

To elucidate the common pathogenetic processes between COVID-19 and TB, we implemented bioinformatics and systematic research to obtain shared pathways and molecular biomarkers. Here, the RNA-seq datasets (GSE196822 and GSE126614) are used to extract shared differentially expressed genes (DEGs) of COVID-19 and TB. The common DEGs were used to identify common pathways, hub genes, transcriptional regulatory networks, and potential drugs.

Results

A total of 96 common DEGs were selected for subsequent analyses. Functional enrichment analyses showed that viral genome replication and immune-related pathways collectively contributed to the development and progression of TB and COVID-19. Based on the protein-protein interaction (PPI) network analysis, we identified 10 hub genes, including IFI44L, ISG15, MX1, IFI44, OASL, RSAD2, GBP1, OAS1, IFI6, and HERC5. Subsequently, the transcription factor (TF)-gene interaction and microRNA (miRNA)-gene coregulatory network identified 61 TFs and 29 miRNAs. Notably, we identified 10 potential drugs to treat TB and COVID-19, namely suloctidil, prenylamine, acetohexamide, terfenadine, prochlorperazine, 3'-azido-3'-deoxythymidine, chlorophyllin, etoposide, clioquinol, and propofol.

Conclusion

This research provides novel strategies and valuable references for the treatment of tuberculosis and COVID-19.

Genetically encoded transcriptional plasticity underlies stress adaptation in Mycobacterium tuberculosis

Transcriptional regulation is a critical adaptive mechanism that allows bacteria to respond to changing environments, yet the concept of transcriptional plasticity (TP) remains largely unexplored. In this study, we investigate the genome-wide TP profiles of Mycobacterium tuberculosis (Mtb) genes by analyzing 894 RNA sequencing samples derived from 73 different environmental conditions. Our data reveal that Mtb genes exhibit significant TP variation that correlates with gene function and gene essentiality. We also found that critical genetic features, such as gene length, GC content, and operon size independently impose constraints on TP, beyond trans-regulation. By extending our analysis to include two other Mycobacterium species -- M. smegmatis and M. abscessus -- we demonstrate a striking conservation of the TP landscape. This study provides a comprehensive understanding of the TP exhibited by mycobacteria genes, shedding light on this significant, yet understudied, genetic feature encoded in bacterial genomes.

Aldehyde accumulation in Mycobacterium tuberculosis with defective proteasomal degradation results in copper sensitivity

Mycobacterium tuberculosis is a major human pathogen and the causative agent of tuberculosis disease. M. tuberculosis is able to persist in the face of host-derived antimicrobial molecules nitric oxide (NO) and copper (Cu). However, M. tuberculosis with defective proteasome activity is highly sensitive to NO and Cu, making the proteasome an attractive target for drug development. Previous work linked NO susceptibility with the accumulation of para-hydroxybenzaldehyde (pHBA) in M. tuberculosis mutants with defective proteasomal degradation. In this study, we found that pHBA accumulation was also responsible for Cu sensitivity in these strains. We showed that exogenous addition of pHBA to wild-type M. tuberculosis cultures sensitized bacteria to Cu to a degree similar to that of a proteasomal degradation mutant. We determined that pHBA reduced the production and function of critical Cu resistance proteins of the regulated in copper repressor (RicR) regulon. Furthermore, we extended these Cu-sensitizing effects to an aldehyde that M. tuberculosis may face within the macrophage. Collectively, this study is the first to mechanistically propose how aldehydes can render M. tuberculosis susceptible to an existing host defense and could support a broader role for aldehydes in controlling M. tuberculosis infections. IMPORTANCE M. tuberculosis is a leading cause of death by a single infectious agent, causing 1.5 million deaths annually. An effective vaccine for M. tuberculosis infections is currently lacking, and prior infection does not typically provide robust immunity to subsequent infections. Nonetheless, immunocompetent humans can control M. tuberculosis infections for decades. For these reasons, a clear understanding of how mammalian immunity inhibits mycobacterial growth is warranted. In this study, we show aldehydes can increase M. tuberculosis susceptibility to copper, an established antibacterial metal used by immune cells to control M. tuberculosis and other microbes. Given that activated macrophages produce increased amounts of aldehydes during infection, we propose host-derived aldehydes may help control bacterial infections, making aldehydes a previously unappreciated antimicrobial defense.

Rv1453 is associated with clofazimine resistance in Mycobacterium tuberculosis

Clofazimine (CFZ) has been repurposed for treating tuberculosis (TB), especially multidrug-resistant tuberculosis (MDR-TB). However, the mechanisms of resistance to clofazimine are poorly understood. We previously reported a mutation located in the intergenic region of Rv1453 that was linked to resistance to CFZ and demonstrated that an Rv1453 knockout resulted in an increased minimum inhibitory concentration (MIC) of CFZ. The current study aims to go back and describe in detail how the mutation was identified and further explore its association with CFZ resistance by testing additional 30 isolates. We investigated MICs of clofazimine against 100 clinical strains isolated from MDR-TB patients by microplate alamarBlue assay. Whole-genome sequencing (WGS) was performed on 11 clofazimine-resistant and 7 clofazimine-susceptible strains, including H37Rv. Among the 11 clofazimine-resistant mutants subjected to WGS, the rate of mutation in the intergenic region of the Rv1453 gene was 55% (6/11) in clofazimine-resistant strains. Among another 30 clofazimine-resistant clinical isolates, 27 had mutations in the intergenic region of the Rv1453 gene. A mutation in the Rv1453 gene associated with clofazimine resistance was identified, which shed light on the mechanisms of action and resistance of clofazimine. IMPORTANCE Tuberculosis (TB) is an infectious disease caused by the bacterium Mycobacterium tuberculosis, especially the emergence of multidrug-resistant tuberculosis (MDR-TB) brings great distress to humans. Clofazimine (CFZ) plays an important role in the treatment of MDR-TB. To understand the underlying mechanism of clofazimine resistance better, in this study, we review and detail the findings of the mutation of intergenic region of Rv1453 and find additional evidence that this mutation is related to clofazimine resistance in 30 additional isolates. The significance of our research is to contribute to a comprehensive understanding of clofazimine-resistant mechanisms, which is critical for reducing the emergence of resistance and for anti-TB drug discovery.

The Rip1 intramembrane protease contributes to iron and zinc homeostasis in Mycobacterium tuberculosis

Mycobacterium tuberculosis is exposed to a variety of stresses during a chronic infection, as the immune system simultaneously produces bactericidal compounds and starves the pathogen of essential nutrients. The intramembrane protease, Rip1, plays an important role in the adaptation to these stresses, at least partially by the cleavage of membrane-bound transcriptional regulators. Although Rip1 is known to be critical for surviving copper intoxication and nitric oxide exposure, these stresses do not fully account for the regulatory protein's essentiality during infection. In this work, we demonstrate that Rip1 is also necessary for growth in low-iron and low-zinc conditions, similar to those imposed by the immune system. Using a newly generated library of sigma factor mutants, we show that the known regulatory target of Rip1, SigL, shares this defect. Transcriptional profiling under iron-limiting conditions supported the coordinated activity of Rip1 and SigL and demonstrated that the loss of these proteins produces an exaggerated iron starvation response. These observations demonstrate that Rip1 coordinates several aspects of metal homeostasis and suggest that a Rip1- and SigL-dependent pathway is necessary to thrive in the iron-deficient environments encountered during infection. IMPORTANCE Metal homeostasis represents a critical point of interaction between the mammalian immune system and potential pathogens. While the host attempts to intoxicate microbes with high concentrations of copper or starve the invader of iron and zinc, successful pathogens have acquired mechanisms to overcome these defenses. Our work identifies a regulatory pathway consisting of the Rip1 intramembrane protease and the sigma factor, SigL, that is essential for the important human pathogen, Mycobacterium tuberculosis, to grow in low-iron or low-zinc conditions such as those encountered during infection. In conjunction with Rip1's known role in resisting copper toxicity, our work implicates this protein as a critical integration point that coordinates the multiple metal homeostatic systems required for this pathogen to survive in host tissue.

Genetically encoded transcriptional plasticity underlies stress adaptation in Mycobacterium tuberculosis

Transcriptional regulation is a critical adaptive mechanism that allows bacteria to respond to changing environments, yet the concept of transcriptional plasticity (TP) remains largely unexplored. In this study, we investigate the genome-wide TP profiles of Mycobacterium tuberculosis (Mtb) genes by analyzing 894 RNA sequencing samples derived from 73 different environmental conditions. Our data reveal that Mtb genes exhibit significant TP variation that correlates with gene function and gene essentiality. We also found that critical genetic features, such as gene length, GC content, and operon size independently impose constraints on TP, beyond trans-regulation. By extending our analysis to include two other Mycobacterium species -- M. smegmatis and M. abscessus -- we demonstrate a striking conservation of the TP landscape. This study provides a comprehensive understanding of the TP exhibited by mycobacteria genes, shedding light on this significant, yet understudied, genetic feature encoded in bacterial genomes.

ResR/McdR-regulated protein translation machinery contributes to drug resilience in Mycobacterium tuberculosis

Survival response of the human tuberculosis pathogen, Mycobacterium tuberculosis (Mtb) to a diverse environmental cues is governed through its versatile transcription regulatory mechanisms with the help of a large pool of transcription regulators (TRs). Rv1830 is one such conserved TR, which remains uncharacterized in Mtb. It was named as McdR based on an effect on cell division upon its overexpression in Mycobacterium smegmatis. Recently, it has been implicated in antibiotic resilience in Mtb and reannotated as ResR. While Rv1830 affects cell division by modulating the expression of M. smegmatis whiB2, the underlying cause of its essentiality and regulation of drug resilience in Mtb is yet to be deciphered. Here we show that ResR/McdR, encoded by ERDMAN_2020 in virulent Mtb Erdman, is pivotal for bacterial proliferation and crucial metabolic activities. Importantly, ResR/McdR directly regulates ribosomal gene expression and protein synthesis, requiring distinct disordered N-terminal sequence. Compared to control, bacteria depleted with resR/mcdR exhibit delayed recovery post-antibiotic treatment. A similar effect upon knockdown of rplN operon genes further implicates ResR/McdR-regulated protein translation machinery in attributing drug resilience in Mtb. Overall, findings from this study suggest that chemical inhibitors of ResR/McdR may be proven effective as adjunctive therapy for shortening the duration of TB treatment.

DNA structural properties of DNA binding sites for 21 transcription factors in the mycobacterial genome

Mycobacterium tuberculosis, the causative agent of tuberculosis, has evolved over time into a multidrug resistance strain that poses a serious global pandemic health threat. The ability to survive and remain dormant within the host macrophage relies on multiple transcription factors contributing to virulence. To date, very limited structural insights from crystallographic and NMR studies are available for TFs and TF-DNA binding events. Understanding the role of DNA structure in TF binding is critical to deciphering MTB pathogenicity and has yet to be resolved at the genome scale. In this work, we analyzed the compositional and conformational preference of 21 mycobacterial TFs, evident at their DNA binding sites, in local and global scales. Results suggest that most TFs prefer binding to genomic regions characterized by unique DNA structural signatures, namely, high electrostatic potential, narrow minor grooves, high propeller twist, helical twist, intrinsic curvature, and DNA rigidity compared to the flanking sequences. Additionally, preference for specific trinucleotide motifs, with clear periodic signals of tetranucleotide motifs, are observed in the vicinity of the TF-DNA interactions. Altogether, our study reports nuanced DNA shape and structural preferences of 21 TFs.

The Mycobacterium tuberculosis protein O-phosphorylation landscape

Bacterial phosphosignalling has been synonymous with two-component systems and their histidine kinases, but many bacteria, including Mycobacterium tuberculosis (Mtb), also code for Ser/Thr protein kinases (STPKs). STPKs are the main phosphosignalling enzymes in eukaryotes but the full extent of phosphorylation on protein Ser/Thr and Tyr (O-phosphorylation) in bacteria is untested. Here we explored the global signalling capacity of the STPKs in Mtb using a panel of STPK loss-of-function and overexpression strains combined with mass spectrometry-based phosphoproteomics. A deep phosphoproteome with >14,000 unique phosphosites shows that O-phosphorylation in Mtb is a vastly underexplored protein modification that affects >80% of the proteome and extensively interfaces with the transcriptional machinery. Mtb O-phosphorylation gives rise to an expansive, distributed and cooperative network of a complexity that has not previously been seen in bacteria and that is on par with eukaryotic phosphosignalling networks. A resource of >3,700 high-confidence direct substrate-STPK interactions and their transcriptional effects provides signalling context for >80% of Mtb proteins and allows the prediction and assembly of signalling pathways for mycobacterial physiology.

Application of combined CRISPR screening for genetic and chemical-genetic interaction profiling in Mycobacterium tuberculosis

CRISPR screening, including CRISPR interference (CRISPRi) and CRISPR-knockout (CRISPR-KO) screening, has become a powerful technology in the genetic screening of eukaryotes. In contrast with eukaryotes, CRISPR-KO screening has not yet been applied to functional genomics studies in bacteria. Here, we constructed genome-scale CRISPR-KO and also CRISPRi libraries in Mycobacterium tuberculosis (Mtb). We first examined these libraries to identify genes essential for Mtb viability. Subsequent screening identified dozens of genes associated with resistance/susceptibility to the antitubercular drug bedaquiline (BDQ). Genetic and chemical validation of the screening results suggested that it provided a valuable resource to investigate mechanisms of action underlying the effects of BDQ and to identify chemical-genetic synergies that can be used to optimize tuberculosis therapy. In summary, our results demonstrate the potential for efficient genome-wide CRISPR-KO screening in bacteria and establish a combined CRISPR screening approach for high-throughput investigation of genetic and chemical-genetic interactions in Mtb.

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.

PrrA modulates Mycobacterium tuberculosis response to multiple environmental cues and is critically regulated by serine/threonine protein kinases

The ability of Mycobacterium tuberculosis (Mtb) to adapt to its surrounding environment is critical for the bacterium to successfully colonize its host. Transcriptional changes are a vital mechanism by which Mtb responds to key environmental signals experienced, such as pH, chloride (Cl^-), nitric oxide (NO), and hypoxia. However, much remains unknown regarding how Mtb coordinates its response to the disparate signals seen during infection. Utilizing a transcription factor (TF) overexpression plasmid library in combination with a pH/Cl^--responsive luciferase reporter, we identified the essential TF, PrrA, part of the PrrAB two-component system, as a TF involved in modulation of Mtb response to pH and Cl^-. Further studies revealed that PrrA also affected Mtb response to NO and hypoxia, with prrA overexpression dampening induction of NO and hypoxia-responsive genes. PrrA is phosphorylated not just by its cognate sensor histidine kinase PrrB, but also by serine/threonine protein kinases (STPKs) at a second distinct site. Strikingly, a STPK-phosphoablative PrrA variant was significantly dampened in its response to NO versus wild type Mtb, disrupted in its ability to adaptively enter a non-replicative state upon extended NO exposure, and attenuated for in vivo colonization. Together, our results reveal PrrA as an important regulator of Mtb response to multiple environmental signals, and uncover a critical role of STPK regulation of PrrA in its function.

Vital to successful host colonization by Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis, is the bacterium’s ability to respond and adapt to changes in its local environment during infection. Here, we discover that the essential transcription factor PrrA, part of the PrrAB two-component system (TCS), modulates Mtb response to four important environmental cues encountered within the host: pH, chloride, nitric oxide, and hypoxia. PrrA acts as a rheostat, adjusting the amplitude of Mtb gene expression changes upon bacterial exposure to each of the four environmental signals. Further, we reveal a critical impact of serine/threonine protein kinases (STPKs) on PrrA function, with prevention of STPK phosphorylation of PrrA disrupting adaptive response of Mtb to growth-inhibiting cues and attenuating the bacterium’s ability to colonize its host. Our work uncovers PrrA as a regulator with broad impact across environmental signals, and highlights how two regulatory systems, TCSs and STPKs, critically interact in coordinating Mtb response to environmental cues.

Expanding the search for small-molecule antibacterials by multidimensional profiling

New techniques for systematic profiling of small-molecule effects can enhance traditional growth inhibition screens for antibiotic discovery and change how we search for new antibacterial agents. Computational models that integrate physicochemical compound properties with their phenotypic and molecular downstream effects can not only predict efficacy of molecules yet to be tested, but also reveal unprecedented insights on compound modes of action (MoAs). The unbiased characterization of compounds that themselves are not growth inhibitory but exhibit diverse MoAs, can expand antibacterial strategies beyond direct inhibition of core essential functions. Early and systematic functional annotation of compound libraries thus paves the way to new models in the selection of lead antimicrobial compounds. In this Review, we discuss how multidimensional small-molecule profiling and the ever-increasing computing power are accelerating the discovery of unconventional antibacterials capable of bypassing resistance and exploiting synergies with established antibacterial treatments and with protective host mechanisms.

Mycobacterium tuberculosis Transcription Factor EmbR Regulates the Expression of Key Virulence Factors That Aid in Ex Vivo and In Vivo Survival

Mycobacterium tuberculosis encodes ~200 transcription factors that modulate gene expression under different microenvironments in the host. Even though high-throughput chromatin immunoprecipitation sequencing and transcriptome sequencing (RNA-seq) studies have identified the regulatory network for ~80% of transcription factors, many transcription factors remain uncharacterized. EmbR is one such transcription factor whose in vivo regulon and biological function are yet to be elucidated. Previous in vitro studies suggested that phosphorylation of EmbR by PknH upregulates the embCAB operon. Using a gene replacement mutant of embR, we investigated its role in modulating cellular morphology, antibiotic resistance, and survival in the host. Contrary to the prevailing hypothesis, under normal growth conditions, EmbR is neither phosphorylated nor impacted by ethambutol resistance through the regulation of the embCAB operon. The embR deletion mutant displayed attenuated M. tuberculosis survival in vivo. RNA-seq analysis suggested that EmbR regulates operons involved in the secretion pathway, lipid metabolism, virulence, and hypoxia, including well-known hypoxia-inducible genes devS and hspX. Lipidome analysis revealed that EmbR modulates levels of all lysophospholipids, several phospholipids, and M. tuberculosis-specific lipids, which is more pronounced under hypoxic conditions. We found that the EmbR mutant is hypersusceptible to hypoxic stress, and RNA sequencing performed under hypoxic conditions indicated that EmbR majorly regulates genes involved in response to acidic pH, hypoxia, and fatty acid metabolism. We observed condition-specific phosphorylation of EmbR, which contributes to EmbR-mediated transcription of several essential genes, ensuring enhanced survival. Collectively, the study establishes EmbR as a key modulator of hypoxic response that facilitates mycobacterial survival in the host.

Multiplexed Strain Phenotyping Defines Consequences of Genetic Diversity in Mycobacterium tuberculosis for Infection and Vaccination Outcomes

There is growing evidence that genetic diversity in Mycobacterium tuberculosis, the causative agent of tuberculosis, contributes to the outcomes of infection and public health interventions, such as vaccination. Epidemiological studies suggest that among the phylogeographic lineages of M. tuberculosis, strains belonging to a sublineage of Lineage 2 (mL2) are associated with concerning clinical features, including hypervirulence, treatment failure, and vaccine escape. The global expansion and increasing prevalence of this sublineage has been attributed to the selective advantage conferred by these characteristics, yet confounding host and environmental factors make it difficult to identify the bacterial determinants driving these associations in human studies. Here, we developed a molecular barcoding strategy to facilitate high-throughput, experimental phenotyping of M. tuberculosis clinical isolates. This approach allowed us to characterize growth dynamics for a panel of genetically diverse M. tuberculosis strains during infection and after vaccination in the mouse model. We found that mL2 strains exhibit distinct growth dynamics in vivo and are resistant to the immune protection conferred by Bacillus Calmette-Guerin (BCG) vaccination. The latter finding corroborates epidemiological observations and demonstrates that mycobacterial features contribute to vaccine efficacy. To investigate the genetic and biological basis of mL2 strains’ distinctive phenotypes, we performed variant analysis, transcriptional studies, and genome-wide transposon sequencing. We identified functional genetic changes across multiple stress and host response pathways in a representative mL2 strain that are associated with variants in regulatory genes. These adaptive changes may underlie the distinct clinical characteristics and epidemiological success of this lineage.

IMPORTANCE Tuberculosis, caused by the bacterium Mycobacterium tuberculosis, is a remarkably heterogeneous disease, a feature that complicates clinical care and public health interventions. The contributions of pathogen genetic diversity to this heterogeneity are uncertain, in part due to the challenges of experimentally manipulating M. tuberculosis, a slow-growing, biosafety level 3 organism. To overcome these challenges, we applied a molecular barcoding strategy to a panel of M. tuberculosis clinical isolates. This novel application of barcoding permitted the high-throughput characterization of M. tuberculosis strain growth dynamics and vaccine resistance in the mouse model of infection. Integrating these results with genomic analyses, we uncover bacterial pathways that contribute to infection outcomes, suggesting targets for improved therapeutics and vaccines.

Multiple genetic paths including massive gene amplification allow Mycobacterium tuberculosis to overcome loss of ESX-3 secretion system substrates

The Mycobacterium tuberculosis (Mtb) ESX-3 type VII secretion system plays a critical role in iron acquisition. Infection of mice with highly attenuated Mtb deletion mutants lacking esxG or esxH, genes encoding key ESX-3 substrates, unexpectedly yielded suppressor mutants with restored capacity to grow in vivo and in vitro in the absence of iron supplementation. Whole-genome sequencing identified two mechanisms of suppression, the disruption of a transcriptional repressor that regulates expression of an ESX-3 paralogous region encoding EsxR and EsxS, and a massive 38- to 60-fold gene amplification of this same region. These data are significant because they reveal a previously unrecognized iron acquisition regulon and inform mechanisms of Mtb chromosome evolution.

Mycobacterium tuberculosis (Mtb) possesses five type VII secretion systems (T7SS), virulence determinants that include the secretion apparatus and associated secretion substrates. Mtb strains deleted for the genes encoding substrates of the ESX-3 T7SS, esxG or esxH, require iron supplementation for in vitro growth and are highly attenuated in vivo. In a subset of infected mice, suppressor mutants of esxG or esxH deletions were isolated, which enabled growth to high titers or restored virulence. Suppression was conferred by mechanisms that cause overexpression of an ESX-3 paralogous region that lacks genes for the secretion apparatus but encodes EsxR and EsxS, apparent ESX-3 orphan substrates that functionally compensate for the lack of EsxG or EsxH. The mechanisms include the disruption of a transcriptional repressor and a massive 38- to 60-fold gene amplification. These data identify an iron acquisition regulon, provide insight into T7SS, and reveal a mechanism of Mtb chromosome evolution involving “accordion-type” amplification.

MadR mediates acyl CoA-dependent regulation of mycolic acid desaturation in mycobacteria

Our studies show that the mycolic acid desaturase regulator (MadR) acts as a molecular switch, controlling the desaturation and biosynthesis of mycolic acids, key lipids of the cell envelopes of mycobacteria. MadR works by a distinct mechanism wherein it binds various acyl-coenzyme As (aceyl-CoAs), but only saturated acyl-CoAs relieve DNA binding and repression. This suggests a unique mechanism that involves sensing of acyl-CoA pools as a checkpoint for coordinating mycolic acid remodeling and biosynthesis in response to cell surface perturbation. Our findings further our understanding of how mycobacteria control cell wall composition in response to stress across various environments ranging from soil to an intracellular niche in infected macrophages, with implications for understanding strategies for pathogenesis in the tubercle bacillus.

Mycobacterium tuberculosis has a lipid-rich cell envelope that is remodeled throughout infection to enable adaptation within the host. Few transcriptional regulators have been characterized that coordinate synthesis of mycolic acids, the major cell wall lipids of mycobacteria. Here, we show that the mycolic acid desaturase regulator (MadR), a transcriptional repressor of the mycolate desaturase genes desA1 and desA2, controls mycolic acid desaturation and biosynthesis in response to cell envelope stress. A madR-null mutant of M. smegmatis exhibited traits of an impaired cell wall with an altered outer mycomembrane, accumulation of a desaturated α-mycolate, susceptibility to antimycobacterials, and cell surface disruption. Transcriptomic profiling showed that enriched lipid metabolism genes that were significantly down-regulated upon madR deletion included acyl-coenzyme A (aceyl-CoA) dehydrogenases, implicating it in the indirect control of β-oxidation pathways. Electromobility shift assays and binding affinities suggest a unique acyl-CoA pool–sensing mechanism, whereby MadR is able to bind a range of acyl-CoAs, including those with unsaturated as well as saturated acyl chains. MadR repression of desA1/desA2 is relieved upon binding of saturated acyl-CoAs of chain length C₁₆ to C₂₄, while no impact is observed upon binding of shorter chain and unsaturated acyl-CoAs. We propose this mechanism of regulation as distinct to other mycolic acid and fatty acid synthesis regulators and place MadR as the key regulatory checkpoint that coordinates mycolic acid remodeling during infection in response to host-derived cell surface perturbation.

Proteome Profiling of Mycobacterium tuberculosis Cells Exposed to Nitrosative Stress

Reactive nitrogen species (RNS) are secreted by human cells in response to infection by Mycobacterium tuberculosis (Mtb). Although RNS can kill Mtb under some circumstances, Mtb can adapt and survive in the presence of RNS by a process that involves modulation of gene expression. Previous studies focused primarily on stress-related changes in the Mtb transcriptome. This study unveils changes in the Mtb proteome in response to a sub-lethal dose of nitric oxide (NO) over several hours of exposure. Proteins were identified using liquid chromatography coupled with electrospray ionization mass spectrometry (LC-MS/MS). A total of 2911 Mtb proteins were identified, of which 581 were differentially abundant (DA) after exposure to NO in at least one of the four time points (30 min, 2 h, 6 h, and 20 h). The proteomic response to NO was marked by two phases, with few DA proteins in the early phase and a multitude of DA proteins in the later phase. The efflux pump Rv1687 stood out as being the only protein more abundant at all the time points and might play a role in the early protection of Mtb against nitrosative stress. These changes appeared to be compensatory in nature, contributing to iron homeostasis, energy metabolism, and other stress responses. This study thereby provides new insights into the response of Mtb to NO at the level of proteomics.

Quantitative prediction of conditional vulnerabilities in regulatory and metabolic networks using PRIME

The ability of Mycobacterium tuberculosis (Mtb) to adopt heterogeneous physiological states underlies its success in evading the immune system and tolerating antibiotic killing. Drug tolerant phenotypes are a major reason why the tuberculosis (TB) mortality rate is so high, with over 1.8 million deaths annually. To develop new TB therapeutics that better treat the infection (faster and more completely), a systems-level approach is needed to reveal the complexity of network-based adaptations of Mtb. Here, we report a new predictive model called PRIME (Phenotype of Regulatory influences Integrated with Metabolism and Environment) to uncover environment-specific vulnerabilities within the regulatory and metabolic networks of Mtb. Through extensive performance evaluations using genome-wide fitness screens, we demonstrate that PRIME makes mechanistically accurate predictions of context-specific vulnerabilities within the integrated regulatory and metabolic networks of Mtb, accurately rank-ordering targets for potentiating treatment with frontline drugs.

Multiform antimicrobial resistance from a metabolic mutation

A critical challenge for microbiology and medicine is how to cure infections by bacteria that survive antibiotic treatment by persistence or tolerance. Seeking mechanisms behind such high survival, we developed a forward-genetic method for efficient isolation of high-survival mutants in any culturable bacterial species. We found that perturbation of an essential biosynthetic pathway (arginine biosynthesis) in a mycobacterium generated three distinct forms of resistance to diverse antibiotics, each mediated by induction of WhiB7: high persistence and tolerance to kanamycin, high survival upon exposure to rifampicin, and minimum inhibitory concentration-shifted resistance to clarithromycin. As little as one base change in a gene that encodes, a metabolic pathway component conferred multiple forms of resistance to multiple antibiotics with different targets. This extraordinary resilience may help explain how substerilizing exposure to one antibiotic in a regimen can induce resistance to others and invites development of drugs targeting the mediator of multiform resistance, WhiB7.

Intricate Genetic Programs Controlling Dormancy in Mycobacterium tuberculosis

Mycobacterium tuberculosis (MTB) displays the remarkable ability to transition in and out of dormancy, a hallmark of the pathogen’s capacity to evade the immune system and exploit susceptible individuals. Uncovering the gene regulatory programs that underlie the phenotypic shifts in MTB during disease latency and reactivation has posed a challenge. We develop an experimental system to precisely control dissolved oxygen levels in MTB cultures in order to capture the transcriptional events that unfold as MTB transitions into and out of hypoxia-induced dormancy. Using a comprehensive genome-wide transcription factor binding map and insights from network topology analysis, we identify regulatory circuits that deterministically drive sequential transitions across six transcriptionally and functionally distinct states encompassing more than three-fifths of the MTB genome. The architecture of the genetic programs explains the transcriptional dynamics underlying synchronous entry of cells into a dormant state that is primed to infect the host upon encountering favorable conditions.

Mycobacterium tuberculosis (MTB) persists within the host by counteracting disparate stressors including hypoxia. Peterson et al. report a transcriptional program that coordinates sequential state transitions to drive MTB in and out of hypoxia-induced dormancy. Among varied properties, this program encodes advanced preparedness to infect the host in favorable conditions.

Dissecting the RecA-(In)dependent Response to Mitomycin C in Mycobacterium tuberculosis Using Transcriptional Profiling and Proteomics Analyses

Mycobacteria exploit at least two independent global systems in response to DNA damage: the LexA/RecA-dependent SOS response and the PafBC-regulated pathway. Intracellular pathogens, such as Mycobacterium tuberculosis, are exposed to oxidative and nitrosative stress during the course of infection while residing inside host macrophages. The current understanding of RecA-independent responses to DNA damage is based on the saprophytic model of Mycobacterium smegmatis, a free-living and nonpathogenic mycobacterium. The aim of the present study was to identify elements of RecA-independent responses to DNA damage in pathogenic intracellular mycobacteria. With the help of global transcriptional profiling, we were able to dissect RecA-dependent and RecA-independent pathways. We profiled the DNA damage responses of an M. tuberculosis strain lacking the recA gene, a strain with an undetectable level of the PafBC regulatory system, and a strain with both systems tuned down simultaneously. RNA-Seq profiling was correlated with the evaluation of cell survival in response to DNA damage to estimate the relevance of each system to the overall sensitivity to genotoxic agents. We also carried out whole-cell proteomics analysis of the M. tuberculosis strains in response to mitomycin C. This approach highlighted that LexA, a well-defined key element of the SOS system, is proteolytically inactivated during RecA-dependent DNA repair, which we found to be transcriptionally repressed in response to DNA-damaging agents in the absence of RecA. Proteomics profiling revealed that AlkB was significantly overproduced in the ΔrecA pafBC^{CRISPRi/dCas9} strain and that Holliday junction resolvase RuvX was a DNA damage response factor that was significantly upregulated regardless of the presence of functional RecA and PafBC systems, thus falling into a third category of DNA damage factors: RecA- and PafBC-independent. While invisible to the mass spectrometer, the genes encoding alkA, dnaB, and dnaE2 were significantly overexpressed in the ΔrecA pafBC^{CRISPRi/dCas9} strain at the transcript level.

Single START-domain protein Mtsp17 is involved in transcriptional regulation in Mycobacterium smegmatis

Tuberculosis caused by the pathogen Mycobacterium tuberculosis (MTB), remains a significant threat to global health. Elucidating the mechanisms of essential MTB genes provides an important theoretical basis for drug exploitation. Gene mtsp17 is essential and is conserved in the Mycobacterium genus. Although Mtsp17 has a structure closely resembling typical steroidogenic acute regulatory protein-related lipid transfer (START) family proteins, its biological function is different. This study characterizes the transcriptomes of Mycobacterium smegmatis to explore the consequences of mtsp17 downregulation on gene expression. Suppression of the mtsp17 gene resulted in significant down-regulation of 3% and upregulation of 1% of all protein-coding genes. Expression of desA1, an essential gene involved in mycolic acid synthesis, and the anti-SigF antagonist MSMEG_0586 were down-regulated in the conditional Mtsp17 knockout mutant and up-regulated in the Mtsp17 over-expression strain. Trends in the changes of 70 of the 79 differentially expressed genes (Log2 fold change > 1.5) in the conditional Mtsp17 knockout strain were the same as in the SigF knockout strain. Our data suggest that Mtsp17 is likely an activator of desA1 and Mtsp17 regulates the SigF regulon by SigF regulatory pathways through the anti-SigF antagonist MSMEG_0586. Our findings indicate the role of Mtsp17 may be in transcriptional regulation, provide new insights into the molecular mechanisms of START family proteins, and uncover a new node in the regulatory network of mycobacteria.

A Role for Mycobacterium tuberculosis Sigma Factor C in Copper Nutritional Immunity

Sigma factor C (SigC) contributes to Mycobacterium tuberculosis virulence in various animal models, but the stress response coordinated by this transcription factor was undefined. The results presented here indicate that SigC prevents copper starvation. Whole genome expression studies demonstrate short-term (4-h) induction of sigC, controlled from a tetracycline-inducible promoter, upregulates ctpB and genes in the nonribosomal peptide synthase (nrp) operon. These genes are expressed at higher levels after 48-h sigC induction, but also elevated are genes encoding copper-responsive regulator RicR and RicR-regulated copper toxicity response operon genes rv0846-rv0850, suggesting prolonged sigC induction results in excessive copper uptake. No growth and global transcriptional differences are observed between a sigC null mutant relative to its parent strain in 7H9 medium. In a copper-deficient medium, however, growth of the sigC deletion strain lags the parent, and 40 genes (including those in the nrp operon) are differentially expressed. Copper supplementation reverses the growth defect and silences most transcriptional differences. Together, these data support SigC as a transcriptional regulator of copper acquisition when the metal is scarce. Attenuation of sigC mutants in severe combined immunodeficient mice is consistent with an inability to overcome innate host defenses that sequester copper ions to deprive invading microbes of this essential micronutrient.

A snapshot of translation in Mycobacterium tuberculosis during exponential growth and nutrient starvation revealed by ribosome profiling

Mycobacterium tuberculosis, which causes tuberculosis, can undergo prolonged periods of non-replicating persistence in the host. The mechanisms underlying this are not fully understood, but translational regulation is thought to play a role. A large proportion of mRNA transcripts expressed in M. tuberculosis lack canonical bacterial translation initiation signals, but little is known about the implications of this for fine-tuning of translation. Here, we perform ribosome profiling to characterize the translational landscape of M. tuberculosis under conditions of exponential growth and nutrient starvation. Our data reveal robust, widespread translation of non-canonical transcripts and point toward different translation initiation mechanisms compared to canonical Shine-Dalgarno transcripts. During nutrient starvation, patterns of ribosome recruitment vary, suggesting that regulation of translation in this pathogen is more complex than originally thought. Our data represent a rich resource for others seeking to understand translational regulation in bacterial pathogens.

Human tuberculosis and Mycobacterium tuberculosis complex: A review on genetic diversity, pathogenesis and omics approaches in host biomarkers discovery

Mycobacterium tuberculosis complex (MTBC) refers to a group of mycobacteria encompassing nine members of closely related species that causes tuberculosis in animals and humans. Among the nine members, Mycobacterium tuberculosis (M. tuberculosis) remains the main causative agent for human tuberculosis that results in high mortality and morbidity globally. In general, MTBC species are low in diversity but exhibit distinctive biological differences and phenotypes among different MTBC lineages. MTBC species are likely to have evolved from a common ancestor through insertions/deletions processes resulting in species speciation with different degrees of pathogenicity. The pathogenesis of human tuberculosis is complex and remains poorly understood. It involves multi-interactions or evolutionary co-options between host factors and bacterial determinants for survival of the MTBC. Granuloma formation as a protection or survival mechanism in hosts by MTBC remains controversial. Additionally, MTBC species are capable of modulating host immune response and have adopted several mechanisms to evade from host immune attack in order to survive in humans. On the other hand, current diagnostic tools for human tuberculosis are inadequate and have several shortcomings. Numerous studies have suggested the potential of host biomarkers in early diagnosis of tuberculosis, in disease differentiation and in treatment monitoring. "Multi-omics" approaches provide holistic views to dissect the association of MTBC species with humans and offer great advantages in host biomarkers discovery. Thus, in this review, we seek to understand how the genetic variations in MTBC lead to species speciation with different pathogenicity. Furthermore, we also discuss how the host and bacterial players contribute to the pathogenesis of human tuberculosis. Lastly, we provide an overview of the journey of "omics" approaches in host biomarkers discovery in human tuberculosis and provide some interesting insights on the challenges and directions of "omics" approaches in host biomarkers innovation and clinical implementation.

Ultra-low Dose Aerosol Infection of Mice with Mycobacterium tuberculosis More Closely Models Human Tuberculosis

Tuberculosis (TB) is a heterogeneous disease manifesting in a subset of individuals infected with aerosolized Mycobacterium tuberculosis (Mtb). Unlike human TB, murine infection results in uniformly high lung bacterial burdens and poorly organized granulomas. To develop a TB model that more closely resembles human disease, we infected mice with an ultra-low dose (ULD) of between 1-3 founding bacteria, reflecting a physiologic inoculum. ULD-infected mice exhibited highly heterogeneous bacterial burdens, well-circumscribed granulomas that shared features with human granulomas, and prolonged Mtb containment with unilateral pulmonary infection in some mice. We identified blood RNA signatures in mice infected with an ULD or a conventional Mtb dose (50-100 CFU) that correlated with lung bacterial burdens and predicted Mtb infection outcomes across species, including risk of progression to active TB in humans. Overall, these findings highlight the potential of the murine TB model and show that ULD infection recapitulates key features of human TB.

The genetic proteome: Using genetics to inform the proteome of mycobacterial pathogens

Mycobacterial pathogens pose a sustained threat to human health. There is a critical need for new diagnostics, therapeutics, and vaccines targeting both tuberculous and nontuberculous mycobacterial species. Understanding the basic mechanisms used by diverse mycobacterial species to cause disease will facilitate efforts to design new approaches toward detection, treatment, and prevention of mycobacterial disease. Molecular, genetic, and biochemical approaches have been widely employed to define fundamental aspects of mycobacterial physiology and virulence. The recent expansion of genetic tools in mycobacteria has further increased the accessibility of forward genetic approaches. Proteomics has also emerged as a powerful approach to further our understanding of diverse mycobacterial species. Detection of large numbers of proteins and their modifications from complex mixtures of mycobacterial proteins is now routine, with efforts of quantification of these datasets becoming more robust. In this review, we discuss the “genetic proteome,” how the power of genetics, molecular biology, and biochemistry informs and amplifies the quality of subsequent analytical approaches and maximizes the potential of hypothesis-driven mycobacterial research. Published proteomics datasets can be used for hypothesis generation and effective post hoc supplementation to experimental data. Overall, we highlight how the integration of proteomics, genetic, molecular, and biochemical approaches can be employed successfully to define fundamental aspects of mycobacterial pathobiology.

DNA Isolation from Mycobacteria

A vast array of molecular biology tools have been developed to investigate the Mycobacterium tuberculosis genome since the advent of its successful sequencing in 1998. These tools, such as quantitative and end point polymerase chain reaction, chromatin immunoprecipitation, and whole genome sequencing, require genomic DNA extracted from lysed mycobacteria. There are numerous methods described in the literature using mechanical, enzymatic, or chemical means to lyse cells and extract genomic DNA to varying degrees of purity. Here, we describe appropriate methods for genomic DNA isolation from solid or liquid cultures from both M. tuberculosis and nontuberculous mycobacteria.

Antitoxin autoregulation of M. tuberculosis toxin-antitoxin expression through negative cooperativity arising from multiple inverted repeat sequences

Toxin-antitoxin systems play key roles in bacterial adaptation, including protection from antibiotic assault and infection by bacteriophages. The type IV toxin-antitoxin system AbiE encodes a DUF1814 nucleotidyltransferase-like toxin, and a two-domain antitoxin. In Streptococcus agalactiae, the antitoxin AbiEi negatively autoregulates abiE expression through positively co-operative binding to inverted repeats within the promoter. The human pathogen Mycobacterium tuberculosis encodes four DUF1814 putative toxins, two of which have antitoxins homologous to AbiEi. One such M. tuberculosis antitoxin, named Rv2827c, is required for growth and whilst the structure has previously been solved, the mode of regulation is unknown. To complete the gaps in our understanding, we first solved the structure of S. agalactiae AbiEi to 1.83 Å resolution for comparison with M. tuberculosis Rv2827c. AbiEi contains an N-terminal DNA binding domain and C-terminal antitoxicity domain, with bilateral faces of opposing charge. The overall AbiEi fold is similar to Rv2827c, though smaller, and with a 65° difference in C-terminal domain orientation. We further demonstrate that, like AbiEi, Rv2827c can autoregulate toxin-antitoxin operon expression. In contrast with AbiEi, the Prv2827c promoter contains two sets of inverted repeats, which bind Rv2827c with differing affinities depending on the sequence consensus. Surprisingly, Rv2827c bound with negative co-operativity to the full Prv2827c promoter, demonstrating an unexpectedly complex form of transcriptional regulation.

Metabolic Fate of Human Immunoactive Sterols in Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) infection is among top ten causes of death worldwide, and the number of drug-resistant strains is increasing. The direct interception of human immune signaling molecules by Mtb remains elusive, limiting drug discovery. Oxysterols and secosteroids regulate both innate and adaptive immune responses. Here we report a functional, structural, and bioinformatics study of Mtb enzymes initiating cholesterol catabolism and demonstrated their interrelation with human immunity. We show that these enzymes metabolize human immune oxysterol messengers. Rv2266 - the most potent among them - can also metabolize vitamin D3 (VD3) derivatives. High-resolution structures show common patterns of sterols binding and reveal a site for oxidative attack during catalysis. Finally, we designed a compound that binds and inhibits three studied proteins. The compound shows activity against Mtb H37Rv residing in macrophages. Our findings contribute to molecular understanding of suppression of immunity and suggest that Mtb has its own transformation system resembling the human phase I drug-metabolizing system.

The Expanding Computational Toolbox for Engineering Microbial Phenotypes at the Genome Scale

Microbial strains are being engineered for an increasingly diverse array of applications, from chemical production to human health. While traditional engineering disciplines are driven by predictive design tools, these tools have been difficult to build for biological design due to the complexity of biological systems and many unknowns of their quantitative behavior. However, due to many recent advances, the gap between design in biology and other engineering fields is closing. In this work, we discuss promising areas of development of computational tools for engineering microbial strains. We define five frontiers of active research: (1) Constraint-based modeling and metabolic network reconstruction, (2) Kinetics and thermodynamic modeling, (3) Protein structure analysis, (4) Genome sequence analysis, and (5) Regulatory network analysis. Experimental and machine learning drivers have enabled these methods to improve by leaps and bounds in both scope and accuracy. Modern strain design projects will require these tools to be comprehensively applied to the entire cell and efficiently integrated within a single workflow. We expect that these frontiers, enabled by the ongoing revolution of big data science, will drive forward more advanced and powerful strain engineering strategies.

Transcriptional regulator-induced phenotype screen reveals drug potentiators in Mycobacterium tuberculosis

Transposon-based strategies provide a powerful and unbiased way to study bacterial stress response^1–8, but these approaches cannot fully capture the complexities of network-based behavior. Here, we present a network-based genetic screening approach: the Transcriptional Regulator Induced Phenotype (TRIP) screen, which we used to identify previously uncharacterized network adaptations of Mycobacterium tuberculosis (Mtb) to the first-line anti-TB drug isoniazid (INH). We found regulators that alter INH susceptibility when induced, several of which could not be identified by standard gene disruption approaches. We then focused on a specific regulator, mce3R, which potentiated INH activity when induced. We compared mce3R-regulated genes with baseline INH transcriptional responses and implicated the gene ctpD (Rv1469) as a putative INH effector. Evaluating a ctpD disruption mutant demonstrated a previously unknown role for this gene in INH susceptibility. Integrating TRIP screening with network information can uncover sophisticated molecular response programs.

Arrayed CRISPRi and quantitative imaging describe the morphotypic landscape of essential mycobacterial genes

Mycobacterium tuberculosis possesses a large number of genes of unknown or predicted function, undermining fundamental understanding of pathogenicity and drug susceptibility. To address this challenge, we developed a high-throughput functional genomics approach combining inducible CRISPR-interference and image-based analyses of morphological features and sub-cellular chromosomal localizations in the related non-pathogen, M. smegmatis. Applying automated imaging and analysis to 263 essential gene knockdown mutants in an arrayed library, we derive robust, quantitative descriptions of bacillary morphologies consequent on gene silencing. Leveraging statistical-learning, we demonstrate that functionally related genes cluster by morphotypic similarity and that this information can be used to inform investigations of gene function. Exploiting this observation, we infer the existence of a mycobacterial restriction-modification system, and identify filamentation as a defining mycobacterial response to histidine starvation. Our results support the application of large-scale image-based analyses for mycobacterial functional genomics, simultaneously establishing the utility of this approach for drug mechanism-of-action studies.

Cell-Cycle-Associated Expression Patterns Predict Gene Function in Mycobacteria

Although the major events in prokaryotic cell cycle progression are likely to be coordinated with transcriptional and metabolic changes, these processes remain poorly characterized. Unlike many rapidly growing bacteria, DNA replication and cell division are temporally resolved in mycobacteria, making these slow-growing organisms a potentially useful system to investigate the prokaryotic cell cycle. To determine whether cell-cycle-dependent gene regulation occurs in mycobacteria, we characterized the temporal changes in the transcriptome of synchronously replicating populations of Mycobacterium tuberculosis (Mtb). By enriching for genes that display a sinusoidal expression pattern, we discover 485 genes that oscillate with a period consistent with the cell cycle. During cytokinesis, the timing of gene induction could be used to predict the timing of gene function, as mRNA abundance was found to correlate with the order in which proteins were recruited to the developing septum. Similarly, the expression pattern of primary metabolic genes could be used to predict the relative importance of these pathways for different cell cycle processes. Pyrimidine synthetic genes peaked during DNA replication, and their depletion caused a filamentation phenotype that phenocopied defects in this process. In contrast, the inosine monophasphate dehydrogenase dedicated to guanosine synthesis, GuaB2, displayed the opposite expression pattern and its depletion perturbed septation. Together, these data imply obligate coordination between primary metabolism and cell division and identify periodically regulated genes that can be related to specific cell biological functions.

RegX3 Activates whiB3 Under Acid Stress and Subverts Lysosomal Trafficking of Mycobacterium tuberculosis in a WhiB3-Dependent Manner

Two-component systems (TCSs) are central to the ability of Mycobacterium tuberculosis to respond to stress. One such paired TCS is SenX3-RegX3, which responds to phosphate starvation. Here we show that RegX3 is required for M. tuberculosis to withstand low pH, one of the challenges encountered by the bacterium in the host environment, and that RegX3 activates the cytosolic redox sensor WhiB3 to launch an appropriate response to acid stress. We show that the whiB3 promoter of M. tuberculosis harbors a RegX3 binding motif. Electrophoretic mobility shift assays (EMSAs) show that phosphorylated RegX3 (RegX3-P) (but not its unphosphorylated counterpart) binds to this motif, whereas a DNA binding mutant, RegX3 (K204A) fails to do so. Mutation of the putative RegX3 binding motif on the whiB3 promoter, abrogates the binding of RegX3-P. The significance of this binding is established by demonstrating that the expression of whiB3 is significantly attenuated under phosphate starvation or under acid stress in the regX3-inactivated mutant, ΔregX3. Green fluorescent protein (GFP)-based reporter assays further confirm the requirement of RegX3 for the activation of the whiB3 promoter. The compromised survival of ΔregX3 under acid stress and its increased trafficking to the lysosomal compartment are reversed upon complementation with either regX3 or whiB3, suggesting that RegX3 exerts its effects in a WhiB3-dependent manner. Finally, using an in vitro granuloma model, we show that granuloma formation is compromised in the absence of regX3, but restored upon complementation with either regX3 or whiB3. Our findings provide insight into an important role of RegX3 in the network that regulates the survival of M. tuberculosis under acid stress similar to that encountered in its intracellular niche. Our results argue strongly in favor of a role of the RegX3-WhiB3 axis in establishment of M. tuberculosis infection.

Redefining genetic essentiality in Mycobacterium tuberculosis

Sequencing transposon mutant libraries have been pivotal in annotating essential and non-essential genes in bacteria. This is particularly very helpful in the case of Mycobacterium tuberculosis with a large part of its genome without known function. It is not known whether there are any variations in the essentiality states as a function of optimal growth in the absence of any selection pressure. We here grow a high-density mutant library of M. tuberculosis through serial cultures and monitor the temporal fluctuations in insertion frequencies across all TA dinucleotides in the genome. Genes that cause morphological and physiological heterogeneity or enable metabolic bypass were found to gradually lose insertions, while genes comprising the toxin-antitoxin systems were found to get enriched with insertions during growth in nutrient replete conditions. High levels of fluctuations were observed in genes involved in cell wall and cell processes, intermediary metabolism, and genes involved in virulence, suggesting new modes of adaptation undertaken by the mutants. We also report the essentiality status of several newly annotated genetic features.

Entropy of a bacterial stress response is a generalizable predictor for fitness and antibiotic sensitivity

Current approaches explore bacterial genes that change transcriptionally upon stress exposure as diagnostics to predict antibiotic sensitivity. However, transcriptional changes are often specific to a species or antibiotic, limiting implementation to known settings only. While a generalizable approach, predicting bacterial fitness independent of strain, species or type of stress, would eliminate such limitations, it is unclear whether a stress-response can be universally captured. By generating a multi-stress and species RNA-Seq and experimental evolution dataset, we highlight the strengths and limitations of existing gene-panel based methods. Subsequently, we build a generalizable method around the observation that global transcriptional disorder seems to be a common, low-fitness, stress response. We quantify this disorder using entropy, which is a specific measure of randomness, and find that in low fitness cases increasing entropy and transcriptional disorder results from a loss of regulatory gene-dependencies. Using entropy as a single feature, we show that fitness and quantitative antibiotic sensitivity predictions can be made that generalize well beyond training data. Furthermore, we validate entropy-based predictions in 7 species under antibiotic and non-antibiotic conditions. By demonstrating the feasibility of universal predictions of bacterial fitness, this work establishes the fundamentals for potentially new approaches in infectious disease diagnostics.

Ectopic Expression of Rv0023 Mediates Isoniazid/Ethionamide Tolerance via Altering NADH/NAD+ Levels in Mycobacterium smegmatis

Tuberculosis (TB) caused by Mycobacterium tuberculosis (Mtb) accounts for nearly 1.2 million deaths per annum worldwide. Due to the emergence of multidrug-resistant (MDR) Mtb strains, TB, a curable and avertable disease, remains one of the leading causes of morbidity and mortality. Isoniazid (INH) is a first-line anti-TB drug while ethionamide (ETH) is used as a second-line anti-TB drug. INH and ETH resistance develop through a network of genes involved in various biosynthetic pathways. In this study, we identified Rv0023, an Mtb protein belonging to the xenobiotic response element (XRE) family of transcription regulators, which has a role in generating higher tolerance toward INH and ETH in Mycobacterium smegmatis (Msmeg). Overexpression of Rv0023 in Msmeg leads to the development of INH- and ETH-tolerant strains. The strains expressing Rv0023 have a higher ratio of NADH/NAD⁺, and this physiological event is known to play a crucial role in the development of INH/ETH co-resistance in Msmeg. Gene expression analysis of some target genes revealed reduction in the expression of the ndh gene, but no direct interaction was observed between Rv0023 and the ndh promoter region. Rv0023 is divergently expressed to Rv0022c (whiB5) and we observed a direct interaction between the recombinant Rv0023 protein with the upstream region of Rv0022c, confirmed using reporter constructs of Msmeg. However, we found no indication that this interaction might play a role in the development of INH/ETH drug tolerance.

Evaluation of specificity determinants in Mycobacterium tuberculosis σ/anti-σ factor interactions

Extra Cytoplasmic Function (ECF) σ factor/regulatory protein (anti-σ factor) pairs govern environment mediated changes in gene expression in bacteria. The release of the ECF σ factor from an inactive σ/anti-σ factor complex is triggered by specific environmental stimuli. The free σ factor then associates with the RNA polymerase and drives the expression of genes in its target regulon. Multiple ECF σ/anti-σ pairs ensure calibrated changes in the expression profile by correlating diverse environmental stimuli with changes in the intracellular levels of different ECF σ factors. Specificity in σ/anti-σ factor interaction is thus essential for accurate signal transduction. Here we describe experiments to evaluate interactions between different M. tuberculosis σ and anti-σ proteins in vitro. The interaction parameters suggest that cross-talk between non-cognate σ/anti-σ pairs is likely. The sequence and conformational determinants that govern interaction specificity in a σ/anti-σ complex are not immediately evident due to substantial structural conservation. Sequence-structure analysis of all σ/anti-σ pairs suggest that conserved residues are not the primary determinants of σ/anti-σ interactions-a finding that suggests a potential route to set tolerance limits in interaction specificity. Non-specific σ/anti-σ interactions are likely to be biologically significant as it can contribute to heterogeneity in cellular responses in a bacterial population under less stringent requirements. This finding is relevant for synthetic biology approaches to engineer bacteria using σ/anti-σ transcription initiation modules for diverse applications in biotechnology.

Protein kinase B controls Mycobacterium tuberculosis growth via phosphorylation of the transcriptional regulator Lsr2 at threonine 112

Mycobacterium tuberculosis (Mtb) is able to persist in the body through months of multi‐drug therapy. Mycobacteria possess a wide range of regulatory proteins, including the protein kinase B (PknB) which controls peptidoglycan biosynthesis during growth. Here, we observed that depletion of PknB resulted in specific transcriptional changes that are likely caused by reduced phosphorylation of the H‐NS‐like regulator Lsr2 at threonine 112. The activity of PknB towards this phosphosite was confirmed with purified proteins, and this site was required for adaptation of Mtb to hypoxic conditions, and growth on solid media. Like H‐NS, Lsr2 binds DNA in sequence‐dependent and non‐specific modes. PknB phosphorylation of Lsr2 reduced DNA binding, measured by fluorescence anisotropy and electrophoretic mobility shift assays, and our NMR structure of phosphomimetic T112D Lsr2 suggests that this may be due to increased dynamics of the DNA‐binding domain. Conversely, the phosphoablative T112A Lsr2 had increased binding to certain DNA sites in ChIP‐sequencing, and Mtb containing this variant showed transcriptional changes that correspond with the change in DNA binding. In summary, PknB controls Mtb growth and adaptations to the changing host environment by phosphorylating the global transcriptional regulator Lsr2.

Phase variation in Mycobacterium tuberculosis glpK produces transiently heritable drug tolerance

The length and complexity of tuberculosis (TB) therapy, as well as the propensity of Mycobacterium tuberculosis to develop drug resistance, are major barriers to global TB control efforts. M. tuberculosis is known to have the ability to enter into a drug-tolerant state, which may explain many of these impediments to TB treatment. We have identified a mechanism of genetically encoded but rapidly reversible drug tolerance in M. tuberculosis caused by transient frameshift mutations in a homopolymeric tract (HT) of 7 cytosines (7C) in the glpK gene. Inactivating frameshift mutations associated with the 7C HT in glpK produce small colonies that exhibit heritable multidrug increases in minimal inhibitory concentrations and decreases in drug-dependent killing; however, reversion back to a fully drug-susceptible large-colony phenotype occurs rapidly through the introduction of additional insertions or deletions in the same glpK HT region. These reversible frameshift mutations in the 7C HT of M. tuberculosis glpK occur in clinical isolates, accumulate in M. tuberculosis-infected mice with further accumulation during drug treatment, and exhibit a reversible transcriptional profile including induction of dosR and sigH and repression of kstR regulons, similar to that observed in other in vitro models of M. tuberculosis tolerance. These results suggest that GlpK phase variation may contribute to drug tolerance, treatment failure, and relapse in human TB. Drugs effective against phase-variant M. tuberculosis may hasten TB treatment and improve cure rates.

Cross-talk between the three furA orthologs in Mycobacterium smegmatis and the contribution to isoniazid resistance

The ferric uptake regulator A (FurA) plays an essential role in responding to oxidative stress in mycobacteria. The genome of Mycobacterium smegmatis harbours three FurA orthologs; however, the potential cross-talk and contribution to drug resistance of different furA operon remain underdetermined. In this study, we characterized the cross-regulation and effect in drug resistance of these orthologs from M. smegmatis. Cross-binding of FurA protein to furA promoter was observed. The binding of FurA1 to furA3p and FurA2 to furA1p or furA3p is even more pronounced than their self-binding. The three FurA proteins are all functional at repressing the expression of the peroxidase enzyme katG1/katG2 in vivo. When overexpressing any of the furA orthologs in M. smegmatis, the bacteria become more resistant to isoniazid (INH). This pattern is consistent with that in Mycobacterium bovis. However, the knockdown of furA does not affect the INH sensitivity. This is the first report of cross-talk and contribution to drug resistance of all three furA orthologs in M. smegmatis.

Genome-wide mutational biases fuel transcriptional diversity in the Mycobacterium tuberculosis complex

The Mycobacterium tuberculosis complex (MTBC) members display different host-specificities and virulence phenotypes. Here, we have performed a comprehensive RNAseq and methylome analysis of the main clades of the MTBC and discovered unique transcriptional profiles. The majority of genes differentially expressed between the clades encode proteins involved in host interaction and metabolic functions. A significant fraction of changes in gene expression can be explained by positive selection on single mutations that either create or disrupt transcriptional start sites (TSS). Furthermore, we show that clinical strains have different methyltransferases inactivated and thus different methylation patterns. Under the tested conditions, differential methylation has a minor direct role on transcriptomic differences between strains. However, disruption of a methyltransferase in one clinical strain revealed important expression differences suggesting indirect mechanisms of expression regulation. Our study demonstrates that variation in transcriptional profiles are mainly due to TSS mutations and have likely evolved due to differences in host characteristics.

The crystal structure of AcrR from Mycobacterium tuberculosis reveals a one-component transcriptional regulation mechanism

Transcriptional regulator proteins are closely involved in essential survival strategies in bacteria. AcrR is a one-component allosteric repressor of the genes associated with lipid transport and antibiotic resistance. When fatty acid ligands bind to the C-terminal ligand-binding cavity of AcrR, a conformational change in the N-terminal operator-binding region of AcrR is triggered, which releases the repressed DNA and initiates transcription. This paper focuses on the structural transition mechanism of AcrR of Mycobacterium tuberculosis upon DNA and ligand binding. AcrR loses its structural integrity upon ligand-mediated structural alteration and bends toward the promoter DNA in a more compact form, initiating a rotational motion. Our functional characterization of AcrR and description of the ligand- and DNA-recognition mechanism may facilitate the discovery of new therapies for tuberculosis.