Statistical analysis of genetic interactions in Tn-Seq data

Genome-Wide CRISPRi Screening of Key Genes for Recombinant Protein Expression in Bacillus Subtilis

Bacillus subtilis is an industrially important microorganism that is often used as a microbial cell factory for the production of recombinant proteins due to its food safety, rapid growth, and powerful secretory capacity. However, the lack of data on functional genes related to recombinant protein production has hindered the further development of B. subtilis cell factories. Here, a strategy combining genome-wide CRISPRi screening and targeted CRISPRa activation to enhance recombinant protein expression is proposed. First, a CRISPRi library covering a total of 4225 coding genes (99.7%) in the B. subtilis genome and built the corresponding high-throughput screening methods is constructed. Twelve key genes for recombinant protein expression are identified, including targets without relevant functional annotations. Meanwhile, the transcription of recombinant protein genes by CRISPRa is up-regulated. These screened or selected genes can be easily applied to metabolic engineering by constructing sgRNA arrays. The relationship between differential pathways and recombinant protein expression in engineered strains by transcriptome analysis is also revealed. High-density fermentation and generalisability validation results prove the reliability of the strategy. This method can be extended to other industrial hosts to support functional gene annotation and the design of novel cell factories.

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.

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.

Activity-based annotation: the emergence of systems biochemistry

Current tools to annotate protein function have failed to keep pace with the speed of DNA sequencing and exponentially growing number of proteins of unknown function (PUFs). A major contributing factor to this mismatch is the historical lack of high-throughput methods to experimentally determine biochemical activity. Activity-based methods, such as activity-based metabolite and protein profiling, are emerging as new approaches for unbiased, global, biochemical annotation of protein function. In this review, we highlight recent experimental, activity-based approaches that offer new opportunities to determine protein function in a biologically agnostic and systems-level manner.

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

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

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

Model-based identification of conditionally-essential genes from transposon-insertion sequencing data

The understanding of bacterial gene function has been greatly enhanced by recent advancements in the deep sequencing of microbial genomes. Transposon insertion sequencing methods combines next-generation sequencing techniques with transposon mutagenesis for the exploration of the essentiality of genes under different environmental conditions. We propose a model-based method that uses regularized negative binomial regression to estimate the change in transposon insertions attributable to gene-environment changes in this genetic interaction study without transformations or uniform normalization. An empirical Bayes model for estimating the local false discovery rate combines unique and total count information to test for genes that show a statistically significant change in transposon counts. When applied to RB-TnSeq (randomized barcode transposon sequencing) and Tn-seq (transposon sequencing) libraries made in strains of Caulobacter crescentus using both total and unique count data the model was able to identify a set of conditionally beneficial or conditionally detrimental genes for each target condition that shed light on their functions and roles during various stress conditions.

Analysis of Gene Essentiality from TnSeq Data Using Transit

TnSeq, or sequencing of transposon insertion libraries, has proven to be a valuable method for probing the functions of genes in a wide range of bacteria. TnSeq has found many applications for studying genes involved in core functions (such as cell division or metabolism), stress response, virulence, etc., as well as to identify potential drug targets. Two of the most commonly used transposons in practice are Himar1, which inserts randomly at TA dinucleotides, and Tn5, which can insert more broadly throughout the genome. These insertions cause putative gene function disruption, and clones with insertions in genes that cannot tolerate disruption (in a given condition) are eliminated from the population. Deep sequencing can be used to efficiently profile the surviving members, with insertions in genes that can be inferred to be non-essential. Data from TnSeq experiments (i.e. transposon insertion counts at specific genomic locations) is inherently noisy, making rigorous statistical analysis (e.g. quantifying significance) challenging. In this chapter, we describe Transit, a Python-based software package for analyzing TnSeq data that combines a variety of data processing tools, quality assessment methods, and analytical algorithms for identifying essential (or conditionally essential) genes.

Loss of YhcB results in dysregulation of coordinated peptidoglycan, LPS and phospholipid synthesis during Escherichia coli cell growth

The cell envelope is essential for viability in all domains of life. It retains enzymes and substrates within a confined space while providing a protective barrier to the external environment. Destabilising the envelope of bacterial pathogens is a common strategy employed by antimicrobial treatment. However, even in one of the best studied organisms, Escherichia coli, there remain gaps in our understanding of how the synthesis of the successive layers of the cell envelope are coordinated during growth and cell division. Here, we used a whole-genome phenotypic screen to identify mutants with a defective cell envelope. We report that loss of yhcB, a conserved gene of unknown function, results in loss of envelope stability, increased cell permeability and dysregulated control of cell size. Using whole genome transposon mutagenesis strategies, we report the comprehensive genetic interaction network of yhcB, revealing all genes with a synthetic negative and a synthetic positive relationship. These genes include those previously reported to have a role in cell envelope biogenesis. Surprisingly, we identified genes previously annotated as essential that became non-essential in a ΔyhcB background. Subsequent analyses suggest that YhcB functions at the junction of several envelope biosynthetic pathways coordinating the spatiotemporal growth of the cell, highlighting YhcB as an as yet unexplored antimicrobial target.

All life depends on a cell envelope to enclose the chemical reactions that make life possible. But how do cell envelopes grow? How each component of the cell envelope is incorporated into the envelope at the correct amount, in the correct place, and at the correct time, to prevent cell death, has been a long-standing question in bacteriology. Using a unique combination of high throughput chemical genetic screens we identified yhcB, a conserved gene of unknown function, required for the maintenance of cell envelope integrity in Escherichia coli. Loss of YhcB results in aberrant cell size driven by the production of excess membrane phospholipids. Subsequent molecular and biochemical analyses suggest YhcB influences the spatiotemporal biogenesis of LPS, peptidoglycan and membrane phospholipids. Our data indicate YhcB is a key regulator of cell envelope growth in Gram-negative bacteria playing a crucial role in coordinating cell width, elongation, and division to maintain cell envelope integrity.

Gene Fitness of Azotobacter vinelandii under Diazotrophic Growth

Azotobacter vinelandii is a nitrogen-fixing free-living soil microbe that has been studied for decades in relation to biological nitrogen fixation (BNF). It is highly amenable to genetic manipulation, helping to unravel the intricate importance of different proteins involved in the process of BNF, including the biosynthesis of cofactors that are essential to assembling the complex metal cofactors that catalyze the difficult reaction of nitrogen fixation. Additionally, A. vinelandii accomplishes this feat while growing as an obligate aerobe, differentiating it from many of the nitrogen-fixing bacteria that are associated with plant roots. The ability to function in the presence of oxygen makes A. vinelandii suitable for application in various potential biotechnological schemes. In this study, we employed transposon sequencing (Tn-seq) to measure the fitness defects associated with disruptions of various genes under nitrogen-fixing dependent growth, versus growth with extraneously provided urea as a nitrogen source. The results allowed us to probe the importance of more than 3,800 genes, revealing that many genes previously believed to be important, can be successfully disrupted without impacting cellular fitness. IMPORTANCE These results provide insights into the functional redundancy in A. vinelandii, while also providing a direct measure of fitness for specific genes associated with the process of BNF. These results will serve as a valuable reference tool in future studies to uncover the mechanisms that govern this process.

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.

Role of LmeA, a Mycobacterial Periplasmic Protein, in Maintaining the Mannosyltransferase MptA and Its Product Lipomannan under Stress

Mycobacteria differentially regulate the cellular amounts of lipoglycans in response to environmental changes, but the molecular mechanisms of this regulation remain unknown. Here, we demonstrate that cellular lipoarabinomannan (LAM) levels rapidly decline under two stress conditions, stationary growth phase and nutrient starvation, while the levels of another related lipoglycan, lipomannan (LM), stay relatively constant. The persistence of LM under stress correlated with the maintenance of two key mannosyltransferases, MptA and MptC, in the LM biosynthetic pathway. We further showed that the stress exposures lead to the upregulation of lmeA gene expression and that the periplasmic protein LmeA plays a key role in maintaining the enzyme MptA and its product LM under stress conditions. These findings reveal new aspects of how lipoglycan biosynthesis is regulated under stress conditions in mycobacteria.

The mycobacterial cell envelope has a diderm structure, composed of an outer mycomembrane, an arabinogalactan-peptidoglycan cell wall, a periplasm, and an inner membrane. Lipomannan (LM) and lipoarabinomannan (LAM) are structural and immunomodulatory components of this cell envelope. LM/LAM biosynthesis involves a number of mannosyltransferases and acyltransferases, and MptA is an α1,6-mannosyltransferase involved in the final extension of the mannan chain. Recently, we reported the periplasmic protein LmeA being involved in the maturation of the mannan backbone in Mycobacterium smegmatis. Here, we examined the role of LmeA under stress conditions. We found that lmeA transcription was upregulated under two stress conditions: stationary growth phase and nutrient starvation. Under both conditions, LAM was decreased, but LM was relatively stable, suggesting that maintaining the cellular level of LM under stress is important. Surprisingly, the protein levels of MptA were decreased in an lmeA deletion (ΔlmeA) mutant under both stress conditions. The transcript levels of mptA in the ΔlmeA mutant were similar to or even higher than those in the wild type, indicating that the decrease of MptA protein was a posttranscriptional event. The ΔlmeA mutant was unable to maintain the cellular level of LM under stress, consistent with the decrease in MptA. Even during active growth, overexpression of LmeA led the cells to produce more LM and become more resistant to several antibiotics. Altogether, our study reveals the roles of LmeA in the homeostasis of the MptA mannosyltransferase, particularly under stress conditions, ensuring the stable expression of LM and the maintenance of cell envelope integrity.

IMPORTANCE Mycobacteria differentially regulate the cellular amounts of lipoglycans in response to environmental changes, but the molecular mechanisms of this regulation remain unknown. Here, we demonstrate that cellular lipoarabinomannan (LAM) levels rapidly decline under two stress conditions, stationary growth phase and nutrient starvation, while the levels of another related lipoglycan, lipomannan (LM), stay relatively constant. The persistence of LM under stress correlated with the maintenance of two key mannosyltransferases, MptA and MptC, in the LM biosynthetic pathway. We further showed that the stress exposures lead to the upregulation of lmeA gene expression and that the periplasmic protein LmeA plays a key role in maintaining the enzyme MptA and its product LM under stress conditions. These findings reveal new aspects of how lipoglycan biosynthesis is regulated under stress conditions in mycobacteria.

A decade of advances in transposon-insertion sequencing

It has been 10 years since the introduction of modern transposon-insertion sequencing (TIS) methods, which combine genome-wide transposon mutagenesis with high-throughput sequencing to estimate the fitness contribution or essentiality of each genetic component in a bacterial genome. Four TIS variations were published in 2009: transposon sequencing (Tn-Seq), transposon-directed insertion site sequencing (TraDIS), insertion sequencing (INSeq) and high-throughput insertion tracking by deep sequencing (HITS). TIS has since become an important tool for molecular microbiologists, being one of the few genome-wide techniques that directly links phenotype to genotype and ultimately can assign gene function. In this Review, we discuss the recent applications of TIS to answer overarching biological questions. We explore emerging and multidisciplinary methods that build on TIS, with an eye towards future applications.

In this Review, several experts discuss progress in the decade since the development of transposon-based approaches for bacterial genetic screens. They describe how advances in both experimental technologies and analytical strategies are resulting in insights into diverse biological processes.

Selection or drift: The population biology underlying transposon insertion sequencing experiments

Transposon insertion sequencing methods such as Tn-seq revolutionized microbiology by allowing the identification of genomic loci that are critical for viability in a specific environment on a genome-wide scale. While powerful, transposon insertion sequencing suffers from limited reproducibility when different analysis methods are compared. From the perspective of population biology, this may be explained by changes in mutant frequency due to chance (drift) rather than differential fitness (selection). Here, we develop a mathematical model of the population biology of transposon insertion sequencing experiments, i.e. the changes in size and composition of the transposon-mutagenized population during the experiment. We use this model to investigate mutagenesis, the growth of the mutant library, and its passage through bottlenecks. Specifically, we study how these processes can lead to extinction of individual mutants depending on their fitness and the distribution of fitness effects (DFE) of the entire mutant population. We find that in typical in vitro experiments few mutants with high fitness go extinct. However, bottlenecks of a size that is common in animal infection models lead to so much random extinction that a large number of viable mutants would be misclassified. While mutants with low fitness are more likely to be lost during the experiment, mutants with intermediate fitness are expected to be much more abundant and can constitute a large proportion of detected hits, i.e. false positives. Thus, incorporating the DFEs of randomly generated mutations in the analysis may improve the reproducibility of transposon insertion experiments, especially when strong bottlenecks are encountered.

Global phenotypic profiling identifies a conserved actinobacterial cofactor for a bifunctional PBP-type cell wall synthase

Members of the Corynebacterineae suborder of Actinobacteria have a unique cell surface architecture and, unlike most well-studied bacteria, grow by tip-extension. To investigate the distinct morphogenic mechanisms shared by these organisms, we performed a genome-wide phenotypic profiling analysis using Corynebacterium glutamicum as a model. A high-density transposon mutagenized library was challenged with a panel of antibiotics and other stresses. The fitness of mutants in each gene under each condition was then assessed by transposon-sequencing. Clustering of the resulting phenotypic fingerprints revealed a role for several genes of previously unknown function in surface biogenesis. Further analysis identified CofA (Cgp_0016) as an interaction partner of the peptidoglycan synthase PBP1a that promotes its stable accumulation at sites of polar growth. The related Mycobacterium tuberculosis proteins were also found to interact, highlighting the utility of our dataset for uncovering conserved principles of morphogenesis for this clinically relevant bacterial suborder.

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.

In vitro exploration of the Xanthomonas hortorum pv. vitians genome using transposon insertion sequencing and comparative genomics to discriminate between core and contextual essential genes

The essential genome of a bacterium encompasses core genes associated with basic cellular processes and conditionally essential genes dependent upon environmental conditions or the genetic context. Comprehensive knowledge of those gene sets allows for a better understanding of fundamental bacterial biology and offers new perspectives for antimicrobial drug research against detrimental bacteria such as pathogens. We investigated the essential genome of Xanthomonas hortorum pv. vitians, a gammaproteobacterial plant pathogen of lettuce (Lactuca sativa L.) which belongs to the plant-pathogen reservoir genus Xanthomonas and is affiliated to the family Xanthomonadaceae. No practical means of disease control or prevention against this pathogen is currently available, and its molecular biology is virtually unknown. To reach a comprehensive overview of the essential genome of X. hortorum pv. vitians LM16734, we developed a mixed approach combining high-quality full genome sequencing, saturated transposon insertion sequencing (Tn-Seq) in optimal growth conditions, and coupled computational analyses such as comparative genomics, synteny assessment and phylogenomics. Among the 370 essential loci identified by Tn-Seq, a majority was bound to critical cell processes conserved across bacteria. The remaining genes were either related to specific ecological features of Xanthomonas or Xanthomonadaceae species, or acquired through horizontal gene transfer of mobile genetic elements and associated with ancestral parasitic gene behaviour and bacterial defence systems. Our study sheds new light on our usual concepts about gene essentiality and is pioneering in the molecular and genomic study of X. hortorum pv. vitians.

A guide to Mycobacterium mutagenesis

The genus Mycobacterium includes several pathogens that cause severe disease in humans, like Mycobacterium tuberculosis (M. tb), the infectious agent causing tuberculosis. Genetic tools to engineer mycobacterial genomes, in a targeted or random fashion, have provided opportunities to investigate M. tb infection and pathogenesis. Furthermore, they have allowed the identification and validation of potential targets for the diagnosis, prevention, and treatment of tuberculosis. This review describes the various methods that are available for the generation of mutants in Mycobacterium species, focusing specifically on tools for altering slow-growing mycobacteria from the M. tb complex. Among others, it incorporates the recent new molecular biological technologies (e.g. ORBIT) to rapidly and/or genome-wide comprehensively obtain targeted mutants in mycobacteria. As such, this review can be used as a guide to select the appropriate genetic tools to generate mycobacterial mutants of interest, which can be used as tools to aid understanding of M. tb infection or to help developing TB intervention strategies.

SB Driver Analysis: a Sleeping Beauty cancer driver analysis framework for identifying and prioritizing experimentally actionable oncogenes and tumor suppressors

Cancer driver prioritization for functional analysis of potential actionable therapeutic targets is a significant challenge. Meta-analyses of mutated genes across different human cancer types for driver prioritization has reaffirmed the role of major players in cancer, including KRAS, TP53 and EGFR, but has had limited success in prioritizing genes with non-recurrent mutations in specific cancer types. Sleeping Beauty (SB) insertional mutagenesis is a powerful experimental gene discovery framework to define driver genes in mouse models of human cancers. Meta-analyses of SB datasets across multiple tumor types is a potentially informative approach to prioritize drivers, and complements efforts in human cancers. Here, we report the development of SB Driver Analysis, an in-silico method for defining cancer driver genes that positively contribute to tumor initiation and progression from population-level SB insertion data sets. We demonstrate that SB Driver Analysis computationally prioritizes drivers and defines distinct driver classes from end-stage tumors that predict their putative functions during tumorigenesis. SB Driver Analysis greatly enhances our ability to analyze, interpret and prioritize drivers from SB cancer datasets and will continue to substantially increase our understanding of the genetic basis of cancer.

Methods for Tn-Seq Analysis in Acinetobacter baumannii

Transposon insertion sequencing (Tn-seq) is a powerful method for identifying genes required for virtually any growth or survival trait in bacteria. The technology employs next-generation DNA sequencing to identify and quantify the relative abundances of individual transposon mutants within complex pools of such mutants. When applied to pools of thousands to millions of random transposon mutants grown under selective pressure, the technique can rapidly identify, at genome scale, the mutants and corresponding genes negatively or positively selected. This chapter presents core protocols for Tn-seq analysis of Acinetobacter baumannii: generation of a high-saturation random transposon mutant pool, and isolation and sequencing of transposon-genome junctions from such a pool for identifying and quantifying the individual mutants. With these tools, the researcher can address diverse biological questions by carrying out selective growth of a mutant pool followed by Tn-seq analysis to identify genotype-phenotype associations.

Unsupervised Learning Approach for Comparing Multiple Transposon Insertion Sequencing Studies

Transposon insertion sequencing (TIS) is a widely used technique for conducting genome-scale forward genetic screens in bacteria. However, few methods enable comparison of TIS data across multiple replicates of a screen or across independent screens, including screens performed in different organisms. Here, we introduce a post hoc analytic framework, comparative TIS (CompTIS), which utilizes unsupervised learning to enable meta-analysis of multiple TIS data sets. CompTIS first implements screen-level principal-component analysis (PCA) and clustering to identify variation between the TIS screens. This initial screen-level analysis facilitates the selection of related screens for additional analyses, reveals the relatedness of complex environments based on growth phenotypes measured by TIS, and provides a useful quality control step. Subsequently, PCA is performed on genes to identify loci whose corresponding mutants lead to concordant/discordant phenotypes across all or in a subset of screens. We used CompTIS to analyze published intestinal colonization TIS data sets from two vibrio species. Gene-level analyses identified both pan-vibrio genes required for intestinal colonization and conserved genes that displayed species-specific requirements. CompTIS is applicable to virtually any combination of TIS screens and can be implemented without regard to either the number of screens or the methods used for upstream data analysis.IMPORTANCE Forward genetic screens are powerful tools for functional genomics. The comparison of similar forward genetic screens performed in different organisms enables the identification of genes with similar or different phenotypes across organisms. Transposon insertion sequencing is a widely used method for conducting genome-scale forward genetic screens in bacteria, yet few bioinformatic approaches have been developed to compare the results of screen replicates and different screens conducted across species or strains. Here, we used principal-component analysis (PCA) and hierarchical clustering, two unsupervised learning approaches, to analyze the relatedness of multiple in vivo screens of pathogenic vibrios. This analytic framework reveals both shared pan-vibrio requirements for intestinal colonization and strain-specific dependencies. Our findings suggest that PCA-based analytics will be a straightforward widely applicable approach for comparing diverse transposon insertion sequencing screens.

High-throughput interaction screens illuminate the role of c-di-AMP in cyanobacterial nighttime survival

The broadly conserved signaling nucleotide cyclic di-adenosine monophosphate (c-di-AMP) is essential for viability in most bacteria where it has been studied. However, characterization of the cellular functions and metabolism of c-di-AMP has largely been confined to the class Bacilli, limiting our functional understanding of the molecule among diverse phyla. We identified the cyclase responsible for c-di-AMP synthesis and characterized the molecule’s role in survival of darkness in the model photosynthetic cyanobacterium Synechococcus elongatus PCC 7942. In addition to the use of traditional genetic, biochemical, and proteomic approaches, we developed a high-throughput genetic interaction screen (IRB-Seq) to determine pathways where the signaling nucleotide is active. We found that in S. elongatus c-di-AMP is produced by an enzyme of the diadenylate cyclase family, CdaA, which was previously unexplored experimentally. A cdaA-null mutant experiences increased oxidative stress and death during the nighttime portion of day-night cycles, in which potassium transport is implicated. These findings suggest that c-di-AMP is biologically active in cyanobacteria and has non-canonical roles in the phylum including oxidative stress management and day-night survival. The pipeline and analysis tools for IRB-Seq developed for this study constitute a quantitative high-throughput approach for studying genetic interactions.

Cyclic di-adenosine monophosphate (c-di-AMP) is a molecule that has significant roles in many microorganisms. This work shows the existence of c-di-AMP for the first time in photosynthetic microorganisms, cyanobacteria, and demonstrates its role in survival during the light-to-dark shifts that occur in day-night cycles. Despite the obvious importance of adaptation to these daily cycles for organisms that are fundamentally reliant on light, such as cyanobacteria, understanding of diurnal physiology is lacking because most cyanobacterial research is conducted during growth in constant light. To identify other players in c-di-AMP’s function we developed a low-cost and efficient method for finding interactions between genes. The technique combines one mutation, in this case for the gene that encodes the enzyme for synthesis of c-di-AMP, with thousands of other individual mutations to find pairwise interactions that affect fitness of the resulting mutants. Mutants are tagged with DNA barcodes to allow their survival to be easily tracked in a population of cells. The method enables us to place the function of c-di-AMP within the context of pathways previously known to be involved in day-night survival. Taken together, this work expands the known roles of c-di-AMP, improves our understanding of cyanobacterial survival in day-night cycles, and presents an improved approach for determining genetic interactions.

When Dicty Met Myco, a (Not So) Romantic Story about One Amoeba and Its Intracellular Pathogen

In recent years, Dictyostelium discoideum has become an important model organism to study the cell biology of professional phagocytes. This amoeba not only shares many molecular features with mammalian macrophages, but most of its fundamental signal transduction pathways are conserved in humans. The broad range of existing genetic and biochemical tools, together with its suitability for cell culture and live microscopy, make D. discoideum an ideal and versatile laboratory organism. In this review, we focus on the use of D. discoideum as a phagocyte model for the study of mycobacterial infections, in particular Mycobacterium marinum. We look in detail at the intracellular cycle of M. marinum, from its uptake by D. discoideum to its active or passive egress into the extracellular medium. In addition, we describe the molecular mechanisms that both the mycobacterial invader and the amoeboid host have developed to fight against each other, and compare and contrast with those developed by mammalian phagocytes. Finally, we introduce the methods and specific tools that have been used so far to monitor the D. discoideum-M. marinum interaction.