Organisms from all domains of life use gene regulation networks to control cell growth, identity, function, and responses to environmental challenges. Although accurate global regulatory models would provide critical evolutionary and functional insights, they remain incomplete, even for the best studied organisms. Efforts to build comprehensive networks are confounded by challenges including network scale, degree of connectivity, complexity of organism–environment interactions, and difficulty of estimating the activity of regulatory factors. Taking advantage of the large number of known regulatory interactions in Bacillus subtilis and two transcriptomics datasets (including one with 38 separate experiments collected specifically for this study), we use a new combination of network component analysis and model selection to simultaneously estimate transcription factor activities and learn a substantially expanded transcriptional regulatory network for this bacterium. In total, we predict 2,258 novel regulatory interactions and recall 74% of the previously known interactions. We obtained experimental support for 391 (out of 635 evaluated) novel regulatory edges (62% accuracy), thus significantly increasing our understanding of various cell processes, such as spore formation.

Exponentially growing Bacillus subtilis cultures are epigenetically differentiated into two subpopulations in which cells are either ON or OFF for σD-dependent gene expression: a pattern suggestive of bistability. The gene encoding σD, sigD, is part of the 31-gene fla/che operon where its location at the 3′ end, 25 kb away from the strong Pfla/che promoter, determines its expression level relative to a threshold. Here we show that addition of a single extra copy of the slrA gene in the chromosome inhibited σD-dependent gene expression. SlrA together with SinR and SlrR reduced sigD transcript by potentiating a distance-dependent decrease in fla/che operon transcript abundance that was not mediated by changes in expression from the Pfla/che promoter. Consistent with acting upstream of σD, SlrA/SinR/SlrR could be bypassed by artificial ectopic expression of sigD that hysteretically maintained for 20 generations by engaging the sigD gene at the native locus. SlrA/SinR/SlrR was also bypassed by increasing fla/che transcription and resulted a hypersensitive output in flagellin expression. Thus, flagellin gene expression demonstrated hypersensitivity and hysteresis and we conclude that σD-dependent gene expression is bistable.

Summary Biofilm formation in Bacillus subtilis requires expression of the eps and tapA-sipW-tasA operons to synthesize the extracellular matrix components, extracellular polysaccharide and TasA amyloid proteins, respectively. Expression of both operons is inhibited by the DNA-binding protein master regulator of biofilm formation SinR and activated by the protein RemA. Here we show that RemA is a DNA-binding protein that binds to multiple sites upstream of the promoters of both operons and is both necessary and sufficient for transcriptional activation in vivo and in vitro. We further show that SinR negatively regulates eps operon expression by occluding RemA binding and thus for the Peps promoter SinR functions as an anti-activator. Finally, transcriptional profiling indicated that RemA was primarily a regulator of the extracellular matrix genes, but it also activated genes involved in osmoprotection, leading to the identification of another direct target, the opuA operon.

We describe an algorithm, multi-species cMonkey, for the simultaneous biclustering of heterogeneous multiple-species data collections and apply the algorithm to a group of bacteria containing Bacillus subtilis, Bacillus anthracis, and Listeria monocytogenes. The algorithm reveals evolutionary insights into the surprisingly high degree of conservation of regulatory modules across these three species and allows data and insights from well-studied organisms to complement the analysis of related but less well studied organisms.

Gene expression during spore development in Bacillus subtilis is controlled by cell type-specific RNA polymerase sigma factors. σFand σE control early stages of development in the forespore and the mother cell, respectively. When, at an intermediate stage in development, the mother cell engulfs the forespore, σF is replaced by σG and σE is replaced by σK. The anti-sigma factor CsfB is produced under the control of σF and binds to and inhibits the auto-regulatory σG, but not σF. A position in region 2.1, occupied by an asparagine in σG and by a glutamate in οF, is sufficient for CsfB discrimination of the two sigmas, and allows it to delay the early to late switch in forespore gene expression. We now show that following engulfment completion, csfB is switched on in the mother cell under the control of σK and that CsfB binds to and inhibits σE but not σK, possibly to facilitate the switch from early to late gene expression. We show that a position in region 2.3 occupied by a conserved asparagine in σE and by a conserved glutamate in σK suffices for discrimination by CsfB. We also show that CsfB prevents activation of σG in the mother cell and the premature σG-dependent activation of σK. Thus, CsfB establishes negative feedback loops that curtail the activity of σE and prevent the ectopic activation of σG in the mother cell. The capacity of CsfB to directly block σE activity may also explain how CsfB plays a role as one of the several mechanisms that prevent σE activation in the forespore. Thus the capacity of CsfB to differentiate between the highly similar σF/σG and σE/σK pairs allows it to rinforce the cell-type specificity of these sigma factors and the transition from early to late development in B. subtilis, and possibly in all sporeformers that encode a CsfB orthologue.