A mechanistic stochastic framework for regulating bacterial cell division

A generalized adder for cell size homeostasis: Effects on stochastic clonal proliferation

Measurements of cell size dynamics have revealed phenomenological principles by which individual cells control their size across diverse organisms. One of the emerging paradigms of cell size homeostasis is the adder, where the cell cycle duration is established such that the cell size increase from birth to division is independent of the newborn cell size. We provide a mechanistic formulation of the adder, considering that cell size follows any arbitrary nonexponential growth law. Our results show that the main requirement to obtain an adder regardless of the growth law (the time derivative of cell size) is that cell cycle regulators are produced at a rate proportional to the growth law, and cell division is triggered when these molecules reach a prescribed threshold level. Among the implications of this generalized adder, we investigate fluctuations in the proliferation of single-cell-derived colonies. Considering exponential cell size growth, random fluctuations in clonal size show a transient increase and then eventually decay to zero over time (i.e., clonal populations become asymptotically more similar). In contrast, several forms of nonexponential cell size dynamics (with adder-based cell size control) yield qualitatively different results: clonal size fluctuations monotonically increase over time, reaching a nonzero value. These results characterize the interplay between cell size homeostasis at the single-cell level and clonal proliferation at the population level, explaining the broad fluctuations in clonal sizes seen in barcoded human cell lines.

Sequestration of gene products by decoys enhances precision in the timing of intracellular events

Expressed gene products often interact ubiquitously with binding sites at nucleic acids and macromolecular complexes, known as decoys. The binding of transcription factors (TFs) to decoys can be crucial in controlling the stochastic dynamics of gene expression. Here, we explore the impact of decoys on the timing of intracellular events, as captured by the time taken for the levels of a given TF to reach a critical threshold level, known as the first passage time (FPT). Although nonlinearity introduced by binding makes exact mathematical analysis challenging, employing suitable approximations and reformulating FPT in terms of an alternative variable, we analytically assess the impact of decoys. The stability of the decoy-bound TFs against degradation impacts FPT statistics crucially. Decoys reduce noise in FPT, and stable decoy-bound TFs offer greater timing precision with less expression cost than their unstable counterparts. Interestingly, when both bound and free TFs decay at the same rate, decoy binding does not directly alter FPT noise. We verify these results by performing exact stochastic simulations. These results have important implications for the precise temporal scheduling of events involved in the functioning of biomolecular clocks, development processes, cell-cycle control, and cell-size homeostasis.

Bacterial cell size modulation along the growth curve across nutrient conditions

Under stable growth conditions, bacteria maintain cell size homeostasis through coordinated elongation and division. However, fluctuations in nutrient availability result in dynamic regulation of the target cell size. Using microscopy imaging and mathematical modelling, we examine how bacterial cell volume changes over the growth curve in response to nutrient conditions. We find that two rod-shaped bacteria, Escherichia coli and Salmonella enterica, exhibit similar cell volume distributions in stationary phase cultures irrespective of growth media. Cell resuspension in rich media results in a transient peak with a five-fold increase in cell volume ≈ 2h after resuspension. This maximum cell volume, which depends on nutrient composition, subsequently decreases to the stationary phase cell size. Continuous nutrient supply sustains the maximum volume. In poor nutrient conditions, cell volume shows minimal changes over the growth curve, but a markedly decreased cell width compared to other conditions. The observed cell volume dynamics translate into non-monotonic dynamics in the ratio between biomass (optical density) and cell number (colony-forming units), highlighting their non-linear relationship. Our findings support a heuristic model comparing modulation of cell division relative to growth across nutrient conditions and providing novel insight into the mechanisms of cell size control under dynamic environmental conditions.

A Generalized Adder mechanism for Cell Size Homeostasis: Implications for Stochastic Dynamics of Clonal Proliferation

Measurements of cell size dynamics have revealed phenomenological principles by which individual cells control their size across diverse organisms. One of the emerging paradigms of cell size homeostasis is the adder, where the cell cycle duration is established such that the cell size increase from birth to division is independent of the newborn cell size. We provide a mechanistic formulation of the adder considering that cell size follows any arbitrary non-exponential growth law. Our results show that the main requirement to obtain an adder regardless of the growth law (the time derivative of cell size) is that cell cycle regulators are produced at a rate proportional to the growth law and cell division is triggered when these molecules reach a prescribed threshold level. Among the implications of this generalized adder, we investigate fluctuations in the proliferation of single-cell derived colonies. Considering exponential cell size growth, random fluctuations in clonal size show a transient increase and then eventually decay to zero over time (i.e., clonal populations become asymptotically more similar). In contrast, several forms of non-exponential cell size dynamics (with adder-based cell size control) yield qualitatively different results: clonal size fluctuations monotonically increase over time reaching a non-zero value. These results characterize the interplay between cell size homeostasis at the single-cell level and clonal proliferation at the population level, explaining the broad fluctuations in clonal sizes seen in barcoded human cell lines.

Cell identity revealed by precise cell cycle state mapping links data modalities

Several methods for cell cycle inference from sequencing data exist and are widely adopted. In contrast, methods for classification of cell cycle state from imaging data are scarce. We have for the first time integrated sequencing and imaging derived cell cycle pseudo-times for assigning 449 imaged cells to 693 sequenced cells at an average resolution of 3.4 and 2.4 cells for sequencing and imaging data respectively. Data integration revealed thousands of pathways and organelle features that are correlated with each other, including several previously known interactions and novel associations. The ability to assign the transcriptome state of a profiled cell to its closest living relative, which is still actively growing and expanding opens the door for genotype-phenotype mapping at single cell resolution forward in time.

Mechanisms of cell size regulation in slow-growing Escherichia coli cells: discriminating models beyond the adder

Under ideal conditions, Escherichia coli cells divide after adding a fixed cell size, a strategy known as the adder. This concept applies to various microbes and is often explained as the division that occurs after a certain number of stages, associated with the accumulation of precursor proteins at a rate proportional to cell size. However, under poor media conditions, E. coli cells exhibit a different size regulation. They are smaller and follow a sizer-like division strategy where the added size is inversely proportional to the size at birth. We explore three potential causes for this deviation: degradation of the precursor protein and two models where the propensity for accumulation depends on the cell size: a nonlinear accumulation rate, and accumulation starting at a threshold size termed the commitment size. These models fit the mean trends but predict different distributions given the birth size. To quantify the precision of the models to explain the data, we used the Akaike information criterion and compared them to open datasets of slow-growing E. coli cells in different media. We found that none of the models alone can consistently explain the data. However, the degradation model better explains the division strategy when cells are larger, whereas size-related models (power-law and commitment size) account for smaller cells. Our methodology proposes a data-based method in which different mechanisms can be tested systematically.

Feedback between stochastic gene networks and population dynamics enables cellular decision-making

Phenotypic selection occurs when genetically identical cells are subject to different reproductive abilities due to cellular noise. Such noise arises from fluctuations in reactions synthesizing proteins and plays a crucial role in how cells make decisions and respond to stress or drugs. We propose a general stochastic agent-based model for growing populations capturing the feedback between gene expression and cell division dynamics. We devise a finite state projection approach to analyze gene expression and division distributions and infer selection from single-cell data in mother machines and lineage trees. We use the theory to quantify selection in multi-stable gene expression networks and elucidate that the trade-off between phenotypic switching and selection enables robust decision-making essential for synthetic circuits and developmental lineage decisions. Using live-cell data, we demonstrate that combining theory and inference provides quantitative insights into bet-hedging-like response to DNA damage and adaptation during antibiotic exposure in Escherichia coli.

Effects of Molecular Noise on Cell Size Control

Cells employ control strategies to maintain a stable size. Dividing at a target size (the "sizer" strategy) is thought to produce the tightest size distribution. However, this result follows from phenomenological models that ignore the molecular mechanisms required to implement the strategy. Here we investigate a simple mechanistic model for exponentially growing cells whose division is triggered at a molecular abundance threshold. We find that size noise inherits the molecular noise and is consequently minimized not by the sizer but by the "adder" strategy, where a cell divides after adding a target amount to its birth size. We derive a lower bound on size noise that agrees with publicly available data from six microfluidic studies on Escherichia coli bacteria.

Energy allocation theory for bacterial growth control in and out of steady state

Efficient allocation of energy resources to key physiological functions allows living organisms to grow and thrive in diverse environments and adapt to a wide range of perturbations. To quantitatively understand how unicellular organisms utilize their energy resources in response to changes in growth environment, we introduce a theory of dynamic energy allocation which describes cellular growth dynamics based on partitioning of metabolizable energy into key physiological functions: growth, division, cell shape regulation, energy storage and loss through dissipation. By optimizing the energy flux for growth, we develop the equations governing the time evolution of cell morphology and growth rate in diverse environments. The resulting model accurately captures experimentally observed dependencies of bacterial cell size on growth rate, superlinear scaling of metabolic rate with cell size, and predicts nutrient-dependent trade-offs between energy expended for growth, division, and shape maintenance. By calibrating model parameters with available experimental data for the model organism E. coli, our model is capable of describing bacterial growth control in dynamic conditions, particularly during nutrient shifts and osmotic shocks. The model captures these perturbations with minimal added complexity and our unified approach predicts the driving factors behind a wide range of observed morphological and growth phenomena.

Coupling Cell Size Regulation and Proliferation Dynamics of C. glutamicum Reveals Cell Division Based on Surface Area

Single cells actively coordinate growth and division to regulate their size, yet how this size homeostasis at the single-cell level propagates over multiple generations to impact clonal expansion remains fundamentally unexplored. Classical timer models for cell proliferation (where the duration of the cell cycle is an independent variable) predict that the stochastic variation in colony size will increase monotonically over time. In stark contrast, implementing size control according to adder strategy (where on average a fixed size added from cell birth to division) leads to colony size variations that eventually decay to zero. While these results assume a fixed size of the colony-initiating progenitor cell, further analysis reveals that the magnitude of the intercolony variation in population number is sensitive to heterogeneity in the initial cell size. We validate these predictions by tracking the growth of isogenic microcolonies of Corynebacterium glutamicum in microfluidic chambers. Approximating their cell shape to a capsule, we observe that the degree of random variability in cell size is different depending on whether the cell size is quantified as per length, surface area, or volume, but size control remains an adder regardless of these size metrics. A comparison of the observed variability in the colony population with the predictions suggests that proliferation matches better with a cell division based on the cell surface. In summary, our integrated mathematical-experimental approach bridges the paradigms of single-cell size regulation and clonal expansion at the population levels. This innovative approach provides elucidation of the mechanisms of size homeostasis from the stochastic dynamics of colony size for rod-shaped microbes.

Mechanisms of Cell Size Regulation in Slow-Growing Escherichia coli Cells: Discriminating Models Beyond the Adder

Under ideal conditions, Escherichia coli cells divide after adding a fixed cell size, a strategy known as the adder. This concept applies to various microbes and is often explained as the division that occurs after a certain number of stages, associated with the accumulation of precursor proteins at a rate proportional to cell size. However, under poor media conditions, E. coli cells exhibit a different size regulation. They are smaller and follow a sizer-like division strategy where the added size is inversely proportional to the size at birth. We explore three potential causes for this deviation: precursor protein degradation, nonlinear accumulation rate, and a threshold size termed the commitment size. These models fit mean trends but predict different distributions given the birth size. To validate these models, we used the Akaike information criterion and compared them to open datasets of slow-growing E. coli cells in different media. the degradation model could explain the division strategy for media where cells are larger, while the commitment size model could account for smaller cells. The power-law model, finally, better fits the data at intermediate regimes.

PyEcoLib: a python library for simulating stochastic cell size dynamics

Recently, there has been an increasing need for tools to simulate cell size regulation due to important applications in cell proliferation and gene expression. However, implementing the simulation usually presents some difficulties, as the division has a cycle-dependent occurrence rate. In this article, we gather a recent theoretical framework inPyEcoLib, a python-based library to simulate the stochastic dynamics of the size of bacterial cells. This library can simulate cell size trajectories with an arbitrarily small sampling period. In addition, this simulator can include stochastic variables, such as the cell size at the beginning of the experiment, the cycle duration timing, the growth rate, and the splitting position. Furthermore, from a population perspective, the user can choose between tracking a single lineage or all cells in a colony. They can also simulate the most common division strategies (adder, timer, and sizer) using the division rate formalism and numerical methods. As an example of PyecoLib applications, we explain how to couple size dynamics with gene expression predicting, from simulations, how the noise in protein levels increases by increasing the noise in division timing, the noise in growth rate and the noise in cell splitting position. The simplicity of this library and its transparency about the underlying theoretical framework yield the inclusion of cell size stochasticity in complex models of gene expression.

Dynamic proteome trade-offs regulate bacterial cell size and growth in fluctuating nutrient environments

Bacteria dynamically regulate cell size and growth to thrive in changing environments. While previous studies have characterized bacterial growth physiology at steady-state, a quantitative understanding of bacterial physiology in time-varying environments is lacking. Here we develop a quantitative theory connecting bacterial growth and division rates to proteome allocation in time-varying nutrient environments. In such environments, cell size and growth are regulated by trade-offs between prioritization of biomass accumulation or division, resulting in decoupling of single-cell growth rate from population growth rate. Specifically, bacteria transiently prioritize biomass accumulation over production of division machinery during nutrient upshifts, while prioritizing division over growth during downshifts. When subjected to pulsatile nutrient concentration, we find that bacteria exhibit a transient memory of previous metabolic states due to the slow dynamics of proteome reallocation. This allows for faster adaptation to previously seen environments and results in division control which is dependent on the time-profile of fluctuations.

Beyond the average: An updated framework for understanding the relationship between cell growth, DNA replication, and division in a bacterial system

Our understanding of the bacterial cell cycle is framed largely by population-based experiments that focus on the behavior of idealized average cells. Most famously, the contributions of Cooper and Helmstetter help to contextualize the phenomenon of overlapping replication cycles observed in rapidly growing bacteria. Despite the undeniable value of these approaches, their necessary reliance on the behavior of idealized average cells masks the stochasticity inherent in single-cell growth and physiology and limits their mechanistic value. To bridge this gap, we propose an updated and agnostic framework, informed by extant single-cell data, that quantitatively accounts for stochastic variations in single-cell dynamics and the impact of medium composition on cell growth and cell cycle progression. In this framework, stochastic timers sensitive to medium composition impact the relationship between cell cycle events, accounting for observed differences in the relationship between cell cycle events in slow- and fast-growing cells. We conclude with a roadmap for potential application of this framework to longstanding open questions in the bacterial cell cycle field.

Cellular resource allocation strategies for cell size and shape control in bacteria

Bacteria are highly adaptive microorganisms that thrive in a wide range of growth conditions via changes in cell morphologies and macromolecular composition. How bacterial morphologies are regulated in diverse environmental conditions is a long-standing question. Regulation of cell size and shape implies control mechanisms that couple the growth and division of bacteria to their cellular environment and macromolecular composition. In the past decade, simple quantitative laws have emerged that connect cell growth to proteomic composition and the nutrient availability. However, the relationships between cell size, shape, and growth physiology remain challenging to disentangle and unifying models are lacking. In this review, we focus on regulatory models of cell size control that reveal the connections between bacterial cell morphology and growth physiology. In particular, we discuss how changes in nutrient conditions and translational perturbations regulate the cell size, growth rate, and proteome composition. Integrating quantitative models with experimental data, we identify the physiological principles of bacterial size regulation, and discuss the optimization strategies of cellular resource allocation for size control.

Patterns of interdivision time correlations reveal hidden cell cycle factors

The time taken for cells to complete a round of cell division is a stochastic process controlled, in part, by intracellular factors. These factors can be inherited across cellular generations which gives rise to, often non-intuitive, correlation patterns in cell cycle timing between cells of different family relationships on lineage trees. Here, we formulate a framework of hidden inherited factors affecting the cell cycle that unifies known cell cycle control models and reveals three distinct interdivision time correlation patterns: aperiodic, alternator, and oscillator. We use Bayesian inference with single-cell datasets of cell division in bacteria, mammalian and cancer cells, to identify the inheritance motifs that underlie these datasets. From our inference, we find that interdivision time correlation patterns do not identify a single cell cycle model but generally admit a broad posterior distribution of possible mechanisms. Despite this unidentifiability, we observe that the inferred patterns reveal interpretable inheritance dynamics and hidden rhythmicity of cell cycle factors. This reveals that cell cycle factors are commonly driven by circadian rhythms, but their period may differ in cancer. Our quantitative analysis thus reveals that correlation patterns are an emergent phenomenon that impact cell proliferation and these patterns may be altered in disease.

Effects of antibiotics on bacterial cell morphology and their physiological origins

Characterizing the physiological response of bacterial cells to antibiotic treatment is crucial for the design of antibacterial therapies and for understanding the mechanisms of antibiotic resistance. While the effects of antibiotics are commonly characterized by their minimum inhibitory concentrations or the minimum bactericidal concentrations, the effects of antibiotics on cell morphology and physiology are less well characterized. Recent technological advances in single-cell studies of bacterial physiology have revealed how different antibiotic drugs affect the physiological state of the cell, including growth rate, cell size and shape, and macromolecular composition. Here, we review recent quantitative studies on bacterial physiology that characterize the effects of antibiotics on bacterial cell morphology and physiological parameters. In particular, we present quantitative data on how different antibiotic targets modulate cellular shape metrics including surface area, volume, surface-to-volume ratio, and the aspect ratio. Using recently developed quantitative models, we relate cell shape changes to alterations in the physiological state of the cell, characterized by changes in the rates of cell growth, protein synthesis and proteome composition. Our analysis suggests that antibiotics induce distinct morphological changes depending on their cellular targets, which may have important implications for the regulation of cellular fitness under stress.

Two different cell-cycle processes determine the timing of cell division in Escherichia coli

Cells must control the cell cycle to ensure that key processes are brought to completion. In Escherichia coli, it is controversial whether cell division is tied to chromosome replication or to a replication-independent inter-division process. A recent model suggests instead that both processes may limit cell division with comparable odds in single cells. Here, we tested this possibility experimentally by monitoring single-cell division and replication over multiple generations at slow growth. We then perturbed cell width, causing an increase of the time between replication termination and division. As a consequence, replication became decreasingly limiting for cell division, while correlations between birth and division and between subsequent replication-initiation events were maintained. Our experiments support the hypothesis that both chromosome replication and a replication-independent inter-division process can limit cell division: the two processes have balanced contributions in non-perturbed cells, while our width perturbations increase the odds of the replication-independent process being limiting.

Exponential trajectories, cell size fluctuations, and the adder property in bacteria follow from simple chemical dynamics and division control

Experiments on steady-state bacterial cultures have uncovered several quantitative regularities at the system level. These include, first, the exponential growth of cell size with time and the balanced growth of intracellular chemicals between cell birth and division, which are puzzling given the nonlinear and decentralized chemical dynamics in the cell. We model a cell as a set of chemical populations undergoing nonlinear mass action kinetics in a container whose volume is a linear function of the chemical populations. This turns out to be a special class of dynamical systems that generically has attractors in which all populations grow exponentially with time at the same rate. This explains exponential balanced growth of bacterial cells without invoking any regulatory mechanisms and suggests that this could be a robust property of protocells as well. Second, we consider the hypothesis that cells commit themselves to division when a certain internal chemical population reaches a threshold of N molecules. We show that this hypothesis leads to a simple explanation of some of the variability observed across cells in a bacterial culture. In particular, it reproduces the adder property of cell size fluctuations observed recently in E. coli; the observed correlations among interdivision time, birth volume, and added volume in a generation; and the observed scale of the fluctuations (CV ≈ 10-30%) when N is between 10 and 100. Third, upon including a suitable regulatory mechanism that optimizes the growth rate of the cell, the model reproduces the observed bacterial growth laws including the dependence of the growth rate and ribosomal protein fraction on the medium. Thus, the models provide a framework for unifying diverse aspects of bacterial growth physiology under one roof. They also suggest new questions for experimental and theoretical enquiry.

Cell size distribution of lineage data: analytic results and parameter inference

Recent advances in single-cell technologies have enabled time-resolved measurements of the cell size over several cell cycles. These data encode information on how cells correct size aberrations so that they do not grow abnormally large or small. Here, we formulate a piecewise deterministic Markov model describing the evolution of the cell size over many generations, for all three cell size homeostasis strategies (timer, sizer, and adder). The model is solved to obtain an analytical expression for the non-Gaussian cell size distribution in a cell lineage; the theory is used to understand how the shape of the distribution is influenced by the parameters controlling the dynamics of the cell cycle and by the choice of cell tracking protocol. The theoretical cell size distribution is found to provide an excellent match to the experimental cell size distribution of E. coli lineage data collected under various growth conditions.

A bacterial size law revealed by a coarse-grained model of cell physiology

Universal observations in Biology are sometimes described as “laws”. In E. coli, experimental studies performed over the past six decades have revealed major growth laws relating ribosomal mass fraction and cell size to the growth rate. Because they formalize complex emerging principles in biology, growth laws have been instrumental in shaping our understanding of bacterial physiology. Here, we discovered a novel size law that connects cell size to the inverse of the metabolic proteome mass fraction and the active fraction of ribosomes. We used a simple whole-cell coarse-grained model of cell physiology that combines the proteome allocation theory and the structural model of cell division. This integrated model captures all available experimental data connecting the cell proteome composition, ribosome activity, division size and growth rate in response to nutrient quality, antibiotic treatment and increased protein burden. Finally, a stochastic extension of the model explains non-trivial correlations observed in single cell experiments including the adder principle. This work provides a simple and robust theoretical framework for studying the fundamental principles of cell size determination in unicellular organisms.

Bacteria respond to environmental changes by adjusting their molecular composition, cell size and growth rate. This plasticity is thought to result from years of evolution and to be at least in part optimal for bacterial physiology. Over the past decades, quantitative studies of bacterial growth have revealed simple phenomenological relationships, called “growth laws”, which link cell size and cell composition to the growth rate. Simplified mathematical models of cell physiology are useful tools to gain quantitative understanding of the molecular mechanisms that underlie growth laws. For instance, these models helped explaining how optimal allocation of cellular resource to physiological processes and pathways governs the cell molecular composition in response to specific environmental conditions. In this study, we have extended and integrated existing mathematical models and used experimental data from several recent studies to understand the co-regulation of cell composition, cell size and the cellular growth rate. The model predictions uncovered a novel “size law” that links cell size to the levels of metabolic proteins and the fraction of active ribosomes present in the cell. This work provides a useful theoretical tool and a quantitative basis for understanding mechanistically bacterial physiology as a function of external conditions.

Factors That Affect the Enlargement of Bacterial Protoplasts and Spheroplasts

Cell enlargement is essential for the microinjection of various substances into bacterial cells. The cell wall (peptidoglycan) inhibits cell enlargement. Thus, bacterial protoplasts/spheroplasts are used for enlargement because they lack cell wall. Though bacterial species that are capable of gene manipulation are limited, procedure for bacterial cell enlargement does not involve any gene manipulation technique. In order to prevent cell wall resynthesis during enlargement of protoplasts/spheroplasts, incubation media are supplemented with inhibitors of peptidoglycan biosynthesis such as penicillin. Moreover, metal ion composition in the incubation medium affects the properties of the plasma membrane. Therefore, in order to generate enlarged cells that are suitable for microinjection, metal ion composition in the medium should be considered. Experiment of bacterial protoplast or spheroplast enlargement is useful for studies on bacterial plasma membrane biosynthesis. In this paper, we have summarized the factors that influence bacterial cell enlargement.

Modeling homeostasis mechanisms that set the target cell size

How organisms maintain cell size homeostasis is a long-standing problem that remains unresolved, especially in multicellular organisms. Recent experiments in diverse animal cell types demonstrate that within a cell population, cellular proliferation is low for small and large cells, but high at intermediate sizes. Here we use mathematical models to explore size-control strategies that drive such a non-monotonic profile resulting in the proliferation capacity being maximized at a target cell size. Our analysis reveals that most models of size control yield proliferation capacities that vary monotonically with cell size, and non-monotonicity requires two key mechanisms: (1) the growth rate decreases with increasing size for excessively large cells; and (2) cell division occurs as per the Adder model (division is triggered upon adding a fixed size from birth), or a Sizer-Adder combination. Consistent with theory, Jurkat T cell growth rates increase with size for small cells, but decrease with size for large cells. In summary, our models show that regulation of both growth and cell-division timing is necessary for size control in animal cells, and this joint mechanism leads to a target cell size where cellular proliferation capacity is maximized.

Extended Robust Boolean Network of Budding Yeast Cell Cycle

Background:

How to explore the dynamics of transition probabilities between phases of budding yeast cell cycle (BYCC) network based on the dynamics of protein activities that control this network? How to identify the robust structure of protein interactions of BYCC Boolean network (BN)? Budding yeast allows scientists to put experiments into effect in order to discover the intracellular cell cycle regulating structures which are well simulated by mathematical modeling.

Methods:

We extended an available deterministic BN of proteins responsible for the cell cycle to a Markov chain model containing apoptosis besides G1, S, G2, M, and stationary G1. Using genetic algorithm (GA), we estimated the kinetic parameters of the extended BN model so that the subsequent transition probabilities derived using Markov chain model of cell states as normal cell cycle becomes the maximum while the structure of chemical interactions of extended BN of cell cycle becomes more stable.

Results:

Using kinetic parameters optimized by GA, the probability of the subsequent transitions between cell cycle phases is maximized. The relative basin size of stationary G1 increased from 86% to 96.48% while the number of attractors decreased from 7 in the original model to 5 in the extended one. Hence, an increase in the robustness of the system has been achieved.

Conclusion:

The structure of interacting proteins in cell cycle network affects its robustness and probabilities of transitions between different cell cycle phases. Markov chain and BN are good approaches to study the stability and dynamics of the cell cycle network.

Optimum Threshold Minimizes Noise in Timing of Intracellular Events

How the noisy expression of regulatory proteins affects timing of intracellular events is an intriguing fundamental problem that influences diverse cellular processes. Here we use the bacteriophage λ to study event timing in individual cells where cell lysis is the result of expression and accumulation of a single protein (holin) in the Escherichia coli cell membrane up to a critical threshold level. Site-directed mutagenesis of the holin gene generated phage variants that vary in their lysis times from 30 to 190 min. Observation of the lysis times of single cells reveals an intriguing finding-the noise in lysis timing first decreases with increasing lysis time to reach a minimum and then sharply increases at longer lysis times. A mathematical model with stochastic expression of holin together with dilution from cell growth was sufficient to explain the non-monotonic noise profile and identify holin accumulation thresholds that generate precision in lysis timing.

General quantitative relations linking cell growth and the cell cycle in Escherichia coli

Growth laws emerging from studies of cell populations provide essential constraints on the global mechanisms that coordinate cell growth^1-3. The foundation of bacterial cell cycle studies relies on two interconnected dogmas that were proposed more than 50 years ago-the Schaechter-Maaloe-Kjeldgaard growth law that relates cell mass to growth rate¹ and Donachie's hypothesis of a growth-rate-independent initiation mass⁴. These dogmas spurred many efforts to understand their molecular bases and physiological consequences^5-14. Although they are generally accepted in the fast-growth regime, that is, for doubling times below 1 h, extension of these dogmas to the slow-growth regime has not been consistently achieved. Here, through a quantitative physiological study of Escherichia coli cell cycles over an extensive range of growth rates, we report that neither dogma holds in either the slow- or fast-growth regime. In their stead, linear relations between the cell mass and the rate of chromosome replication-segregation were found across the range of growth rates. These relations led us to propose an integral-threshold model in which the cell cycle is controlled by a licensing process, the rate of which is related in a simple way to chromosomal dynamics. These results provide a quantitative basis for predictive understanding of cell growth-cell cycle relationships.

Unification of cell division control strategies through continuous rate models

Recent experiments support the adder model for E. coli division control. This model posits that bacteria grow, on average, a fixed size before division. It also predicts decorrelation between the noise in the added size and the size at birth. Here we develop a theory based on stochastic hybrid systems which could explain the main division strategies, including not only the adder strategy but the whole range from sizer to timer. We use experiments to explore the division control of E. coli growing with glycerol as carbon source. In this medium, the division strategy is sizerlike, which means that the added size decreases with the size at birth. We found, as our theory predicts, that in a sizerlike strategy the mean added size decreases with the size at birth while the noise in added size increases. We discuss possible molecular mechanisms underlying this strategy and propose a general model that encompasses the different division strategies.

Surface-to-volume scaling and aspect ratio preservation in rod-shaped bacteria

Rod-shaped bacterial cells can readily adapt their lengths and widths in response to environmental changes. While many recent studies have focused on the mechanisms underlying bacterial cell size control, it remains largely unknown how the coupling between cell length and width results in robust control of rod-like bacterial shapes. In this study we uncover a conserved surface-to-volume scaling relation in Escherichia coli and other rod-shaped bacteria, resulting from the preservation of cell aspect ratio. To explain the mechanistic origin of aspect-ratio control, we propose a quantitative model for the coupling between bacterial cell elongation and the accumulation of an essential division protein, FtsZ. This model reveals a mechanism for why bacterial aspect ratio is independent of cell size and growth conditions, and predicts cell morphological changes in response to nutrient perturbations, antibiotics, MreB or FtsZ depletion, in quantitative agreement with experimental data.

Burst statistics in an early biofilm quorum sensing model: the role of spatial colony-growth heterogeneity

Quorum-sensing bacteria in a growing colony of cells send out signalling molecules (so-called “autoinducers”) and themselves sense the autoinducer concentration in their vicinity. Once—due to increased local cell density inside a “cluster” of the growing colony—the concentration of autoinducers exceeds a threshold value, cells in this clusters get “induced” into a communal, multi-cell biofilm-forming mode in a cluster-wide burst event. We analyse quantitatively the influence of spatial disorder, the local heterogeneity of the spatial distribution of cells in the colony, and additional physical parameters such as the autoinducer signal range on the induction dynamics of the cell colony. Spatial inhomogeneity with higher local cell concentrations in clusters leads to earlier but more localised induction events, while homogeneous distributions lead to comparatively delayed but more concerted induction of the cell colony, and, thus, a behaviour close to the mean-field dynamics. We quantify the induction dynamics with quantifiers such as the time series of induction events and burst sizes, the grouping into induction families, and the mean autoinducer concentration levels. Consequences for different scenarios of biofilm growth are discussed, providing possible cues for biofilm control in both health care and biotechnology.

Modulation of first-passage time for gene expression via asymmetric cell division

How to balance the size of exponentially growing cells has always been a focus of biologists. Recent experiments have uncovered that the cell is divided into two daughter cells only when the level of time-keeper protein reaches a fixed threshold and cell division in prokaryote is not completely symmetric. The timing of cell division is essentially random because gene expression is stochastic, but cells seen to manage to have precise timing of cell division events. Although the inter-cellular variability of gene expression has attracted much attention, the randomness of event timing has been rarely studied. In our analysis, the timing of cell division is formulated as the first-passage time (denoted by FPT) for time-keeper protein’s level to cross a critical threshold firstly, we derive exact analytical formulae for the mean and noise of FPT based on stochastic gene expression model with asymmetric cell division. The results of numerical simulation show that the regulatory factors (division rate, newborn cell size, exponential growth rate and threshold) have significant influence on the mean and noise of FPT. We also show that both the increase of division rate and newborn cell size could reduce the mean of FPT and increase the noise of FPT, the larger the exponential growth rate is, the smaller the mean and noise of FPT will be; and the larger the threshold value is, the higher the mean of FPT is and the lower the noise is. In addition, compared with symmetric division, asymmetric division can reduce the mean of FPT and improve the noise of FPT. In summary, our results provide insight into the relationship between regulatory factors and FPT and reveal that asymmetric division is an effective mechanism to shorten the mean of FPT.

Mechanistic Origin of Cell-Size Control and Homeostasis in Bacteria

Evolutionarily divergent bacteria share a common phenomenological strategy for cell-size homeostasis under steady-state conditions. In the presence of inherent physiological stochasticity, cells following this "adder" principle gradually return to their steady-state size by adding a constant volume between birth and division, regardless of their size at birth. However, the mechanism of the adder has been unknown despite intense efforts. In this work, we show that the adder is a direct consequence of two general processes in biology: (1) threshold-accumulation of initiators and precursors required for cell division to a respective fixed number-and (2) balanced biosynthesis-maintenance of their production proportional to volume growth. This mechanism is naturally robust to static growth inhibition but also allows us to "reprogram" cell-size homeostasis in a quantitatively predictive manner in both Gram-negative Escherichia coli and Gram-positive Bacillus subtilis. By generating dynamic oscillations in the concentration of the division protein FtsZ, we were able to oscillate cell size at division and systematically break the adder. In contrast, periodic induction of replication initiator protein DnaA caused oscillations in cell size at initiation but did not alter division size or the adder. Finally, we were able to restore the adder phenotype in slow-growing E. coli, the only known steady-state growth condition wherein E. coli significantly deviates from the adder, by repressing active degradation of division proteins. Together, these results show that cell division and replication initiation are independently controlled at the gene-expression level and that division processes exclusively drive cell-size homeostasis in bacteria. VIDEO ABSTRACT.

Evidence that the human cell cycle is a series of uncoupled, memoryless phases

The cell cycle is canonically described as a series of four consecutive phases: G1, S, G2, and M. In single cells, the duration of each phase varies, but the quantitative laws that govern phase durations are not well understood. Using time-lapse microscopy, we found that each phase duration follows an Erlang distribution and is statistically independent from other phases. We challenged this observation by perturbing phase durations through oncogene activation, inhibition of DNA synthesis, reduced temperature, and DNA damage. Despite large changes in durations in cell populations, phase durations remained uncoupled in individual cells. These results suggested that the independence of phase durations may arise from a large number of molecular factors that each exerts a minor influence on the rate of cell cycle progression. We tested this model by experimentally forcing phase coupling through inhibition of cyclin-dependent kinase 2 (CDK2) or overexpression of cyclin D. Our work provides an explanation for the historical observation that phase durations are both inherited and independent and suggests how cell cycle progression may be altered in disease states.

Controlling organism size by regulating constituent cell numbers

How living cells employ counting mechanisms to regulate their numbers or density is a long-standing problem in developmental biology that ties directly with organism or tissue size. Diverse cells types have been shown to regulate their numbers via secretion of factors in the extracellular space. These factors act as a proxy for the number of cells and function to reduce cellular proliferation rates creating a negative feedback. It is desirable that the production rate of such factors be kept as low as possible to minimize energy costs and detection by predators. Here we formulate a stochastic model of cell proliferation with feedback control via a secreted extracellular factor. Our results show that while low levels of feedback minimizes random fluctuations in cell numbers around a given set point, high levels of feedback amplify Poisson fluctuations in secreted-factor copy numbers. This trade-off results in an optimal feedback strength, and sets a fundamental limit to noise suppression in cell numbers with short-lived factors providing more efficient noise buffering. We further expand the model to consider external disturbances in key physiological parameters, such as, proliferation and factor synthesis rates. Intriguingly, while negative feedback effectively mitigates disturbances in the proliferation rate, it amplifies disturbances in the synthesis rate. In summary, these results provide unique insights into the functioning of feedback-based counting mechanisms, and apply to organisms ranging from unicellular prokaryotes and eukaryotes to human cells.

Homeostasis of protein and mRNA concentrations in growing cells

Many experiments show that the numbers of mRNA and protein are proportional to the cell volume in growing cells. However, models of stochastic gene expression often assume constant transcription rate per gene and constant translation rate per mRNA, which are incompatible with these experiments. Here, we construct a minimal gene expression model to fill this gap. Assuming ribosomes and RNA polymerases are limiting in gene expression, we show that the numbers of proteins and mRNAs both grow exponentially during the cell cycle and that the concentrations of all mRNAs and proteins achieve cellular homeostasis; the competition between genes for the RNA polymerases makes the transcription rate independent of the genome number. Furthermore, by extending the model to situations in which DNA (mRNA) can be saturated by RNA polymerases (ribosomes) and becomes limiting, we predict a transition from exponential to linear growth of cell volume as the protein-to-DNA ratio increases.

Modulation of first-passage time for bursty gene expression via random signals

The stochastic nature of cell-specific signal molecules (such as transcription factor, ribosome, etc.) and the intrinsic stochastic nature of gene expression process result in cell-to-cell variations at protein levels. Increasing experimental evidences suggest that cell phenotypic variations often depend on the accumulation of some special proteins. Hence, a natural and fundamental question is: How does input signal affect the timing of protein count up to a given threshold? To this end, we study effects of input signal on the first-passage time (FPT), the time at which the number of proteins crosses a given threshold. Input signal is distinguished into two types: constant input signal and random input signal, regulating only burst frequency (or burst size) of gene expression. Firstly, we derive analytical formulae for FPT moments in each case of constant signal regulation and random signal regulation. Then, we find that random input signal tends to increases the mean and noise of FPT compared with constant input signal. Finally, we observe that different regulation ways of random signal have different effects on FPT, that is, burst size modulation tends to decrease the mean of FPT and increase the noise of FPT compared with burst frequency modulation. Our findings imply a fundamental mechanism that random fluctuating environment may prolong FPT. This can provide theoretical guidance for studies of some cellular key events such as latency of HIV and lysis time of bacteriophage λ. In conclusion, our results reveal impacts of external signal on FPT and aid understanding the regulation mechanism of gene expression.

A Markovian Approach towards Bacterial Size Control and Homeostasis in Anomalous Growth Processes

Regardless of the progress achieved during recent years, the mechanisms coupling growth and division to attain cell size homeostasis in bacterial populations are still not well understood. In particular, there is a gap of knowledge about the mechanisms controlling anomalous growth events that are ubiquitous even in wild-type phenotypes. Thus, when cells exceed the doubling size the divisome dynamics sets a characteristic length scale that suggests a sizer property. Yet, it has been recently shown that the size at birth and the size increment still satisfy an adder-like correlation. Herein we propose a Markov chain model, that we complement with computational and experimental approaches, to clarify this issue. In this context, we show that classifying cells as a function of the characteristic size set by the divisome dynamics provides a compelling framework to understand size convergence, growth, and division at the large length scale, including the adaptation to, and rescue from, filamentation processes. Our results reveal the independence of size homeostasis on the division pattern of long cells and help to reconcile sizer concepts at the single cell level with an adder-like behavior at a population level.

Modeling Cell Size Regulation: From Single-Cell-Level Statistics to Molecular Mechanisms and Population-Level Effects

Most microorganisms regulate their cell size. In this article, we review some of the mathematical formulations of the problem of cell size regulation. We focus on coarse-grained stochastic models and the statistics that they generate. We review the biologically relevant insights obtained from these models. We then describe cell cycle regulation and its molecular implementations, protein number regulation, and population growth, all in relation to size regulation. Finally, we discuss several future directions for developing understanding beyond phenomenological models of cell size regulation.

Cell size control and gene expression homeostasis in single-cells

Growth of a cell and its subsequent division into daughters is a fundamental aspect of all cellular living systems. During these processes, how do individual cells correct size aberrations so that they do not grow abnormally large or small? How do cells ensure that the concentration of essential gene products are maintained at desired levels, in spite of dynamic/stochastic changes in cell size during growth and division? Both these questions have fascinated researchers for over a century. We review how advances in singe-cell technologies and measurements are providing unique insights into these questions across organisms from prokaryotes to human cells. More specifically, diverse strategies based on timing of cell-cycle events, regulating growth, and number of daughters are employed to maintain cell size homeostasis. Interestingly, size homeostasis often results in size optimality - proliferation of individual cells in a population is maximized at an optimal cell size. We further discuss how size-dependent expression or gene-replication timing can buffer concentration of a gene product from cell-to-cell size variations within a population. Finally, we speculate on an intriguing hypothesis that specific size control strategies may have evolved as a consequence of gene-product concentration homeostasis.

Learning from Noise: How Observing Stochasticity May Aid Microbiology

For many decades, the wedding of quantitative data with mathematical modeling has been fruitful, leading to important biological insights. Here, we review some of the ongoing efforts to gain insights into problems in microbiology - and, in particular, cell-cycle progression and its regulation - through observation and quantitative analysis of the natural fluctuations in the system. We first illustrate this idea by reviewing a classic example in microbiology - the Luria-Delbrück experiment - and discussing how, in that case, useful information was obtained by looking beyond the mean outcome of the experiment, but instead paying attention to the variability between replicates of the experiment. We then switch gears to the contemporary problem of cell cycle progression and discuss in more detail how insights into cell size regulation and, when relevant, coupling between the cell cycle and the circadian clock, can be gained by studying the natural fluctuations in the system and their statistical properties. We end with a more general discussion of how (in this context) the correct level of phenomenological model should be chosen, as well as some of the pitfalls associated with this type of analysis. Throughout this review the emphasis is not on providing details of the experimental setups or technical details of the models used, but rather, in fleshing out the conceptual structure of this particular approach to the problem. For this reason, we choose to illustrate the framework on a rather broad range of problems, and on organisms from all domains of life, to emphasize the commonality of the ideas and analysis used (as well as their differences).

Archaeal cells share common size control with bacteria despite noisier growth and division

In nature, microorganisms exhibit different volumes spanning six orders of magnitude ¹ . Despite their capability to create different sizes, a clonal population in a given environment maintains a uniform size across individual cells. Recent studies in eukaryotic and bacterial organisms showed that this homogeneity in cell size can be accomplished by growing a constant size between two cell cycle events (that is, the adder model ^2-6 ). Demonstration of the adder model led to the hypothesis that this phenomenon is a consequence of convergent evolution. Given that archaeal cells share characteristics with both bacteria and eukaryotes, we investigated whether and how archaeal cells exhibit control over cell size. To this end, we developed a soft-lithography method of growing the archaeal cells to enable quantitative time-lapse imaging and single-cell analysis, which would be useful for other microorganisms. Using this method, we demonstrated that Halobacterium salinarum, a hypersaline-adapted archaeal organism, grows exponentially at the single-cell level and maintains a narrow-size distribution by adding a constant length between cell division events. Interestingly, the archaeal cells exhibited greater variability in cell division placement and exponential growth rate across individual cells in a population relative to those observed in Escherichia coli ^6-9 . Here, we present a theoretical framework that explains how these larger fluctuations in archaeal cell cycle events contribute to cell size variability and control.

The Effects of Stochasticity at the Single-Cell Level and Cell Size Control on the Population Growth

Establishing a quantitative connection between the population growth rate and the generation times of single cells is a prerequisite for understanding evolutionary dynamics of microbes. However, existing theories fail to account for the experimentally observed correlations between mother-daughter generation times that are unavoidable when cell size is controlled for, which is essentially always the case. Here, we study population-level growth in the presence of cell size control and corroborate our theory using experimental measurements of single-cell growth rates. We derive a closed formula for the population growth rate and demonstrate that it only depends on the single-cell growth rate variability, not other sources of stochasticity. Our work provides an evolutionary rationale for the narrow growth rate distributions often observed in nature: when single-cell growth rates are less variable but have a fixed mean, the population will exhibit an enhanced population growth rate as long as the correlations between the mother and daughter cells' growth rates are not too strong.

Analysis of Noise Mechanisms in Cell-Size Control

At the single-cell level, noise arises from multiple sources, such as inherent stochasticity of biomolecular processes, random partitioning of resources at division, and fluctuations in cellular growth rates. How these diverse noise mechanisms combine to drive variations in cell size within an isoclonal population is not well understood. Here, we investigate the contributions of different noise sources in well-known paradigms of cell-size control, such as adder (division occurs after adding a fixed size from birth), sizer (division occurs after reaching a size threshold), and timer (division occurs after a fixed time from birth). Analysis reveals that variation in cell size is most sensitive to errors in partitioning of volume among daughter cells, and not surprisingly, this process is well regulated among microbes. Moreover, depending on the dominant noise mechanism, different size-control strategies (or a combination of them) provide efficient buffering of size variations. We further explore mixer models of size control, where a timer phase precedes/follows an adder, as has been proposed in Caulobacter crescentus. Although mixing a timer and an adder can sometimes attenuate size variations, it invariably leads to higher-order moments growing unboundedly over time. This results in a power-law distribution for the cell size, with an exponent that depends inversely on the noise in the timer phase. Consistent with theory, we find evidence of power-law statistics in the tail of C. crescentus cell-size distribution, although there is a discrepancy between the observed power-law exponent and that predicted from the noise parameters. The discrepancy, however, is removed after data reveal that the size added by individual newborns in the adder phase itself exhibits power-law statistics. Taken together, this study provides key insights into the role of noise mechanisms in size homeostasis, and suggests an inextricable link between timer-based models of size control and heavy-tailed cell-size distributions.

Evidence that Listeria innocua modulates its membrane's stored curvature elastic stress, but not fluidity, through the cell cycle

This paper reports that the abundances of endogenous cardiolipin and phosphatidylethanolamine halve during elongation of the Gram-positive bacterium Listeria innocua. The lyotropic phase behaviour of model lipid systems that describe these modulations in lipid composition indicate that the average stored curvature elastic stress of the membrane is reduced on elongation of the cell, while the fluidity appears to be maintained. These findings suggest that phospholipid metabolism is linked to the cell cycle and that changes in membrane composition can facilitate passage to the succeding stage of the cell cycle. This therefore suggests a means by which bacteria can manage the physical properties of their membranes through the cell cycle.

Cell size control - a mechanism for maintaining fitness and function

The maintenance of cell size homeostasis has been studied for years in different cellular systems. With the focus on 'what regulates cell size', the question 'why cell size needs to be maintained' has been largely overlooked. Recent evidence indicates that animal cells exhibit nonlinear cell size dependent growth rates and mitochondrial metabolism, which are maximal in intermediate sized cells within each cell population. Increases in intracellular distances and changes in the relative cell surface area impose biophysical limitations on cells, which can explain why growth and metabolic rates are maximal in a specific cell size range. Consistently, aberrant increases in cell size, for example through polyploidy, are typically disadvantageous to cellular metabolism, fitness and functionality. Accordingly, cellular hypertrophy can potentially predispose to or worsen metabolic diseases. We propose that cell size control may have emerged as a guardian of cellular fitness and metabolic activity.

Mycobacteria Modify Their Cell Size Control under Sub-Optimal Carbon Sources

The decision to divide is the most important one that any cell must make. Recent single cell studies suggest that most bacteria follow an “adder” model of cell size control, incorporating a fixed amount of cell wall material before dividing. Mycobacteria, including the causative agent of tuberculosis Mycobacterium tuberculosis, are known to divide asymmetrically resulting in heterogeneity in growth rate, doubling time, and other growth characteristics in daughter cells. The interplay between asymmetric cell division and adder size control has not been extensively investigated. Moreover, the impact of changes in the environment on growth rate and cell size control have not been addressed for mycobacteria. Here, we utilize time-lapse microscopy coupled with microfluidics to track live Mycobacterium smegmatis cells as they grow and divide over multiple generations, under a variety of growth conditions. We demonstrate that, under optimal conditions, M. smegmatis cells robustly follow the adder principle, with constant added length per generation independent of birth size, growth rate, and inherited pole age. However, the nature of the carbon source induces deviations from the adder model in a manner that is dependent on pole age. Understanding how mycobacteria maintain cell size homoeostasis may provide crucial targets for the development of drugs for the treatment of tuberculosis, which remains a leading cause of global mortality.

Long-term microfluidic tracking of coccoid cyanobacterial cells reveals robust control of division timing

Background

Cyanobacteria are important agents in global carbon and nitrogen cycling and hold great promise for biotechnological applications. Model organisms such as Synechocystis sp. and Synechococcus sp. have advanced our understanding of photosynthetic capacity and circadian behavior, mostly using population-level measurements in which the behavior of individuals cannot be monitored. Synechocystis sp. cells are small and divide slowly, requiring long-term experiments to track single cells. Thus, the cumulative effects of drift over long periods can cause difficulties in monitoring and quantifying cell growth and division dynamics.

Results

To overcome this challenge, we enhanced a microfluidic cell-culture device and developed an image analysis pipeline for robust lineage reconstruction. This allowed simultaneous tracking of many cells over multiple generations, and revealed that cells expand exponentially throughout their cell cycle. Generation times were highly correlated for sister cells, but not between mother and daughter cells. Relationships between birth size, division size, and generation time indicated that cell-size control was inconsistent with the “sizer” rule, where division timing is based on cell size, or the “timer” rule, where division occurs after a fixed time interval. Instead, single cell growth statistics were most consistent with the “adder” rule, in which division occurs after a constant increment in cell volume. Cells exposed to light-dark cycles exhibited growth and division only during the light period; dark phases pause but do not disrupt cell-cycle control.

Conclusions

Our analyses revealed that the “adder” model can explain both the growth-related statistics of single Synechocystis cells and the correlation between sister cell generation times. We also observed rapid phenotypic response to light-dark transitions at the single cell level, highlighting the critical role of light in cyanobacterial cell-cycle control. Our findings suggest that by monitoring the growth kinetics of individual cells we can build testable models of circadian control of the cell cycle in cyanobacteria.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-016-0344-4) contains supplementary material, which is available to authorized users.

Effects of cell-cycle-dependent expression on random fluctuations in protein levels

Expression of many genes varies as a cell transitions through different cell-cycle stages. How coupling between stochastic expression and cell cycle impacts cell-to-cell variability (noise) in the level of protein is not well understood. We analyse a model where a stable protein is synthesized in random bursts, and the frequency with which bursts occur varies within the cell cycle. Formulae quantifying the extent of fluctuations in the protein copy number are derived and decomposed into components arising from the cell cycle and stochastic processes. The latter stochastic component represents contributions from bursty expression and errors incurred during partitioning of molecules between daughter cells. These formulae reveal an interesting trade-off: cell-cycle dependencies that amplify the noise contribution from bursty expression also attenuate the contribution from partitioning errors. We investigate the existence of optimum strategies for coupling expression to the cell cycle that minimize the stochastic component. Intriguingly, results show that a zero production rate throughout the cell cycle, with expression only occurring just before cell division, minimizes noise from bursty expression for a fixed mean protein level. By contrast, the optimal strategy in the case of partitioning errors is to make the protein just after cell division. We provide examples of regulatory proteins that are expressed only towards the end of the cell cycle, and argue that such strategies enhance robustness of cell-cycle decisions to the intrinsic stochasticity of gene expression.

Intercellular Variability in Protein Levels from Stochastic Expression and Noisy Cell Cycle Processes

Inside individual cells, expression of genes is inherently stochastic and manifests as cell-to-cell variability or noise in protein copy numbers. Since proteins half-lives can be comparable to the cell-cycle length, randomness in cell-division times generates additional intercellular variability in protein levels. Moreover, as many mRNA/protein species are expressed at low-copy numbers, errors incurred in partitioning of molecules between two daughter cells are significant. We derive analytical formulas for the total noise in protein levels when the cell-cycle duration follows a general class of probability distributions. Using a novel hybrid approach the total noise is decomposed into components arising from i) stochastic expression; ii) partitioning errors at the time of cell division and iii) random cell-division events. These formulas reveal that random cell-division times not only generate additional extrinsic noise, but also critically affect the mean protein copy numbers and intrinsic noise components. Counter intuitively, in some parameter regimes, noise in protein levels can decrease as cell-division times become more stochastic. Computations are extended to consider genome duplication, where transcription rate is increased at a random point in the cell cycle. We systematically investigate how the timing of genome duplication influences different protein noise components. Intriguingly, results show that noise contribution from stochastic expression is minimized at an optimal genome-duplication time. Our theoretical results motivate new experimental methods for decomposing protein noise levels from synchronized and asynchronized single-cell expression data. Characterizing the contributions of individual noise mechanisms will lead to precise estimates of gene expression parameters and techniques for altering stochasticity to change phenotype of individual cells.

Inside individual cells, gene products often occur at low molecular counts and are subject to considerable stochastic fluctuations (noise) in copy numbers over time. An important consequence of noisy expression is that the level of a protein can vary considerably even among genetically identical cells exposed to the same environment. Such non-genetic phenotypic heterogeneity is physiologically relevant and critically influences diverse cellular processes. In addition to noise sources inherent in gene product synthesis, recent experimental studies have uncovered additional noise mechanisms that critically effect expression. For example, the time within the cell cycle when a gene duplicates, and the time taken to complete cell cycle are governed by random processes. The key contribution of this work is development of novel mathematical results quantifying how cell cycle-related noise sources combine with stochastic expression to drive intercellular variability in protein molecular counts. Derived formulas lead to many counterintuitive results, such as increasing randomness in the timing of cell division can lower noise in the level of a protein. Finally, these results inform experimental strategies to systematically dissect the contributions of different noise sources in the expression of a gene of interest.