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

Bacterial cell volume regulation and the importance of cyclic di-AMP

SUMMARYNucleotide-derived second messengers are present in all domains of life. In prokaryotes, most of their functionality is associated with general lifestyle and metabolic adaptations, often in response to environmental fluctuations of physical parameters. In the last two decades, cyclic di-AMP has emerged as an important signaling nucleotide in many prokaryotic lineages, including Firmicutes, Actinobacteria, and Cyanobacteria. Its importance is highlighted by the fact that both the lack and overproduction of cyclic di-AMP affect viability of prokaryotes that utilize cyclic di-AMP, and that it generates a strong innate immune response in eukaryotes. In bacteria that produce the second messenger, most molecular targets of cyclic di-AMP are associated with cell volume control. Besides, other evidence links the second messenger to cell wall remodeling, DNA damage repair, sporulation, central metabolism, and the regulation of glycogen turnover. In this review, we take a biochemical, quantitative approach to address the main cellular processes that are directly regulated by cyclic di-AMP and show that these processes are very connected and require regulation of a similar set of proteins to which cyclic di-AMP binds. Altogether, we argue that cyclic di-AMP is a master regulator of cell volume and that other cellular processes can be connected with cyclic di-AMP through this core function. We further highlight important directions in which the cyclic di-AMP field has to develop to gain a full understanding of the cyclic di-AMP signaling network and why some processes are regulated, while others are not.

Bacillus subtilis stress-associated mutagenesis and developmental DNA repair

SUMMARYThe metabolic conditions that prevail during bacterial growth have evolved with the faithful operation of repair systems that recognize and eliminate DNA lesions caused by intracellular and exogenous agents. This idea is supported by the low rate of spontaneous mutations (10-9) that occur in replicating cells, maintaining genome integrity. In contrast, when growth and/or replication cease, bacteria frequently process DNA lesions in an error-prone manner. DNA repairs provide cells with the tools needed for maintaining homeostasis during stressful conditions and depend on the developmental context in which repair events occur. Thus, different physiological scenarios can be anticipated. In nutritionally stressed bacteria, different components of the base excision repair pathway may process damaged DNA in an error-prone approach, promoting genetic variability. Interestingly, suppressing the mismatch repair machinery and activating specific DNA glycosylases promote stationary-phase mutations. Current evidence also suggests that in resting cells, coupling repair processes to actively transcribed genes may promote multiple genetic transactions that are advantageous for stressed cells. DNA repair during sporulation is of interest as a model to understand how transcriptional processes influence the formation of mutations in conditions where replication is halted. Current reports indicate that transcriptional coupling repair-dependent and -independent processes operate in differentiating cells to process spontaneous and induced DNA damage and that error-prone synthesis of DNA is involved in these events. These and other noncanonical ways of DNA repair that contribute to mutagenesis, survival, and evolution are reviewed in this manuscript.

Single-cell heterogeneity in ribosome content and the consequences for the growth laws

Previous work has suggested that the ribosome content of a cell is optimized to maximize growth given the nutrient availability. The resulting correlation between ribosome number and growth rate appears to be independent of the rate limiting nutrient and has been reported in many organisms. The robustness and universality of this observation has given it the classification of a "growth law." These laws have had powerful impacts on many biological disciplines. They have fueled predictions about how organisms evolve to maximize reproduction, and informed models about how cells regulate growth. Due to methodological limitations, this growth law has rarely been studied at the level of individual cells. While populations of fast-growing cells tend to have more ribosomes than populations of slow-growing cells, it is unclear if individual cells tightly regulate their ribosome content to match their environment. Here, we use recent ground-breaking single-cell RNA sequencing techniques to study this growth law at the single-cell level in two different microbes, S. cerevisiae (a single-celled yeast and eukaryote) and B. subtilis (a bacterium and prokaryote). In both species, we find enormous variation in the ribosomal content of single cells that is not predictive of growth rate. Fast-growing populations include cells showing transcriptional signatures of slow growth and stress, as do cells with the highest ribosome content we survey. Broadening our focus to the levels of non-ribosomal transcripts reveals subpopulations of cells in unique transcriptional states suggestive of divergent growth strategies. These results suggest that single-cell ribosome levels are not finely tuned to match population growth rates or nutrient availability, at least not in a way that can be captured by a unifying law that applies to all cell types. Overall, this work encourages the expansion of these "laws" and other models that predict how growth rates are regulated or how they evolve to consider single-cell heterogeneity.

DarA-the central processing unit for the integration of osmotic with potassium and amino acid homeostasis in Bacillus subtilis

Cyclic di-adenosine monophosphate (c-di-AMP) is a second messenger involved in diverse metabolic processes including osmolyte uptake, cell wall homeostasis, as well as antibiotic and heat resistance. This study investigates the role of the c-di-AMP receptor protein DarA in the osmotic stress response in Bacillus subtilis. Through a series of experiments, we demonstrate that DarA plays a central role in the cellular response to osmotic fluctuations. Our findings show that DarA becomes essential under extreme potassium limitation as well as upon salt stress, highlighting its significance in mediating osmotic stress adaptation. Suppressor screens with darA mutants reveal compensatory mechanisms involving the accumulation of osmoprotectants, particularly potassium and citrulline. Mutations affecting various metabolic pathways, including the citric acid cycle as well as glutamate and arginine biosynthesis, indicate a complex interplay between the osmotic stress response and metabolic regulation. In addition, the growth defects of the darA mutant during potassium starvation and salt stress in a strain lacking the high-affinity potassium uptake systems KimA and KtrAB can be rescued by increased affinity of the remaining potassium channel KtrCD or by increased expression of ktrD, thus resulting in increased potassium uptake. Finally, the darA mutant can respond to salt stress by the increased expression of MleN , which can export sodium ions.IMPORTANCEEnvironmental bacteria are exposed to rapidly changing osmotic conditions making an effective adaptation to these changes crucial for the survival of the cells. In Gram-positive bacteria, the second messenger cyclic di-AMP plays a key role in this adaptation by controlling (i) the influx of physiologically compatible organic osmolytes and (ii) the biosynthesis of such osmolytes. In several bacteria, cyclic di-adenosine monophosphate (c-di-AMP) can bind to a signal transduction protein, called DarA, in Bacillus subtilis. So far, no function for DarA has been discovered in any organism. We have identified osmotically challenging conditions that make DarA essential and have identified suppressor mutations that help the bacteria to adapt to those conditions. Our results indicate that DarA is a central component in the integration of osmotic stress with the synthesis of compatible amino acid osmolytes and with the homeostasis of potassium, the first response to osmotic stress.

Lipid phase separation impairs membrane thickness sensing by the Bacillus subtilis sensor kinase DesK

Membrane fluidity and thickness have emerged as crucial factors for the activity of and resistance to several antimicrobials. However, the lack of tools to study membrane fluidity and, in particular, thickness in living bacteria limits our understanding of this interplay. The Bacillus subtilis histidine kinase/phosphatase DesK is a molecular sensor that directly detects membrane thickness. It controls activity of DesR, which regulates expression of the lipid desaturase Des, known for its role in cold adaptation and daptomycin susceptibility. We hypothesized that this property could be exploited to develop biosensors and reporters for antibiotic-induced changes in membrane fluidity and thickness. To test this, we designed three assays based on the des system: activation of the Pdes promoter as reporter for membrane thickening, localization of DesK-GFP(green-fluorescent protein) as proxy for rigidified membrane domains, and antibiotic sensitivity of des, desK, and desR deletion mutants as readout for the importance of membrane rigidification/thickening under the tested condition. While we could not confirm the suitability of the des system as reporter for antibiotic-induced changes in membrane thickness, we did observe that des expression is only activated by mild temperature shocks, likely due to partitioning of the sensor DesK into fluid membrane domains upon phase separation, precluding effective thickness sensing under harsh cold shock and antibiotic stress conditions. Similarly, we did not observe any sensitivity of the deletion mutants to either temperature or antibiotic stress, raising the question to what extent the des system contributes to fluidity adaptation under these conditions.

IMPORTANCE

The B. subtilis des system is a prime model for direct molecular membrane thickness sensor and, as such, has been well studied in vitro. Our study shows that our understanding of its function in vivo and its importance under temperature and antibiotic stress is still very limited. Specifically, our results suggest that (i) the des system senses very subtle membrane fluidity changes that escape detection by established fluidity reporters like laurdan; (ii) membrane thickness sensing by DesK is impaired by phase separation due to partitioning of the protein into the fluid phase; and (iii) fluidity adaptations by Des are too subtle to elicit growth defects under rigidifying conditions, raising the question of how much the des system contributes to adaptation of overall membrane fluidity.

Enhanced surface colonisation and competition during bacterial adaptation to a fungus

Bacterial-fungal interactions influence microbial community performance of most ecosystems and elicit specific microbial behaviours, including stimulating specialised metabolite production. Here, we use a co-culture experimental evolution approach to investigate bacterial adaptation to the presence of a fungus, using a simple model of bacterial-fungal interactions encompassing the bacterium Bacillus subtilis and the fungus Aspergillus niger. We find in one evolving population that B. subtilis was selected for enhanced production of the lipopeptide surfactin and accelerated surface spreading ability, leading to inhibition of fungal expansion and acidification of the environment. These phenotypes were explained by specific mutations in the DegS-DegU two-component system. In the presence of surfactin, fungal hyphae exhibited bulging cells with delocalised secretory vesicles possibly provoking an RlmA-dependent cell wall stress. Thus, our results indicate that the presence of the fungus selects for increased surfactin production, which inhibits fungal growth and facilitates the competitive success of the bacterium.

Bipartite networks represent causality better than simple networks: evidence, algorithms, and applications

A network, whose nodes are genes and whose directed edges represent positive or negative influences of a regulatory gene and its targets, is often used as a representation of causality. To infer a network, researchers often develop a machine learning model and then evaluate the model based on its match with experimentally verified "gold standard" edges. The desired result of such a model is a network that may extend the gold standard edges. Since networks are a form of visual representation, one can compare their utility with architectural or machine blueprints. Blueprints are clearly useful because they provide precise guidance to builders in construction. If the primary role of gene regulatory networks is to characterize causality, then such networks should be good tools of prediction because prediction is the actionable benefit of knowing causality. But are they? In this paper, we compare prediction quality based on "gold standard" regulatory edges from previous experimental work with non-linear models inferred from time series data across four different species. We show that the same non-linear machine learning models have better predictive performance, with improvements from 5.3% to 25.3% in terms of the reduction in the root mean square error (RMSE) compared with the same models based on the gold standard edges. Having established that networks fail to characterize causality properly, we suggest that causality research should focus on four goals: (i) predictive accuracy; (ii) a parsimonious enumeration of predictive regulatory genes for each target gene g; (iii) the identification of disjoint sets of predictive regulatory genes for each target g of roughly equal accuracy; and (iv) the construction of a bipartite network (whose node types are genes and models) representation of causality. We provide algorithms for all goals.

Bacillus subtilis remains translationally active after CRISPRi-mediated replication initiation arrest

Initiation of bacterial DNA replication takes place at the origin of replication (oriC), a region characterized by the presence of multiple DnaA boxes that serve as the binding sites for the master initiator protein DnaA. This process is tightly controlled by modulation of the availability or activity of DnaA and oriC during development or stress conditions. Here, we aimed to uncover the physiological and molecular consequences of stopping replication in the model bacterium Bacillus subtilis. We successfully arrested replication in B. subtilis by employing a clustered regularly interspaced short palindromic repeats interference (CRISPRi) approach to specifically target the key DnaA boxes 6 and 7, preventing DnaA binding to oriC. In this way, other functions of DnaA, such as a transcriptional regulator, were not significantly affected. When replication initiation was halted by this specific artificial and early blockage, we observed that non-replicating cells continued translation and cell growth, and the initial replication arrest did not induce global stress conditions such as the SOS response.IMPORTANCEAlthough bacteria constantly replicate under laboratory conditions, natural environments expose them to various stresses such as lack of nutrients, high salinity, and pH changes, which can trigger non-replicating states. These states can enable bacteria to (i) become tolerant to antibiotics (persisters), (ii) remain inactive in specific niches for an extended period (dormancy), and (iii) adjust to hostile environments. Non-replicating states have also been studied because of the possibility of repurposing energy for the production of additional metabolites or proteins. Using clustered regularly interspaced short palindromic repeats interference (CRISPRi) targeting bacterial replication initiation sequences, we were able to successfully control replication initiation in Bacillus subtilis. This precise approach makes it possible to study non-replicating phenotypes, contributing to a better understanding of bacterial adaptive strategies.

Control of three-carbon amino acid homeostasis by promiscuous importers and exporters in Bacillus subtilis: role of the "sleeping beauty" amino acid exporters

The Gram-positive model bacterium Bacillus subtilis can acquire amino acids by import, de novo biosynthesis, or degradation of proteins and peptides. The accumulation of several amino acids inhibits the growth of B. subtilis, probably due to misincorporation into cellular macromolecules such as proteins or peptidoglycan or due to interference with other amino acid biosynthetic pathways. Here, we studied the adaptation of B. subtilis to toxic concentrations of the three-carbon amino acids L-alanine, β-alanine, and 2,3-diaminopropionic acid, as well as the two-carbon amino acid glycine. Resistance to the non-proteinogenic amino acid β-alanine, which is a precursor for coenzyme A biosynthesis, is achieved by mutations that either activate a cryptic amino acid exporter, AexA (previously YdeD), or inactivate the amino acid importers AimA, AimB (previously YbxG), and BcaP. The aexA gene is very poorly expressed under most conditions studied. However, mutations affecting the transcription factor AerA (previously YdeC) can result in strong constitutive aexA expression. AexA is the first characterized member of a group of amino acid exporters in B. subtilis, which are all very poorly expressed. Therefore, we suggest to call this group "sleeping beauty amino acid exporters." 2,3-Diaminopropionic acid can also be exported by AexA, and this amino acid also seems to be a natural substrate of AerA/AexA, as it can cause a slight but significant induction of aexA expression, and AexA also provides some natural resistance toward 2,3-diaminopropionic acid. Moreover, our work shows how low-specificity amino acid transporters contribute to amino acid homeostasis in B. subtilis.IMPORTANCEEven though Bacillus subtilis is one of the most-studied bacteria, amino acid homeostasis in this organism is not fully understood. We have identified import and export systems for the C2 and C3 amino acids. Our work demonstrates that the responsible amino acid permeases contribute in a rather promiscuitive way to amino acid uptake. In addition, we have discovered AexA, the first member of a group of very poorly expressed amino acid exporters in B. subtilis that we call "sleeping beauty amino acid exporters." The expression of these transporters is typically triggered by mutations in corresponding regulator genes that are acquired upon exposure to toxic amino acids. These exporters are ubiquitous in all domains of life. It is tempting to speculate that many of them are not expressed until the cells experience selective pressure by toxic compounds, and they protect the cells from rare but potentially dangerous encounters with such compounds.

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.

From the outer space to the inner cell: deconvoluting the complexity of Bacillus subtilis disulfide stress responses by redox state and absolute abundance quantification of extracellular, membrane, and cytosolic proteins

Understanding cellular mechanisms of stress management relies on omics data as a valuable resource. However, the lack of absolute quantitative data on protein abundances remains a significant limitation, particularly when comparing protein abundances across different cell compartments. In this study, we aimed to gain deeper insights into the proteomic responses of the Gram-positive model bacterium Bacillus subtilis to disulfide stress. We determined proteome-wide absolute abundances, focusing on different sub-cellular locations (cytosol and membrane) as well as the extracellular medium, and combined these data with redox state determination. To quantify secreted proteins in the culture medium, we developed a simple and straightforward protocol for the absolute quantification of extracellular proteins in bacteria. We concentrated extracellular proteins, which are highly diluted in the medium, using StrataClean beads along with a set of standard proteins to determine the extent of the concentration step. The resulting data set provides new insights into protein abundances in different sub-cellular compartments and the extracellular medium, along with a comprehensive proteome-wide redox state determination. Our study offers a quantitative understanding of disulfide stress management, protein production, and secretion in B. subtilis.

IMPORTANCE

Stress responses play a crucial role in bacterial survival and adaptation. The ability to quantitatively measure protein abundances and redox states in different cellular compartments and the extracellular environment is essential for understanding stress management mechanisms. In this study, we addressed the knowledge gap regarding absolute quantification of extracellular proteins and compared protein concentrations in various sub-cellular locations and in the extracellular medium under disulfide stress conditions. Our findings provide valuable insights into the protein production and secretion dynamics of B. subtilis, shedding light on its stress response strategies. Furthermore, the developed protocol for absolute quantification of extracellular proteins in bacteria presents a practical and efficient approach for future studies in the field. Overall, this research contributes to the quantitative understanding of stress management mechanisms and protein dynamics in B. subtilis, which can be used to enhance bacterial stress tolerance and protein-based biotechnological applications.

A historical sequence deletion in a commonly used Bacillus subtilis chromosome integration vector generates undetected loss-of-function mutations

Since the 1980s, chromosome-integration vectors have been used as a core method of engineering Bacillus subtilis. One of the most frequently used vector backbones contains chromosomally derived regions that direct homologous recombination into the amyE locus. Here, we report a gap in the homology region inherited from the original amyE integration vector, leading to erroneous recombination in a subset of transformants and a loss-of-function mutation in the downstream gene. Internal to the homology arm that spans the 3' portion of amyE and the downstream gene ldh, an unintentional 227 bp deletion generates two crossover events. The major event yields the intended genotype, but the minor event, occurring in ~10 % of colonies, results in a truncation of ldh, which encodes lactate dehydrogenase. Although both types of colonies test positive for amyE disruption by starch plating, the potential defect in fermentative metabolism may be left undetected and confound the results of subsequent experiments.

Genome-wide annotation of transcript boundaries using bacterial Rend-seq datasets

Accurate annotation to single-nucleotide resolution of the transcribed regions in genomes is key to optimally analyse RNA-seq data, understand regulatory events and for the design of experiments. However, currently most genome annotations provided by GenBank generally lack information about untranslated regions. Additionally, information regarding genomic locations of non-coding RNAs, such as sRNAs, or anti-sense RNAs is frequently missing. To provide such information, diverse RNA-seq technologies, such as Rend-seq, have been developed and applied to many bacterial species. However, incorporating this vast amount of information into annotation files has been limited and is bioinformatically challenging, resulting in UTRs and other non-coding elements being overlooked or misrepresented. To overcome this problem, we present pyRAP (python Rend-seq Annotation Pipeline), a software package that analyses Rend-seq datasets to accurately resolve transcript boundaries genome-wide. We report the use of pyRAP to find novel transcripts, transcript isoforms, and RNase-dependent sRNA processing events. In Bacillus subtilis we uncovered 63 novel transcripts and provide genomic coordinates with single-nucleotide resolution for 2218 5'UTRs, 1864 3'UTRs and 161 non-coding RNAs. In Escherichia coli, we report 117 novel transcripts, 2429 5'UTRs, 1619 3'UTRs and 91 non-coding RNAs, and in Staphylococcus aureus, 16 novel transcripts, 664 5'UTRs, 696 3'UTRs, and 81 non-coding RNAs. Finally, we use pyRAP to produce updated annotation files for B. subtilis 168, E. coli K-12 MG1655, and S. aureus 8325 for use in the wider microbial genomics research community.

Bacillus subtilis as a host for natural product discovery and engineering of biosynthetic gene clusters

Covering: up to October 2023Many bioactive natural products are synthesized by microorganisms that are either difficult or impossible to cultivate under laboratory conditions, or that produce only small amounts of the desired compound. By transferring biosynthetic gene clusters (BGCs) into alternative host organisms that are more easily cultured and engineered, larger quantities can be obtained and new analogues with potentially improved biological activity or other desirable properties can be generated. Moreover, expression of cryptic BGCs in a suitable host can facilitate the identification and characterization of novel natural products. Heterologous expression therefore represents a valuable tool for natural product discovery and engineering as it allows the study and manipulation of their biosynthetic pathways in a controlled setting, enabling innovative applications. Bacillus is a genus of Gram-positive bacteria that is widely used in industrial biotechnology as a host for the production of proteins from diverse origins, including enzymes and vaccines. However, despite numerous successful examples, Bacillus species remain underexploited as heterologous hosts for the expression of natural product BGCs. Here, we review important advantages that Bacillus species offer as expression hosts, such as high secretion capacity, natural competence for DNA uptake, and the increasing availability of a wide range of genetic tools for gene expression and strain engineering. We evaluate different strain optimization strategies and other critical factors that have improved the success and efficiency of heterologous natural product biosynthesis in B. subtilis. Finally, future perspectives for using B. subtilis as a heterologous host are discussed, identifying research gaps and promising areas that require further exploration.

8-OxoG-Dependent Regulation of Global Protein Responses Leads to Mutagenesis and Stress Survival in Bacillus subtilis

The guanine oxidized (GO) system of Bacillus subtilis, composed of the YtkD (MutT), MutM and MutY proteins, counteracts the cytotoxic and genotoxic effects of the oxidized nucleobase 8-OxoG. Here, we report that in growing B. subtilis cells, the genetic inactivation of GO system potentiated mutagenesis (HPM), and subsequent hyperresistance, contributes to the damaging effects of hydrogen peroxide (H2O2) (HPHR). The mechanism(s) that connect the accumulation of the mutagenic lesion 8-OxoG with the ability of B. subtilis to evolve and survive the noxious effects of oxidative stress were dissected. Genetic and biochemical evidence indicated that the synthesis of KatA was exacerbated, in a PerR-independent manner, and the transcriptional coupling repair factor, Mfd, contributed to HPHR and HPM of the ΔGO strain. Moreover, these phenotypes are associated with wider pleiotropic effects, as revealed by a global proteome analysis. The inactivation of the GO system results in the upregulated production of KatA, and it reprograms the synthesis of the proteins involved in distinct types of cellular stress; this has a direct impact on (i) cysteine catabolism, (ii) the synthesis of iron-sulfur clusters, (iii) the reorganization of cell wall architecture, (iv) the activation of AhpC/AhpF-independent organic peroxide resistance, and (v) increased resistance to transcription-acting antibiotics. Therefore, to contend with the cytotoxic and genotoxic effects derived from the accumulation of 8-OxoG, B. subtilis activates the synthesis of proteins belonging to transcriptional regulons that respond to a wide, diverse range of cell stressors.

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.

Control of asparagine homeostasis in Bacillus subtilis: identification of promiscuous amino acid importers and exporters

The Gram-positive model bacterium B. subtilis is able to import all proteinogenic amino acids from the environment as well as to synthesize them. However, the players involved in the acquisition of asparagine have not yet been identified for this bacterium. In this work, we used d-asparagine as a toxic analog of l-asparagine to identify asparagine transporters. This revealed that d- but not l-asparagine is taken up by the malate/lactate antiporter MleN. Specific strains that are sensitive to the presence of l-asparagine due to the lack of the second messenger cyclic di-AMP or due to the intracellular accumulation of this amino acid were used to isolate and characterize suppressor mutants that were resistant to the presence of otherwise growth-inhibiting concentrations of l-asparagine. These screens identified the broad-spectrum amino acid importers AimA and BcaP as responsible for the acquisition of l-asparagine. The amino acid exporter AzlCD allows detoxification of l-asparagine in addition to 4-azaleucine and histidine. This work supports the idea that amino acids are often transported by promiscuous importers and exporters. However, our work also shows that even stereo-enantiomeric amino acids do not necessarily use the same transport systems.IMPORTANCETransport of amino acid is a poorly studied function in many bacteria, including the model organism Bacillus subtilis. The identification of transporters is hampered by the redundancy of transport systems for most amino acids as well as by the poor specificity of the transporters. Here, we apply several strategies to use the growth-inhibitive effect of many amino acids under defined conditions to isolate suppressor mutants that exhibit either reduced uptake or enhanced export of asparagine, resulting in the identification of uptake and export systems for l-asparagine. The approaches used here may be useful for the identification of transporters for other amino acids both in B. subtilis and in other bacteria.

Bacillus subtilis Induces Human Beta Defensin-2 Through its Lipoproteins in Human Intestinal Epithelial Cells

Human intestinal epithelial cells (IECs) play an important role in maintaining gut homeostasis by producing antimicrobial peptides (AMPs). Bacillus subtilis, a commensal bacterium, is considered a probiotic. Although its protective effects on intestinal health are widely reported, the key component of B. subtilis responsible for its beneficial effects remains elusive. In this study, we tried to identify the key molecules responsible for B. subtilis-induced AMPs and their molecular mechanisms in a human IEC line, Caco-2. B. subtilis increased human beta defensin (HBD)-2 mRNA expression in a dose- and time-dependent manner. Among the B. subtilis microbe-associated molecular patterns, lipoprotein (LPP) substantially increased the mRNA expression and protein production of HBD-2, whereas lipoteichoic acid and peptidoglycan did not show such effects. Those results were confirmed in primary human IECs. In addition, both LPP recognition and HBD-2 secretion mainly took place on the apical side of fully differentiated and polarized Caco-2 cells through Toll-like receptor 2-mediated JNK/p38 MAP kinase/AP-1 and NF-κB pathways. HBD-2 efficiently inhibited the growth of the intestinal pathogens Staphylococcus aureus and Bacillus cereus. Furthermore, LPPs pre-incubated with lipase or proteinase K decreased LPP-induced HBD-2 expression, suggesting that the lipid and protein moieties of LPP are crucial for HBD-2 expression. Q Exactive Plus mass spectrometry identified 35 B. subtilis LPP candidates within the LPP preparation, and most of them were ABC transporters. Taken together, these results suggest that B. subtilis promotes HBD-2 secretion in human IECs mainly with its LPPs, which might enhance the protection from intestinal pathogens.

The previously uncharacterized RnpM (YlxR) protein modulates the activity of ribonuclease P in Bacillus subtilis in vitro

Even though Bacillus subtilis is one of the most studied organisms, no function has been identified for about 20% of its proteins. Among these unknown proteins are several RNA- and ribosome-binding proteins suggesting that they exert functions in cellular information processing. In this work, we have investigated the RNA-binding protein YlxR. This protein is widely conserved in bacteria and strongly constitutively expressed in B. subtilis suggesting an important function. We have identified the RNA subunit of the essential RNase P as the binding partner of YlxR. The main activity of RNase P is the processing of 5' ends of pre-tRNAs. In vitro processing assays demonstrated that the presence of YlxR results in reduced RNase P activity. Chemical cross-linking studies followed by in silico docking analysis and experiments with site-directed mutant proteins suggest that YlxR binds to the region of the RNase P RNA that is important for binding and cleavage of the pre-tRNA substrate. We conclude that the YlxR protein is a novel interaction partner of the RNA subunit of RNase P that serves to finetune RNase P activity to ensure appropriate amounts of mature tRNAs for translation. We rename the YlxR protein RnpM for RNase P modulator.

Structural basis of ribosomal 30S subunit degradation by RNase R

Protein synthesis is a major energy-consuming process of the cell that requires the controlled production1-3 and turnover4,5 of ribosomes. Although the past few years have seen major advances in our understanding of ribosome biogenesis, structural insight into the degradation of ribosomes has been lacking. Here we present native structures of two distinct small ribosomal 30S subunit degradation intermediates associated with the 3' to 5' exonuclease ribonuclease R (RNase R). The structures reveal that RNase R binds at first to the 30S platform to facilitate the degradation of the functionally important anti-Shine-Dalgarno sequence and the decoding-site helix 44. RNase R then encounters a roadblock when it reaches the neck region of the 30S subunit, and this is overcome by a major structural rearrangement of the 30S head, involving the loss of ribosomal proteins. RNase R parallels this movement and relocates to the decoding site by using its N-terminal helix-turn-helix domain as an anchor. In vitro degradation assays suggest that head rearrangement poses a major kinetic barrier for RNase R, but also indicate that the enzyme alone is sufficient for complete degradation of 30S subunits. Collectively, our results provide a mechanistic basis for the degradation of 30S mediated by RNase R, and reveal that RNase R targets orphaned 30S subunits using a dynamic mechanism involving an anchored switching of binding sites.

A new family of bacterial ribosome hibernation factors

To conserve energy during starvation and stress, many organisms use hibernation factor proteins to inhibit protein synthesis and protect their ribosomes from damage1,2. In bacteria, two families of hibernation factors have been described, but the low conservation of these proteins and the huge diversity of species, habitats and environmental stressors have confounded their discovery3-6. Here, by combining cryogenic electron microscopy, genetics and biochemistry, we identify Balon, a new hibernation factor in the cold-adapted bacterium Psychrobacter urativorans. We show that Balon is a distant homologue of the archaeo-eukaryotic translation factor aeRF1 and is found in 20% of representative bacteria. During cold shock or stationary phase, Balon occupies the ribosomal A site in both vacant and actively translating ribosomes in complex with EF-Tu, highlighting an unexpected role for EF-Tu in the cellular stress response. Unlike typical A-site substrates, Balon binds to ribosomes in an mRNA-independent manner, initiating a new mode of ribosome hibernation that can commence while ribosomes are still engaged in protein synthesis. Our work suggests that Balon-EF-Tu-regulated ribosome hibernation is a ubiquitous bacterial stress-response mechanism, and we demonstrate that putative Balon homologues in Mycobacteria bind to ribosomes in a similar fashion. This finding calls for a revision of the current model of ribosome hibernation inferred from common model organisms and holds numerous implications for how we understand and study ribosome hibernation.

Beyond blast: enabling microbiologists to better extract literature, taxonomic distributions and gene neighbourhood information for protein families

Capturing the published corpus of information on all members of a given protein family should be an essential step in any study focusing on specific members of that family. Using a previously gathered dataset of more than 280 references mentioning a member of the DUF34 (NIF3/Ngg1-interacting Factor 3) family, we evaluated the efficiency of different databases and search tools, and devised a workflow that experimentalists can use to capture the most information published on members of a protein family in the least amount of time. To complement this workflow, web-based platforms allowing for the exploration of protein family members across sequenced genomes or for the analysis of gene neighbourhood information were reviewed for their versatility and ease of use. Recommendations that can be used for experimentalist users, as well as educators, are provided and integrated within a customized, publicly accessible Wiki.

Targeting the Bottlenecks in Levan Biosynthesis Pathway in Bacillus subtilis and Strain Optimization by Computational Modeling and Omics Integration

Levan is a fructan polymer with many industrial applications such as the formulation of hydrogels, drug delivery, and wound healing, among others. To this end, metabolic systems engineering is a valuable method to improve the yield of a specific metabolite in a wide range of bacterial and eukaryotic organisms. In this study, we report a systems biology approach integrating genomics data for the Bacillus subtilis model, wherein the metabolic pathway for levan biosynthesis is unpacked. We analyzed a revised genome-scale enzyme-constrained metabolic model (ecGEM) and performed simulations to increase levan biopolymer production capacity in B. subtilis. We used the model ec_iYO844_lvn to (1) identify the essential genes and bottlenecks in levan production, and (2) specifically design an engineered B. subtilis strain capable of producing higher levan yields. The FBA and FVA analysis showed the maximal growth rate of the organism up to 0.624 hr-1 at 20 mmol gDw-1 hr-1 of sucrose intake. Gene knockout analyses were performed to identify gene knockout targets to increase the levan flux in B. subtilis. Importantly, we found that the pgk and ctaD genes are the two target genes for the knockout. The perturbation of these two genes has flux gains for levan production reactions with 1.3- and 1.4-fold the relative flux span in the mutant strains, respectively, compared to the wild type. In all, this work identifies the bottlenecks in the production of levan and possible ways to overcome them. Our results provide deeper insights on the bacterium's physiology and new avenues for strain engineering.

Competition for iron shapes metabolic antagonism between Bacillus subtilis and Pseudomonas marginalis

Siderophores have long been implicated in sociomicrobiology as determinants of bacterial interrelations. For plant-associated genera, like Bacillus and Pseudomonas, siderophores are well known for their biocontrol functions. Here, we explored the functional role of the Bacillus subtilis siderophore bacillibactin (BB) in an antagonistic interaction with Pseudomonas marginalis. The presence of BB strongly influenced the outcome of the interaction in an iron-dependent manner. The BB producer B. subtilis restricts colony spreading of P. marginalis by repressing the transcription of histidine kinase-encoding gene gacS, thereby abolishing production of secondary metabolites such as pyoverdine and viscosin. By contrast, lack of BB restricted B. subtilis colony growth. To explore the specificity of the antagonism, we cocultured B. subtilis with a collection of fluorescent Pseudomonas spp. and found that the Bacillus-Pseudomonas interaction is conserved, expanding our understanding of the interplay between two of the most well-studied genera of soil bacteria.

A historical sequence deletion in a commonly used Bacillus subtilis chromosome integration vector generates undetected loss-of-function mutations

Since the 1980s, chromosome-integration vectors have been used as a core method of engineering Bacillus subtilis. One of the most frequently used vector backbones contains chromosomally derived regions that direct homologous recombination into the amyE locus. Here, we report a gap in the homology region inherited from the original amyE integration vector, leading to erroneous recombination in a subset of transformants and a loss-of-function mutation in the downstream gene. Internal to the homology arm that spans the 3' portion of amyE and the downstream gene ldh, an unintentional 227-bp deletion generates two crossover events. The major event yields the intended genotype, but the minor event, occurring in ~10% of colonies, results in a truncation of ldh, which encodes lactate dehydrogenase. Although both types of colonies test positive for amyE disruption by starch plating, the potential defect in fermentative metabolism may be left undetected and confound the results of subsequent experiments.

Let There Be Light: Genome Reduction Enables Bacillus subtilis to Produce Disulfide-Bonded Gaussia Luciferase

Bacillus subtilis is a major workhorse for enzyme production in industrially relevant quantities. Compared to mammalian-based expression systems, B. subtilis presents intrinsic advantages, such as high growth rates, high space-time yield, unique protein secretion capabilities, and low maintenance costs. However, B. subtilis shows clear limitations in the production of biopharmaceuticals, especially proteins from eukaryotic origin that contain multiple disulfide bonds. In the present study, we deployed genome minimization, signal peptide screening, and coexpression of recombinant thiol oxidases as strategies to improve the ability of B. subtilis to secrete proteins with multiple disulfide bonds. Different genome-reduced strains served as the chassis for expressing the model protein Gaussia Luciferase (GLuc), which contains five disulfide bonds. These chassis lack extracellular proteases, prophages, and key sporulation genes. Importantly, compared to the reference strain with a full-size genome, the best-performing genome-minimized strain achieved over 3000-fold increased secretion of active GLuc while growing to lower cell densities. Our results show that high-level GLuc secretion relates, at least in part, to the absence of major extracellular proteases. In addition, we show that the thiol-disulfide oxidoreductase requirements for disulfide bonding have changed upon genome reduction. Altogether, our results highlight genome-engineered Bacillus strains as promising expression platforms for proteins with multiple disulfide bonds.

Bacillus subtilis NDmed, a model strain for biofilm genetic studies

The Bacillus subtilis strain NDmed was isolated from an endoscope washer-disinfector in a medical environment. NDmed can form complex macrocolonies with highly wrinkled architectural structures on solid medium. In static liquid culture, it produces thick pellicles at the interface with air as well as remarkable highly protruding ''beanstalk-like'' submerged biofilm structures at the solid surface. Since these mucoid submerged structures are hyper-resistant to biocides, NDmed has the ability to protect pathogens embedded in mixed-species biofilms by sheltering them from the action of these agents. Additionally, this non-domesticated and highly biofilm forming strain has the propensity of being genetically manipulated. Due to all these properties, the NDmed strain becomes a valuable model for the study of B. subtilis biofilms. This review focuses on several studies performed with NDmed that have highlighted the sophisticated genetic dynamics at play during B. subtilis biofilm formation. Further studies in project using modern molecular tools of advanced technologies with this strain, will allow to deepen our knowledge on the emerging properties of multicellular bacterial life.

Simultaneous spatiotemporal transcriptomics and microscopy of Bacillus subtilis swarm development reveal cooperation across generations

Development of microbial communities is a complex multiscale phenomenon with wide-ranging biomedical and ecological implications. How biological and physical processes determine emergent spatial structures in microbial communities remains poorly understood due to a lack of simultaneous measurements of gene expression and cellular behaviour in space and time. Here we combined live-cell microscopy with a robotic arm for spatiotemporal sampling, which enabled us to simultaneously acquire phenotypic imaging data and spatiotemporal transcriptomes during Bacillus subtilis swarm development. Quantitative characterization of the spatiotemporal gene expression patterns revealed correlations with cellular and collective properties, and phenotypic subpopulations. By integrating these data with spatiotemporal metabolome measurements, we discovered a spatiotemporal cross-feeding mechanism fuelling swarm development: during their migration, earlier generations deposit metabolites which are consumed by later generations that swarm across the same location. These results highlight the importance of spatiotemporal effects during the emergence of phenotypic subpopulations and their interactions in bacterial communities.

Migration of surface-associated microbial communities in spaceflight habitats

Astronauts are spending longer periods locked up in ships or stations for scientific and exploration spatial missions. The International Space Station (ISS) has been inhabited continuously for more than 20 years and the duration of space stays by crews could lengthen with the objectives of human presence on the moon and Mars. If the environment of these space habitats is designed for the comfort of astronauts, it is also conducive to other forms of life such as embarked microorganisms. The latter, most often associated with surfaces in the form of biofilm, have been implicated in significant degradation of the functionality of pieces of equipment in space habitats. The most recent research suggests that microgravity could increase the persistence, resistance and virulence of pathogenic microorganisms detected in these communities, endangering the health of astronauts and potentially jeopardizing long-duration manned missions. In this review, we describe the mechanisms and dynamics of installation and propagation of these microbial communities associated with surfaces (spatial migration), as well as long-term processes of adaptation and evolution in these extreme environments (phenotypic and genetic migration), with special reference to human health. We also discuss the means of control envisaged to allow a lasting cohabitation between these vibrant microscopic passengers and the astronauts.

Direct comparison of spatial transcriptional heterogeneity across diverse Bacillus subtilis biofilm communities

Bacillus subtilis can form various types of spatially organised communities on surfaces, such as colonies, pellicles and submerged biofilms. These communities share similarities and differences, and phenotypic heterogeneity has been reported for each type of community. Here, we studied spatial transcriptional heterogeneity across the three types of surface-associated communities. Using RNA-seq analysis of different regions or populations for each community type, we identified genes that are specifically expressed within each selected population. We constructed fluorescent transcriptional fusions for 17 of these genes, and observed their expression in submerged biofilms using time-lapse confocal laser scanning microscopy (CLSM). We found mosaic expression patterns for some genes; in particular, we observed spatially segregated cells displaying opposite regulation of carbon metabolism genes (gapA and gapB), indicative of distinct glycolytic or gluconeogenic regimes coexisting in the same biofilm region. Overall, our study provides a direct comparison of spatial transcriptional heterogeneity, at different scales, for the three main models of B. subtilis surface-associated communities.

ListiWiki: A database for the foodborne pathogen Listeria monocytogenes

Listeria monocytogenes is a Gram positive foodborne pathogen that regularly causes outbreaks of systemic infectious diseases. The bacterium maintains a facultative intracellular lifestyle; it thrives under a variety of environmental conditions and is able to infect human host cells. L. monocytogenes is genetically tractable and therefore has become an attractive model system to study the mechanisms employed by facultative intracellular bacteria to invade eukaryotic cells and to replicate in their cytoplasm. Besides its importance for basic research, L. monocytogenes also serves as a paradigmatic pathogen in genomic epidemiology, where the relative stability of its genome facilitates successful outbreak detection and elucidation of transmission chains in genomic pathogen surveillance systems. In both terms, it is necessary to keep the annotation of the L. monocytogenes genome up to date. Therefore, we have created the database ListiWiki (http://listiwiki.uni-goettingen.de/) which stores comprehensive information on the widely used L. monocytogenes reference strain EDG-e. ListiWiki is designed to collect information on genes, proteins and RNAs and their relevant functional characteristics, but also further information such as mutant phenotypes, available biological material, and publications. In its present form, ListiWiki combines the most recent annotation of the EDG-e genome with published data on gene essentiality, gene expression and subcellular protein localization. ListiWiki also predicts protein-protein interactions networks based on protein homology to Bacillus subtilis proteins, for which detailed interaction maps have been compiled in the sibling database SubtiWiki. Furthermore, crystallographic information of proteins is made accessible through integration of Protein Structure Database codes and AlphaFold structure predictions. ListiWiki is an easy-to-use web interface that has been developed with a focus on an intuitive access to all information. Use of ListiWiki is free of charge and its content can be edited by all members of the scientific community after registration. In our labs, ListiWiki has already become an important and easy to use tool to quickly access genome annotation details that we can keep updated with advancing knowledge. It also might be useful to promote the comprehensive understanding of the physiology and virulence of an important human pathogen.

Co-evolution at protein-protein interfaces guides inference of stoichiometry of oligomeric protein complexes by de novo structure prediction

The quaternary structure with specific stoichiometry is pivotal to the specific function of protein complexes. However, determining the structure of many protein complexes experimentally remains a major bottleneck. Structural bioinformatics approaches, such as the deep learning algorithm Alphafold2-multimer (AF2-multimer), leverage the co-evolution of amino acids and sequence-structure relationships for accurate de novo structure and contact prediction. Pseudo-likelihood maximization direct coupling analysis (plmDCA) has been used to detect co-evolving residue pairs by statistical modeling. Here, we provide evidence that combining both methods can be used for de novo prediction of the quaternary structure and stoichiometry of a protein complex. We achieve this by augmenting the existing AF2-multimer confidence metrics with an interpretable score to identify the complex with an optimal fraction of native contacts of co-evolving residue pairs at intermolecular interfaces. We use this strategy to predict the quaternary structure and non-trivial stoichiometries of Bacillus subtilis spore germination protein complexes with unknown structures. Co-evolution at intermolecular interfaces may therefore synergize with AI-based de novo quaternary structure prediction of structurally uncharacterized bacterial protein complexes.

An exhaustive multiple knockout approach to understanding cell wall hydrolase function in Bacillus subtilis

IMPORTANCE

In order to grow, bacterial cells must both create and break down their cell wall. The enzymes that are responsible for these processes are the target of some of our best antibiotics. Our understanding of the proteins that break down the wall- cell wall hydrolases-has been limited by redundancy among the large number of hydrolases many bacteria contain. To solve this problem, we identified 42 cell wall hydrolases in Bacillus subtilis and created a strain lacking 40 of them. We show that cells can survive using only a single cell wall hydrolase; this means that to understand the growth of B. subtilis in standard laboratory conditions, it is only necessary to study a very limited number of proteins, simplifying the problem substantially. We additionally show that the ∆40 strain is a research tool to characterize hydrolases, using it to identify three "helper" hydrolases that act in certain stress conditions.

The many roles of cyclic di-AMP to control the physiology of Bacillus subtilis

The dinucleotide cyclic di-AMP (c-di-AMP) is synthesized as a second messenger in the Gram-positive model bacterium Bacillus subtilis as well as in many bacteria and archaea. Bacillus subtilis possesses three diadenylate cyclases and two phosphodiesterases that synthesize and degrade the molecule, respectively. Among the second messengers, c-di-AMP is unique since it is essential for B. subtilis on the one hand but toxic upon accumulation on the other. This role as an "essential poison" is related to the function of c-di-AMP in the control of potassium homeostasis. C-di-AMP inhibits the expression and activity of potassium uptake systems by binding to riboswitches and transporters and activates the activity of potassium exporters. In this way, c-di-AMP allows the adjustment of uptake and export systems to achieve a balanced intracellular potassium concentration. C-di-AMP also binds to two dedicated signal transduction proteins, DarA and DarB. Both proteins seem to interact with other proteins in their apo state, i.e. in the absence of c-di-AMP. For DarB, the (p)ppGpp synthetase/hydrolase Rel and the pyruvate carboxylase PycA have been identified as targets. The interactions trigger the synthesis of the alarmone (p)ppGpp and of the acceptor molecule for the citric acid cycle, oxaloacetate, respectively. In the absence of c-di-AMP, many amino acids inhibit the growth of B. subtilis. This feature can be used to identify novel players in amino acid homeostasis. In this review, we discuss the different functions of c-di-AMP and their physiological relevance.

Fibration symmetry uncovers minimal regulatory networks for logical computation in bacteria

Symmetry principles have proven important in physics, deep learning and geometry, allowing for the reduction of complicated systems to simpler, more comprehensible models that preserve the system's features of interest. Biological systems often show a high level of complexity and consist of a high number of interacting parts. Using symmetry fibrations, the relevant symmetries for biological 'message-passing' networks, we reduced the gene regulatory networks of E. coli and B. subtilis bacteria in a way that preserves information flow and highlights the computational capabilities of the network. Nodes that share isomorphic input trees are grouped into equivalence classes called fibers, whereby genes that receive signals with the same 'history' belong to one fiber and synchronize. We further reduce the networks to its computational core by removing "dangling ends" via k-core decomposition. The computational core of the network consists of a few strongly connected components in which signals can cycle while signals are transmitted between these "information vortices" in a linear feed-forward manner. These components are in charge of decision making in the bacterial cell by employing a series of genetic toggle-switch circuits that store memory, and oscillator circuits. These circuits act as the central computation machine of the network, whose output signals then spread to the rest of the network.

A review on surfactin: molecular regulation of biosynthesis

Surfactin has many biological activities, such as inhibiting plant diseases, resisting bacteria, fungi, viruses, tumors, mycoplasma, anti-adhesion, etc. It has great application potential in agricultural biological control, clinical medical treatment, environmental treatment and other fields. However, the low yield has been the bottleneck of its popularization and application. It is very important to understand the synthesis route and control strategy of surfactin to improve its yield and purity. In this paper, based on the biosynthetic pathway and regulatory factors of surfactin, its biosynthesis regulation strategy was comprehensively summarized, involving enhancement of endogenous and exogenous precursor supply, modification of the synthesis pathway of lipid chain and peptide chain, improvement of secretion and efflux, and manipulation some global regulatory factors, such as Spo0A, AbrB, ComQXP, phrCSF, etc. to directly or indirectly stimulate surfactin synthesis. And the current production and separation and purification process of surfactin are briefly described. This review also provides a scientific reference for promoting surfactin production and its applications in various fields.

High-throughput functional genomics: A (myco)bacterial perspective

Advances in sequencing technologies have enabled unprecedented insights into bacterial genome composition and dynamics. However, the disconnect between the rapid acquisition of genomic data and the (much slower) confirmation of inferred genetic function threatens to widen unless techniques for fast, high-throughput functional validation can be applied at scale. This applies equally to Mycobacterium tuberculosis, the leading infectious cause of death globally and a pathogen whose genome, despite being among the first to be sequenced two decades ago, still contains many genes of unknown function. Here, we summarize the evolution of bacterial high-throughput functional genomics, focusing primarily on transposon (Tn)-based mutagenesis and the construction of arrayed mutant libraries in diverse bacterial systems. We also consider the contributions of CRISPR interference as a transformative technique for probing bacterial gene function at scale. Throughout, we situate our analysis within the context of functional genomics of mycobacteria, focusing specifically on the potential to yield insights into M. tuberculosis pathogenicity and vulnerabilities for new drug and regimen development. Finally, we offer suggestions for future approaches that might be usefully applied in elucidating the complex cellular biology of this major human pathogen.

Glucose controls manganese homeostasis through transcription factors regulating known and newly identified manganese transporter genes in Bacillus subtilis

Mn2+ is an essential nutrient whose concentration is tightly controlled in bacteria. In Bacillus subtilis, the Mn2+-activated transcription factor MntR controls Mn2+ transporter genes. However, factors regulating intracellular Mn2+ concentration are incompletely understood. Here, we found that glucose addition induces an increase in intracellular Mn2+ concentration. We determined this upshift was mediated by glucose induction of the major Mn2+ importer gene mntH by the transcription factor AhrC, which is known to be involved in arginine metabolism and to be indirectly induced by glucose. In addition, we identified novel AhrC-regulated genes encoding the Mn2+ importer YcsG and the ABC-type exporter YknUV. We found the expression of these genes was also regulated by glucose and contributes to the glucose induction of Mn2+ concentrations. ycsG expression is regulated by MntR as well. Furthermore, we analyzed the interaction of AhrC and MntR with the promoter driving ycsG expression and examined the Mn2+-dependent induction of this promoter to identify the transcription factors responsible for the Mn2+ induction. RNA-Seq revealed that disruption of ahrC and mntR affected the expression of 502 and 478 genes, respectively (false discovery rate, <0.001, log2[fold change] ≥ |2|. The AhrC- and/or MntR-dependent expression of twenty promoters was confirmed by LacZ analysis, and AhrC or MntR binding to some of these promoters was observed via EMSA. The finding that glucose promotes an increase in intracellular Mn2+ levels without changes in extracellular Mn2+ concentrations is reasonable for the bacterium, as intracellular Mn2+ is required for enzymes and pathways mediating glucose metabolism.

Plasmid - Chromosome interplay in natural and non-natural hosts: global transcription study of three Bacillus cereus group strains carrying pCER270 plasmid

The Bacillus cereus group comprises genetically related Gram-positive spore-forming bacteria that colonize a wide range of ecological niches and hosts. Despite their high degree of genome conservation, extrachromosomal genetic material diverges between these species. The discriminating properties of the B. cereus group strains are mainly due to plasmid-borne toxins, reflecting the importance of horizontal gene transfers in bacterial evolution and species definition. To investigate how a newly acquired megaplasmid can impact the transcriptome of its host, we transferred the pCER270 from the emetic B. cereus strains to phylogenetically distant B. cereus group strains. RNA-sequencing experiments allowed us to determine the transcriptional influence of the plasmid on host gene expression and the impact of the host genomic background on the pCER270 gene expression. Our results show a transcriptional cross-regulation between the megaplasmid and the host genome. pCER270 impacted carbohydrate metabolism and sporulation genes expression, with a higher effect in the natural host of the plasmid, suggesting a role of the plasmid in the adaptation of the carrying strain to its environment. In addition, the host genomes also modulated the expression of pCER270 genes. Altogether, these results provide an example of the involvement of megaplasmids in the emergence of new pathogenic strains.

DNA damage repair proteins across the Tree of Life

Genome maintenance is orchestrated by a highly regulated DNA damage response with specific DNA repair pathways. Here, we investigate the phylogenetic diversity in the recognition and repair of three well-established DNA lesions, primarily repaired by base excision repair (BER) and ribonucleotide excision repair (RER): (1) 8-oxoguanine, (2) abasic site, and (3) incorporated ribonucleotide in DNA in 11 species: Escherichia coli, Bacillus subtilis, Halobacterium salinarum, Trypanosoma brucei, Tetrahymena thermophila, Saccharomyces cerevisiae, Schizosaccharomyces pombe, Caenorhabditis elegans, Homo sapiens, Arabidopsis thaliana, and Zea mays. Using quantitative mass spectrometry, we identified 337 binding proteins across these species. Of these proteins, 99 were previously characterized to be involved in DNA repair. Through orthology, network, and domain analysis, we linked 44 previously unconnected proteins to DNA repair. Our study presents a resource for future study of the crosstalk and evolutionary conservation of DNA damage repair across all domains of life.

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.

PAE viewer: a webserver for the interactive visualization of the predicted aligned error for multimer structure predictions and crosslinks

The development of AlphaFold for protein structure prediction has opened a new era in structural biology. This is even more the case for AlphaFold-Multimer for the prediction of protein complexes. The interpretation of these predictions has become more important than ever, but it is difficult for the non-specialist. While an evaluation of the prediction quality is provided for monomeric protein predictions by the AlphaFold Protein Structure Database, such a tool is missing for predicted complex structures. Here, we present the PAE Viewer webserver (http://www.subtiwiki.uni-goettingen.de/v4/paeViewerDemo), an online tool for the integrated visualization of predicted protein complexes using a 3D structure display combined with an interactive representation of the Predicted Aligned Error (PAE). This metric allows an estimation of the quality of the prediction. Importantly, our webserver also allows the integration of experimental cross-linking data which helps to interpret the reliability of the structure predictions. With the PAE Viewer, the user obtains a unique online tool which for the first time allows the intuitive evaluation of the PAE for protein complex structure predictions with integrated crosslinks.

Valorization of soybean pulp for sustainable α-ketoisocaproate production using engineered Bacillus subtilis whole-cell biocatalyst

The disposal of soybean pulp (okara) (∼14 M tons annually) represents a global concern. α-ketoisocaproate (KIC) is an intrinsic l-leucine metabolite boosting mammalian muscle growth and has great potential in animal husbandry. However, the use of pure l-leucine (5000 USD/kg) for KIC (22 USD/kg) bioproduction is cost-prohibitive in practice, while okara rich in l-leucine (10%) could serve as an economical alternative. Following the concept of a circular bioeconomy, we managed to develop a cost-efficient platform to valorize okara into KIC. In this study, proteolytic Bacillus subtilis strain 168 capable of utilizing okara as a comprehensive substrate was employed as the whole-cell biocatalyst for KIC bioproduction. First, we elucidated the function of genes involved in KIC downstream metabolism in strain 168, including those encoding 2-oxoisovalerate dehydrogenase (bkdAA), 2-oxoisovalerate decarboxylase (bkdAB), enoyl-CoA hydratase (fadB), and bifunctional enoyl-CoA hydratase/3-hydroxyacyl-CoA dehydrogenase (fadN). Among those KIC downstream metabolizing mutants of strain 168, the 2-oxoisovalerate decarboxylase gene knockout strain (ΔbkdAB) was found to have a better accumulation of KIC. To further improve the KIC yield, a soluble l-amino acid deaminase (LAAD) from Proteus vulgaris was heterologously expressed in the ΔbkdAB strain and a ∼50% conversion of total l-leucine contained in okara was catalyzed into KIC, along with a ∼50% reduction of CO₂ emission compared to the wild-type cultures. Altogether, this renovated biocatalytic system provides an alternative platform to valorize okara for producing value-added chemicals in an eco-friendly manner.

A Sporulation-Independent Way of Life for Bacillus thuringiensis in the Late Stages of an Infection

The formation of endospores has been considered the unique survival and transmission mode of sporulating Firmicutes due to the exceptional resistance and persistence of this bacterial form. However, nonsporulated bacteria (Spo^-) were reported at the early stages following the death of a host infected with Bacillus thuringiensis, an entomopathogenic sporulating bacterium. Here, we investigated the characteristics of the bacterial population in the late stages of an infection in the B. thuringiensis/Galleria mellonella infection model. Using fluorescent reporters and molecular markers coupled to flow cytometry, we demonstrated that the Spo^- cells persist and constitute about half of the population 2 weeks post-infection (p.i.). Protein synthesis and growth recovery assays indicated that they are in a metabolically slowed-down state. These bacteria were extremely resistant to the insect cadaver environment, which did not support growth of in vitro-grown vegetative cells and spores. A transcriptomic analysis of this subpopulation at 7 days p.i. revealed a signature profile of this state, and the expression analysis of individual genes at the cell level showed that more bacteria mount an oxidative stress response as their survival time increases, in agreement with the increase of the free radical level in the host cadaver and in the number of reactive oxygen species (ROS)-producing bacteria. Altogether, these data show for the first time that nonsporulated bacteria are able to survive for a prolonged period of time in the context of an infection and indicate that they engage in a profound adaptation process that leads to their persistence in the host cadaver. IMPORTANCE Bacillus thuringiensis is an entomopathogenic bacterium widely used as a biopesticide. It belongs to the Bacillus cereus group, comprising the foodborne pathogen B. cereus sensu stricto and the anthrax agent Bacillus anthracis. Like other Firmicutes when they encounter harsh conditions, these Gram-positive bacteria can form dormant cells called spores. Due to its highly resistant nature, the spore was considered the unique mode of long-term survival, eclipsing any other form of persistence. Breaking this paradigm, we observed that B. thuringiensis was able to persist in its host cadaver in a nonsporulated form for at least 14 days. Our results show that these bacteria survived in the cadaver environment, which proved hostile for actively growing bacteria by engaging in a profound adaptation process. Studying this facet of the life cycle of a sporulating bacterium provides new fundamental knowledge and might lead to the development of strategies to combat sporulating pathogenic species.

Structural Analysis of Bacillus subtilis Sigma Factors

Bacteria use an array of sigma factors to regulate gene expression during different stages of their life cycles. Full-length, atomic-level structures of sigma factors have been challenging to obtain experimentally as a result of their many regions of intrinsic disorder. AlphaFold has now supplied plausible full-length models for most sigma factors. Here we discuss the current understanding of the structures and functions of sigma factors in the model organism, Bacillus subtilis, and present an X-ray crystal structure of a region of B. subtilis SigE, a sigma factor that plays a critical role in the developmental process of spore formation.

Protein complexes in cells by AI-assisted structural proteomics

Accurately modeling the structures of proteins and their complexes using artificial intelligence is revolutionizing molecular biology. Experimental data enable a candidate-based approach to systematically model novel protein assemblies. Here, we use a combination of in-cell crosslinking mass spectrometry and co-fractionation mass spectrometry (CoFrac-MS) to identify protein-protein interactions in the model Gram-positive bacterium Bacillus subtilis. We show that crosslinking interactions prior to cell lysis reveals protein interactions that are often lost upon cell lysis. We predict the structures of these protein interactions and others in the SubtiWiki database with AlphaFold-Multimer and, after controlling for the false-positive rate of the predictions, we propose novel structural models of 153 dimeric and 14 trimeric protein assemblies. Crosslinking MS data independently validates the AlphaFold predictions and scoring. We report and validate novel interactors of central cellular machineries that include the ribosome, RNA polymerase, and pyruvate dehydrogenase, assigning function to several uncharacterized proteins. Our approach uncovers protein-protein interactions inside intact cells, provides structural insight into their interaction interfaces, and is applicable to genetically intractable organisms, including pathogenic bacteria.

A model industrial workhorse: Bacillus subtilis strain 168 and its genome after a quarter of a century

The vast majority of genomic sequences are automatically annotated using various software programs. The accuracy of these annotations depends heavily on the very few manual annotation efforts that combine verified experimental data with genomic sequences from model organisms. Here, we summarize the updated functional annotation of Bacillus subtilis strain 168, a quarter century after its genome sequence was first made public. Since the last such effort 5 years ago, 1168 genetic functions have been updated, allowing the construction of a new metabolic model of this organism of environmental and industrial interest. The emphasis in this review is on new metabolic insights, the role of metals in metabolism and macromolecule biosynthesis, functions involved in biofilm formation, features controlling cell growth, and finally, protein agents that allow class discrimination, thus allowing maintenance management, and accuracy of all cell processes. New 'genomic objects' and an extensive updated literature review have been included for the sequence, now available at the International Nucleotide Sequence Database Collaboration (INSDC: AccNum AL009126.4).

Understudied proteins and understudied functions in the model bacterium Bacillus subtilis-A major challenge in current research

Model organisms such as the Gram-positive bacterium Bacillus subtilis have been studied intensively for decades. However, even for model organisms, no function has been identified for about one fourth of all proteins. It has recently been realized that such understudied proteins as well as poorly studied functions set a limitation to our understanding of the requirements for cellular life, and the Understudied Proteins Initiative has been launched. Of poorly studied proteins, those that are strongly expressed are likely to be important to the cell and should therefore be considered high priority in further studies. Since the functional analysis of unknown proteins can be extremely laborious, a minimal knowledge is required prior to targeted functional studies. In this review, we discuss strategies to obtain such a minimal annotation, for example, from global interaction, expression, or localization studies. We present a set of 41 highly expressed and poorly studied proteins of B. subtilis. Several of these proteins are thought or known to bind RNA and/or the ribosome, some may control the metabolism of B. subtilis, and another subset of particularly small proteins may act as regulatory elements to control the expression of downstream genes. Moreover, we discuss the challenges of poorly studied functions with a focus on RNA-binding proteins, amino acid transport, and the control of metabolic homeostasis. The identification of the functions of the selected proteins not only will strongly advance our knowledge on B. subtilis, but also on other organisms since many of the proteins are conserved in many groups of bacteria.

Insights into the Roles of Lipoteichoic Acids and MprF in Bacillus subtilis

Gram-positive bacterial cells are protected from the environment by a cell envelope that is comprised of a thick layer of peptidoglycan that maintains cell shape and teichoic acid polymers whose biological function remains unclear. In Bacillus subtilis, the loss of all class A penicillin-binding proteins (aPBPs), which function in peptidoglycan synthesis, is conditionally lethal. Here, we show that this lethality is associated with an alteration of lipoteichoic acids (LTAs) and the accumulation of the major autolysin LytE in the cell wall. Our analysis provides further evidence that the length and abundance of LTAs act to regulate the cellular level and activity of autolytic enzymes, specifically LytE. Importantly, we identify a novel function for the aminoacyl-phosphatidylglycerol synthase MprF in the modulation of LTA biosynthesis in both B. subtilis and Staphylococcus aureus. This finding has implications for our understanding of antimicrobial resistance (particularly to daptomycin) in clinically relevant bacteria and the involvement of MprF in the virulence of pathogens such as methicillin-resistant S. aureus (MRSA). IMPORTANCE In Gram-positive bacteria such as Bacillus subtilis and Staphylococcus aureus, the cell envelope is a structure that protects the cells from the environment but is also dynamic in that it must be modified in a controlled way to allow cell growth. In this study, we show that lipoteichoic acids (LTAs), which are anionic polymers attached to the membrane, have a direct role in modulating the cellular abundance of cell wall-degrading enzymes. We also find that the apparent length of the LTA is modulated by the activity of the enzyme MprF, previously implicated in modifications of the cell membrane leading to resistance to antimicrobial peptides. These findings are important contributions to our understanding of how bacteria balance cell wall synthesis and degradation to permit controlled growth and division. These results also have implications for the interpretation of antibiotic resistance, particularly for the clinical treatment of MRSA infections.