Combining Shigella Tn-seq data with gold-standard E. coli gene deletion data suggests rare transitions between essential and non-essential gene functionality

Essential Genes of Vibrio anguillarum and Other Vibrio spp. Guide the Development of New Drugs and Vaccines

Essential genes in bacterial pathogens are potential drug targets and vaccine candidates because disrupting their function is lethal. The development of new antibiotics, in addition to effective prevention measures such as vaccination, contributes to addressing the global problem of bacterial antibiotic resistance. The aim of this present study was to determine the essential genes of Vibrio anguillarum, a bacterial pathogen of aquatic animals, as a means to identify putative targets for novel drugs and to assist the prioritisation of potential vaccine candidates. Essential genes were characterised by a Tn-seq approach using the TnSC189 mariner transposon to construct a library of 52,662 insertion mutants. In total, 329 essential genes were identified, with 34.7% found within the core genome of this species; each of these genes represents a strong potential drug target. Seven essential gene products were predicted to reside in the cell membrane or be released extracellularly, thus serving as putative vaccine candidates. Comparison to essential gene data from five other studies of Vibrio species revealed 13 proteins to be conserved across the studies, while 25 genes were specific to V. anguillarum and not found to be essential in the other Vibrio spp. This study provides new information on the essential genes of Vibrio species and the methodology may be applied to other pathogens to guide the development of new drugs and vaccines, which will assist efforts to counter antibiotic resistance.

How sequence populations persist inside bacterial genomes

Compared to their eukaryotic counterparts, bacterial genomes are small and contain extremely tightly packed genes. Repetitive sequences are rare but not completely absent. One of the most common repeat families is REPINs. REPINs can replicate in the host genome and form populations that persist for millions of years. Here, we model the interactions of these intragenomic sequence populations with the bacterial host. We first confirm well-established results, in the presence and absence of horizontal gene transfer (hgt) sequence populations either expand until they drive the host to extinction or the sequence population gets purged from the genome. We then show that a sequence population can be stably maintained, when each individual sequence provides a benefit that decreases with increasing sequence population size. Maintaining a sequence population of stable size also requires the replication of the sequence population to be costly to the host, otherwise the sequence population size will increase indefinitely. Surprisingly, in regimes with high hgt rates, the benefit conferred by the sequence population does not have to exceed the damage it causes to its host. Our analyses provide a plausible scenario for the persistence of sequence populations in bacterial genomes. We also hypothesize a limited biologically relevant parameter range for the provided benefit, which can be tested in future experiments.

The impact of genetic diversity on gene essentiality within the Escherichia coli species

Bacteria from the same species can differ widely in their gene content. In Escherichia coli, the set of genes shared by all strains, known as the core genome, represents about half the number of genes present in any strain. Although recent advances in bacterial genomics have unravelled genes required for fitness in various experimental conditions, most studies have focused on single model strains. As a result, the impact of the species' genetic diversity on core processes of the bacterial cell remains largely under-investigated. Here, we have developed a CRISPR interference platform for high-throughput gene repression that is compatible with most E. coli isolates and closely related species. We have applied it to assess the importance of ~3,400 nearly ubiquitous genes in three growth conditions in 18 representative E. coli strains spanning most common phylogroups and lifestyles of the species. Our screens revealed extensive variations in gene essentiality between strains and conditions. Investigation of the genetic determinants for these variations highlighted the importance of epistatic interactions with mobile genetic elements. In particular, we have shown how prophage-encoded defence systems against phage infection can trigger the essentiality of persistent genes that are usually non-essential. This study provides broad insights into the evolvability of gene essentiality and argues for the importance of studying various isolates from the same species under diverse conditions.

Ecology, Structure, and Evolution of Shigella Phages

Numerous bacteriophages-viruses of bacteria, also known as phages-have been described for hundreds of bacterial species. The Gram-negative Shigella species are close relatives of Escherichia coli, yet relatively few previously described phages appear to exclusively infect this genus. Recent efforts to isolate Shigella phages have indicated these viruses are surprisingly abundant in the environment and have distinct genomic and structural properties. In addition, at least one model system used for experimental evolution studies has revealed a unique mechanism for developing faster infection cycles. Differences between these bacteriophages and other well-described model systems may mirror differences between their hosts' ecology and defense mechanisms. In this review, we discuss the history of Shigella phages and recent developments in their isolation and characterization and the structural information available for three model systems, Sf6, Sf14, and HRP29; we also provide an overview of potential selective pressures guiding both Shigella phage and host evolution.

A Mini-ISY100 Transposon Delivery System Effective in γ Proteobacteria

Transposons are invaluable biological tools for the genetic manipulation of microorganisms. ISY100 from Synechocystis sp. PCC6803 is a member of the Tc1/mariner/IS630 superfamily, and is characterized by high transposition efficiency and a strong preference for TA target sequences. In this paper, we describe the design and application of a mini-ISY100 suicide vector for the in vivo creation of stable random transposon insertion libraries. The system was successfully applied in seven species belonging to four different orders of γ proteobacteria. In all cases, delivery using conjugation consistently showed the highest transposition efficiency compared to chemical transformation or electroporation. We determined the frequency of transposon insertions in all the species and proved the utility of the system by identifying genes involved in colony coloration in Shewanella oneidensis. The ease and the efficiency of the protocol developed here allow the creation of complete knock-out libraries in an extensive range of host microorganisms in less than a week with no requirement for preparatory modification.

Tn-Core: A Toolbox for Integrating Tn-seq Gene Essentiality Data and Constraint-Based Metabolic Modeling

The design of synthetic cells requires a detailed understanding of the relevance of genes and gene networks underlying complex cellular phenotypes. Transposon-sequencing (Tn-seq) and constraint-based metabolic modeling can be used to probe the core genetic and metabolic networks underlying a biological process. Integrating these highly complementary experimental and in silico approaches has the potential to yield a highly comprehensive understanding of the core networks of a cell. Specifically, it can facilitate the interpretation of Tn-seq data sets and identify gaps in the data that could hinder the engineering of the cellular system, while also providing refined models for the accurate predictions of cellular metabolism. Here, we present Tn-Core, the first easy-to-use computational pipeline specifically designed for integrating Tn-seq data with metabolic modeling, prepared for use by both experimental and computational biologists. Tn-Core is a MATLAB toolbox that contains several custom functions, and it is built upon existing functions within the COBRA Toolbox and the TIGER Toolbox. Tn-Core takes as input a genome-scale metabolic model, Tn-seq data, and optionally RNA-seq data, and returns: (i) a context-specific core metabolic model; (ii) an evaluation of redundancies within core metabolic pathways, and optionally (iii) a refined genome-scale metabolic model. A simple, user-friendly workflow, requiring limited knowledge of metabolic modeling, is provided that allows users to run the analyses and export the data as easy-to-explore files of value to both experimental and computational biologists. We demonstrate the utility of Tn-Core using Sinorhizobium meliloti, Pseudomonas aeruginosa, and Rhodobacter sphaeroides genome-scale metabolic reconstructions as case studies.

Robustness encoded across essential and accessory replicons of the ecologically versatile bacterium Sinorhizobium meliloti

Bacterial genome evolution is characterized by gains, losses, and rearrangements of functional genetic segments. The extent to which large-scale genomic alterations influence genotype-phenotype relationships has not been investigated in a high-throughput manner. In the symbiotic soil bacterium Sinorhizobium meliloti, the genome is composed of a chromosome and two large extrachromosomal replicons (pSymA and pSymB, which together constitute 45% of the genome). Massively parallel transposon insertion sequencing (Tn-seq) was employed to evaluate the contributions of chromosomal genes to growth fitness in both the presence and absence of these extrachromosomal replicons. Ten percent of chromosomal genes from diverse functional categories are shown to genetically interact with pSymA and pSymB. These results demonstrate the pervasive robustness provided by the extrachromosomal replicons, which is further supported by constraint-based metabolic modeling. A comprehensive picture of core S. meliloti metabolism was generated through a Tn-seq-guided in silico metabolic network reconstruction, producing a core network encompassing 726 genes. This integrated approach facilitated functional assignments for previously uncharacterized genes, while also revealing that Tn-seq alone missed over a quarter of wild-type metabolism. This work highlights the many functional dependencies and epistatic relationships that may arise between bacterial replicons and across a genome, while also demonstrating how Tn-seq and metabolic modeling can be used together to yield insights not obtainable by either method alone.

S. meliloti, which has traditionally facilitated ground-breaking insights into symbiotic communication, is also emerging as an excellent model for studying the evolution of functional relationships between bacterial chromosomes and anciently acquired accessory replicons. Multi-replicon genome architecture is present in ~ 10% of presently sequenced bacterial genomes. The S. meliloti genome is composed of three circular replicons, two of which are dispensable even though they encompass nearly half of the protein-coding genes in this organism. The construction of strains lacking these replicons has enabled a straightforward, genome-wide analysis of interactions between the chromosome and the non-essential replicons, revealing extensive functional cooperation between these genomic components. This analysis enabled a substantial refinement of a metabolic network model for S. meliloti. The integration of massively parallel genotype-phenotype screening with in silico metabolic reconstruction has enhanced our understanding of metabolic network structure as it relates to genome evolution in S. meliloti, and exemplifies an approach that may be productively applied to other taxa. The combined experimental and computational approach employed here further provides unique insights into the pervasive genetic interactions that may exist within large bacterial genomes.

The Essential Genome of Escherichia coli K-12

Transposon-directed insertion site sequencing (TraDIS) is a high-throughput method coupling transposon mutagenesis with short-fragment DNA sequencing. It is commonly used to identify essential genes. Single gene deletion libraries are considered the gold standard for identifying essential genes. Currently, the TraDIS method has not been benchmarked against such libraries, and therefore, it remains unclear whether the two methodologies are comparable. To address this, a high-density transposon library was constructed in Escherichia coli K-12. Essential genes predicted from sequencing of this library were compared to existing essential gene databases. To decrease false-positive identification of essential genes, statistical data analysis included corrections for both gene length and genome length. Through this analysis, new essential genes and genes previously incorrectly designated essential were identified. We show that manual analysis of TraDIS data reveals novel features that would not have been detected by statistical analysis alone. Examples include short essential regions within genes, orientation-dependent effects, and fine-resolution identification of genome and protein features. Recognition of these insertion profiles in transposon mutagenesis data sets will assist genome annotation of less well characterized genomes and provides new insights into bacterial physiology and biochemistry.

Incentives to define lists of genes that are essential for bacterial survival include the identification of potential targets for antibacterial drug development, genes required for rapid growth for exploitation in biotechnology, and discovery of new biochemical pathways. To identify essential genes in Escherichia coli, we constructed a transposon mutant library of unprecedented density. Initial automated analysis of the resulting data revealed many discrepancies compared to the literature. We now report more extensive statistical analysis supported by both literature searches and detailed inspection of high-density TraDIS sequencing data for each putative essential gene for the E. coli model laboratory organism. This paper is important because it provides a better understanding of the essential genes of E. coli, reveals the limitations of relying on automated analysis alone, and provides a new standard for the analysis of TraDIS data.

A Comprehensive Overview of Online Resources to Identify and Predict Bacterial Essential Genes

Genes critical for the survival or reproduction of an organism in certain circumstances are classified as essential genes. Essential genes play a significant role in deciphering the survival mechanism of life. They may be greatly applied to pharmaceutics and synthetic biology. The continuous progress of experimental method for essential gene identification has accelerated the accumulation of gene essentiality data which facilitates the study of essential genes in silico. In this article, we present some available online resources related to gene essentiality, including bioinformatic software tools for transposon sequencing (Tn-seq) analysis, essential gene databases and online services to predict bacterial essential genes. We review several computational approaches that have been used to predict essential genes, and summarize the features used for gene essentiality prediction. In addition, we evaluate the available online bacterial essential gene prediction servers based on the experimentally validated essential gene sets of 30 bacteria from DEG. This article is intended to be a quick reference guide for the microbiologists interested in the essential genes.

Creation of a Dense Transposon Insertion Library Using Bacterial Conjugation in Enterobacterial Strains Such As Escherichia Coli or Shigella flexneri

Transposon mutagenesis is a method that allows gene disruption via the random genomic insertion of a piece of DNA called a transposon. The protocol below outlines a method for high efficiency transfer between bacterial strains of a plasmid harboring a transposon containing a kanamycin resistance marker. The plasmid-borne transposase is encoded by a variant tnp gene that inserts the transposon into the genome of the recipient strain with very low insertional bias. This method thus allows the creation of large mutant libraries in which transposons have been inserted into unique genomic positions in a recipient strain of either Escherichia coli or Shigella flexneri bacteria. By using bacterial conjugation, as opposed to other methods such as electroporation or chemical transformation, large libraries with hundreds of thousands of unique clones can be created. This yields high-density insertion libraries, with insertions occurring as frequently as every 4-6 base pairs in non-essential genes. This method is superior to other methods as it allows for an inexpensive, easy to use, and high efficiency method for the creation of a dense transposon insertion library. The transposon library can be used in downstream applications such as transposon sequencing (Tn-Seq), to infer genetic interaction networks, or more simply, in mutational (forward genetic) screens.