Noise in a phosphorelay drives stochastic entry into sporulation in Bacillus subtilis

Altruistic feeding and cell-cell signaling during bacterial differentiation actively enhance phenotypic heterogeneity

Starvation triggers bacterial spore formation, a committed differentiation program that transforms a vegetative cell into a dormant spore. Cells in a population enter sporulation nonuniformly to secure against the possibility that favorable growth conditions, which put sporulation-committed cells at a disadvantage, may resume. This heterogeneous behavior is initiated by a passive mechanism: stochastic activation of a master transcriptional regulator. Here, we identify a cell-cell communication pathway containing the proteins ShfA (YabQ) and ShfP (YvnB) that actively promotes phenotypic heterogeneity, wherein Bacillus subtilis cells that start sporulating early use a calcineurin-like phosphoesterase to release glycerol, which simultaneously acts as a signaling molecule and a nutrient to delay nonsporulating cells from entering sporulation. This produced a more diverse population that was better poised to exploit a sudden influx of nutrients compared to those generating heterogeneity via stochastic gene expression alone. Although conflict systems are prevalent among microbes, genetically encoded cooperative behavior in unicellular organisms can evidently also boost inclusive fitness.

Shaping of microbial phenotypes by trade-offs

Growth rate maximization is an important fitness strategy for microbes. However, the wide distribution of slow-growing oligotrophic microbes in ecosystems suggests that rapid growth is often not favored across ecological environments. In many circumstances, there exist trade-offs between growth and other important traits (e.g., adaptability and survival) due to physiological and proteome constraints. Investments on alternative traits could compromise growth rate and microbes need to adopt bet-hedging strategies to improve fitness in fluctuating environments. Here we review the mechanistic role of trade-offs in controlling bacterial growth and further highlight its ecological implications in driving the emergences of many important ecological phenomena such as co-existence, population heterogeneity and oligotrophic/copiotrophic lifestyles.

Contributions of transcriptional noise to leukaemia evolution: KAT2A as a case-study

Transcriptional noise is proposed to participate in cell fate changes, but contributions to mammalian cell differentiation systems, including cancer, remain associative. Cancer evolution is driven by genetic variability, with modulatory or contributory participation of epigenetic variants. Accumulation of epigenetic variants enhances transcriptional noise, which can facilitate cancer cell fate transitions. Acute myeloid leukaemia (AML) is an aggressive cancer with strong epigenetic dependencies, characterized by blocked differentiation. It constitutes an attractive model to probe links between transcriptional noise and malignant cell fate regulation. Gcn5/KAT2A is a classical epigenetic transcriptional noise regulator. Its loss increases transcriptional noise and modifies cell fates in stem and AML cells. By reviewing the analysis of KAT2A-depleted pre-leukaemia and leukaemia models, I discuss that the net result of transcriptional noise is diversification of cell fates secondary to alternative transcriptional programmes. Cellular diversification can enable or hinder AML progression, respectively, by differentiation of cell types responsive to mutations, or by maladaptation of leukaemia stem cells. KAT2A-dependent noise-responsive genes participate in ribosome biogenesis and KAT2A loss destabilizes translational activity. I discuss putative contributions of perturbed translation to AML biology, and propose KAT2A loss as a model for mechanistic integration of transcriptional and translational control of noise and fate decisions. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Tools and methods for high-throughput single-cell imaging with the mother machine

Despite much progress, image processing remains a significant bottleneck for high-throughput analysis of microscopy data. One popular platform for single-cell time-lapse imaging is the mother machine, which enables long-term tracking of microbial cells under precisely controlled growth conditions. While several mother machine image analysis pipelines have been developed in the past several years, adoption by a non-expert audience remains a challenge. To fill this gap, we implemented our own software, MM3, as a plugin for the multidimensional image viewer napari. napari-MM3 is a complete and modular image analysis pipeline for mother machine data, which takes advantage of the high-level interactivity of napari. Here, we give an overview of napari-MM3 and test it against several well-designed and widely used image analysis pipelines, including BACMMAN and DeLTA. Researchers often analyze mother machine data with custom scripts using varied image analysis methods, but a quantitative comparison of the output of different pipelines has been lacking. To this end, we show that key single-cell physiological parameter correlations and distributions are robust to the choice of analysis method. However, we also find that small changes in thresholding parameters can systematically alter parameters extracted from single-cell imaging experiments. Moreover, we explicitly show that in deep learning-based segmentation, 'what you put is what you get' (WYPIWYG) - that is, pixel-level variation in training data for cell segmentation can propagate to the model output and bias spatial and temporal measurements. Finally, while the primary purpose of this work is to introduce the image analysis software that we have developed over the last decade in our lab, we also provide information for those who want to implement mother machine-based high-throughput imaging and analysis methods in their research.

Dissecting cell heterogeneities in bacterial biofilms and their implications for antibiotic tolerance

Bacterial biofilms consist of large, self-formed aggregates where resident bacteria can exhibit very different physiological states and phenotypes. This heterogeneity of cell types is crucial for many structural and functional emergent properties of biofilms. Consequently, it becomes essential to understand what drives cells to differentiate and how they achieve it within the three-dimensional landscape of the biofilms. Here, we discuss recent advances in comprehending two forms of cell heterogeneity that, while recognized to coexist within biofilms, have proven challenging to distinguish. These two forms include cell heterogeneity arising as a consequence of bacteria physiologically responding to resource gradients formed across the biofilms and cell-to-cell phenotypic heterogeneity, which emerges locally within biofilm subzones among neighboring bacteria due to stochastic variations in gene expression. We describe the defining features and concepts related to both forms of cell heterogeneity and discuss their implications, with a particular focus on antibiotic tolerance.

Altruistic feeding and cell-cell signaling during bacterial differentiation actively enhance phenotypic heterogeneity

Starvation triggers bacterial spore formation, a committed differentiation program that transforms a vegetative cell into a dormant spore. Cells in a population enter sporulation non-uniformly to secure against the possibility that favorable growth conditions, which puts sporulation-committed cells at a disadvantage, may resume. This heterogeneous behavior is initiated by a passive mechanism: stochastic activation of a master transcriptional regulator. Here, we identify a cell-cell communication pathway that actively promotes phenotypic heterogeneity, wherein Bacillus subtilis cells that start sporulating early utilize a calcineurin-like phosphoesterase to release glycerol, which simultaneously acts as a signaling molecule and a nutrient to delay non-sporulating cells from entering sporulation. This produced a more diverse population that was better poised to exploit a sudden influx of nutrients compared to those generating heterogeneity via stochastic gene expression alone. Although conflict systems are prevalent among microbes, genetically encoded cooperative behavior in unicellular organisms can evidently also boost inclusive fitness.

Tools and methods for high-throughput single-cell imaging with the mother machine

Despite much progress, image processing remains a significant bottleneck for high-throughput analysis of microscopy data. One popular platform for single-cell time-lapse imaging is the mother machine, which enables long-term tracking of microbial cells under precisely controlled growth conditions. While several mother machine image analysis pipelines have been developed in the past several years, adoption by a non-expert audience remains a challenge. To fill this gap, we implemented our own software, MM3, as a plugin for the multidimensional image viewer napari. napari-MM3 is a complete and modular image analysis pipeline for mother machine data, which takes advantage of the high-level interactivity of napari. Here, we give an overview of napari-MM3 and test it against several well-designed and widely-used image analysis pipelines, including BACMMAN and DeLTA. Researchers often analyze mother machine data with custom scripts using varied image analysis methods, but a quantitative comparison of the output of different pipelines has been lacking. To this end, we show that key single-cell physiological parameter correlations and distributions are robust to the choice of analysis method. However, we also find that small changes in thresholding parameters can systematically alter parameters extracted from single-cell imaging experiments. Moreover, we explicitly show that in deep learning based segmentation, "what you put is what you get" (WYPIWYG) - i.e., pixel-level variation in training data for cell segmentation can propagate to the model output and bias spatial and temporal measurements. Finally, while the primary purpose of this work is to introduce the image analysis software that we have developed over the last decade in our lab, we also provide information for those who want to implement mother-machine-based high-throughput imaging and analysis methods in their research.

Artificial activation of both σH and Spo0A in Bacillus subtilis enforced initiation of spore development at the vegetatively growing phase

When Bacillus subtilis cells face environmental deterioration, such as exhaustion of nutrients and an increase in cell density, they form spores. It is known that phosphorylation of Spo0A and activation of σH are key events at the initiation of sporulation. However, the initiation of sporulation is an extremely complicated process, and the relationship between these two events remains to be elucidated. To determine the minimum requirements for triggering sporulation initiation, we attempted to induce cell sporulation at the log phase, regardless of nutrients and cell density. In rich media such as Luria-Bertani (LB) medium, the cells of B. subtilis do not sporulate efficiently, possibly because of excess nutrition. When the amount of xylose in the LB medium was limited, σH -dependent transcription of the strain, in which sigA was under the control of the xylose-inducible promoter, was induced, and the frequency of sporulation was elevated according to the decreased level of σA. We also employed a fusion of sad67, which codes for an active form of Spo0A, and the IPTG-inducible promoter. The combination of lowered σA expression and activated Spo0A allowed the cells in the log phase to stop growing and rush into spore development. This observation of enforced initiation of sporulation in the mutant strain was detected even in the presence of the wild-type strain, suggesting that only intracellular events initiate and fulfill spore development regardless of extracellular conditions. Under natural sporulation conditions, the amount of σA did not change drastically throughout growth. Mechanisms that sequester σA from the core RNA polymerase and help σH to become active exist, but this has not yet been elucidated.

Bacterial phenotypic heterogeneity through the lens of single-cell RNA sequencing

Bacterial transcription is not monolithic. Microbes exist in a wide variety of cell states that help them adapt to their environment, acquire and produce essential nutrients, and engage in both competition and cooperation with their neighbors. While we typically think of bacterial adaptation as a group behavior, where all cells respond in unison, there is often a mixture of phenotypic responses within a bacterial population, where distinct cell types arise. A primary phenomenon driving these distinct cell states is transcriptional heterogeneity. Given that bacterial mRNA transcripts are extremely short-lived compared to eukaryotes, their transcriptional state is closely associated with their physiology, and thus the transcriptome of a bacterial cell acts as a snapshot of the behavior of that bacterium. Therefore, the application of single-cell transcriptomics to microbial populations will provide novel insight into cellular differentiation and bacterial ecology. In this review, we provide an overview of transcriptional heterogeneity in microbial systems, discuss the findings already provided by single-cell approaches, and plot new avenues of inquiry in transcriptional regulation, cellular biology, and mechanisms of heterogeneity that are made possible when microbial communities are analyzed at single-cell resolution.

The motility-matrix production switch in Bacillus subtilis-a modeling perspective

Phenotype switching can be triggered by external stimuli and by intrinsic stochasticity. Here, we focus on the motility-matrix production switch in Bacillus subtilis. We use modeling to describe the SinR-SlrR bistable switch and its regulation by SinI and to distinguish different sources of stochasticity. Our simulations indicate that intrinsic fluctuations in the synthesis of SinI are insufficient to drive spontaneous switching and suggest that switching is triggered by upstream noise from the Spo0A phosphorelay. IMPORTANCE The switch from motility to matrix production is the first step toward biofilm formation and, thus, to multicellular behavior in Bacillus subtilis. The transition is governed by a bistable switch based on the interplay of the regulators SinR and SlrR, while SinI transmits upstream signals to that switch. Quantitative modeling can be used to study the switching dynamics. Here, we build such a model step by step to describe the dynamics of the switch and its regulation and to study how spontaneous switching is triggered by upstream noise from the Spo0A phosphorelay.

Identification of potential RapJ hits as sporulation pathway inducer candidates in Bacillus coagulans via structure-based virtual screening and molecular dynamics simulation studies

BACKGROUND

The bacterium Bacillus coagulans has attracted interest because of its ability to produce spores and advantageous probiotic traits, such as facilitating food digestion in the intestine, managing some disorders, and controlling the symbiotic microbiota. Spore-forming probiotic bacteria are especially important in the probiotic industry compared to non-spore-forming bacteria due to their stability during production and high resistance to adverse factors such as stomach acid. When spore-forming bacteria are exposed to environmental stresses, they enter the sporulation pathway to survive. This pathway is activated by the final phosphorylation of the master regulator of spore response, Spo0A, and upon achieving the phosphorylation threshold. Spo0A is indirectly inhibited by some enzymes of the aspartate response regulator phosphatase (Rap) family, such as RapJ. RapJ is one of the most important Rap enzymes in the sporogenesis pathway, which is naturally inhibited by the pentapeptides.

METHODS

This study used structure-based virtual screening and molecular dynamics (MD) simulation studies to find potential RapJ hits that could induce the sporulation pathway. The crystal structures of RapJ complexed with pentapeptide clearly elucidated their interactions with the enzyme active site.

RESULTS

Based on the binding compartment, through molecular docking, MD simulation, hydrogen bonds, and binding-free energy calculations, a series of novel hits against RapJ named tandutinib, infigratinib, sitravatinib, linifanib, epertinib, surufatinib, and acarbose were identified. Among these compounds, acarbose obtained the highest score, especially in terms of the number of hydrogen bonds, which plays a major role in stabilizing RapJ-ligand complexes, and also according to the occupancy percentages of hydrogen bonds, its hydrogen bonds were more stable during the simulation time. Consequently, acarbose is probably the most suitable hit for RapJ enzyme. Notably, experimental validation is crucial to confirm the effectiveness of the selected ligands.

Microfluidics for long-term single-cell time-lapse microscopy: Advances and applications

Cells are inherently dynamic, whether they are responding to environmental conditions or simply at equilibrium, with biomolecules constantly being made and destroyed. Due to their small volumes, the chemical reactions inside cells are stochastic, such that genetically identical cells display heterogeneous behaviors and gene expression profiles. Studying these dynamic processes is challenging, but the development of microfluidic methods enabling the tracking of individual prokaryotic cells with microscopy over long time periods under controlled growth conditions has led to many discoveries. This review focuses on the recent developments of one such microfluidic device nicknamed the mother machine. We overview the original device design, experimental setup, and challenges associated with this platform. We then describe recent methods for analyzing experiments using automated image segmentation and tracking. We further discuss modifications to the experimental setup that allow for time-varying environmental control, replicating batch culture conditions, cell screening based on their dynamic behaviors, and to accommodate a variety of microbial species. Finally, this review highlights the discoveries enabled by this technology in diverse fields, such as cell-size control, genetic mutations, cellular aging, and synthetic biology.

Spores-on-a-chip: new frontiers for spore research

In recent years, microfluidic technologies have become widespread in biological science. However, the suitability of this technique for understanding different aspects of spore research has hardly been considered. Herein, we review recent developments in 'spores-on-a-chip' technologies, highlighting how they could be exploited to drive new frontiers in spore research.

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.

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

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

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

Proceedings from the 3rd International Conference on Microbiome Engineering

The human microbiome has been inextricably linked to multiple facets of human physiology. From an engineering standpoint, the ability to precisely control the composition and activity of the microbiome holds great promise for furthering our understanding of disease etiology and for new avenues of therapeutic and diagnostic agents. While the field of microbiome research is still in its infancy, growing engineering efforts are emerging to enable new studies in the microbiome and to rapidly translate these findings to microbiome‐based interventions. At the 3rd International Conference on Microbiome Engineering, leading experts in the field presented state‐of‐the‐art work in microbiome engineering, discussing probiotics, prebiotics, engineered microbes, microbially derived biomolecules, and bacteriophage.

Learning from Nature: Bacterial Spores as a Target for Current Technologies in Medicine (Review)

The capability of some representatives of Clostridium spp. and Bacillus spp. genera to form spores in extreme external conditions long ago became a subject of medico-biological investigations. Bacterial spores represent dormant cellular forms of gram-positive bacteria possessing a high potential of stability and the capability to endure extreme conditions of their habitat. Owing to these properties, bacterial spores are recognized as the most stable systems on the planet, and spore-forming microorganisms became widely spread in various ecosystems.

Spore-forming bacteria have been attracted increased interest for years due to their epidemiological danger. Bacterial spores may be in the quiescent state for dozens or hundreds of years but after they appear in the favorable conditions of a human or animal organism, they turn into vegetative forms causing an infectious process. The greatest threat among the pathogenic spore-forming bacteria is posed by the causative agents of anthrax (B. anthracis), food toxicoinfection (B. cereus), pseudomembranous colitis (C. difficile), botulism (C. botulinum), gas gangrene (C. perfringens).

For the effective prevention of severe infectious diseases first of all it is necessary to study the molecular structure of bacterial spores and the biochemical mechanisms of sporulation and to develop innovative methods of detection and disinfection of dormant cells. There is another side of the problem: the necessity to investigate exo- and endospores from the standpoint of obtaining similar artificially synthesized models in order to use them in the latest medical technologies for the development of thermostable vaccines, delivery of biologically active substances to the tissues and intracellular structures. In recent years, bacterial spores have become an interesting object for the exploration from the point of view of a new paradigm of unicellular microbiology in order to study microbial heterogeneity by means of the modern analytical tools.

Effects of orphan histidine kinases on clostridial sporulation progression and metabolism

Solventogenesis and sporulation of clostridia are the main responsive adaptations to the acidic environment during acetone-butanol-ethanol (ABE) fermentation. It was hypothesized that five orphan histidine kinases (HKs) including Cac3319, Cac0323, Cac0903, Cac2730, and Cac0437 determined the cell fates between sporulation and solventogenesis. In this study, the comparative genomic analysis revealed that a mutation in cac0437 appeared to contribute to the nonsporulating feature of ATCC 55025. Hence, the individual and interactive roles of five HKs in regulating cell growth, metabolism, and sporulation were investigated. The fermentation results of mutants with different HK expression levels suggested that cac3319 and cac0437 played critical roles in regulating sporulation and acids and butanol biosynthesis. Morphological analysis revealed that cac3319 knockout abolished sporulation (Stage 0) whereas cac3319 overexpression promoted spore development (Stage VII), and cac0437 knockout initiated but blocked sporulation before Stage II, indicating the progression of sporulation was altered through engineering HKs. By combinatorial HKs knockout, the interactive effects between two different HKs were investigated. This study elucidated the regulatory roles of HKs in clostridial differentiation and demonstrated that HK engineering can be effectively used to control sporulation and enhance butanol biosynthesis.

Bacillus subtilis: a bacterium for all seasons

Microbes inhabit ever-changing environments in which the conditions for propagation, whether they be favorable or unfavorable, are often transient. Plants and animals are able to cope with changing circumstances in their external environment by maintaining homeostasis for internal cells, tissues and organs. Microbes do not have this luxury. Instead, they have evolved mechanisms for adapting to change by being versatile. In some cases, these adaptations involve the expression of a limited set of genes that exploit the availability of a food source or provide protection against chemical, radiation or thermal stress. In other cases, however, microbes adapt in a more elaborate manner by entering a specialized state that enables them to exploit a particular niche or protects them against environmental extremes. Entry into such a state can occur as a direct response to an external cue or stochastically (by chance) as part of a bet-hedging strategy, or sometimes a combination of the two. An example of a versatile microbe that exhibits a wide variety of states is the bacterium Bacillus subtilis, the subject of this Primer. Many of the states exhibited by B. subtilis are similar to states observed in other bacteria. What is special about B. subtilis is the unusually rich repertoire of alternative states exhibited by one bacterium, enabling it to cope with a wide range of environmental challenges.

Cell-to-Cell Heterogeneity in Trypanosomes

Recent developments in single-cell and single-molecule techniques have revealed surprising levels of heterogeneity among isogenic cells. These advances have transformed the study of cell-to-cell heterogeneity into a major area of biomedical research, revealing that it can confer essential advantages, such as priming populations of unicellular organisms for future environmental stresses. Protozoan parasites, such as trypanosomes, face multiple and often hostile environments, and to survive, they undergo multiple changes, including changes in morphology, gene expression, and metabolism. But why does only a subset of proliferative cells differentiate to the next life cycle stage? Why do only some bloodstream parasites undergo antigenic switching while others stably express one variant surface glycoprotein? And why do some parasites invade an organ while others remain in the bloodstream? Building on extensive research performed in bacteria, here we suggest that biological noise can contribute to the fitness of eukaryotic pathogens and discuss the importance of cell-to-cell heterogeneity in trypanosome infections. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Valuing what happens: a biogenic approach to valence and (potentially) affect

Valence is half of the pair of properties that constitute core affect, the foundation of emotion. But what is valence, and where is it found in the natural world? Currently, this question cannot be answered. The idea that emotion is the body's way of driving the organism to secure its survival, thriving and reproduction runs like a leitmotif from the pathfinding work of Antonio Damasio through four book-length neuroscientific accounts of emotion recently published by the field's leading practitioners. Yet while Damasio concluded 20 years ago that the homeostasis-affect linkage is rooted in unicellular life, no agreement exists about whether even non-human animals with brains experience emotions. Simple neural animals-those less brainy than bees, fruit flies and other charismatic invertebrates-are not even on the radar of contemporary affective research, to say nothing of aneural organisms. This near-sightedness has effectively denied the most productive method available for getting a grip on highly complex biological processes to a scientific domain whose importance for understanding biological decision-making cannot be underestimated. Valence arguably is the fulcrum around which the dance of life revolves. Without the ability to discriminate advantage from harm, life very quickly comes to an end. In this paper, we review the concept of valence, where it came from, the work it does in current leading theories of emotion, and some of the odd features revealed via experiment. We present a biologically grounded framework for investigating valence in any organism and sketch a preliminary pathway to a computational model. This article is part of the theme issue 'Basal cognition: conceptual tools and the view from the single cell'.

Microbial single-cell RNA sequencing by split-pool barcoding

Single-cell RNA sequencing (scRNA-seq) has become an essential tool for characterizing gene expression in eukaryotes, but current methods are incompatible with bacteria. Here, we introduce microSPLiT (microbial split-pool ligation transcriptomics), a high-throughput scRNA-seq method for Gram-negative and Gram-positive bacteria that can resolve heterogeneous transcriptional states. We applied microSPLiT to >25,000 Bacillus subtilis cells sampled at different growth stages, creating an atlas of changes in metabolism and lifestyle. We retrieved detailed gene expression profiles associated with known, but rare, states such as competence and prophage induction and also identified unexpected gene expression states, including the heterogeneous activation of a niche metabolic pathway in a subpopulation of cells. MicroSPLiT paves the way to high-throughput analysis of gene expression in bacterial communities that are otherwise not amenable to single-cell analysis, such as natural microbiota.

Machine learning uncovers independently regulated modules in the Bacillus subtilis transcriptome

The transcriptional regulatory network (TRN) of Bacillus subtilis coordinates cellular functions of fundamental interest, including metabolism, biofilm formation, and sporulation. Here, we use unsupervised machine learning to modularize the transcriptome and quantitatively describe regulatory activity under diverse conditions, creating an unbiased summary of gene expression. We obtain 83 independently modulated gene sets that explain most of the variance in expression and demonstrate that 76% of them represent the effects of known regulators. The TRN structure and its condition-dependent activity uncover putative or recently discovered roles for at least five regulons, such as a relationship between histidine utilization and quorum sensing. The TRN also facilitates quantification of population-level sporulation states. As this TRN covers the majority of the transcriptome and concisely characterizes the global expression state, it could inform research on nearly every aspect of transcriptional regulation in B. subtilis.

The systems-level regulatory structure underlying gene expression in bacteria can be inferred using machine learning algorithms. Here we show this structure for Bacillus subtilis, present five hypotheses gleaned from it, and analyse the process of sporulation from its perspective.

Efficient and cost-effective bacterial mRNA sequencing from low input samples through ribosomal RNA depletion

BACKGROUND

RNA sequencing is a powerful approach to quantify the genome-wide distribution of mRNA molecules in a population to gain deeper understanding of cellular functions and phenotypes. However, unlike eukaryotic cells, mRNA sequencing of bacterial samples is more challenging due to the absence of a poly-A tail that typically enables efficient capture and enrichment of mRNA from the abundant rRNA molecules in a cell. Moreover, bacterial cells frequently contain 100-fold lower quantities of RNA compared to mammalian cells, which further complicates mRNA sequencing from non-cultivable and non-model bacterial species. To overcome these limitations, we report EMBR-seq (Enrichment of mRNA by Blocked rRNA), a method that efficiently depletes 5S, 16S and 23S rRNA using blocking primers to prevent their amplification.

RESULTS

EMBR-seq results in 90% of the sequenced RNA molecules from an E. coli culture deriving from mRNA. We demonstrate that this increased efficiency provides a deeper view of the transcriptome without introducing technical amplification-induced biases. Moreover, compared to recent methods that employ a large array of oligonucleotides to deplete rRNA, EMBR-seq uses a single or a few oligonucleotides per rRNA, thereby making this new technology significantly more cost-effective, especially when applied to varied bacterial species. Finally, compared to existing commercial kits for bacterial rRNA depletion, we show that EMBR-seq can be used to successfully quantify the transcriptome from more than 500-fold lower starting total RNA.

CONCLUSIONS

EMBR-seq provides an efficient and cost-effective approach to quantify global gene expression profiles from low input bacterial samples.

Functional roles of microbial cell-to-cell heterogeneity and emerging technologies for analysis and control

Clonal cell populations often display significant cell-to-cell phenotypic heterogeneity, even when maintained under constant external conditions. This variability can result from the inherently stochastic nature of transcription and translation processes, which leads to varying numbers of transcripts and proteins per cell. Here, we showcase studies that reveal links between stochastic cellular events and biological functions in isogenic microbial populations. Then, we highlight emerging tools from engineering, computation, and synthetic and molecular biology that enable precise measurement, control, and analysis of gene expression noise in microorganisms. The capabilities offered by this sophisticated toolbox will shape future directions in the field and generate insight into the behavior of living systems at the single-cell level.

Abhyankar W, N. Swarge B, van der Wel N, Kramer G, Brul S, J. de Koning L. Artificial Sporulation Induction (ASI) by kinA Overexpression Affects the Proteomes and Properties of Bacillus subtilis Spores

To facilitate more accurate spore proteomic analysis, the current study focuses on inducing homogeneous sporulation by overexpressing kinA and assesses the effect of synchronized sporulation initiation on spore resistance, structures, the germination behavior at single-spore level and the proteome. The results indicate that, in our set up, the sporulation by overexpressing kinA can generate a spore yield of 70% within 8 h. The procedure increases spore wet heat resistance and thickness of the spore coat and cortex layers, whilst delaying the time to spore phase-darkening and burst after addition of germinant. The proteome analysis reveals that the upregulated proteins in the kinA induced spores, compared to spores without kinA induction, as well as the 'wildtype' spores, are mostly involved in spore formation. The downregulated proteins mostly belong to the categories of coping with stress, carbon and nitrogen metabolism, as well as the regulation of sporulation. Thus, while kinA overexpression enhances synchronicity in sporulation initiation, it also has profound effects on the central equilibrium of spore formation and spore germination, through modulation of the spore molecular composition and stress resistance physiology.

Stochastic pulsing of gene expression enables the generation of spatial patterns in Bacillus subtilis biofilms

Stochastic pulsing of gene expression can generate phenotypic diversity in a genetically identical population of cells, but it is unclear whether it has a role in the development of multicellular systems. Here, we show how stochastic pulsing of gene expression enables spatial patterns to form in a model multicellular system, Bacillus subtilis bacterial biofilms. We use quantitative microscopy and time-lapse imaging to observe pulses in the activity of the general stress response sigma factor σ^B in individual cells during biofilm development. Both σ^B and sporulation activity increase in a gradient, peaking at the top of the biofilm, even though σ^B represses sporulation. As predicted by a simple mathematical model, increasing σ^B expression shifts the peak of sporulation to the middle of the biofilm. Our results demonstrate how stochastic pulsing of gene expression can play a key role in pattern formation during biofilm development.

Stochastic pulsing of gene expression can generate phenotypic diversity in a genetically identical population of cells. Here, the authors show that stochastic pulsing in the expression of a sigma factor enables the formation of spatial patterns in a multicellular system, Bacillus subtilis bacterial biofilms.

Stochastic Expression of Sae-Dependent Virulence Genes during Staphylococcus aureus Biofilm Development Is Dependent on SaeS

The intricate process of biofilm formation in the human pathogen Staphylococcus aureus involves distinct stages during which a complex mixture of matrix molecules is produced and modified throughout the developmental cycle. Early in biofilm development, a subpopulation of cells detaches from its substrate in an event termed "exodus" that is mediated by SaePQRS-dependent stochastic expression of a secreted staphylococcal nuclease, which degrades extracellular DNA within the matrix, causing the release of cells and subsequently allowing for the formation of metabolically heterogenous microcolonies. Since the SaePQRS regulatory system is involved in the transcriptional control of multiple S. aureus virulence factors, the expression of several additional virulence genes was examined within a developing biofilm by introducing fluorescent gene reporter plasmids into wild-type S. aureus and isogenic regulatory mutants and growing these strains in a microfluidic system that supplies the bacteria with a constant flow of media while simultaneously imaging developing biofilms in 5-min intervals. This study demonstrated that multiple virulence genes, including nuc, were expressed stochastically within a specialized subpopulation of cells in nascent biofilms. We demonstrated that virulence genes regulated by SaePQRS were stochastically expressed in nearly all strains examined whereas Agr-regulated genes were expressed more homogenously within maturing microcolonies. The commonly used Newman strain contains a variant of SaeS (SaeS^P) that confers constitutive kinase activity to the protein and caused this strain to lack the stochastic expression pattern observed in other strain backgrounds. Importantly, repair of the SaeS^P allele resulting in reversion to the well-conserved SaeS ^L allele found in other strains restored stochastic expression in this strain.IMPORTANCE Staphylococcus aureus is an important human pathogen capable of colonizing diverse tissue types and inducing severe disease in both immunocompromised and otherwise healthy individuals. Biofilm infections caused by this bacterial species are of particular concern because of their persistence, even in the face of intensive therapeutic intervention. The results of the current study demonstrate the stochastic nature of Sae-mediated virulence gene expression in S. aureus and indicate that this regulatory system may function as a "bistable switch" in a manner similar to that seen with regulators controlling competence gene expression in Bacillus subtilis and persister cell formation in Escherichia coli The results of this study provide a new perspective on the complex mechanisms utilized by S. aureus during the establishment of infections.

Sporulation and Germination in Clostridial Pathogens

As obligate anaerobes, clostridial pathogens depend on their metabolically dormant, oxygen-tolerant spore form to transmit disease. However, the molecular mechanisms by which those spores germinate to initiate infection and then form new spores to transmit infection remain poorly understood. While sporulation and germination have been well characterized in Bacillus subtilis and Bacillus anthracis, striking differences in the regulation of these processes have been observed between the bacilli and the clostridia, with even some conserved proteins exhibiting differences in their requirements and functions. Here, we review our current understanding of how clostridial pathogens, specifically Clostridium perfringens, Clostridium botulinum, and Clostridioides difficile, induce sporulation in response to environmental cues, assemble resistant spores, and germinate metabolically dormant spores in response to environmental cues. We also discuss the direct relationship between toxin production and spore formation in these pathogens.

High-throughput detection and tracking of cells and intracellular spots in mother machine experiments

The analysis of bacteria at the single-cell level is essential to characterization of processes in which cellular heterogeneity plays an important role. BACMMAN (bacteria mother machine analysis) is a software allowing fast and reliable automated image analysis of high-throughput 2D or 3D time-series images from experiments using the 'mother machine', a very popular microfluidic device allowing biological processes in bacteria to be investigated at the single-cell level. Here, we describe how to use some of the BACMMAN features, including (i) segmentation and tracking of bacteria and intracellular fluorescent spots, (ii) visualization and editing of the results, (iii) configuration of the image-processing pipeline for different datasets and (iv) BACMMAN coupling to data analysis software for visualization and analysis of data subsets with specific properties. Among software specifically dedicated to the analysis of mother machine data, only BACMMAN allows segmentation and tracking of both bacteria and intracellular spots. For a single position, single channel with 1,000 frames (2-GB dataset), image processing takes ~6 min on a regular computer. Numerous implemented algorithms, easy configuration and high modularity ensure wide applicability of the BACMMAN software.

Challenges in measuring and understanding biological noise

Biochemical reactions are intrinsically stochastic, leading to variation in the production of mRNAs and proteins within cells. In the scientific literature, this source of variation is typically referred to as 'noise'. The observed variability in molecular phenotypes arises from a combination of processes that amplify and attenuate noise. Our ability to quantify cell-to-cell variability in numerous biological contexts has been revolutionized by recent advances in single-cell technology, from imaging approaches through to 'omics' strategies. However, defining, accurately measuring and disentangling the stochastic and deterministic components of cell-to-cell variability is challenging. In this Review, we discuss the sources, impact and function of molecular phenotypic variability and highlight future directions to understand its role.

Microbial life cycles link global modularity in regulation to mosaic evolution

Microbes are exposed to changing environments, to which they can respond by adopting various lifestyles such as swimming, colony formation or dormancy. These lifestyles are often studied in isolation, thereby giving a fragmented view of the life cycle as a whole. Here, we study lifestyles in the context of this whole. We first use machine learning to reconstruct the expression changes underlying life cycle progression in the bacterium Bacillus subtilis, based on hundreds of previously acquired expression profiles. This yields a timeline that reveals the modular organization of the life cycle. By analysing over 380 Bacillales genomes, we then show that life cycle modularity gives rise to mosaic evolution in which life stages such as motility and sporulation are conserved and lost as discrete units. We postulate that this mosaic conservation pattern results from habitat changes that make these life stages obsolete or detrimental. Indeed, when evolving eight distinct Bacillales strains and species under laboratory conditions that favour colony growth, we observe rapid and parallel losses of the sporulation life stage across species, induced by mutations that affect the same global regulator. We conclude that a life cycle perspective is pivotal to understanding the causes and consequences of modularity in both regulation and evolution.

Comparative genomics and proteogenomics highlight key molecular players involved in Frankia sporulation

Sporulation is a microbial adaptive strategy to resist inhospitable conditions for vegetative growth and to disperse to colonise more favourable environments. This microbial trait is widespread in Actinobacteria. Among them, Frankia strains are able to differentiate sporangia in pure culture, while others can sporulate even when in symbiosis with sporulation occurring within host cells. The molecular determinants controlling Frankia sporulation have not been yet described. In order to highlight, for the first time, the molecular players potentially involved in Frankia sporulation, we conducted (i) a comparison of protein contents between Frankia spores and hyphae and (ii) a comparative genomic analysis of Frankia proteomes with sporulating and non-sporulating Actinobacteria. Among the main results, glycogen-metabolism related proteins, as well as oxidative stress response and protease-like proteins were overdetected in hyphae, recalling lytic processes that allow Streptomyces cells to erect sporogenic hyphae. Several genes encoding transcriptional regulators, including GntR-like, appeared up-regulated in spores, as well as tyrosinase, suggesting their potential role in mature spore metabolism. Finally, our results highlighted new proteins potentially involved in Frankia sporulation, including a pyrophosphate-energized proton pump and YaaT, described as involved in the phosphorelay allowing sporulation in Bacillus subtilis, leading us to discuss the role of a phosphorelay in Frankia sporulation.

Stress-induced protein aggregates shape population heterogeneity in bacteria

The concept of phenotypic heterogeneity preparing a subpopulation of isogenic cells to better cope with anticipated stresses has been well established. However, less is known about how stress itself can drive subsequent cellular individualization in clonal populations. In this perspective, we focus on the impact of stress-induced cellular protein aggregates, and how their segregation and disaggregation can act as a deterministic incentive for heterogeneity in the population emerging from a stressed ancestor.

Metal removal and morphological changes of B. megaterium in the presence of a spent catalyst

Spent catalysts represent an environmental concern, mainly due to their elevated metal content. Although conventional treatment methods for spent catalysts are available, they generate large volumes of potentially harmful wastes and gaseous emissions. To overcome the environmental impact, biotechnological approaches are currently being explored and developed. Thus, the current study assayed the capability of Bacillus megaterium strain MNSH1-9K-1 to remove Al, Ni, V and Ti contained in the spent catalyst coded as ECAT-TL-II. To this end, B. megaterium MNSH1-9K-1 growth and metal uptake abilities in the presence of ECAT-TL-II spent catalyst at 15% (wt/vol) pulp density were evaluated in modified Starkey medium at 37 °C and 200 rpm. The results presented here show B. megaterium resistance capability to the high-metal content residue, and its Al, V and Ni removal ability, in 1,059.15 ± 197.28 mg kg^-1 of Al, 43.39 ± 24.13 mg kg^-1 of V and 0.58 ± 0.00 mg kg^-1 of Ni, corresponding to the 0.79%, 1.63% and 0.46% of each metal content, respectively, while no Ti removal was detected. Besides, it was observed that the sporulation process took place in B. megaterium cells in the presence of the spent catalyst. The results shown in this study suggest the potential of the strain MNSH1-9K-1 for the removal of metals contained in high-metal content residues, contributing also to the knowledge of the metal resistance and removal abilities of B. megaterium in the presence of a spent catalyst, and how morphological cell changes may be occurring while metal removal is taking place.

In-depth analysis of Bacillus subtilis proteome identifies new ORFs and traces the evolutionary history of modified proteins

Bacillus subtilis is a sporulating Gram-positive bacterium widely used in basic research and biotechnology. Despite being one of the best-characterized bacterial model organism, recent proteomics studies identified only about 50% of its theoretical protein count. Here we combined several hundred MS measurements to obtain a comprehensive map of the proteome, phosphoproteome and acetylome of B. subtilis grown at 37 °C in minimal medium. We covered 75% of the theoretical proteome (3,159 proteins), detected 1,085 phosphorylation and 4,893 lysine acetylation sites and performed a systematic bioinformatic characterization of the obtained data. A subset of analyzed MS files allowed us to reconstruct a network of Hanks-type protein kinases, Ser/Thr/Tyr phosphatases and their substrates. We applied genomic phylostratigraphy to gauge the evolutionary age of B. subtilis protein classes and revealed that protein modifications were present on the oldest bacterial proteins. Finally, we performed a proteogenomic analysis by mapping all MS spectra onto a six-frame translation of B. subtilis genome and found evidence for 19 novel ORFs. We provide the most extensive overview of the proteome and post-translational modifications for B. subtilis to date, with insights into functional annotation and evolutionary aspects of the B. subtilis genome.

Single-cell Microfluidic Analysis of Bacillus subtilis

Microfluidic technology overcomes many of the limitations to traditional analytical methods in microbiology. Unlike bulk-culture methods, it offers single-cell resolution and long observation times spanning hundreds of generations; unlike agarose pad-based microscopy, it has uniform growth conditions that can be tightly controlled. Because the continuous flow of growth medium isolates the cells in a microfluidic device from unpredictable variations in the local chemical environment caused by cell growth and metabolism, authentic changes in gene expression and cell growth in response to specific stimuli can be more confidently observed. Bacillus subtilis is used here as a model bacterial species to demonstrate a "mother machine"-type method for cellular analysis. We show how to construct and plumb a microfluidic device, load it with cells, initiate microscopic imaging, and expose cells to a stimulus by switching from one growth medium to another. A stress-responsive reporter is used as an example to reveal the type of data that may be obtained by this method. We also briefly discuss further applications of this method for other types of experiments, such as analysis of bacterial sporulation.

Transcription termination factor Rho and microbial phenotypic heterogeneity

Populations of genetically identical microorganisms exhibit high degree of cell-to-cell phenotypic diversity even when grown in uniform environmental conditions. Heterogeneity is a genetically determined trait, which ensures bacterial adaptation and survival in the ever changing environmental conditions. Fluctuations in gene expression (noise) at the level of transcription initiation largely contribute to cell-to-cell variability within population. Not surprisingly, the analyses of the mechanisms driving phenotypic heterogeneity are mainly focused on the activity of promoters and transcriptional factors. Less attention is currently given to a role of intrinsic and factor-dependent transcription terminators. Here, we discuss recent advances in understanding the regulatory role of the multi-functional transcription termination factor Rho, the major inhibitor of pervasive transcription in bacteria and the emerging global regulator of gene expression. We propose that termination activity of Rho might be among the mechanisms by which cells manage the intensity of transcriptional noise, thus affecting population heterogeneity.

Noise in a phosphorelay drives stochastic entry into sporulation in Bacillus subtilis

Entry into sporulation in Bacillus subtilis is governed by a phosphorelay in which phosphoryl groups from a histidine kinase are successively transferred via relay proteins to the response regulator Spo0A. Spo0A~P, in turn, sets in motion events that lead to asymmetric division and activation of the cell-specific transcription factor σ^F, a hallmark for entry into sporulation. Here, we have used a microfluidics-based platform to investigate the activation of Spo0A and σ^F in individual cells held under constant, sporulation-inducing conditions. The principal conclusions were that: (i) activation of σ^F occurs with an approximately constant probability after adaptation to conditions of nutrient limitation; (ii) activation of σ^F is tightly correlated with, and preceded by, Spo0A~P reaching a high threshold level; (iii) activation of Spo0A takes place abruptly just prior to asymmetric division; and (iv) the primary source of noise in the activation of Spo0A is the phosphorelay. We propose that cells exhibit a constant probability of attaining a high threshold level of Spo0A~P due to fluctuations in the flux of phosphoryl groups through the phosphorelay.