Post-transcriptionally generated cell heterogeneity regulates biofilm formation in Bacillus subtilis

Regulation of biofilm gene expression by DNA replication in Bacillus subtilis

Bacillus subtilis relies on biofilms for survival in harsh environments. Extracellular polymeric substance (EPS) is a crucial component of biofilms, yet the dynamics of EPS production in single cells remain elusive. To unveil the modulation of EPS synthesis, we built a minimal network model comprising the SinI-SinR-SlrR module, Spo0A, and EPS. Stochastic simulations revealed that antagonistic interplay between SinI and SinR enables EPS production in bursts. SlrR widens these bursts and increases their frequency by stabilizing SinR-SlrR complexes and depleting free SinR. DNA replication and chromosomal positioning of key genes dictate pulsatile changes in the slrR:sinR gene dosage ratio (gr) and Spo0A-P levels, each promoting EPS production in distinct phases of the cell cycle. As the cell cycle lengthens with nutrient stress, the duty cycle of gr pulsing decreases, whereas the amplitude of Spo0A-P pulses elevates. This coordinated response facilitates keeping a constant proportion of EPS-secreting cells within colonies across diverse nutrient conditions. Our results suggest that bacteria may 'encode' eps expression through strategic chromosomal organization. This work illuminates how stochastic protein interactions, gene copy number imbalance, and cell-cycle dynamics orchestrate EPS synthesis, offering a deeper understanding of biofilm formation.

The Slowdown of Growth Rate Controls the Single-Cell Distribution of Biofilm Matrix Production via an SinI-SinR-SlrR Network

In Bacillus subtilis, master regulator Spo0A controls several cell-differentiation pathways. Under moderate starvation, phosphorylated Spo0A (Spo0A~P) induces biofilm formation by indirectly activating genes controlling matrix production in a subpopulation of cells via an SinI-SinR-SlrR network. Under severe starvation, Spo0A~P induces sporulation by directly and indirectly regulating sporulation gene expression. However, what determines the heterogeneity of individual cell fates is not fully understood. In particular, it is still unclear why, despite being controlled by a single master regulator, biofilm matrix production and sporulation seem mutually exclusive on a single-cell level. In this work, with mathematical modeling, we showed that the fluctuations in the growth rate and the intrinsic noise amplified by the bistability in the SinI-SinR-SlrR network could explain the single-cell distribution of matrix production. Moreover, we predicted an incoherent feed-forward loop; the decrease in the cellular growth rate first activates matrix production by increasing in Spo0A phosphorylation level but then represses it via changing the relative concentrations of SinR and SlrR. Experimental data provide evidence to support model predictions. In particular, we demonstrate how the degree to which matrix production and sporulation appear mutually exclusive is affected by genetic perturbations. IMPORTANCE The mechanisms of cell-fate decisions are fundamental to our understanding of multicellular organisms and bacterial communities. However, even for the best-studied model systems we still lack a complete picture of how phenotypic heterogeneity of genetically identical cells is controlled. Here, using B. subtilis as a model system, we employ a combination of mathematical modeling and experiments to explain the population-level dynamics and single-cell level heterogeneity of matrix gene expression. The results demonstrate how the two cell fates, biofilm matrix production and sporulation, can appear mutually exclusive without explicitly inhibiting one another. Such a mechanism could be used in a wide range of other biological systems.

The Biofilm Regulatory Network from Bacillus subtilis: A Structure-Function Analysis

Bacterial biofilms are notorious for their ability to protect bacteria from environmental challenges, most importantly the action of antibiotics. Bacillus subtilis is an extensively studied model organism used to understand the process of biofilm formation. A complex network of principal regulatory proteins including Spo0A, AbrB, AbbA, Abh, SinR, SinI, SlrR, and RemA, work in concert to transition B. subtilis from the free-swimming planktonic state to the biofilm state. In this review, we explore, connect, and summarize decades worth of structural and biochemical studies that have elucidated this protein signaling network. Since structure dictates function, unraveling aspects of protein molecular mechanisms will allow us to devise ways to exploit critical features of the biofilm regulatory pathway, such as possible therapeutic intervention. This review pools our current knowledge base of B. subtilis biofilm regulatory proteins and highlights potential therapeutic intervention points.

Bacillus subtilis Histidine Kinase KinC Activates Biofilm Formation by Controlling Heterogeneity of Single-Cell Responses

In Bacillus subtilis, biofilm and sporulation pathways are both controlled by a master regulator, Spo0A, which is activated by phosphorylation via a phosphorelay-a cascade of phosphotransfer reactions commencing with autophosphorylation of histidine kinases KinA, KinB, KinC, KinD, and KinE. However, it is unclear how the kinases, despite acting via the same regulator, Spo0A, differentially regulate downstream pathways, i.e., how KinA mainly activates sporulation genes and KinC mainly activates biofilm genes. In this work, we found that KinC also downregulates sporulation genes, suggesting that KinC has a negative effect on Spo0A activity. To explain this effect, with a mathematical model of the phosphorelay, we revealed that unlike KinA, which always activates Spo0A, KinC has distinct effects on Spo0A at different growth stages: during fast growth, KinC acts as a phosphate source and activates Spo0A, whereas during slow growth, KinC becomes a phosphate sink and contributes to decreasing Spo0A activity. However, under these conditions, KinC can still increase the population-mean biofilm matrix production activity. In a population, individual cells grow at different rates, and KinC would increase the Spo0A activity in the fast-growing cells but reduce the Spo0A activity in the slow-growing cells. This mechanism reduces single-cell heterogeneity of Spo0A activity, thereby increasing the fraction of cells that activate biofilm matrix production. Thus, KinC activates biofilm formation by controlling the fraction of cells activating biofilm gene expression. IMPORTANCE In many bacterial and eukaryotic systems, multiple cell fate decisions are activated by a single master regulator. Typically, the activities of the regulators are controlled posttranslationally in response to different environmental stimuli. The mechanisms underlying the ability of these regulators to control multiple outcomes are not understood in many systems. By investigating the regulation of Bacillus subtilis master regulator Spo0A, we show that sensor kinases can use a novel mechanism to control cell fate decisions. By acting as a phosphate source or sink, kinases can interact with one another and provide accurate regulation of the phosphorylation level. Moreover, this mechanism affects the cell-to-cell heterogeneity of the transcription factor activity and eventually determines the fraction of different cell types in the population. These results demonstrate the importance of intercellular heterogeneity for understanding the effects of genetic perturbations on cell fate decisions. Such effects can be applicable to a wide range of cellular systems.

Glucose-Mediated Protein Arginine Phosphorylation/Dephosphorylation Regulates ylxR Encoding Nucleoid-Associated Protein and Cell Growth in Bacillus subtilis

Glucose is the most favorable carbon source for many bacteria, and these bacteria have several glucose-responsive networks. We proposed new glucose responsive system, which includes protein acetylation and probable translation control through TsaEBD, which is a tRNA modification enzyme required for the synthesis of threonylcarbamoyl adenosine (t⁶A)-tRNA. The system also includes nucleoid-associated protein YlxR, regulating more than 400 genes including many metabolic genes and the ylxR-containing operon driven by the PylxS promoter is induced by glucose. Thus, transposon mutagenesis was performed for searching regulatory factors for PylxS expression. As a result, ywlE was identified. The McsB kinase phosphorylates arginine (Arg) residues of proteins and the YwlE phosphatase counteracts against McsB through Arg-dephosphorylation. Phosphorylated Arg has been known to function as a tag for ClpCP-dependent protein degradation. The previous analysis identified TsaD as an Arg-phosphorylated protein. Our results showed that the McsB/YwlE system regulates PylxS expression through ClpCP-mediated protein degradation of TsaD. In addition, we observed that glucose induced ywlE expression and repressed mcsB expression. It was concluded that these phenomena would cause glucose induction (GI) of PylxS, based on the Western blot analyses of TsaD-FLAG. These observations and the previous those that many glycolytic enzymes are Arg-phosphorylated suggested that the McsB/YwlE system might be involved in cell growth in glucose-containing medium. We observed that the disruption of mcsB and ywlE resulted in an increase of cell mass and delayed growth, respectively, in semi-synthetic medium. These results provide us broader insights to the physiological roles of the McsB/YwlE system and protein Arg-phosphorylation.

Bacillus subtilis Nucleoid-Associated Protein YlxR Is Involved in Bimodal Expression of the Fructoselysine Utilization Operon (frlBONMD-yurJ) Promoter

Bacteria must survive harsh environmental fluctuations at times and have evolved several strategies. “Collective” behaviors have been identified due to recent progress in single-cell analysis. Since most bacteria exist as single cells, bacterial populations are often considered clonal. However, accumulated evidence suggests this is not the case. Gene expression and protein expression are often not homogeneous, resulting in phenotypic heterogeneity. In extreme cases, this leads to bistability, the existence of two stable states. In many cases, expression of key master regulators is bimodal via positive feedback loops causing bimodal expression of the target genes. We observed bimodal expression of metabolic genes for alternative carbon sources. Expression profiles of the frlBONMD-yurJ operon driven by the frlB promoter (PfrlB), which encodes degradation enzymes and a transporter for amino sugars including fructoselysine, were investigated using transcriptional lacZ and gfp, and translational fluorescence reporter mCherry fusions. Disruption effects of genes encoding CodY, FrlR, RNaseY, and nucleoid-associated protein YlxR, four known regulatory factors for PfrlB, were examined for expression of each fusion construct. Expression of PfrlB-gfp and PfrlB-mCherry, which were located at amyE and its original locus, respectively, was bimodal; and disruption of ylxR resulted in the disappearance of the clear bimodal expression pattern in flow cytometric analyses. This suggested a role for YlxR on the bimodal expression of PfrlB. The data indicated that YlxR acted on the bimodal expression of PfrlB through both transcription and translation. YlxR regulates many genes, including those related to translation, supporting the above notion. Depletion of RNaseY abolished heterogenous expression of transcriptional PfrlB-gfp but not bimodal expression of translational PfrlB-mCherry, suggesting the role of RNaseY in regulation of the operon through mRNA stability control and regulatory mechanism for PfrlB-mCherry at the translational level. Based on these results, we discuss the meaning and possible cause of bimodal PfrlB expression.

Bacillus subtilis YlxR, Which Is Involved in Glucose-Responsive Metabolic Changes, Regulates Expression of tsaD for Protein Quality Control of Pyruvate Dehydrogenase

Glucose is the most favorable carbon source for many bacteria, which have several glucose-responsive gene networks. Recently, we found that in Bacillus subtilis glucose induces the expression of the extracellular sigma factor genes sigX and sigM through the acetylation of CshA (RNA helicase), which associates with RNA polymerase (RNAP). We performed a transposon mutagenesis screen for mutants with no glucose induction (GI) of sigX-lacZ. While screening for such mutants, we recently found that the GI of sigX/M involves YlxR, a nucleoid-associated protein (NAP) that regulates nearly 400 genes, including metabolic genes. It has been shown that acetylated CshA positively regulates expression of ylxR-containing operon. Here, we report additional mutations in yqfO or tsaD required for the GI of sigX. YqfO contains a universally conserved domain with unknown function. YqfO and YlxR were found to regulate expression of the tsaEBD-containing operon. Mutational analysis using lacZ fusions revealed the adenine-rich cis-element for YlxR. TsaD is a component of the TsaEBD enzyme required for the synthesis of threonylcarbamoyl adenosine (t6A). The t6A modification of tRNA is universal across the three domains of life. Western blot analysis showed that the tsaD mutation in the presence of glucose reduced levels of soluble PdhA, PdhB, and PdhD, which are subunits of the pyruvate dehydrogenase complex (PDHc). This resulted in severely defective PDHc function and thus reduced concentrations of cellular acetyl-CoA, a reaction product of PDHc and plausible source for CshA acetylation. Thus, we discuss a suggested glucose-responsive system (GRS) involving self-reinforcing CshA acetylation. This self-reinforcing pathway may contribute to the maintenance of the acetyl-CoA pool for protein acetylation.

Newly Identified Nucleoid-Associated-Like Protein YlxR Regulates Metabolic Gene Expression in Bacillus subtilis

Expression of genes encoding NAPs is often temporally regulated. According to results from single-cell analysis, the ylxR gene is induced by glucose and expressed in a bistable mode. These characteristics have not previously been reported for NAP gene expression. Transcriptional profiling of the ylxR disruptant revealed a change in the expression levels of approximately 400 genes, including genes for synthesis of 12 amino acids and 4 nucleotides, in addition to the SigX/M regulons. Thus, YlxR is a critical regulator of glucose response in B. subtilis.

Glucose is the most favorable carbon source for the majority of bacteria, which have several glucose-responsive gene networks. Recently, we found that in Bacillus subtilis, glucose induces expression of the extracellular sigma factor genes sigX/M. To explore the factors affecting this phenomenon, we performed a transposon mutagenesis screen for mutants with no glucose induction (GI) of sigX-lacZ and identified ylxR. YlxR is widely conserved in eubacteria. Further analysis revealed that ylxR is induced by glucose addition. In vitro DNA-binding and cytological studies suggested that YlxR is a nucleoid-associated protein (NAP) in B. subtilis. In many cases, NAPs influence transcription, recombination, and genome stability. Thus, we performed transcriptome sequencing (RNA-Seq) analysis to evaluate the impact of ylxR disruption on the transcriptome in the presence of glucose and observed that YlxR has a profound impact on metabolic gene expression in addition to that of four sigma factor genes. The wide fluctuations of gene expression may result in abolition of GI of sigX/M in the ylxR disruptant.

IMPORTANCE Expression of genes encoding NAPs is often temporally regulated. According to results from single-cell analysis, the ylxR gene is induced by glucose and expressed in a bistable mode. These characteristics have not previously been reported for NAP gene expression. Transcriptional profiling of the ylxR disruptant revealed a change in the expression levels of approximately 400 genes, including genes for synthesis of 12 amino acids and 4 nucleotides, in addition to the SigX/M regulons. Thus, YlxR is a critical regulator of glucose response in B. subtilis.

Dual Regulation of Bacillus subtilis kinB Gene Encoding a Sporulation Trigger by SinR through Transcription Repression and Positive Stringent Transcription Control

It is known that transcription of kinB encoding a trigger for Bacillus subtilis sporulation is under repression by SinR, a master repressor of biofilm formation, and under positive stringent transcription control depending on the adenine species at the transcription initiation nucleotide (nt). Deletion and base substitution analyses of the kinB promoter (P_kinB) region using lacZ fusions indicated that either a 5-nt deletion (Δ5, nt -61/-57, +1 is the transcription initiation nt) or the substitution of G at nt -45 with A (G-45A) relieved kinB repression. Thus, we found a pair of SinR-binding consensus sequences (GTTCTYT; Y is T or C) in an inverted orientation (SinR-1) between nt -57/-42, which is most likely a SinR-binding site for kinB repression. This relief from SinR repression likely requires SinI, an antagonist of SinR. Surprisingly, we found that SinR is essential for positive stringent transcription control of P_kinB. Electrophoretic mobility shift assay (EMSA) analysis indicated that SinR bound not only to SinR-1 but also to SinR-2 (nt -29/-8) consisting of another pair of SinR consensus sequences in a tandem repeat arrangement; the two sequences partially overlap the ‘-35’ and ‘-10’ regions of P_kinB. Introduction of base substitutions (T-27C C-26T) in the upstream consensus sequence of SinR-2 affected positive stringent transcription control of P_kinB, suggesting that SinR binding to SinR-2 likely causes this positive control. EMSA also implied that RNA polymerase and SinR are possibly bound together to SinR-2 to form a transcription initiation complex for kinB transcription. Thus, it was suggested in this work that derepression of kinB from SinR repression by SinI induced by Spo0A∼P and occurrence of SinR-dependent positive stringent transcription control of kinB might induce effective sporulation cooperatively, implying an intimate interplay by stringent response, sporulation, and biofilm formation.

Glucose Induces ECF Sigma Factor Genes, sigX and sigM, Independent of Cognate Anti-sigma Factors through Acetylation of CshA in Bacillus subtilis

Extracytoplasmic function (ECF) σ factors have roles related to cell envelope and/or cell membrane functions, in addition to other cellular functions. Without cell-surface stresses, ECF σ factors are sequestered by the cognate anti-σ factor, leading to inactivation and the resultant repression of regulons due to the inhibition of transcription of their own genes. Bacillus subtilis has seven ECF σ factors including σ^X and σ^M that transcribe their own structural genes. Here, we report that glucose addition to the medium induced sigX and sigM transcription independent of their anti-σ factors. This induction was dependent on an intracellular acetyl-CoA pool. Transposon mutagenesis searching for the mutants showing no induction of sigX and sigM revealed that the cshA gene encoding DEAD-box RNA helicase is required for gene induction. Global analysis of the acetylome in B. subtilis showed CshA has two acetylated lysine residues. We found that in a cshA mutant with acetylation-abolishing K to R exchange mutations, glucose induction of sigX and sigM was abolished and that glucose addition stimulated acetylation of CshA in the wild type strain. Thus, we present a model wherein glucose addition results in a larger acetyl-CoA pool, probably leading to increased levels of acetylated CshA. CshA is known to associate with RNA polymerase (RNAP), and thus RNAP with acetylated CshA could stimulate the autoregulation of sigX and sigM. This is a unique model showing a functional link between nutritional signals and the basal transcription machinery.

Second Messenger Signaling in Bacillus subtilis: Accumulation of Cyclic di-AMP Inhibits Biofilm Formation

The Gram-positive model organism Bacillus subtilis produces the essential second messenger signaling nucleotide cyclic di-AMP. In B. subtilis and other bacteria, c-di-AMP has been implicated in diverse functions such as control of metabolism, cell division and cell wall synthesis, and potassium transport. To enhance our understanding of the multiple functions of this second messenger, we have studied the consequences of c-di-AMP accumulation at a global level by a transcriptome analysis. C-di-AMP accumulation affected the expression of about 700 genes, among them the two major operons required for biofilm formation. The expression of both operons was severely reduced both in the laboratory and a non-domesticated strain upon accumulation of c-di-AMP. In excellent agreement, the corresponding strain was unable to form complex colonies. In B. subtilis, the transcription factor SinR controls the expression of biofilm genes by binding to their promoter regions resulting in transcription repression. Inactivation of the sinR gene restored biofilm formation even at high intracellular c-di-AMP concentrations suggesting that the second messenger acts upstream of SinR in the signal transduction pathway. As c-di-AMP accumulation did not affect the intracellular levels of SinR, we conclude that the nucleotide affects the activity of SinR.