Reversible and noisy progression towards a commitment point enables adaptable and reliable cellular decision-making

Stem cell therapy for heart failure in the clinics: new perspectives in the era of precision medicine and artificial intelligence

Stem/progenitor cells have been widely evaluated as a promising therapeutic option for heart failure (HF). Numerous clinical trials with stem/progenitor cell-based therapy (SCT) for HF have demonstrated encouraging results, but not without limitations or discrepancies. Recent technological advancements in multiomics, bioinformatics, precision medicine, artificial intelligence (AI), and machine learning (ML) provide new approaches and insights for stem cell research and therapeutic development. Integration of these new technologies into stem/progenitor cell therapy for HF may help address: 1) the technical challenges to obtain reliable and high-quality therapeutic precursor cells, 2) the discrepancies between preclinical and clinical studies, and 3) the personalized selection of optimal therapeutic cell types/populations for individual patients in the context of precision medicine. This review summarizes the current status of SCT for HF in clinics and provides new perspectives on the development of computation-aided SCT in the era of precision medicine and AI/ML.

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.

A circadian clock in a nonphotosynthetic prokaryote

Circadian clocks create a 24-hour temporal structure, which allows organisms to occupy a niche formed by time rather than space. They are pervasive throughout nature, yet they remain unexpectedly unexplored and uncharacterized in nonphotosynthetic bacteria. Here, we identify in Bacillus subtilis circadian rhythms sharing the canonical properties of circadian clocks: free-running period, entrainment, and temperature compensation. We show that gene expression in B. subtilis can be synchronized in 24-hour light or temperature cycles and exhibit phase-specific characteristics of entrainment. Upon release to constant dark and temperature conditions, bacterial biofilm populations have temperature-compensated free-running oscillations with a period close to 24 hours. Our work opens the field of circadian clocks in the free-living, nonphotosynthetic prokaryotes, bringing considerable potential for impact upon biomedicine, ecology, and industrial processes.

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

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

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

Forecasting cell fate during antibiotic exposure using stochastic gene expression

Antibiotic killing does not occur at a single, precise time for all cells within a population. Variability in time to death can be caused by stochastic expression of genes, resulting in differences in endogenous stress-resistance levels between individual cells in a population. Here we investigate whether single-cell differences in gene expression prior to antibiotic exposure are related to cell survival times after antibiotic exposure for a range of genes of diverse function. We quantified the time to death of single cells under antibiotic exposure in combination with expression of reporters. For some reporters, including genes involved in stress response and cellular processes like metabolism, the time to cell death had a strong relationship with the initial expression level of the genes. Our results highlight the single-cell level non-uniformity of antibiotic killing and also provide examples of key genes where cell-to-cell variation in expression is strongly linked to extended durations of antibiotic survival.

Time Series Analysis of the Bacillus subtilis Sporulation Network Reveals Low Dimensional Chaotic Dynamics

Chaotic behavior refers to a behavior which, albeit irregular, is generated by an underlying deterministic process. Therefore, a chaotic behavior is potentially controllable. This possibility becomes practically amenable especially when chaos is shown to be low-dimensional, i.e., to be attributable to a small fraction of the total systems components. In this case, indeed, including the major drivers of chaos in a system into the modeling approach allows us to improve predictability of the systems dynamics. Here, we analyzed the numerical simulations of an accurate ordinary differential equation model of the gene network regulating sporulation initiation in Bacillus subtilis to explore whether the non-linearity underlying time series data is due to low-dimensional chaos. Low-dimensional chaos is expectedly common in systems with few degrees of freedom, but rare in systems with many degrees of freedom such as the B. subtilis sporulation network. The estimation of a number of indices, which reflect the chaotic nature of a system, indicates that the dynamics of this network is affected by deterministic chaos. The neat separation between the indices obtained from the time series simulated from the model and those obtained from time series generated by Gaussian white and colored noise confirmed that the B. subtilis sporulation network dynamics is affected by low dimensional chaos rather than by noise. Furthermore, our analysis identifies the principal driver of the networks chaotic dynamics to be sporulation initiation phosphotransferase B (Spo0B). We then analyzed the parameters and the phase space of the system to characterize the instability points of the network dynamics, and, in turn, to identify the ranges of values of Spo0B and of the other drivers of the chaotic dynamics, for which the whole system is highly sensitive to minimal perturbation. In summary, we described an unappreciated source of complexity in the B. subtilis sporulation network by gathering evidence for the chaotic behavior of the system, and by suggesting candidate molecules driving chaos in the system. The results of our chaos analysis can increase our understanding of the intricacies of the regulatory network under analysis, and suggest experimental work to refine our behavior of the mechanisms underlying B. subtilis sporulation initiation control.

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.

Visualizing retinoic acid morphogen gradients

Morphogens were originally defined as secreted signaling molecules that diffuse from local sources to form concentration gradients, which specify multiple cell fates. More recently morphogen gradients have been shown to incorporate a range of mechanisms including short-range signal activation, transcriptional/translational feedback, and temporal windows of target gene induction. Many critical cell-cell signals implicated in both embryonic development and disease, such as Wnt, fibroblast growth factor (Fgf), hedgehog (Hh), transforming growth factor beta (TGFb), and retinoic acid (RA), are thought to act as morphogens, but key information on signal propagation and ligand distribution has been lacking for most. The zebrafish provides unique advantages for genetics and imaging to address gradients during early embryonic stages when morphogens help establish major body axes. This has been particularly informative for RA, where RA response elements (RAREs) driving fluorescent reporters as well as Fluorescence Resonance Energy Transfer (FRET) reporters of receptor binding have provided evidence for gradients, as well as regulatory mechanisms that attenuate noise and enhance gradient robustness in vivo. Here we summarize available tools in zebrafish and discuss their utility for studying dynamic regulation of RA morphogen gradients, through combined experimental and computational approaches.

Regulatory mechanisms link phenotypic plasticity to evolvability

Organisms have a remarkable capacity to respond to environmental change. They can either respond directly, by means of phenotypic plasticity, or they can slowly adapt through evolution. Yet, how phenotypic plasticity links to evolutionary adaptability is largely unknown. Current studies of plasticity tend to adopt a phenomenological reaction norm (RN) approach, which neglects the mechanisms underlying plasticity. Focusing on a concrete question - the optimal timing of bacterial sporulation - we here also consider a mechanistic approach, the evolution of a gene regulatory network (GRN) underlying plasticity. Using individual-based simulations, we compare the RN and GRN approach and find a number of striking differences. Most importantly, the GRN model results in a much higher diversity of responsive strategies than the RN model. We show that each of the evolved strategies is pre-adapted to a unique set of unseen environmental conditions. The regulatory mechanisms that control plasticity therefore critically link phenotypic plasticity to the adaptive potential of biological populations.

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.

Fluctuations, or “noise”, in the response of a system is usually thought to be harmful. However, it is becoming increasingly clear that in single-celled organisms, noise can sometimes help cells survive. This is because noise can enhance the diversity of responses within a cell population. In this study, we identify a novel benefit of noise in the competence response of a population of Bacillus subtilis bacteria, where competence is the ability of bacteria to take in DNA from their environment when under stress. We use computational modeling and experiments to show that noise increases the range of stress levels for which these bacteria exhibit a highly dynamic response, meaning that they are neither unresponsive, nor permanently in the competent state. Since a dynamic response is thought to be optimal for survival, this study suggests that noise is exploited to increase the fitness of the bacterial population.

Coordination of cell decisions and promotion of phenotypic diversity in B. subtilis via pulsed behavior of the phosphorelay

The phosphorelay of Bacillus subtilis, a kinase cascade that activates master regulator Spo0A ~ P in response to starvation signals, is the core of a large network controlling the cell's decision to differentiate into sporulation and other phenotypes. This article reviews recent advances in understanding the origins and purposes of the complex dynamical behavior of the phosphorelay, which pulses with peaks of activity coordinated with the cell cycle. The transient imbalance in the expression of two critical genes caused by their strategic placement at opposing ends of the chromosome proved to be the key for this pulsed behavior. Feedback control loops in the phosphorelay use these pulses to implement a timer mechanism, which creates several windows of opportunity for phenotypic transitions over multiple generations. This strategy allows the cell to coordinate multiple differentiation programs in a decision process that fosters phenotypic diversity and adapts to current conditions.

New Elements To Consider When Modeling the Hazards Associated with Botulinum Neurotoxin in Food

Botulinum neurotoxins (BoNTs) produced by the anaerobic bacterium Clostridium botulinum are the most potent biological substances known to mankind. BoNTs are the agents responsible for botulism, a rare condition affecting the neuromuscular junction and causing a spectrum of diseases ranging from mild cranial nerve palsies to acute respiratory failure and death. BoNTs are a potential biowarfare threat and a public health hazard, since outbreaks of foodborne botulism are caused by the ingestion of preformed BoNTs in food. Currently, mathematical models relating to the hazards associated with C. botulinum, which are largely empirical, make major contributions to botulinum risk assessment. Evaluated using statistical techniques, these models simulate the response of the bacterium to environmental conditions. Though empirical models have been successfully incorporated into risk assessments to support food safety decision making, this process includes significant uncertainties so that relevant decision making is frequently conservative and inflexible. Progression involves encoding into the models cellular processes at a molecular level, especially the details of the genetic and molecular machinery. This addition drives the connection between biological mechanisms and botulism risk assessment and hazard management strategies. This review brings together elements currently described in the literature that will be useful in building quantitative models of C. botulinum neurotoxin production. Subsequently, it outlines how the established form of modeling could be extended to include these new elements. Ultimately, this can offer further contributions to risk assessments to support food safety decision making.

Chance and Necessity in Bacillus subtilis Development

Bacillus subtilis is an important model bacterium for the study of developmental adaptations that enhance survival in the face of fluctuating environmental challenges. These adaptations include sporulation, biofilm formation, motility, cannibalism, and competence. Remarkably, not all the cells in a given population exhibit the same response. The choice of fate by individual cells is random but is also governed by complex signal transduction pathways and cross talk mechanisms that reinforce decisions once made. The interplay of stochastic and deterministic mechanisms governing the selection of developmental fate on the single-cell level is discussed in this article.

A survey of HK, HPt, and RR domains and their organization in two-component systems and phosphorelay proteins of organisms with fully sequenced genomes

Two Component Systems and Phosphorelays (TCS/PR) are environmental signal transduction cascades in prokaryotes and, less frequently, in eukaryotes. The internal domain organization of proteins and the topology of TCS/PR cascades play an important role in shaping the responses of the circuits. It is thus important to maintain updated censuses of TCS/PR proteins in order to identify the various topologies used by nature and enable a systematic study of the dynamics associated with those topologies. To create such a census, we analyzed the proteomes of 7,609 organisms from all domains of life with fully sequenced and annotated genomes. To begin, we survey each proteome searching for proteins containing domains that are associated with internal signal transmission within TCS/PR: Histidine Kinase (HK), Response Regulator (RR) and Histidine Phosphotranfer (HPt) domains, and analyze how these domains are arranged in the individual proteins. Then, we find all types of operon organization and calculate how much more likely are proteins that contain TCS/PR domains to be coded by neighboring genes than one would expect from the genome background of each organism. Finally, we analyze if the fusion of domains into single TCS/PR proteins is more frequently observed than one might expect from the background of each proteome. We find 50 alternative ways in which the HK, HPt, and RR domains are observed to organize into single proteins. In prokaryotes, TCS/PR coding genes tend to be clustered in operons. 90% of all proteins identified in this study contain just one of the three domains, while 8% of the remaining proteins combine one copy of an HK, a RR, and/or an HPt domain. In eukaryotes, 25% of all TCS/PR proteins have more than one domain. These results might have implications for how signals are internally transmitted within TCS/PR cascades. These implications could explain the selection of the various designs in alternative circumstances.

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.

How long can bacteriophage λ change its mind?

A key event in the lifecycle of a temperate bacteriophage is the choice between lysis and lysogeny upon infection of a susceptible host cell. In a recent paper, we showed that a prolonged period exists after the decision to lysogenize, during which bacteriophage λ can abandon the initial decision, and instead develop lytically, as a response to the accumulation of the late lytic regulatory protein Q. Here, we present evidence that expression of Q does not induce replication of λ DNA, suggesting that the DNA to be packaged into the resulting phage progeny was already present at the time of the initial decision to lysogenize. We summarize our findings in a working model of the key determinants of the duration of the post-decision period during which it is possible for the infected cell to switch from the lysogeny decision to successful lytic development.

MetDFBA: incorporating time-resolved metabolomics measurements into dynamic flux balance analysis

This paper presents MetDFBA, a new approach incorporating experimental metabolomics time-series into constraint-based modeling. The method can be used for hypothesis testing and predicting dynamic flux profiles.

Computational modelling and analysis of the molecular network regulating sporulation initiation in Bacillus subtilis

BACKGROUND

Bacterial spores are important contaminants in food, and the spore forming bacteria are often implicated in food safety and food quality considerations. Spore formation is a complex developmental process involving the expression of more than 500 genes over the course of 6 to 8 hrs. The process culminates in the formation of resting cells capable of resisting environmental extremes and remaining dormant for long periods of time, germinating when conditions promote further vegetative growth. Experimental observations of sporulation and germination are problematic and time consuming so that reliable models are an invaluable asset in terms of prediction and risk assessment. In this report we develop a model which assists in the interpretation of sporulation dynamics.

RESULTS

This paper defines and analyses a mathematical model for the network regulating Bacillus subtilis sporulation initiation, from sensing of sporulation signals down to the activation of the early genes under control of the master regulator Spo0A. Our model summarises and extends other published modelling studies, by allowing the user to execute sporulation initiation in a scenario where Isopropyl β-D-1-thiogalactopyranoside (IPTG) is used as an artificial sporulation initiator as well as in modelling the induction of sporulation in wild-type cells. The analysis of the model results and the comparison with experimental data indicate that the model is good at predicting inducible responses to sporulation signals. However, the model is unable to reproduce experimentally observed accumulation of phosphorelay sporulation proteins in wild type B. subtilis. This model also highlights that the phosphorelay sub-component, which relays the signals detected by the sensor kinases to the master regulator Spo0A, is crucial in determining the response dynamics of the system.

CONCLUSION

We show that there is a complex connectivity between the phosphorelay features and the master regulatory Spo0A. Additional we discovered that the experimentally observed regulation of the phosphotransferase Spo0B for wild-type B. subtilis may be playing an important role in the network which suggests that modelling of sporulation initiation may require additional experimental support.

Advancing stem cell therapy from bench to bedside: lessons from drug therapies

The inadequacy of existing therapeutic tools together with the paucity of organ donors have always led medical researchers to innovate the current treatment methods or to discover new ways to cure disease. Emergence of cell-based therapies has provided a new framework through which it has given the human world a new hope. Though relatively a new concept, the pace of advancement clearly reveals the significant role that stem cells will ultimately play in the near future. However, there are numerous uncertainties that are prevailing against the present setting of clinical trials related to stem cells: like the best route of cell administration, appropriate dosage, duration and several other applications. A better knowledge of these factors can substantially improve the effectiveness of disease cure or organ repair using this latest therapeutic tool. From a certain perspective, it could be argued that by considering certain proven clinical concepts and experience from synthetic drug system, we could improve the overall efficacy of cell-based therapies. In the past, studies on synthetic drug therapies and their clinical trials have shown that all the aforementioned factors have critical ascendancy over its therapeutic outcomes. Therefore, based on the knowledge gained from synthetic drug delivery systems, we hypothesize that by employing many of the clinical approaches from synthetic drug therapies to this new regenerative therapeutic tool, the efficacy of stem cell-based therapies can also be improved.

Commitment to lysogeny is preceded by a prolonged period of sensitivity to the late lytic regulator Q in bacteriophage λ

A key event in development is the irreversible commitment to a particular cell fate, which may be concurrent with or delayed with respect to the initial cell fate decision. In this work, we use the paradigmatic bacteriophage λ lysis-lysogeny decision circuit to study the timing of commitment. The lysis-lysogeny decision is made based on the expression trajectory of CII. The chosen developmental strategy is manifested by repression of the pR and pL promoters by CI (lysogeny) or by antitermination of late gene expression by Q (lysis). We found that expression of Q in trans from a plasmid at the time of infection resulted in a uniform lytic decision. Furthermore, expression of Q up to 50 min after infection results in lysis of the majority of cells which initially chose lysogenic development. In contrast, expression of Q in cells containing a single chromosomal prophage had no effect on cell growth, indicating commitment to lysogeny. Notably, if the prophage was present in 10 plasmid-borne copies, Q expression resulted in lytic development, suggesting that the cellular phage chromosome number is the critical determinant of the timing of lysogenic commitment. Based on our results, we conclude that (i) the lysogenic decision made by the CI-Cro switch soon after infection can be overruled by ectopic Q expression at least for a time equivalent to one phage life cycle, (ii) the presence of multiple λ chromosomes is a prerequisite for a successful Q-mediated switch from lysogenic to lytic development, and (iii) phage chromosomes within the same cell can reach different decisions.

Nutrient-regulated proteolysis of MrpC halts expression of genes important for commitment to sporulation during Myxococcus xanthus development

Starved Myxococcus xanthus cells glide to aggregation centers and form fruiting bodies in which rod-shaped cells differentiate into ovoid spores. Commitment to development was investigated by adding nutrients at specific times after starvation and determining whether development halted or proceeded. At 24 h poststarvation, some rod-shaped cells were committed to subsequent shape change and to becoming sonication-resistant spores, but nutrients caused partial disaggregation of fruiting bodies. By 30 h poststarvation, 10-fold more cells were committed to becoming sonication-resistant spores, and compact fruiting bodies persisted after nutrient addition. During the critical period of commitment around 24 to 30 h poststarvation, the transcription factors MrpC and FruA cooperatively regulate genes important for sporulation. FruA responds to short-range C-signaling, which increases as cells form fruiting bodies. MrpC was found to be highly sensitive to nutrient-regulated proteolysis both before and during the critical period of commitment to sporulation. The rapid turnover of MrpC upon nutrient addition to developing cells halted expression of the dev operon, which is important for sporulation. Regulated proteolysis of MrpC appeared to involve ATP-independent metalloprotease activity and may provide a mechanism for monitoring whether starvation persists and halting commitment to sporulation if nutrients reappear.

Tunable stochastic pulsing in the Escherichia coli multiple antibiotic resistance network from interlinked positive and negative feedback loops

Cells live in uncertain, dynamic environments and have many mechanisms for sensing and responding to changes in their surroundings. However, sudden fluctuations in the environment can be catastrophic to a population if it relies solely on sensory responses, which have a delay associated with them. Cells can reconcile these effects by using a tunable stochastic response, where in the absence of a stressor they create phenotypic diversity within an isogenic population, but use a deterministic response when stressors are sensed. Here, we develop a stochastic model of the multiple antibiotic resistance network of Escherichia coli and show that it can produce tunable stochastic pulses in the activator MarA. In particular, we show that a combination of interlinked positive and negative feedback loops plays an important role in setting the dynamics of the stochastic pulses. Negative feedback produces a pulsatile response that is tunable, while positive feedback serves to amplify the effect. Our simulations show that the uninduced native network is in a parameter regime that is of low cost to the cell (taxing resistance mechanisms are expressed infrequently) and also elevated noise strength (phenotypic variability is high). The stochastic pulsing can be tuned by MarA induction such that variability is decreased once stresses are sensed, avoiding the detrimental effects of noise when an optimal MarA concentration is needed. We further show that variability in the expression of MarA can act as a bet hedging mechanism, allowing for survival in time-varying stress environments, however this effect is tunable to allow for a fully induced, deterministic response in the presence of a stressor.

Modeling gene expression in time and space

Cell populations rarely exhibit gene-expression profiles that are homogeneous in time and space. In the temporal domain, dynamical behaviors such as oscillations and pulses of protein production pervade cell biology, underlying phenomena as diverse as circadian rhythmicity, cell cycle control, stress and damage responses, and stem-cell pluripotency. In multicellular populations, spatial heterogeneities are crucial for decision making and development, among many other functions. Cells need to exquisitely coordinate this temporal and spatial variation to survive. Although the spatiotemporal character of gene expression is challenging to quantify experimentally at the level of individual cells, it is beneficial from the modeling viewpoint, because it provides strong constraints that can be probed by theoretically analyzing mathematical models of candidate gene and protein circuits. Here, we review recent examples of temporal dynamics and spatial patterning in gene expression to show how modeling such phenomenology can help us unravel the molecular mechanisms of cellular function.

Dynamics and precision in retinoic acid morphogen gradients

Retinoic acid (RA) regulates many cellular behaviors during embryonic development and adult homeostasis. Like other morphogens, RA forms gradients through the use of localized sources and sinks, feedback, and interactions with other signals; this has been particularly well studied in the context of hindbrain segmentation in vertebrate embryos. Yet, as a small lipophilic molecule derived from a dietary source-vitamin A-RA differs markedly from better-studied polypeptide morphogens in its mechanisms of transport, signaling, and removal. Computational models suggest that the distinctive features of RA gradients make them particularly robust to large perturbations. Such features include combined positive and negative feedback effects via intracellular fatty acid binding proteins and RA-degrading enzymes. Here, we discuss how these features, together with feedback interactions among RA target genes, help enable RA to specify multiple, accurate pattern elements in the developing hindbrain, despite operating in an environment of high cellular and biochemical uncertainty and noise.

Ultrasensitivity of the Bacillus subtilis sporulation decision

Starving Bacillus subtilis cells execute a gene expression program resulting in the formation of stress-resistant spores. Sporulation master regulator, Spo0A, is activated by a phosphorelay and controls the expression of a multitude of genes, including the forespore-specific sigma factor σ(F) and the mother cell-specific sigma factor σ(E). Identification of the system-level mechanism of the sporulation decision is hindered by a lack of direct control over Spo0A activity. This limitation can be overcome by using a synthetic system in which Spo0A activation is controlled by inducing expression of phosphorelay kinase KinA. This induction results in a switch-like increase in the number of sporulating cells at a threshold of KinA. Using a combination of mathematical modeling and single-cell microscopy, we investigate the origin and physiological significance of this ultrasensitive threshold. The results indicate that the phosphorelay is unable to achieve a sufficiently fast and ultrasensitive response via its positive feedback architecture, suggesting that the sporulation decision is made downstream. In contrast, activation of σ(F) in the forespore and of σ(E) in the mother cell compartments occurs via a cascade of coherent feed-forward loops, and thereby can produce fast and ultrasensitive responses as a result of KinA induction. Unlike σ(F) activation, σ(E) activation in the mother cell compartment only occurs above the KinA threshold, resulting in completion of sporulation. Thus, ultrasensitive σ(E) activation explains the KinA threshold for sporulation induction. We therefore infer that under uncertain conditions, cells initiate sporulation but postpone making the sporulation decision to average stochastic fluctuations and to achieve a robust population response.

Dynamical principles of cell-cycle arrest: reversible, irreversible, and mixed strategies

Living cells often alternate between proliferating and nonproliferating states as part of individual or collective strategies to adapt to complex and changing environments. To this aim, they have evolved a biochemical regulatory network enabling them to switch between cell-division cycles (i.e., oscillatory state) and cell-cycle arrests (i.e., steady state) in response to extracellular cues. This can be achieved by means of a variety of bifurcation mechanisms that potentially give rise to qualitatively distinct cell-cycle arrest properties. In this paper, we study the dynamics of a minimal biochemical network model in which a cell-division oscillator and a differentiation switch mutually antagonize. We identify the existence of three biologically plausible bifurcation scenarios organized around a codimension-four swallowtail-homoclinic singularity. As a result, the model exhibits a broad repertoire of cell-cycle arrest properties in terms of reversibility of these arrests, tunability of interdivision time, and ability to track time-varying signals. This dynamic versatility would explain the diversity of cell-cycle arrest strategies developed in different living species and functional contexts.

Strategic cell-cycle regulatory features that provide mammalian cells with tunable G1 length and reversible G1 arrest

Transitions between consecutive phases of the eukaryotic cell cycle are driven by the catalytic activity of selected sets of cyclin-dependent kinases (Cdks). Yet, their occurrence and precise timing is tightly scheduled by a variety of means including Cdk association with inhibitory/adaptor proteins (CKIs). Here we focus on the regulation of G1-phase duration by the end of which cells of multicelled organisms must decide whether to enter S phase or halt, and eventually then, differentiate, senesce or die to obey the homeostatic rules of their host. In mammalian cells, entry in and progression through G1 phase involve sequential phosphorylation and inactivation of the retinoblastoma Rb proteins, first, by cyclin D-Cdk4,6 with the help of CKIs of the Cip/Kip family and, next, by the cyclin E-Cdk2 complexes that are negatively regulated by Cip/Kip proteins. Using a dynamical modeling approach, we show that the very way how the Rb and Cip/Kip regulatory modules interact differentially with cyclin D-Cdk4,6 and cyclin E-Cdk2 provides to mammalian cells a powerful means to achieve an exquisitely-sensitive control of G1-phase duration and fully reversible G1 arrests. Consistently, corruption of either one of these two modules precludes G1 phase elongation and is able to convert G1 arrests from reversible to irreversible. This study unveils fundamental design principles of mammalian G1-phase regulation that are likely to confer to mammalian cells the ability to faithfully control the occurrence and timing of their division process in various conditions.

Noise drives sharpening of gene expression boundaries in the zebrafish hindbrain

Morphogens provide positional information for spatial patterns of gene expression during development. However, stochastic effects such as local fluctuations in morphogen concentration and noise in signal transduction make it difficult for cells to respond to their positions accurately enough to generate sharp boundaries between gene expression domains. During development of rhombomeres in the zebrafish hindbrain, the morphogen retinoic acid (RA) induces expression of hoxb1a in rhombomere 4 (r4) and krox20 in r3 and r5. Fluorescent in situ hybridization reveals rough edges around these gene expression domains, in which cells co-express hoxb1a and krox20 on either side of the boundary, and these sharpen within a few hours. Computational analysis of spatial stochastic models shows, surprisingly, that noise in hoxb1a/krox20 expression actually promotes sharpening of boundaries between adjacent segments. In particular, fluctuations in RA initially induce a rough boundary that requires noise in hoxb1a/krox20 expression to sharpen. This finding suggests a novel noise attenuation mechanism that relies on intracellular noise to induce switching and coordinate cellular decisions during developmental patterning.