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

Membrane Localization of RNase Y Is Important for Global Gene Expression in Bacillus subtilis

RNase Y is a key endoribonuclease that regulates global mRNA turnover and processing in Bacillus subtilis and likely many other bacteria. This enzyme is anchored to the cell membrane, creating a pseudo-compartmentalization that aligns with its role in initiating the decay of mRNAs primarily translated at the cell periphery. However, the reasons behind and the consequences of RNase Y's membrane attachment remain largely unknown. In our study, we examined a strain expressing wild-type levels of a cytoplasmic form of RNase Y from its chromosomal locus. This strain exhibits a slow-growth phenotype, similar to that of an RNase Y null mutant. Genome-wide data reveal a significant impact on the expression of hundreds of genes. While certain RNA substrates clearly depend on RNase Y's membrane attachment, others do not. We observed no correlation between mRNA stabilization in the mutant strains and the cellular location or function of the encoded proteins. Interestingly, the Y-complex, a specificity factor for RNase Y, also appears also recognize the cytoplasmic form of the enzyme, restoring wild-type levels of the corresponding transcripts. We propose that membrane attachment of RNase Y is crucial for its functional interaction with many coding and non-coding RNAs, limiting the cleavage of specific substrates, and potentially avoiding unfavorable competition with other ribonucleases like RNase J, which shares a similar evolutionarily conserved cleavage specificity.

The Construction of an Environmentally Friendly Super-Secreting Strain of Bacillus subtilis through Systematic Modulation of Its Secretory Pathway Using the CRISPR-Cas9 System

Achieving commercially significant yields of recombinant proteins in Bacillus subtilis requires the optimization of its protein production pathway, including transcription, translation, folding, and secretion. Therefore, in this study, our aim was to maximize the secretion of a reporter α-amylase by overcoming potential bottlenecks within the secretion process one by one, using a clustered regularly interspaced short palindromic repeat-Cas9 (CRISPR-Cas9) system. The strength of single and tandem promoters was evaluated by measuring the relative α-amylase activity of AmyQ integrated into the B. subtilis chromosome. Once a suitable promoter was selected, the expression levels of amyQ were upregulated through the iterative integration of up to six gene copies, thus boosting the α-amylase activity 20.9-fold in comparison with the strain harboring a single amyQ gene copy. Next, α-amylase secretion was further improved to a 26.4-fold increase through the overexpression of the extracellular chaperone PrsA and the signal peptide peptidase SppA. When the final expression strain was cultivated in a 3 L fermentor for 90 h, the AmyQ production was enhanced 57.9-fold. The proposed strategy allows for the development of robust marker-free plasmid-less super-secreting B. subtilis strains with industrial relevance.

Large language model based framework for automated extraction of genetic interactions from unstructured data

Extracting biological interactions from published literature helps us understand complex biological systems, accelerate research, and support decision-making in drug or treatment development. Despite efforts to automate the extraction of biological relations using text mining tools and machine learning pipelines, manual curation continues to serve as the gold standard. However, the rapidly increasing volume of literature pertaining to biological relations poses challenges in its manual curation and refinement. These challenges are further compounded because only a small fraction of the published literature is relevant to biological relation extraction, and the embedded sentences of relevant sections have complex structures, which can lead to incorrect inference of relationships. To overcome these challenges, we propose GIX, an automated and robust Gene Interaction Extraction framework, based on pre-trained Large Language models fine-tuned through extensive evaluations on various gene/protein interaction corpora including LLL and RegulonDB. GIX identifies relevant publications with minimal keywords, optimises sentence selection to reduce computational overhead, simplifies sentence structure while preserving meaning, and provides a confidence factor indicating the reliability of extracted relations. GIX's Stage-2 relation extraction method performed well on benchmark protein/gene interaction datasets, assessed using 10-fold cross-validation, surpassing state-of-the-art approaches. We demonstrated that the proposed method, although fully automated, performs as well as manual relation extraction, with enhanced robustness. We also observed GIX's capability to augment existing datasets with new sentences, incorporating newly discovered biological terms and processes. Further, we demonstrated GIX's real-world applicability in inferring E. coli gene circuits.

Two distinct regulatory systems control pulcherrimin biosynthesis in Bacillus subtilis

Regulation of transcription is a fundamental process that allows bacteria to respond to external stimuli with appropriate timing and magnitude of response. In the soil bacterium Bacillus subtilis, transcriptional regulation is at the core of developmental processes needed for cell survival. Gene expression in cells transitioning from exponential phase to stationary phase is under the control of a group of transcription factors called transition state regulators (TSRs). TSRs influence numerous developmental processes including the decision between biofilm formation and motility, genetic competence, and sporulation, but the extent to which TSRs influence bacterial physiology remains to be fully elucidated. Here, we demonstrate two TSRs, ScoC and AbrB, along with the MarR-family transcription factor PchR negatively regulate production of the iron chelator pulcherrimin in B. subtilis. Genetic analysis of the relationship between the three transcription factors indicate that all are necessary to limit pulcherrimin production during exponential phase and influence the rate and total amount of pulcherrimin produced. Similarly, expression of the pulcherrimin biosynthesis gene yvmC was found to be under control of ScoC, AbrB, and PchR and correlated with the amount of pulcherrimin produced by each background. Lastly, our in vitro data indicate a weak direct role for ScoC in controlling pulcherrimin production along with AbrB and PchR. The layered regulation by two distinct regulatory systems underscores the important role for pulcherrimin in B. subtilis physiology.

PMF-GRN: a variational inference approach to single-cell gene regulatory network inference using probabilistic matrix factorization

Inferring gene regulatory networks (GRNs) from single-cell data is challenging due to heuristic limitations. Existing methods also lack estimates of uncertainty. Here we present Probabilistic Matrix Factorization for Gene Regulatory Network Inference (PMF-GRN). Using single-cell expression data, PMF-GRN infers latent factors capturing transcription factor activity and regulatory relationships. Using variational inference allows hyperparameter search for principled model selection and direct comparison to other generative models. We extensively test and benchmark our method using real single-cell datasets and synthetic data. We show that PMF-GRN infers GRNs more accurately than current state-of-the-art single-cell GRN inference methods, offering well-calibrated uncertainty estimates.

MOB rules: Antibiotic Exposure Reprograms Metabolism to Mobilize Bacillus subtilis in Competitive Interactions

Antibiotics have dose-dependent effects on exposed bacteria. The medicinal use of antibiotics relies on their growth-inhibitory activities at sufficient concentrations. At subinhibitory concentrations, exposure effects vary widely among different antibiotics and bacteria. Bacillus subtilis responds to bacteriostatic translation inhibitors by mobilizing a population of cells (MOB-Mobilized Bacillus) to spread across a surface. How B. subtilis regulates the antibiotic-induced mobilization is not known. In this study, we used chloramphenicol to identify regulatory functions that B. subtilis requires to coordinate cell mobilization following subinhibitory exposure. We measured changes in gene expression and metabolism and mapped the results to a network of regulatory proteins that direct the mobile response. Our data reveal that several transcriptional regulators coordinately control the reprogramming of metabolism to support mobilization. The network regulates changes in glycolysis, nucleotide metabolism, and amino acid metabolism that are signature features of the mobilized population. Among the hundreds of genes with changing expression, we identified two, pdhA and pucA, where the magnitudes of their changes in expression, and in the abundance of associated metabolites, reveal hallmark metabolic features of the mobilized population. Using reporters of pdhA and pucA expression, we visualized the separation of major branches of metabolism in different regions of the mobilized population. Our results reveal a regulated response to chloramphenicol exposure that enables a population of bacteria in different metabolic states to mount a coordinated mobile response.

Improving gene set enrichment analysis (GSEA) by using regulation directionality

To infer the biological meaning from transcriptome data, it is useful to focus on genes that are regulated by the same regulator, i.e., regulons. Unfortunately, current gene set enrichment analysis (GSEA) tools do not consider whether a gene is activated or repressed by a regulator. This distinction is crucial when analyzing regulons since a regulator can work as an activator of certain genes and as a repressor of other genes, yet both sets of genes belong to the same regulon. Therefore, simply averaging expression differences of the genes of such a regulon will not properly reflect the activity of the regulator. What makes it more complicated is the fact that many genes are regulated by different transcription factors, and current transcriptome analysis tools are unable to indicate which regulator is most likely responsible for the observed expression difference of a gene. To address these challenges, we developed the gene set enrichment analysis program GINtool. Additional features of GINtool are novel graphical representations to facilitate the visualization of gene set analyses of transcriptome data, the possibility to include functional categories as gene sets for analysis, and the option to analyze expression differences within operons, which is useful when analyzing prokaryotic transcriptome and also proteome data.IMPORTANCEMeasuring the activity of all genes in cells is a common way to elucidate the function and regulation of genes. These transcriptome analyses produce large amounts of data since genomes contain thousands of genes. The analysis of these large data sets is challenging. Therefore, we developed a new software tool called GINtool that can facilitate the analysis of transcriptome data by using prior knowledge of gene sets controlled by the same regulator, the so-called regulons. An important novelty of GINtool is that it can take into account the directionality of gene regulation in these analyses, i.e., whether a gene is activated or repressed, which is crucial to assess whether a regulon or functional category is affected. GINtool also includes new graphical methods to facilitate the visual inspection of regulation events in transcriptome data sets. These and additional analysis methods included in GINtool make it a powerful software tool to analyze transcriptome data.

Response of Bacillus velezensis 83 to interaction with Colletotrichum gloeosporioides resembles a Greek phalanx-style formation: A stress resistant phenotype with antibiosis capacity

Plant growth-promoting rhizobacteria, such as Bacillus spp., establish beneficial associations with plants and may inhibit the growth of phytopathogenic fungi. However, these bacteria are subject to multiple biotic stimuli from their competitors, causing stress and modifying their development. This work is a study of an in vitro interaction between two model microorganisms of socioeconomic relevance, using population dynamics and transcriptomic approaches. Co-cultures of Bacillus velezensis 83 with the phytopathogenic fungus Colletotrichum gloeosporioides 09 were performed to evaluate the metabolic response of the bacteria under conditions of non-nutritional limitation. The bacterial response was associated with the induction of a stress-resistant phenotype, characterized by a lower specific growth rate, but with antimicrobial production capacity. About 12% of co-cultured B. velezensis 83 coding sequences were differentially expressed, including the up-regulation of the general stress response (sigB regulon), and the down-regulation of alternative carbon sources catabolism (glucose preference). Defense strategies in B. velezensis are a determining factor in order to preserve the long-term viability of its population. Mostly, the presence of the fungus does not affect the expression of antibiosis genes, except for those corresponding to surfactin/bacillomycin D production. Indeed, the up-regulation of antibiosis genes expression is associated with bacterial growth, regardless of the presence of the fungus. This behavior in B. velezensis 83 resembles the strategy used by the classical Greek phalanx formation: by sacrificing growth rate and metabolic versatility, resources can be redistributed to defense (stress resistant phenotype) while maintaining the attack (antibiosis capacity). The presented results are the first characterization of the molecular phenotype at the transcriptome level of a biological control agent under biotic stress caused by a phytopathogen without nutrient limitation.

Two Distinct Regulatory Systems Control Pulcherrimin Biosynthesis in Bacillus subtilis

Regulation of transcription is a fundamental process that allows bacteria to respond to external stimuli with appropriate timing and magnitude of response. In the soil bacterium Bacillus subtilis, transcriptional regulation is at the core of developmental processes needed for cell survival. Gene expression in cells transitioning from exponential phase to stationary phase is under the control of a group of transcription factors called transition state regulators (TSRs). TSRs influence numerous developmental processes including the decision between biofilm formation and motility, genetic competence, and sporulation, but the extent to which TSRs influence bacterial physiology remains to be fully elucidated. Here, we demonstrate two TSRs, ScoC and AbrB, along with the MerR-family transcription factor PchR negatively regulate production of the iron chelator pulcherrimin in B. subtilis. Genetic analysis of the relationship between the three transcription factors indicate that all are necessary to limit pulcherrimin production during exponential phase and influence the rate and total amount of pulcherrimin produced. Similarly, expression of the pulcherrimin biosynthesis gene yvmC was found to be under control of ScoC, AbrB, and PchR and correlated with the amount of pulcherrimin produced by each background. Lastly, our in vitro data indicate a weak direct role for ScoC in controlling pulcherrimin production along with AbrB and PchR. The layered regulation by two distinct regulatory systems underscores the important, and somewhat enigmatic, role for pulcherrimin in B. subtilis physiology.

Arbitrium communication controls phage lysogeny through non-lethal modulation of a host toxin-antitoxin defence system

Temperate Bacillus phages often utilize arbitrium communication to control lysis/lysogeny decisions, but the mechanisms by which this control is exerted remains largely unknown. Here we find that the arbitrium system of Bacillus subtilis phage ϕ3T modulates the host-encoded MazEF toxin-antitoxin system to this aim. Upon infection, the MazF ribonuclease is activated by three phage genes. At low arbitrium signal concentrations, MazF is inactivated by two phage-encoded MazE homologues: the arbitrium-controlled AimX and the later-expressed YosL proteins. At high signal, MazF remains active, promoting lysogeny without harming the bacterial host. MazF cleavage sites are enriched on transcripts of phage lytic genes but absent from the phage repressor in ϕ3T and other Spβ-like phages. Combined with low activation levels of MazF during infections, this pattern explains the phage-specific effect. Our results show how a bacterial toxin-antitoxin system has been co-opted by a phage to control lysis/lysogeny decisions without compromising host viability.

Induction of the CtsR regulon improves Xylanase production in Bacillus subtilis

BACKGROUND

The bacterium Bacillus subtilis is extensively used for the commercial production of enzymes due to its efficient protein secretion capacity. However, the efficiency of secretion varies greatly between enzymes, and despite many years of research, optimization of enzyme production is still largely a matter of trial-and-error. Genome-wide transcriptome analysis seems a useful tool to identify relevant secretion bottlenecks, yet to this day, only a limited number of transcriptome studies have been published that focus on enzyme secretion in B. subtilis. Here, we examined the effect of high-level expression of the commercially important enzyme endo-1,4-β-xylanase XynA on the B. subtilis transcriptome using RNA-seq.

RESULTS

Using the novel gene-set analysis tool GINtool, we found a reduced activity of the CtsR regulon when XynA was overproduced. This regulon comprises several protein chaperone genes, including clpC, clpE and clpX, and is controlled by transcriptional repression. CtsR levels are directly controlled by regulated proteolysis, involving ClpC and its cognate protease ClpP. When we abolished this negative feedback, by inactivating the repressor CtsR, the XynA production increased by 25%.

CONCLUSIONS

Overproduction of enzymes can reduce the pool of Clp protein chaperones in B. subtilis, presumably due to negative feedback regulation. Breaking this feedback can improve enzyme production yields. Considering the conserved nature of Clp chaperones and their regulation, this method might benefit high-yield enzyme production in other organisms.

Defining early steps in Bacillus subtilis biofilm biosynthesis

IMPORTANCE

Biofilms are the communal way of life that microbes adopt to increase survival. Key to our ability to systematically promote or ablate biofilm formation is a detailed understanding of the biofilm matrix macromolecules. Here, we identify the first two essential steps in the Bacillus subtilis biofilm matrix exopolysaccharide (EPS) synthesis pathway. Together, our studies and approaches provide the foundation for the sequential characterization of the steps in EPS biosynthesis, using prior steps to enable chemoenzymatic synthesis of the undecaprenyl diphosphate-linked glycan substrates.

Trace impurities in sodium phosphate influences the physiological activity of Escherichia coli in M9 minimal medium

In the field of applied microbiology, reproducibility and experimental variability are important factors that influence both basic research as well as process development for industrial applications. Experimental reproducibility and accuracy depend not only on culture conditions such as temperature and aeration but also on raw materials and procedures used for media preparation. The M9 minimal medium is one of the most common synthetic media for culturing Escherichia coli and other bacteria. This synthetic medium can be used to observe and evaluate the physiological activity of microbes under minimal nutritional requirements and determine the limiting factor for the desired phenotype. Although one of the advantages using the M9 medium is that its composition can be modulated, it is difficult to control presence of trace components and impurities from the reagents for preparing this medium. Herein, we showed that trace ingredients present in the reagents used for M9 media preparation affect the bacterial physiological activities (e.g., cell growth, substrate consumption, and byproduct formation). Additionally, we systematically identified the trace ingredient that influenced phenotypic differences. Our results showed that the selection of reagents and accuracy during reagent preparation is important for experimental reproducibility in the field of bio-engineering and systems biology focused on the systematic and continuous development of biomolecular systems (e.g., biorefinery, metabolic engineering, and synthetic biology).

Pleiotropic hubs drive bacterial surface competition through parallel changes in colony composition and expansion

Bacteria commonly adhere to surfaces where they compete for both space and resources. Despite the importance of surface growth, it remains largely elusive how bacteria evolve on surfaces. We previously performed an evolution experiment where we evolved distinct Bacilli populations under a selective regime that favored colony spreading. In just a few weeks, colonies of Bacillus subtilis showed strongly advanced expansion rates, increasing their radius 2.5-fold relative to that of the ancestor. Here, we investigate what drives their rapid evolution by performing a uniquely detailed analysis of the evolutionary changes in colony development. We find mutations in diverse global regulators, RicT, RNAse Y, and LexA, with strikingly similar pleiotropic effects: They lower the rate of sporulation and simultaneously facilitate colony expansion by either reducing extracellular polysaccharide production or by promoting filamentous growth. Combining both high-throughput flow cytometry and gene expression profiling, we show that regulatory mutations lead to highly reproducible and parallel changes in global gene expression, affecting approximately 45% of all genes. This parallelism results from the coordinated manner by which regulators change activity both during colony development-in the transition from vegetative growth to dormancy-and over evolutionary time. This coordinated activity can however also break down, leading to evolutionary divergence. Altogether, we show how global regulators function as major pleiotropic hubs that drive rapid surface adaptation by mediating parallel changes in both colony composition and expansion, thereby massively reshaping gene expression.

Major distinctions between the two oligopeptide permease systems of Bacillus subtilis with respect to signaling, development and evolutionary divergence

Oligopeptide-permeases (Opps) are used by bacteria to import short peptides. In addition to their metabolic benefit, imported short peptides are used in many Gram-positive bacteria as signalling molecules of the RRNPP super-family of quorum-sensing systems, making Opps an integral part of cell–cell communication. In some Gram-positive bacteria there exist multiple Opps and the relative importance of those to RRNPP quorum sensing are not fully clear. Specifically, in Bacillus subtilis , the Gram-positive model species, there exist two homologous oligopeptide permeases named Opp and App. Previous work showed that the App system is mutated in lab strain 168 and its recovery partially complements an Opp mutation for several developmental processes. Yet, the nature of the impact of App on signalling and development in wild-type strains, where both permeases are active was not studied. Here we re-examine the impact of the two permease systems. We find that App has a minor contribution to biofilm formation, surfactin production and phage infection compared to the effect of Opp. This reduced effect is also reflected in its lower ability to import the signals of four different Rap-Phr RRNPP systems. Further analysis of the App system revealed that, unlike Opp, some App genes have undergone horizontal transfer, resulting in two distinct divergent alleles of this system in B. subtilis strains. We found that both alleles were substantially better adapted than the Opp system to import an exogenous RRNPP signal of the Bacillus cereus group PlcR-PapR system. In summary, we find that the App system has only a minor role in signalling but may still be crucial for the import of other peptides.

Tuning transcription factor DegU for developing extracellular protease overproducer in Bacillus pumilus

BACKGROUND

Global transcription machinery engineering (gTME) is an effective approach employed in strain engineering to rewire gene expression and reshape cellular metabolic fluxes at the transcriptional level.

RESULTS

In this study, we utilized gTME to engineer the positive transcription factor, DegU, in the regulation network of major alkaline protease, AprE, in Bacillus pumilus. To validate its functionality when incorporated into the chromosome, we performed several experiments. First, three negative transcription factors, SinR, Hpr, and AbrB, were deleted to promote AprE synthesis. Second, several hyper-active DegU mutants, designated as DegU(hy), were selected using the fluorescence colorimetric method with the host of the Bacillus subtilis ΔdegSU mutant. Third, we integrated a screened degU(L113F) sequence into the chromosome of the Δhpr mutant of B. pumilus SCU11 to replace the original degU gene using a CRISPR/Cas9 system. Finally, based on transcriptomic and molecular dynamic analysis, we interpreted the possible mechanism of high-yielding and found that the strain produced alkaline proteases 2.7 times higher than that of the control strain (B. pumilus SCU11) in LB medium.

CONCLUSION

Our findings serve as a proof-of-concept that tuning the global regulator is feasible and crucial for improving the production performance of B. pumilus. Additionally, our study established a paradigm for gene function research in strains that are difficult to handle.

The genome of Bacillus tequilensis EA-CB0015 sheds light into its epiphytic lifestyle and potential as a biocontrol agent

Different Bacillus species have successfully been used as biopesticides against a broad range of plant pathogens. Among these, Bacillus tequilensis EA-CB0015 has shown to efficiently control Black sigatoka disease in banana plants, presumably by mechanisms of adaptation that involve modifying the phyllosphere environment. Here, we report the complete genome of strain EA-CB0015, its precise taxonomic identity, and determined key genetic features that may contribute to its effective biocontrol of plant pathogens. We found that B. tequilensis EA-CB0015 harbors a singular 4 Mb circular chromosome, with 3,951 protein-coding sequences. Multi-locus sequence analysis (MLSA) and average nucleotide identity (ANI) analysis classified strain EA-CB0015 as B. tequilensis. Encoded within its genome are biosynthetic gene clusters (BGCs) for surfactin, iturin, plipastatin, bacillibactin, bacilysin, subtilosin A, sporulation killing factor, and other natural products that may facilitate inter-microbial warfare. Genes for indole-acetic acid (IAA) synthesis, the use of diverse carbon sources, and a multicellular lifestyle involving motility, biofilm formation, quorum sensing, competence, and sporulation suggest EA-CB0015 is adept at colonizing plant surfaces. Defensive mechanisms to survive invading viral infections and preserve genome integrity include putative type I and type II restriction modification (RM) and toxin/antitoxin (TA) systems. The presence of bacteriophage sequences, genomic islands, transposable elements, virulence factors, and antibiotic resistance genes indicate prior occurrences of genetic exchange. Altogether, the genome of EA-CB0015 supports its function as a biocontrol agent against phytopathogens and suggest it has adapted to thrive within phyllosphere environments.

Defining Early Steps in B. subtilis Biofilm Biosynthesis

The Bacillus subtilis extracellular biofilm matrix includes an exopolysaccharide that is critical for the architecture and function of the community. To date, our understanding of the biosynthetic machinery and the molecular composition of the exopolysaccharide of B. subtilis remains unclear and incomplete. This report presents synergistic biochemical and genetic studies built from a foundation of comparative sequence analyses targeted at elucidating the activities of the first two membrane-committed steps in the exopolysaccharide biosynthetic pathway. By taking this approach, we determined the nucleotide sugar donor and lipid-linked acceptor substrates for the first two enzymes in the B. subtilis biofilm exopolysaccharide biosynthetic pathway. EpsL catalyzes the first phosphoglycosyl transferase step using UDP-di- N -acetyl bacillosamine as phospho-sugar donor. EpsD is a GT-B fold glycosyl transferase that facilitates the second step in the pathway that utilizes the product of EpsL as an acceptor substrate and UDP- N -acetyl glucosamine as the sugar donor. Thus, the study defines the first two monosaccharides at the reducing end of the growing exopolysaccharide unit. In doing so we provide the first evidence of the presence of bacillosamine in an exopolysaccharide synthesized by a Gram-positive bacterium.

IMPORTANCE

Biofilms are the communal way of life that microbes adopt to increase survival. Key to our ability to systematically promote or ablate biofilm formation is a detailed understanding of the biofilm matrix macromolecules. Here we identify the first two essential steps in the Bacillus subtilis biofilm matrix exopolysaccharide synthesis pathway. Together our studies and approaches provide the foundation for the sequential characterization of the steps in exopolysaccharide biosynthesis, using prior steps to enable chemoenzymatic synthesis of the undecaprenol diphosphate-linked glycan substrates.

SigB modulates expression of novel SigB regulon members via Bc1009 in non-stressed and heat-stressed cells revealing its alternative roles in Bacillus cereus

BACKGROUND

The Bacillus cereus Sigma B (SigB) dependent general stress response is activated via the two-component RsbKY system, which involves a phosphate transfer from RsbK to RsbY. It has been hypothesized that the Hpr-like phosphocarrier protein (Bc1009) encoded by bc1009 in the SigB gene cluster may play a role in this transfer, thereby acting as a regulator of SigB activation. Alternatively, Bc1009 may be involved in the activation of a subset of SigB regulon members.

RESULTS

We first investigated the potential role of bc1009 to act as a SigB regulator but ruled out this possibility as the deletion of bc1009 did not affect the expression of sigB and other SigB gene cluster members. The SigB-dependent functions of Bc1009 were further examined in B. cereus ATCC14579 via comparative proteome profiling (backed up by transcriptomics) of wt, Δbc1009 and ΔsigB deletion mutants under heat stress at 42 °C. This revealed 284 proteins displaying SigB-dependent alterations in protein expression levels in heat-stressed cells, including a subgroup of 138 proteins for which alterations were also Bc1009-dependent. Next to proteins with roles in stress defense, newly identified SigB and Bc1009-dependent proteins have roles in cell motility, signal transduction, transcription, cell wall biogenesis, and amino acid transport and metabolism. Analysis of lethal stress survival at 50 °C after pre-adaptation at 42 °C showed intermediate survival efficacy of Δbc1009 cells, highest survival of wt, and lowest survival of ΔsigB cells, respectively. Additional comparative proteome analysis of non-stressed wt and mutant cells at 30 °C revealed 96 proteins with SigB and Bc1009-dependent differences in levels: 51 were also identified under heat stress, and 45 showed significant differential expression at 30 °C. This includes proteins with roles in carbohydrate/ion transport and metabolism. Overlapping functions at 30 °C and 42 °C included proteins involved in motility, and ΔsigB and Δbc1009 cells showed reduced motility compared to wt cells in swimming assays at both temperatures.

CONCLUSION

Our results extend the B. cereus SigB regulon to > 300 members, with a novel role of SigB-dependent Bc1009 in the activation of a subregulon of  > 180 members, conceivably via interactions with other transcriptional regulatory networks.

Termination factor Rho mediates transcriptional reprogramming of Bacillus subtilis stationary phase

Transcription termination factor Rho is known for its ubiquitous role in suppression of pervasive, mostly antisense, transcription. In the model Gram-positive bacterium Bacillus subtilis, de-repression of pervasive transcription by inactivation of rho revealed the role of Rho in the regulation of post-exponential differentiation programs. To identify other aspects of the regulatory role of Rho during adaptation to starvation, we have constructed a B. subtilis strain (Rho+) that expresses rho at a relatively stable high level in order to compensate for its decrease in the wild-type cells entering stationary phase. The RNAseq analysis of Rho+, WT and Δrho strains (expression profiles can be visualized at http://genoscapist.migale.inrae.fr/seb_rho/) shows that Rho over-production enhances the termination efficiency of Rho-sensitive terminators, thus reducing transcriptional read-through and antisense transcription genome-wide. Moreover, the Rho+ strain exhibits global alterations of sense transcription with the most significant changes observed for the AbrB, CodY, and stringent response regulons, forming the pathways governing the transition to stationary phase. Subsequent physiological analyses demonstrated that maintaining rho expression at a stable elevated level modifies stationary phase-specific physiology of B. subtilis cells, weakens stringent response, and thereby negatively affects the cellular adaptation to nutrient limitations and other stresses, and blocks the development of genetic competence and sporulation. These results highlight the Rho-specific termination of transcription as a novel element controlling stationary phase. The release of this control by decreasing Rho levels during the transition to stationary phase appears crucial for the functionality of complex gene networks ensuring B. subtilis survival in stationary phase.

Resonant learning in scale-free networks

Large networks of interconnected components, such as genes or machines, can coordinate complex behavioral dynamics. One outstanding question has been to identify the design principles that allow such networks to learn new behaviors. Here, we use Boolean networks as prototypes to demonstrate how periodic activation of network hubs provides a network-level advantage in evolutionary learning. Surprisingly, we find that a network can simultaneously learn distinct target functions upon distinct hub oscillations. We term this emergent property resonant learning, as the new selected dynamical behaviors depend on the choice of the period of the hub oscillations. Furthermore, this procedure accelerates the learning of new behaviors by an order of magnitude faster than without oscillations. While it is well-established that modular network architecture can be selected through evolutionary learning to produce different network behaviors, forced hub oscillations emerge as an alternative evolutionary learning strategy for which network modularity is not necessarily required.

Prediction and validation of novel SigB regulon members in Bacillus subtilis and regulon structure comparison to Bacillales members

BACKGROUND

Sigma factor B (SigB) is the central regulator of the general stress response in Bacillus subtilis and regulates a group of genes in response to various stressors, known as the SigB regulon members. Genes that are directly regulated by SigB contain a promotor binding motif (PBM) with a previously identified consensus sequence.

RESULTS

In this study, refined SigB PBMs were derived and different spacer compositions and lengths (N₁₂-N₁₇) were taken into account. These were used to identify putative SigB-regulated genes in the B. subtilis genome, revealing 255 genes: 99 had been described in the literature and 156 genes were newly identified, increasing the number of SigB putative regulon members (with and without a SigB PBM) to > 500 in B. subtilis. The 255 genes were assigned to five categories (I-V) based on their similarity to the original SigB consensus sequences. The functionalities of selected representatives per category were assessed using promoter-reporter fusions in wt and ΔsigB mutants upon exposure to heat, ethanol, and salt stress. The activity of the P_rsbV (I) positive control was induced upon exposure to all three stressors. P_ytoQ (II) showed SigB-dependent activity only upon exposure to ethanol, whereas P_pucI (II) with a N₁₇ spacer and P_ylaL (III) with a N₁₆ spacer showed mild induction regardless of heat/ethanol/salt stress. P_ywzA (III) and P_yaaI (IV) displayed ethanol-specific SigB-dependent activities despite a lower-level conserved - 10 binding motif. P_gtaB (V) was SigB-induced under ethanol and salt stress while lacking a conserved - 10 binding region. The activities of P_ygaO and P_ykaA (III) did not show evident changes under the conditions tested despite having a SigB PBM that highly resembled the consensus. The identified extended SigB regulon candidates in B. subtilis are mainly involved in coping with stress but are also engaged in other cellular processes. Orthologs of SigB regulon candidates with SigB PBMs were identified in other Bacillales genomes, but not all showed a SigB PBM. Additionally, genes involved in the integration of stress signals to activate SigB were predicted in these genomes, indicating that SigB signaling and regulon genes are species-specific.

CONCLUSION

The entire SigB regulatory network is sophisticated and not yet fully understood even for the well-characterized organism B. subtilis 168. Knowledge and information gained in this study can be used in further SigB studies to uncover a complete picture of the role of SigB in B. subtilis and other species.

Bacillus subtilis forms twisted cells with cell wall integrity defects upon removal of the molecular chaperones DnaK and trigger factor

The protein homeostasis network ensures a proper balance between synthesis, folding, and degradation of all cellular proteins. DnaK and trigger factor (TF) are ubiquitous bacterial molecular chaperones that assist in protein folding, as well as preventing protein misfolding and aggregation. In Escherichia coli, DnaK and TF possess partially overlapping functions. Their combined depletion results in proteostasis collapse and is synthetically lethal at temperatures above 30°C. To increase our understanding on how proteostasis is maintained in Gram-positive bacteria, we have investigated the physiological effects of deleting dnaK and tig (encoding for DnaK and TF) in Bacillus subtilis. We show that combined deletion of dnaK and tig in B. subtilis is non-lethal, but causes a severe pleiotropic phenotype, including an aberrant twisted and filamentous cell morphology, as well as decreased tolerance to heat and to cell wall active antibiotics and hydrolytic enzymes, indicative of defects in cell wall integrity. In addition, cells lacking DnaK and TF have a much smaller colony size due to defects in motility. Despite these physiological changes, we observed no major compromises in important cellular processes such as cell growth, FtsZ localization and division and only moderate defects in spore formation. Finally, through suppressor analyses, we found that the wild-type cell shape can be partially restored by mutations in genes involved in metabolism or in other diverse cellular processes.

PRO-Simat: Protein network simulation and design tool

PRO-Simat is a simulation tool for analysing protein interaction networks, their dynamic change and pathway engineering. It provides GO enrichment, KEGG pathway analyses, and network visualisation from an integrated database of more than 8 million protein-protein interactions across 32 model organisms and the human proteome. We integrated dynamical network simulation using the Jimena framework, which quickly and efficiently simulates Boolean genetic regulatory networks. It enables simulation outputs with in-depth analysis of the type, strength, duration and pathway of the protein interactions on the website. Furthermore, the user can efficiently edit and analyse the effect of network modifications and engineering experiments. In case studies, applications of PRO-Simat are demonstrated: (i) understanding mutually exclusive differentiation pathways in Bacillus subtilis, (ii) making Vaccinia virus oncolytic by switching on its viral replication mainly in cancer cells and triggering cancer cell apoptosis and (iii) optogenetic control of nucleotide processing protein networks to operate DNA storage. Multilevel communication between components is critical for efficient network switching, as demonstrated by a general census on prokaryotic and eukaryotic networks and comparing design with synthetic networks using PRO-Simat. The tool is available at https://prosimat.heinzelab.de/ as a web-based query server.

Predicting stress response and improved protein overproduction in Bacillus subtilis

Bacillus subtilis is a well-characterized microorganism and a model for the study of Gram-positive bacteria. The bacterium can produce proteins at high densities and yields, which has made it valuable for industrial bioproduction. Like other cell factories, metabolic modeling of B. subtilis has discovered ways to optimize its metabolism toward various applications. The first genome-scale metabolic model (M-model) of B. subtilis was published more than a decade ago and has been applied extensively to understand metabolism, to predict growth phenotypes, and served as a template to reconstruct models for other Gram-positive bacteria. However, M-models are ill-suited to simulate the production and secretion of proteins as well as their proteomic response to stress. Thus, a new generation of metabolic models, known as metabolism and gene expression models (ME-models), has been initiated. Here, we describe the reconstruction and validation of a ME model of B. subtilis, iJT964-ME. This model achieved higher performance scores on the prediction of gene essentiality as compared to the M-model. We successfully validated the model by integrating physiological and omics data associated with gene expression responses to ethanol and salt stress. The model further identified the mechanism by which tryptophan synthesis is upregulated under ethanol stress. Further, we employed iJT964-ME to predict amylase production rates under two different growth conditions. We analyzed these flux distributions and identified key metabolic pathways that permitted the increase in amylase production. Models like iJT964-ME enable the study of proteomic response to stress and the illustrate the potential for optimizing protein production in bacteria.

The B. subtilis Rok protein is an atypical H-NS-like protein irresponsive to physico-chemical cues

Nucleoid-associated proteins (NAPs) play a central role in chromosome organization and environment-responsive transcription regulation. The Bacillus subtilis-encoded NAP Rok binds preferentially AT-rich regions of the genome, which often contain genes of foreign origin that are silenced by Rok binding. Additionally, Rok plays a role in chromosome architecture by binding in genomic clusters and promoting chromosomal loop formation. Based on this, Rok was proposed to be a functional homolog of E. coli H-NS. However, it is largely unclear how Rok binds DNA, how it represses transcription and whether Rok mediates environment-responsive gene regulation. Here, we investigated Rok's DNA binding properties and the effects of physico-chemical conditions thereon. We demonstrate that Rok is a DNA bridging protein similar to prototypical H-NS-like proteins. However, unlike these proteins, the DNA bridging ability of Rok is not affected by changes in physico-chemical conditions. The DNA binding properties of the Rok interaction partner sRok are affected by salt concentration. This suggests that in a minority of Bacillus strains Rok activity can be modulated by sRok, and thus respond indirectly to environmental stimuli. Despite several functional similarities, the absence of a direct response to physico-chemical changes establishes Rok as disparate member of the H-NS family.

Genome Mining Reveals High Biosynthetic Potential of Biocontrol Agent Bacillus velezensis B.BV10

The present study demonstrates the biocontrol potential of a plant growth-promoting bacterial strain using three different approaches: (i) an in vitro evaluation of antagonistic activity against important phytopathogenic fungi; (ii) an evaluation under greenhouse conditions with strawberry plants to assess the control of gray mold; and (iii) an in silico whole genome sequence mining to assign genetic features such as gene clusters or isolated genes to the strain activity. The in vitro assay showed that the B.BV10 strain presented antagonistic activity, inhibiting the mycelial growth in all the phytopathogenic fungi evaluated. The application of the Bacillus velezensis strain B.BV10 under greenhouse conditions reduced the presence of Botrytis cinerea and increased the mean fruit biomass. The genome of B.BV10 was estimated at 3,917,533 bp, with a GC content of 46.6% and 4088 coding DNA sequences, and was identified as B. velezensis. Biosynthetic gene clusters related to the synthesis of the molecules with antifungal activity were found in its genome. Genes related to the regulation/formation of biofilms, motility, and the important properties for the rhizospheric colonization were also found in the genome. The current study offers a comprehensive understanding of the genomic architecture and control activity of phytopathogenic fungi by the B. velezensis strain B.BV10 that may substantiate the industrialization of this strain in the future.

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.

Phage-Encoded Sigma Factors Alter Bacterial Dormancy

By entering a reversible state of reduced metabolic activity, dormant microorganisms are able to tolerate suboptimal conditions that would otherwise reduce their fitness. Dormancy may also benefit bacteria by serving as a refuge from parasitic infections. Here, we focus on dormancy in the Bacillota, where endospore development is transcriptionally regulated by the expression of sigma factors. A disruption of this process could influence the survivorship or reproduction of phages that infect spore-forming hosts with implications for coevolutionary dynamics. We characterized the distribution of sigma factors in over 4,000 genomes of diverse phages capable of infecting hosts that span the bacterial domain. From this, we identified homologs of sporulation-specific sigma factors in phages that infect spore-forming hosts. Unlike sigma factors required for phage reproduction, we provide evidence that sporulation-like sigma factors are nonessential for lytic infection. However, when expressed in the spore-forming Bacillus subtilis, some of these phage-derived sigma factors can activate the bacterial sporulation gene network and lead to a reduction in spore yield. Our findings suggest that the acquisition of host-like transcriptional regulators may allow phages to manipulate a complex and ancient trait in one of the most abundant cell types on Earth. IMPORTANCE As obligate parasites, phages exert strong top-down pressure on host populations with eco-evolutionary implications for community dynamics and ecosystem functioning. The process of phage infection, however, is constrained by bottom-up processes that influence the energetic and nutritional status of susceptible hosts. Many phages have acquired auxiliary genes from bacteria, which can be used to exploit host metabolism with consequences for phage fitness. In this study, we demonstrate that phages infecting spore-forming bacteria carry homologs of sigma factors, which their hosts use to orchestrate gene expression during spore development. By tapping into regulatory gene networks, phages may manipulate the physiology and survival strategies of nongrowing bacteria in ways that influence host-parasite coevolution.

A Comparative Analysis of the Core Proteomes within and among the Bacillus subtilis and Bacillus cereus Evolutionary Groups Reveals the Patterns of Lineage- and Species-Specific Adaptations

By integrating phylogenomic and comparative analyses of 1104 high-quality genome sequences, we identify the core proteins and the lineage-specific fingerprint proteins of the various evolutionary clusters (clades/groups/species) of the Bacillus genus. As fingerprints, we denote those core proteins of a certain lineage that are present only in that particular lineage and absent in any other Bacillus lineage. Thus, these lineage-specific fingerprints are expected to be involved in particular adaptations of that lineage. Intriguingly, with a few notable exceptions, the majority of the Bacillus species demonstrate a rather low number of species-specific fingerprints, with the majority of them being of unknown function. Therefore, species-specific adaptations are mostly attributed to highly unstable (in evolutionary terms) accessory proteomes and possibly to changes at the gene regulation level. A series of comparative analyses consistently demonstrated that the progenitor of the Cereus Clade underwent an extensive genomic expansion of chromosomal protein-coding genes. In addition, the majority (76–82%) of the B. subtilis proteins that are essential or play a significant role in sporulation have close homologs in most species of both the Subtilis and the Cereus Clades. Finally, the identification of lineage-specific fingerprints by this study may allow for the future development of highly specific vaccines, therapeutic molecules, or rapid and low-cost molecular tests for species identification.

Synthesis and mechanism-of-action of a novel synthetic antibiotic based on a dendritic system with bow-tie topology

Over the course of the last decades, the continuous exposure of bacteria to antibiotics-at least in parts due to misprescription, misuse, and misdosing-has led to the widespread development of antimicrobial resistances. This development poses a threat to the available medication in losing their effectiveness in treating bacterial infections. On the drug development side, only minor advances have been made to bring forward novel therapeutics. In addition to increasing the efforts and approaches of tapping the natural sources of new antibiotics, synthetic approaches to developing novel antimicrobials are being pursued. In this study, BDTL049 was rationally designed using knowledge based on the properties of natural antibiotics. BDTL049 is a carbosilane dendritic system with bow-tie type topology, which has antimicrobial activity at concentrations comparable to clinically established natural antibiotics. In this report, we describe its mechanism of action on the Gram-positive model organism Bacillus subtilis. Exposure to BDTL049 resulted in a complex transcriptional response, which pointed toward disturbance of the cell envelope homeostasis accompanied by disruption of other central cellular processes of bacterial metabolism as the primary targets of BDTL049 treatment. By applying a combination of whole-cell biosensors, molecular staining, and voltage sensitive dyes, we demonstrate that the mode of action of BDTL049 comprises membrane depolarization concomitant with pore formation. As a result, this new molecule kills Gram-positive bacteria within minutes. Since BDTL049 attacks bacterial cells at different targets simultaneously, this might decrease the chances for the development of bacterial resistances, thereby making it a promising candidate for a future antimicrobial agent.

The impact of PrsA over-expression on the Bacillus subtilis transcriptome during fed-batch fermentation of alpha-amylase production

The production of the alpha-amylase (AMY) enzyme in Bacillus subtilis at a high rate leads to the accumulation of unfolded AMY, which causes secretion stress. The over-expression of the PrsA chaperone aids enzyme folding and reduces stress. To identify affected pathways and potential mechanisms involved in the reduced growth, we analyzed the transcriptomic differences during fed-batch fermentation between a PrsA over-expressing strain and control in a time-series RNA-seq experiment. We observe transcription in 542 unannotated regions, of which 234 had significant changes in expression levels between the samples. Moreover, 1,791 protein-coding sequences, 80 non-coding genes, and 20 riboswitches overlapping UTR regions of coding genes had significant changes in expression. We identified putatively regulated biological processes via gene-set over-representation analysis of the differentially expressed genes; overall, the analysis suggests that the PrsA over-expression affects ATP biosynthesis activity, amino acid metabolism, and cell wall stability. The investigation of the protein interaction network points to a potential impact on cell motility signaling. We discuss the impact of these highlighted mechanisms for reducing secretion stress or detrimental aspects of PrsA over-expression during AMY production.

Experimental evolution of Bacillus subtilis on Arabidopsis thaliana roots reveals fast adaptation and improved root colonization

Bacillus subtilis is known to promote plant growth and protect plants against disease. B. subtilis rapidly adapts to Arabidopsis thaliana root colonization, as evidenced by improved root colonizers already after 12 consecutive transfers between seedlings in a hydroponic setup. Re-sequencing of single evolved isolates and endpoint populations revealed mutations in genes related to different bacterial traits, in accordance with evolved isolates displaying increased root colonization associated with robust biofilm formation in response to the plant polysaccharide xylan and impaired motility. Interestingly, evolved isolates suffered a fitness disadvantage in a non-selective environment, demonstrating an evolutionary cost of adaptation to the plant root. Finally, increased root colonization by an evolved isolate was also demonstrated in the presence of resident soil microbes. Our findings highlight how a plant growth-promoting rhizobacterium rapidly adapts to an ecologically relevant environment and reveal evolutionary consequences that are fundamental to consider when evolving strains for biocontrol purposes.

•

Bacillus subtilis shows fast adaptation to Arabidopsis thaliana roots in a hydroponic setup

•

Evolved isolates exhibit robust biofilms in response to xylan and impaired motility

•

Adaptation to A. thaliana roots is accompanied by an evolutionary cost

•

An evolved isolate shows higher root colonization in the presence of soil bacteria

Effects of 4-Br-A23187 on Bacillus subtilis cells and unilamellar vesicles reveal it to be a potent copper ionophore

Ionophores are small molecules or peptides that transport metal ions across biological membranes. Their transport capabilities are typically characterized in vitro using vesicles and single ion species. It is difficult to infer from these data which effects ionophores have on living cells in a complex environment (e.g. culture medium), since net ion movement is influenced by many factors including ion composition of the medium, concentration gradients, pH gradient, and protein-mediated transport processes across the membrane. To gain insights into the antibacterial mechanism of action of the semisynthetic polyether ionophore 4-Br-A23187, known to efficiently transport zinc and manganese in vitro, we investigated its effects on the gram-positive model organism Bacillus subtilis. In addition to monitoring cellular ion concentrations, the physiological impact of treatment was assessed on the proteome level. 4-Br-A23187 treatment resulted in an increase in intracellular copper levels, the extent of which depended on the copper concentration of the medium. Effects of copper accumulation mirrored by the proteomic response included oxidative stress, disturbance of proteostasis, metal and sulfur homeostasis. The antibiotic effect of 4-Br-A23187 is further aggravated by a decrease in intracellular manganese and magnesium. A liposome model confirmed that 4-Br-A23187 acts as copper ionophore in vitro. This article is protected by copyright. All rights reserved.

The Acetyltransferase RibT From Bacillus subtilis Affects in vivo Dynamics of the Multimeric Heavy Riboflavin Synthase Complex

Flavins are ubiquitous molecules in life as they serve as important enzyme cofactors. In the Gram-positive, soil-dwelling bacterium Bacillus subtilis, four well-characterized gene products (the enzymes RibDG, RibE, RibAB, and RibH) catalyze the biosynthesis of riboflavin (RF) from guanosine-triphosphate (GTP) and ribulose-5-phosphate (R5P). The corresponding genes form an operon together with the gene ribT (ribDG-E-AB-H-T), wherein the function of this terminal gene remained enigmatic. RibT has been structurally characterized as a GCN5-like acetyltransferase (GNAT), however, with unidentified target molecules. Bacterial two-hybrid system revealed interactions between RibT, RibH, and RibE, forming the heavy RF synthase complex. Applying single particle tracking (SPT), we found that confined (sub)diffusion of RibT is largely dependent on interacting RibE and, to a lesser degree, on interacting RibH. By induced expression of otherwise low-expressed ribT from an ectopic locus, we observed a decrease in the subpopulation considered to represent capsids of the heavy RF synthase and an increase in the subpopulation thought to represent pentamers of RibH, pointing to a putative role for RibT in capsid disassembly. Complementarily, either deletion of ribT or mutation of a key residue from RibH (K29) suspected to be the substrate of RibT for acetylation leads to increased levels of subpopulations considered as capsids of RibH-mVenus (RibH-mV) in comparison to wild-type (wt)-like cells. Thus, we provide evidence for an indirect involvement of RibT in RF biosynthesis by a putative capsid disassembling mechanism considered to involve acetylation of RibH residue K29 at the three-fold symmetry axis of 60-mer capsids.

Surfactin Shows Relatively Low Antimicrobial Activity against Bacillus subtilis and Other Bacterial Model Organisms in the Absence of Synergistic Metabolites

Surfactin is described as a powerful biosurfactant and is natively produced by Bacillus subtilis in notable quantities. Among other industrially relevant characteristics, antimicrobial properties have been attributed to surfactin-producing Bacillus isolates. To investigate this property, stress approaches were carried out with biotechnologically established strains of Corynebacterium glutamicum, Bacillus subtilis, Escherichia coli and Pseudomonas putida with the highest possible amounts of surfactin. Contrary to the popular opinion, the highest growth-reducing effects were detectable in B. subtilis and E. coli after surfactin treatment of 100 g/L with 35 and 33%, respectively, while P. putida showed no growth-specific response. In contrast, other antimicrobial biosurfactants, like rhamnolipids and sophorolipids, showed significantly stronger effects on bacterial growth. Since the addition of high amounts of surfactin in defined mineral salt medium reduced the cell growth of B. subtilis by about 40%, the initial stress response at the protein level was analyzed by mass spectrometry, showing induction of stress proteins under control of alternative sigma factors σB and σW as well as the activation of LiaRS two-component system. Overall, although surfactin is associated with antimicrobial properties, relatively low growth-reducing effects could be demonstrated after the surfactin addition, challenging the general claim of the antimicrobial properties of surfactin.

High-Throughput Time-Lapse Fluorescence Microscopy Screening for Heterogeneously Expressed Genes in Bacillus subtilis

Elucidating phenotypic heterogeneity in clonal bacterial populations is important for both the fundamental understanding of bacterial behavior and the synthetic engineering of bacteria in biotechnology. In this study, we present and validate a high-throughput and high-resolution time-lapse fluorescence microscopy-based strategy to easily and systematically screen for heterogeneously expressed genes in the Bacillus subtilis model bacterium. This screen allows detection of expression patterns at high spatial and temporal resolution, which often escape detection by other approaches, and can readily be extrapolated to other bacteria. A proof-of-concept screening in B. subtilis revealed both recognized and yet unrecognized heterogeneously expressed genes, thereby validating the approach.

IMPORTANCE Differential gene expression among isogenic siblings often leads to phenotypic heterogeneity and the emergence of complex social behavior and functional capacities within clonal bacterial populations. Despite the importance of such features for both the fundamental understanding and synthetic engineering of bacterial behavior, approaches to systematically map such population heterogeneity are scarce. In this context, we have elaborated a new time-lapse fluorescence microscopy-based strategy to easily and systematically screen for such heterogeneously expressed genes in bacteria with high resolution and throughput. A proof-of-concept screening in the Bacillus subtilis model bacterium revealed both recognized and yet unrecognized heterogeneously expressed genes, thereby validating our approach.

Discovery and initial characterization of YloC, a novel endoribonuclease in Bacillus subtilis

The Bacillus subtilis genome is predicted to encode numerous ribonucleases, including four 3' exoribonucleases that have been characterized to some extent. A strain containing gene knockouts of all four known 3' exoribonucleases is viable, suggesting that one or more additional RNases remain to be discovered. A protein extract from the quadruple RNase mutant strain was fractionated and RNase activity was followed, resulting in the identification of an enzyme activity catalyzed by the YloC protein. YloC is an endoribonuclease and is a member of the highly conserved "YicC family" of proteins that is widespread in bacteria. YloC is a metal-dependent enzyme that catalyzes the cleavage of single-stranded RNA, preferentially at U residues, and exists in an oligomeric form, most likely a hexamer. As such, YloC shares some characteristics with the SARS-CoV Nsp15 endoribonuclease. While the in vivo function of YloC in B. subtilis is yet to be determined, YloC was found to act similarly to YicC in an Escherichia coli in vivo assay that assesses decay of the small RNA, RyhB. Thus, YloC may play a role in small RNA regulation.

Prediction of whole-cell transcriptional response with machine learning

MOTIVATION

Applications in synthetic and systems biology can benefit from measuring whole-cell response to biochemical perturbations. Execution of experiments to cover all possible combinations of perturbations is infeasible. In this paper, we present the host response model (HRM), a machine learning approach that maps response of single perturbations to transcriptional response of the combination of perturbations.

RESULTS

The HRM combines high-throughput sequencing with machine learning to infer links between experimental context, prior knowledge of cell regulatory networks, and RNASeq data to predict a gene's dysregulation. We find that the HRM can predict the directionality of dysregulation to a combination of inducers with an accuracy of >90% using data from single inducers. We further find that the use of prior, known cell regulatory networks doubles the predictive performance of the HRM (an R2 from 0.3 to 0.65). The model was validated in two organisms, Escherichia coli and Bacillus subtilis, using new experiments conducted after training. Finally, while the HRM is trained with gene expression data, the direct prediction of differential expression makes it possible to also conduct enrichment analyses using its predictions. We show that the HRM can accurately classify >95% of the pathway regulations. The HRM reduces the number of RNASeq experiments needed as responses can be tested in silico prior to the experiment.

AVAILABILITY AND IMPLEMENTATION

The HRM software and tutorial are available at https://github.com/sd2e/CDM and the configurable differential expression analysis tools and tutorials are available at https://github.com/SD2E/omics_tools.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Redirected Stress Responses in a Genome-Minimized 'midiBacillus' Strain with Enhanced Capacity for Protein Secretion

Genome engineering offers the possibility to create completely novel cell factories with enhanced properties for biotechnological applications. In recent years, genome minimization was extensively explored in the Gram-positive bacterial cell factory Bacillus subtilis, where up to 42% of the genome encoding dispensable functions was removed. Such studies showed that some strains with minimized genomes gained beneficial features, especially for secretory protein production. However, strains with the most minimal genomes displayed growth defects. This focused our attention on strains with less extensive genomic deletions that display close-to-wild-type growth properties while retaining the acquired beneficial traits in secretory protein production. A strain of this category is B. subtilis IIG-Bs27-47-24, here referred to as midiBacillus, which lacks 30.95% of the parental genome. To date, it was unknown how the altered genomic configuration of midiBacillus impacts cell physiology in general, and protein secretion in particular. The present study bridges this knowledge gap through comparative quantitative proteome analyses with focus on protein secretion. Interestingly, the results show that the secretion stress responses of midiBacillus, as elicited by high-level expression of the immunodominant staphylococcal antigen A, are completely different from secretion stress responses that occur in the parental strain 168. We further show that midiBacillus has an increased capacity for translation and that a variety of critical Sec secretion machinery components is present at elevated levels. Altogether, our observations demonstrate that high-level protein secretion has different consequences for wild-type and genome-engineered Bacillus strains, dictated by the altered genomic and proteomic configurations.

IMPORTANCE Our present study showcases a genome-minimized nonpathogenic bacterium, the so-called midiBacillus, as a chassis for the development of future industrial strains that serve in the production of high-value difficult-to-produce proteins. In particular, we explain how midiBacillus, which lacks about one-third of the original genome, effectively secretes a protein of the major human pathogen Staphylococcus aureus that cannot be produced by the parental Bacillus subtilis strain. This is important, because the secreted S. aureus protein is exemplary for a range of targets that can be implemented in future antistaphylococcal immunotherapies. Accordingly, we anticipate that midiBacillus chassis will contribute to the development of vaccines that protect both humans and livestock against diseases caused by S. aureus, a bacterial pathogen that is increasingly difficult to fight with antibiotics, because it has accumulated resistances to essentially all antibiotics that are currently in clinical practice.

A High-Content Microscopy Screening Identifies New Genes Involved in Cell Width Control in Bacillus subtilis

How cells control their shape and size is a fundamental question of biology. In most bacteria, cell shape is imposed by the peptidoglycan (PG) polymeric meshwork that surrounds the cell. Thus, bacterial cell morphogenesis results from the coordinated action of the proteins assembling and degrading the PG shell. Remarkably, during steady-state growth, most bacteria maintain a defined shape along generations, suggesting that error-proof mechanisms tightly control the process. In the rod-shaped model for the Gram-positive bacterium Bacillus subtilis, the average cell length varies as a function of the growth rate, but the cell diameter remains constant throughout the cell cycle and across growth conditions. Here, in an attempt to shed light on the cellular circuits controlling bacterial cell width, we developed a screen to identify genetic determinants of cell width in B. subtilis. Using high-content screening (HCS) fluorescence microscopy and semiautomated measurement of single-cell dimensions, we screened a library of ∼4,000 single knockout mutants. We identified 13 mutations significantly altering cell diameter, in genes that belong to several functional groups. In particular, our results indicate that metabolism plays a major role in cell width control in B. subtilis.

IMPORTANCE Bacterial shape is primarily dictated by the external cell wall, a vital structure that, as such, is the target of countless antibiotics. Our understanding of how bacteria synthesize and maintain this structure is therefore a cardinal question for both basic and applied research. Bacteria usually multiply from generation to generation while maintaining their progenies with rigorously identical shapes. This implies that the bacterial cells constantly monitor and maintain a set of parameters to ensure this perpetuation. Here, our study uses a large-scale microscopy approach to identify at the whole-genome level, in a model bacterium, the genes involved in the control of one of the most tightly controlled cellular parameters, the cell width.

Glucolipids and lipoteichoic acids affect the activity of SigI, an alternative sigma factor, and WalKR, an essential two-component system, in Bacillus subtilis

In a Bacillus subtilis ugtP mutant lacking glucolipids, SigI was activated in the log phase, and the activation of SigI in the mutant was suppressed by the expression of native ugtP. By contrast, SigI was inhibited in a yfnI mutant lacking one of the lipoteichoic acid (LTA) synthase genes, and the inhibition was suppressed by the expression of yfnI. A series of mutation analyses of the sigI promoter revealed that the two WalR binding sites were involved in the increase of P_sigI -lacZ activity in the ugtP mutant and decrease of the lacZ activity in the yfnI mutant. Transcription from the SigI recognition sequence was enhanced in the ugtP mutant, whereas yfnI disruption inhibited the transcription from the SigA recognition sequence in the sigI promoter. We found that not only SigI but also WalKR, the essential two-component system, was activated in the ugtP mutant and inhibited in the yfnI mutant. The walK mutants with activated WalR exhibited abnormal morphology, but this phenotype was suppressed by the addition of MgSO₄ . We conclude that glucolipids and LTA are key compounds in the maintenance of normal cell surface structure in B. subtilis.

Identification and characterization of a non-canonical menaquinone-linked formate dehydrogenase

The Molybdenum/Tungsten-bispyranopterin guanine dinucleotides (Mo/W-bisPGD) family of Formate Dehydrogenases (FDHs) plays roles in several metabolic pathways ranging from carbon fixation to energy harvesting owing to their reaction with a wide variety of redox partners. Indeed, this metabolic plasticity results from the diverse structures, cofactor content, and substrates employed by partner subunits interacting with the catalytic hub. Here, we unveiled two non-canonical FDHs in Bacillus subtilis which are organized into two-subunit complexes with unique features, ForCE1 and ForCE2. We show that the ForC catalytic subunit interacts with an unprecedented partner subunit, ForE, and that its amino acid sequence within the active site deviates from the consensus residues typically associated with FDH activity, as a histidine residue is naturally substituted with a glutamine. The ForE essential subunit mediates the utilization of menaquinone as an electron acceptor as shown by the formate:menadione oxidoreductase activity of both enzymes, their copurification with menaquinone, and the distinctive detection of a protein-bound neutral menasemiquinone radical by multifrequency electron paramagnetic resonance (EPR) experiments on the purified enzymes. Moreover, EPR characterization of both FDHs reveals the presence of several [Fe-S] clusters with distinct relaxation properties and a weakly anisotropic Mo(V) EPR signature, consistent with the characteristic Mo/bisPGD cofactor of this enzyme family. Altogether, this work enlarges our knowledge of the FDH family by identifying a non-canonical FDH, which differs in terms of architecture, amino acid conservation around the Mo cofactor, and reactivity.

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.

An Organogold Compound as Potential Antimicrobial Agent against Drug-Resistant Bacteria: Initial Mechanistic Insights

The rise of antimicrobial resistance has necessitated novel strategies to efficiently combat pathogenic bacteria. Metal-based compounds have been proven as a possible alternative to classical organic drugs. Here, we have assessed the antibacterial activity of seven gold complexes of different families. One compound, a cyclometalated Au(III) C^N complex, showed activity against Gram-positive bacteria, including multi-drug resistant clinical strains. The mechanism of action of this compound was studied in Bacillus subtilis. Overall, the studies point towards a complex mode of antibacterial action, which does not include induction of oxidative stress or cell membrane damage. A number of genes related to metal transport and homeostasis were upregulated upon short treatment of the cells with gold compound. Toxicity tests conducted on precision-cut mouse tissue slices ex vivo revealed that the organogold compound is poorly toxic to mouse liver and kidney tissues, and may thus, be treated as an antibacterial drug candidate.

Inducible Population Quality Control of Engineered Bacillus subtilis for Improved N-Acetylneuraminic Acid Biosynthesis

Biosynthesis by microorganisms using renewable feedstocks is an important approach for realizing sustainable chemical manufacturing. However, cell-to-cell variation in biosynthesis capability during fermentation restricts the robustness and efficiency of bioproduction, hampering the industrialization of biosynthesis. Herein, we developed an inducible population quality control system (iPopQC) for dynamically modulating the producing and nonproducing subpopulations of engineered Bacillus subtilis, which was constructed via inducible promoter- and metabolite-responsive biosensor-based genetic circuit for regulating essential genes. Moreover, iPopQC achieved a 1.97-fold increase in N-acetylneuraminic acid (NeuAc) titer by enriching producing cell subpopulation during cultivation, representing 52% higher than that of previous PopQC. Strains with double-output iPopQC cocoupling the expression of double essential genes with NeuAc production improved production robustness further, retaining NeuAc production throughout 96 h of fermentation, upon which the strains cocoupling one essential gene expression with NeuAc production abolished the production ability.

Designing efficient genetic code expansion in Bacillus subtilis to gain biological insights

Bacillus subtilis is a model gram-positive bacterium, commonly used to explore questions across bacterial cell biology and for industrial uses. To enable greater understanding and control of proteins in B. subtilis, here we report broad and efficient genetic code expansion in B. subtilis by incorporating 20 distinct non-standard amino acids within proteins using 3 different families of genetic code expansion systems and two choices of codons. We use these systems to achieve click-labelling, photo-crosslinking, and translational titration. These tools allow us to demonstrate differences between E. coli and B. subtilis stop codon suppression, validate a predicted protein-protein binding interface, and begin to interrogate properties underlying bacterial cytokinesis by precisely modulating cell division dynamics in vivo. We expect that the establishment of this simple and easily accessible chemical biology system in B. subtilis will help uncover an abundance of biological insights and aid genetic code expansion in other organisms.

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.