Global assessment of genomic regions required for growth in Mycobacterium tuberculosis

TnSeq identifies genetic requirements of Mycobacterium tuberculosis for survival under vaccine-induced immunity

Mycobacterium tuberculosis (Mtb), the etiologic agent of tuberculosis (TB), remains a persistent global health challenge due to the lack of an effective vaccine. The only licensed TB vaccine, Bacille Calmette-Guerin (BCG), is a live attenuated strain of Mycobacterium bovis that protects young children from severe disease but fails to provide protection through adulthood. It is unclear why BCG provides incomplete protection despite inducing a robust Th1 immune response. We set out to interrogate mycobacterial determinants of vaccine escape using a functional genomics approach, TnSeq, to define bacterial genes required for survival in mice vaccinated with BCG, the live attenuated Mtb vaccine strain, ΔLprG, and in mice with Mtb immunity conferred by prior infection. We find that critical virulence genes associated with acute infection and exponential growth are less essential in hosts with adaptive immunity, including genes encoding the Esx-1 and Mce1 systems. Genetic requirements for Mtb growth in vaccinated and previously Mtb-infected hosts mirror the genetic requirements reported for bacteria under in vitro conditions that reflect aspects of the adaptive immune response. Across distinct immunization conditions, differences in genetic requirements between live attenuated vaccines and vaccination routes are observed, suggesting that different immunization strategies impose distinct bacterial stressors. Collectively, these data support the idea that Mtb requires genes that enable stress adaptation and growth arrest upon encountering the restrictive host environment induced by the adaptive immune response. We demonstrate that TnSeq can be used to understand the bacterial genetic requirements for survival in vaccinated hosts across pre-clinical live attenuated vaccines and therefore may be applied to other vaccine modalities. Understanding how Mtb survives vaccine-induced immunity has the potential to inform the development of new vaccines or adjuvant therapies.

Interred mechanisms of resistance and host immune evasion revealed through network-connectivity analysis of M. tuberculosis complex graph pangenome

Mycobacterium tuberculosis complex successfully adapts to environmental pressures through mechanisms of rapid adaptation which remain poorly understood despite knowledge gained through decades of research. In this study, we used 110 reference-quality, complete de novo assembled, long-read sequenced clinical genomes to study patterns of structural adaptation through a graph-based pangenome analysis, elucidating rarely studied mechanisms that enable enhanced clinical phenotypes offering a novel perspective to the species' adaptation. Across isolates, we identified a pangenome of 4,325 genes (3,767 core and 558 accessory), revealing 290 novel genes, and a substantially more complete account of difficult-to-sequence esx/pe/pgrs/ppe genes. Seventy-four percent of core genes were deemed non-essential in vitro, 38% of which support the pathogen's survival in vivo, suggesting a need to broaden current perspectives on essentiality. Through information-theoretic analysis, we reveal the ppe genes that contribute most to the species' diversity-several with known consequences for antigenic variation and immune evasion. Construction of a graph pangenome revealed topological variations that implicate genes known to modulate host immunity (Rv0071-73, Rv2817c, cas2), defense against phages/viruses (cas2, csm6, and Rv2817c-2821c), and others associated with host tissue colonization. Here, the prominent trehalose transport pathway stands out for its involvement in caseous granuloma catabolism and the development of post-primary disease. We show paralogous duplications of genes implicated in bedaquiline (mmpL5 in all L1 isolates) and ethambutol (embC-A) resistance, with a paralogous duplication of its regulator (embR) in 96 isolates. We provide hypotheses for novel mechanisms of immune evasion and antibiotic resistance through gene dosing that can escape detection by molecular diagnostics.IMPORTANCEM. tuberculosis complex (MTBC) has killed over a billion people in the past 200 years alone and continues to kill nearly 1.5 million annually. The pathogen has a versatile ability to diversify under immune and drug pressure and survive, even becoming antibiotic persistent or resistant in the face of harsh chemotherapy. For proper diagnosis and design of an appropriate treatment regimen, a full understanding of this diversification and its clinical consequences is desperately needed. A mechanism of diversification that is rarely studied systematically is MTBC's ability to structurally change its genome. In this article, we have de novo assembled 110 clinical genomes (the largest de novo assembled set to date) and performed a pangenomic analysis. Our pangenome provides structural variation-based hypotheses for novel mechanisms of immune evasion and antibiotic resistance through gene dosing that can compromise molecular diagnostics and lead to further emergence of antibiotic resistance.

CRISPR/Cas12a-mediated gene silencing across diverse functional genes demonstrates single gene-specific spacer efficacy in Mycobacterium smegmatis

BACKGROUND

Tuberculosis, a persistent global health threat, necessitates a comprehensive understanding of the genes and pathways crucial for the survival and virulence of the causative pathogen, Mycobacterium tuberculosis. Working with M. tuberculosis (M.tb) presents significant challenges; therefore, the use of M. smegmatis as a surrogate system for conducting genetic studies of M.tb has proven to be highly valuable. Development of novel genetic tools to probe cellular processes accelerates the progress in the field of drug development and also helps in understanding the basic physiology of the bacterium.

RESULTS

This study reports the successful implementation and evaluation of the CRISPR-Cas12a system for gene repression in Mycobacterium smegmatis, a surrogate for M. tuberculosis. We engineered a Cas12a-based CRISPR interference (CRISPRi) system and assessed its functionality. Targeting 45 genes with a single sgRNA per gene, we achieved efficient gene repression, leading to marked phenotypic changes. Each knockdown strain was evaluated individually for growth phenotypes, and a comparison of the results with the reported essential gene library probed with dCas9 demonstrated congruous results across diverse gene categories. The study shows that CRISPR/Cas12a system can be effectively utilised with a single gene specific target for efficient silencing of the gene and highlights the importance of subsequent growth assays required to evaluate the vulnerability of targeted gene silencing.

CONCLUSION

Our findings reveal the robustness and versatility of the dCas12a-CRISPRi system in M. smegmatis, providing a valuable tool for functional genomics research. This work showcases the potential of the dCas12a-CRISPRi system in investigating essential genes, enabling a deeper understanding of the biology and potential therapeutic targets in mycobacterium species.

Widespread loss-of-function mutations implicating preexisting resistance to new or repurposed anti-tuberculosis drugs

BACKGROUND

Five New or Repurposed Drugs (NRDs) were approved in the last decade for treatment of multi-drug resistant tuberculosis: bedaquiline, clofazimine, linezolid, delamanid, and pretomanid. Unfortunately, resistance to these drugs emerged faster than anticipated, potentially due to preexisting resistance in naïve strains. Previous investigations into the rapid emergence have mostly included short variants. For the first time, we utilize de novo-assembled genomes, and systematically include Structural Variations (SV) and heterogeneity to comprehensively study this rapid emergence. We show high prevalence of preexisting resistance, identify novel markers of resistance, and lay the foundation for preventing preexisting resistance in future drug development.

METHODS

First, a systematic literature review revealed 313 NRD resistance variants in 13 genes. Next, 409 globally diverse clinical isolates collected prior to the drugs' programmatic use (308 were multidrug resistant, 106 had de novo assembled genomes) were utilized to study the 13 genes comprehensively for conventional, structural, and heterogeneous variants.

FINDINGS

We identified 5 previously reported and 67 novel putative NRD resistance variants. These variants were 2 promoter mutations (in 8/409 isolates), 13 frameshifts (21/409), 6 SVs (9/409), 35 heterogeneous frameshifts (32/409) and 11 heterogeneous SVs (12/106). Delamanid and pretomanid resistance mutations were most prevalent (48/409), while linezolid resistance mutations were least prevalent (8/409).

INTERPRETATION

Preexisting mutations implicated in resistance to at least one NRD was highly prevalent (85/409, 21 %). This was mostly caused by loss-of-function mutations in genes responsible for prodrug activation and efflux pump regulation. These preexisting mutations may have emerged through a bet-hedging strategy, or through cross-resistance with non-tuberculosis drugs such as metronidazole. Future drugs that could be resisted through loss-of-function in non-essential genes may suffer from preexisting resistance. The methods used here for comprehensive preexisting resistance assessment (especially SVs and heterogeneity) may mitigate this risk during early-stage drug development.

Factors governing attachment of Rhizobium leguminosarum to legume roots at acid, neutral, and alkaline pHs

Rhizobial attachment to host legume roots is the first physical interaction of bacteria and plants in symbiotic nitrogen fixation. The pH-dependent primary attachment of Rhizobium leguminosarum biovar viciae 3841 to Pisum sativum (pea) roots was investigated by genome-wide insertion sequencing, luminescence-based attachment assays, and proteomic analysis. Under acid, neutral, or alkaline pH, a total of 115 genes are needed for primary attachment under one or more environmental pH, with 22 genes required for all. These include components of cell surfaces and membranes, together with enzymes that construct and modify them. Mechanisms of dealing with stress also play a part; however, exact requirements vary depending on environmental pH. RNASeq showed that knocking out the two transcriptional regulators required for attachment causes massive changes in the bacterial cell surface. Approximately half of the 54 proteins required for attachment at pH 7.0 have a role in the later stages of nodule formation. We found no evidence for a single rhicadhesin responsible for alkaline attachment, although sonicated cell surface fractions inhibited root attachment at alkaline pH. Our results demonstrate the complexity of primary root attachment and illustrate the diversity of mechanisms involved.

IMPORTANCE

The first step by which bacteria interact with plant roots is by attachment. In this study, we use a combination of insertion sequencing and biochemical analysis to determine how bacteria attach to pea roots and how this is influenced by pH. We identify several key adhesins, which are molecules that enable bacteria to stick to roots. This includes a novel filamentous hemagglutinin which is needed at all pHs for attachment. Overall, 115 proteins are required for attachment at one or more pHs.

Cyclic AMP binding to a universal stress protein in Mycobacterium tuberculosis is essential for viability

Mycobacterial genomes encode multiple adenylyl cyclases and cAMP effector proteins, underscoring the diverse ways these bacteria utilize cAMP. We identified universal stress proteins, Rv1636 and MSMEG_3811 in Mycobacterium tuberculosis and Mycobacterium smegmatis, respectively, as abundantly expressed, novel cAMP-binding proteins. Rv1636 is secreted via the SecA2 secretion system in M. tuberculosis but is not directly responsible for the efflux of cAMP from the cell. In slow-growing mycobacteria, intrabacterial concentrations of Rv1636 were equivalent to the concentrations of cAMP present in the cell. In contrast, levels of intrabacterial MSMEG_3811 in M. smegmatis were lower than that of cAMP and therefore, overexpression of Rv1636 increased levels of "bound" cAMP. While msmeg_3811 could be readily deleted from the genome of M. smegmatis, we found that the rv1636 gene is essential for the viability of M. tuberculosis and is dependent on the cAMP-binding ability of Rv1636. Therefore, Rv1636 may function to regulate cAMP signaling by direct sequestration of the second messenger. This is the first evidence of a "sponge" for any second messenger in bacterial signaling that would allow mycobacterial cells to regulate the available intrabacterial "free" pool of cAMP.

Temporal genome-wide fitness analysis of Mycobacterium marinum during infection reveals the genetic requirement for virulence and survival in amoebae and microglial cells

Tuberculosis remains the most pervasive infectious disease and the recent emergence of drug-resistant strains emphasizes the need for more efficient drug treatments. A key feature of pathogenesis, conserved between the human pathogen Mycobacterium tuberculosis and the model pathogen Mycobacterium marinum, is the metabolic switch to lipid catabolism and altered expression of virulence genes at different stages of infection. This study aims to identify genes involved in sustaining viable intracellular infection. We applied transposon sequencing (Tn-Seq) to M. marinum, an unbiased genome-wide strategy combining saturation insertional mutagenesis and high-throughput sequencing. This approach allowed us to identify the localization and relative abundance of insertions in pools of transposon mutants. Gene essentiality and fitness cost of mutations were quantitatively compared between in vitro growth and different stages of infection in two evolutionary distinct phagocytes, the amoeba Dictyostelium discoideum and the murine BV2 microglial cells. In the M. marinum genome, 57% of TA sites were disrupted and 568 genes (10.2%) were essential, which is comparable to previous Tn-Seq studies on M. tuberculosis and M. bovis. Major pathways involved in the survival of M. marinum during infection of D. discoideum are related to DNA damage repair, lipid and vitamin metabolism, the type VII secretion system (T7SS) ESX-1, and the Mce1 lipid transport system. These pathways, except Mce1 and some glycolytic enzymes, were similarly affected in BV2 cells. These differences suggest subtly distinct nutrient availability or requirement in different host cells despite the known predominant use of lipids in both amoeba and microglial cells.IMPORTANCEThe emergence of biochemically and genetically tractable host model organisms for infection studies holds the promise to accelerate the pace of discoveries related to the evolution of innate immunity and the dissection of conserved mechanisms of cell-autonomous defenses. Here, we have used the genetically and biochemically tractable infection model system Dictyostelium discoideum/Mycobacterium marinum to apply a genome-wide transposon-sequencing experimental strategy to reveal comprehensively which mutations confer a fitness advantage or disadvantage during infection and compare these to a similar experiment performed using the murine microglial BV2 cells as host for M. marinum to identify conservation of virulence pathways between hosts.

Proline Dehydrogenase and Pyrroline 5 Carboxylate Dehydrogenase from Mycobacterium tuberculosis: Evidence for Substrate Channeling

In Mycobacterium tuberculosis, proline dehydrogenase (PruB) and ∆1-pyrroline-5-carboxylate (P5C) dehydrogenase (PruA) are monofunctional enzymes that catalyze proline oxidation to glutamate via the intermediates P5C and L-glutamate-γ-semialdehyde. Both enzymes are essential for the replication of pathogenic M. tuberculosis. Highly active enzymes were expressed and purified using a Mycobacterium smegmatis expression system. The purified enzymes were characterized using natural substrates and chemically synthesized analogs. The structural requirements of the quinone electron acceptor were examined. PruB displayed activity with all tested lipoquinone analogs (naphthoquinone or benzoquinone). In PruB assays utilizing analogs of the native naphthoquinone [MK-9 (II-H2)] specificity constants Kcat/Km were an order of magnitude greater for the menaquinone analogs than the benzoquinone analogs. In addition, mycobacterial PruA was enzymatically characterized for the first time using exogenous chemically synthesized P5C. A Km value of 120 ± 0.015 µM was determined for P5C, while the Km value for NAD+ was shown to be 33 ± 4.3 µM. Furthermore, proline competitively inhibited PruA activity and coupled enzyme assays, suggesting that the recombinant purified monofunctional PruB and PruA enzymes of M. tuberculosis channel substrate likely increase metabolic flux and protect the bacterium from methylglyoxal toxicity.

High-throughput functional genomics: A (myco)bacterial perspective

Advances in sequencing technologies have enabled unprecedented insights into bacterial genome composition and dynamics. However, the disconnect between the rapid acquisition of genomic data and the (much slower) confirmation of inferred genetic function threatens to widen unless techniques for fast, high-throughput functional validation can be applied at scale. This applies equally to Mycobacterium tuberculosis, the leading infectious cause of death globally and a pathogen whose genome, despite being among the first to be sequenced two decades ago, still contains many genes of unknown function. Here, we summarize the evolution of bacterial high-throughput functional genomics, focusing primarily on transposon (Tn)-based mutagenesis and the construction of arrayed mutant libraries in diverse bacterial systems. We also consider the contributions of CRISPR interference as a transformative technique for probing bacterial gene function at scale. Throughout, we situate our analysis within the context of functional genomics of mycobacteria, focusing specifically on the potential to yield insights into M. tuberculosis pathogenicity and vulnerabilities for new drug and regimen development. Finally, we offer suggestions for future approaches that might be usefully applied in elucidating the complex cellular biology of this major human pathogen.

Mycobacterial Genetic Technologies for Probing the Host-Pathogen Microenvironment

Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis, is one of the oldest and most successful pathogens in the world. Diverse selective pressures encountered within host cells have directed the evolution of unique phenotypic traits, resulting in the remarkable evolutionary success of this largely obligate pathogen. Despite centuries of study, the genetic repertoire utilized by Mtb to drive virulence and host immune evasion remains to be fully understood. Various genetic approaches have been and continue to be developed to tackle the challenges of functional gene annotation and validation in an intractable organism such as Mtb. In vitro and ex vivo systems remain the primary approaches to generate and confirm hypotheses that drive a general understanding of mycobacteria biology. However, it remains of great importance to characterize genetic requirements for successful infection within a host system as in vitro and ex vivo studies fail to fully replicate the complex microenvironment experienced by Mtb. In this review, we evaluate the employment of the mycobacterial genetic toolkit to probe the host-pathogen interface by surveying the current state of mycobacterial genetic studies within host systems, with a major focus on the murine model. Specifically, we discuss the different ways that these tools have been utilized to examine various aspects of infection, including bacterial survival/virulence, bacterial evasion of host immunity, and development of novel antibacterial/vaccine strategies.

Identification of essential genes in Coxiella burnetii

Coxiella burnetii is an intracellular pathogen responsible for causing Q fever in humans, a disease with varied presentations ranging from a mild flu-like sickness to a debilitating illness that can result in endocarditis. The intracellular lifestyle of C. burnetii is unique, residing in an acidic phagolysosome-like compartment within host cells. An understanding of the core molecular biology of C. burnetii will greatly increase our understanding of C. burnetii growth, survival and pathogenesis. We used transposon-directed insertion site sequencing (TraDIS) to reveal C. burnetii Nine Mile Phase II genes fundamental for growth and in vitro survival. Screening a transposon library containing >10 000 unique transposon mutants revealed 512 predicted essential genes. Essential routes of synthesis were identified for the mevalonate pathway, as well as peptidoglycan and biotin synthesis. Some essential genes identified (e.g. predicted type IV secretion system effector genes) are typically considered to be associated with C. burnetii virulence, a caveat concerning the axenic media used in the study. Investigation into the conservation of the essential genes identified revealed that 78 % are conserved across all C. burnetii strains sequenced to date, which probably play critical functions. This is the first report of a whole genome transposon screen in C. burnetii that has been undertaken for the identification of essential genes.

Construction of high-density transposon mutant library of Staphylococcus aureus using bacteriophage ϕ11

Transposon mutant libraries are an important resource to study bacterial metabolism and pathogenesis. The fitness analysis of mutants in the libraries under various growth conditions provides important clues to study the physiology and biogenesis of structural components of a bacterial cell. A transposon library in conjunction with next-generation sequencing techniques, collectively named transposon sequencing (Tn-seq), enables high-throughput genome profiling and synthetic lethality analysis. Tn-seq has also been used to identify essential genes and to study the mode of action of antibacterials. To construct a high-density transposon mutant library, an efficient delivery system for transposition in a model bacterium is essential. Here, I describe a detailed protocol for generating a high-density phage-based transposon mutant library in a Staphylococcus aureus strain, and this protocol is readily applicable to other S. aureus strains including USA300 and MW2.

Genetic Engineering in Mycobacteria

Abstract

Genetic tools for targeted modification of the mycobacterial genome contribute to the understanding of the physiology and virulence mechanisms of mycobacteria. Human and animal pathogens, such as the Mycobacterium tuberculosis complex, which causes tuberculosis, and M. leprae, which causes leprosy, are of particular importance. Genetic research opens up novel opportunities to identify and validate new targets for antibacterial drugs and to develop improved vaccines. Although mycobacteria are difficult to work with due to their slow growth rate and a limited possibility to transfer genetic information, significant progress has been made in developing genetic engineering methods for mycobacteria. The review considers the main approaches to changing the mycobacterial genome in a targeted manner, including homologous and site-specific recombination and use of the CRISPR/Cas system.

Comparative genome analysis of mycobacteria focusing on tRNA and non-coding RNA

BACKGROUND

The Mycobacterium genus encompasses at least 192 named species, many of which cause severe diseases such as tuberculosis. Non-tuberculosis mycobacteria (NTM) can also infect humans and animals. Some are of emerging concern because they show high resistance to commonly used antibiotics while others are used and evaluated in bioremediation or included in anticancer vaccines.

RESULTS

We provide the genome sequences for 114 mycobacterial type strains and together with 130 available mycobacterial genomes we generated a phylogenetic tree based on 387 core genes and supported by average nucleotide identity (ANI) data. The 244 genome sequences cover most of the species constituting the Mycobacterium genus. The genome sizes ranged from 3.2 to 8.1 Mb with an average of 5.7 Mb, and we identified 14 new plasmids. Moreover, mycobacterial genomes consisted of phage-like sequences ranging between 0 and 4.64% dependent on mycobacteria while the number of IS elements varied between 1 and 290. Our data also revealed that, depending on the mycobacteria, the number of tRNA and non-coding (nc) RNA genes differ and that their positions on the chromosome varied. We identified a conserved core set of 12 ncRNAs, 43 tRNAs and 18 aminoacyl-tRNA synthetases among mycobacteria.

CONCLUSIONS

Phages, IS elements, tRNA and ncRNAs appear to have contributed to the evolution of the Mycobacterium genus where several tRNA and ncRNA genes have been horizontally transferred. On the basis of our phylogenetic analysis, we identified several isolates of unnamed species as new mycobacterial species or strains of known mycobacteria. The predicted number of coding sequences correlates with genome size while the number of tRNA, rRNA and ncRNA genes does not. Together these findings expand our insight into the evolution of the Mycobacterium genus and as such they establish a platform to understand mycobacterial pathogenicity, their evolution, antibiotic resistance/tolerance as well as the function and evolution of ncRNA among mycobacteria.

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

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

Mycobacterium tuberculosis KasA as a drug target: Structure-based inhibitor design

Recent studies have reported the β-ketoacyl-acyl carrier protein KasA as a druggable target for Mycobacterium tuberculosis. This review summarizes the current status of major classes of KasA inhibitors with an emphasis on significant contributions from structure-based design methods leveraging X-ray crystal structures of KasA alone and in complex with inhibitors. The issues addressed within each inhibitor class are discussed while detailing the characterized interactions with KasA and structure-activity relationships. A critical analysis of these findings should lay the foundation for new KasA inhibitors to study the basic biology of M. tuberculosis and to form the basis of new antitubercular molecules of clinical significance with activity against drug-sensitive and drug-resistant infections.

Uncovering interactions between mycobacterial respiratory complexes to target drug-resistant Mycobacterium tuberculosis

Mycobacterium tuberculosis remains a leading cause of infectious disease morbidity and mortality for which new drug combination therapies are needed. Mycobacterial bioenergetics has emerged as a promising space for the development of novel therapeutics. Further to this, unique combinations of respiratory inhibitors have been shown to have synergistic or synthetic lethal interactions, suggesting that combinations of bioenergetic inhibitors could drastically shorten treatment times. Realizing the full potential of this unique target space requires an understanding of which combinations of respiratory complexes, when inhibited, have the strongest interactions and potential in a clinical setting. In this review, we discuss (i) chemical-interaction, (ii) genetic-interaction and (iii) chemical-genetic interaction studies to explore the consequences of inhibiting multiple mycobacterial respiratory components. We provide potential mechanisms to describe the basis for the strongest interactions. Finally, whilst we place an emphasis on interactions that occur with existing bioenergetic inhibitors, by highlighting interactions that occur with alternative respiratory components we envision that this information will provide a rational to further explore alternative proteins as potential drug targets and as part of unique drug combinations.

Molecular Connectivity between Extracytoplasmic Sigma Factors and PhoP Accounts for Coupled Mycobacterial Stress Response

Mycobacterium tuberculosis encounters numerous stress conditions within the host, but how it is able to mount a coupled stress response remains unknown. Growing evidence suggests that under acidic pH, M. tuberculosis modulates redox homeostasis. In an attempt to dissect the mechanistic details of responses to multiple stress conditions, here we studied the significance of connectivity of extracytoplasmic sigma factors with PhoP. We show that PhoP impacts the mycothiol redox state, and the H37Rv ΔphoP deletion mutant strain displays a significantly higher susceptibility to redox stress than the wild-type bacilli. To probe how the two regulators PhoP and redox-active sigma factor SigH contribute to redox homeostasis, we show that SigH controls expression of redox-active thioredoxin genes, a major mycobacterial antioxidant system, and under redox stress, SigH, but not PhoP, is recruited at the target promoters. Consistent with these results, interaction between PhoP and SigH fails to impact redox-dependent gene expression. This is in striking contrast to our previous results showing PhoP-dependent SigE recruitment within acid-inducible mycobacterial promoters to maintain pH homeostasis. Our subsequent results demonstrate reduced PhoP-SigH interaction in the presence of diamide and enhanced PhoP-SigE interaction under low pH. These contrasting results uncover the underlying mechanism of the mycobacterial adaptive program, coupling low pH with maintenance of redox homeostasis. IMPORTANCE M. tuberculosis encounters reductive stress under acidic pH. To investigate the mechanism of coupled stress response, we show that PhoP plays a major role in mycobacterial redox stress response. We observed a strong correlation of phoP-dependent redox-active expression of thioredoxin genes, a major mycobacterial antioxidant system. Further probing of functioning of regulators revealed that while PhoP controls pH homeostasis via its interaction with SigE, direct recruitment of SigH, but not PhoP-SigH interaction, controls expression of thioredoxin genes. These strikingly contrasting results showing enhanced PhoP-SigE interaction under acidic pH and reduced PhoP-SigH interaction under redox conditions uncover the underlying novel mechanism of the mycobacterial adaptive program, coupling low pH with maintenance of redox homeostasis.

Gradients in gene essentiality reshape antibacterial research

Essential genes encode the processes that are necessary for life. Until recently, commonly applied binary classifications left no space between essential and non-essential genes. In this review, we frame bacterial gene essentiality in the context of genetic networks. We explore how the quantitative properties of gene essentiality are influenced by the nature of the encoded process, environmental conditions and genetic background, including a strain's distinct evolutionary history. The covered topics have important consequences for antibacterials, which inhibit essential processes. We argue that the quantitative properties of essentiality can thus be used to prioritize antibacterial cellular targets and desired spectrum of activity in specific infection settings. We summarize our points with a case study on the core essential genome of the cystic fibrosis pathobiome and highlight avenues for targeted antibacterial development.

Understanding the Genetic Diversity of Mycobacterium africanum Using Phylogenetics and Population Genomics Approaches

A total of two lineages of Mycobacterium tuberculosis var. africanum (Maf), L5 and L6, which are members of the Mycobacterium tuberculosis complex (MTBC), are responsible for causing tuberculosis in West Africa. Regions of difference (RDs) are usually used for delineation of MTBC. With increased data availability, single nucleotide polymorphisms (SNPs) promise to provide better resolution. Publicly available 380 Maf samples were analyzed for identification of “core-cluster-specific-SNPs,” while additional 270 samples were used for validation. RD-based methods were used for lineage-assignment, wherein 31 samples remained unidentified. The genetic diversity of Maf was estimated based on genome-wide SNPs using phylogeny and population genomics approaches. Lineage-based clustering (L5 and L6) was observed in the whole genome phylogeny with distinct sub-clusters. Population stratification using both model-based and de novo approaches supported the same observations. L6 was further delineated into three sub-lineages (L6.1–L6.3), whereas L5 was grouped as L5.1 and L5.2 based on the occurrence of RD711. L5.1 and L5.2 were further divided into two (L5.1.1 and L5.1.2) and four (L5.2.1–L5.2.4) sub-clusters, respectively. Unassigned samples could be assigned to definite lineages/sub-lineages based on clustering observed in phylogeny along with high-confidence posterior membership scores obtained during population stratification. Based on the (sub)-clusters delineated, “core-cluster-specific-SNPs” were derived. Synonymous SNPs (137 in L5 and 128 in L6) were identified as biomarkers and used for validation. Few of the cluster-specific missense variants in L5 and L6 belong to the central carbohydrate metabolism pathway which include His6Tyr (Rv0946c), Glu255Ala (Rv1131), Ala309Gly (Rv2454c), Val425Ala and Ser112Ala (Rv1127c), Gly198Ala (Rv3293) and Ile137Val (Rv0363c), Thr421Ala (Rv0896), Arg442His (Rv1248c), Thr218Ile (Rv1122), and Ser381Leu (Rv1449c), hinting at the differential growth attenuation. Genes harboring multiple (sub)-lineage-specific “core-cluster” SNPs such as Lys117Asn, Val447Met, and Ala455Val (Rv0066c; icd2) present across L6, L6.1, and L5, respectively, hinting at the association of these SNPs with selective advantage or host-adaptation. Cluster-specific SNPs serve as additional markers along with RD-regions for Maf delineation. The identified SNPs have the potential to provide insights into the genotype–phenotype correlation and clues for endemicity of Maf in the African population.

A Feedback Regulatory Loop Containing McdR and WhiB2 Controls Cell Division and DNA Repair in Mycobacteria

Cell division must be coordinated with DNA repair, which is strictly regulated in response to different drugs and environmental stresses in bacteria. However, the mechanisms by which mycobacteria orchestrate these two processes remain largely uncharacterized. Here, we report a regulatory loop between two essential mycobacterial regulators, McdR (Rv1830) and WhiB2, in coordinating the processes of cell division and DNA repair. McdR inhibits cell division-associated whiB2 expression by binding to the AATnACAnnnnTGTnATT motif in the promoter region. Furthermore, McdR overexpression simultaneously activates imuAB and dnaE2 expression to promote error-prone DNA repair, which facilitates genetic adaptation to stress conditions. Through a feedback mechanism, WhiB2 activates mcdR expression by binding to the cGACACGc motif in the promoter region. Importantly, analyses of mutations in clinical Mycobacterium tuberculosis strains indicate that disruption of this McdR-WhiB2 feedback regulatory loop influences expression of both cell growth- and DNA repair-associated genes, which further supports the contribution of McdR-WhiB2 regulatory loop in regulating mycobacterial cell growth and drug resistance. This highly conserved feedback regulatory loop provides fresh insight into the link between mycobacterial cell growth control and stress responses. IMPORTANCE Drug-resistant M. tuberculosis poses a threat to the control and prevention of tuberculosis (TB) worldwide. Thus, there is a need to identify the mechanisms enabling M. tuberculosis to adapt and grow under drug-induced stress. Rv1830 has been shown to be associated with drug resistance in M. tuberculosis, but its mechanisms have not yet been elucidated. Here, we reveal a regulatory role of Rv1830, which coordinates cell division and DNA repair in mycobacteria, and rename it McdR (mycobacterial cell division regulator). An increase in McdR levels represses the expression of cell division-associated whiB2 but activates the DNA repair-associated, error-prone enzymes ImuA/B and DnaE2, which in turn facilitates adaptation to stress responses and drug resistance. Furthermore, WhiB2 activates the transcription of mcdR to form a conserved regulatory loop. These data provide new insights into the mechanisms controlling mycobacterial cell growth and stress responses.

Identification of putative essential protein domains from high-density transposon insertion sequencing

A first clue to gene function can be obtained by examining whether a gene is required for life in certain standard conditions, that is, whether a gene is essential. In bacteria, essential genes are usually identified by high-density transposon mutagenesis followed by sequencing of insertion sites (Tn-seq). These studies assign the term "essential" to whole genes rather than the protein domain sequences that encode the essential functions. However, genes can code for multiple protein domains that evolve their functions independently. Therefore, when essential genes code for more than one protein domain, only one of them could be essential. In this study, we defined this subset of genes as "essential domain-containing" (EDC) genes. Using a Tn-seq data set built-in Burkholderia cenocepacia K56-2, we developed an in silico pipeline to identify EDC genes and the essential protein domains they encode. We found forty candidate EDC genes and demonstrated growth defect phenotypes using CRISPR interference (CRISPRi). This analysis included two knockdowns of genes encoding the protein domains of unknown function DUF2213 and DUF4148. These putative essential domains are conserved in more than two hundred bacterial species, including human and plant pathogens. Together, our study suggests that essentiality should be assigned to individual protein domains rather than genes, contributing to a first functional characterization of protein domains of unknown function.

An essential periplasmic protein coordinates lipid trafficking and is required for asymmetric polar growth in mycobacteria

Mycobacteria, including the human pathogen Mycobacterium tuberculosis, grow by inserting new cell wall material at their poles. This process and that of division are asymmetric, producing a phenotypically heterogeneous population of cells that respond non-uniformly to stress (Aldridge et al., 2012; Rego et al., 2017). Surprisingly, deletion of a single gene - lamA - leads to more symmetry, and to a population of cells that is more uniformly killed by antibiotics (Rego et al., 2017). How does LamA create asymmetry? Here, using a combination of quantitative time-lapse imaging, bacterial genetics, and lipid profiling, we find that LamA recruits essential proteins involved in cell wall synthesis to one side of the cell - the old pole. One of these proteins, MSMEG_0317, here renamed PgfA, was of unknown function. We show that PgfA is a periplasmic protein that interacts with MmpL3, an essential transporter that flips mycolic acids in the form of trehalose monomycolate (TMM), across the plasma membrane. PgfA interacts with a TMM analog suggesting a direct role in TMM transport. Yet our data point to a broader function as well, as cells with altered PgfA levels have differences in the abundance of other lipids and are differentially reliant on those lipids for survival. Overexpression of PgfA, but not MmpL3, restores growth at the old poles in cells missing lamA. Together, our results suggest that PgfA is a key determinant of polar growth and cell envelope composition in mycobacteria, and that the LamA-mediated recruitment of this protein to one side of the cell is a required step in the establishment of cellular asymmetry.

Probing Differences in Gene Essentiality Between the Human and Animal Adapted Lineages of the Mycobacterium tuberculosis Complex Using TnSeq

Members of the Mycobacterium tuberculosis complex (MTBC) show distinct host adaptations, preferences and phenotypes despite being >99% identical at the nucleic acid level. Previous studies have explored gene expression changes between the members, however few studies have probed differences in gene essentiality. To better understand the functional impacts of the nucleic acid differences between Mycobacterium bovis and Mycobacterium tuberculosis, we used the Mycomar T7 phagemid delivery system to generate whole genome transposon libraries in laboratory strains of both species and compared the essentiality status of genes during growth under identical in vitro conditions. Libraries contained insertions in 54% of possible TA sites in M. bovis and 40% of those present in M. tuberculosis, achieving similar saturation levels to those previously reported for the MTBC. The distributions of essentiality across the functional categories were similar in both species. 527 genes were found to be essential in M. bovis whereas 477 genes were essential in M. tuberculosis and 370 essential genes were common in both species. CRISPRi was successfully utilised in both species to determine the impacts of silencing genes including wag31, a gene involved in peptidoglycan synthesis and Rv2182c/Mb2204c, a gene involved in glycerophospholipid metabolism. We observed species specific differences in the response to gene silencing, with the inhibition of expression of Mb2204c in M. bovis showing significantly less growth impact than silencing its orthologue (Rv2182c) in M. tuberculosis. Given that glycerophospholipid metabolism is a validated pathway for antimicrobials, our observations suggest that target vulnerability in the animal adapted lineages cannot be assumed to be the same as the human counterpart. This is of relevance for zoonotic tuberculosis as it implies that the development of antimicrobials targeting the human adapted lineage might not necessarily be effective against the animal adapted lineage. The generation of a transposon library and the first reported utilisation of CRISPRi in M. bovis will enable the use of these tools to further probe the genetic basis of survival under disease relevant conditions.

Reprogramming Mycobacterium tuberculosis CRISPR System for Gene Editing and Genome-wide RNA Interference Screening

Mycobacterium tuberculosis is the causative agent of tuberculosis (TB), which is still the leading cause of mortality from a single infectious disease worldwide. The development of novel anti-TB drugs and vaccines is severely hampered by the complicated and time-consuming genetic manipulation techniques for M. tuberculosis. Here, we harnessed an endogenous type III-A CRISPR/Cas10 system of M. tuberculosis for efficient gene editing and RNA interference (RNAi). This simple and easy method only needs to transform a single mini-CRISPR array plasmid, thus avoiding the introduction of exogenous protein and minimizing proteotoxicity. We demonstrated that M. tuberculosis genes can be efficiently and specifically knocked in/out by this system as confirmed by DNA high-throughput sequencing. This system was further applied to single- and multiple-gene RNAi. Moreover, we successfully performed genome-wide RNAi screening to identify M. tuberculosis genes regulating in vitro and intracellular growth. This system can be extensively used for exploring the functional genomics of M. tuberculosis and facilitate the development of novel anti-TB drugs and vaccines.

A large-scale whole-genome sequencing analysis reveals false positives of bacterial essential genes

Essential genes are crucial for bacterial viability and represent attractive targets for novel anti-pathogen drug discovery. However, essential genes determined by the transposon insertion sequencing (Tn-seq) approach often contain many false positives. We hypothesized that some of those false positives are genes that are actually deleted from the genome, so they do not present any transposon insertion in the course of Tn-seq analysis. Based on this assumption, we performed a large-scale whole-genome sequencing analysis for the bacterium of interest. Our analysis revealed that some "essential genes" are indeed removed from the analyzed bacterial genomes. Since these genes were kicked out by bacteria, they should not be defined as essential. Our work showed that gene deletion is one of the false positive sources of essentiality determination, which is apparently underestimated in previous studies. We suggest subtracting the genome backgrounds before the evaluation of Tn-seq, and created a list of false positive gene essentiality as a reference for the downstream application. KEY POINTS: • Discovery of false positives of essential genes defined previously through the analyses of a large scale of whole-genome sequencing data • These false positives are the results of gene deletions in the studied genomes • Sequencing the target genome before Tn-seq analysis is of importance while some studies neglected it.

Rv0180c contributes to Mycobacterium tuberculosis cell shape and to infectivity in mice and macrophages

Mycobacterium tuberculosis, the main causative agent of human tuberculosis, is transmitted from person to person via small droplets containing very few bacteria. Optimizing the chance to seed in the lungs is therefore a major adaptation to favor survival and dissemination in the human population. Here we used TnSeq to identify genes important for the early events leading to bacterial seeding in the lungs. Beside several genes encoding known virulence factors, we found three new candidates not previously described: rv0180c, rv1779c and rv1592c. We focused on the gene, rv0180c, of unknown function. First, we found that deletion of rv0180c in M. tuberculosis substantially reduced the initiation of infection in the lungs of mice. Next, we established that Rv0180c enhances entry into macrophages through the use of complement-receptor 3 (CR3), a major phagocytic receptor for M. tuberculosis. Silencing CR3 or blocking the CR3 lectin site abolished the difference in entry between the wild-type parental strain and the Δrv0180c::km mutant. However, we detected no difference in the production of both CR3-known carbohydrate ligands (glucan, arabinomannan, mannan), CR3-modulating lipids (phthiocerol dimycocerosate), or proteins in the capsule of the Δrv0180c::km mutant in comparison to the wild-type or complemented strains. By contrast, we established that Rv0180c contributes to the functionality of the bacterial cell envelope regarding resistance to toxic molecule attack and cell shape. This alteration of bacterial shape could impair the engagement of membrane receptors that M. tuberculosis uses to invade host cells, and open a new perspective on the modulation of bacterial infectivity.

The epidemic efficiency of tuberculosis bacilli is determined by their capacity to transmit via aerosol. Currently, the bacterial functions that favor Mycobacterium tuberculosis seeding in the lung of naïve host remain mostly unknown. Here we implemented a genome-wide approach to identify M. tuberculosis mutants deficient for seeding and early replication in the lung of mice. In addition to genes known to encode virulence factors, we identified three genes not previously described. We used complementary approaches to characterize the phenotype of a M. tuberculosis mutant with insertion within the rv0180c gene. We found that this mutant is impaired for seeding in the lung of mice and for invasion and replication in human macrophages. In macrophages, the defect relies on a lack of engagement of CR3 receptor. Although we did not detect any difference between the wild type strain and the rv0180c mutant with regard to potential CR3-ligand, we found that the bacterial cell envelope is altered in the rv0180c mutant. Our study provides new insight into bacterial genes required for early interaction of M. tuberculosis with the host and perspective to understand the bacterial functions enhancing infectivity.

Total Synthesis of Tetrahydrolipstatin, Its Derivatives, and Evaluation of Their Ability to Potentiate Multiple Antibiotic Classes against Mycobacterium Species

Tetrahydrolipstatin (THL, 1a) has been shown to inhibit both mammalian and bacterial α/β hydrolases. In the case of bacterial systems, THL is a known inhibitor of several Mycobacterium tuberculosis hydrolases involved in mycomembrane biosynthesis. Herein we report a highly efficient eight-step asymmetric synthesis of THL using a route that allows modification of the THL α-chain substituent to afford compounds 1a through 1e. The key transformation in the synthesis was use of a (TPP)CrCl/Co2(CO)8-catalyzed regioselective and stereospecific carbonylation on an advanced epoxide intermediate to yield a trans-β-lactone. These compounds are modest inhibitors of Ag85A and Ag85C, two α/β hydrolases of M. tuberculosis involved in the biosynthesis of the mycomembrane. Among these compounds, 10d showed the highest inhibitory effect on Ag85A (34 ± 22 μM) and Ag85C (66 ± 8 μM), and its X-ray structure was solved in complex with Ag85C to 2.5 Å resolution. In contrast, compound 1e exhibited the best-in-class MICs of 50 μM (25 μg/mL) and 16 μM (8.4 μg/mL) against M. smegmatis and M. tuberculosis H37Ra, respectively, using a microtiter assay plate. Combination of 1e with 13 well-established antibiotics synergistically enhanced the potency of few of these antibiotics in M. smegmatis and M. tuberculosis H37Ra. Compound 1e applied at concentrations 4-fold lower than its MIC enhanced the MIC of the synergistic antibiotic by 2-256-fold. In addition to observing synergy with first-line drugs, rifamycin and isoniazid, the MIC of vancomycin against M. tuberculosis H37Ra was 65 μg/mL; however, the MIC was lowered to 0.25 μg/mL in the presence of 2.1 μg/mL 1e demonstrating the potential of targeting mycobacterial hydrolases involved in mycomembrane and peptidoglycan biosynthesis.

Mutant resources for functional genomics in Dictyostelium discoideum using REMI-seq technology

Background

Genomes can be sequenced with relative ease, but ascribing gene function remains a major challenge. Genetically tractable model systems are crucial to meet this challenge. One powerful model is the social amoeba Dictyostelium discoideum, a eukaryotic microbe widely used to study diverse questions in the cell, developmental and evolutionary biology.

Results

We describe REMI-seq, an adaptation of Tn-seq, which allows high throughput, en masse, and quantitative identification of the genomic site of insertion of a drug resistance marker after restriction enzyme-mediated integration. We use REMI-seq to develop tools which greatly enhance the efficiency with which the sequence, transcriptome or proteome variation can be linked to phenotype in D. discoideum. These comprise (1) a near genome-wide resource of individual mutants and (2) a defined pool of ‘barcoded’ mutants to allow large-scale parallel phenotypic analyses. These resources are freely available and easily accessible through the REMI-seq website that also provides comprehensive guidance and pipelines for data analysis. We demonstrate that integrating these resources allows novel regulators of cell migration, phagocytosis and macropinocytosis to be rapidly identified.

Conclusions

We present methods and resources, generated using REMI-seq, for high throughput gene function analysis in a key model system.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01108-y.

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.

Structure-based identification of natural compound inhibitor against M. tuberculosis thioredoxin reductase: insight from molecular docking and dynamics simulation

Antioxidant systems of M. tuberculosis (Mtb) play an important role in providing resistance in the hostile environment of mononuclear phagocytes. Thioredoxin system is a known antioxidant system that consists of three copies of thioredoxins (Trxs) and a single copy of thioredoxin reductase (TrxR). TrxR has been validated as an essential gene known to be involved in the reduction of peroxides, dinitrobenzenes and hydroperoxides, and is crucial in maintaining the survival of Mtb in macrophages. Recently, it has been demonstrated to be a druggable target. In this study, molecular docking was applied to screen more than 20,000 natural compounds from the Traditional Chinese Medicine database. Theoretical calculation of ΔGbinding by the Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) methods indicated two top-hit compounds that bind with a high affinity to the allosteric site, consisting of a hinge region, of TrxR. Further, stability and binding analysis of both compounds were carried out with molecular dynamics simulation. An analysis of conformational variation by principal component analysis (PCA) and protein contact network (PCN) uncovered the conformational changes in the compound-bound forms of protein. The NADPH domain formed many new interactions with the FAD domain in the compound-bound form, signifying that the binding may render an effect on the protein structure and function. Our results suggest that these two compounds could potentially be used for structure-based lead inhibitors against TrxR. The inhibitor selected as lead compound will be used further as a scaffold to optimize as novel anti-tuberculosis therapeutic.Communicated by Ramaswamy H. Sarma.

CRISPR Interference of Adenylate Cyclases from Mycobacterium tuberculosis

This work describes a modification of the pRH2521 vector of the pRH2502/pRH2521 system for CRISPR-dCas9-mediated RNA interference. The modification enabled an increase in the cloning efficiency of guide RNA spacers. The ability of the modified pRH2502/pRH2521 system to suppress the transcription of certain genes was evaluated with the use of genes of Mycobacterium tuberculosis adenylate cyclases. The results revealed the limitations of the pRH2502/pRH2521 system for CRISPR interference associated with the probability of the detection of a protospacer adjacent motif (PAM) in the gene promoter region.

Profiling Pretomanid as a Therapeutic Option for TB Infection: Evidence to Date

Tuberculosis (TB) is the most deadly infectious disease globally. Although most individuals achieve a cure, a substantial portion develop multi-drug resistant TB which is exceedingly difficult to treat, and the number of effective agents is dwindling. Development of new anti-tubercular medications is imperative to combat existing drug resistance and accelerate global eradication of TB. Pretomanid (PA-824) represents one of the newest drug classes (ie, nitroimidazooxazines) approved in 2019 by the United States Food and Drug Administration as part of a multi-drug regimen (with bedaquiline and linezolid, BPaL) and recommended by the World Health Organization (WHO) to treat extensively-resistant (XR-TB) and multi-drug resistant tuberculosis (MDR-TB). Approval was granted through the FDA's Limited Population Pathway for Antibacterial and Antifungal Drugs, which accelerates approval for antimicrobial drugs used to treat life-threatening or serious infections in a limited population with unmet need. This review details the pharmacology, efficacy, and safety of this new agent and describes evidence to date for its role in the treatment of drug resistant TB including published, ongoing, and planned studies.

Metabolic bifunctionality of Rv0812 couples folate and peptidoglycan biosynthesis in Mycobacterium tuberculosis

Comparative sequence analysis has enabled the annotation of millions of genes from organisms across the evolutionary tree. However, this approach has inherently biased the annotation of phylogenetically ubiquitous, rather than species-specific, functions. The ecologically unusual pathogen Mycobacterium tuberculosis (Mtb) has evolved in humans as its sole reservoir and emerged as the leading bacterial cause of death worldwide. However, the physiological factors that define Mtb’s pathogenicity are poorly understood. Here, we report the structure and function of a protein that is required for optimal in vitro fitness and bears homology to two distinct enzymes, Rv0812. Despite diversification of related orthologues into biochemically distinct enzyme families, rv0812 encodes a single active site with aminodeoxychorismate lyase and D–amino acid transaminase activities. The mutual exclusivity of substrate occupancy in this active site mediates coupling between nucleic acid and cell wall biosynthesis, prioritizing PABA over D-Ala/D-Glu biosynthesis. This bifunctionality reveals a novel, enzymatically encoded fail-safe mechanism that may help Mtb and other bacteria couple replication and division.

Pan and Core Genome Analysis of 183 Mycobacterium tuberculosis Strains Revealed a High Inter-Species Diversity among the Human Adapted Strains

Tuberculosis (TB) is an airborne communicable disease with high morbidity and mortality rates, especially in developing countries. The causal agents of TB belong to the complex Mycobacterium tuberculosis (MTBc), which is composed of different human and animal TB associated species. Some animal associated species have zoonotic potential and add to the burden of TB management. The BCG ("Bacillus Calmette-Guérin") vaccine is widely used for the prevention against TB, but its use is limited in immunocompromised patients and animals due to the adverse effects and disseminated life-threatening complications. In this study, we aimed to carry out a comparative genome analysis between the human adapted species including BCG vaccine strains to identify and pinpoint the conserved genes related to the virulence across all the species, which could add a new value for vaccine development. For this purpose, the sequences of 183 Mycobacterium tuberculosis (MTB) strains were retrieved from the freely available WGS dataset at NCBI. The species included: 168 sensu stricto MTB species with other human MTB complex associated strains: M. tuberculosis var. africanum (3), M. tuberculosis var. bovis (2 draft genomes) and 10 BCG species, which enabled the analysis of core genome which contains the conserved genes and some virulence factor determinants. Further, a phylogenetic tree was constructed including the genomes of human (183); animals MTB adapted strains (6) and the environmental Mycobacterium strain "M. canettii". Our results showed that the core genome consists of 1166 conserved genes among these species, which represents a small portion of the pangenome (7036 genes). The remaining genes in the pangenome (5870) are accessory genes, adding a high inter-species diversity. Further, the core genome includes several virulence-associated genes and this could explain the rare infectiousness potential of some attenuated vaccine strains in some patients. This study reveals that low number of conserved genes in human adapted MTBc species and high inter-species diversity of the pan-genome could be considered for vaccine candidate development.

Targeting amino acid metabolism of Mycobacterium tuberculosis for developing inhibitors to curtail its survival

Tuberculosis caused by the bacterium, Mycobacterium tuberculosis (Mtb), continues to remain one of the most devastating infectious diseases afflicting humans. Although there are several drugs for treating tuberculosis available currently, the emergence of the drug resistant forms of this pathogen has made its treatment and eradication a challenging task. While the replication machinery, protein synthesis and cell wall biogenesis of Mtb have been targeted often for anti-tubercular drug development a number of essential metabolic pathways crucial to its survival have received relatively less attention. In this context a number of amino acid biosynthesis pathways have recently been shown to be essential for the survival and pathogenesis of Mtb. Many of these pathways and or their key enzymes homologs are absent in humans hence they could be harnessed for anti-tubercular drug development. In this review, we describe comprehensively the amino acid metabolic pathways essential in Mtb and the key enzymes involved therein that are being investigated for developing inhibitors that compromise the survival and pathogenesis caused by this pathogen.

DEG 15, an update of the Database of Essential Genes that includes built-in analysis tools

Essential genes refer to genes that are required by an organism to survive under specific conditions. Studies of the minimal-gene-set for bacteria have elucidated fundamental cellular processes that sustain life. The past five years have seen a significant progress in identifying human essential genes, primarily due to the successful use of CRISPR/Cas9 in various types of human cells. DEG 15, a new release of the Database of Essential Genes (www.essentialgene.org), has provided major advancements, compared to DEG 10. Specifically, the number of eukaryotic essential genes has increased by more than fourfold, and that of prokaryotic ones has more than doubled. Of note, the human essential-gene number has increased by more than tenfold. Moreover, we have developed built-in analysis modules by which users can perform various analyses, such as essential-gene distributions between bacterial leading and lagging strands, sub-cellular localization distribution, enrichment analysis of gene ontology and KEGG pathways, and generation of Venn diagrams to compare and contrast gene sets between experiments. Additionally, the database offers customizable BLAST tools for performing species- and experiment-specific BLAST searches. Therefore, DEG comprehensively harbors updated human-curated essential-gene records among prokaryotes and eukaryotes with built-in tools to enhance essential-gene analysis.

CRISPR Interference (CRISPRi) for Targeted Gene Silencing in Mycobacteria

The genetic basis for Mycobacterium tuberculosis pathogenesis is incompletely understood. One reason for this knowledge gap is the relative difficulty of genetic manipulation of M. tuberculosis. To close this gap, we recently developed a robust CRISPR interference (CRISPRi) platform for programmable gene silencing in mycobacteria. In this chapter, we: (1) discuss some of the advantages and disadvantages of CRISPRi relative to more traditional genetic approaches; and (2) provide a protocol for the application of CRISPRi to reduce transcription of target genes in mycobacteria.

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.

RNA Identification of PRIME Cells Predicting Rheumatoid Arthritis Flares

BACKGROUND

Rheumatoid arthritis, like many inflammatory diseases, is characterized by episodes of quiescence and exacerbation (flares). The molecular events leading to flares are unknown.

METHODS

We established a clinical and technical protocol for repeated home collection of blood in patients with rheumatoid arthritis to allow for longitudinal RNA sequencing (RNA-seq). Specimens were obtained from 364 time points during eight flares over a period of 4 years in our index patient, as well as from 235 time points during flares in three additional patients. We identified transcripts that were differentially expressed before flares and compared these with data from synovial single-cell RNA-seq. Flow cytometry and sorted-blood-cell RNA-seq in additional patients were used to validate the findings.

RESULTS

Consistent changes were observed in blood transcriptional profiles 1 to 2 weeks before a rheumatoid arthritis flare. B-cell activation was followed by expansion of circulating CD45-CD31-PDPN+ preinflammatory mesenchymal, or PRIME, cells in the blood from patients with rheumatoid arthritis; these cells shared features of inflammatory synovial fibroblasts. Levels of circulating PRIME cells decreased during flares in all 4 patients, and flow cytometry and sorted-cell RNA-seq confirmed the presence of PRIME cells in 19 additional patients with rheumatoid arthritis.

CONCLUSIONS

Longitudinal genomic analysis of rheumatoid arthritis flares revealed PRIME cells in the blood during the period before a flare and suggested a model in which these cells become activated by B cells in the weeks before a flare and subsequently migrate out of the blood into the synovium. (Funded by the National Institutes of Health and others.).

Mycobacterium smegmatis HtrA Blocks the Toxic Activity of a Putative Cell Wall Amidase

Mycobacterium tuberculosis, the causative agent of tuberculosis, withstands diverse environmental stresses in the host. The periplasmic protease HtrA is required only to survive extreme conditions in most bacteria but is predicted to be essential for normal growth in mycobacteria. We confirm that HtrA is indeed essential in Mycobacterium smegmatis and interacts with another essential protein of unknown function, LppZ. However, the loss of any of three unlinked genes, including those encoding Ami3, a peptidoglycan muramidase, and Pmt, a mannosyltransferase, suppresses the essentiality of both HtrA and LppZ, indicating the functional relevance of these genes' protein products. Our data indicate that HtrA-LppZ is required to counteract the accumulation of active Ami3, which is toxic under the stabilizing influence of Pmt-based mannosylation. This suggests that HtrA-LppZ blocks the toxicity of a cell wall enzyme to maintain mycobacterial homeostasis.

Role of the unique, non-essential phosphatidylglycerol::prolipoprotein diacylglyceryl transferase (Lgt) in Corynebacterium glutamicum

Bacterial lipoproteins are secreted proteins that are post-translationally lipidated. Following synthesis, preprolipoproteins are transported through the cytoplasmic membrane via the Sec or Tat translocon. As they exit the transport machinery, they are recognized by a phosphatidylglycerol::prolipoprotein diacylglyceryl transferase (Lgt), which converts them to prolipoproteins by adding a diacylglyceryl group to the sulfhydryl side chain of the invariant Cys₊₁ residue. Lipoprotein signal peptidase (LspA or signal peptidase II) subsequently cleaves the signal peptide, liberating the α-amino group of Cys₊₁, which can eventually be further modified. Here, we identified the lgt and lspA genes from Corynebacterium glutamicum and found that they are unique but not essential. We found that Lgt is necessary for the acylation and membrane anchoring of two model lipoproteins expressed in this species: MusE, a C. glutamicum maltose-binding lipoprotein, and LppX, a Mycobacterium tuberculosis lipoprotein. However, Lgt is not required for these proteins' signal peptide cleavage, or for LppX glycosylation. Taken together, these data show that in C. glutamicum the association of some lipoproteins with membranes through the covalent attachment of a lipid moiety is not essential for further post-translational modification.

A systematic evaluation of Mycobacterium tuberculosis Genome-Scale Metabolic Networks

Metabolism underpins the pathogenic strategy of the causative agent of TB, Mycobacterium tuberculosis (Mtb), and therefore metabolic pathways have recently re-emerged as attractive drug targets. A powerful approach to study Mtb metabolism as a whole, rather than just individual enzymatic components, is to use a systems biology framework, such as a Genome-Scale Metabolic Network (GSMN) that allows the dynamic interactions of all the components of metabolism to be interrogated together. Several GSMNs networks have been constructed for Mtb and used to study the complex relationship between the Mtb genotype and its phenotype. However, the utility of this approach is hampered by the existence of multiple models, each with varying properties and performances. Here we systematically evaluate eight recently published metabolic models of Mtb-H37Rv to facilitate model choice. The best performing models, sMtb2018 and iEK1011, were refined and improved for use in future studies by the TB research community.

Growth promotion in Corynebacterium glutamicum by overexpression of the NCgl2986 gene encoding a protein homologous to peptidoglycan amidases

We previously reported the extracellular production of antibody fragment Fab by Corynebacterium glutamicum. In the course of searching for genes which improve the secretion efficiency of Fab, we coincidentally found that the final growth increased significantly when the NCgl2986 gene encoding an amidase-like protein was overexpressed. This effect was observed when cells were grown on the production medium MMTG, which contains high concentrations of glucose and neutralizing agent CaCO₃, but not on MMTG without CaCO₃ or Lennox medium. Not only turbidity but also dry cell weight was increased by NCgl2986 overexpression, although the growth rate was not affected. It was recently reported that the Mycobacterium tuberculosis homolog Rv3915 functions as an activator of MurA protein, which catalyzes the initial step of peptidoglycan synthesis. Growth promotion was also observed when the MurA protein was overproduced. His-tagged NCgl2986 protein was purified, but its peptidoglycan hydrolyzing activity could not be detected. These results suggest that NCgl2986 promotes cell growth by activating the peptidoglycan synthetic pathway.

Genome-wide identification of essential genes in Mycobacterium intracellulare by transposon sequencing - Implication for metabolic remodeling

The global incidence of the human nontuberculous mycobacteria (NTM) disease is rapidly increasing. However, knowledge of gene essentiality under optimal growth conditions and conditions relevant to the natural ecology of NTM, such as hypoxia, is lacking. In this study, we utilized transposon sequencing to comprehensively identify genes essential for growth in Mycobacterium intracellulare. Of 5126 genes of M. intracellulare ATCC13950, 506 genes were identified as essential genes, of which 280 and 158 genes were shared with essential genes of M. tuberculosis and M. marinum, respectively. The shared genes included target genes of existing antituberculous drugs including SQ109, which targets the trehalose monomycolate transporter MmpL3. From 175 genes showing decreased fitness as conditionally essential under hypoxia, preferential carbohydrate metabolism including gluconeogenesis, glyoxylate cycle and succinate production was suggested under hypoxia. Virulence-associated genes including proteasome system and mycothiol redox system were also identified as conditionally essential under hypoxia, which was further supported by the higher effective suppression of bacterial growth under hypoxia compared to aerobic conditions in the presence of these inhibitors. This study has comprehensively identified functions essential for growth of M. intracellulare under conditions relevant to the host environment. These findings provide critical functional genomic information for drug discovery.

Transposon mutagenesis in Mycobacterium kansasii links a small RNA gene to colony morphology and biofilm formation and identifies 9,885 intragenic insertions that do not compromise colony outgrowth

Mycobacterium kansasii (Mk) is a resilient opportunistic human pathogen that causes tuberculosis‐like chronic pulmonary disease and mortality stemming from comorbidities and treatment failure. The standard treatment of Mk infections requires costly, long‐term, multidrug courses with adverse side effects. The emergence of drug‐resistant isolates further complicates the already challenging drug therapy regimens and threatens to compromise the future control of Mk infections. Despite the increasingly recognized global burden of Mk infections, the biology of this opportunistic pathogen remains essentially unexplored. In particular, studies reporting gene function or generation of defined mutants are scarce. Moreover, no transposon (Tn) mutagenesis tool has been validated for use in Mk, a situation limiting the repertoire of genetic approaches available to accelerate the dissection of gene function and the generation of gene knockout mutants in this poorly characterized pathogen. In this study, we validated the functionality of a powerful Tn mutagenesis tool in Mk and used this tool in conjunction with a forward genetic screen to establish a previously unrecognized role of a conserved mycobacterial small RNA gene of unknown function in colony morphology features and biofilm formation. We also combined Tn mutagenesis with next‐generation sequencing to identify 12,071 Tn insertions that do not compromise viability in vitro. Finally, we demonstrated the susceptibility of the Galleria mellonella larva to Mk, setting the stage for further exploration of this simple and economical infection model system to the study of this pathogen.

Comprehensive analysis of iron utilization by Mycobacterium tuberculosis

Iron is essential for nearly all bacterial pathogens, including Mycobacterium tuberculosis (Mtb), but is severely limited in the human host. To meet its iron needs, Mtb secretes siderophores, small molecules with high affinity for iron, and takes up iron-loaded mycobactins (MBT) and carboxymycobactins (cMBT), from the environment. Mtb is also capable of utilizing heme and hemoglobin which contain more than 70% of the iron in the human body. However, many components of these iron acquisition pathways are still unknown. In this study, a high-density transposon mutagenesis coupled with deep sequencing (TnSeq) showed that Mtb exhibits nearly opposite requirements for 165 genes in the presence of heme and hemoglobin versus MBT and cMBT as iron sources. The ESX-3 secretion system was assessed as essential for siderophore-mediated iron uptake and, surprisingly, also for heme utilization by Mtb. Predictions derived from the TnSeq analysis were validated by growth experiments with isogenic Mtb mutants. These results showed that (i) the efflux pump MmpL5 plays a dominant role in siderophore secretion, (ii) the Rv2047c protein is essential for growth of Mtb in the presence of mycobactin, and (iii) the transcriptional repressor Zur is required for heme utilization by Mtb. The novel genetic determinants of iron utilization revealed in this study will stimulate further experiments in this important area of Mtb physiology.

Comparative genomics shows differences in the electron transport and carbon metabolic pathways of Mycobacterium africanum relative to Mycobacterium tuberculosis and suggests an adaptation to low oxygen tension

The geographically restricted Mycobacterium africanum lineages (MAF) are primarily found in West Africa, where they account for a significant proportion of tuberculosis. Despite this phenomenon, little is known about the co-evolution of these ancient lineages with West Africans. MAF and M. tuberculosis sensu stricto lineages (MTB) differ in their clinical, in vitro and in vivo characteristics for reasons not fully understood. Therefore, we compared genomes of 289 MAF and 205 MTB clinical isolates from the 6 main human-adapted M. tuberculosis complex lineages, for mutations in their Electron Transport Chain and Central Carbon Metabolic pathway in order to explain these metabolic differences. Furthermore, we determined, in silico, whether each mutation could affect the function of genes encoding enzymes in these pathways. We found more mutations with the potential to affect enzymes in these pathways in MAF lineages compared to MTB lineages. We also found that similar mutations occurred in these pathways between MAF and some MTB lineages. Generally, our findings show further differences between MAF and MTB lineages that may have contributed to the MAF clinical and growth phenotype and indicate potential adaptation of MAF lineages to a distinct ecological niche, which we suggest includes areas characterized by low oxygen tension.

Global Assessment of Mycobacterium avium subsp. hominissuis Genetic Requirement for Growth and Virulence

Nontuberculous mycobacterial infections caused by the opportunistic pathogen Mycobacterium avium subsp. hominissuis (MAH) are currently receiving renewed attention due to increased incidence combined with difficult treatment. Insights into the disease-causing mechanisms of this species have been hampered by difficulties in genetic manipulation of the bacteria. Here, we identified and sequenced a highly transformable, virulent MAH clinical isolate susceptible to high-density transposon mutagenesis, facilitating global gene disruption and subsequent investigation of MAH gene function. By transposon insertion sequencing (TnSeq) of this strain, we defined the MAH genome-wide genetic requirement for virulence and in vitro growth and organized ∼3,500 identified transposon mutants for hypothesis-driven research. The majority (96%) of the genes we identified as essential for MAH in vitro had a mutual ortholog in the related and highly virulent Mycobacterium tuberculosis (Mtb). However, passaging our library through a mouse model of infection revealed a substantial number (54% of total hits) of novel virulence genes. More than 97% of the MAH virulence genes had a mutual ortholog in Mtb Finally, we validated novel genes required for successful MAH infection: one encoding a probable major facilitator superfamily (MFS) transporter and another encoding a hypothetical protein located in the immediate vicinity of six other identified virulence genes. In summary, we provide new, fundamental insights into the underlying genetic requirement of MAH for growth and host infection.IMPORTANCE Pulmonary disease caused by nontuberculous mycobacteria is increasing worldwide. The majority of these infections are caused by the Mycobacterium avium complex (MAC), whereof >90% are due to Mycobacterium avium subsp. hominissuis (MAH). Treatment of MAH infections is currently difficult, with a combination of antibiotics given for at least 12 months. To control MAH by improved therapy, prevention, and diagnostics, we need to understand the underlying mechanisms of infection. Here, we provide crucial insights into MAH's global genetic requirements for growth and infection. We find that the vast majority of genes required for MAH growth and virulence (96% and 97%, respectively) have mutual orthologs in the tuberculosis-causing pathogen M. tuberculosis (Mtb). However, we also find growth and virulence genes specific to MAC species. Finally, we validate novel mycobacterial virulence factors that might serve as future drug targets for MAH-specific treatment or translate to broader treatment of related mycobacterial diseases.

Statistical analysis of variability in TnSeq data across conditions using zero-inflated negative binomial regression

Background

Deep sequencing of transposon mutant libraries (or TnSeq) is a powerful method for probing essentiality of genomic loci under different environmental conditions. Various analytical methods have been described for identifying conditionally essential genes whose tolerance for insertions varies between two conditions. However, for large-scale experiments involving many conditions, a method is needed for identifying genes that exhibit significant variability in insertions across multiple conditions.

Results

In this paper, we introduce a novel statistical method for identifying genes with significant variability of insertion counts across multiple conditions based on Zero-Inflated Negative Binomial (ZINB) regression. Using likelihood ratio tests, we show that the ZINB distribution fits TnSeq data better than either ANOVA or a Negative Binomial (in a generalized linear model). We use ZINB regression to identify genes required for infection of M. tuberculosis H37Rv in C57BL/6 mice. We also use ZINB to perform a analysis of genes conditionally essential in H37Rv cultures exposed to multiple antibiotics.

Conclusions

Our results show that, not only does ZINB generally identify most of the genes found by pairwise resampling (and vastly out-performs ANOVA), but it also identifies additional genes where variability is detectable only when the magnitudes of insertion counts are treated separately from local differences in saturation, as in the ZINB model.