Genome-wide gene expression tuning reveals diverse vulnerabilities of M. tuberculosis.

Antibacterial agents target the products of essential genes but rarely achieve complete target inhibition. Thus, the all-or-none definition of essentiality afforded by traditional genetic approaches fails to discern the most attractive bacterial targets: those whose incomplete inhibition results in major fitness costs. In contrast, gene "vulnerability" is a continuous, quantifiable trait that relates the magnitude of gene inhibition to the effect on bacterial fitness. We developed a CRISPR interference-based functional genomics method to systematically titrate gene expression in Mycobacterium tuberculosis (Mtb) and monitor fitness outcomes. We identified highly vulnerable genes in various processes, including novel targets unexplored for drug discovery. Equally important, we identified invulnerable essential genes, potentially explaining failed drug discovery efforts. Comparison of vulnerability between the reference and a hypervirulent Mtb isolate revealed incomplete conservation of vulnerability and that differential vulnerability can predict differential antibacterial susceptibility. Our results quantitatively redefine essential bacterial processes and identify high-value targets for drug development.

Incomplete transcripts dominate the Mycobacterium tuberculosis transcriptome.

Mycobacterium tuberculosis (Mtb) is a bacterial pathogen that causes tuberculosis (TB), an infectious disease that is responsible for major health and economic costs worldwide1. Mtb encounters diverse environments during its life cycle and responds to these changes largely by reprogramming its transcriptional output2. However, the mechanisms of Mtb transcription and how they are regulated remain poorly understood. Here we use a sequencing method that simultaneously determines both termini of individual RNA molecules in bacterial cells3 to profile the Mtb transcriptome at high resolution. Unexpectedly, we find that most Mtb transcripts are incomplete, with their 5' ends aligned at transcription start sites and 3' ends located 200-500 nucleotides downstream. We show that these short RNAs are mainly associated with paused RNA polymerases (RNAPs) rather than being products of premature termination. We further show that the high propensity of Mtb RNAP to pause early in transcription relies on the binding of the σ-factor. Finally, we show that a translating ribosome promotes transcription elongation, revealing a potential role for transcription-translation coupling in controlling Mtb gene expression. In sum, our findings depict a mycobacterial transcriptome that prominently features incomplete transcripts resulting from RNAP pausing. We propose that the pausing phase constitutes an important transcriptional checkpoint in Mtb that allows the bacterium to adapt to environmental changes and could be exploited for TB therapeutics.

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.

Mutations in rv0678 Confer Low-Level Resistance to Benzothiazinone DprE1 Inhibitors in Mycobacterium tuberculosis.

Tuberculosis (TB) is the leading cause of death from any bacterial infection, causing 1.5 million deaths worldwide each year. Due to the emergence of drug-resistant strains of Mycobacterium tuberculosis (Mtb) there have been significant efforts aimed at developing novel drugs to treat TB. One promising drug target in Mtb is the arabinogalactan biosynthetic enzyme DprE1, and there have been over a dozen unique chemical scaffolds identified which inhibit the activity of this protein. Among the most promising lead compounds are the benzothiazinones BTZ043 and PBTZ169, both of which are currently in or have completed phase IIa clinical trials. Due to the potential clinical utility of these drugs, we sought to identify potential synergistic interactions and new mechanisms of resistance using a genome-scale CRISPRi chemical-genetic screen with PBTZ169. We found that knockdown of rv0678, the negative regulator of the mmpS5/L5 drug efflux pump, confers resistance to PBTZ169. Mutations in rv0678 are the most common form of resistance to bedaquiline and there is already abundant evidence of these mutations emerging in bedaquiline-treated patients. We confirmed that rv0678 mutations from clinical isolates confer low level cross-resistance to BTZ043 and PBTZ169. While it is yet unclear whether rv0678 mutations would render benzothiazinones ineffective in treating TB, these results highlight the importance of monitoring for clinically prevalent rv0678 mutations during ongoing BTZ043 and PBTZ169 clinical trials.

Correction: TnSeq of Mycobacterium tuberculosis clinical isolates reveals strain-specific antibiotic liabilities.

[This corrects the article DOI: 10.1371/journal.ppat.1006939.].

Small RNA profiling in Mycobacterium tuberculosis identifies MrsI as necessary for an anticipatory iron sparing response.

One key to the success of Mycobacterium tuberculosis as a pathogen is its ability to reside in the hostile environment of the human macrophage. Bacteria adapt to stress through a variety of mechanisms, including the use of small regulatory RNAs (sRNAs), which posttranscriptionally regulate bacterial gene expression. However, very little is currently known about mycobacterial sRNA-mediated riboregulation. To date, mycobacterial sRNA discovery has been performed primarily in log-phase growth, and no direct interaction between any mycobacterial sRNA and its targets has been validated. Here, we performed large-scale sRNA discovery and expression profiling in M. tuberculosis during exposure to five pathogenically relevant stresses. From these data, we identified a subset of sRNAs that are highly induced in multiple stress conditions. We focused on one of these sRNAs, ncRv11846, here renamed mycobacterial regulatory sRNA in iron (MrsI). We characterized the regulon of MrsI and showed in mycobacteria that it regulates one of its targets, bfrA, through a direct binding interaction. MrsI mediates an iron-sparing response that is required for optimal survival of M. tuberculosis under iron-limiting conditions. However, MrsI is induced by multiple host-like stressors, which appear to trigger MrsI as part of an anticipatory response to impending iron deprivation in the macrophage environment.

TnSeq of Mycobacterium tuberculosis clinical isolates reveals strain-specific antibiotic liabilities.

Once considered a phenotypically monomorphic bacterium, there is a growing body of work demonstrating heterogeneity among Mycobacterium tuberculosis (Mtb) strains in clinically relevant characteristics, including virulence and response to antibiotics. However, the genetic and molecular basis for most phenotypic differences among Mtb strains remains unknown. To investigate the basis of strain variation in Mtb, we performed genome-wide transposon mutagenesis coupled with next-generation sequencing (TnSeq) for a panel of Mtb clinical isolates and the reference strain H37Rv to compare genetic requirements for in vitro growth across these strains. We developed an analytic approach to identify quantitative differences in genetic requirements between these genetically diverse strains, which vary in genomic structure and gene content. Using this methodology, we found differences between strains in their requirements for genes involved in fundamental cellular processes, including redox homeostasis and central carbon metabolism. Among the genes with differential requirements were katG, which encodes the activator of the first-line antitubercular agent isoniazid, and glcB, which encodes malate synthase, the target of a novel small-molecule inhibitor. Differences among strains in their requirement for katG and glcB predicted differences in their response to these antimicrobial agents. Importantly, these strain-specific differences in antibiotic response could not be predicted by genetic variants identified through whole genome sequencing or by gene expression analysis. Our results provide novel insight into the basis of variation among Mtb strains and demonstrate that TnSeq is a scalable method to predict clinically important phenotypic differences among Mtb strains.

Cdc15 integrates Tem1 GTPase-mediated spatial signals with Polo kinase-mediated temporal cues to activate mitotic exit.

In budding yeast, a Ras-like GTPase signaling cascade known as the mitotic exit network (MEN) promotes exit from mitosis. To ensure the accurate execution of mitosis, MEN activity is coordinated with other cellular events and restricted to anaphase. The MEN GTPase Tem1 has been assumed to be the central switch in MEN regulation. We show here that during an unperturbed cell cycle, restricting MEN activity to anaphase can occur in a Tem1 GTPase-independent manner. We found that the anaphase-specific activation of the MEN in the absence of Tem1 is controlled by the Polo kinase Cdc5. We further show that both Tem1 and Cdc5 are required to recruit the MEN kinase Cdc15 to spindle pole bodies, which is both necessary and sufficient to induce MEN signaling. Thus, Cdc15 functions as a coincidence detector of two essential cell cycle oscillators: the Polo kinase Cdc5 synthesis/degradation cycle and the Tem1 G-protein cycle. The Cdc15-dependent integration of these temporal (Cdc5 and Tem1 activity) and spatial (Tem1 activity) signals ensures that exit from mitosis occurs only after proper genome partitioning.

DNA Replication Fidelity in the Mycobacterium tuberculosis Complex.

Mycobacterium tuberculosis is genetically isolated, with no evidence for horizontal gene transfer or the acquisition of episomal genetic information in the modern evolution of strains of the Mycobacterium tuberculosis complex. When considered in the context of the specific features of the disease M. tuberculosis causes (e.g., transmission via cough aerosol, replication within professional phagocytes, subclinical persistence, and stimulation of a destructive immune pathology), this implies that to understand the mechanisms ensuring preservation of genomic integrity in infecting mycobacterial populations is to understand the source of genetic variation, including the emergence of microdiverse sub-populations that may be linked to the acquisition of drug resistance. In this chapter, we focus on mechanisms involved in maintaining DNA replication fidelity in M. tuberculosis, and consider the potential to target components of the DNA replication machinery as part of novel therapeutic regimens designed to curb the emerging threat of drug-resistance.

Precise genome modification in the crop species Zea mays using zinc-finger nucleases.

Agricultural biotechnology is limited by the inefficiencies of conventional random mutagenesis and transgenesis. Because targeted genome modification in plants has been intractable, plant trait engineering remains a laborious, time-consuming and unpredictable undertaking. Here we report a broadly applicable, versatile solution to this problem: the use of designed zinc-finger nucleases (ZFNs) that induce a double-stranded break at their target locus. We describe the use of ZFNs to modify endogenous loci in plants of the crop species Zea mays. We show that simultaneous expression of ZFNs and delivery of a simple heterologous donor molecule leads to precise targeted addition of an herbicide-tolerance gene at the intended locus in a significant number of isolated events. ZFN-modified maize plants faithfully transmit these genetic changes to the next generation. Insertional disruption of one target locus, IPK1, results in both herbicide tolerance and the expected alteration of the inositol phosphate profile in developing seeds. ZFNs can be used in any plant species amenable to DNA delivery; our results therefore establish a new strategy for plant genetic manipulation in basic science and agricultural applications.

Weak links: Advancing target-based drug discovery by identifying the most vulnerable targets.

Mycobacterium tuberculosis remains the most common infectious killer worldwide despite decades of antitubercular drug development. Effectively controlling the tuberculosis (TB) pandemic will require innovation in drug discovery. In this review, we provide a brief overview of the two main approaches to discovering new TB drugs-phenotypic screens and target-based drug discovery-and outline some of the limitations of each method. We then explore recent advances in genetic tools that aim to overcome some of these limitations. In particular, we highlight a novel metric to prioritize essential targets, termed vulnerability. Stratifying targets based on their vulnerability presents new opportunities for future target-based drug discovery campaigns.

Beyond antibiotic resistance: The whiB7 transcription factor coordinates an adaptive response to alanine starvation in mycobacteria.

Pathogenic mycobacteria are a significant cause of morbidity and mortality worldwide. The conserved whiB7 stress response reduces the effectiveness of antibiotic therapy by activating several intrinsic antibiotic resistance mechanisms. Despite our comprehensive biochemical understanding of WhiB7, the complex set of signals that induce whiB7 expression remain less clear. We employed a reporter-based, genome-wide CRISPRi epistasis screen to identify a diverse set of 150 mycobacterial genes whose inhibition results in constitutive whiB7 expression. We show that whiB7 expression is determined by the amino acid composition of the 5' regulatory uORF, thereby allowing whiB7 to sense amino acid starvation. Although deprivation of many amino acids can induce whiB7, whiB7 specifically coordinates an adaptive response to alanine starvation by engaging in a feedback loop with the alanine biosynthetic enzyme, aspC. These findings describe a metabolic function for whiB7 and help explain its evolutionary conservation across mycobacterial species occupying diverse ecological niches.

Functional genomics, proteomics, and regulatory DNA analysis in isogenic settings using zinc finger nuclease-driven transgenesis into a safe harbor locus in the human genome.

Isogenic settings are routine in model organisms, yet remain elusive for genetic experiments on human cells. We describe the use of designed zinc finger nucleases (ZFNs) for efficient transgenesis without drug selection into the PPP1R12C gene, a "safe harbor" locus known as AAVS1. ZFNs enable targeted transgenesis at a frequency of up to 15% following transient transfection of both transformed and primary human cells, including fibroblasts and hES cells. When added to this locus, transgenes such as expression cassettes for shRNAs, small-molecule-responsive cDNA expression cassettes, and reporter constructs, exhibit consistent expression and sustained function over 50 cell generations. By avoiding random integration and drug selection, this method allows bona fide isogenic settings for high-throughput functional genomics, proteomics, and regulatory DNA analysis in essentially any transformed human cell type and in primary cells.

Incomplete transcripts dominate the Mycobacterium tuberculosis transcriptome.

Mycobacterium tuberculosis (Mtb) is a bacterial pathogen that causes tuberculosis, an infectious disease that inflicts major health and economic costs around the world 1 . Mtb encounters a diversity of environments during its lifecycle, and responds to these changing environments by reprogramming its transcriptional output 2 . However, the transcriptomic features of Mtb remain poorly characterized. In this work, we comprehensively profile the Mtb transcriptome using the SEnd-seq method that simultaneously captures the 5' and 3' ends of RNA 3 . Surprisingly, we find that the RNA coverage for most of the Mtb transcription units display a gradual drop-off within a 200-500 nucleotide window downstream of the transcription start site, yielding a massive number of incomplete transcripts with heterogeneous 3' ends. We further show that the accumulation of these short RNAs is mainly due to the intrinsically low processivity of the Mtb transcription machinery rather than trans-acting factors such as Rho. Finally, we demonstrate that transcription-translation coupling plays a critical role in generating full-length protein-coding transcripts in Mtb. In sum, our results depict a mycobacterial transcriptome that is dominated by incomplete RNA products, suggesting a distinctive set of transcriptional regulatory mechanisms that could be exploited for new therapeutics.

Beyond antibiotic resistance: the whiB7 transcription factor coordinates an adaptive response to alanine starvation in mycobacteria.

Pathogenic mycobacteria are a significant cause of morbidity and mortality worldwide. These bacteria are highly intrinsically drug resistant, making infections challenging to treat. The conserved whiB7 stress response is a key contributor to mycobacterial intrinsic drug resistance. Although we have a comprehensive structural and biochemical understanding of WhiB7, the complex set of signals that activate whiB7 expression remain less clear. It is believed that whiB7 expression is triggered by translational stalling in an upstream open reading frame (uORF) within the whiB7 5' leader, leading to antitermination and transcription into the downstream whiB7 ORF. To define the signals that activate whiB7, we employed a genome-wide CRISPRi epistasis screen and identified a diverse set of 150 mycobacterial genes whose inhibition results in constitutive whiB7 activation. Many of these genes encode amino acid biosynthetic enzymes, tRNAs, and tRNA synthetases, consistent with the proposed mechanism for whiB7 activation by translational stalling in the uORF. We show that the ability of the whiB7 5' regulatory region to sense amino acid starvation is determined by the coding sequence of the uORF. The uORF shows considerable sequence variation among different mycobacterial species, but it is universally and specifically enriched for alanine. Providing a potential rationalization for this enrichment, we find that while deprivation of many amino acids can activate whiB7 expression, whiB7 specifically coordinates an adaptive response to alanine starvation by engaging in a feedback loop with the alanine biosynthetic enzyme, aspC. Our results provide a holistic understanding of the biological pathways that influence whiB7 activation and reveal an extended role for the whiB7 pathway in mycobacterial physiology, beyond its canonical function in antibiotic resistance. These results have important implications for the design of combination drug treatments to avoid whiB7 activation, as well as help explain the conservation of this stress response across a wide range of pathogenic and environmental mycobacteria.

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 essential M. tuberculosis Clp protease is functionally asymmetric in vivo.

The Clp protease system is a promising, noncanonical drug target against Mycobacterium tuberculosis (Mtb). Unlike in Escherichia coli, the Mtb Clp protease consists of two distinct proteolytic subunits, ClpP1 and ClpP2, which hydrolyze substrates delivered by the chaperones ClpX and ClpC1. While biochemical approaches uncovered unique aspects of Mtb Clp enzymology, its essentiality complicates in vivo studies. To address this gap, we leveraged new genetic tools to mechanistically interrogate the in vivo essentiality of the Mtb Clp protease. While validating some aspects of the biochemical model, we unexpectedly found that only the proteolytic activity of ClpP1, but not of ClpP2, is essential for substrate degradation and Mtb growth and infection. Our observations not only support a revised model of Mtb Clp biology, where ClpP2 scaffolds chaperone binding while ClpP1 provides the essential proteolytic activity of the complex; they also have important implications for the ongoing development of inhibitors toward this emerging therapeutic target.

CRISPR Interference Reveals That All-Trans-Retinoic Acid Promotes Macrophage Control of Mycobacterium tuberculosis by Limiting Bacterial Access to Cholesterol and Propionyl Coenzyme A.

Macrophages are a protective replicative niche for Mycobacterium tuberculosis (Mtb) but can kill the infecting bacterium when appropriately activated. To identify mechanisms of clearance, we compared levels of bacterial restriction by human macrophages after treatment with 26 compounds, including some currently in clinical trials for tuberculosis. All-trans-retinoic acid (ATRA), an active metabolite of vitamin A, drove the greatest increase in Mtb control. Bacterial clearance was transcriptionally and functionally associated with changes in macrophage cholesterol trafficking and lipid metabolism. To determine how these macrophage changes affected bacterial control, we performed the first Mtb CRISPR interference screen in an infection model, identifying Mtb genes specifically required to survive in ATRA-activated macrophages. These data showed that ATRA treatment starves Mtb of cholesterol and the downstream metabolite propionyl coenzyme A (propionyl-CoA). Supplementation with sources of propionyl-CoA, including cholesterol, abrogated the restrictive effect of ATRA. This work demonstrates that targeting the coupled metabolism of Mtb and the macrophage improves control of infection and that it is possible to genetically map the mode of bacterial death using CRISPR interference. IMPORTANCE Tuberculosis, caused by the bacterium Mycobacterium tuberculosis, is a leading cause of death due to infectious disease. Improving the immune response to tuberculosis holds promise for fighting the disease but is limited by our lack of knowledge as to how the immune system kills M. tuberculosis. Our research identifies a potent way to make relevant immune cells more effective at fighting M. tuberculosis and then uses paired human and bacterial genomic methods to determine the mechanism of that improved bacterial clearance.

Iron limitation in M. tuberculosis has broad impact on central carbon metabolism.

Mycobacterium tuberculosis (Mtb), the cause of the human pulmonary disease tuberculosis (TB), contributes to approximately 1.5 million deaths every year. Prior work has established that lipids are actively catabolized by Mtb in vivo and fulfill major roles in Mtb physiology and pathogenesis. We conducted a high-throughput screen to identify inhibitors of Mtb survival in its host macrophage. One of the hit compounds identified in this screen, sAEL057, demonstrates highest activity on Mtb growth in conditions where cholesterol was the primary carbon source. Transcriptional and functional data indicate that sAEL057 limits Mtb's access to iron by acting as an iron chelator. Furthermore, pharmacological and genetic inhibition of iron acquisition results in dysregulation of cholesterol catabolism, revealing a previously unappreciated linkage between these pathways. Characterization of sAEL057's mode of action argues that Mtb's metabolic regulation reveals vulnerabilities in those pathways that impact central carbon metabolism.

CRISPRi chemical genetics and comparative genomics identify genes mediating drug potency in Mycobacterium tuberculosis.

Mycobacterium tuberculosis (Mtb) infection is notoriously difficult to treat. Treatment efficacy is limited by Mtb's intrinsic drug resistance, as well as its ability to evolve acquired resistance to all antituberculars in clinical use. A deeper understanding of the bacterial pathways that influence drug efficacy could facilitate the development of more effective therapies, identify new mechanisms of acquired resistance, and reveal overlooked therapeutic opportunities. Here we developed a CRISPR interference chemical-genetics platform to titrate the expression of Mtb genes and quantify bacterial fitness in the presence of different drugs. We discovered diverse mechanisms of intrinsic drug resistance, unveiling hundreds of potential targets for synergistic drug combinations. Combining chemical genetics with comparative genomics of Mtb clinical isolates, we further identified several previously unknown mechanisms of acquired drug resistance, one of which is associated with a multidrug-resistant tuberculosis outbreak in South America. Lastly, we found that the intrinsic resistance factor whiB7 was inactivated in an entire Mtb sublineage endemic to Southeast Asia, presenting an opportunity to potentially repurpose the macrolide antibiotic clarithromycin to treat tuberculosis. This chemical-genetic map provides a rich resource to understand drug efficacy in Mtb and guide future tuberculosis drug development and treatment.