[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Advertisement

  • Loading metrics

Open Access

Peer-reviewed

Research Article

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

  • Mohammad Soltani,

    Affiliation Electrical and Computer Engineering Department, University of Delaware, Newark, Delaware, United States of America

    ⨯
  • Cesar A. Vargas-Garcia,

    Affiliation Electrical and Computer Engineering Department, University of Delaware, Newark, Delaware, United States of America

    ⨯
  • Duarte Antunes,

    Affiliation Mechanical Engineering Department, Eindhoven University of Technology, Eindhoven, Netherlands

    ⨯
  • Abhyudai Singh

    absingh@udel.udu

    Affiliations Electrical and Computer Engineering Department, University of Delaware, Newark, Delaware, United States of America, Biomedical Engineering Department, University of Delaware, Newark, Delaware, United States of America, Mathematical Sciences Department, University of Delaware, Newark, Delaware, United States of America, Center for Bioinformatics and Computational Biology, University of Delaware, Newark, Delaware, United States of America

    ⨯

Abstract

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.

Author Summary

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.

Citation: Soltani M, Vargas-Garcia CA, Antunes D, Singh A (2016) Intercellular Variability in Protein Levels from Stochastic Expression and Noisy Cell Cycle Processes. PLoS Comput Biol 12(8): e1004972. https://doi.org/10.1371/journal.pcbi.1004972

Editor: Suckjoon Jun, UCSD, UNITED STATES

Received: September 8, 2015; Accepted: July 29, 2016; Published: August 18, 2016

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

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: This work is supported by the National Science Foundation Grant DMS-1312926. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

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

Introduction

The level of a protein can deviate considerably from cell-to-cell, in spite of the fact that cells are genetically-identical and are in the same extracellular environment [1–3]. This intercellular variation or noise in protein counts has been implicated in diverse processes such as corrupting functioning of gene networks [4–6], driving probabilistic cell-fate decisions [7–12], buffering cell populations from hostile changes in the environment [13–16], and causing clonal cells to respond differently to the same stimulus [17–19]. An important source of noise driving random fluctuations in protein levels is stochastic gene expression due to the inherent probabilistic nature of biochemical processes [20–23]. Recent experimental studies have uncovered additional noise sources that affect protein copy numbers. For example, the time take to complete cell cycle (i.e., time between two successive cell-division events) has been observed to be stochastic across organisms [24–32]. Moreover, given that many proteins/mRNAs are present inside cells at low-copy numbers, errors incurred in partitioning of molecules between two daughter cells are significant [33–35]. Finally, the time at which a particular gene of interest is duplicated can also vary between cells [36, 37]. We investigate how such noise sources in the cell-cycle process combine with stochastic gene expression to generate intercellular variability in protein copy numbers (Fig 1).

thumbnail
Download:
Fig 1. Sample trajectory of the protein level in a single cell with different sources of noise.

Stochastically expressed proteins accumulate within the cell at a certain rate. At a random point in the cell cycle, gene duplication results in an increase in production rate. Stochastic cell-division events lead to random partitioning of protein molecules between two daughter cells with each cell receiving, on average, half the number of proteins in the mother cell just before division. The steady-state protein copy number distribution obtained from a large number of trajectories is shown on the right. The total noise in the protein level, as measured by the squared coefficient of variation (CV2) can be broken into contributions from individual noise mechanisms.

https://doi.org/10.1371/journal.pcbi.1004972.g001

Prior studies that quantify the effects of cell division on the protein noise level have been restricted to specific cases. For example, noise computations have been done in stochastic gene expression models, where cell divisions occur at deterministic time intervals [33, 38, 39]. Recently, we have analyzed a deterministic model of gene expression with random cell-division events [40]. Building up on this work, we formulate a mathematical model that couples stochastic expression of a stable protein with random cell-division events that follow a general class of probability distributions. Moreover, at the time of cell division, proteins are randomly partitioned between two daughter cells based on a framework that allows the partitioning errors to be higher or lower than as predicted by binomial partitioning. For this class of models, we derive an exact analytical formula for the protein noise level as quantified by the steady-state squared Coefficient of Variation (CV2). This formula is further decomposed into individual components representing contributions from different noise sources. A systematic investigation of this formula leads to novel insights, such as identification of regimes where increasing randomness in the timing of cell-division events decreases the protein noise level.

Next, we extend the above model to include genome-duplication events that increase the gene’s transcription rate [36, 41]. To our knowledge, this is the first study integrating randomness in the genome-duplication process with stochastic gene expression. An exact formula for the protein noise level is derived for this extended model and used to investigate how the timing of duplication affects different noise components. Counter intuitively, results show that doubling of the transcription rate within the cell cycle can lead to smaller fluctuations in protein levels as compared to a constant transcription rate through out the cell cycle. Finally, we discuss how formulas obtained in this study can be used to infer parameters and characterize the gene expression process from single-cell studies.

Methods

Coupling gene expression to cell division

We consider the standard model of stochastic gene expression [42, 43], where mRNAs are transcribed at exponentially distributed time intervals from a constitutive gene with rate kx. For the time being, we exclude genome duplication and the transcription rate is fixed throughout the cell cycle. Assuming short-lived mRNAs, each transcription event results in a burst of proteins [43–45]. The corresponding jump in protein levels is shown as (1) where x(t) is the protein population count in the mother cell at time t, B is a random burst size drawn from a positively-valued distribution and represents the number of protein molecules synthesized in a single-mRNA lifetime. Motivated by observations in E. coli and mammalian cells, where many proteins have half-lives considerably longer than the cell-doubling time, we assume a stable protein with no active degradation [46–48]. Thus, proteins accumulate within the cell till the time of cell division, at which point they are randomly partitioned between two daughter cells.

Let cell division events occur at times ts, s ∈ {1, 2, …}. The cell-cycle time (2) follows an arbitrary positively-valued probability distribution with the following mean and squared coefficient of variation (CV2) (3) where 〈.〉 denotes expected value through out this paper. The random change in x(t) during cell division is given by (4) where x(ts) denotes the protein levels in the mother cell just before division and x+(ts) denotes the protein levels in one of the daughter cells just after division. Conditioned on x(ts), x+(ts) is assumed to have the following statistics (5) The first equation implies symmetric partitioning, i.e., on average each of the daughter cells inherits half the number protein molecules just before division. The second equation in Eq (5) describes the variance of x+(ts) and quantifies the error in partitioning of molecules through the non-negative parameter α. For example, α = 0 represents deterministic partitioning where x+(ts) = x(ts)/2 with probability equal to one. A more realistic model for partitioning is each molecule having an equal probability of being in the each daughter cell [49–51]. This results in a binomial distribution for x+(ts) (6) and corresponds to α = 1 in Eq (5). Interestingly, recent studies have shown that partitioning of proteins that form clusters or multimers can result in α > 1 in Eq (5), i.e., partitioning errors are much higher than as predicted by the binomial distribution [33, 39]. In contrast, if molecules push each other to opposite poles of the cell, then the partitioning errors will be smaller than as predicted by Eq (6) and α < 1.

The model with all the different noise mechanisms (stochastic expression; random cell-division events and partitioning errors) is illustrated in Fig 2A and referred to as the full model. We also introduce two additional hybrid models [52, 53], where protein production and partitioning are considered in their deterministic limit (Fig 2B and 2C). Note that unlike the full model, where x(t) takes non-negative integer values, x(t) is continuous in the hybrid models. We will use these hybrid models for decomposing the protein noise level obtained from the full model into individual components representing contributions from different noise sources.

thumbnail
Download:
Fig 2. Stochastic models of gene expression with cell division.

Arrows denote stochastic events that change the protein level by discrete jumps as shown in Eqs (1) and (4). The differential equation within the circle represents the time evolution of x(t) in between events. A) Model with all the different sources of noise: proteins are expressed in stochastic bursts, cell division occurs at random times, and molecules are partitioned between the two daughter cells based on Eq (5). The trivial dynamics signifies that the protein level is constant in-between stochastic events. B) Hybrid model where randomness in cell-division events is the only source of noise. Protein production is modeled deterministic through a differential equation and partitioning errors are absent, i.e., α = 0 in Eq (5). C) Hybrid model where noise comes from both cell-division events and partitioning errors. Protein production is considered to be deterministic as in Fig 2B. Since x(t) is continuous here, x+(ts) has a positively-valued continuous distribution with same mean and variance as in Eq (5)

https://doi.org/10.1371/journal.pcbi.1004972.g002

Modeling the cell-cycle time

In order to quantify the steady-state protein mean and noise, we need to define the stochastic process that governs the timing of cell division. Variations in the duration of cell cycle can result from a variety of factors, such as cell physiology, growth rate, cell size and expression of genes that affect cell-cycle time such as FtsZ [24–32]. Given these complexities, we take a phenomenological approach to modeling cell-cycle time, and assume it to be an independent and identically distributed random variable that is drawn from a mixture of Erlang distributions (also known as phase-type distribution). The motivation for choosing this distribution is two fold:

  1. Mixture of Erlang distributions can be represented via a continuous-time Markov chain, allowing mathematical tractability in terms of deriving/analyzing time evolution of moments.
  2. This class of distribution is fairly general, in the sense that, any positively-valued distribution with CV ≤ 1 can be modeled via a mixture of Erlang distributions [54].

Consider a mixture of n Erlang distributions with mixing probabilities pi, i = {1, …, n}. Recall that an Erlang distribution of order i is the distribution of the sum of i independent and identical exponential random variables. The cell-cycle time is assumed to have an Erlang distribution of order i with probability pi and can be represented by a continuous-time Markov chain with states Gij, j = {1, …, i}, i = {1, …, n} (Fig 3). Let Bernoulli random variables gij = 1 if the system resides in state Gij and 0 otherwise. The probability of transition Gij → Gi(j+1) in the next infinitesimal time interval [t, t + dt) is given by ikgij dt, implying that the time spent in each state Gij is exponentially distributed with mean 1/ik. To summarize, at the start of cell cycle, a state Gi1, i = {1, …, n} is chosen with probability pi and cell division occurs after transitioning through i exponentially distributed steps. Based on this formulation, the probability of a cell-division event occurring and a new cell cycle obtained from an Erlang distribution of size i starting in the next time interval [t, t + dt) is given by , and whenever the event occurs, the protein level changes as per Eq (4). Finally, the mean, the squared coefficient of variation, and the skewness of the cell-cycle time in terms of the Markov chain parameters are given by (7) [55], where 〈T3〉 is the third-order moment of the cell-cycle time. An important property of this class of distributions is that increasing also makes the distribution highly skewed, because from Eq (7) both the CV and skewness are linear combinations of pi, albeit with different linear coefficients that decrease with i. Considering that , the only way to increase is by increasing smaller-index components and decreasing larger-index components of the distribution (i.e. increasing pi and decreasing pj, where i < j). Since higher values of i are more penalized in the skewness equation, this would correspond to making the distribution more positively skewed. Hence high values of also means high values of skewness, thus occurrences of longer cell cycles are more probable. As we will shortly see, this property leads to mean protein levels being dependent on .

thumbnail
Download:
Fig 3. A continuous-time Markov chain model for the cell-cycle time.

Left: The cell-cycle time is assumed to follow a mixture of Erlang distributions. At the start of cell cycle, a state Gi1, i = {1, …, n} is chosen with probability pi. The cell cycle transitions through states Gij, j = {1, …, i} residing for an exponentially distributed time with mean 1/ik in each state. Cell division occurs after exit from Gii and the above process is repeated.

https://doi.org/10.1371/journal.pcbi.1004972.g003

Results

Computing the average number of protein molecules

All the models shown in Fig 2 are identical in terms of finding 〈x(t)〉 and in principle any one of them could have been used. We choose to analyze the full model illustrated in Fig 2A. Time evolution of the statistical moments of x(t) can be obtained from the Kolmogorov forward equations corresponding to the full model in Fig 2A combined with the cell-division process described in Fig 3. We refer the reader to [52, 56, 57] for an introduction to moment dynamics for stochastic and hybrid systems. Analysis in Appendix A in S1 Text shows (8) Note that the time-derivative of the mean protein level (first-order moment) is unclosed, in the sense that, it depends on the second-order moment 〈xgij〉. Typically, approximate closure methods are used to solve moments in such cases [52, 57–62]. However, the fact that gij is binary can be exploited to automatically close moment dynamics. In particular, since gij ∈ {0, 1} (9) for any non-negative integer m. Moreover, as only a single state gij can be 1 at any time (10) Using Eqs (9) and (10), the time evolutions of 〈gij〉 and 〈xgij〉 are obtained as (11) (12) (13) (14) and only depend on 〈gij〉 and 〈xgij〉 (see Appendix A in S1 Text). Thus, Eqs (8) and (11)–(14) constitute a closed system of linear differential equations from which moments can be computed exactly.

To obtain an analytical formula for the average number of proteins, we start by performing a steady-state analysis of Eq (8) that yields (15) where denotes the expected value in the limit t → ∞. Using Eq (15), is determined from Eq (13), and then all moments are obtained recursively by performing a steady-state analysis of Eq (14) for j = {2, …, i}. This analysis results in (16) Using Eqs (7), (16) and the fact that we obtain the following expression for the mean protein level (17) It is important to point that Eq (17) holds irrespective of the complexity, i.e., the number of states Gij used in the phase-type distribution to approximate the cell-cycle time distribution. As expected, increases linearly with the average cell-cycle time duration 〈T〉 with longer cell cycles resulting in more accumulation of proteins. Consistent with previous findings, Eq (17) shows that the mean protein level is also affected by the randomness in the cell-cycle times [40, 63]. For example, reduces by 25% as T changes from being exponentially distributed to periodic for fixed 〈T〉. Next, we determine the noise in protein copy numbers, as quantified by the squared coefficient of variation.

Computing the protein noise level

Recall that the full model introduced in Fig 2A has three distinct noise mechanisms. Our strategy for computing the protein noise level is to first analyze the model with a single noise source, and then consider models with two and three sources. As shown below, this approach provides a systematic dissection of the protein noise level into components representing contributions from different mechanisms.

Contribution from randomness in cell-cycle times.

We begin with the model shown in Fig 2B, where noise comes from a single source—random cell-division events. For this model, the time evolution of the second-order moment of the protein copy number is obtained as (18) and depends on third-order moments 〈x2 gjj〉 (see Appendix B in S1 Text). Using the approach introduced earlier for obtaining the mean protein level, we close moment equations by writing the time evolution of moments 〈x2 gij〉. Using Eqs (9) and (10) (19) (20) Note that the moment dynamics for 〈x〉 and 〈xgij〉 obtained in the previous section (Eqs (8), (13) and (14)) are identical for all the models in Fig 2, irrespective of whether the noise mechanism is modeled deterministically or stochastically. Eqs (8), (11)–(14) and (18)–(20) represent a closed set of linear differential equations and their steady-state analysis yields (21) From Eq (21) (22) (23) Using Eq (22) and the mean protein count quantified in Eq (17), we obtain the following squared coefficient of variation (24) where represents the noise contribution from random cell-division events. Since cell division is a global event that affects expression of all genes, this noise contribution can also be referred to as extrinsic noise [49, 64–67]. In reality, there would be other sources of extrinsic noise, such as, fluctuations in the gene-expression machinery that we have ignored in this analysis.

Note that as T approaches a delta distribution, i.e., cell divisions occur at fixed time intervals. We discuss simplifications of Eq (24) in various limits. For example, if the time taken to complete cell cycle is lognormally distributed, then (25) and extrinsic noise monotonically increases with . If fluctuations in T around 〈T〉 are small, then using Taylor series (26) Substituting Eq (26) in Eq (24) and ignoring and higher order terms yields (27) where the first term is the extrinsic noise for and the second term is the additional noise due to random cell-division events.

Contribution from partitioning errors.

Next, we consider the model illustrated in Fig 2C with both random cell-division events and partitioning of protein between the two daughter cells. Thus, the protein noise level here represents the contribution from both these sources. Analysis in Appendix C in S1 Text shows that the time evolution of 〈x2〉 and 〈x2 gij〉 are given by (28) (29) (30) Note that Eqs (28) and (29) are slightly different from their counterparts obtained in the previous section (Eqs (18) and (19)) with additional terms that depend on α, where α quantifies the degree of partitioning error as defined in Eq (5). As expected, Eqs (28) and (29) reduce to Eqs (18) and (19) when α = 0 (i.e., deterministic partitioning). Computing by performing a steady-state analysis of Eqs (28)–(30) and using a similar approach as in Eq (22) we obtain (31) Finding CV2 of the protein level and subtracting the extrinsic noise found in Eq (24) yields (32) where represents the contribution of partitioning errors to the protein noise level. Intriguingly, while increases with α, it decrease with . Thus, as cell-division times become more random for a fixed 〈T〉 and , the noise contribution from partitioning errors decrease. It turns out that this dependence of on CVT is a direct result of the second equation in Eq (5), where stochasticity in the partitioning process increases linearly with x(ts), the number of protein molecules just before division. Based on Eq (17), one needs to reduce kx or 〈B〉 to maintain a fixed for increasing randomness in cell-division times. Since the average number of protein molecules just before division is 2kx〈B〉〈T〉 (see Appendix D in S1 Text), a reduction in kx or 〈B〉 results in a lower number of protein molecules before division, and hence, lesser noise from partitioning as per Eq (5) and a smaller . This reasoning is supported by the fact that if we redefine the noise in the partitioning process to make it independent of x(ts), i.e. modify Eq (5) as (33) then the noise contribution from partitioning errors is given by (34) and the dependency of on CVT disappears (Appendix D in S1 Text).

Contribution from stochastic expression.

Finally, we consider the full model in Fig 2A with all the three different noise sources. For this model, moment dynamics is obtained as (see Appendix E in S1 Text) (35) (36) (37)

Compared to Eqs (28)–(30), (35)–(37) have additional terms of the form kx〈B2〉, where 〈B2〉 is the second-order moment of the protein burst size in Eq (1). Performing an identical analysis as before we obtain (38) which yields the following total protein noise level (39) that can be decomposed into three terms. The first term represents the contribution from random cell-division events and is given by Eq (24). The second term is the contribution from partitioning errors determined in the previous section (partitioning noise), and the final term is the additional noise representing the contribution from stochastic expression (production noise). A common approach to study gene expression noise is to decompose it into intrinsic and extrinsic components. These components are obtained experimentally using the dual-color assay that measures the correlation in the expression of two identical copies of the gene [49]. As per this definition, represents the extrinsic noise as random cell-division events are common to all genes and makes expression levels more correlated in individual cells. In contrast, the contributions from noisy production and partitioning represent the intrinsic noise as they are specific to an individual gene and make expression levels less correlated.

An interesting observation from Eq (39) is that has opposite effects on and (for fixed mean protein level). While monotonically decreases with increasing , increases with . Thus, if is small and α is large, then the noise contributed from partitioning dominates the total noise, and making cell-cycle duration more random will reduce the total noise. However, since both and are monotonically increasing functions of , the total noise will begin to increase with once these noise sources become dominant. It turns out that in certain cases the intrinsic noise becomes invariant of . For example, when B = 1 with probability one, i.e., proteins are synthesized one at a time at exponentially distributed time intervals and α = 1 (binomial partitioning) (40) In this limit the intrinsic noise is always 1/Mean irrespective of the cell-cycle time distribution T [33]. Note that the average number of proteins itself depends on T as shown in Eq (17). Another important limit is , in which case Eq (39) reduces to (41) and is similar to the result obtained in [38] for deterministic cell-division times and binomial partitioning.

Fig 4 shows how different protein noise components change as a function of the mean protein level as the gene’s transcription rate kx is modulated. The extrinsic noise is primarily determined by the distribution of the cell-cycle time and is completely independent of the mean. In contrast, both and scale inversely with the mean, albeit with different scaling factors (Fig 4). This observation is particularly important since many single-cell studies in E. coli, yeast and mammalian cells have found the protein noise levels to scale inversely with the mean across different genes [68–71]. Based on this scaling it is often assumed that the observed cell-to-cell variability in protein copy numbers is a result of stochastic expression. However, as our results show, noise generated thorough partitioning errors is also consistent with these experimental observations and it may be impossible to distinguish between these two noise mechanisms based on protein CV2 versus mean plots unless α is known.

thumbnail
Download:
Fig 4. Scaling of noise as a function of the mean protein level for different mechanisms.

The contribution of random cell-division events to the noise in protein copy numbers (extrinsic noise) is invariant of the mean. In contrast, contributions from partitioning errors at the time of cell division (partitioning noise) and stochastic expression (production noise) scale inversely with the mean. The scaling factors are shown as a function of the protein random burst size B, noise in cell-cycle time () and magnitude of partitioning errors quantified by α (see Eq (5)). With increasing mean level the total noise first decreases and then reaches a baseline that corresponds to extrinsic noise. For this plot, B is assumed to be geometrically-distributed with mean 〈B〉 = 1.5, and α = 1 (i.e., binomial partitioning).

https://doi.org/10.1371/journal.pcbi.1004972.g004

Quantifying the effects of gene duplication on protein noise

The full model introduced in Fig 2 assumes that the transcription rate (i.e., the protein burst arrival rate) is constant throughout the cell cycle. This model is now extended to incorporate gene duplication during cell cycle, which increases the burst arrival (transcription) rate by f times (Fig 5). Note that due to gene dosage compensation, doubling the genome does not always correspond to f = 2 [72–74]. If f > 1, then accumulation of proteins will be bilinear as illustrated in Fig 1. As before, we again take a phenomenological approach to model the timing of gene duplication. The cell-cycle time T is divided into two intervals: time from the start of cell cycle to gene duplication (T1), and time from gene duplication to cell division (T2). T1 and T2 are independent random variables, each drawn from a mixture of Erlang distributions (see Fig. B in S1 Text). The mean cell-cycle duration and its noise can be expressed as (42) where denotes the squared coefficient of variation of the random variable X. An important variable in this formulation is β, which represents the average time of gene duplication normalized by the mean cell-cycle time. Thus, β values close to 0 (1) imply that the gene is duplicated early (late) in the cell-cycle process. Moreover, the noise in the gene-duplication time is controlled via .

thumbnail
Download:
Fig 5. Model illustrating stochastic expression together with random gene-duplication and cell-division events.

At the start of cell cycle, protein production occurs in stochastic bursts with rate kx. Genome duplication occurs at a random point T1 within the cell cycle and increases the burst arrival rate to fkx (f > 1). Cell division occurs after time T2 from genome duplication, at which point the burst arrival rate reverts back to kx and proteins are randomly partitioned between cells based on Eq (4).

https://doi.org/10.1371/journal.pcbi.1004972.g005

We refer the reader to Appendix F in S1 Text for a detailed analysis of the model in Fig 5 and only present the main results on the protein mean and noise levels. The steady-state mean protein count is given by (43) and decreases with β, i.e., a gene that duplicates early has on average, more number of proteins. When β = 1, then the transcription rate is kx throughout the cell cycle and we recover the mean protein level obtained in Eq (17). Similarly, when β = 0 the transcription rate is fkx and we obtain f times of the amount as in Eq (17). As per our earlier observation, more randomness in the timing of genome duplication and cell division (i.e., higher and values) increases .

Our analysis shows that the total protein noise level can be decomposed into three components (44) where is the extrinsic noise from random genome-duplication/cell-division events, and the sum of the contributions from partitioning errors () and stochastic expression () is the intrinsic noise. We refer the reader to Appendix F in S1 Text for noise formulas for any f, and only present formulas for f = 2 here. In this case, the intrinsic noise is obtained as (45)

Note that for β = 0 and 1, we recover the intrinsic noise level in Eq (39) from Eq (45). Interestingly, for B = 1 with probability 1 and α = 1, the intrinsic noise is always 1/Mean irrespective of the values chosen for , and β. For high precision in the timing of cell-cycle events (CVT1 → 0, CVT2 → 0) (46) where mean protein level is given by (47) We investigate how different noise components in Eq (46) vary with β as the mean protein level is held fixed by changing kx. Fig 6 shows that follows a U-shaped profile with the optima occurring at and the corresponding minimum value being ≈ 5% lower than its value at β = 0. An implication of this result is that if stochastic expression is the dominant noise source, then gene duplication can result in slightly lower protein noise levels. In contrast to , has a maxima at which is ≈ 6% higher than its value at β = 0 (Fig 6). Analysis in Appendix F5 in S1 Text reveals that and follow the same qualitative shapes as in Fig 6 for any and . Interestingly, when , the maximum and minimum values of and always occur at albeit with different optimal values than Fig 6 (see Fig. C in S1 Text). For example, if (i.e., exponentially distributed T1 and T2), then the maximum value of is 20% higher and the minimum value of is 10% lower than their respective value for β = 0. Given that the effect of changing β on and is small and antagonistic, the overall affect of genome duplication on intrinsic noise may be minimal and hard to detect experimentally.

thumbnail
Download:
Fig 6. Contributions from different noise sources as a function of the timing of genome duplication for .

Different noise components in Eq (46) are plotted as a function of β, which represents the fraction of time within the cell cycle at which gene duplication occurs. The mean protein level is held constant by simultaneously changing the transcription rate kx. Noise levels are normalized by their respective value at β = 0. The noise contribution from partitioning errors is maximized at β ≈ 0.6. In contrast, the contribution from stochastic expression is minimum at β ≈ 0.6. The extrinsic noise contribution from random gene-duplication and cell-division events is maximum at β ≈ 0.2 and minimum at β ≈ 0.8. For this plot, the mean of the protein is 170 molecules per cell; and the bursts are geometrically distributed with 〈B〉 = 10.

https://doi.org/10.1371/journal.pcbi.1004972.g006

While the above analysis is for a stable protein, a natural question to ask is how would these results change for an unstable protein? Consider an unstable protein with half-life considerably shorter than the cell-cycle duration. This rapid turnover ensures that the protein level equilibrates instantaneously after cell-division and gene-duplication events. Let γx denote the protein decay rate. Then, the mean protein level before and after genome duplication is and , respectively. Note that in the limit of large γx there is no noise contribution from partitioning errors since errors incurred at the time of cell division would be instantaneously corrected. The extrinsic noise, which can be interpreted as the protein noise level for deterministic protein production and decay is obtained as (for analysis on general f see Appendix G in S1 Text) (48) which is similar to noise level reported in [75]. When β = 0 or 1, the transcription rate and the protein level are constant within the cell cycle and . Moreover, is maximized at β = 2/3 with a value of 1/8. Thus, in contrast to a stable protein, extrinsic noise in an unstable protein is strongly dependent on the timing of gene duplication. Next, consider the intrinsic noise component. Analysis in Appendix G in S1 Text shows that the noise contribution from random protein production and decay is (49) While the mean protein level is strongly dependent on β, the intrinsic noise Fano factor is independent of it. Thus, similar to what was observed for a stable protein, the intrinsic noise in an unstable protein is invariant of β for a fixed .

Discussion

We have investigated a model of protein expression in bursts coupled to discrete gene-duplication and cell-division events. The novelty of our modeling framework lies in describing the size of protein bursts B, the time between cell birth and gene duplication T1, the time between gene duplication and cell division T2, and partitioning of molecules during cell division through general statistical distributions. Exact formulas connecting the protein mean and noise levels to these underlying distributions were derived. Furthermore, the protein noise level, as measured by the squared coefficient of variation, was decomposed into three components representing contributions from gene-duplication/cell-division events, stochastic expression and random partitioning. While the first component is independent of the mean protein level, the other two components are inversely proportional to it. Some important insights are as follows:

  • The mean protein level is affected by both the first and second-order moments of T1 and T2. In particular, randomness in these times (for a fixed mean) increases the average protein count. This increase can be attributed to the fact that increasing cell-cycle time variations leads to positively skewed distributions, making longer cell cycles (and hence higher protein accumulation) more likely.
  • Random gene-duplication/cell-division events create an extrinsic noise term which is completely determined by moments of T1 and T2 up to order three. Interestingly, noise in the timing of these events also critically affects the intrinsic noise contributions from stochastic expression and partitioning. Hence, ignoring the effect of cell-cycle time variations, may lead to erroneous estimation of intrinsic noise.
  • The noise contribution from partitioning errors decreases with increasing randomness in T1 and T2. Thus, if is sufficiently small and α is large compared to B in Eq (45), increasing noise in the timing of cell-cycle events decreases the total noise level.

A key limitation of our approach is to model timing of gene-duplication/cell-division events through independent random variables. There is always non-zero correlation in the cell-cycle durations of mother and daughter cells [76–78]. Moreover, in the same cell cycle, times T1 and T2 could be dependent [79]. While our assumption on independence of timing maybe unrealistic, it played an important role in deriving exact analytical formulas for the protein mean and noise levels. We have used Monte Carlo simulations to investigate scenarios where T1 and T2, or successive cell-division events, have memory and are dependent random variables (see Appendix H in S1 Text). Our analysis reveals that the results presented in Figs 4 and 6 hold even when the assumption of independent timing is perturbed over biologically meaningful parameters.

Effect of gene duplication on noise level

In this first-of-its-kind study, we have investigated how discrete f-fold changes in the transcription rate due to gene duplication affect the intercellular variability in protein levels. Not surprisingly, the timing of genome duplication strongly affects the mean protein level— can change up to f folds depending on whether the gene duplicates early (β = 0) or late (β = 1) in the cell cycle. Results show that genome duplication has counter intuitive effects on the protein noise level (Fig 6). For example, if stochastic expression is the dominant source of noise, then doubling of transcription due to duplication results in lower noise, as compared to constant transcription throughout the cell cycle. This is because for the same mean protein level, there are more burst (transcription) events in the case of genome duplication (f = 2) than constant transcription (f = 1). For example, consider deterministic timing () and gene duplication in the middle of the cell cycle (β = 0.5). Then, for the case β = 1, there are on average kx〈T〉 burst events per cell cycle. For the same , there are 1.05kx〈T〉 production events in the case of gene duplication (β = 0.5). This slight increase in the number of transcription events leads to better averaging of bursty protein synthesis and lower noise levels. Overall, the effect of β on different noise component is quite modest: as β varies, noise components deviate at maximum ≈ 20% from their values at β = 0 (Fig 6). These results are in contrast to the case of an unstable protein, where noise from the cell-cycle process is strongly dependent on β as shown in Eq (48).

Noise in synchronized cells

The mathematical framework introduced for modeling timing of cell division can be easily used to compute noise in synchronized cells. For example, let the cell-cycle duration be an Erlang distribution with shape parameter n and rate parameter nk (i.e., pn = 1 in Fig 3), which can be biologically interpreted as cells moving through n cell-cycle stages Gn1, Gn2, …, Gnn. Statistical moments conditioned on the cell-cycle stage Gnj can be obtained using (50) Using Eq (50) and moments obtained from Eqs (16) and (35)–(37), yields the following conditional mean (51) which increases with cell-cycle stage (i.e., higher values of j). The protein noise level given that cells are in stage Gnj is given by (52) Note that if n is large then the first term, which represents the noise contribution from the cell-cycle process, is negligible and can be dropped. Interesting, while the noise contribution from partitioning errors decreases with cell-cycle stage, the noise contribution from stochastic expression increases with j. Moreover, for B = 1 with probability 1 and α = 1, the intrinsic noise is always 1/Mean irrespective of j. Assuming high n, the noise at cell birth (j = 1) and division (j = n) are obtained as (53) (54) respectively. Thus, measurements of Eqs (53) and (54) by synchronizing cells (or by using cell size as a proxy for cell-cycle stage) can be used to quantify α and 〈B2〉/〈B〉, providing a novel way to separate these noise contributions. Next, we discuss how noise in asynchronous cell can be used to quantify these parameters.

Parameter inference in asynchronous cells

Simple models of bursty expression and decay predict the distribution of protein levels to be negative binomial (or gamma distributed in the continuous framework) [80, 81]. These distributions are characterized by two parameter—the burst arrival rate kx and the average burst size 〈B〉, which can be estimated from measured protein mean and noise levels. This method has been used for estimating kx and 〈B〉 across different genes in E. coli [47, 82]. Our detailed model that takes into account partitioning errors predicts (ignoring gene-duplication effects) (55) Using and a geometrically distributed B [50, 83–85], Eq (55) reduces to (56) Given measurements of intrinsic noise and the mean protein level, 〈B〉 can be estimated from Eq (56) assuming α = 1 (i.e., binomial partitioning). Once 〈B〉 is known, kx is obtained from the mean protein level given by Eq (17). Since for many genes 〈B〉 ≈ 0.5–5 [47], the contribution of the first term in Eq (56) is significant, and ignoring it could lead to overestimation of 〈B〉. Overestimation would be even more severe if α happen to be much higher than 1, as would be the case for proteins that form aggregates or multimers [33]. One approach to estimate both 〈B〉 and α is to measure intrinsic noise changes in response to perturbing 〈B〉 by, for example, changing the mRNA translation rate through mutations in the ribosomal-binding sites (RBS). Consider a hypothetical scenario where the Fano Factor (intrinsic noise times the mean level) is 6. Let mutations in the RBS reduces by 50%, implying a 50% reduction in 〈B〉. If the Fano factor changes from 6 to 4 due to this mutation, then 〈B〉 = 3.6 and 〈α〉 = 3.25.

Our recent work has shown that higher-order statistics of protein levels (i.e., skewness and kurtosis) or transient changes in protein noise levels in response to blocking transcription provide additional information for discriminating between noise mechanisms [86, 87]. Up till now these studies have ignored noise sources in the cell-cycle process. It remains to be seen if such methods can be used for separating the noise contributions of partitioning errors and stochastic expression to reliably estimate 〈B〉 and α.

Integrating cell size and promoter switching

An important limitation of our modeling approach is that it does not take into account the size of growing cells. Recent experimental studies have provided important insights into the regulatory mechanisms controlling cell size [88–91]. More specifically, studies in E. coli and yeast argue for an “adder” model, where cell-cycle timing is controlled so as to add a constant volume between cell birth and division [78, 91–93]. Assuming exponential growth, this implies that the time taken to complete cell cycle is negatively correlated with cell size at birth. More importantly, cell size also affects gene expression—in mammalian cells transcription rates linearly increase with the cell size [94]. Thus, as cells become bigger they also produce more mRNAs to ensure gene product concentrations remains more or less constant. A main direction of future work would be to explicitly include cell size with size-dependent expression and timing of cell division determined by the adder model. This formulation will for the first time, allow simultaneous investigation of stochasticity in cell size, protein molecular count and concentration.

Our study ignores genetic promoter switching between active and inactive states, which has been shown to be a major source of noise in the expression of genes across organisms [95–104]. Taking into account promote switching is particularly important for genome-duplication studies, where doubling the number of gene copies could lead to more efficient averaging of promoter fluctuations. Another direction of future work will be to incorporate this additional noise source into the modeling framework and investigate its contribution as a function of gene-duplication timing.

Supporting Information

S1 Text. Supplementary Material.

Additional proofs and simulations.

https://doi.org/10.1371/journal.pcbi.1004972.s001

(PDF)

Author Contributions

  1. Conceived and designed the experiments: MS DA AS.
  2. Performed the experiments: MS CAVG.
  3. Analyzed the data: MS CAVG.
  4. Wrote the paper: MS AS.

References

  1. 1. Blake WJ, Kaern M, Cantor CR, Collins JJ. Noise in eukaryotic gene expression. Nature. 2003;422:633–637. pmid:12687005
  2. 2. Raser JM, O’Shea EK. Noise in Gene Expression: Origins, Consequences, and Control. Science. 2005;309:2010–2013. pmid:16179466
  3. 3. Neuert G, Munsky B, Tan RZ, Teytelman L, Khammash M, van Oudenaarden A. Systematic Identification of Signal-Activated Stochastic Gene Regulation. Science. 2013;339:584–587. pmid:23372015
  4. 4. Libby E, Perkins TJ, Swain PS. Noisy information processing through transcriptional regulation. Proceedings of the National Academy of Sciences. 2007;104:7151–7156.
  5. 5. Fraser HB, Hirsh AE, Giaever G, Kumm J, Eisen MB. Noise Minimization in Eukaryotic Gene Expression. PLOS Biology. 2004;2:e137. pmid:15124029
  6. 6. Lehner B. Selection to minimise noise in living systems and its implications for the evolution of gene expression. Molecular Systems Biology. 2008;4:170. pmid:18319722
  7. 7. Losick R, Desplan C. Stochasticity and Cell Fate. Science. 2008;320:65–68. pmid:18388284
  8. 8. Arkin AP, Ross J, McAdams HH. Stochastic Kinetic Analysis of Developmental Pathway Bifurcation in Phage λ-Infected Escherichia coli Cells. Genetics. 1998;149:1633–1648. pmid:9691025
  9. 9. Weinberger L, Burnett J, Toettcher J, Arkin A, Schaffer D. Stochastic Gene Expression in a Lentiviral Positive-Feedback Loop: HIV-1 Tat Fluctuations Drive Phenotypic Diversity. Cell. 2005;122:169–182. pmid:16051143
  10. 10. Weinberger LS, Dar RD, Simpson ML. Transient-mediated fate determination in a transcriptional circuit of HIV. Nature Genetics. 2008;40:466–470. pmid:18344999
  11. 11. Singh A, Weinberger LS. Stochastic gene expression as a molecular switch for viral latency. Current Opinion in Microbiology. 2009;12:460–466. pmid:19595626
  12. 12. Dar RD, Hosmane NN, Arkin MR, Siliciano RF, Weinberger LS. Screening for noise in gene expression identifies drug synergies. Science. 2014;344:1392–1396. pmid:24903562
  13. 13. Eldar A, Elowitz MB. Functional roles for noise in genetic circuits. Nature. 2010;467:167–173. pmid:20829787
  14. 14. Veening JW, Smits WK, Kuipers OP. Bistability, Epigenetics, and Bet-Hedging in Bacteria. Annual Review of Microbiology. 2008;62:193–210. pmid:18537474
  15. 15. Kussell E, Leibler S. Phenotypic diversity, population growth, and information in fluctuating environments. Science. 2005;309:2075–2078. pmid:16123265
  16. 16. Balaban NQ, Merrin J, Chait R, Kowalik L, Leibler S. Bacterial persistence as a phenotypic switch. Science. 2004;305:1622–1625. pmid:15308767
  17. 17. Sánchez-Romero MA, Casadesús J. Contribution of phenotypic heterogeneity to adaptive antibiotic resistance. Proceedings of the National Academy of Sciences. 2014;111:355–360.
  18. 18. Neildez-Nguyen TMA, Parisot A, Vignal C, Rameau P, Stockholm D, Picot J, et al. Epigenetic gene expression noise and phenotypic diversification of clonal cell populations. Differentiation. 2008;76:33–40. pmid:17825084
  19. 19. Paldi A. Stochastic gene expression during cell differentiation: order from disorder? Cellular and Molecular Life Sciences. 2003;60:1775–1778. pmid:14523542
  20. 20. Raj A, van Oudenaarden A. Nature, nurture, or chance: stochastic gene expression and its consequences. Cell. 2008;135:216–226. pmid:18957198
  21. 21. Kaern M, Elston TC, Blake WJ, Collins JJ. Stochasticity in gene expression: from theories to phenotypes. Nature Reviews Genetics. 2005;6:451–464. pmid:15883588
  22. 22. Magklara A, Lomvardas S. Stochastic gene expression in mammals: lessons from olfaction. Trends in Cell Biology. 2013;23:449–456. pmid:23689023
  23. 23. Munsky B, Trinh B, Khammash M. Listening to the noise: random fluctuations reveal gene network parameters. Molecular systems biology. 2009;5:318. pmid:19888213
  24. 24. Wang P, Robert L, Pelletier J, Dang WL, Taddei F, Wright A, et al. Robust growth of Escherichia coli. Current biology. 2010;20:1099–1103. pmid:20537537
  25. 25. Lambert G, Kussell E. Quantifying Selective Pressures Driving Bacterial Evolution Using Lineage Analysis. Physical Review X. 2015;5:011016. pmid:26213639
  26. 26. Tsukanov R, Reshes G, Carmon G, Fischer-Friedrich E, Gov NS, Fishov I, et al. Timing of Z-ring localization in Escherichia coli. Physical Biology. 2011;8:066003. pmid:22015938
  27. 27. Reshes G, Vanounou S, Fishov I, Feingold M. Cell shape dynamics in Escherichia coli. Biophysical Journal. 2008;94:251–264. pmid:17766333
  28. 28. Reshes G, Vanounou S, Fishov I, Feingold M. Timing the start of division in E. coli: a single-cell study. Physical Biology. 2008;5:046001. pmid:18997273
  29. 29. Roeder AH, Chickarmane V, Chuna A, Obara B, Manjunath BS, Meyerowitz EM. Variability in the control of cell division underlies sepal epidermal patterning in Arabidopsis thaliana. PLOS Biology. 2010;8:e1000367. pmid:20485493
  30. 30. Zilman A, Ganusov VV, Perelson AS. Stochastic models of lymphocyte proliferation and death. PLOS ONE. 2010;5:e12775. pmid:20941358
  31. 31. Hawkins ED, Markham JF, McGuinness LP, Hodgkin PD. A single-cell pedigree analysis of alternative stochastic lymphocyte fates. Proceedings of the National Academy of Sciences. 2009;106:13457–13462.
  32. 32. Stukalin EB, Aifuwa I, Kim JS, Wirtz D, Sun SX. Age-dependent stochastic models for understanding population fluctuations in continuously cultured cells. Journal of The Royal Society Interface. 2013;101.
  33. 33. Huh D, Paulsson J. Random partitioning of molecules at cell division. Proceedings of the National Academy of Sciences. 2011;108:15004–15009.
  34. 34. Gonze D. Modeling the effect of cell division on genetic oscillators. Journal of Theoretical Biology. 2013;325:22–33. pmid:23434891
  35. 35. Lloyd-Price J, Tran H, Ribeiro AS. Dynamics of small genetic circuits subject to stochastic partitioning in cell division. Journal of Theoretical Biology. 2014;356:11–19. pmid:24768865
  36. 36. Zopf CJ, Quinn K, Zeidman J, Maheshri N. Cell-Cycle Dependence of Transcription Dominates Noise in Gene Expression. PLOS Computational Biology. 2013;9:e1003161. pmid:23935476
  37. 37. Narula J, Kuchina A, Lee DyD, Fujita M, Süel GM, Igoshin OA. Chromosomal Arrangement of Phosphorelay Genes Couples Sporulation and DNA Replication. Cell. 2015;162:328–337. pmid:26165942
  38. 38. Schwabe A, Bruggeman FJ. Contributions of Cell Growth and Biochemical Reactions to Nongenetic Variability of Cells. Biophysical Journal. 2014;107:301–313. pmid:25028872
  39. 39. Huh D, Paulsson J. Non-genetic heterogeneity from stochastic partitioning at cell division. Nature Genetics. 2011;43:95–100. pmid:21186354
  40. 40. Antunes D, Singh A. Quantifying gene expression variability arising from randomness in cell division times. Journal of Mathematical Biology. 2015; 71: 437–463. pmid:25182129
  41. 41. Yu J, Xiao J, Ren X, Lao K, Xie XS. Probing Gene Expression in Live Cells, One Protein Molecule at a Time. Science. 2006;311:1600–1603. pmid:16543458
  42. 42. Paulsson J. Model of stochastic gene expression. Physics of Life Reviews. 2005;2:157–175.
  43. 43. Shahrezaei V, Swain PS. Analytical distributions for stochastic gene expression. Proceedings of the National Academy of Sciences. 2008;105:17256–17261.
  44. 44. Singh A, Hespanha JP. Optimal Feedback Strength for Noise Suppression in Autoregulatory Gene Networks. Biophysical Journal. 2009;96:4013–4023. pmid:19450473
  45. 45. Jia T, Kulkarni RV. Intrinsic Noise in Stochastic Models of Gene Expression with Molecular Memory and Bursting. Physical Review Letters. 2011; 106: 058102. pmid:21405439
  46. 46. Alon U. An Introduction to Systems Biology: Design Principles of Biological Circuits. Chapman and Hall/CRC; 2006.
    • 47. Taniguchi Y, Choi PJ, Li GW, Chen H, Babu M, Hearn J, et al. Quantifying E. coli proteome and transcriptome with single-molecule sensitivity in single cells. Science. 2010;329:533–538. pmid:20671182
    • 48. Schwanhausser B, Busse D, Li N, Dittmar G, Schuchhardt J, Wolf J, et al. Global quantification of mammalian gene expression control. Nature. 2011;473:337–342. pmid:21593866
    • 49. Swain PS, Elowitz MB, Siggia ED. Intrinsic and extrinsic contributions to stochasticity in gene expression. Proceedings of the National Academy of Sciences. 2002;99:12795–12800.
    • 50. Berg OG. A model for the statistical fluctuations of protein numbers in a microbial population. Journal of Theoretical Biology. 1978;71:587–603. pmid:96307
    • 51. Rigney DR. Stochastic model of constitutive protein levels in growing and dividing bacterial cells. Journal of Theoretical Biology. 1979;76:453–480. pmid:439915
    • 52. Singh A, Hespanha JP. Stochastic hybrid systems for studying biochemical processes. Philosophical Transactions of the Royal Society of London A: Mathematical, Physical and Engineering Sciences. 2010;368:4995–5011.
    • 53. Daigle BJ, Soltani M, Petzold LR, Singh A. Inferring single-cell gene expression mechanisms using stochastic simulation. Bioinformatics. 2015;31:1428–1435. pmid:25573914
    • 54. Lagershausen S. Performance Analysis of Closed Queueing Networks. Lecture Notes in Economics and Mathematical Systems. Springer; 2013.
      • 55. Buchholz P, Kriege J, Felko I. Input Modeling with Phase-Type Distributions and Markov Models. Springer; 2014.
        • 56. Hespanha JP, Singh A. Stochastic Models for Chemically Reacting Systems Using Polynomial Stochastic Hybrid Systems. International Journal of Robust and Nonlinear Control. 2005;15:669–689.
        • 57. Singh A, Hespanha JP. Approximate Moment Dynamics for Chemically Reacting Systems. IEEE Transactions on Automatic Control. 2011;56:414–418.
        • 58. Gomez-Uribe CA, Verghese GC. Mass Fluctuation Kinetics: Capturing Stochastic Effects in Systems of Chemical Reactions through Coupled Mean-Variance Computations. Journal of Chemical Physics. 2007;126:024109. pmid:17228945
        • 59. Lee CH, Kim K, Kim P. A moment closure method for stochastic reaction networks. Journal of Chemical Physics. 2009;130:134107. pmid:19355717
        • 60. Goutsias J. Classical versus stochastic kinetics modeling of biochemical reaction systems. Biophysical Journal. 2007;92:2350–2365. pmid:17218456
        • 61. Gillespie CS. Moment-closure approximations for mass-action models. IET systems biology. 2009;3:52–58. pmid:19154084
        • 62. Soltani M, Vargas-Garcia CA, Singh A. Conditional moment closure schemes for studying stochastic dynamics of genetic circuits. IEEE Transactions on Biomedical Systems and Circuits. 2015;9:518–526.
        • 63. Wang H, Yuan Z, Liu P, Zhou T. Division time-based amplifiers for stochastic gene expression. Molecular BioSystems. 2015;11:2417–2428. pmid:26178011
        • 64. Hilfinger A, Paulsson J. Separating intrinsic from extrinsic fluctuations in dynamic biological systems. Proceedings of the National Academy of Sciences. 2011;108:12167–12172.
        • 65. Singh A, Soltani M. Quantifying Intrinsic and Extrinsic Variability in Stochastic Gene Expression Models. PLOS ONE. 2013;8:e84301. pmid:24391934
        • 66. Shahrezaei V, Ollivier JF, Swain PS. Colored extrinsic fluctuations and stochastic gene expression. Molecular Systems Biology. 2008;4:196. pmid:18463620
        • 67. Scott M, Ingalls B, Kaern M. Estimations of intrinsic and extrinsic noise in models of nonlinear genetic networks. Chaos. 2006;16:026107. pmid:16822039
        • 68. Ozbudak EM, Thattai M, Kurtser I, Grossman AD, van Oudenaarden A. Regulation of noise in the expression of a single gene. Nature Genetics. 2002;31:69–73. pmid:11967532
        • 69. Newman JRS, Ghaemmaghami S, Ihmels J, Breslow DK, Noble M, DeRisi JL, et al. Single-cell proteomic analysis of S. cerevisiae reveals the architecture of biological noise. Nature. 2006;441:840–846.
        • 70. Singh A, Razooky B, Cox CD, Simpson ML, Weinberger LS. Transcriptional Bursting from the HIV-1 Promoter Is a Significant Source of Stochastic Noise in HIV-1 Gene Expression. Biophysical Journal. 2010;98:L32–L34. pmid:20409455
        • 71. Bar-Even A, Paulsson J, Maheshri N, Carmi M, O’Shea E, Pilpel Y, et al. Noise in protein expression scales with natural protein abundance. Nature Genetics. 2006;38:636–643. pmid:16715097
        • 72. Yunger S, Rosenfeld L, Garini Y, Shav-Tal Y. Single-allele analysis of transcription kinetics in living mammalian cells. Nature Methods. 2010;7:631–633. pmid:20639867
        • 73. Lott SE, Villalta JE, Schroth GP, Luo S, Tonkin LA, Eisen MB. Noncanonical Compensation of Zygotic X Transcription in Early Drosophila melanogaster Development Revealed through Single-Embryo RNA-Seq. PLOS Biology. 2011;9:e1000590. pmid:21346796
        • 74. Walker N, Nghe P, Tans SJ. Generation and filtering of gene expression noise by the bacterial cell cycle. BMC Biology. 2016;14:1–10.
        • 75. Keren L, van Dijk D, Weingarten-Gabbay S, Davidi D, Jona G, Weinberger A, et al. Noise in gene expression is coupled to growth rate. Genome Research. 2015;25:1893–1902. pmid:26355006
        • 76. Siegal-Gaskins D, Crosson S. Tightly Regulated and Heritable Division Control in Single Bacterial Cells. Biophysical Journal. 2008;95:2063–2072. pmid:18469083
        • 77. Cerulus B, New AM, Pougach K, Verstrepen KJ. Noise and Epigenetic Inheritance of Single-Cell Division Times Influence Population Fitness. Current Biology. 2016; 26: 1138–1147. pmid:27068419
        • 78. Taheri-Araghi S, Bradde S, Sauls JT, Hill NS, Levin PA, Paulsson J, et al. Cell-Size Control and Homeostasis in Bacteria. Current Biology. 2015;25:385–391. pmid:25544609
        • 79. Adiciptaningrum A, Osella M, Moolman MC, Lagomarsino MC, Tans SJ. Stochasticity and homeostasis in the E. coli replication and division cycle. Scientific Reports. 2015;5:18261. pmid:26671779
        • 80. Friedman N, Cai L, Xie XS. Linking stochastic dynamics to population distribution: an analytical framework of gene expression. Physical Review Letters. 2006;97:168302. pmid:17155441
        • 81. Paulsson J, Ehrenberg M. Random Signal Fluctuations Can Reduce Random Fluctuations in Regulated Components of Chemical Regulatory Networks. Physical Review Letters. 2000; 84:5447–5450. pmid:10990965
        • 82. Sherman MS, Cohen BA. A Computational Framework for Analyzing Stochasticity in Gene Expression. PLOS Computational Biology. 2014;10:e1003596. pmid:24811315
        • 83. Golding I, Paulsson J, Zawilski SM, Cox EC. Real-time kinetics of gene activity in individual bacteria. Cell. 2005;123:1025–1036. pmid:16360033
        • 84. McAdams HH, Arkin A. Stochastic mechanisms in gene expression. Proceedings of the National Academy of Sciences. 1997; 94:814–819. pmid:9023339
        • 85. Cai L, Friedman N, Xie XS. Stochastic protein expression in individual cells at the single molecule level. Nature. 2006;440:358–362. pmid:16541077
        • 86. Kumar N, Singh A, Kulkarni RV. Transcriptional bursting in gene expression: analytical results for general stochastic models. PLOS Computational Biology. 2015;11:e1004292. pmid:26474290
        • 87. Singh A. Transient Changes in Intercellular Protein Variability Identify Sources of Noise in Gene Expression. Biophysical Journal. 2014;107:2214–2220. pmid:25418106
        • 88. Osella M, Nugent E, Lagomarsino MC. Concerted control of Escherichia coli cell division. Proceedings of the National Academy of Sciences. 2014;111:3431–3435.
        • 89. Robert L, Hoffmann M, Krell N, Aymerich S, Robert J, Doumic M. Division in Escherichia coli is triggered by a size-sensing rather than a timing mechanism. BMC Biology. 2014;12:17. pmid:24580833
        • 90. Kafri R, Levy J, Ginzberg MB, Oh S, Lahav G, Kirschner MW. Dynamics extracted from fixed cells reveal feedback linking cell growth to cell cycle. Nature. 2013;494:480–483. pmid:23446419
        • 91. Ghusinga KR, Vargas-Garcia CA, Singh A. A mechanistic stochastic framework for regulating bacterial cell division. Scientific Reports. 2016; 6: 30229.
        • 92. Amir A. Cell Size Regulation in Bacteria. Physical Review Letters. 2014;112:208102.
        • 93. Campos M, Surovtsev IV, Kato S, Paintdakhi A, Beltran B, Ebmeier SE, et al. A Constant Size Extension Drives Bacterial Cell Size Homeostasis. Cell. 2014;159:1433–1446. pmid:25480302
        • 94. Padovan-Merhar O, Nair GP, Biaesch AG, Mayer A, Scarfone S, Foley SW, et al. Single Mammalian Cells Compensate for Differences in Cellular Volume and DNA Copy Number through Independent Global Transcriptiona Mechanisms. Molecular Cell. 2015;58:339–352. pmid:25866248
        • 95. Suter DM, Molina N, Gatfield D, Schneider K, Schibler U, Naef F. Mammalian genes are transcribed with widely different bursting kinetics. Science. 2011;332:472–474. pmid:21415320
        • 96. Brown CR, Mao C, Falkovskaia E, Jurica MS, Boeger H. Linking Stochastic Fluctuations in Chromatin Structure and Gene Expression. PLOS Biology. 2013;11:e1001621. pmid:23940458
        • 97. Raj A, Peskin CS, Tranchina D, Vargas DY, Tyagi S. Stochastic mRNA synthesis in mammalian cells. PLOS Biology. 2006;4:e309. pmid:17048983
        • 98. Hornung G, Bar-Ziv R, Rosin D, Tokuriki N, Tawfik DS, Oren M, et al. Noise-mean relationship in mutated promoters. Genome Research. 2012;22:2409–2417. pmid:22820945
        • 99. Singh A, Razooky BS, Dar RD, Weinberger LS. Dynamics of protein noise can distinguish between alternate sources of gene-expression variability. Molecular Systems Biology. 2012;8:607. pmid:22929617
        • 100. Dar RD, Razooky BS, Singh A, Trimeloni T, McCollum J, Cox C, et al. Transcriptional burst frequency and burst size are equally modulated across the human genome. Proceedings of the National Academy of Sciences. 2012;109:17454–17459.
        • 101. Corrigan AM, Chubb JR. Regulation of Transcriptional Bursting by a Naturally Oscillating Signal. Current Biology. 2014;24:205–211. pmid:24388853
        • 102. Bothma JP, Garcia HG, Esposito E, Schlissel G, Gregor T, Levine M. Dynamic regulation of eve stripe 2 expression reveals transcriptional bursts in living Drosophila embryos. Proceedings of the National Academy of Sciences. 2014;111:10598–10603.
        • 103. Chubb JR, Trcek T, Shenoy SM, Singer RH. Transcriptional Pulsing of a Developmental Gene. Current Biology. 2006;16:1018–1025. pmid:16713960
        • 104. Chong S, Chen C, Ge H, Xie XS. Mechanism of Transcriptional Bursting in Bacteria. Cell. 2014;158:314–326. pmid:25036631