Reversible and Noisy Progression towards a Commitment Point Enables Adaptable and Reliable Cellular Decision-Making: e1002273

Abstract

Cells must make reliable decisions under fluctuating extracellular conditions, but also be flexible enough to adapt to such changes. How cells reconcile these seemingly contradictory requirements through the dynamics of cellular decision-making is poorly understood. To study this issue we quantitatively measured gene expression and protein localization in single cells of the model organism Bacillus subtilis during the progression to spore formation. We found that sporulation proceeded through noisy and reversible steps towards an irreversible, all-or-none commitment point. Specifically, we observed cell-autonomous and spontaneous bursts of gene expression and transient protein localization events during sporulation. Based on these measurements we developed mathematical population models to investigate how the degree of reversibility affects cellular decision-making. In particular, we evaluated the effect of reversibility on the 1) reliability in the progression to sporulation, and 2) adaptability under changing extracellular stress conditions. Results show that reversible progression allows cells to remain responsive to long-term environmental fluctuations. In contrast, the irreversible commitment point supports reliable execution of cell fate choice that is robust against short-term reductions in stress. This combination of opposite dynamic behaviors (reversible and irreversible) thus maximizes both adaptable and reliable decision-making over a broad range of changes in environmental conditions. These results suggest that decision-making systems might employ a general hybrid strategy to cope with unpredictably fluctuating environmental conditions.

Abstract

Cells must make reliable decisions under fluctuating extracellular conditions, but also be flexible enough to adapt to such changes. How cells reconcile these seemingly contradictory requirements through the dynamics of cellular decision-making is poorly understood. To study this issue we quantitatively measured gene expression and protein localization in single cells of the model organism Bacillus subtilis during the progression to spore formation. We found that sporulation proceeded through noisy and reversible steps towards an irreversible, all-or-none commitment point. Specifically, we observed cell-autonomous and spontaneous bursts of gene expression and transient protein localization events during sporulation. Based on these measurements we developed mathematical population models to investigate how the degree of reversibility affects cellular decision-making. In particular, we evaluated the effect of reversibility on the 1) reliability in the progression to sporulation, and 2) adaptability under changing extracellular stress conditions. Results show that reversible progression allows cells to remain responsive to long-term environmental fluctuations. In contrast, the irreversible commitment point supports reliable execution of cell fate choice that is robust against short-term reductions in stress. This combination of opposite dynamic behaviors (reversible and irreversible) thus maximizes both adaptable and reliable decision-making over a broad range of changes in environmental conditions. These results suggest that decision-making systems might employ a general hybrid strategy to cope with unpredictably fluctuating environmental conditions.

Author Summary

Cells must continuously make decisions in response to changes in their environment. These decisions must be irreversible, to prevent cells from reverting back to unfit cellular states, but also be flexible, to allow cells to go back to their previous state upon environmental changes. Using single-cell time-lapse fluorescence microscopy, we show that these seemingly contradictory properties coexist in Bacillus subtilis cells during their progression to spore formation. We suggest, on the basis of a mathematical population model, that reversible progression towards the irreversible decision to sporulate optimizes respectively adaptability and reliability of decision-making over a broad range of changes in environmental conditions.

Citation: Kuchina A, Espinar L, Garcia-Ojalvo J, Süel GM (2011) Reversible and Noisy Progression towards a Commitment Point Enables Adaptable and Reliable Cellular Decision-Making. PLoS Comput Biol 7(11): e1002273. doi:10.1371/journal.pcbi.1002273

Editor: Markus W. Covert, Stanford University, United States of America

Received: July 23, 2011; Accepted: September 29, 2011; Published: November 10, 2011

Copyright: © 2011 Kuchina et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Funding: GMS acknowledges support by NIH Grant No. NIGMS RO1 GM088428, Welch Foundation (Grant No. I-1674), and James S. McDonnell Foundation (Grant No. 220020141). JGO acknowledges financial support from the Ministerio de Ciencia e Innovación (Spain, Project No. FIS2009-13360) and the ICREA Academia programme. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Cellular decision-making underlies many biological processes such as multipotent differentiation, where cells commit to one of several distinct fates. Such cell fate choice must permit individual cells to reach a decision even in fluctuating extracellular environments [1]. At the same time, cells must also be able to adapt their cell fate choice to changes in these conditions. It is unclear how individual cells reconcile these opposing requirements of decisiveness and adaptability during decision-making. Decisive cellular differentiation mechanisms have been proposed to combine ultra-sensitivity and positive feedback to generate an irreversible and all-or-none cell fate choice such as those observed during Xenopus oocyte maturation [2] and yeast mating decision [3]. However, individual cells with irreversible responses can lack the flexibility to respond proportionally to changing environments, since even small changes can trigger irreversible responses. In contrast, progression of cellular differentiation through reversible intermediate states permits flexibility and proportional responses to environmental changes. For example, multipotent differentiation of hematopoietic stem cells is a stepwise process with numerous reversible intermediate states that allows cells to gradually adapt to changes in extracellular signals [4], [5], [6], [7], [8], [9], [10]. Despite these recent insights, how multipotent differentiation systems reach a decisive cell fate choice while maintaining the ability to respond to changes in the environment is largely unknown.

To understand cellular decision-making it is critical to determine the single-cell dynamics underlying the progression to cell fate choice. However, these dynamics are poorly characterized in multipotent differentiation systems ranging from bacteria to mammalian stem cells. Simultaneous measurement of multiple components of a differentiation program in the same cell can reveal the dynamics of cellular decision-making underlying multipotent differentiation. The soil bacterium Bacillus subtilis serves as an ideal model system in which the dynamics of multiple genes within a differentiation circuit are simultaneously measurable in single cells [11], [12], [13]. In stressful environments the majority of B. subtilis cells form spores that survive environmental extremes [14], [15]. The sporulation program has been well characterized genetically and multiple stages of sporulation have been described [16], [17], [18], [19], [20]. However, despite these important insights, how individual cells proceed to spore formation and thus the dynamics of the sporulation program in single cells has not been determined.

To uncover cell fate choice dynamics in B. subtilis, we simultaneously measured the activities of multiple sporulation circuit elements with single-cell resolution. We found that individual cells progress gradually to spore formation through reversible activities of early sporulation components. These measurements also allowed us to confirm the previously identified commitment point of cellular decision-making at which cells irreversibly proceed to complete spore formation [21]. More importantly, single cells analysis revealed the precise all-or-none dynamics of this decision-making point that was obscured by variability at the single-cell level in population measurements. Modeling of alternative sporulation dynamics showed that the combination of reversible and irreversible dynamics employed by B. subtilis can represent a general strategy to maximize reliable and yet adaptable cellular decision-making over a broad range of randomly fluctuating environmental conditions.

Results

Single cells of Bacillus subtilis progress reversibly and gradually to spore formation

First, we established the single-cell dynamics of the progression toward spore formation by measuring the temporal profile of four sporulation components (Fig. 1A). Initiation of sporulation is controlled by a multicomponent phosphorelay including two phosphotransferases, Spo0F and Spo0B, and a transcription factor Spo0A [14], [19], [22], [23], [24], [25]. Spo0F protein senses and integrates a variety of inputs and as a result becomes phosphorylated [23], [24]. The phosphate group is subsequently transferred through Spo0B to Spo0A, a master regulator of sporulation, which upon phosphorylation directly controls expression of over 120 genes [20], [26], [27], including spo0F (in what constitutes a feedback mechanism [28], [29], [30]). Among these is the gene for SpoIIE, a protein phosphatase that is required for the activation of the forespore specific transcription factor (σF) and whose localization to the asymmetric septum is the first morphological marker for forespore formation [31], [32], [33], [34], [35]. Activation of σF in the forespore switches on the expression of SpoIIR, which in turn leads to the activation of σE, a transcription factor specific to the mother cell that switches on the expression of a large number of late sporulation genes [36], [37], [38]. At this point, the cell becomes irreversibly committed to sporulation [21], [37], [38], [39]. Using fluorescent reporter constructs, we quantified the activities of the Spo0A (Pspo0A), Spo0F (Pspo0F) and SpoIIR (PspoIIR) promoters (Fig. 1A, top two rows). We also visualized SpoIIE protein localization by utilizing a functional translational fusion to SpoIIE (Fig. 1A, bottom row). Each sporulation reporter was simultaneously measured in combination with the PspoIIR reporter. Measurement of these overlapping pair-wise combinations of reporters allowed us to establish a relative temporal profile for multiple sporulation steps with single-cell resolution (Fig. 1B and Supporting Fig. S1) that was consistent with the genetically established hierarchy within the sporulation circuit [14], [18], [22], [40], [41], [42] (Fig. 1C).

[Figure omitted, see PDF]

Figure 1. Reversibility of sporulation component activities during the progression to spore formation.

(A) Filmstrips of phase and fluorescence image overlays of sporulating B. subtilis cells expressing pair-wise combinations of fluorescent proteins. For all panels, CFP fluorescence expressed from the PspoIIR promoter is colored in green. In the same cells, YFP is expressed from the Pspo0A promoter (upper panel, blue), from the Pspo0F promoter (middle panel, red), and as a translational fusion to SpoIIE protein (lower panel, orange). Time is indicated in hours. (B) Mean quantitative time traces of sporulation reporters during typical sporulation events. Dynamics of the reporters were obtained from strains expressing pair-wise combinations of indicated sporulation reporters (strains 0A-IIR, n = 33; 0F-IIR, n = 30; IIE-IIR, n = 28; strain definitions and genetic background can be found in Supporting Text S2 and Supporting Table S1) and aligned in time with respect to activation of the common PspoIIR reporter (defined by PspoIIR fluorescence higher than 70% of maximum PspoIIR fluorescence at sporulation). Activities of all reporters were measured as mean fluorescence intensities in single cells. All traces have been normalized for amplitude. Dashed lines above and below the mean curves indicate standard deviation (SD). See also Supporting Fig. S1. (C) Scheme representing the temporal sequence of cellular states during the progression to sporulation, summarizing data from the plot shown in (B).

Single cell measurements of sporulation reporters revealed that progression to spore formation is comprised of reversible steps. During the early progression to spore formation, individual cells exhibited bursts of Pspo0A gene expression that did not result in spore formation (Fig. 2A). Activation of Spo0A during B. subtilis sporulation is known to be heterogeneous among single cells [43], [44]. However, spontaneous bursts of gene expression at single-cell level during sporulation have not been described to date. Bursts in gene expression were not limited to Pspo0A, but were also observed for Pspo0F (Fig. 2B). Given that Pspo0F is, as mentioned above, transcriptionally activated by phosphorylated Spo0A, the dynamics of that promoter is reporting the activity of Spo0A. Consequently, the observed bursts in Pspo0F indicate that not only the expression, but also the activity of the Spo0A master regulator exhibits bursts in single cells. In contrast, bursts of gene expression were not observed for the late stage sporulation reporter PspoIIR (Fig. 2A-C): after the sharp signal increase observed for PspoIIR, a spore is always formed (and thus the time traces shown in the plots cannot be continued). Gene expression bursts can introduce variability and reversibility during the early stages of the progression to spore formation, whereas the later stages do not appear to be subject to such stochasticity. Additionally, we observed reversible protein localization of SpoIIE in approximately 2±1% (SEM) of cells (Fig. 2C). Specifically, we observed transient localization events of the SpoIIE-YFP fusion protein to the asymmetric septum that did not give rise to spore formation. In these cells, SpoIIE either switched its localization between opposite poles, or completely delocalized and cells continued with cell division (Fig. 2C). Similar bursts of promoter activity and protein localization in single cells have also been reported in other systems, and have been attributed to the stochastic and reversible nature of the underlying biochemical reactions [43], [45], [46], [47], [48], [49], [50].

[Figure omitted, see PDF]

Figure 2. Smoothed quantitative time traces of normalized fluorescence signals reporting sporulation component activities in single cells.

(A) 0A(5x)-IIR strain with Pspo0A-CFP (blue) and PspoIIR-YFP (green). (B) 0F-IIR strain with Pspo0F-YFP (red) and PspoIIR-CFP (green). (C) IIE-IIR strain showing localization of fusion protein SpoIIE-YFP (orange) and PspoIIR-CFP (green). For SpoIIE-YFP, local fluorescence measurements near the cell poles are shown. On top of each panel, cell images are shown at sample time points. The measured cell is indicated with a yellow outline.

The single-cell reversibility observed here is distinct from what has been reported in population-level studies which showed that sporulation can be aborted upon transfer from stress to rich media conditions [21]. In contrast, the single-cell reversibility discussed here appears to be a cell autonomous behavior that occurs randomly without requiring a change in environmental conditions. Such cell autonomous and random sampling of heterogeneous behavior in single cells has also been observed in other systems where it may provide a fitness advantage, such as the generation of antibiotic resistant Escherichia coli persister cells [51] and the survival of Saccharomyces cerevisiae to cellular stress [50]. Together, our data reveal that the early steps of progression towards spore formation are subject to stochastic reversibility at the single-cell level, suggesting an adaptable progression to spore formation.

Bacillus subtilis commit to sporulation in an all-or-none fashion

Detailed analysis of single cell dynamics revealed that the irreversible commitment for spore formation is executed within a narrow time window. Specifically, we find that PspoIIR activation (which as mentioned above, denotes the sporulation commitment point) is switch-like (Fig. 2A-C, green line). This timing was quantified by using the morphological appearance of an actual spore (indicated by a phase bright spot in our single-cell movies) as a reference time point from which to measure commitment time (Fig. 3A). Analysis of single-cell data shows that the temporal distance between PspoIIR activation and the formation of the morphologically visible forespore is tightly distributed (CV = 0.2, and 15% of median cell cycle duration) (Fig. 3B). These results show the precise timing of the irreversible decision and the completion of the spore formation process. Therefore, single cell measurements of sporulation dynamics exposed a precise cell intrinsic decision point in time that revealed switch-like dynamics of the commitment to spore formation. The temporal precision of this decision-making point was previously concealed in population measurements (see Fig. 1B) by the cell-to-cell variability in the reversible progression to spore formation.

[Figure omitted, see PDF]

Figure 3. Precise timing of PspoIIR activation prior to spore formation.

(A) Filmstrips showing fluorescence (top) and phase (bottom) images of a sample sporulating cell (strain 0A-IIR) activating PspoIIR expression in the forespore (green). Phase bright spore appearance is indicated by an arrow. (B) Histogram presenting, in green, the distribution of time between PspoIIR activation and phase bright spore formation (n = 31, sample events highlighted by yellow boxes in the filmstrips in panel (A)), compared to sample cell cycle durations, in light grey, measured in the same cells under movie-making conditions. PspoIIR activation is measured as >60% of maximum PspoIIR fluorescence at sporulation.

Mathematical modeling reveals advantage of B. subtilis cellular decision-making dynamics under alternating stress conditions

Our measurements suggest that the decision to sporulate in B. subtilis is made in an irreversible all-or-none manner, following a reversible and gradual progression toward this decision point. Therefore, sporulation combines reversible and irreversible dynamics, two seemingly opposed decision mechanisms. In order to examine the potential advantages of this hybrid mechanism, we compared the response of three simplified models of cellular decision-making to a variable-stress environment. These models describe the dynamics of a population of cells progressing toward sporulation under stress, in terms of the number of cells existing at any given time in a given state along the progression, beginning with the vegetative state and ending in the spore state. The dynamics of cell populations in all these states is given by a set of coupled ordinary differential equations that are linear, and thus can be solved exactly (Supporting Text S1). The three models, shown schematically in Fig. 4A and described in detail in the Supporting Information (Supporting Text S1 and Supporting Fig. S2), involve either 1) a purely irreversible/all-or-none, 2) a purely reversible/gradual, or 3) a "hybrid" process that takes cells from their initial vegetative state to their final spore state. In the irreversible-only scheme, cells decide in a single step whether or not to sporulate, and the decision is irreversible. The reversible-only model takes cells gradually toward sporulation through multiple reversible intermediate states, without any irreversible commitment step taking place along the process. Only the ultimate transition to the spore state, taking place after the decision, is irreversible. Finally, the third "hybrid" model features the actual sporulation dynamics identified here for B. subtilis that combines a gradual progression through reversible intermediate states with an irreversible all-or-none decision prior to the spore transition.

[Figure omitted, see PDF]

Figure 4. Mathematical comparison of distinct cell-fate choice dynamics under alternating environmental conditions.

(A) Schematic representation of the three models. Irreversible-only, top panel: Decision to sporulate is made by cells in a single step. The model comprises three cellular states: vegetative (V), decided (D) and spores (S). Reversible-only, middle panel: In this model vegetative cells (V) proceed toward sporulation through two intermediate states I1 and I2. All transitions, including the transition to the decided state (D), are reversible. Hybrid, lower panel: This model combines characteristics from the two mechanisms described above. In this scenario, cells initially progress toward sporulation through intermediate states I1 and I2 with reversible transitions, with the irreversible final transition to the decided state (D). For all models the transition rates between sporulation stages, as well as growth and death (Ø) rates, are coupled to the stress level. For details, see Supporting Text S1. (B) Relative fitness of the hybrid model with respect to the irreversible-only and reversible-only models described in Fig. 3A. The blue and red lines correspond to the ratio of the total cell populations (summed over all cell states) of the hybrid model and the irreversible-only and reversible-only models, respectively. A random dichotomic variation of the environment is applied to all models, with the relative average duration of high versus low stress phases increasing in the x axis. The two insets depict sample profiles of stress for relative average durations of high stress phase equal to 0.2 and 0.8. The shaded areas in light blue and light red represent the maximum variation resulting from simultaneous random changes in all non-zero parameters of up to 20% around their baseline values given in Supporting Table S2. (C) Average fitness of the hybrid model with respect to the irreversible-only and reversible-only models obtained from a random sampling of environmental conditions. The inset shows short sample series of the average stress levels used to study the response of the three models to random stress conditions. Relative stress durations where chosen from a uniform random distribution in the range [0, 1]. The blue and red bars represent the ratio of the total cell population of the hybrid model to that of irreversible-only and reversible-only, respectively, averaged over 100 values of the relative average duration of high versus low stress scanned in (B). For each stress level, the models were integrated over 30 hours, in the same conditions as in (B). The error bars represent the maximum variation resulting from parameter changes up to 20% their baseline values, as in (B).

We subjected the three models described above to a random variation of the environment, with alternation of a high and a low level of stress (insets in Fig. 4B and Supporting Fig. S3). For all models, the rates of progression toward/back from sporulation, growth and death rates are all coupled to the level of stress (Supporting Text S1). We then systematically varied the ratio of high stress to total cycle duration, which was kept constant. The relative fitness of the hybrid model with respect to the irreversible-only and reversible-only models was measured as the ratio of total cells in the population between the models (Fig. 4B, blue and red lines respectively).

Results of this analysis show that the survival rate of the different models depends on the stress profile. None of the models dominates over the entire range of environmental stress conditions, but the hybrid model performs best overall. For short phases of high stress, the hybrid outgrows the irreversible-only model (left half of main plot in Fig. 4B, blue line), since the latter is driven irreversibly to sporulation even for short stress periods. Under these conditions, the hybrid and reversible-only models perform similarly (red line) since in both models reversible progression delays sporulation. However, in the opposite limit where high stress durations approach the total duration (right half of main plot in Fig. 4B), the irreversible-only model outgrows the hybrid model, since responsiveness is a disadvantage. For such prolonged high stress durations, spores are at an advantage because non-spore cells have a higher death than growth rate under stress. In this limit, the hybrid model is in turn more reliable than the reversible-only model, since the reversible-only model is driven away from sporulation even by short intervals of rich phase. Parameter sensitivity analysis showed that these results hold over a range of parameter values (Fig. 4B, light blue and light red regions). Taken together, these results show that the hybrid model outperforms both the irreversible-only and reversible-only schemes over a broader range of randomly alternating stress profiles, given its responsiveness to the long-term recovery from stress and its reliability during short-term release of stress conditions.

This prediction was tested directly by subjecting all three models to random changes in the fraction of time spent in high stress. In these extended simulations cells were exposed to 100 random samplings of environmental stress ratios (part of an environmental profile is shown in Fig. 4C, inset). We then calculated the ratio of the total cell population of the hybrid model to that of the irreversible-only and reversible-only models averaged over the entire simulation, defining the relative fitness of each model (Fig. 4C). Our data demonstrate that the hybrid model has an overall higher fitness over both the irreversible-only and reversible-only models when populations are subjected to a randomly changing range of environmental conditions. Therefore, the combination of gradual progression toward an all-or-none decision in B. subtilis, as represented in the hybrid model and observed experimentally, enables the system to cope with a broader range of unpredictable stress conditions.

Discussion

Single cell measurements of sporulation components allowed us to precisely determine both the temporal progression of the earliest sporulation events and the consecutive switch-like dynamics of the commitment point. The temporal precision and switch-like nature of the decision point would have been concealed in population measurements that describe average behaviors of cells and cannot account for cell-to-cell variability. Similarly, a recent study of the post-infection decision in bacteriophage lambda to undergo lysis or lysogeny showed how cell-cell variability obscured the precision of cell fate choice [52]. Therefore, stochastic fluctuations and noise observed at the single-cell level can conceal the temporal precision of numerous cellular processes. Our approach based on measuring multiple gene circuit components simultaneously in the same cell can identify cell intrinsic temporal reference points that can accurately establish the relative timing of cellular processes.

Multipotent differentiation in B. subtilis appears to combine the opposing dynamic regimes of reversible/gradual and irreversible/all-or-none behavior to reconcile the seemingly contradictory requirements of adaptable and reliable cellular decision-making. The time required for cells to reach spore formation is highly variable, allowing cells to generate a broad distribution of wait times prior to cell fate choice. This extends the period during which B. subtilis cells are known to be responsive to environmental changes. Therefore, variability in sporulation progression can be beneficial biologically. These findings thus extend previous results regarding the advantageous role of stochasticity in enhancing survival under uncertain environments [13], [50], [53], [54], [55], [56]. In contrast, the actual decision to sporulate is governed by irreversible and all-or-none dynamics, which provides reliable execution of the sporulation program. Our modeling suggests that this combination of gradual and all-or-none dynamics during cell fate choice allows B. subtilis to successfully survive under a broad range of alternating environmental stress profiles. In mammalian cells, the decision for apoptosis has been described as a slow progression towards a fast decision, suggesting that a hybrid strategy similar to that described here might be employed [57]. Therefore, this strategy of combining opposing dynamic behaviors may be common to various other biological processes, and may represent a general mechanism for decision-making under unpredictably changing environments. Models for decision-making under uncertain conditions have been developed and applied to numerous unrelated complex systems, such as finance [58]. Perhaps a stochastically reversible progression to an irreversible all-or-none switch as observed here in cells, may also serve as a strategy for adaptable and reliable decision-making in other complex systems that are subject to unpredictable conditions.

Materials and Methods

Strain construction

Bacillus subtilis strains used in the study are isogenic to wild-type B. subtilis PY79 strain and are listed in Supporting Table S1. Promoter - fluorescent proteins fusions were generated using fusion polymerase chain reaction and cloned into B. subtilis chromosomal integration vectors following standard protocols. The following B. subtilis chromosomal integration vectors were used: pSac-Cm, integrating into the sacA locus (constructed by R. Middleton and obtained from the Bacillus Genetic Stock Center) and pLD30 integrating into the amyE locus (kind gift from Jonathan Dworkin, Columbia University). We have also utilized the bifunctional cloning plasmid pHP13 carrying the replication origin of the cryptic B. Subtilis plasmid pTA1060 (5 copies per genome) [59]. Standard one-step transformation protocol was followed to transform B. subtilis with these constructs.

Growth and imaging conditions

For imaging, B. subtilis cells were grown at 37°C in LB with the following concentrations of appropriate antibiotic: 5 µg/ml chloramphenicol, 5 µg/ml neomycin, 5 µg/ml erythromycin or 100 µg/ml spectinomycin. After reaching OD 1.8, the cells were resuspended in 0.5 volume of Resuspension Media (RM; see Supporting Text S2, Materials and Methods) supplemented with 0.02% glucose. The cells were incubated at 37°C for 1 hour, then diluted 10-fold in RM and applied onto a 1.5% low-melting agarose pad placed into a coverslip-bottom Willco dish for imaging.

Time-lapse microscopy

Growth of microcolonies was observed with fluorescence time-lapse microscopy at 37°C with an Olympus IX-81 inverted microscope with a motorized stage (ASI) and an incubation chamber. Image sets were acquired every 20 min with a Hamamatsu ORCA-ER camera. The imaging time has been optimized in order to prevent phototoxicity [11]. Custom Visual Basic software in combination with the Image Pro Plus (Media Cybernetics) was used to automate image acquisition and microscope control.

Image analysis

Time-lapse movie data analysis was performed by custom software developed with MATLAB image processing and statistics toolboxes (The Mathworks) described in [60] and [11].

Mathematical modeling

Extended description of the models and methods is available in Supporting Text S1, Mathematical Modeling.

Supporting Information

Figure S1.

Single-cell traces of sporulating cells expressing pair-wise combinations of fluorescent sporulation markers. In panels (A) through (C), quantitative time traces of sporulation reporters during typical sporulation events are shown in lighter colors, with mean trace shown on top in bright color and dashed lines indicating standard deviation (SD). Dynamics of the reporters in each panel were obtained from strains expressing pair-wise combinations of indicated sporulation reporters: (A) strain 0A-IIR, n = 33, (B) strain 0F-IIR, n = 30, (C) strain IIE-IIR, n = 28. All traces were aligned with respect to PspoIIR activation (in green) defined as >70% of fluorescence intensity relative to maximum intensity at sporulation. Data from these panels were combined to produce Fig. 1B.

(PDF)

Figure S2.

Three alternative models of sporulation progression.

(PDF)

Figure S3.

Stress modulation.

(PDF)

Table S1.

Strain definitions and genetic background. In the first column, promoters expressing fluorescent proteins are abbreviated as follows: "0A", Pspo0A; "0F", Pspo0F; IIE, spoIIE; "IIR", PspoIIR. (5x) denotes average copy number of plasmid pHP13 (this strain was used where indicated in place of 0A-IIR for increased signal-to-noise ratio).

(PDF)

Table S2.

Parameter values in the rich and poor phases.

(PDF)

Text S1.

Mathematical modeling.

(PDF)

Text S2.

Materials and methods.

(PDF)

Author Contributions

Conceived and designed the experiments: AK LE JGO GMS. Performed the experiments: AK LE. Analyzed the data: AK LE JGO GMS. Contributed reagents/materials/analysis tools: AK LE JGO GMS. Wrote the paper: AK LE JGO GMS.

References

Perkins TJ, Swain PS (2009) Strategies for cellular decision-making. Mol Syst Biol 5: 326. Find this article online

Ferrell JE Jr, Machleder EM (1998) The biochemical basis of an all-or-none cell fate switch in Xenopus oocytes. Science 280: 895-898. Find this article online

Malleshaiah MK, Shahrezaei V, Swain PS, Michnick SW (2010) The scaffold protein Ste5 directly controls a switch-like mating decision in yeast. Nature 465: 101-105. Find this article online

Bertrand JY, Traver D (2009) Hematopoietic cell development in the zebrafish embryo. Curr Opin Hematol 16: 243-248. Find this article online

Sanosaka T, Namihira M, Nakashima K (2009) Epigenetic mechanisms in sequential differentiation of neural stem cells. Epigenetics 4: 89-92. Find this article online

Parra M (2009) Epigenetic events during B lymphocyte development. Epigenetics 4: 462-468. Find this article online

Malhotra S, Kincade PW (2009) Wnt-related molecules and signaling pathway equilibrium in hematopoiesis. Cell Stem Cell 4: 27-36. Find this article online

Quesenberry PJ, Aliotta JM (2008) The paradoxical dynamism of marrow stem cells: considerations of stem cells, niches, and microvesicles. Stem Cell Rev 4: 137-147. Find this article online

Prindull GA, Fibach E (2007) Are postnatal hemangioblasts generated by dedifferentiation from committed hematopoietic stem cells? Exp Hematol 35: 691-701. Find this article online

Moldovan NI (2005) Functional adaptation: the key to plasticity of cardiovascular "stem" cells? Stem Cells Dev 14: 111-121. Find this article online

Suel GM, Garcia-Ojalvo J, Liberman LM, Elowitz MB (2006) An excitable gene regulatory circuit induces transient cellular differentiation. Nature 440: 545-550. Find this article online

Suel GM, Kulkarni RP, Dworkin J, Garcia-Ojalvo J, Elowitz MB (2007) Tunability and noise dependence in differentiation dynamics. Science 315: 1716-1719. Find this article online

Cagatay T, Turcotte M, Elowitz MB, Garcia-Ojalvo J, Suel GM (2009) Architecture-dependent noise discriminates functionally analogous differentiation circuits. Cell 139: 512-522. Find this article online

Piggot PJ, Hilbert DW (2004) Sporulation of Bacillus subtilis. Curr Opin Microbiol 7: 579-586. Find this article online

Errington J (1993) Bacillus subtilis sporulation: regulation of gene expression and control of morphogenesis. Microbiol Rev 57: 1-33. Find this article online

Schultz D, Wolynes PG, Ben Jacob E, Onuchic JN (2009) Deciding fate in adverse times: sporulation and competence in Bacillus subtilis. Proc Natl Acad Sci U S A 106: 21027-21034. Find this article online

Veening JW, Smits WK, Hamoen LW, Kuipers OP (2006) Single cell analysis of gene expression patterns of competence development and initiation of sporulation in Bacillus subtilis grown on chemically defined media. J Appl Microbiol 101: 531-541. Find this article online

Hilbert DW, Piggot PJ (2004) Compartmentalization of gene expression during Bacillus subtilis spore formation. Microbiol Mol Biol Rev 68: 234-262. Find this article online

Fujita M, Losick R (2005) Evidence that entry into sporulation in Bacillus subtilis is governed by a gradual increase in the level and activity of the master regulator Spo0A. Genes Dev 19: 2236-2244. Find this article online

Fujita M, Gonzalez-Pastor JE, Losick R (2005) High- and low-threshold genes in the Spo0A regulon of Bacillus subtilis. J Bacteriol 187: 1357-1368. Find this article online

Dworkin J, Losick R (2005) Developmental commitment in a bacterium. Cell 121: 401-409. Find this article online

Stragier P, Losick R (1996) Molecular genetics of sporulation in Bacillus subtilis. Annu Rev Genet 30: 297-241. Find this article online

Burbulys D, Trach KA, Hoch JA (1991) Initiation of sporulation in B. subtilis is controlled by a multicomponent phosphorelay. Cell 64: 545-552. Find this article online

Hoch JA (1993) The phosphorelay signal transduction pathway in the initiation of Bacillus subtilis sporulation. J Cell Biochem 51: 55-61. Find this article online

Hoch JA (1993) Regulation of the phosphorelay and the initiation of sporulation in Bacillus subtilis. Annu Rev Microbiol 47: 441-465. Find this article online

Molle V, Fujita M, Jensen ST, Eichenberger P, Gonzalez-Pastor JE, et al. (2003) The Spo0A regulon of Bacillus subtilis. Mol Microbiol 50: 1683-1701. Find this article online

Tzeng YL, Zhou XZ, Hoch JA (1998) Phosphorylation of the Spo0B response regulator phosphotransferase of the phosphorelay initiating development in Bacillus subtilis. J Biol Chem 273: 23849-23855. Find this article online

Lewandoski M, Dubnau E, Smith I (1986) Transcriptional regulation of the spo0F gene of Bacillus subtilis. J Bacteriol 168: 870-877. Find this article online

Strauch MA, Wu JJ, Jonas RH, Hoch JA (1993) A positive feedback loop controls transcription of the spoOF gene, a component of the sporulation phosphorelay in Bacillus subtilis. Mol Microbiol 7: 967-974. Find this article online

Bai U, Lewandoski M, Dubnau E, Smith I (1990) Temporal regulation of the Bacillus subtilis early sporulation gene spo0F. J Bacteriol 172: 5432-5439. Find this article online

Arigoni F, Duncan L, Alper S, Losick R, Stragier P (1996) SpoIIE governs the phosphorylation state of a protein regulating transcription factor sigma F during sporulation in Bacillus subtilis. Proc Natl Acad Sci U S A 93: 3238-3242. Find this article online

Arigoni F, Pogliano K, Webb CD, Stragier P, Losick R (1995) Localization of protein implicated in establishment of cell type to sites of asymmetric division. Science 270: 637-640. Find this article online

Barak I, Behari J, Olmedo G, Guzman P, Brown DP, et al. (1996) Structure and function of the Bacillus SpoIIE protein and its localization to sites of sporulation septum assembly. Mol Microbiol 19: 1047-1060. Find this article online

Duncan L, Alper S, Arigoni F, Losick R, Stragier P (1995) Activation of cell-specific transcription by a serine phosphatase at the site of asymmetric division. Science 270: 641-644. Find this article online

Khvorova A, Zhang L, Higgins ML, Piggot PJ (1998) The spoIIE locus is involved in the Spo0A-dependent switch in the location of FtsZ rings in Bacillus subtilis. J Bacteriol 180: 1256-1260. Find this article online

Eichenberger P, Fujita M, Jensen ST, Conlon EM, Rudner DZ, et al. (2004) The program of gene transcription for a single differentiating cell type during sporulation in Bacillus subtilis. PLoS Biol 2: e328. Find this article online

Karow ML, Glaser P, Piggot PJ (1995) Identification of a gene, spoIIR, that links the activation of sigma E to the transcriptional activity of sigma F during sporulation in Bacillus subtilis. Proc Natl Acad Sci U S A 92: 2012-2016. Find this article online

Londono-Vallejo JA, Stragier P (1995) Cell-cell signaling pathway activating a developmental transcription factor in Bacillus subtilis. Genes Dev 9: 503-508. Find this article online

Hofmeister AE, Londono-Vallejo A, Harry E, Stragier P, Losick R (1995) Extracellular signal protein triggering the proteolytic activation of a developmental transcription factor in B. subtilis. Cell 83: 219-226. Find this article online

Losick R, Stragier P (1992) Crisscross regulation of cell-type-specific gene expression during development in B. subtilis. Nature 355: 601-604. Find this article online

Phillips ZE, Strauch MA (2002) Bacillus subtilis sporulation and stationary phase gene expression. Cell Mol Life Sci 59: 392-402. Find this article online

Campo N, Marquis KA, Rudner DZ (2008) SpoIIQ anchors membrane proteins on both sides of the sporulation septum in Bacillus subtilis. J Biol Chem 283: 4975-4982. Find this article online

de Jong IG, Veening JW, Kuipers OP (2010) Heterochronic phosphorelay gene expression as a source of heterogeneity in Bacillus subtilis spore formation. J Bacteriol 192: 2053-2067. Find this article online

Chastanet A, Vitkup D, Yuan GC, Norman TM, Liu JS, et al. (2010) Broadly heterogeneous activation of the master regulator for sporulation in Bacillus subtilis. Proc Natl Acad Sci U S A 107: 8486-8491. Find this article online

Cai L, Friedman N, Xie XS (2006) Stochastic protein expression in individual cells at the single molecule level. Nature 440: 358-362. Find this article online

Dubnau D, Losick R (2006) Bistability in bacteria. Mol Microbiol 61: 564-572. Find this article online

Zhang XP, Liu F, Cheng Z, Wang W (2009) Cell fate decision mediated by p53 pulses. Proc Natl Acad Sci U S A 106: 12245-12250. Find this article online

Schultz D, Ben Jacob E, Onuchic JN, Wolynes PG (2007) Molecular level stochastic model for competence cycles in Bacillus subtilis. Proc Natl Acad Sci U S A 104: 17582-17587. Find this article online

Schultz D, Onuchic JN, Wolynes PG (2007) Understanding stochastic simulations of the smallest genetic networks. J Chem Phys 126: 245102. Find this article online

Blake WJ, Balazsi G, Kohanski MA, Isaacs FJ, Murphy KF, et al. (2006) Phenotypic consequences of promoter-mediated transcriptional noise. Mol Cell 24: 853-865. Find this article online

Balaban NQ, Merrin J, Chait R, Kowalik L, Leibler S (2004) Bacterial persistence as a phenotypic switch. Science 305: 1622-1625. Find this article online

Zeng L, Skinner SO, Zong C, Sippy J, Feiss M, et al. (2010) Decision making at a subcellular level determines the outcome of bacteriophage infection. Cell 141: 682-691. Find this article online

Thattai M, van Oudenaarden A (2004) Stochastic gene expression in fluctuating environments. Genetics 167: 523-530. Find this article online

Kussell E, Leibler S (2005) Phenotypic diversity, population growth, and information in fluctuating environments. Science 309: 2075-2078. Find this article online

Wolf DM, Vazirani VV, Arkin AP (2005) Diversity in times of adversity: probabilistic strategies in microbial survival games. J Theor Biol 234: 227-253. Find this article online

Wolf DM, Vazirani VV, Arkin AP (2005) A microbial modified prisoner's dilemma game: how frequency-dependent selection can lead to random phase variation. J Theor Biol 234: 255-262. Find this article online

Choi HS, Han S, Yokota H, Cho KH (2007) Coupled positive feedbacks provoke slow induction plus fast switching in apoptosis. FEBS Lett 581: 2684-2690. Find this article online

Black F, Myron S (1973) The Pricing of Options and Corporate Liabilities. J Polit Econ 81: 637-654. Find this article online

Haima P, Bron S, Venema G (1987) The effect of restriction on shotgun cloning and plasmid stability in Bacillus subtilis Marburg. Mol Gen Genet 209: 335-342. Find this article online

Rosenfeld N, Young JW, Alon U, Swain PS, Elowitz MB (2005) Gene regulation at the single-cell level. Science 307: 1962-1965. Find this article online