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.

Magnesium Flux Modulates Ribosomes to Increase Bacterial Survival.

Bacteria exhibit cell-to-cell variability in their resilience to stress, for example, following antibiotic exposure. Higher resilience is typically ascribed to "dormant" non-growing cellular states. Here, by measuring membrane potential dynamics of Bacillus subtilis cells, we show that actively growing bacteria can cope with ribosome-targeting antibiotics through an alternative mechanism based on ion flux modulation. Specifically, we observed two types of cellular behavior: growth-defective cells exhibited a mathematically predicted transient increase in membrane potential (hyperpolarization), followed by cell death, whereas growing cells lacked hyperpolarization events and showed elevated survival. Using structural perturbations of the ribosome and proteomic analysis, we uncovered that stress resilience arises from magnesium influx, which prevents hyperpolarization. Thus, ion flux modulation provides a distinct mechanism to cope with ribosomal stress. These results suggest new approaches to increase the effectiveness of ribosome-targeting antibiotics and reveal an intriguing connection between ribosomes and the membrane potential, two fundamental properties of cells.

SnapShot: Electrochemical Communication in Biofilms.

The role of electricity in biological systems was first appreciated through electrical stimulation experiments performed by Luigi Galvani in the 18th century. These pioneering experiments demonstrated that the behavior of living tissues is governed by the flow of electrochemical species-an insight that gave rise to the modern field of electrophysiology. Since then, electrophysiology has largely remained a bastion of neuroscience. However, exciting recent developments have demonstrated that even simple bacteria residing in communities use electrochemical communication to coordinate population-level behaviors. These recent works are defining the emerging field of bacterial biofilm electrophysiology. To view this SnapShot, open or download the PDF.

Species-Independent Attraction to Biofilms through Electrical Signaling.

Bacteria residing within biofilm communities can coordinate their behavior through cell-to-cell signaling. However, it remains unclear if these signals can also influence the behavior of distant cells that are not part of the community. Using a microfluidic approach, we find that potassium ion channel-mediated electrical signaling generated by a Bacillus subtilis biofilm can attract distant cells. Integration of experiments and mathematical modeling indicates that extracellular potassium emitted from the biofilm alters the membrane potential of distant cells, thereby directing their motility. This electrically mediated attraction appears to be a generic mechanism that enables cross-species interactions, as Pseudomonas aeruginosa cells also become attracted to the electrical signal released by the B. subtilis biofilm. Cells within a biofilm community can thus not only coordinate their own behavior but also influence the behavior of diverse bacteria at a distance through long-range electrical signaling. PAPERCLIP.

Chromosomal Arrangement of Phosphorelay Genes Couples Sporulation and DNA Replication.

Genes encoding proteins in a common regulatory network are frequently located close to one another on the chromosome to facilitate co-regulation or couple gene expression to growth rate. Contrasting with these observations, here, we demonstrate a functional role for the arrangement of Bacillus subtilis sporulation network genes on opposite sides of the chromosome. We show that the arrangement of two sporulation network genes, one located close to the origin and the other close to the terminus, leads to a transient gene dosage imbalance during chromosome replication. This imbalance is detected by the sporulation network to produce cell-cycle coordinated pulses of the sporulation master regulator Spo0A∼P. This pulsed response allows cells to decide between sporulation and continued vegetative growth during each cell cycle spent in starvation. The simplicity of this coordination mechanism suggests that it may be widely applicable in a variety of gene regulatory and stress-response settings. VIDEO ABSTRACT.

Identification of F-actin as the dynamic hub in a microbial-induced GTPase polarity circuit.

Polarity in mammalian cells emerges from the assembly of signaling molecules into extensive biochemical interaction networks. Despite their complexity, bacterial pathogens have evolved parsimonious mechanisms to hijack these systems. Here, we develop a tractable experimental and theoretical model to uncover fundamental operating principles, in both mammalian cell polarity and bacterial pathogenesis. Using synthetic derivatives of the enteropathogenic Escherichia coli guanine-nucleotide exchange factor (GEF) Map, we discover that Cdc42 GTPase signal transduction is controlled by the interaction between Map and F-actin. Mathematical modeling reveals how actin dynamics coupled to a Map-dependent positive feedback loop spontaneously polarizes Cdc42 on the plasma membrane. By rewiring the pathogenic signaling circuit to operate through ß-integrin stimulation, we further show how Cdc42 is polarized in response to an extracellular spatial cue. Thus, a molecular pathway of polarity is proposed, centered on the interaction between GEFs and F-actin, which is likely to function in diverse biological systems.

Architecture-dependent noise discriminates functionally analogous differentiation circuits.

Gene regulatory circuits with different architectures (patterns of regulatory interactions) can generate similar dynamics. This raises the question of why a particular circuit architecture is selected to implement a given cellular process. To investigate this problem, we compared the Bacillus subtilis circuit that regulates differentiation into the competence state to an engineered circuit with an alternative architecture (SynEx) in silico and in vivo. Time-lapse microscopy measurements showed that SynEx cells generated competence dynamics similar to native cells and reconstituted the physiology of differentiation. However, architectural differences between the circuits altered the dynamic distribution of stochastic fluctuations (noise) during circuit operation. This distinction in noise causes functional differences between the circuits by selectively controlling the timing of competence episodes and response of the system to various DNA concentrations. These results reveal a tradeoff between temporal precision and physiological response range that is controlled by distinct noise characteristics of alternative circuit architectures.

Ion channels enable electrical communication in bacterial communities.

The study of bacterial ion channels has provided fundamental insights into the structural basis of neuronal signalling; however, the native role of ion channels in bacteria has remained elusive. Here we show that ion channels conduct long-range electrical signals within bacterial biofilm communities through spatially propagating waves of potassium. These waves result from a positive feedback loop, in which a metabolic trigger induces release of intracellular potassium, which in turn depolarizes neighbouring cells. Propagating through the biofilm, this wave of depolarization coordinates metabolic states among cells in the interior and periphery of the biofilm. Deletion of the potassium channel abolishes this response. As predicted by a mathematical model, we further show that spatial propagation can be hindered by specific genetic perturbations to potassium channel gating. Together, these results demonstrate a function for ion channels in bacterial biofilms, and provide a prokaryotic paradigm for active, long-range electrical signalling in cellular communities.

Metabolic co-dependence gives rise to collective oscillations within biofilms.

Cells that reside within a community can cooperate and also compete with each other for resources. It remains unclear how these opposing interactions are resolved at the population level. Here we investigate such an internal conflict within a microbial (Bacillus subtilis) biofilm community: cells in the biofilm periphery not only protect interior cells from external attack but also starve them through nutrient consumption. We discover that this conflict between protection and starvation is resolved through emergence of long-range metabolic co-dependence between peripheral and interior cells. As a result, biofilm growth halts periodically, increasing nutrient availability for the sheltered interior cells. We show that this collective oscillation in biofilm growth benefits the community in the event of a chemical attack. These findings indicate that oscillations support population-level conflict resolution by coordinating competing metabolic demands in space and time, suggesting new strategies to control biofilm growth.

Bistable emergence of oscillations in growing Bacillus subtilis biofilms.

Biofilm communities of Bacillus subtilis bacteria have recently been shown to exhibit collective growth-rate oscillations mediated by electrochemical signaling to cope with nutrient starvation. These oscillations emerge once the colony reaches a large enough number of cells. However, it remains unclear whether the amplitude of the oscillations, and thus their effectiveness, builds up over time gradually or if they can emerge instantly with a nonzero amplitude. Here we address this question by combining microfluidics-based time-lapse microscopy experiments with a minimal theoretical description of the system in the form of a delay-differential equation model. Analytical and numerical methods reveal that oscillations arise through a subcritical Hopf bifurcation, which enables instant high-amplitude oscillations. Consequently, the model predicts a bistable regime where an oscillating and a nonoscillating attractor coexist in phase space. We experimentally validate this prediction by showing that oscillations can be triggered by perturbing the media conditions, provided the biofilm size lies within an appropriate range. The model also predicts that the minimum size at which oscillations start decreases with stress, a fact that we also verify experimentally. Taken together, our results show that collective oscillations in cell populations can emerge suddenly with nonzero amplitude via a discontinuous transition.

Spiral Wave Propagation in Communities with Spatially Correlated Heterogeneity.

Many multicellular communities propagate signals in a directed manner via excitable waves. Cell-to-cell heterogeneity is a ubiquitous feature of multicellular communities, but the effects of heterogeneity on wave propagation are still unclear. Here, we use a minimal FitzHugh-Nagumo-type model to investigate excitable wave propagation in a two-dimensional heterogeneous community. The model shows three dynamic regimes in which waves either propagate directionally, die out, or spiral indefinitely, and we characterize how these regimes depend on the heterogeneity parameters. We find that in some parameter regimes, spatial correlations in the heterogeneity enhance directional propagation and suppress spiraling. However, in other regimes, spatial correlations promote spiraling, a surprising feature that we explain by demonstrating that these spirals form by a second, distinct mechanism. Finally, we characterize the dynamics using techniques from percolation theory. Despite the fact that percolation theory does not completely describe the dynamics quantitatively because it neglects the details of the excitable propagation, we find that it accounts for the transitions between the dynamic regimes and the general dependency of the spiral period on the heterogeneity and thus provides important insights. Our results reveal that the spatial structure of cell-to-cell heterogeneity can have important consequences for signal propagation in cellular communities.

Statistics of correlated percolation in a bacterial community.

Signal propagation over long distances is a ubiquitous feature of multicellular communities, but cell-to-cell variability can cause propagation to be highly heterogeneous. Simple models of signal propagation in heterogenous media, such as percolation theory, can potentially provide a quantitative understanding of these processes, but it is unclear whether these simple models properly capture the complexities of multicellular systems. We recently discovered that in biofilms of the bacterium Bacillus subtilis, the propagation of an electrical signal is statistically consistent with percolation theory, and yet it is reasonable to suspect that key features of this system go beyond the simple assumptions of basic percolation theory. Indeed, we find here that the probability for a cell to signal is not independent from other cells as assumed in percolation theory, but instead is correlated with its nearby neighbors. We develop a mechanistic model, in which correlated signaling emerges from cell division, phenotypic inheritance, and cell displacement, that reproduces the experimentally observed correlations. We find that the correlations do not significantly affect the spatial statistics, which we rationalize using a renormalization argument. Moreover, the fraction of signaling cells is not constant in space, as assumed in percolation theory, but instead varies within and across biofilms. We find that this feature lowers the fraction of signaling cells at which one observes the characteristic power-law statistics of cluster sizes, consistent with our experimental results. We validate the model using a mutant biofilm whose signaling probability decays along the propagation direction. Our results reveal key statistical features of a correlated signaling process in a multicellular community. More broadly, our results identify extensions to percolation theory that do or do not alter its predictions and may be more appropriate for biological systems.

Slowdown of growth controls cellular differentiation.

How can changes in growth rate affect the regulatory networks behavior and the outcomes of cellular differentiation? We address this question by focusing on starvation response in sporulating Bacillus subtilis We show that the activity of sporulation master regulator Spo0A increases with decreasing cellular growth rate. Using a mathematical model of the phosphorelay-the network controlling Spo0A-we predict that this increase in Spo0A activity can be explained by the phosphorelay protein accumulation and lengthening of the period between chromosomal replication events caused by growth slowdown. As a result, only cells growing slower than a certain rate reach threshold Spo0A activity necessary for sporulation. This growth threshold model accurately predicts cell fates and explains the distribution of sporulation deferral times. We confirm our predictions experimentally and show that the concentration rather than activity of phosphorelay proteins is affected by the growth slowdown. We conclude that sensing the growth rates enables cells to indirectly detect starvation without the need for evaluating specific stress signals.

Noise Expands the Response Range of the Bacillus subtilis Competence Circuit.

Gene regulatory circuits must contend with intrinsic noise that arises due to finite numbers of proteins. While some circuits act to reduce this noise, others appear to exploit it. A striking example is the competence circuit in Bacillus subtilis, which exhibits much larger noise in the duration of its competence events than a synthetically constructed analog that performs the same function. Here, using stochastic modeling and fluorescence microscopy, we show that this larger noise allows cells to exit terminal phenotypic states, which expands the range of stress levels to which cells are responsive and leads to phenotypic heterogeneity at the population level. This is an important example of how noise confers a functional benefit in a genetic decision-making circuit.

Inverse Gillespie for inferring stochastic reaction mechanisms from intermittent samples.

Gillespie stochastic simulation is used extensively to investigate stochastic phenomena in many fields, ranging from chemistry to biology to ecology. The inverse problem, however, has remained largely unsolved: How to reconstruct the underlying reactions de novo from sparse observations. A key challenge is that often only aggregate concentrations, proportional to the population numbers, are observable intermittently. We discovered that under specific assumptions, the set of relative population updates in phase space forms a convex polytope whose vertices are indicative of the dominant underlying reactions. We demonstrate the validity of this simple principle by reconstructing stochastic models (reaction structure plus propensities) from a variety of simulated and experimental systems, where hundreds and even thousands of reactions may be occurring in between observations. In some cases, the inferred models provide mechanistic insight. This principle can lead to the understanding of a broad range of phenomena, from molecular biology to population ecology.

Circuit-level input integration in bacterial gene regulation.

Gene regulatory circuits can receive multiple simultaneous inputs, which can enter the system through different locations. It is thus necessary to establish how these genetic circuits integrate multiple inputs as a function of their relative entry points. Here, we use the dynamic circuit regulating competence for DNA uptake in Bacillus subtilis as a model system to investigate this issue. Specifically, we map the response of single cells in vivo to a combination of (i) a chemical signal controlling the constitutive expression of key competence genes, and (ii) a genetic perturbation in the form of copy number variation of one of these genes, which mimics the level of stress signals sensed by the bacteria. Quantitative time-lapse fluorescence microscopy shows that a variety of dynamical behaviors can be reached by the combination of the two inputs. Additionally, the integration depends strongly on the relative locations where the two perturbations enter the circuit. Specifically, when the two inputs act upon different circuit elements, their integration generates novel dynamical behavior, whereas inputs affecting the same element do not. An in silico bidimensional bifurcation analysis of a mathematical model of the circuit offers good quantitative agreement with the experimental observations, and sheds light on the dynamical mechanisms leading to the different integrated responses exhibited by the gene regulatory circuit.

Localized cell death focuses mechanical forces during 3D patterning in a biofilm.

From microbial biofilm communities to multicellular organisms, 3D macroscopic structures develop through poorly understood interplay between cellular processes and mechanical forces. Investigating wrinkled biofilms of Bacillus subtilis, we discovered a pattern of localized cell death that spatially focuses mechanical forces, and thereby initiates wrinkle formation. Deletion of genes implicated in biofilm development, together with mathematical modeling, revealed that ECM production underlies the localization of cell death. Simultaneously with cell death, we quantitatively measured mechanical stiffness and movement in WT and mutant biofilms. Results suggest that localized cell death provides an outlet for lateral compressive forces, thereby promoting vertical mechanical buckling, which subsequently leads to wrinkle formation. Guided by these findings, we were able to generate artificial wrinkle patterns within biofilms. Formation of 3D structures facilitated by cell death may underlie self-organization in other developmental systems, and could enable engineering of macroscopic structures from cell populations.

Synthetic biological networks.

Despite their obvious relationship and overlap, the field of physics is blessed with many insightful laws, while such laws are sadly absent in biology. Here we aim to discuss how the rise of a more recent field known as synthetic biology may allow us to more directly test hypotheses regarding the possible design principles of natural biological networks and systems. In particular, this review focuses on synthetic gene regulatory networks engineered to perform specific functions or exhibit particular dynamic behaviors. Advances in synthetic biology may set the stage to uncover the relationship of potential biological principles to those developed in physics.

Biological role of noise encoded in a genetic network motif.

Genetic circuits that regulate distinct cellular processes can differ in their wiring pattern of interactions (architecture) and susceptibility to stochastic fluctuations (noise). Whether the link between circuit architecture and noise is of biological importance remains, however, poorly understood. To investigate this problem, we performed a computational study of gene expression noise for all possible circuit architectures of feed-forward loop (FFL) motifs. Results revealed that FFL architectures fall into two categories depending on whether their ON (stimulated) or OFF (unstimulated) steady states exhibit noise. To explore the biological importance of this difference in noise behavior, we analyzed 858 documented FFLs in Escherichia coli that were divided into 39 functional categories. The majority of FFLs were found to regulate two subsets of functional categories. Interestingly, these two functional categories associated with FFLs of opposite noise behaviors. This opposite noise preference revealed two noise-based strategies to cope with environmental constraints where cellular responses are either initiated or terminated stochastically to allow probabilistic sampling of alternative states. FFLs may thus be selected for their architecture-dependent noise behavior, revealing a biological role for noise that is encoded in gene circuit architectures.

Temporal competition between differentiation programs determines cell fate choice.

Multipotent differentiation, where cells adopt one of several possible fates, occurs in diverse systems ranging from bacteria to mammals. This decision-making process is driven by multiple differentiation programs that operate simultaneously in the cell. How these programs interact to govern cell fate choice is poorly understood. To investigate this issue, we simultaneously measured activities of the competing sporulation and competence programs in single Bacillus subtilis cells. This approach revealed that these competing differentiation programs progress independently without cross-regulation before the decision point. Cells seem to arrive at a fate choice through differences in the relative timing between the two programs. To test this proposed dynamic mechanism, we altered the relative timing by engineering artificial cross-regulation between the sporulation and competence circuits. Results suggest a simple model that does not require a checkpoint or intricate cross-regulation before cellular decision-making. Rather, cell fate choice appears to be the outcome of a 'molecular race' between differentiation programs that compete in time, providing a simple dynamic mechanism for decision-making.