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

High-throughput single-cell transcriptomics of bacteria using combinatorial barcoding

Microbial split-pool ligation transcriptomics (microSPLiT) is a high-throughput single-cell RNA sequencing method for bacteria. With four combinatorial barcoding rounds, microSPLiT can profile transcriptional states in hundreds of thousands of Gram-negative and Gram-positive bacteria in a single experiment without specialized equipment. As bacterial samples are fixed and permeabilized before barcoding, they can be collected and stored ahead of time. During the first barcoding round, the fixed and permeabilized bacteria are distributed into a 96-well plate, where their transcripts are reverse transcribed into cDNA and labeled with the first well-specific barcode inside the cells. The cells are mixed and redistributed two more times into new 96-well plates, where the second and third barcodes are appended to the cDNA via in-cell ligation reactions. Finally, the cells are mixed and divided into aliquot sub-libraries, which can be stored until future use or prepared for sequencing with the addition of a fourth barcode. It takes 4 days to generate sequencing-ready libraries, including 1 day for collection and overnight fixation of samples. The standard plate setup enables single-cell transcriptional profiling of up to 1 million bacterial cells and up to 96 samples in a single barcoding experiment, with the possibility of expansion by adding barcoding rounds. The protocol requires experience in basic molecular biology techniques, handling of bacterial samples and preparation of DNA libraries for next-generation sequencing. It can be performed by experienced undergraduate or graduate students. Data analysis requires access to computing resources, familiarity with Unix command line and basic experience with Python or R.

Information content and optimization of self-organized developmental systems

A key feature of many developmental systems is their ability to self-organize spatial patterns of functionally distinct cell fates. To ensure proper biological function, such patterns must be established reproducibly, by controlling and even harnessing intrinsic and extrinsic fluctuations. While the relevant molecular processes are increasingly well understood, we lack a principled framework to quantify the performance of such stochastic self-organizing systems. To that end, we introduce an information-theoretic measure for self-organized fate specification during embryonic development. We show that the proposed measure assesses the total information content of fate patterns and decomposes it into interpretable contributions corresponding to the positional and correlational information. By optimizing the proposed measure, our framework provides a normative theory for developmental circuits, which we demonstrate on lateral inhibition, cell type proportioning, and reaction-diffusion models of self-organization. This paves a way toward a classification of developmental systems based on a common information-theoretic language, thereby organizing the zoo of implicated chemical and mechanical signaling processes.

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.

Perception of butenolides by Bacillus subtilis via the α/β hydrolase RsbQ

The regulation of behavioral and developmental decisions by small molecules is common to all domains of life. In plants, strigolactones and karrikins are butenolide growth regulators that influence several aspects of plant growth and development, as well as interactions with symbiotic fungi.1,2,3 DWARF14 (D14) and KARRIKIN INSENSITIVE2 (KAI2) are homologous enzyme-receptors that perceive strigolactones and karrikins, respectively, and that require hydrolase activity to effect signal transduction.4,5,6,7 RsbQ, a homolog of D14 and KAI2 from the gram-positive bacterium Bacillus subtilis, regulates growth responses to nutritional stress via the alternative transcription factor SigmaB (σB).8,9 However, the molecular function of RsbQ is unknown. Here, we show that RsbQ perceives butenolide compounds that are bioactive in plants. RsbQ is thermally destabilized by the synthetic strigolactone GR24 and its desmethyl butenolide equivalent dGR24. We show that, like D14 and KAI2, RsbQ is a functional butenolide hydrolase that undergoes covalent modification of the catalytic histidine residue. Exogenous application of both GR24 and dGR24 inhibited the endogenous signaling function of RsbQ in vivo, with dGR24 being 10-fold more potent. Application of dGR24 to B. subtilis phenocopied loss-of-function rsbQ mutations and led to a significant downregulation of σB-regulated transcripts. We also discovered that exogenous butenolides promoted the transition from planktonic to biofilm growth. Our results suggest that butenolides may serve as inter-kingdom signaling compounds between plants and bacteria to help shape rhizosphere communities.

The σB alternative sigma factor circuit modulates noise to generate different types of pulsing dynamics

Single-cell approaches are revealing a high degree of heterogeneity, or noise, in gene expression in isogenic bacteria. How gene circuits modulate this noise in gene expression to generate robust output dynamics is unclear. Here we use the Bacillus subtilis alternative sigma factor σB as a model system for understanding the role of noise in generating circuit output dynamics. σB controls the general stress response in B. subtilis and is activated by a range of energy and environmental stresses. Recent single-cell studies have revealed that the circuit can generate two distinct outputs, stochastic pulsing and a single pulse response, but the conditions under which each response is generated are under debate. We implement a stochastic mathematical model of the σB circuit to investigate this and find that the system's core circuit can generate both response types. This is despite one response (stochastic pulsing) being stochastic in nature, and the other (single response pulse) being deterministic. We demonstrate that the main determinant for whichever response is generated is the degree with which the input pathway activates the core circuit, although the noise properties of the input pathway also biases the system towards one or the other type of output. Thus, our work shows how stochastic modelling can reveal the mechanisms behind non-intuitive gene circuit output dynamics.

Spatial transcriptome uncovers rich coordination of metabolism in E. coli K12 biofilm

Microbial communities often display region-specific properties, which give rise to complex interactions and emergent behaviors that are critical to the homeostasis and stress response of the communities. However, systems-level understanding of these properties still remains elusive. In this study, we established RAINBOW-seq and profiled the transcriptome of Escherichia coli biofilm communities with high spatial resolution and high gene coverage. We uncovered three modes of community-level coordination, including cross-regional resource allocation, local cycling and feedback signaling, which were mediated by strengthened transmembrane transport and spatially specific activation of metabolism. As a consequence of such coordination, the nutrient-limited region of the community maintained an unexpectedly high level of metabolism, enabling it to express many signaling genes and functionally unknown genes with potential sociality functions. Our work provides an extended understanding of the metabolic interplay in biofilms and presents a new approach of investigating complex interactions in bacterial communities on the systems level.

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.

Advances in linking single-cell bacterial stress response to population-level survival

Stress response mechanisms can allow bacteria to survive a myriad of challenges, including nutrient changes, antibiotic encounters, and antagonistic interactions with other microbes. Expression of these stress response pathways, in addition to other cell features such as growth rate and metabolic state, can be heterogeneous across cells and over time. Collectively, these single-cell-level phenotypes contribute to an overall population-level response to stress. These include diversifying actions, which can be used to enable bet-hedging, and coordinated actions, such as biofilm production, horizontal gene transfer, and cross-feeding. Here, we highlight recent results and emerging technologies focused on both single-cell and population-level responses to stressors, and we draw connections about the combined impact of these effects on survival of bacterial communities.

Navigating Environmental Transitions: the Role of Phenotypic Variation in Bacterial Responses

The ability of bacteria to respond to changes in their environment is critical to their survival, allowing them to withstand stress, form complex communities, and induce virulence responses during host infection. A remarkable feature of many of these bacterial responses is that they are often variable across individual cells, despite occurring in an isogenic population exposed to a homogeneous environmental change, a phenomenon known as phenotypic heterogeneity. Phenotypic heterogeneity can enable bet-hedging or division of labor strategies that allow bacteria to survive fluctuating conditions. Investigating the significance of phenotypic heterogeneity in environmental transitions requires dynamic, single-cell data. Technical advances in quantitative single-cell measurements, imaging, and microfluidics have led to a surge of publications on this topic. Here, we review recent discoveries on single-cell bacterial responses to environmental transitions of various origins and complexities, from simple diauxic shifts to community behaviors in biofilm formation to virulence regulation during infection. We describe how these studies firmly establish that this form of heterogeneity is prevalent and a conserved mechanism by which bacteria cope with fluctuating conditions. We end with an outline of current challenges and future directions for the field. While it remains challenging to predict how an individual bacterium will respond to a given environmental input, we anticipate that capturing the dynamics of the process will begin to resolve this and facilitate rational perturbation of environmental responses for therapeutic and bioengineering purposes.

Mechanisms driving spatial distribution of residents in colony biofilms: an interdisciplinary perspective

Biofilms are consortia of microorganisms that form collectives through the excretion of extracellular matrix compounds. The importance of biofilms in biological, industrial and medical settings has long been recognized due to their emergent properties and impact on surrounding environments. In laboratory situations, one commonly used approach to study biofilm formation mechanisms is the colony biofilm assay, in which cell communities grow on solid-gas interfaces on agar plates after the deposition of a population of founder cells. The residents of a colony biofilm can self-organize to form intricate spatial distributions. The assay is ideally suited to coupling with mathematical modelling due to the ability to extract a wide range of metrics. In this review, we highlight how interdisciplinary approaches have provided deep insights into mechanisms causing the emergence of these spatial distributions from well-mixed inocula.

To biofilm or not to biofilm

A new model helps to predict under which conditions a species of bacteria will switch to a static lifestyle.

A Noninvasive Method for Time-Lapse Imaging of Microbial Interactions and Colony Dynamics

Complex interactions between microbial populations can greatly affect the overall properties of a microbial community, sometimes leading to cooperation and mutually beneficial coexistence, or competition and the death or displacement of organisms or subpopulations. Interactions between different biofilm populations are highly relevant in diverse scientific areas, from antimicrobial resistance to microbial ecology. The utilization of modern microscopic techniques has provided a new and interesting insight into how bacteria interact at the cellular level to form and maintain microbial biofilms. However, our ability to follow complex intraspecies and interspecies interactions in vivo at the microscopic level has remained somewhat limited. Here, we detailed BacLive, a novel noninvasive method for tracking bacterial growth and biofilm dynamics using high-resolution fluorescence microscopy and an associated ImageJ processing macro (https://github.com/BacLive) for easier data handling and image analysis. Finally, we provided examples of how BacLive can be used in the analysis of complex bacterial communities.

IMPORTANCE Communication and interactions between single cells are continuously defining the structure and composition of microbial communities temporally and spatially. Methods routinely used to study these communities at the cellular level rely on sample manipulation which makes microscopic time-lapse experiments impossible. BacLive was conceived as a method for the noninvasive study of the formation and development of bacterial communities, such as biofilms, and the formation dynamics of specialized subpopulations in time-lapse experiments at a colony level. In addition, we developed a tool to simplify the processing and analysis of the data generated by this method.

Nutrient Gradients Mediate Complex Colony-Level Antibiotic Responses in Structured Microbial Populations

Antibiotic treatments often fail to eliminate bacterial populations due to heterogeneity in how individual cells respond to the drug. In structured bacterial populations such as biofilms, bacterial metabolism and environmental transport processes lead to an emergent phenotypic structure and self-generated nutrient gradients toward the interior of the colony, which can affect cell growth, gene expression and susceptibility to the drug. Even in single cells, survival depends on a dynamic interplay between the drug's action and the expression of resistance genes. How expression of resistance is coordinated across populations in the presence of such spatiotemporal environmental coupling remains elusive. Using a custom microfluidic device, we observe the response of spatially extended microcolonies of tetracycline-resistant E. coli to precisely defined dynamic drug regimens. We find an intricate interplay between drug-induced changes in cell growth and growth-dependent expression of resistance genes, resulting in the redistribution of metabolites and the reorganization of growth patterns. This dynamic environmental feedback affects the regulation of drug resistance differently across the colony, generating dynamic phenotypic structures that maintain colony growth during exposure to high drug concentrations and increase population-level resistance to subsequent exposures. A mathematical model linking metabolism and the regulation of gene expression is able to capture the main features of spatiotemporal colony dynamics. Uncovering the fundamental principles that govern collective mechanisms of antibiotic resistance in spatially extended populations will allow the design of optimal drug regimens to counteract them.

A segmentation clock patterns cellular differentiation in a bacterial biofilm

Contrary to multicellular organisms that display segmentation during development, communities of unicellular organisms are believed to be devoid of such sophisticated patterning. Unexpectedly, we find that the gene expression underlying the nitrogen stress response of a developing Bacillus subtilis biofilm becomes organized into a ring-like pattern. Mathematical modeling and genetic probing of the underlying circuit indicate that this patterning is generated by a clock and wavefront mechanism, similar to that driving vertebrate somitogenesis. We experimentally validated this hypothesis by showing that predicted nutrient conditions can even lead to multiple concentric rings, resembling segments. We additionally confirmed that this patterning mechanism is driven by cell-autonomous oscillations. Importantly, we show that the clock and wavefront process also spatially patterns sporulation within the biofilm. Together, these findings reveal a biofilm segmentation clock that organizes cellular differentiation in space and time, thereby challenging the paradigm that such patterning mechanisms are exclusive to plant and animal development.

Emergent Biological Endurance Depends on Extracellular Matrix Composition of Three-Dimensionally Printed Escherichia coli Biofilms

Biofilms are three-dimensional (3D) bacterial communities that exhibit a highly self-organized nature in terms of their composition and complex architecture. Bacteria in biofilms display emergent biological properties, such as resistance to antimicrobials and disinfectants that the individual planktonic cells lack. Bacterial biofilms possess specialized architectural features including unique extracellular matrix compositions and a distinct spatially patterned arrangement of cells and matrix components within the biofilm. It is unclear which of these architectural elements of bacterial biofilms lead to the development of their emergent biological properties. Here, we report a 3D printing-based technique for studying the emergent resistance behaviors of Escherichia coli biofilms as a function of their architecture. Cellulose and curli are the major extracellular-matrix components in E. coli biofilms. We show that 3D-printed biofilms expressing either curli alone or both curli and cellulose in their extracellular matrices show higher resistance to exposure against disinfectants than 3D prints expressing either cellulose alone or no biofilm-matrix components. The 3D-printed biofilms expressing cellulose and/or curli also show thicker anaerobic zones than nonbiofilm-forming E. coli 3D prints. Thus, the matrix composition plays a crucial role in the emergent spatial patterning and biological endurance of 3D-printed biofilms. In contrast, initial spatial distribution of bacterial density or curli-producing cells does not have an effect on biofilm resistance phenotypes. Further, these 3D-printed biofilms could be reversibly attached to different surfaces (bacterial cellulose, glass, and polystyrene) and display resistance to physical distortions by retaining their shape and structure. This physical robustness highlights their potential in applications including bioremediation, protective coatings against pathogens on medical devices, or wastewater treatment, among many others. This new understanding of the emergent behavior of bacterial biofilms could aid in the development of novel engineered living materials using synthetic biology and materials science approaches.

Bacillus subtilis biofilm formation and social interactions

Biofilm formation is a process in which microbial cells aggregate to form collectives that are embedded in a self-produced extracellular matrix. Bacillus subtilis is a Gram-positive bacterium that is used to dissect the mechanisms controlling matrix production and the subsequent transition from a motile planktonic cell state to a sessile biofilm state. The collective nature of life in a biofilm allows emergent properties to manifest, and B. subtilis biofilms are linked with novel industrial uses as well as probiotic and biocontrol processes. In this Review, we outline the molecular details of the biofilm matrix and the regulatory pathways and external factors that control its production. We explore the beneficial outcomes associated with biofilms. Finally, we highlight major advances in our understanding of concepts of microbial evolution and community behaviour that have resulted from studies of the innate heterogeneity of biofilms.

Single-Cell RNA Sequencing Elucidates the Structure and Organization of Microbial Communities

Clonal bacterial populations exhibit various forms of heterogeneity, including co-occurrence of cells with different morphological traits, biochemical properties, and gene expression profiles. This heterogeneity is prevalent in a variety of environments. For example, the productivity of large-scale industrial fermentations and virulence of infectious diseases are shaped by cell population heterogeneity and have a direct impact on human life. Due to the need and importance to better understand this heterogeneity, multiple methods of examining single-cell heterogeneity have been developed. Traditionally, fluorescent reporters or probes are used to examine a specific gene of interest, providing a useful but inherently biased approach. In contrast, single-cell RNA sequencing (scRNA-seq) is an agnostic approach to examine heterogeneity and has been successfully applied to eukaryotic cells. Unfortunately, current extensively utilized methods of eukaryotic scRNA-seq present difficulties when applied to bacteria. Specifically, bacteria have a cell wall which makes eukaryotic lysis methods incompatible, bacterial mRNA has a shorter half-life and lower copy numbers, and isolating an individual bacterial species from a mixed community is difficult. Recent work has demonstrated that these technical hurdles can be overcome, providing valuable insight into factors influencing microbial heterogeneity. This perspective describes the emerging microbial scRNA-seq toolkit. We outline the benefit of these new tools in elucidating numerous scientific questions in microbiological studies and offer insight about the possible rules that govern the segregation of traits in individual microbial cells.

Swarming bacteria undergo localized dynamic phase transition to form stress-induced biofilms

Self-organized multicellular behaviors enable cells to adapt and tolerate stressors to a greater degree than isolated cells. However, whether and how cellular communities alter their collective behaviors adaptively upon exposure to stress is largely unclear. Here, we investigate this question using Bacillus subtilis, a model system for bacterial multicellularity. We discover that, upon exposure to a spatial gradient of kanamycin, swarming bacteria activate matrix genes and transit to biofilms. The initial stage of this transition is underpinned by a stress-induced multilayer formation, emerging from a biophysical mechanism reminiscent of motility-induced phase separation (MIPS). The physical nature of the process suggests that stressors which suppress the expansion of swarms would induce biofilm formation. Indeed, a simple physical barrier also induces a swarm-to-biofilm transition. Based on the gained insight, we propose a strategy of antibiotic treatment to inhibit the transition from swarms to biofilms by targeting the localized phase transition.

Community effects allow bioelectrical reprogramming of cell membrane potentials in multicellular aggregates: Model simulations

Bioelectrical patterns are established by spatiotemporal correlations of cell membrane potentials at the multicellular level, being crucial to development, regeneration, and tumorigenesis. We have conducted multicellular simulations on bioelectrical community effects and intercellular coupling in multicellular aggregates. The simulations aim at establishing under which conditions a local heterogeneity consisting of a small patch of cells can be stabilized against a large aggregate of surrounding identical cells which are in a different bioelectrical state. In this way, instructive bioelectrical information can be persistently encoded in spatiotemporal patterns of separated domains with different cell polarization states. The multicellular community effects obtained are regulated both at the single-cell and intercellular levels, and emerge from a delicate balance between the degrees of intercellular coupling in: (i) the small patch, (ii) the surrounding bulk, and (iii) the interface that separates these two regions. The model is experimentally motivated and consists of two generic voltage-gated ion channels that attempt to establish the depolarized and polarized cell states together with coupling conductances whose individual and intercellular different states permit a dynamic multicellular connectivity. The simulations suggest that community effects may allow the reprogramming of single-cell bioelectrical states, in agreement with recent experimental data. A better understanding of the resulting electrical regionalization can assist the electroceutical correction of abnormally depolarized regions initiated in the bulk of normal tissues as well as suggest new biophysical mechanisms for the establishment of target patterns in multicellular engineering.

The Stress-Responsive Alternative Sigma Factor SigB of Bacillus subtilis and Its Relatives: An Old Friend With New Functions

Alternative sigma factors have led the core RNA polymerase (RNAP) to recognize different sets of promoters to those recognized by the housekeeping sigma A-directed RNAP. This change in RNAP promoter selectivity allows a rapid and flexible reformulation of the genetic program to face environmental and metabolic stimuli that could compromise bacterial fitness. The model bacterium Bacillus subtilis constitutes a matchless living system in the study of the role of alternative sigma factors in gene regulation and physiology. SigB from B. subtilis was the first alternative sigma factor described in bacteria. Studies of SigB during the last 40 years have shown that it controls a genetic universe of more than 150 genes playing crucial roles in stress response, adaption, and survival. Activation of SigB relies on three separate pathways that specifically respond to energy, environmental, and low temperature stresses. SigB homologs, present in other Gram-positive bacteria, also play important roles in virulence against mammals. Interestingly, during recent years, other unexpected B. subtilis responses were found to be controlled by SigB. In particular, SigB controls the efficiencies of spore and biofilm formation, two important features that play critical roles in adaptation and survival in planktonic and sessile B. subtilis communities. In B. subtilis, SigB induces the expression of the Spo0E aspartyl-phosphatase, which is responsible for the blockage of sporulation initiation. The upregulated activity of Spo0E connects the two predominant adaptive pathways (i.e., sporulation and stress response) present in B. subtilis. In addition, the RsbP serine-phosphatase, belonging to the energy stress arm of the SigB regulatory cascade, controls the expression of the key transcription factor SinR to decide whether cells residing in the biofilm remain in and maintain biofilm growth or scape to colonize new niches through biofilm dispersal. SigB also intervenes in the recognition of and response to surrounding microorganisms, a new SigB role that could have an agronomic impact. SigB is induced when B. subtilis is confronted with phytopathogenic fungi (e.g., Fusarium verticillioides) and halts fungal growth to the benefit of plant growth. In this article, we update and review literature on the different regulatory networks that control the activation of SigB and the new roles that have been described the recent years.

The σB-Mediated General Stress Response of Listeria monocytogenes: Life and Death Decision Making in a Pathogen

Sensing and responding to environmental cues is critical for the adaptability and success of the food-borne bacterial pathogen Listeria monocytogenes. A supramolecular multi-protein complex known as the stressosome, which acts as a stress sensing hub, is responsible for orchestrating the activation of a signal transduction pathway resulting in the activation of σ^B, the sigma factor that controls the general stress response (GSR). When σ^B is released from the anti-sigma factor RsbW, a rapid up-regulation of the large σ^B regulon, comprised of ≥ 300 genes, ensures that cells respond appropriately to the new environmental conditions. A diversity of stresses including low pH, high osmolarity, and blue light are known to be sensed by the stressosome, resulting in a generalized increase in stress resistance. Appropriate activation of the stressosome and deployment of σ^B are critical to fitness as there is a trade-off between growth and stress protection when the GSR is deployed. We review the recent developments in this field and describe an up-to-date model of how this sensory organelle might integrate environmental signals to produce an appropriate activation of the GSR. Some of the outstanding questions and challenges in this fascinating field are also discussed.