A vast collection of microbial genes that are toxic to bacteria

Computational evaluation of cellular metabolic costs successfully predicts genes whose expression is deleterious

Gene suppression and overexpression are both fundamental tools in linking genotype to phenotype in model organisms. Computational methods have proven invaluable in studying and predicting the deleterious effects of gene deletions, and yet parallel computational methods for overexpression are still lacking. Here, we present Expression-Dependent Gene Effects (EDGE), an in silico method that can predict the deleterious effects resulting from overexpression of either native or foreign metabolic genes. We first test and validate EDGE's predictive power in bacteria through a combination of small-scale growth experiments that we performed and analysis of extant large-scale datasets. Second, a broad cross-species analysis, ranging from microorganisms to multiple plant and human tissues, shows that genes that EDGE predicts to be deleterious when overexpressed are indeed typically down-regulated. This reflects a universal selection force keeping the expression of potentially deleterious genes in check. Third, EDGE-based analysis shows that cancer genetic reprogramming specifically suppresses genes whose overexpression impedes proliferation. The magnitude of this suppression is large enough to enable an almost perfect distinction between normal and cancerous tissues based solely on EDGE results. We expect EDGE to advance our understanding of human pathologies associated with up-regulation of particular transcripts and to facilitate the utilization of gene overexpression in metabolic engineering.

Comprehensive analysis of the HEPN superfamily: identification of novel roles in intra-genomic conflicts, defense, pathogenesis and RNA processing

Background

The major role of enzymatic toxins that target nucleic acids in biological conflicts at all levels has become increasingly apparent thanks in large part to the advances of comparative genomics. Typically, toxins evolve rapidly hampering the identification of these proteins by sequence analysis. Here we analyze an unexpectedly widespread superfamily of toxin domains most of which possess RNase activity.

Results

The HEPN superfamily is comprised of all α-helical domains that were first identified as being associated with DNA polymerase β-type nucleotidyltransferases in prokaryotes and animal Sacsin proteins. Using sensitive sequence and structure comparison methods, we vastly extend the HEPN superfamily by identifying numerous novel families and by detecting diverged HEPN domains in several known protein families. The new HEPN families include the RNase LS and LsoA catalytic domains, KEN domains (e.g. RNaseL and Ire1) and the RNase domains of RloC and PrrC. The majority of HEPN domains contain conserved motifs that constitute a metal-independent endoRNase active site. Some HEPN domains lacking this motif probably function as non-catalytic RNA-binding domains, such as in the case of the mannitol repressor MtlR. Our analysis shows that HEPN domains function as toxins that are shared by numerous systems implicated in intra-genomic, inter-genomic and intra-organismal conflicts across the three domains of cellular life. In prokaryotes HEPN domains are essential components of numerous toxin-antitoxin (TA) and abortive infection (Abi) systems and in addition are tightly associated with many restriction-modification (R-M) and CRISPR-Cas systems, and occasionally with other defense systems such as Pgl and Ter. We present evidence of multiple modes of action of HEPN domains in these systems, which include direct attack on viral RNAs (e.g. LsoA and RNase LS) in conjunction with other RNase domains (e.g. a novel RNase H fold domain, NamA), suicidal or dormancy-inducing attack on self RNAs (RM systems and possibly CRISPR-Cas systems), and suicidal attack coupled with direct interaction with phage components (Abi systems). These findings are compatible with the hypothesis on coupling of pathogen-targeting (immunity) and self-directed (programmed cell death and dormancy induction) responses in the evolution of robust antiviral strategies. We propose that altruistic cell suicide mediated by HEPN domains and other functionally similar RNases was essential for the evolution of kin and group selection and cell cooperation. HEPN domains were repeatedly acquired by eukaryotes and incorporated into several core functions such as endonucleolytic processing of the 5.8S-25S/28S rRNA precursor (Las1), a novel ER membrane-associated RNA degradation system (C6orf70), sensing of unprocessed transcripts at the nuclear periphery (Swt1). Multiple lines of evidence suggest that, similar to prokaryotes, HEPN proteins were recruited to antiviral, antitransposon, apoptotic systems or RNA-level response to unfolded proteins (Sacsin and KEN domains) in several groups of eukaryotes.

Conclusions

Extensive sequence and structure comparisons reveal unexpectedly broad presence of the HEPN domain in an enormous variety of defense and stress response systems across the tree of life. In addition, HEPN domains have been recruited to perform essential functions, in particular in eukaryotic rRNA processing. These findings are expected to stimulate experiments that could shed light on diverse cellular processes across the three domains of life.

Reviewers

This article was reviewed by Martijn Huynen, Igor Zhulin and Nick Grishin

Characterizing and measuring bias in sequence data

Background

DNA sequencing technologies deviate from the ideal uniform distribution of reads. These biases impair scientific and medical applications. Accordingly, we have developed computational methods for discovering, describing and measuring bias.

Results

We applied these methods to the Illumina, Ion Torrent, Pacific Biosciences and Complete Genomics sequencing platforms, using data from human and from a set of microbes with diverse base compositions. As in previous work, library construction conditions significantly influence sequencing bias. Pacific Biosciences coverage levels are the least biased, followed by Illumina, although all technologies exhibit error-rate biases in high- and low-GC regions and at long homopolymer runs. The GC-rich regions prone to low coverage include a number of human promoters, so we therefore catalog 1,000 that were exceptionally resistant to sequencing. Our results indicate that combining data from two technologies can reduce coverage bias if the biases in the component technologies are complementary and of similar magnitude. Analysis of Illumina data representing 120-fold coverage of a well-studied human sample reveals that 0.20% of the autosomal genome was covered at less than 10% of the genome-wide average. Excluding locations that were similar to known bias motifs or likely due to sample-reference variations left only 0.045% of the autosomal genome with unexplained poor coverage.

Conclusions

The assays presented in this paper provide a comprehensive view of sequencing bias, which can be used to drive laboratory improvements and to monitor production processes. Development guided by these assays should result in improved genome assemblies and better coverage of biologically important loci.

Genomic and proteomic analyses of the terminally redundant genome of the Pseudomonas aeruginosa phage PaP1: establishment of genus PaP1-like phages

We isolated and characterized a new Pseudomonas aeruginosa myovirus named PaP1. The morphology of this phage was visualized by electron microscopy and its genome sequence and ends were determined. Finally, genomic and proteomic analyses were performed. PaP1 has an icosahedral head with an apex diameter of 68–70 nm and a contractile tail with a length of 138–140 nm. The PaP1 genome is a linear dsDNA molecule containing 91,715 base pairs (bp) with a G+C content of 49.36% and 12 tRNA genes. A strategy to identify the genome ends of PaP1 was designed. The genome has a 1190 bp terminal redundancy. PaP1 has 157 open reading frames (ORFs). Of these, 143 proteins are homologs of known proteins, but only 38 could be functionally identified. Sodium dodecyl sulfate-polyacrylamide gel electrophoresis and high-performance liquid chromatography-mass spectrometry allowed identification of 12 ORFs as structural protein coding genes within the PaP1 genome. Comparative genomic analysis indicated that the Pseudomonas aeruginosa phage PaP1, JG004, PAK_P1 and vB_PaeM_C2-10_Ab1 share great similarity. Besides their similar biological characteristics, the phages contain 123 core genes and have very close phylogenetic relationships, which distinguish them from other known phage genera. We therefore propose that these four phages be classified as PaP1-like phages, a new phage genus of Myoviridae that infects Pseudomonas aeruginosa.

YgfX (CptA) is a multimeric membrane protein that interacts with the succinate dehydrogenase assembly factor SdhE (YgfY)

Serratia sp. strain ATCC 39006 produces the red-pigmented antibiotic prodigiosin. Prodigiosin biosynthesis is regulated by a complex hierarchy that includes the uncharacterized protein YgfX (DUF1434). The ygfX gene is co-transcribed with sdhE, an FAD assembly factor essential for the flavinylation and activation of the SdhA subunit of succinate dehydrogenase (SDH), a central enzyme in the tricarboxylic acid cycle and electron transport chain. The sdhEygfX operon is highly conserved within the Enterobacteriaceae, suggesting that SdhE and YgfX function together. We performed an extensive mutagenesis to gain molecular insights into the uncharacterized protein YgfX, and have investigated the relationship between YgfX and SdhE. YgfX localized to the membrane, interacted with itself, forming dimers or larger multimers, and interacted with SdhE. The transmembrane helices of YgfX were critical for protein function and the formation of YgfX multimers. Site-directed mutagenesis of residues conserved in DUF1434 proteins revealed a periplasmic tryptophan and a cytoplasmic aspartate that were crucial for YgfX activity. Both of these amino acids were required for the formation of YgfX multimers and interactions with SdhE but not membrane localization. Multiple cell division proteins were identified as putative interaction partners of YgfX and overexpression of YgfX had effects on cell morphology. These findings represent an important step in understanding the function of DUF1434 proteins. In contrast to a recent report, we found no evidence that YgfX and SdhE form a toxin-antitoxin system. In summary, YgfX functions as a multimeric membrane-bound protein that interacts with SdhE, an important FAD assembly factor that controls SDH activity.

Development of a Fur-dependent and tightly regulated expression system in Escherichia coli for toxic protein synthesis

Background

There is a continuous demanding for tightly regulated prokaryotic expression systems, which allow functional synthesis of toxic proteins in Escherichia coli for bioscience or biotechnology application. However, most of the current promoter options either are tightly repressed only with low protein production levels, or produce substantial protein but lacking of the necessary repression to avoid mutations initiated by leaky expression in the absence of inducer. The aim of this study was to develop a tightly regulated, relatively high-efficient expression vector in E. coli based on the principle of iron uptake system.

Results

By using GFP as reporter, P_fhuA with the highest relative fluorescence units, but leaky expression, was screened from 23 iron-regulated promoter candidates. P_fhuA was repressed by ferric uptake regulator (Fur)-Fe²⁺ complex binding to Fur box locating at the promoter sequence. Otherwise, P_fhuA was activated without Fur-Fe²⁺ binding in the absence of iron. In order to improve the tightness of P_fhuA regulation for toxic gene expression, Fur box in promoter sequence and fur expression were refined through five different approaches. Eventually, through substituting E. coli consensus Fur box for original one of P_fhuA, the induction ratio of modified P_fhuA (named P_fhuA1) was improved from 3 to 101. Under the control of P_fhuA1, strong toxic gene E was successfully expressed in high, middle, low copy-number vectors, and other two toxic proteins, Gef and MazF were functionally synthesized without E. coli death before induction.

Conclusions

The features of easy control, tight regulation and relatively high efficiency were combined in the newly engineered P_fhuA1. Under this promoter, the toxic genes E, gef and mazF were functionally expressed in E. coli induced by iron chelator in a tightly controllable way. This study provides a tightly regulated expression system that might enable the stable cloning, and functional synthesis of toxic proteins for their function study, bacterial programmed cell death in biological containment system and bacterial vector vaccine development.

Discovery of functional toxin/antitoxin systems in bacteria by shotgun cloning

Toxin-antitoxin (TA) modules, composed of a toxic protein and a counteracting antitoxin, play important roles in bacterial physiology. We examined the experimental insertion of 1.5 million genes from 388 microbial genomes into an Escherichia coli host using more than 8.5 million random clones. This revealed hundreds of genes (toxins) that could only be cloned when the neighboring gene (antitoxin) was present on the same clone. Clustering of these genes revealed TA families widespread in bacterial genomes, some of which deviate from the classical characteristics previously described for such modules. Introduction of these genes into E. coli validated that the toxin toxicity is mitigated by the antitoxin. Infection experiments with T7 phage showed that two of the new modules can provide resistance against phage. Moreover, our experiments revealed an "antidefense" protein in phage T7 that neutralizes phage resistance. Our results expose active fronts in the arms race between bacteria and phage.

Comparative genomics of defense systems in archaea and bacteria

Our knowledge of prokaryotic defense systems has vastly expanded as the result of comparative genomic analysis, followed by experimental validation. This expansion is both quantitative, including the discovery of diverse new examples of known types of defense systems, such as restriction-modification or toxin-antitoxin systems, and qualitative, including the discovery of fundamentally new defense mechanisms, such as the CRISPR-Cas immunity system. Large-scale statistical analysis reveals that the distribution of different defense systems in bacterial and archaeal taxa is non-uniform, with four groups of organisms distinguishable with respect to the overall abundance and the balance between specific types of defense systems. The genes encoding defense system components in bacterial and archaea typically cluster in defense islands. In addition to genes encoding known defense systems, these islands contain numerous uncharacterized genes, which are candidates for new types of defense systems. The tight association of the genes encoding immunity systems and dormancy- or cell death-inducing defense systems in prokaryotic genomes suggests that these two major types of defense are functionally coupled, providing for effective protection at the population level.

Production of bulk chemicals via novel metabolic pathways in microorganisms

Metabolic engineering has been playing important roles in developing high performance microorganisms capable of producing various chemicals and materials from renewable biomass in a sustainable manner. Synthetic and systems biology are also contributing significantly to the creation of novel pathways and the whole cell-wide optimization of metabolic performance, respectively. In order to expand the spectrum of chemicals that can be produced biotechnologically, it is necessary to broaden the metabolic capacities of microorganisms. Expanding the metabolic pathways for biosynthesizing the target chemicals requires not only the enumeration of a series of known enzymes, but also the identification of biochemical gaps whose corresponding enzymes might not actually exist in nature; this issue is the focus of this paper. First, pathway prediction tools, effectively combining reactions that lead to the production of a target chemical, are analyzed in terms of logics representing chemical information, and designing and ranking the proposed metabolic pathways. Then, several approaches for potentially filling in the gaps of the novel metabolic pathway are suggested along with relevant examples, including the use of promiscuous enzymes that flexibly utilize different substrates, design of novel enzymes for non-natural reactions, and exploration of hypothetical proteins. Finally, strain optimization by systems metabolic engineering in the context of novel metabolic pathways constructed is briefly described. It is hoped that this review paper will provide logical ways of efficiently utilizing 'big' biological data to design and develop novel metabolic pathways for the production of various bulk chemicals that are currently produced from fossil resources.

Assessing the effectiveness of functional genetic screens for the identification of bioactive metabolites

A common limitation for the identification of novel activities from functional (meta) genomic screens is the low number of active clones detected relative to the number of clones screened. Here we demonstrate that constructing libraries with strains known to produce bioactives can greatly enhance the screening efficiency, by increasing the “hit-rate” and unmasking multiple activities from the same bacterial source.

Live virus-free or die: coupling of antivirus immunity and programmed suicide or dormancy in prokaryotes

Background

The virus-host arms race is a major theater for evolutionary innovation. Archaea and bacteria have evolved diverse, elaborate antivirus defense systems that function on two general principles: i) immune systems that discriminate self DNA from nonself DNA and specifically destroy the foreign, in particular viral, genomes, whereas the host genome is protected, or ii) programmed cell suicide or dormancy induced by infection.

Presentation of the hypothesis

Almost all genomic loci encoding immunity systems such as CRISPR-Cas, restriction-modification and DNA phosphorothioation also encompass suicide genes, in particular those encoding known and predicted toxin nucleases, which do not appear to be directly involved in immunity. In contrast, the immunity systems do not appear to encode antitoxins found in typical toxin-antitoxin systems. This raises the possibility that components of the immunity system themselves act as reversible inhibitors of the associated toxin proteins or domains as has been demonstrated for the Escherichia coli anticodon nuclease PrrC that interacts with the PrrI restriction-modification system. We hypothesize that coupling of diverse immunity and suicide/dormancy systems in prokaryotes evolved under selective pressure to provide robustness to the antivirus response. We further propose that the involvement of suicide/dormancy systems in the coupled antivirus response could take two distinct forms:

1) induction of a dormancy-like state in the infected cell to ‘buy time’ for activation of adaptive immunity; 2) suicide or dormancy as the final recourse to prevent viral spread triggered by the failure of immunity.

Testing the hypothesis

This hypothesis entails many experimentally testable predictions. Specifically, we predict that Cas2 protein present in all cas operons is a mRNA-cleaving nuclease (interferase) that might be activated at an early stage of virus infection to enable incorporation of virus-specific spacers into the CRISPR locus or to trigger cell suicide when the immune function of CRISPR-Cas systems fails. Similarly, toxin-like activity is predicted for components of numerous other defense loci.

Implications of the hypothesis

The hypothesis implies that antivirus response in prokaryotes involves key decision-making steps at which the cell chooses the path to follow by sensing the course of virus infection.

Reviewers

This article was reviewed by Arcady Mushegian, Etienne Joly and Nick Grishin. For complete reviews, go to the Reviewers’ reports section.

Successes and failures in modular genetic engineering

Synthetic biology relies on engineering concepts such as abstraction, standardization, and decoupling to develop systems that address environmental, clinical, and industrial needs. Recent advances in applying modular design to system development have enabled creation of increasingly complex systems. However, several challenges to module and system development remain, including syntactic errors, semantic errors, parameter mismatches, contextual sensitivity, noise and evolution, and load and stress. To combat these challenges, researchers should develop a framework for describing and reasoning about biological information, design systems with modularity in mind, and investigate how to predictively describe the diverse sources and consequences of metabolic load and stress.

PanDaTox: a tool for accelerated metabolic engineering

Metabolic engineering is often facilitated by cloning of genes encoding enzymes from various heterologous organisms into E. coli. Such engineering efforts are frequently hampered by foreign genes that are toxic to the E. coli host. We have developed PanDaTox (www.weizmann.ac.il/pandatox), a web-based resource that provides experimental toxicity information for more than 1.5 million genes from hundreds of different microbial genomes. The toxicity predictions, which were extensively experimentally verified, are based on serial cloning of genes into E. coli as part of the Sanger whole genome shotgun sequencing process. PanDaTox can accelerate metabolic engineering projects by allowing researchers to exclude toxic genes from the engineering plan and verify the clonability of selected genes before the actual metabolic engineering experiments are conducted.