Skip to main content
  • Research article
  • Open access
  • Published:

Effects of promoter leakage on dynamics of gene expression

Abstract

Background

Quantitative analysis of simple molecular networks is an important step forward understanding fundamental intracellular processes. As network motifs occurring recurrently in complex biological networks, gene auto-regulatory circuits have been extensively studied but gene expression dynamics remain to be fully understood, e.g., how promoter leakage affects expression noise is unclear.

Results

In this work, we analyze a gene model with auto regulation, where the promoter is assumed to have one active state with highly efficient transcription and one inactive state with very lowly efficient transcription (termed as promoter leakage). We first derive the analytical distribution of gene product, and then analyze effects of promoter leakage on expression dynamics including bursting kinetics. Interestingly, we find that promoter leakage always reduces expression noise and that increasing the leakage rate tends to simplify phenotypes. In addition, higher leakage results in fewer bursts.

Conclusions

Our results reveal the essential role of promoter leakage in controlling expression dynamics and further phenotype. Specifically, promoter leakage is a universal mechanism of reducing expression noise, controlling phenotypes in different environments and making the gene produce generate fewer bursts.

Background

Gene expression dynamics is a lasting issue in Systems Biology and has attracted extensive attention. While recent advances in experimental methods allow direct observations of gene expression levels in individual cells, there is considerable interest in theoretically understanding how different molecular mechanisms of gene expression influence variations in mRNA or protein levels across a population of cells. In fact, quantifying the contributions of different sources of noise using stochastic models of gene expression is an important step towards understanding fundamental cellular processes and cell-to-cell variability.

Many theoretical models of gene expression have been proposed and become more and more subtle, from the initial one-state model [1,2] to the common two-state model [3] to those considering many detailed processes or factors, such as chromatin remodeling [4-6], TATA-box- mediated promoter [7-9], transcription additional re-initiation [10], growth rate [11], copy number variations [12], recruitment of transcription factors [13,14], alternative splicing [15]. In spite of these, few gene models in previous studies considered promoter leakage. Here, by promoter leakage we mean that transcription efficiency at the promoter inactive state is much lower than that at the promoter active state. In fact, it has been experimentally verified that transcription takes place not only at the active state but also at the inactive state of promoter, e.g., different nucleosome protein binding sites can lead to different expression efficiencies, some of which are high whereas the others are very low [16,17]; a basal transcription rate at each open promoter state implies that the promoter has leakage; and transcription can take place at some promoter state with a very low rate due to the pre-initiation complex formed at the TATA box [7]. Some studies have shown that increasing the leakage rate may eliminate bistability [18] whereas decreasing the leakage rate of protein production can lead to persistent oscillations [19,20]. But there has been no systematic study on how promoter leakage affects dynamics of gene expression including the noise in gene product, probability distribution, and bursting kinetics (characterized by burst size and burst frequency). This paper will address these issues by analyzing three cases: the common ON-OFF model, the ON-OFF model with negative regulation and the ON-OFF model with positive regulation, referring to schematic Figure 1(A). We will present a systematic investigation for each case.

Figure 1
figure 1

A schematic diagram for effect of promoter leakage on expression dynamics. (A) schematic description of two-state gene auto-regulatory model with leakage, where gene product is produced not only at the ON state with high efficiency but also at the OFF state with very low efficiency (termed as promoter leakage); (B) representative time series for changes in the number of gene products (left) and effect of leakage on population distribution (right).

To be clear, we will analyze relevant literatures and present the major findings of this study from three different aspects.

First, how the architecture of promoter affects expression noise is a hot issue and has attracted attention of many scholars [21,22]. Here, the promoter architecture mainly includes the pattern of transitions among promoter activity states and exits of transcription [22,23]. Some transcription rates may be very small (this corresponds to promoter leakage). On the other hand, it is known that feedback is a ubiquitous mechanism controlling signals. In general, the effect of feedback on expression noise depends on the type of feedback (positive or negative) and model parameters [24]. It has been shown that direct negative feedback decreases fluctuations and is a ubiquitous mechanism for homoeostatic control [25,26]. A question naturally arises: in the case of feedback loop, how does promoter leakage impact expression noise? By analyzing effects of promoter leakage on gene expression using the above-mentioned three gene models, we find that promoter leakage always reduces expression noise, regardless of models, feedback types and model parameters as well as whether or not the expression level is fixed. This implies an important fact that promoter leakage is a universal mechanism of reducing expression noise.

Second, it is well known that probability distribution not only provides the most complete information on stochastic dynamics of a stochastic model but also can clearly indicate the system’s states in the stochastic sense since a different peak in the distribution corresponds to a different state of the system. In addition, peakedness of probability distribution has biological implication, e.g., bimodal or multimodal gene expression (i.e., mRNA or protein distribution exhibits two or multiple peaks) is a cause of phenotypic diversity in genetically identical cell populations [27]; the amount of phenotypic variation in gene product population can determine fitness by affecting growth rate, robustness and adaptation [27-31]. Increased phenotypic variation can enhance adaptation and growth of cells in fluctuating environments as well as robustness of the population to external stresses [32-34]. It is particularly beneficial to microbial cells that need to adapt efficiently to sudden changes in environmental conditions [35,36]. And decreased phenotypic variation, where every gene product is as close as possible to an optimal level of gene expression that maximizes fitness, is usually advantageous to cells in constant environments [28,37,38]. Given the above facts or results, an unsolved question is how promoter leakage influences probability distribution and phenotypic heterogeneity. Although a theoretical study [39] illustrated the impact of promoter leakage on a probability distribution, the corresponding model did not consider feedback. We find that promoter leakage tends to unimodalize the mRNA or protein distribution regardless of feedback type, implying that promoter leakage may cause phenotypic simplification. In a word, phenotypic diversity is beneficial to survival of cells in a fluctuated environment on the one hand and phenotypic oneness is beneficial to fitness of cells in a constant environment on the other hand. Through our study, we conclude that promoter leakage is a mechanism of controlling phenotypes in different environments: cells would make use of not only the mechanism of promoter leakage-induced diverse phenotypes to adapt changes in environments but also the mechanism of promoter leakage-induced optimal phenotypes to adapt to constant environments.

Third, it has been shown that mRNAs or proteins are produced often in a bursty manner [40-42]. In fact, single molecule measurements have provided evidence for transcriptional or translational bursting [9,43,44]. An accepted view is that bursting dynamics are responsible for the generation of cellular heterogeneity in the response of genetically identical cells to the same stimulus, e.g., the authors in Ref. [34] demonstrated in yeast cells that high levels of cell-to-cell variability, originated from slow promoter state fluctuations, may confer cell colonies with an enhanced probability of survival when subjected to external stresses, such as addition of high concentrations of antibiotic. Regarding bursting dynamics, there have been many studies. For example, for a two-state gene model, it has been shown that the higher the burst frequency (BF) is, the lower is the gene expression noise whereas the larger the burst size (BS) is, the higher is the noise [42,45,46]. More interestingly, R.D. Dar, et al. [42] showed by analyzing 8,000 individual human genomic loci that transcriptional bursting dominates across the human genome, both BF and BS change by chromosomal location, and transcriptional activators alter BF and BS, depending on the expression level of the locus. In addition, G. Hornung, et al. [8] showed that burst size is a promoter-specific property that is relatively robust to sequence mutations but depends strongly on the interaction between the TATA box and promoter nucleosomes. In spite of these studies, it is not clear whether the results obtained in the case of neither feedback nor promoter leakage are still correct in the case of feedback or promoter leakage or both. We will address this issue. By model analysis, we find that a higher leakage rate of promoter produces fewer bursts, regardless of feedback type, referring to Figure 1(B).

Results

Gene model and analytical distribution

In principle, the chemical master equation (CME) provides the most complete model for probabilistic behavior of any biochemical reaction network including gene auto-regulatory circuits. Alternatively, the joint probability density function (or the joint probability distribution) for all the reactive species in a reaction network provides the most complete information on stochastic properties of this system. Therefore, finding probability distribution becomes a common interest in understanding all the possible stochastic properties of a biochemical reaction network. Here, we will derive the analytical expression for the steady-state distribution in a two-state gene auto regulatory model with promoter leakage.

Before presenting our analytical result, let us simply introduce our model to be investigated. Assume that a gene has two activity states (we will use D to represent promoter state): the active one (or the ON state) where DNAs are transcribed into mRNAs that are then translated into proteins (denoted by P), and the inactive one (or the OFF state) where transcription is not nonexistent but is lowly efficient, that is, there is a smaller transcription rate at the inactive state in contrast to that at the active state. The latter case is often called promoter leakage [47,48]. The promoter leakage was neglected in previous studies of gene models [49-52], or equally, the transcription rate at the inactive state was assumed as zero. In addition, we assume that the gene product regulates its expression as a transcription factor, thus forming a feedback loop. This regulation may be positive or negative. In particular, we assume that the transcription factor regulates the gene expression not in a manner of sequestration (by sequestration we mean that dissociation of protein from transcription factor-DNA complex is a slow process, or binding of transcription factor-DNA is strong) [53,54] but in a manner of directly changing transition rates between promoter states (i.e., assuming that association and dissociation of a transcription factor are a very fast process. Specifically, the transcription factor first binds fast to DNA and then dissociate rapidly from the DNA, leading that the transcription factor changes directly the switching rates between promoter states without consuming itself) [55,56]. To simplify our analysis, we further integrate transcription and translation into a single-step process. This simplification has been extensively made [52,57]. In fact, it has been shown that protein’s half-life is in general much longer than that of mRNA’s half-life [52,57], e.g., Shahrezaei and Swain [52] did a survey for ~ 2000 genes in budding yeast and found that the expressions of most of these genes satisfy this condition. A main reason for this consideration is to derive the analytical expression for the probability distribution of gene product or to give the analytical formula for calculating the noise intensity, which in turn can clearly describe our qualitative results. In the Additional file 1, we also investigate slightly more complex gene models, which consider either a slow process for binding and dissociation of a transcription factor or two-step process of transcription and translation. By stochastic simulation, we find that the results using the complex models are basically similar to those using the simplified model, indicating that the latter can well capture the effect of promoter leakage on expression dynamics.

A transcription factor may regulate gene expression in an enhancing or a repressing manner. To derive our analytical distributions in a unified framework, we introduce the following set of biochemical reactions based on the above hypotheses.

$$ \begin{array}{l}{D}_0\overset{\gamma_1}{\to }{D}_1,{D}_1\overset{\gamma_0}{\to }{D}_0,{D}_1+P\overset{f}{\to }{D}_0+P\\ {}{D}_1\overset{\lambda_1}{\to }{D}_1+P,{D}_0\overset{\lambda_0}{\to }{D}_0+P,P\overset{d}{\to}\phi \end{array} $$
(1)

which describe our gene model with positive or negative regulation. Here, D 1 and D 0 represent the ON and the OFF states respectively if λ 1λ 0 and the OFF and the ON states respectively if λ 0λ 1. In the former case, the third reaction in the first row of Eq. (1) describes the negative regulation and λ 1λ 0 are transcription rate and promoter leakage rate respectively, whereas in the latter case, this reaction describes the positive regulation and λ 0λ 1 are transcription rate and promoter leakage rate respectively. Parameters γ 1 and γ 0 are transition rates between the promoter activity states, f represents feedback strength, and d is the degradation rate of gene product.

Then, we establish our mathematical model to be studied. Let P 0(n, t) and P 1(n, t) represent the probability that the gen product has n molecules at time t when the gene is at D 0 and D 1 states, respectively. Then, the discrete CME for the full reaction network takes the form

$$ \frac{\partial }{\partial t}\left(\begin{array}{c}\hfill {P}_0\left(n,t\right)\hfill \\ {}\hfill {P}_1\left(n,t\right)\hfill \end{array}\right)=\left(\begin{array}{cc}\hfill -{\gamma}_1\hfill & \hfill {\gamma}_0+fn\hfill \\ {}\hfill {\gamma}_1\hfill & \hfill -{\gamma}_0-fn\hfill \end{array}\right)\left(\begin{array}{c}\hfill {P}_0\left(n,t\right)\hfill \\ {}\hfill {P}_1\left(n,t\right)\hfill \end{array}\right)+\left(\begin{array}{c}\hfill \left[{\lambda}_0\left({E}^{-1}-I\right)+d\left(E-I\right)\right]{P}_0\left(n,t\right)\hfill \\ {}\hfill \left[{\lambda}_1\left({E}^{-1}-I\right)+d\left(E-I\right)\right]{P}_1\left(n,t\right)\hfill \end{array}\right) $$
(2)

where E is the step operator and I is the identity operator.

Next, we focus on finding the steady-state solution of Eq. (2). The basic idea is first to introduce probability-generating functions, and then to solve a coupled set of ordinary differential equations with respect to these functions. The overall procedure for finding the stationary distribution is technical. Here, we list results only, and the details for derivation are put in the Additional file 2 of this paper.

For convenience, let all the parameters be normalized by the degradation rate d, that is, γ 1/d → γ 1, γ 0/d → γ 0, f/d → f, λ 0/d → λ 0, λ 1/d → λ 1. Then, then analytical probability distribution can be expressed as

$$ P(n)=\frac{gA}{n!}{\displaystyle \sum_{m=0}^n\left(\begin{array}{c}\hfill n\hfill \\ {}\hfill m\hfill \end{array}\right){\lambda}_0^{n-m}}{\left[\left(f+1\right)Q\right]}^m\frac{{\left(\alpha -1\right)}_m}{{\left(\beta -1\right)}_m}{}_1F_1\left(\alpha +m-1,\beta +m-1;-Q\right) $$
(3)

which is the linear superposition of confluent hypergeometric functions. In Eq. (3), λ = λ 1 − λ 0 R = λ − f λ 0, \( Q=\frac{\lambda -f\kern0.1em {\lambda}_0}{{\left(f+1\right)}^2} \), \( g=\frac{\lambda +{\gamma}_1+{\gamma}_0}{\gamma_1}-\frac{R}{\gamma_1\left(f+1\right)} \), \( \beta =1+\frac{\lambda +{\gamma}_0+{\gamma}_1}{f+1}-\frac{R}{{\left(f+1\right)}^2} \) \( \alpha =1+\frac{\lambda {\gamma}_1}{R} \) and \( A={e}^{-{\lambda}_0}{\left[g{}_1F_1\left(\alpha -1,\beta -1;f\kern0.1em Q\right)\right]}^{-1} \). The symbol \( \left(\begin{array}{c}\hfill n\hfill \\ {}\hfill m\hfill \end{array}\right) \) represents the combinational number of choosing m molecule from n molecules, and (c) n is the Pochhammer symbol and is defined as (c) n  = Γ(c + n)/Γ(c). In principle, this analytical distribution gives all the stochastic information about the underlying gene model.

The more useful is that we can give the analytical formula for calculating the noise intensity for the gene product, where by noise intensity we mean that it is the ratio of variance over the square of mean. In fact, note that mean and variance can be calculated according to the following general formulae.

$$ \left\langle n\right\rangle ={G}^{\prime }(1),{\sigma}_n^2={G}^{{\prime\prime} }(1)+{G}^{\prime }(1)-{\left[{G}^{\prime }(1)\right]}^2 $$
(4)

where G(z) = G 0(z) + G 1(z) is the total generating function. Two factorial generating functions G 0 (z) and G 1(z) are analytically given in the Additional file 2. By calculation, we find

$$ {G}^{\prime }(1)=A{e}^{\lambda_0}\left[g{\lambda}_0{}_1F_1\left(\alpha -1,\kern0.5em \beta -1,\kern0.5em f\kern0.1em Q\right)+\lambda {}_1F_1\left(\alpha, \kern0.5em \beta, \kern0.5em f\kern0.1em Q\right)\right] $$
(5)
$$ \begin{array}{l}{G}^{{\prime\prime} }(1)=A{e}^{\lambda_0}\left[g{\lambda}_0^2{}_1F_1\left(\alpha -1,\beta -1,f\kern0.1em Q\right)+2\lambda {\lambda}_0{}_1F_1\left(\alpha, \beta, f\kern0.1em Q\right)\right.\\ {}\begin{array}{cc}\hfill \hfill & \hfill \hfill \end{array}+\frac{\lambda \left(\lambda {\gamma}_1+R\right)}{\lambda +{\gamma}_1+{\gamma}_0+f+1-R/\left(1+f\right)}{}_1F_1\left(\alpha +1,\beta +1,f\kern0.1em Q\right)\Big]\end{array} $$
(6)

Thus, the noise intensity is given according to the formula.

$$ {\eta}_n^2=\frac{{\sigma_n}^2}{\left\langle n\right\rangle {}^2}=\frac{G^{{\prime\prime} }(1)+{G}^{\prime }(1)-{\left[{G}^{\prime }(1)\right]}^2}{{\left[{G}^{\prime }(1)\right]}^2} $$
(7)

Next, we focus on effects of promoter leakage on expression noise as well as burst dynamics. Note that the promoter leakage rate is λ 0 in the case of negative feedback whereas it is λ 1 in the case of positive feedback. As pointed out in the previous section, by promoter leakage we mean that the parameter λ 0 is in general more than zero but much smaller than λ 1 in the former case. Conversely, λ 1 is in general more than zero but much smaller than λ 0 in the latter case. We point out that to demonstrate the more remarkable numerical effects of promoter leakage, we sometimes let the leakage rate be a little smaller but be not much smaller than the normal transcription rate. For convenience, we always set d = 1 in our simulation.

Leakage always attenuates expression noise

In previous studies [49-52], the effect of promoter leakage on gene expression was frequently neglected. Here, we numerically show that the promoter leakage has unneglectable effects on gene expression and in particular on expression noise. More precisely, the promoter leakage always reduces the noise in gene product. The numerical results are shown in Figure 2.

Figure 2
figure 2

Effects of promoter leakage on gene expression noise. (A) The dependence of the noise intensity on the promoter leakage rate in the case that the gene product amount is fixed at a certain value, showing the noise intensity is always a monotonically decreasing function of the leakage rate, regardless of ways to keep the average expression level fixed (e.g., decreasing the transition rate from OFF to ON (γ 1) (solid red line); increasing the transition rate from ON to OFF (γ 0) (black dotted line); increasing the feedback strength (f) (green dashed line); decreasing the maximum transcription rate (λ 1) (blue dash dot line). (B) Results in the case that the mean expression is not fixed, showing that increasing the leakage rate reduces the expression noise, where 4 colored lines correspond to 4 different sets of parameter values: λ 1 = 40, γ 0 = 0.1, γ 1 = 0.1, f = 0.01(green); λ 1 = 40, γ 0 = 0.1, γ 1 = 0.2, f = 0(red); λ 1 = 40, γ 0 = 0.2, γ 1 = 0.1, f = 0(black); λ 1 = 30, γ 0 = 0.1, γ 1 = 0.1, f = 0(blue). This subfigure shows that the conclusion that promoter leakage always reduces noise is independent of model parameters. (C) The noise as a function of the mean for different values of the leakage rate, where the parameter values are the same as those used in Figure 2 (B). The subfigure shows that the larger the promoter leakage rate is, the more is the number of gene product molecules, implying that promoter leakage always reduces expression noise. In Figure 2 (B) and (C), lines represent theoretical results whereas circles represent stochastically simulating results.

In our numerical calculation, we let the leakage rate increase but keep the mean expression of gene product fixed. This constraint implies that two of the system parameters are dependent of each other if the other parameters are fixed. For example, consider the case of λ 1λ 0. If gene product is kept at a fixed level, then increasing leakage rate (λ 0) implies: 1), a decrease in the transition rate from OFF to ON (γ 1); 2), an increase in the transition rate from ON to OFF (γ 0); 3), an increase in the feedback strength (f); or 4), a decrease in the maximum transcription rate λ 1. Specifically, in Figure 2(A), changing λ 0 implies varying γ 1 (solid red line) if we fix λ 1 = 40, γ 0 = 0.1, and f = 0; changing λ 0 implies varying γ 0 (black dotted line) if we fix λ 1 = 40, f = 0, and γ 1 = 0.1; changing λ 0 implies varying f (green dashed line) if we fix λ 1 = 40, γ 0 = 0.1 and γ 1 = 0.1; changing λ 0 implies varying λ 1 (blue dash dot line) if we fix γ 0 = 0.1, γ 1 = 0.1 and f = 0. We point out that this classification of parameter values implies that positive and negative feedbacks can be considered in a unified framework. Therefore, it is unnecessary to distinguish the case of positive feedback from that of negative feedback for simulation.

From Figure 2(A), we observe that in the case that the average expression level is fixed, the noise intensity η n is a monotonically decreasing function of the promoter leakage rate, regardless of ways that the mean is fixed. In other words, the noise intensity always decreases with the increase of the leakage rate. This implies that the promoter leakage plays a role of attenuating the noise in the gene product, no matter what the property of feedback (positive or negative). Figure 2(B) is used to show that the conclusion that the promoter leakage always reduces the expression noise does not depend on the choice of parameter values, thus being qualitatively invariant. In this figure, we choose 4 different sets of parameter values (see the caption of Figure 2(B)) to demonstrate numerical results. The combination of Figures 2(A) and (B) implies that promoter leakage is a mechanism of efficiently reducing expression noise. To check if this mechanism is universal, we also investigate slightly more complex models in the Additional file 1, which consider either strong binding of transcription factors to DNAs or the two-step processes of transcription and translation. By numerical simulation, we find that the qualitative conclusion that promoter leakage always reduces noise still keeps invariant, referring to Additional file 1: Figure S1 (C), Figure S2 (C) and Figure S3 (C).

Intuitively, the larger the promoter leakage rate is, the more is the number of gene-product molecules. This will lead to the reduction of the noise in gene product (see Figure 2(C)). This result is easily seen in the case that the mean level is not fixed but seems not apparent in the case that the mean level is fixed. However, if the mean expression level is fixed, then the above qualitative conclusion seems unrelated to the property or type of feedback. See the following content for interpretation.

Promoter leakage tends to unimodalize distribution

Several examples have shown that the noise in gene expression is a potential mechanism to generate phenotypic heterogeneity [27,49,58]. The phenotypic diversity has been a focus of attention in biology, since the amount of phenotypic variation (also known as gene expression noise) in a cell population can determine fitness by affecting growth rate, robustness and adaptation [28-31]. Population diversity offers an alternate way that cells adapt to randomly fluctuating environments [30]. Increasing phenotypic variation is particularly beneficial to organisms that need to adapt efficiently to sudden changes in chemical composition, local temperature, or illumination [30]. In contrast, decreasing phenotypic variation is usually advantageous in constant environments [28,37,38]. The peaks of gene product distribution are a cause of generating phenotypic diversity in genetically identical cell populations [27], and the gene product noise always decreases with the increase of the leakage rate as shown above, so we naturally consider the relationship between promoter leakage and phenotypic selection. For this, we will investigate the effect of promoter leakage on peakedness of probability distribution.

It has been shown that a two-state gene model can exhibit bimodal distributions if neither promoter leakage nor auto regulation is considered [28,59]. Bimodality can also occur even in the presence of regulation [60,61]. If promoter leakage is considered, however, we find that the situation is different. The numerical results are demonstrated in Figure 3, where we consider three cases: no feedback, which corresponds to f = 0; negative feedback, which implies that D 1 represents the active state but λ 1 corresponds to the normal transcription rate whereas λ 0 to the leakage rate; and positive feedback, which implies that D 0 represents the active state but λ 0 corresponds to the normal transcription rate whereas λ 1 to the leakage rate.

Figure 3
figure 3

Effects of promoter leakage on distribution. (A) no feedback (f = 0): the gene product distribution changes from bimodality to unimodality when the leakage rate increases. The parameter values are set as λ 1 = 40, γ 0 = 0.1, γ 1 = 0.2; (B) negative feedback: only one peak closed to the origin gradually becomes another peak away from the origin with the increase of the leakage rate. Other parameter values are set as λ 1 = 40, γ 0 = 0.1, γ 1 = 0.2, f = 0.1; (C) positive feedback: two peaks gradually become one peak away from the origin with the increase of the leakage rate, where the parameter values are set as λ 0 = 40, γ 0 = 0.1, γ 1 = 0.5, f = 0.1.

From Figure 3(A) where no feedback is considered, we observe that with the increase of the leakage rate from zero to 25, the gene product distribution finally becomes a single peak from initial two peaks. Figure 3(B) shows that in the case of negative feedback, the gene product distribution becomes a single peak away from the origin from one single peak near the origin with the increase of the leakage rate from zero to 10. In contrast, Figure 3(C) shows that in the case of positive feedback, the gene product distribution becomes a single peak away from the origin from two distinct peaks with the increase of the leakage rate from zero to 10. Compared to the case in Figure 3(A), this change in Figure 3(C) is faster.

Thus, we obtain an interesting conclusion, that is, promoter leakage tends to unimodalize the distribution, independently of the property or type of feedback. We point out that neither transcription factor long residency on the DNA nor transcription and translation in different time scales influences this qualitative conclusion. In fact, we have shown in the Additional file 1 that if these factors are considered, then the unimodal distribution of gene product will become more apparent with the increase of the promoter leakage rate, referring to Additional file 1: Figure S1 (B), Figure S2 (B) and Figure S3 (B). Since multimodality is an important source resulting in the diversity of phenotype and since the unimodal distribution implies the singleness of phenotype. While increased phenotypic diversity can enhance adaptation and growth of cells in fluctuating environments [33], so our conclusion provides an important hint, that is, promoter leakage would not be too remarkable in in vivo organisms.

As mentioned in the introduction, diversity of phenotype is beneficial to survival of cells in a fluctuated environment whereas oneness of phenotype is beneficial to fitness of cells in a constant environment [35-38]. Thus, the above results imply that promoter leakage would be a mechanism of effectively controlling phenotype in different environments: cells would make use of not only the mechanism of promoter leakage-induced diverse phenotypes to adapt fluctuated environments but also the mechanism of promoter leakage-induced optimal phenotypes to adapt constant environments.

Promoter leakage can result in fewer bursts

As is well known, mRNAs or proteins are synthesized often in a burst manner [40-42]. Bursting kinetics is commonly characterized by two indices: burst size (BS) and burst frequency (BF). A question is how promoter leakage impacts bursting kinetics. Previous studies did not give a positive answer to this issue although it has been shown that the larger the BS is, the higher is the gene expression noise whereas the higher the BF is, the lower is the noise [42,46,47]. Here, we will show that promoter leakage has unneglectable influences on bursting kinetics, remarkably making the gene product produce fewer bursts.

Before presenting results, let us simply introduce computation formulae associated with bursting kinetics. Recall that in the case of no feedback, the mean BF and the mean BS are calculated according to the following formulae [42]

$$ \left\langle BF\right\rangle =\frac{1}{\tau_{OFF}},\left\langle BS\right\rangle ={k}_{transcription}\cdot {\tau}_{ON} $$
(8)

where τ OFF and τ ON represent the mean time dwelling at OFF and ON states respectively, and k transcription is the transcription rate when the gene is at ON state. First, consider that our model has no feedback, i.e., f = 0. If λ 0 compared to λ 1 is so small that it may be ignored (in this case, D 1 represents the ON state), then we have τ OFF  = 1/γ 1 [62,63], which implies 〈BF〉 = γ 1, and τ ON  = 1/γ 0 [62,63], which implies 〈BS〉 = λ 1/γ 0. Similarly, if λ 1 compared to λ 0 is so small that it may be ignored (in this case, D 0 represents the ON state), then we have τ OFF  = 1/γ 0, which implies 〈BF〉 = γ 0, and τ ON  = 1/γ 1, which implies 〈BS〉 = λ 0/γ 1. Thus, for the former case, we see from Eq. (11) that increasing the leakage rate λ 0 does not change the mean BF nor change the mean BS, referring to Figure 4(A) and (D), but may make the original low expression amount of the gene product have a rise. Moreover, it can be seen from the time series shown in Figure 5 (more precisely, by comparing Figure 5(A) which corresponds to the case of no leakage with Figure 5(D) which corresponds to the case that the leakage rate is 10) that the number of gene product molecules tends to centralize a certain value. The similar phenomena can take place in the latter case.

Figure 4
figure 4

Effects of promoter leakage on bursting kinetics. (A and D) no feedback: promoter leakage does not change burst size nor burst frequency, where the parameter values are the same as those used in Figure 3 (A). (B and E) negative feedback: promoter leakage does not change burst frequency, but reduces burst size, where the parameter values are the same as those used in Figure 3(B). (C and F) positive feedback: promoter leakage does not change burst size but increases burst frequency, where the parameter values are the same as those used in Figure 3(C).

Figure 5
figure 5

Effects of promoter leakage on the time-dependent number of gene product molecules, where the leakage rate is set as either λ 0 = 0 (A, B) or λ 0 = 10 (D, E). (A and D) no feedback, where the parameter values are the same as those used in Figure 3(A); (B and E) negative feedback, where the parameter values are the same as those used in Figure 3(B); (C and F) positive feedback, where the leakage rate λ 1 is set as 10 for (F) or as 0 for (C), and the other parameter values are the same as those used in Figure 3(C).

Next, consider the case that there is feedback. In this case, increasing the promoter leakage rate will change the number of gene product molecules, thus in turn influencing the switching rates between the ON and the OFF states. Furthermore, this will change the dwelling times at two states of the promoter, thus influencing BS and BF. More specifically, in the case of negative feedback, increasing the leakage rate can increase the number of gene product molecules, leading that the effect of negative feedback becomes more remarkable. In other words, increasing the leakage rate will make the transition rate from the ON state to the OFF state become larger, i.e., will make the dwelling time at the ON state become shorter. This will lead to the decrease of BS, referring to Figure 4(B), and further the reduction of the expression noise. In addition, negative feedback does not change the transition rate from the OFF state to the ON state, so does not influence the dwelling time at the OFF state nor change the mean BF, referring to Figure 4(E). Figure 5(B) and (E) further justify the numerical results shown in Figure 4(B) and (E). Similarly, in the case of positive feedback, increasing the promoter leakage rate will increase the number of the gene product molecules, leading that the role of positive feedback becomes more remarkable. Specifically, increasing the leakage rate will make the transition rate from the OFF state to the ON state become larger, i.e., make the dwelling time at the OFF state become smaller. This will lead to the increase of BF, referring to Figure 4(F), and further the reduction of the expression noise. In addition, positive feedback does not change the transition rate from the ON state to the OFF state, so does not influence the dwelling time at the ON state nor change the mean BS, referring to Figure 4(C). The time series shown in Figure 5(C) and (F) further justify the numerical results shown in Figure 4(C) and (F).

Then, we want to know how increasing leakage influences bursting kinetics in the case of keeping the mean level fixed. By numerical analysis, we find that for negative feedback (λ 1λ 0), the mean BS always decreases with the increase of the leakage rate λ 0, regardless of the way that the mean is fixed, but the mean BF decreases only when the transition rate from OFF to ON (γ 1) decreases, and is kept invariant in the other three cases: increasing the transition rate from ON to OFF (γ 0); increasing the feedback strength (f); and decreasing the maximum transcription rate (λ 1), referring to Additional file 1: Figure S4 (A). In the case of positive feedback (λ 1λ 0)), if increasing the leakage rate (λ 1) corresponds to the change of the transition rate from OFF to ON (γ 0) or that of positive feedback intensity (f), then the mean BS does not change but the mean BF decreases; and if increasing the leakage rate (λ 1) corresponds to change of the transition rate from ON to OFF (γ 1) or the maximum transcription rate (λ 0), then the mean BS decreases and the mean BF increases. referring to Additional file 1: Figure S4 (B).

In a word, under the condition that the mean level is fixed, increasing the leakage rate can make the mean BS decrease in most cases and the mean BS invariant in few cases, regardless of the way the mean is fixed and the type of feedback. This implies that increasing promoter leakage tends to generate fewer bursts.

Summarizing the above analysis, we know that in the case of no feedback, promoter leakage does not influence promoter dynamics, but in the presence of feedback, the former influences the latter and this influence is through the way that feedback indirectly impacts bursting kinetics. Specifically, the promoter leakage makes the burst size be reduced and the burst frequency be enlarged, indicating that the promoter leakage makes the gene produce fewer bursts. This conclusion can be also seen from the time series shown in Figure 5. In fact, by comparing Figure 4(D), (E) and (F) with their corresponding Figure 5(D), (E) and (F), we find that fewer bursts are generated due to the effect of promoter leakage.

Discussion

Biochemical reactions associated with gene expression are all essentially single-molecule events and thus stochastic, resulting in substantial randomness in the production of mRNA or protein. This noise can significantly influence the expression levels of gene products and has been identified as a source of cell-to-cell variability. To capture effects of expression noise, many gene models have been proposed, such as those with simple promoter structures [2,7] or with a DNA loop [51] or with a more complex promoter structure [23], and those with auto regulation [7]. Almost these models, however, neglected the effect of promoter leakage on gene expression. Here, we have introduced and analyzed a stochastic model of gene expression, which considers not only promoter activity and regulation but also promoter leakage. Interestingly, we have derived the analytical distribution (seeing Eq. (2)), which can reproduce some known distributions obtained in simplified cases, e.g., with or without auto-regulation, without promoter leakage. More importantly, our results on the noise in gene product indicate that promoter leakage can be taken as a mechanism of attenuating expression noise, implying that previous gene models of no promoter leakage would overestimate expression noise. Similarly, previous estimations on bursting dynamics characterized by burst frequency and burst size would be inaccurate since they did not consider the effect of promoter leakage. In fact, we have shown that promoter leakage can significantly impact bursting dynamics (Figure 4) and expression noise (Figure 2). An intuitive interpretation for this impact is as follows. Promoter leakage can increase the amount of gene product, thus possibly reducing the intrinsic noise of gene product. Meanwhile, it also can influence promoter noise when the gene product as a transcription factor auto-regulates the transition rates between promoter states.

Regarding our model, here we present simple discussions. To simplify our analysis, our first assumption is that the binding and dissociation of a transcription factor is a very fast process, that is, it binds quickly to DNAs and dissociates rapidly from the DNAs. With this assumption, one can view that the transcription factor changes only the transition rates between promoter activity states without consuming itself. The similar assumption has been before made to obtain analytical distributions in gene models [55,56]. On the other hand, experimental studies indicated that the association and disassociation of transcription factors to promoter sites may be a slow process. Moreover, it was theoretically shown that the mode of binding of transcription factors to DNAs can affect the properties of expression noise [64,65]. If a slow process for binding of a transcription factor is introduced to the models studied here and even if transcription and translation processes are considered, then two qualitative conclusions obtained here, that is, increasing the leakage rate reduces noise and makes the gene product distribution be uni-modalized, will still held. For this, we have performed numerical simulation, with results shown in Additional file 1: Figure S1, Figure S2 and Figure S3. These results indicates that our simplified model have well captured effects of promoter leakage on expression dynamics including distribution, noise and bursting dynamics.

Our second assumption is that the gene promoter has one active state and one inactive state. In many cases, however, the promoter may have multiple activity states [23,63]. For example, the PRM promoter of phage lambda in E. coli is regulated by two different TFs binding to two sets of three operators that are brought together by looping out the intervening DNA. As a result, the number of regulatory states of the PRM promoter is up to 128 [66]. In particular, eukaryotic promoter structures would be more complex since they involve nucleosomes competing with or being removed by transcription factors [67]. In spite of this, our qualitative conclusions will not be ruined although the quantitative results would be modified if a more complex promoter structure is considered (data are not shown).

It should be pointed out that gene expression is a complex biochemical process. Except for the factors considered here such as transition between two activity state of the promoter, auto-regulation and promoter leakage, gene expression also involves other factors such as RNA nuclear retention, chromatin remodeling [4-6], combinatorial regulation from many transcription factors [68,69], DNA loop [70,71] and alternative splicing [15] as well as binding of transcription factors [64,72]. However, it is unclear how these factors impact gene expression in the combined case of feedback and promoter leakage.

Finally, our study would have biological implications in fields such as synthetic biology. First, according to our results obtained here, promoter leakage in circuits of stochastic gene expression especially in those with repressive regulation should be controlled and otherwise it would influence functions of the corresponding circuits, e.g., bistability in circuits of positive feedbacks [73], and information processing in circuits with AND gate [74]. See the review article [75] for more details. Thus, introducing promoter leakage to gene circuits to be designed would achieve better design effects. Second, our results can provide a guideline for biologists who design synthesized circuits used to probe for the relationship between bistability and phenotype. Third, our results would imply that promoter leakage is a mechanism of efficiently adjusting phenotypic diversity of in vivo organisms.

Conclusions

Promoter leakage is an unneglectable factor in gene expression and plays a significant role of controlling expression dynamics and phenotypic diversity. Specifically, promoter leakage (1) always reduces expression noise; (2) tends to unimodalize the gene-product distribution; (2) makes the gene produce fewer burst in contrast to the case of no promoter leakage. These results imply that promoter leakage may be taken as a strategy of efficiently controlling cell-to-cell variability.

Methods

Derivation of analysis results

In order to derive the analysis expression of the probability distribution and noise, we transform the chemical master equation into a confluence hypergeometric equation [51,53], using the probability-generating function.

Stochastic simulations

Stochastic simulations of our model studied here have been carried out using the Gillespie algorithm [76].

Abbreviations

BS:

Burst size

BF:

Burst frequency

CME:

Chemical master equation

References

  1. Minoru SHK. A stochastic model for gene induction. J Theor Biol. 1991;153(2):181–94.

    Article  Google Scholar 

  2. Friedman N, Cai L, Xie XS. Linking stochastic dynamics to population distribution: an analytical framework of gene expression. Phys Rev Lett. 2006;97(16):168302.

    Article  PubMed  Google Scholar 

  3. Raj A, Peskin CS, Tranchina D, Vargas DY, Tyagi S. Stochastic mRNA synthesis in mammalian cells. PLoS Biol. 2006;4(10):e309.

    Article  PubMed Central  PubMed  Google Scholar 

  4. Cosma MP, Tanaka T, Nasmyth K. Ordered recruitment of transcription and chromatin remodeling factors to a cell cycle–and developmentally regulated promoter. Cell. 1999;97(3):299–311.

    Article  CAS  PubMed  Google Scholar 

  5. Martinowich K, Hattori D, Wu H, Fouse S, He F, Hu Y, et al. DNA methylation-related chromatin remodeling in activity-dependent BDNF gene regulation. Science. 2003;302(5646):890–3.

    Article  CAS  PubMed  Google Scholar 

  6. Sanchez A, Garcia HG, Jones D, Phillips R, Kondev J. Effect of promoter architecture on the cell-to-cell variability in gene expression. PLoS Comp Biol. 2011;7(3):e1001100.

    Article  CAS  Google Scholar 

  7. Sanchez A, Choubey S, Kondev J. Regulation of noise in gene expression. Annu Rev Biophys. 2013;42:469–91.

    Article  CAS  PubMed  Google Scholar 

  8. Hornung G, Bar-Ziv R, Rosin D, Tokuriki N, Tawfik DS, Oren M, et al. Noise–mean relationship in mutated promoters. Genome Res. 2012;22(12):2409–17.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Newman JR, 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(7059):840–6.

    Article  CAS  PubMed  Google Scholar 

  10. Liu B, Yuan Z, Aihara K, Chen L. Reinitiation enhances reliable transcriptional responses in eukaryotes. JR Soc Interface. 2014;11(97):20140326.

    Article  Google Scholar 

  11. Klumpp S, Zhang Z, Hwa T. Growth rate-dependent global effects on gene expression in bacteria. Cell. 2009;139(7):1366–75.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Mileyko Y, Joh RI, Weitz JS. Small-scale copy number variation and large-scale changes in gene expression. Proc Natl Acad Sci U S A. 2008;105(43):16659–64.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. San Burger A, Walczak AM, Wolynes PG. Influence of decoys on the noise and dynamics of gene expression. Phy Rev E. 2012;86(4):041920.

    Article  Google Scholar 

  14. Burger A, Walczak AM, Wolynes PG. Abduction and asylum in the lives of transcription factors. Proc Natl Acad Sci U S A. 2010;107(9):4016–21.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Wang Q, Zhou T. Alternative-splicing-mediated gene expression. Phy Rev E. 2014;89(1):012713.

    Article  Google Scholar 

  16. Clapier CR, Cairns BR. The biology of chromatin remodeling complexes. Annu Rev Biochem. 2009;78:273–304.

    Article  CAS  PubMed  Google Scholar 

  17. Narlikar GJ, Fan HY, Kingston RE. Cooperation between complexes that regulate chromatin structure and transcription. Cell. 2002;108(4):475–87.

    Article  CAS  PubMed  Google Scholar 

  18. Venturelli OS, El-Samad H, Murray RM. Synergistic dual positive feedback loops established by molecular sequestration generate robust bimodal response. Proc Natl Acad Sci U S A. 2012;109(48):E3324–33.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Yoda M, Ushikubo T, Inoue W, Sasai M. Roles of noise in single and coupled multiple genetic oscillators. J Chem Phys. 2007;126(11):115101.

    Article  PubMed  Google Scholar 

  20. Hasty J, Isaacs F, Dolnik M, McMillen D, Collins JJ. Designer gene networks: towards fundamental cellular control. Chaos. 2001;11(1):207–20.

    Article  CAS  PubMed  Google Scholar 

  21. Pedraza JM, Paulsson J. Effects of molecular memory and bursting on fluctuations in gene expression. Science. 2008;319(5861):339–43.

    Article  CAS  PubMed  Google Scholar 

  22. Carey LB, Van Dijk D, Sloot PM, Kaandorp JA, Segal E. Promoter sequence determines the relationship between expression level and noise. PLoS Biol. 2013;11(4):e1001528.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Zhang J, Zhou T. Promoter-mediated transcriptional dynamics. Biophys J. 2014;106(2):479–88.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Marquez-Lago TT, Stelling J. Counter-intuitive stochastic behavior of simple gene circuits with negative feedback. Biophys J. 2010;98(9):1742–50.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Rosenfeld N, Elowitz MB, Alon U. Negative autoregulation speeds the response times of transcription networks. J Mol Biol. 2002;323(5):785–93.

    Article  CAS  PubMed  Google Scholar 

  26. Paulsson J. Summing up the noise in gene networks. Nature. 2004;427(6973):415–8.

    Article  CAS  PubMed  Google Scholar 

  27. Kærn M, Elston TC, Blake WJ, Collins JJ. Stochasticity in gene expression: from theories to phenotypes. Nat Rev Gen. 2005;6(6):451–64.

    Article  Google Scholar 

  28. Hung M, Chang E, Hussein R, Frazier K, Shin JE, Sagawa S, et al. Modulating the frequency and bias of stochastic switching to control phenotypic variation. Nat Commun. 2014;5:4574.

    Article  CAS  PubMed  Google Scholar 

  29. Fraser D, Kærn M. A chance at survival: gene expression noise and phenotypic diversification strategies. Mol Microbiol. 2009;71(6):1333–40.

    Article  CAS  PubMed  Google Scholar 

  30. Kussell E, Leibler S. Phenotypic diversity, population growth, and information in fluctuating environments. Science. 2005;309(5743):2075–8.

    Article  CAS  PubMed  Google Scholar 

  31. Kussell E, Kishony R, Balaban NQ, Leibler S. Bacterial persistence a model of survival in changing environments. Genetics. 2005;169(4):1807–14.

    Article  PubMed Central  PubMed  Google Scholar 

  32. Acar M, Mettetal JT, van Oudenaarden A. Stochastic switching as a survival strategy in fluctuating environments. Nat Genet. 2008;40(4):471–5.

    Article  CAS  PubMed  Google Scholar 

  33. Balaban NQ, Merrin J, Chait R, Kowalik L, Leibler S. Bacterial persistence as a phenotypic switch. Science. 2004;305(5690):1622–5.

    Article  CAS  PubMed  Google Scholar 

  34. Blake WJ, Balázsi G, Kohanski MA, Isaacs FJ, Murphy KF, Kuang Y, et al. Phenotypic consequences of promoter-mediated transcriptional noise. Mol Cell. 2006;24(6):853–65.

    Article  CAS  PubMed  Google Scholar 

  35. Booth IR. Stress and the single cell: intrapopulation diversity is a mechanism to ensure survival upon exposure to stress. Int J Food Microbiol. 2002;78(1):19–30.

    Article  PubMed  Google Scholar 

  36. Maughan H, Nicholson WL. Stochastic processes influence stationary-phase decisions in Bacillus subtilis. J Bacteriol. 2004;186(7):2212–4.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Lehner B. Selection to minimise noise in living systems and its implications for the evolution of gene expression. Mol Sys Biol. 2008;4(1):170.

    Google Scholar 

  38. Singh A, Vargas CA, Karmakar R. Stochastic analysis and inference of a two-state genetic promoter model. Washington, DC: American Control Conference (ACC); 2013. p. 4563–8.

    Google Scholar 

  39. Singh A, Vargas CA, Karmakar R. Stochastic analysis and inference of a two-state genetic promoter model. In: Proc. of the 2013 Amer. Control Conference, IEEE; 2013. P. 4563–68

  40. Cai L, Friedman N, Xie XS. Stochastic protein expression in individual cells at the single molecule level. Nature. 2006;440(7082):358–62.

    Article  CAS  PubMed  Google Scholar 

  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(5767):1600–3.

    Article  CAS  PubMed  Google Scholar 

  42. Dar RD, Razooky BS, Singh A, Trimeloni TV, McCollum JM, Cox CD, et al. Transcriptional burst frequency and burst size are equally modulated across the human genome. Proc Natl Acad Sci U S A. 2012;109(43):17454–9.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. Golding I, Paulsson J, Zawilski SM, Cox EC. Real-time kinetics of gene activity in individual bacteria. Cell. 2005;123(6):1025–36.

    Article  CAS  PubMed  Google Scholar 

  44. Chubb JR, Trcek T, Shenoy SM, Singer RH. Transcriptional pulsing of a developmental gene. Curr Biol. 2006;16(10):1018–25.

    Article  CAS  PubMed  Google Scholar 

  45. Yella VR, Bansal M. DNA structural features and architecture of promoter regions play a role in gene responsiveness of S. cerevisiae. J Bioinform Comput Biol. 2013;11(06):1343001.

    Article  PubMed  Google Scholar 

  46. Hansen AS. Promoter decoding of transcription factor dynamics involves a trade‐off between noise and control of gene expression. Mol Sys Biol. 2013;9(1):704.

    Article  CAS  Google Scholar 

  47. Baler R, Covington S, Klein DC. Rat arylalkylamine N-acetyltransferase gene: upstream and intronic components of a bipartite promoter. Biol Cell. 1999;91(9):699–705.

    Article  CAS  PubMed  Google Scholar 

  48. Reed CC, Waterhouse A, Kirby S, Kay P, Owens RT, McQuillan DJ, et al. Decorin prevents metastatic spreading of breast cancer. Oncogene. 2004;24(6):1104–10.

    Article  Google Scholar 

  49. McAdams HH, Arkin A. Stochastic mechanisms in gene expression. Proc Natl Acad Sci U S A. 1997;94(3):814–9.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Raser JM, O’Shea EK. Control of stochasticity in eukaryotic gene expression. Science. 2004;304(5678):1811–4.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Zhou T, Zhang J. Analytical results for a multistate gene model. SIAM J Appl Math. 2012;72(3):789–818.

    Article  Google Scholar 

  52. Shahrezaei V, Swain PS. Analytical distributions for stochastic gene expression. Proc Natl Acad Sci U S A. 2008;105(45):17256–61.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  53. Grima R, Schmidt DR, Newman TJ. Steady-state fluctuations of a genetic feedback loop: an exact solution. J Chem Phys. 2012;137(3):035104.

    Article  CAS  PubMed  Google Scholar 

  54. Grönlund A, Lötstedt P, Elf J. Transcription factor binding kinetics constrain noise suppression via negative feedback. Nat Commun. 2013;4:1864.

    Article  PubMed  Google Scholar 

  55. Karpova TS, Kim MJ, Spriet C, Nalley K, Stasevich TJ, Kherrouche Z, et al. Concurrent fast and slow cycling of a transcriptional activator at an endogenous promoter. Science. 2008;319(5862):466–9.

    Article  CAS  PubMed  Google Scholar 

  56. Blossey R, Giuraniuc CV. Mean-field versus stochastic models for transcriptional regulation. Phy Rev E. 2008;78(3):031909.

    Article  CAS  Google Scholar 

  57. Paulsson J. Models of stochastic gene expression. Phys Life Rev. 2005;2(2):157–75.

    Article  Google Scholar 

  58. Ozbudak EM, Thattai M, Kurtser I, Grossman AD, van Oudenaarden A. Regulation of noise in the expression of a single gene. Nat Genet. 2002;31(1):69–73.

    Article  CAS  PubMed  Google Scholar 

  59. Schwabe A, Rybakova KN, Bruggeman FJ. Transcription stochasticity of complex gene regulation models. Biophys J. 2012;103(6):1152–61.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  60. To TL, Maheshri N. Noise can induce bimodality in positive transcriptional feedback loops without bistability. Science. 2010;327(5969):1142–5.

    Article  CAS  PubMed  Google Scholar 

  61. Weinberger LS, Shenk T. An HIV feedback resistor: auto-regulatory circuit deactivator and noise buffer. PLoS Biol. 2006;5(1):e9.

    Article  PubMed Central  Google Scholar 

  62. Sanchez A, Choubey S, Kondev J. Stochastic models of transcription: from single molecules to single cells. Methods. 2013;62(1):13–25.

    Article  CAS  PubMed  Google Scholar 

  63. Halme A, Bumgarner S, Styles C, Fink GR. Genetic and epigenetic regulation of the < i > FLO</i > gene family generates cell-surface variation in yeast. Cell. 2004;116(3):405–15.

    Article  CAS  PubMed  Google Scholar 

  64. Vilar JMG, Saiz L. Systems biophysics of gene expression. Biophys J. 2013;104(12):2574–85.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  65. Saiz L, Vilar JMG. Stochastic dynamics of macromolecular-assembly networks. Mol Sys Biol. 2006;2:0024.

    Google Scholar 

  66. Vilar JMG, Saiz L. CplexA: a mathematica package to study macromolecular-assembly control of gene expression. Bioinformatics. 2010;26(16):2060–1.

    Article  CAS  PubMed  Google Scholar 

  67. Boeger H, Griesenbeck J, Kornberg RD. Nucleosome retention and the stochastic nature of promoter chromatin remodeling for transcription. Cell. 2008;133(4):716–26.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  68. Sanchez A, Kondev J. Transcriptional control of noise in gene expression. Proc Natl Acad Sci U S A. 2008;105(13):5081–6.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  69. William JB, Mads K, Charles RC, Collins JJ. Noise in eukaryotic gene expression. Nature. 2003;422(01546):633–7.

    Google Scholar 

  70. Vilar JMG, Saiz L. Suppression and enhancement of transcriptional noise by DNA looping. Phy Rev E. 2014;89(6):062703.

    Article  Google Scholar 

  71. Hensel Z, Weng X, Lagda AC, Xiao J. Transcription-factor-mediated DNA looping probed by high-resolution, single-molecule imaging in live E. Coli cells. PLoS Biol. 2013;11(6):e1001591.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  72. Stekel DJ, Jenkins DJ. Strong negative self regulation of prokaryotic transcription factors increases the intrinsic noise of protein expression. BMC Syst Biol. 2008;2:6.

    Article  PubMed Central  PubMed  Google Scholar 

  73. Becskei A, Séraphin B, Serrano L. Positive feedback in eukaryotic gene networks cell differentiation by graded to binary response conversion. EMBO J. 2001;20(10):2528–35.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  74. Minaba M, Kato Y. High-yield, zero-leakage expression system with a translational switch using site-specific unnatural amino acid incorporation. Appl Environ Microbiol. 2014;80(5):1718–25.

    Article  PubMed Central  PubMed  Google Scholar 

  75. Chen YY, Galloway KE, Smolke CD. Synthetic biology: advancing biological frontiers by building synthetic systems. Genome Biol. 2012;13(2):240.

    Article  PubMed Central  PubMed  Google Scholar 

  76. Gillespie DT. Exact stochastic simulation of coupled chemical reactions. J Chem Phys. 1977;81(25):2340–61.

    Article  CAS  Google Scholar 

Download references

Acknowledgments

This work was supported by Science and Technology Department No. 2014CB964703 and No. 2010CB945400, and the Natural Science Foundation No. 91230204,No. 11475273, and the New Star of Science and Technology Project of Zhujiang No.2012J2200017, and Specialized Research Fund for the Doctoral Program of Higher Education No. 20110171120045, and the Natural Science Foundation of Hunan Province No.14JJ3134,and the Natural Science Foundation of Hainan Province No.114002 of P. R. China.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Tianshou Zhou.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

LH, ZY and TZ designed the study. LH performed the theoretical derivations of the model and conceived/analysed the qualitative results. LH, ZY, PL and TZ interpreted the results and wrote the paper. All authors read and approved the final manuscript.

Additional files

Additional file 1:

Effects of promoter leakage on expression dynamics in more complex situations. Figure S1. The effect of promoter leakage on expression dynamics in a gene auto-repressing model that assumes the dissociation of transcription factors from genomic binding site is a slow process. Figure S2. The effect of promoter leakage on expression dynamics in a gene auto-enhancing model that assumes the dissociation of transcription factors from genomic binding site is a slow process. Figure S3. Effects of promoter leakage on distribution and noise in a full two-state gene model. Figure S4. The effect of promoter leakage on burst dynamics in the case that amount of the gene product is fixed.

Additional file 2:

Theoretical derivations of probability and noise.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Huang, L., Yuan, Z., Liu, P. et al. Effects of promoter leakage on dynamics of gene expression. BMC Syst Biol 9, 16 (2015). https://0-doi-org.brum.beds.ac.uk/10.1186/s12918-015-0157-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12918-015-0157-z

Keywords