Essential genes in Bacillus subtilis: a re-evaluation after ten years

Discovering genetic modulators of the protein homeostasis system through multilevel analysis

Every protein progresses through a natural lifecycle from birth to maturation to death; this process is coordinated by the protein homeostasis system. Environmental or physiological conditions trigger pathways that maintain the homeostasis of the proteome. An open question is how these pathways are modulated to respond to the many stresses that an organism encounters during its lifetime. To address this question, we tested how the fitness landscape changes in response to environmental and genetic perturbations using directed and massively parallel transposon mutagenesis in Caulobacter crescentus. We developed a general computational pipeline for the analysis of gene-by-environment interactions in transposon mutagenesis experiments. This pipeline uses a combination of general linear models, statistical knockoffs, and a nonparametric Bayesian statistical model to identify essential genetic network components that are shared across environmental perturbations. This analysis allows us to quantify the similarity of proteotoxic environmental perturbations from the perspective of the fitness landscape. We find that essential genes vary more by genetic background than by environmental conditions, with limited overlap among mutant strains targeting different facets of the protein homeostasis system. We also identified 146 unique fitness determinants across different strains, with 19 genes common to at least two strains, showing varying resilience to proteotoxic stresses. Experiments exposing cells to a combination of genetic perturbations and dual environmental stressors show that perturbations that are quantitatively dissimilar from the perspective of the fitness landscape are likely to have a synergistic effect on the growth defect.

Signal recognition particle RNA is critical for genetic competence and virulence of Streptococcus pneumoniae

Streptococcus pneumoniae (pneumococcus) causes a wide range of important human infectious diseases, including pneumonia, pneumonia-derived sepsis, otitis media, and meningitis. Pneumococcus produces numerous secreted proteins that are critical for normal physiology and pathogenesis. The membrane targeting and translocation of these secreted proteins are partly mediated by the signal recognition particle (SRP) complex, which consists of 4.5S small cytoplasmic RNA (ScRNA), and the Ffh, and FtsY proteins. Here, we report that pneumococcal ∆scRNA, ∆ffh, and ∆ftsY mutants were significantly impaired in competence induction, competence pili production, exogenous DNA uptake, and genetic transformation. Also, the ∆scRNA mutant was significantly attenuated in the mouse models of bacteremia and pneumonia. Interestingly, unlike the ∆scRNA, both ∆ffh and ∆ftsY mutants had growth defects on Todd-Hewitt Agar, which were alleviated by the provision of free amino acids or serum. Differences in nutritional requirements between ∆ffh and ∆ftsY vs ∆scRNA suggest that Ffh and FtsY may be partially functional in the absence of ScRNA. Finally, the insertase YidC2, which could functionally rescue some SRP mutations in other streptococcal species, was not essential for pneumococcal genetic transformation. Collectively, these results indicate that ScRNA is crucial for the successful development of genetic competence and virulence in pneumococcus.

IMPORTANCE

Streptococcus pneumoniae (pneumococcus) causes multiple important infectious diseases in humans. The signal recognition particle (SRP) complex, which comprised 4.5S small cytoplasmic RNA (ScRNA), and the Ffh and FtsY proteins, mediates membrane targeting and translocation of secreted proteins in all organisms. However, the role of SRP and ScRNA has not been characterized during the induction of the competence system for genetic transformation and virulence in pneumococcus. By using a combination of genetic, biochemical, proteomic, and imaging approaches, we demonstrated that the SRP complex plays a significant role in membrane targeting of competence system-regulated effectors important for genetic transformation, virulence during bacteremia and pneumonia infections, and nutritional acquisition.

Reduction-to-synthesis: the dominant approach to genome-scale synthetic biology

Advances in systems and synthetic biology have propelled the construction of reduced bacterial genomes. Genome reduction was initially focused on exploring properties of minimal genomes, but more recently it has been deployed as an engineering strategy to enhance strain performance. This review provides the latest updates on reduced genomes, focusing on dual-track approaches of top-down reduction and bottom-up synthesis for their construction. Using cases from studies that are based on established industrial workhorse strains, we discuss the construction of a series of synthetic phenotypes that are candidates for biotechnological applications. Finally, we address the possible uses of reduced genomes for biotechnological applications and the needed future research directions that may ultimately lead to the total synthesis of rationally designed genomes.

The gap gene of Rhizobium etli is required for both free life and symbiosis with common beans

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH or Gap) is a ubiquitous enzyme essential for carbon and energy metabolism in most organisms. Despite its primary role in sugar metabolism, GAPDH is recognized for its involvement in diverse cellular processes, being considered a paradigm among multifunctional/moonlighting proteins. Besides its canonical cytoplasmic location, GAPDH has been detected on cell surfaces or as a secreted protein in prokaryotes, yet little is known about its possible roles in plant symbiotic bacteria. Here we report that Rhizobium etli, a nitrogen-fixing symbiont of common beans, carries a single gap gene responsible for both GAPDH glycolytic and gluconeogenic activities. An active Gap protein is required throughout all stages of the symbiosis between R. etli and its host plant Phaseolus vulgaris. Both glycolytic and gluconeogenic Gap metabolic activities likely contribute to bacterial fitness during early and intermediate stages of the interaction, whereas GAPDH gluconeogenic activity seems critical for nodule invasion and nitrogen fixation. Although the R. etli Gap protein is secreted in a c-di-GMP related manner, no involvement of the R. etli gap gene in c-di-GMP related phenotypes, such as flocculation, biofilm formation or EPS production, was observed. Notably, the R. etli gap gene fully complemented a double gap1/gap2 mutant of Pseudomonas syringae for free life growth, albeit only partially in planta, suggesting potential specific roles for each type of Gap protein. Nevertheless, further research is required to unravel additional functions of the R. etli Gap protein beyond its essential metabolic roles.

Lethal metabolism of Candida albicans respiratory mutants

The destructive impact of fungi in agriculture and animal and human health, coincident with increases in antifungal resistance, underscores the need for new and alternative drug targets to counteract these trends. Cellular metabolism relies on many intermediates with intrinsic toxicity and promiscuous enzymatic activity generates others. Fuller knowledge of these toxic entities and their generation may offer opportunities of antifungal development. From this perspective our observation of media-conditional lethal metabolism in respiratory mutants of the opportunistic fungal pathogen Candida albicans was of interest. C. albicans mutants defective in NADH:ubiquinone oxidoreductase (Complex I of the electron transport chain) exhibit normal growth in synthetic complete medium. In YPD medium, however, the mutants grow normally until early stationary phase whereupon a dramatic loss of viability occurs. Upwards of 90% of cells die over the subsequent four to six hours with a loss of membrane integrity. The extent of cell death was proportional to the amount of BactoPeptone, and to a lesser extent, the amount of yeast extract. YPD medium conditioned by growth of the mutant was toxic to wild-type cells indicating mutant metabolism established a toxic milieu in the media. Conditioned media contained a volatile component that contributed to toxicity, but only in the presence of a component of BactoPeptone. Fractionation experiments revealed purine nucleosides or bases as the synergistic component. GC-mass spectrometry analysis revealed acetal (1,1-diethoxyethane) as the active volatile. This previously unreported and lethal synergistic interaction of acetal and purines suggests a hitherto unrecognized toxic metabolism potentially exploitable in the search for antifungal targets.

Recent advances in the characterization of essential genes and development of a database of essential genes

Over the past few decades, there has been a significant interest in the study of essential genes, which are crucial for the survival of an organism under specific environmental conditions and thus have practical applications in the fields of synthetic biology and medicine. An increasing amount of experimental data on essential genes has been obtained with the continuous development of technological methods. Meanwhile, various computational prediction methods, related databases and web servers have emerged accordingly. To facilitate the study of essential genes, we have established a database of essential genes (DEG), which has become popular with continuous updates to facilitate essential gene feature analysis and prediction, drug and vaccine development, as well as artificial genome design and construction. In this article, we summarized the studies of essential genes, overviewed the relevant databases, and discussed their practical applications. Furthermore, we provided an overview of the main applications of DEG and conducted comprehensive analyses based on its latest version. However, it should be noted that the essential gene is a dynamic concept instead of a binary one, which presents both opportunities and challenges for their future development.

Structural assembly of the bacterial essential interactome

The study of protein interactions in living organisms is fundamental for understanding biological processes and central metabolic pathways. Yet, our knowledge of the bacterial interactome remains limited. Here, we combined gene deletion mutant analysis with deep-learning protein folding using AlphaFold2 to predict the core bacterial essential interactome. We predicted and modeled 1402 interactions between essential proteins in bacteria and generated 146 high-accuracy models. Our analysis reveals previously unknown details about the assembly mechanisms of these complexes, highlighting the importance of specific structural features in their stability and function. Our work provides a framework for predicting the essential interactomes of bacteria and highlight the potential of deep-learning algorithms in advancing our understanding of the complex biology of living organisms. Also, the results presented here offer a promising approach to identify novel antibiotic targets.

Two glyceraldehyde-3-phosphate dehydrogenases with distinctive roles in Pseudomonas syringae pv. tomato DC3000

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH or Gap) is a ubiquitously distributed enzyme that plays an essential role in the glycolytic and gluconeogenic pathways. However, additional roles have been described unrelated to its enzymatic function in diverse organisms, often linked to its presence in the cell surface or as a secreted protein. Despite being a paradigm among multifunctional/moonlighting proteins, little is known about its possible roles in phytopathogenic bacteria. In the present work we have studied three putative gap paralogous genes identified in the genome of Pseudomonas syringae pv. tomato (Pto) DC3000, an important model in molecular plant pathology, with the aim of determining their physiological and possible non-canonical roles in this bacterium and in the plant infection process. We have established that the Gap1 protein has a predominantly glycolytic activity, whereas the NADPH-dependent Gap2 main activity is gluconeogenic. The third paralogue lacks GAPDH activity in Pto but is indispensable for vitamin B6 metabolism and displays erythrose-4-phosphate dehydrogenase activity, thus referred as epd. Both Gap enzymes exhibit distinct functional characteristics depending on the bacterium physiological state, with Gap1 presenting a substantial role in motility, biosurfactant production and biofilm formation. On the other hand, solely Gap2 appears to be essential for growth on tomato plant. Furthermore, Gap1 and Gap2 present a distinctive transcriptional regulation and both have been identified exported outside the cells with different definite media compositions. This serves as compelling evidence of additional roles beyond their central metabolic functions.

Signatures of kin selection in a natural population of the bacteria Bacillus subtilis

Laboratory experiments have suggested that bacteria perform a range of cooperative behaviors, which are favored because they are directed toward relatives (kin selection). However, there is a lack of evidence for cooperation and kin selection in natural bacterial populations. Molecular population genetics offers a promising method to study natural populations because the theory predicts that kin selection will lead to relaxed selection, which will result in increased polymorphism and divergence at cooperative genes. Examining a natural population of Bacillus subtilis, we found consistent evidence that putatively cooperative traits have higher polymorphism and greater divergence than putatively private traits expressed at the same rate. In addition, we were able to eliminate alternative explanations for these patterns and found more deleterious mutations in genes controlling putatively cooperative traits. Overall, our results suggest that cooperation is favored by kin selection, with an average relatedness of r = .79 between interacting individuals.

Bacillus subtilis, a Swiss Army Knife in Science and Biotechnology

Next to Escherichia coli, Bacillus subtilis is the most studied and best understood organism that also serves as a model for many important pathogens. Due to its ability to form heat-resistant spores that can germinate even after very long periods of time, B. subtilis has attracted much scientific interest. Another feature of B. subtilis is its genetic competence, a developmental state in which B. subtilis actively takes up exogenous DNA. This makes B. subtilis amenable to genetic manipulation and investigation. The bacterium was one of the first with a fully sequenced genome, and it has been subject to a wide variety of genome- and proteome-wide studies that give important insights into many aspects of the biology of B. subtilis. Due to its ability to secrete large amounts of proteins and to produce a wide range of commercially interesting compounds, B. subtilis has become a major workhorse in biotechnology. Here, we review the development of important aspects of the research on B. subtilis with a specific focus on its cell biology and biotechnological and practical applications from vitamin production to concrete healing. The intriguing complexity of the developmental programs of B. subtilis, paired with the availability of sophisticated tools for genetic manipulation, positions it at the leading edge for discovering new biological concepts and deepening our understanding of the organization of bacterial cells.

Meddling with Metal Sensors: Fur-Family Proteins as Signaling Hubs

The ferric uptake regulator (Fur) protein is the founding member of the FUR superfamily of metalloregulatory proteins that control metal homeostasis in bacteria. FUR proteins regulate metal homeostasis in response to the binding of iron (Fur), zinc (Zur), manganese (Mur), or nickel (Nur). FUR family proteins are generally dimers in solution, but the DNA-bound complex can involve a single dimer, a dimer-of-dimers, or an extended array of bound protein. Elevated FUR levels due to changes in cell physiology increase DNA occupancy and may also kinetically facilitate protein dissociation. Interactions between FUR proteins and other regulators are commonplace, often including cooperative and competitive DNA-binding interactions within the regulatory region. Further, there are many emerging examples of allosteric regulators that interact directly with FUR family proteins. Here, we focus on newly uncovered examples of allosteric regulation by diverse Fur antagonists (Escherichia coli YdiV/SlyD, Salmonella enterica EIIANtr, Vibrio parahaemolyticus FcrX, Acinetobacter baumannii BlsA, Bacillus subtilis YlaN, and Pseudomonas aeruginosa PacT) as well as one Zur antagonist (Mycobacterium bovis CmtR). Small molecules and metal complexes may also serve as regulatory ligands, with examples including heme binding to Bradyrhizobium japonicum Irr and 2-oxoglutarate binding to Anabaena FurA. How these protein-protein and protein-ligand interactions act in conjunction with regulatory metal ions to facilitate signal integration is an active area of investigation.

Trimming the genomic fat: minimising and re-functionalising genomes using synthetic biology

Naturally evolved organisms typically have large genomes that enable their survival and growth under various conditions. However, the complexity of genomes often precludes our complete understanding of them, and limits the success of biotechnological designs. In contrast, minimal genomes have reduced complexity and therefore improved engineerability, increased biosynthetic capacity through the removal of unnecessary genetic elements, and less recalcitrance to complete characterisation. Here, we review the past and current genome minimisation and re-functionalisation efforts, with an emphasis on the latest advances facilitated by synthetic genomics, and provide a critical appraisal of their potential for industrial applications.

How To Deal with Toxic Amino Acids: the Bipartite AzlCD Complex Exports Histidine in Bacillus subtilis

The Gram-positive model bacterium Bacillus subtilis can use several amino acids as sources of carbon and nitrogen. However, some amino acids inhibit the growth of this bacterium. This amino acid toxicity is often enhanced in strains lacking the second messenger cyclic dimeric adenosine 3',5'-monophosphate (c-di-AMP). We observed that the presence of histidine is also toxic for a B. subtilis strain that lacks all three c-di-AMP synthesizing enzymes. However, suppressor mutants emerged, and whole-genome sequencing revealed mutations in the azlB gene that encode the repressor of the azl operon. This operon encodes an exporter and an importer for branched-chain amino acids. The suppressor mutations result in an overexpression of the azl operon. Deletion of the azlCD genes encoding the branched-chain amino acid exporter restored the toxicity of histidine, indicating that this exporter is required for histidine export and for resistance to otherwise toxic levels of the amino acid. The higher abundance of the amino acid exporter AzlCD increased the extracellular concentration of histidine, thus confirming the new function of AzlCD as a histidine exporter. Unexpectedly, the AzlB-mediated repression of the operon remains active even in the presence of amino acids, suggesting that the expression of the azl operon requires the mutational inactivation of AzlB. IMPORTANCE Amino acids are building blocks for protein biosynthesis in each living cell. However, due to their reactivity and the similarity between several amino acids, they may also be involved in harmful reactions or in noncognate interactions and thus may be toxic. Bacillus subtilis can deal with otherwise toxic histidine by overexpressing the bipartite amino acid exporter AzlCD. Although encoded in an operon that also contains a gene for an amino acid importer, the corresponding genes are not expressed, irrespective of the availability of amino acids in the medium. This suggests that the azl operon is a last resort by which to deal with histidine stress that can be expressed due to the mutational inactivation of the cognate repressor AzlB.

Listeria monocytogenes genes supporting growth under standard laboratory cultivation conditions and during macrophage infection

The Gram-positive bacterium Listeria monocytogenes occurs widespread in the environment and infects humans when ingested along with contaminated food. Such infections are particularly dangerous for risk group patients, for whom they represent a life-threatening disease. To invent novel strategies to control contamination and disease, it is important to identify those cellular processes that maintain pathogen growth inside and outside the host. Here, we have applied transposon insertion sequencing (Tn-Seq) to L. monocytogenes for the identification of such processes on a genome-wide scale. Our approach identified 394 open reading frames that are required for growth under standard laboratory conditions and 42 further genes, which become necessary during intracellular growth in macrophages. Most of these genes encode components of the translation machinery and act in chromosome-related processes, cell division, and biosynthesis of the cellular envelope. Several cofactor biosynthesis pathways and 29 genes with unknown functions are also required for growth, suggesting novel options for the development of antilisterial drugs. Among the genes specifically required during intracellular growth are known virulence factors, genes compensating intracellular auxotrophies, and several cell division genes. Our experiments also highlight the importance of PASTA kinase signaling for general viability and of glycine metabolism and chromosome segregation for efficient intracellular growth of L. monocytogenes.

Bacterial genome reductions: Tools, applications, and challenges

Bacterial cells are widely used to produce value-added products due to their versatility, ease of manipulation, and the abundance of genome engineering tools. However, the efficiency of producing these desired biomolecules is often hindered by the cells’ own metabolism, genetic instability, and the toxicity of the product. To overcome these challenges, genome reductions have been performed, making strains with the potential of serving as chassis for downstream applications. Here we review the current technologies that enable the design and construction of such reduced-genome bacteria as well as the challenges that limit their assembly and applicability. While genomic reductions have shown improvement of many cellular characteristics, a major challenge still exists in constructing these cells efficiently and rapidly. Computational tools have been created in attempts at minimizing the time needed to design these organisms, but gaps still exist in modelling these reductions in silico. Genomic reductions are a promising avenue for improving the production of value-added products, constructing chassis cells, and for uncovering cellular function but are currently limited by their time-consuming construction methods. With improvements to and the creation of novel genome editing tools and in silico models, these approaches could be combined to expedite this process and create more streamlined and efficient cell factories.

EcoliGD: An Online Tool for Designing Escherichia coli Genome

Synthetic biology is an important interdisciplinary field that has emerged in this century, focusing on the rewriting and reprogramming of DNA through the cycles of "design-edit", and so, the cell's own operating system, its genome, is naturally coming into focus. Here, we propose EcoliGD, an online genome design tool with a visual interactive interface and the function of browsing information, as well as the ability to perform insertion, exchange, deletion, and codon replacement operations on the E. coli genome and display the results in real-time. Users can utilize EcoliGD to check various functional characteristic about E. coli genes, to help them build their genomes. Furthermore, we also collected experimentally verified large genomic segments that have been successfully deleted from the genome for users to choose from and simplify the genome. EcoliGD can help recode the entire E. coli genome, providing a novel way to explore the diversity and function of this microorganism. The EcoliGD web tool is available at http://guolab.whu.edu.cn/EcoliGD/.

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.

Update on the Protein Homeostasis Network in Bacillus subtilis

Protein homeostasis is fundamental to cell function and survival. It relies on an interconnected network of processes involving protein synthesis, folding, post-translational modification and degradation as well as regulators of these processes. Here we provide an update on the roles, regulation and subcellular localization of the protein homeostasis machinery in the Gram-positive model organism Bacillus subtilis. We discuss emerging ideas and current research gaps in the field that, if tackled, increase our understanding of how Gram-positive bacteria, including several human pathogens, maintain protein homeostasis and cope with stressful conditions that challenge their survival.

Effectiveness of the Bacillus subtilis genome-reduced strain as an ethanol production host

We investigated the performance of a genome-reduced strain of Bacillus subtilis MGB874, whose 0.87 Mbp of genomic DNA was cumulatively deleted, as an ethanol production host. A recombinant strain A267_EtOH was constructed by introducing the pdc and adhB genes from Zymomonas mobilis, both of which were expressed from an isopropyl-β-d-1-thiogalactopyranoside-inducible spac promoter, into the A267 strain, a tryptophan prototrophic derivative of the MGB874 with disruption of metabolic pathways for producing lactic acid, acetic acid, and acetoin. Focusing on the stationary phase in fed-batch fermentation, 1.6 g L-1 ethanol was produced by the A267_EtOH strain after 144 h. Moreover, its ethanol production further increased by approximately 3.7-fold (5.9 g L-1) at 80 h through replacing the spac promoter for expressing pdc and adhB genes with the lytR promoter and the yield was about 112%. These results indicate that the MGB874 is an effective host for ethanol production during the stationary phase.

The current state of SubtiWiki, the database for the model organism Bacillus subtilis

Bacillus subtilis is a Gram-positive model bacterium with extensive documented annotation. However, with the rise of high-throughput techniques, the amount of complex data being generated every year has been increasing at a fast pace. Thus, having platforms ready to integrate and give a representation to these data becomes a priority. To address it, SubtiWiki (http://subtiwiki.uni-goettingen.de/) was created in 2008 and has been growing in data and viewership ever since. With millions of requests every year, it is the most visited B. subtilis database, providing scientists all over the world with curated information about its genes and proteins, as well as intricate protein–protein interactions, regulatory elements, expression data and metabolic pathways. However, there is still a large portion of annotation to be unveiled for some biological elements. Thus, to facilitate the development of new hypotheses for research, we have added a Homology section covering potential protein homologs in other organisms. Here, we present the recent developments of SubtiWiki and give a guided tour of our database and the current state of the data for this organism.

The Bacillus subtilis Minimal Genome Compendium

To better understand cellular life, it is essential to decipher the contribution of individual components and their interactions. Minimal genomes are an important tool to investigate these interactions. Here, we provide a database of 105 fully annotated genomes of a series of strains with sequential deletion steps of the industrially relevant model bacterium Bacillus subtilis starting with the laboratory wild type strain B. subtilis 168 and ending with B. subtilis PG38, which lacks approximately 40% of the original genome. The annotation is supported by sequencing of key intermediate strains as well as integration of literature knowledge for the annotation of the deletion scars and their potential effects. The strain compendium presented here represents a comprehensive genome library of the entire MiniBacillus project. This resource will facilitate the more effective application of the different strains in basic science as well as in biotechnology.

A pan-genome method to determine core regions of the Bacillus subtilis and Escherichia coli genomes

Background: Synthetic engineering of bacteria to produce industrial products is a burgeoning field of research and application. In order to optimize genome design, designers need to understand which genes are essential, which are optimal for growth, and locations in the genome that will be tolerated by the organism when inserting engineered cassettes.

Methods: We present a pan-genome based method for the identification of core regions in a genome that are strongly conserved at the species level.

Results: We show that the core regions determined by our method contain all or almost all essential genes. This demonstrates the accuracy of our method as essential genes should be core genes. We show that we outperform previous methods by this measure. We also explain why there are exceptions to this rule for our method.

Conclusions: We assert that synthetic engineers should avoid deleting or inserting into these core regions unless they understand and are manipulating the function of the genes in that region. Similarly, if the designer wishes to streamline the genome, non-core regions and in particular low penetrance genes would be good targets for deletion. Care should be taken to remove entire cassettes with similar penetrance of the genes within cassettes as they may harbor toxin/antitoxin genes which need to be removed in tandem. The bioinformatic approach introduced here saves considerable time and effort relative to knockout studies on single isolates of a given species and captures a broad understanding of the conservation of genes that are core to a species.

SynWiki: Functional annotation of the first artificial organism Mycoplasma mycoides JCVI-syn3A

The new field of synthetic biology aims at the creation of artificially designed organisms. A major breakthrough in the field was the generation of the artificial synthetic organism Mycoplasma mycoides JCVI-syn3A. This bacterium possesses only 452 protein-coding genes, the smallest number for any organism that is viable independent of a host cell. However, about one third of the proteins have no known function indicating major gaps in our understanding of simple living cells. To facilitate the investigation of the components of this minimal bacterium, we have generated the database SynWiki (http://synwiki.uni-goettingen.de/). SynWiki is based on a relational database and gives access to published information about the genes and proteins of M. mycoides JCVI-syn3A. To gain a better understanding of the functions of the genes and proteins of the artificial bacteria, protein-protein interactions that may provide clues for the protein functions are included in an interactive manner. SynWiki is an important tool for the synthetic biology community that will support the comprehensive understanding of a minimal cell as well as the functional annotation of so far uncharacterized proteins.

Quasi-essentiality of RNase Y in Bacillus subtilis is caused by its critical role in the control of mRNA homeostasis

RNA turnover is essential in all domains of life. The endonuclease RNase Y (rny) is one of the key components involved in RNA metabolism of the model organism Bacillus subtilis. Essentiality of RNase Y has been a matter of discussion, since deletion of the rny gene is possible, but leads to severe phenotypic effects. In this work, we demonstrate that the rny mutant strain rapidly evolves suppressor mutations to at least partially alleviate these defects. All suppressor mutants had acquired a duplication of an about 60 kb long genomic region encompassing genes for all three core subunits of the RNA polymerase—α, β, β′. When the duplication of the RNA polymerase genes was prevented by relocation of the rpoA gene in the B. subtilis genome, all suppressor mutants carried distinct single point mutations in evolutionary conserved regions of genes coding either for the β or β’ subunits of the RNA polymerase that were not tolerated by wild type bacteria. In vitro transcription assays with the mutated polymerase variants showed a severe decrease in transcription efficiency. Altogether, our results suggest a tight cooperation between RNase Y and the RNA polymerase to establish an optimal RNA homeostasis in B. subtilis cells.

Functional Redundancy and Specialization of the Conserved Cold Shock Proteins in Bacillus subtilis

Many bacteria encode so-called cold shock proteins. These proteins are characterized by a conserved protein domain. Often, the bacteria have multiple cold shock proteins that are expressed either constitutively or at low temperatures. In the Gram-positive model bacterium Bacillussubtilis, two of three cold shock proteins, CspB and CspD, belong to the most abundant proteins suggesting a very important function. To get insights into the role of these highly abundant proteins, we analyzed the phenotypes of single and double mutants, tested the expression of the csp genes and the impact of CspB and CspD on global gene expression in B. subtilis. We demonstrate that the simultaneous loss of both CspB and CspD results in a severe growth defect, in the loss of genetic competence, and the appearance of suppressor mutations. Overexpression of the third cold shock protein CspC could compensate for the loss of CspB and CspD. The transcriptome analysis revealed that the lack of CspB and CspD affects the expression of about 20% of all genes. In several cases, the lack of the cold shock proteins results in an increased read-through at transcription terminators suggesting that CspB and CspD might be involved in the control of transcription termination.

Non-essential ribosomal proteins in bacteria and archaea identified using COGs

Ribosomal proteins (RPs) are highly conserved across the bacterial and archaeal domains. Although many RPs are essential for survival, genome analysis demonstrates the absence of some RP genes in many bacterial and archaeal genomes. Furthermore, global transposon mutagenesis and/or targeted deletion showed that elimination of some RP genes had only a moderate effect on the bacterial growth rate. Here, we systematically analyze the evolutionary conservation of RPs in prokaryotes by compiling the list of the ribosomal genes that are missing from one or more genomes in the recently updated version of the Clusters of Orthologous Genes (COG) database. Some of these absences occurred because the respective genes carried frameshifts, presumably, resulting from sequencing errors, while others were overlooked and not translated during genome annotation. Apart from these annotation errors, we identified multiple genuine losses of RP genes in a variety of bacteria and archaea. Some of these losses are clade-specific, whereas others occur in symbionts and parasites with dramatically reduced genomes. The lists of computationally and experimentally defined non-essential ribosomal genes show a substantial overlap, revealing a common trend in prokaryote ribosome evolution that could be linked to the architecture and assembly of the ribosomes. Thus, RPs that are located at the surface of the ribosome and/or are incorporated at a late stage of ribosome assembly are more likely to be non-essential and to be lost during microbial evolution, particularly, in the course of genome compaction.IMPORTANCEIn many prokaryote genomes, one or more ribosomal protein (RP) genes are missing. Analysis of 1,309 prokaryote genomes included in the COG database shows that only about half of the RPs are universally conserved in bacteria and archaea. In contrast, up to 16 other RPs are missing in some genomes, primarily, tiny (<1 Mb) genomes of host-associated bacteria and archaea. Ten universal and nine archaea-specific ribosomal proteins show clear patterns of lineage-specific gene loss. Most of the RPs that are frequently lost from bacterial genomes are located on the ribosome periphery and are non-essential in Escherichia coli and Bacillus subtilis These results reveal general trends and common constraints in the architecture and evolution of ribosomes in prokaryotes.

A pan-genome method to determine core regions of the Bacillus subtilis and Escherichia coli genomes

Background: Synthetic engineering of bacteria to produce industrial products is a burgeoning field of research and application. In order to optimize genome design, designers need to understand which genes are essential, which are optimal for growth, and locations in the genome that will be tolerated by the organism when inserting engineered cassettes. Methods: We present a pan-genome based method for the identification of core regions in a genome that are strongly conserved at the species level. Results: We show that the core regions determined by our method contain all or almost all essential genes. This demonstrates the accuracy of our method as essential genes should be core genes. We show that we outperform previous methods by this measure. We also explain why there are exceptions to this rule for our method. Conclusions: We assert that synthetic engineers should avoid deleting or inserting into these core regions unless they understand and are manipulating the function of the genes in that region. Similarly, if the designer wishes to streamline the genome, non-core regions and in particular low penetrance genes would be good targets for deletion. Care should be taken to remove entire cassettes with similar penetrance of the genes within cassettes as they may harbor toxin/antitoxin genes which need to be removed in tandem. The bioinformatic approach introduced here saves considerable time and effort relative to knockout studies on single isolates of a given species and captures a broad understanding of the conservation of genes that are core to a species.

Essentiality of c-di-AMP in Bacillus subtilis: Bypassing mutations converge in potassium and glutamate homeostasis

In order to adjust to changing environmental conditions, bacteria use nucleotide second messengers to transduce external signals and translate them into a specific cellular response. Cyclic di-adenosine monophosphate (c-di-AMP) is the only known essential nucleotide second messenger. In addition to the well-established role of this second messenger in the control of potassium homeostasis, we observed that glutamate is as toxic as potassium for a c-di-AMP-free strain of the Gram-positive model bacterium Bacillus subtilis. In this work, we isolated suppressor mutants that allow growth of a c-di-AMP-free strain under these toxic conditions. Characterization of glutamate resistant suppressors revealed that they contain pairs of mutations, in most cases affecting glutamate and potassium homeostasis. Among these mutations, several independent mutations affected a novel glutamate transporter, AimA (Amino acid importer A, formerly YbeC). This protein is the major transporter for glutamate and serine in B. subtilis. Unexpectedly, some of the isolated suppressor mutants could suppress glutamate toxicity by a combination of mutations that affect phospholipid biosynthesis and a specific gain-of-function mutation of a mechanosensitive channel of small conductance (YfkC) resulting in the acquisition of a device for glutamate export. Cultivation of the c-di-AMP-free strain on complex medium was an even greater challenge because the amounts of potassium, glutamate, and other osmolytes are substantially higher than in minimal medium. Suppressor mutants viable on complex medium could only be isolated under anaerobic conditions if one of the two c-di-AMP receptor proteins, DarA or DarB, was absent. Also on complex medium, potassium and osmolyte toxicity are the major bottlenecks for the growth of B. subtilis in the absence of c-di-AMP. Our results indicate that the essentiality of c-di-AMP in B. subtilis is caused by the global impact of the second messenger nucleotide on different aspects of cellular physiology.

Two-Component-System RspA1/A2-Dependent Regulation on Primary Metabolism in Streptomyces albus A30 Cultivated With Glutamate as the Sole Nitrogen Source

In our previous study, a two-component-system (TCS) RspA1/A2 was identified and proven to play a positive role in the regulation of salinomycin (antibiotic) biosynthesis in Streptomyces albus. However, the regulatory mechanism of RspA1/A2 using a carbon source (glucose or acetate) for the cell growth of S. albus is still unclear till present research work. Therefore, in this work, the mechanistic pathway of RspA1/A2 on carbon source metabolism is unveiled. Firstly, this work reports that the response regulator RspA1 gene rspA1 knocked-out mutant ΔrspA1 exhibits lower biomass accumulation and lower glucose consumption rates as compared to the parental strain A30 when cultivated in a defined minimal medium (MM) complemented with 75 mM glutamate. Further, it is demonstrated that the regulation of TCS RspA1/A2 on the phosphoenolpyruvate-pyruvate-oxaloacetate node results in decreasing the intracellular acetyl-CoA pool in mutant ΔrspA1. Subsequently, it was verified that the RspA1 could not only directly interact with the promoter regions of key genes encoding AMP-forming acetyl-CoA synthase (ACS), citrate synthase (CS), and pyruvate dehydrogenase complex (PDH) but also bind promoter regions of the genes pyc, pck, and glpX in gluconeogenesis. In addition, the transcriptomic data analysis showed that pyruvate and glutamate transformations supported robust TCS RspA1/A2-dependent regulation of glucose metabolism, which led to a decreased flux of pyruvate into the TCA cycle and an increased flux of gluconeogenesis pathway in mutant ΔrspA1. Finally, a new transcriptional regulatory network of TCS RspA1/A2 on primary metabolism across central carbon metabolic pathways including the glycolysis pathway, TCA cycle, and gluconeogenesis pathway is proposed.

Advances on systems metabolic engineering of Bacillus subtilis as a chassis cell

The Gram-positive model bacterium Bacillus subtilis, has been broadly applied in various fields because of its low pathogenicity and strong protein secretion ability, as well as its well-developed fermentation technology. B. subtilis is considered as an attractive host in the field of metabolic engineering, in particular for protein expression and secretion, so it has been well studied and applied in genetic engineering. In this review, we discussed why B. subtilis is a good chassis cell for metabolic engineering. We also summarized the latest research progress in systematic biology, synthetic biology and evolution-based engineering of B. subtilis, and showed systemic metabolic engineering expedite the harnessing B. subtilis for bioproduction.

Zinc oxide nanoparticles decorated fluorescent and antibacterial glass fiber pre-filter paper

Zinc oxide nanoparticles (ZnO–NPs) were synthesized and decorated simultaneously onto the glass fiber pre-filter paper (GF paper) by the sonochemical method without using any additional reagents (a ‘Green’ synthesis approach). ZnO–NPs decorated GF paper was characterized by electron, confocal laser scanning and atomic force microscopy, fourier transform infrared and atomic emission spectroscopy, X-ray diffraction, and thermogravimetric analysis etc. Due to the massive void volume space, exceptional dimensional stability, large thickness (790 μm) of the GF paper (unlike other paper materials) and ultrasonic irradiation effects, ZnO–NPs were decorated in the enormous amount (96 mg per paper) without causing any adverse effects on the GF paper. Such a huge amount decoration onto GF paper makes it multifunctional, fluorescencet (orange-pink color, 535–624 nm) under ultra-violet light (360 nm) and antibacterial. The antibacterial activity of the ZnO–NPs decorated GF paper was examined against Gram-positive bacteria Bacillus subtilis 168 and Staphylococcus aureus (MCC 2043, pathogenic). The outcomes from the antibacterial experiments revealed ∼99% (2 log) reduction in the survival of the filtered bacteria (B. subtilis) on the ZnO–NPs decorated GF paper due to the toxicity of ZnO–NPs on bacterial cells like cell shrinkage, cytoplasmic leakage, cell burst, etc. Multifunctional, ZnO–NPs decorated GF paper could be used for fluorescencet and antibacterial paper-based applications.

Dissecting the Functional Contributions of the Intrinsically Disordered C-terminal Tail of Bacillus subtilis FtsZ

FtsZ is a bacterial GTPase that is central to the spatial and temporal control of cell division. It is a filament-forming enzyme that encompasses a well-folded core domain and a disordered C-terminal tail (CTT). The CTT is essential for ensuring proper assembly of the cytokinetic ring, and its deletion leads to mis-localization of FtsZ, aberrant assembly, and cell death. In this work, we dissect the contributions of modules within the disordered CTT to assembly and enzymatic activity of Bacillus subtilis FtsZ (Bs-FtsZ). The CTT features a hypervariable C-terminal linker (CTL) and a conserved C-terminal peptide (CTP). Our in vitro studies show that the CTL weakens the driving forces for forming single-stranded active polymers and suppresses lateral associations of these polymers, whereas the CTP promotes the formation of alternative assemblies. Accordingly, in full-length Bs-FtsZ, the CTL acts as a spacer that spatially separates the CTP sticker from the core, thus ensuring filament formation through core-driven polymerization and lateral associations through CTP-mediated interactions. We also find that the CTL weakens GTP binding while enhancing the catalytic rate, whereas the CTP has opposite effects. The joint contributions of the CTL and CTP make Bs-FtsZ, an enzyme that is only half as efficient as a truncated version that lacks the CTT. Overall, our data suggest that the CTT acts as an auto-regulator of Bs-FtsZ assembly and as an auto-inhibitor of enzymatic activity. Based on our results, we propose hypotheses regarding the hypervariability of CTLs and compare FtsZs to other bacterial proteins with tethered IDRs.

Designing minimal genomes using whole-cell models

In the future, entire genomes tailored to specific functions and environments could be designed using computational tools. However, computational tools for genome design are currently scarce. Here we present algorithms that enable the use of design-simulate-test cycles for genome design, using genome minimisation as a proof-of-concept. Minimal genomes are ideal for this purpose as they have a simple functional assay whether the cell replicates or not. We used the first (and currently only published) whole-cell model for the bacterium Mycoplasma genitalium. Our computational design-simulate-test cycles discovered novel in silico minimal genomes which, if biologically correct, predict in vivo genomes smaller than JCVI-Syn3.0; a bacterium with, currently, the smallest genome that can be grown in pure culture. In the process, we identified 10 low essential genes and produced evidence for at least two Mycoplasma genitalium in silico minimal genomes. This work brings combined computational and laboratory genome engineering a step closer.

Genome engineering will one day benefit from computational tools that can design genomes with desired functions. Here the authors develop computational design-simulate-test algorithms to design minimal genomes based on the whole-cell model of Mycoplasma genitalium.

Multiple integration of the gene ganA into the Bacillus subtilis chromosome for enhanced β-galactosidase production using the CRISPR/Cas9 system

The ganA gene from Bacillus subtilis encoding a β-galactosidase for degradation of the galactomannan was integrated in different loci of the B. subtilis chromosome employing the CRISPR/Cas9 system. Hereby a total of five copies of ganA cassettes in which the ganA gene was fused with the glucitol-promoter were inserted in the recipient chromosome wherein hypothetical, sporulation and protease genes were deleted. The strain with five copies of ganA expression cassette showed a β-galactosidase activity similar to the one with the same gene on a pUB110 derived multi-copy plasmid and under the same regulatory control of the glucitol promoter and GutR activator. The production of β-galactosidase in the strain with the multi-copy plasmid decreased rapidly when growth was performed under induced conditions and without antibiotic selection. In contrast, the strain with the five copies of ganA in the chromosome produced β-galactosidase for at least 40 generations. This demonstrates that the CRISPR/Cas9 system is a valuable and easy tool for constructing stable producer strains. The bigger efforts that are needed for the multiple target gene integration into the chromosome compared to cloning in expression vectors were justified by the higher stability of the target genes and the lack of antibiotic resistance genes.

The conditional nature of gene essentiality

Essential genes are classically defined as required for cellular viability and reproductive success. Despite this deceptively simple definition, several lines of evidence suggest that gene essentiality is instead a conditional trait. Indeed, gene essentiality has been shown to depend on the environmental and genetic context as well as the variable ability of cells to acquire adaptive mutations to survive inactivation of seemingly essential genes. Here, we will discuss these findings and highlight the mechanisms underlying the ability of cells to survive an essential gene deletion. Also, since essential genes are prioritized as targets for anticancer therapy, we discuss emergence of bypass resistance mechanisms toward targeted therapies as the result of the conditional nature of gene essentiality. To identify targets associated to a lower risk of relapse (i.e. the return of cancer following remission), we finally call for a coordinated effort to quantify the variable nature of gene essentiality across species, cell types, and growth conditions.

SubtiWiki in 2018: from genes and proteins to functional network annotation of the model organism Bacillus subtilis

Living cells are made up of individual parts, i.e. the genome, the proteins, the RNA and lipid molecules as well as the metabolites and ions. However, life depends on the functional interaction among these components which is often organized in networks. Here, we present the recent development of SubtiWiki, the integrated database for the model bacterium Bacillus subtilis (http://subtiwiki.uni-goettingen.de/). SubtiWiki is based on a relational database and provides access to published information about the genes and proteins of B. subtilis and about metabolic and regulatory pathways. We have included a network visualization tool that can be used to visualize regulatory as well as protein-protein interaction networks. The resulting interactive graphical presentations allow the user to detect novel associations and thus to develop novel hypotheses that can then be tested experimentally. To facilitate the mobile use of SubtiWiki, we provide enhanced versions of the SubtiWiki App that are available for iOS and Android devices. Importantly, the App allows to link private notes and pictures to the gene/protein pages that can be synchronized on multiple devices. SubtiWiki has become one of the most complete resources of knowledge on a living organism.

Topoisomerase IV can functionally replace all type 1A topoisomerases in Bacillus subtilis

DNA topoisomerases play essential roles in chromosome organization and replication. Most bacteria possess multiple topoisomerases which have specialized functions in the control of DNA supercoiling or in DNA catenation/decatenation during recombination and chromosome segregation. DNA topoisomerase I is required for the relaxation of negatively supercoiled DNA behind the transcribing RNA polymerase. Conflicting results have been reported on the essentiality of the topA gene encoding topoisomerase I in the model bacterium Bacillus subtilis. In this work, we have studied the requirement for topoisomerase I in B. subtilis. All stable topA mutants carried different chromosomal amplifications of the genomic region encompassing the parEC operon encoding topoisomerase IV. Using a fluorescent amplification reporter system we observed that each individual topA mutant had acquired such an amplification. Eventually, the amplifications were replaced by a point mutation in the parEC promoter region which resulted in a fivefold increase of parEC expression. In this strain both type I topoisomerases, encoded by topA and topB, were dispensable. Our results demonstrate that topoisomerase IV at increased expression is necessary and sufficient to take over the function of type 1A topoisomerases.

The essential genome of the crenarchaeal model Sulfolobus islandicus

Sulfolobus islandicus is a model microorganism in the TACK superphylum of the Archaea, a key lineage in the evolutionary history of cells. Here we report a genome-wide identification of the repertoire of genes essential to S. islandicus growth in culture. We confirm previous targeted gene knockouts, uncover the non-essentiality of functions assumed to be essential to the Sulfolobus cell, including the proteinaceous S-layer, and highlight essential genes whose functions are yet to be determined. Phyletic distributions illustrate the potential transitions that may have occurred during the evolution of this archaeal microorganism, and highlight sets of genes that may have been associated with each transition. We use this comparative context as a lens to focus future research on archaea-specific uncharacterized essential genes that may provide valuable insights into the evolutionary history of cells.

Sulfolobus islandicus is a model organism within the TACK superphylum of the Archaea. Here, the authors perform a genome-wide analysis of essential genes in this organism, show that the proteinaceous S-layer is not essential, and explore potential stages of evolution of the essential gene repertoire in Archaea.

Stimulation of Vibrio vulnificus Pyruvate Kinase in the Presence of Glucose to Cope With H2O2 Stress Generated by Its Competitors

The bacterial phosphoenolpyruvate (PEP):carbohydrate phosphotransferase system (PTS) regulates a variety of cellular processes in addition to catalyzing the coupled transport and phosphorylation of carbohydrates. We recently reported that, in the presence of glucose, HPr of the PTS is dephosphorylated and interacts with pyruvate kinase A (PykA) catalyzing the conversion of PEP to pyruvate in Vibrio vulnificus. Here, we show that this interaction enables V. vulnificus to survive H₂O₂ stress by increasing pyruvate production. A pykA deletion mutant was more susceptible to H₂O₂ stress than wild-type V. vulnificus without any decrease in the expression level of catalase, and this sensitivity was rescued by the addition of pyruvate. The H₂O₂ sensitivity difference between wild-type and pykA mutant strains becomes more apparent in the presence of glucose. Fungi isolated from the natural habitat of V. vulnificus retarded the growth of the pykA mutant more severely than the wild-type strain in the presence of glucose by glucose oxidase-dependent generation of H₂O₂. These data suggest that V. vulnificus has evolved to resist the killing action of its fungal competitors by increasing pyruvate production in the presence of glucose.

Changes of DNA topology affect the global transcription landscape and allow rapid growth of a Bacillus subtilis mutant lacking carbon catabolite repression

Bacteria are able to prioritize preferred carbon sources from complex mixtures. This is achieved by the regulatory phenomenon of carbon catabolite repression. To allow the simultaneous utilization of multiple carbon sources and to prevent the time-consuming adaptation to each individual nutrient in biotechnological applications, mutants lacking carbon catabolite repression can be used. However, such mutants often exhibit pleiotropic growth defects. We have isolated and characterized mutations that overcome the growth defect of Bacillus subtilis ccpA mutants lacking the major regulator of catabolite repression, in particular their glutamate auxotrophy. Here we show, that distinct mutations affecting the essential DNA topoisomerase I (TopA) cause glutamate prototrophy of the ccpA mutant. These suppressing variants of the TopA enzyme exhibit increased activity resulting in enhanced relaxation of the DNA. Reduced DNA supercoiling results in enhanced expression of the gltAB operon encoding the biosynthetic glutamate synthase. This is achieved by a significant re-organization of the global transcription network accompanied by re-routing of metabolism, which results in inactivation of the glutamate dehydrogenase. Our results provide a link between DNA topology, the global transcriptional network, and glutamate metabolism and suggest that specific topA mutants may be well suited for biotechnological purposes.

Emerging and evolving concepts in gene essentiality

Gene essentiality is a founding concept of genetics with important implications in both fundamental and applied research. Multiple screens have been performed over the years in bacteria, yeasts, animals and more recently in human cells to identify essential genes. A mounting body of evidence suggests that gene essentiality, rather than being a static and binary property, is both context dependent and evolvable in all kingdoms of life. This concept of a non-absolute nature of gene essentiality changes our fundamental understanding of essential biological processes and could directly affect future treatment strategies for cancer and infectious diseases.

Developing genome-reduced Pseudomonas chlororaphis strains for the production of secondary metabolites

BACKGROUND

The current chassis organisms or various types of cell factories have considerable advantages and disadvantages. Therefore, it is necessary to develop various chassis for an efficient production of different bioproducts from renewable resources. In this context, synthetic biology offers unique potentialities to produce value-added products of interests. Microbial genome reduction and modification are important strategies for constructing cellular chassis and cell factories. Many genome-reduced strains from Escherichia coli, Bacillus subtilis, Corynebacterium glutamicum and Streptomyces, have been widely used for the production of amino acids, organic acids, and some enzymes. Some Pseudomonas strains could serve as good candidates for ideal chassis cells since they grow fast and can produce many valuable metabolites with low nutritional requirements and strong environmental adaptability. Pseudomonas chlororaphis GP72 is a non-pathogenic plant growth-promoting rhizobacterium that possesses capacities of tolerating various environmental stresses and synthesizing many kinds of bioactive compounds with high yield. These include phenazine-1-carboxylic acid (PCA) and 2-hydroxyphenazine (2-OH-PHZ), which exhibit strong bacteriostatic and antifungal activity toward some microbial pathogens.

RESULTS

We depleted 685 kb (10.3% of the genomic sequence) from the chromosome of P. chlororaphis GP72(rpeA-) by a markerless deletion method, which included five secondary metabolic gene clusters and 17 strain-specific regions (525 non-essential genes). Then we characterized the 22 multiple-deletion series (MDS) strains. Growth characteristics, production of phenazines and morphologies were changed greatly in mutants with large-fragment deletions. Some of the genome-reduced P. chlororaphis mutants exhibited more productivity than the parental strain GP72(rpeA-). For example, strain MDS22 had 4.4 times higher production of 2-OH-PHZ (99.1 mg/L) than strain GP72(rpeA-), and the specific 2-OH-PHZ production rate (mmol/g/h) increased 11.5-fold. Also and MDS10 had the highest phenazine production (852.0 mg/L) among all the studied strains with a relatively high specific total phenazine production rate (0.0056 g/g/h).

CONCLUSIONS

In conclusion, P. chlororaphis strains with reduced genome performed better in production of secondary metabolites than the parent strain. The newly developed mutants can be used for the further genetic manipulation to construct chassis cells with the less complex metabolic network, better regulation and more efficient productivity for diverse biotechnological applications.

Prediction and analysis of essential genes using the enrichments of gene ontology and KEGG pathways

Identifying essential genes in a given organism is important for research on their fundamental roles in organism survival. Furthermore, if possible, uncovering the links between core functions or pathways with these essential genes will further help us obtain deep insight into the key roles of these genes. In this study, we investigated the essential and non-essential genes reported in a previous study and extracted gene ontology (GO) terms and biological pathways that are important for the determination of essential genes. Through the enrichment theory of GO and KEGG pathways, we encoded each essential/non-essential gene into a vector in which each component represented the relationship between the gene and one GO term or KEGG pathway. To analyze these relationships, the maximum relevance minimum redundancy (mRMR) was adopted. Then, the incremental feature selection (IFS) and support vector machine (SVM) were employed to extract important GO terms and KEGG pathways. A prediction model was built simultaneously using the extracted GO terms and KEGG pathways, which yielded nearly perfect performance, with a Matthews correlation coefficient of 0.951, for distinguishing essential and non-essential genes. To fully investigate the key factors influencing the fundamental roles of essential genes, the 21 most important GO terms and three KEGG pathways were analyzed in detail. In addition, several genes was provided in this study, which were predicted to be essential genes by our prediction model. We suggest that this study provides more functional and pathway information on the essential genes and provides a new way to investigate related problems.

T7 RNA polymerase-driven inducible cell lysis for DNA transfer from Escherichia coli to Bacillus subtilis

The majority of the good DNA editing techniques have been developed in Escherichia coli; however, Bacillus subtilis is better host for a plethora of synthetic biology and biotechnology applications. Reliable and efficient systems for the transfer of synthetic DNA between E. coli and B. subtilis are therefore of the highest importance. Using synthetic biology approaches, such as streamlined lambda Red recombineering and Gibson Isothermal Assembly, we integrated genetic circuits pT7L123, Repr‐ts‐1 and pLT7pol encoding the lysis genes of bacteriophages MS2, ΦX174 and lambda, the thermosensitive repressor and the T7 RNA polymerase into the E. coli chromosome. In this system, T7 RNA polymerase regulated by the thermosensitive repressor drives the expression of the phage lysis genes. We showed that T7 RNA polymerase significantly increases efficiency of cell lysis and transfer of the plasmid and bacterial artificial chromosome‐encoded DNA from the lysed E. coli into B. subtilis. The T7 RNA polymerase‐driven inducible cell lysis system is suitable for the efficient cell lysis and transfer of the DNA engineered in E. coli to other naturally competent hosts, such as B. subtilis.

The Highly Conserved Asp23 Family Protein YqhY Plays a Role in Lipid Biosynthesis in Bacillus subtilis

In most bacteria, fatty acid biosynthesis is an essential process that must be controlled by the availability of precursors and by the needs of cell division. So far, no mechanisms controlling synthesis of malonyl-coenzyme A (CoA), the committed step in fatty acid synthesis, have been identified in the Gram-positive model bacterium Bacillus subtilis. We have studied the localization and function of two highly expressed proteins of unknown function, YqhY and YloU. Both proteins are members of the conserved and widespread Asp23 family. While the deletion of yloU had no effect, loss of the yqhY gene induced the rapid acquisition of suppressor mutations. The vast majority of these mutations affect subunits of the acetyl-CoA carboxylase (ACCase) complex, the enzyme that catalyzes the formation of malonyl-CoA. Moreover, lack of yqhY is accompanied by the formation of lipophilic clusters in the polar regions of the cells indicating an increased activity of ACCase. Our results suggest that YqhY controls the activity of ACCase and that this control results in inhibition of ACCase activity. Hyperactivity of the enzyme complex in the absence of YqhY does then provoke mutations that cause reduced ACCase activity.

Selection of key sequence-based features for prediction of essential genes in 31 diverse bacterial species

Genes that are indispensable for survival are essential genes. Many features have been proposed for computational prediction of essential genes. In this paper, the least absolute shrinkage and selection operator method was used to screen key sequence-based features related to gene essentiality. To assess the effects, the selected features were used to predict the essential genes from 31 bacterial species based on a support vector machine classifier. For all 31 bacterial objects (21 Gram-negative objects and ten Gram-positive objects), the features in the three datasets were reduced from 57, 59, and 58, to 40, 37, and 38, respectively, without loss of prediction accuracy. Results showed that some features were redundant for gene essentiality, so could be eliminated from future analyses. The selected features contained more complex (or key) biological information for gene essentiality, and could be of use in related research projects, such as gene prediction, synthetic biology, and drug design.

Construction and Analysis of Two Genome-Scale Deletion Libraries for Bacillus subtilis

A systems-level understanding of Gram-positive bacteria is important from both an environmental and health perspective and is most easily obtained when high-quality, validated genomic resources are available. To this end, we constructed two ordered, barcoded, erythromycin-resistance- and kanamycin-resistance-marked single-gene deletion libraries of the Gram-positive model organism, Bacillus subtilis. The libraries comprise 3,968 and 3,970 genes, respectively, and overlap in all but four genes. Using these libraries, we update the set of essential genes known for this organism, provide a comprehensive compendium of B. subtilis auxotrophic genes, and identify genes required for utilizing specific carbon and nitrogen sources, as well as those required for growth at low temperature. We report the identification of enzymes catalyzing several missing steps in amino acid biosynthesis. Finally, we describe a suite of high-throughput phenotyping methodologies and apply them to provide a genome-wide analysis of competence and sporulation. Altogether, we provide versatile resources for studying gene function and pathway and network architecture in Gram-positive bacteria.