Pith. sign in

REVIEW 5 major objections 5 minor 43 references

Effects of restrained degradation on gene expression and regulation

T0 review · 5 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Restrained degradation may explain why gene-expression noise slightly exceeds Poisson statistics.

desk verdict New modeling idea with a genuinely broken central normalization: the Fano-factor claim is undefined as written, but the qualitative picture is plausible and fixable. read the letter →

arxiv 1908.10511 v1 pith:5UCP3TQX submitted 2019-08-28 physics.bio-ph

Yan-Ling Feng , Jing Sun , Yi-Fan Liu , Jian-Gong Ren , Jian-Min Dong This is my paper · ORCID
classification physics.bio-ph MSC 60H1082C3192C40 PACS 87.17.-d87.57.cm05.40.-a
keywords LangevinequationFokker-PlanckgeneexpressionnoiseFanofactorcarryingcapacitydegradationsaturationstochasticresonancemeanfirstpassagetime
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper argues that the degradation of mRNA and proteins can be restrained by a finite carrying capacity of the cellular environment—analogous to logistic growth—and that this 'K effect' changes gene-expression statistics. For constitutive expression, the stationary mRNA distribution gains a cubic term in its exponent, exp(γx³/(3KQ)), making the Fano factor slightly larger than 1 under intrinsic noise alone. The paper proposes that this mechanism, not external noise or time delays, is the probable origin of the slight Fano-factor deviation measured in recent experiments. In a genetic regulatory circuit, the same effect shifts probability toward high transcription-factor levels, can abolish bistability, slows switching between states, and reduces the signal-to-noise ratio of stochastic resonance.

What carries the argument

The load-bearing object is the logistic degradation rate γx(1 − x/K) with carrying capacity K, inserted into a Langevin equation with multiplicative external noise and additive internal noise. The stationary solution of the Fokker-Planck equation then carries a term exp(γx³/(3KQ)) in the intrinsic-noise-only limit; this cubic exponent is what tilts the distribution and pushes the Fano factor above 1. For the regulatory circuit, the same machinery enters the modified potential Φ(x) = −∫ A(x)/G²(x) dx used to compute mean first passage times and the stochastic-resonance signal-to-noise ratio.

What would settle it

Measure single-cell mRNA distributions from a constitutive promoter while titrating the level of degrading enzymes: the K effect predicts Fano factor > 1 that grows as the degradation capacity shrinks, whereas the Poisson null predicts Fano = 1 at every condition. Alternatively, fit Eq. (5) to published single-cell data; if the best-fit K is not positive and finite, or if the Fano factor stays at 1 for all mean copy numbers once known extrinsic noise is subtracted, the K-effect explanation is falsified.

Watch

Extended reading notes

Core claim

The central claim is that replacing the linear degradation rate γx with the logistic form γx(1 − x/K) changes the stationary distribution of constitutive gene expression from a Gaussian (the large-copy-number limit of Poisson) into P(x) = N exp(−(x − α/γ)²/(2Q/γ)) exp(γx³/(3KQ)). The extra factor is positive and grows with x, skewing the distribution to the right and inflating the variance relative to the mean, so the Fano factor becomes slightly greater than 1 for finite K and returns to 1 as K → ∞. On this basis the paper attributes the experimentally observed slight deviation of the Fano factor from 1 (Science 346, 1533 (2014)) to the K effect. For the Smolen-type regulatory system, the K effect shifts the bimodal distribution toward the high state, and for strong restraint (small K) the bistability disappears; the mean first passage time increases as K decreases (especially at weak noise), and the stochastic-resonance signal-to-noise ratio is reduced with its peak shifted to larger noise intensities.

Load-bearing premise

The load-bearing assumption is that degradation really follows the logistic law γx(1 − x/K) and that mRNA copy numbers never reach or exceed K; the paper gives no independent measurement of K, and without restricting the state space to x < K the stationary distribution in Eq. (5) cannot be normalized.

Editorial extensions

If this is right

  • If the K effect is real, the Fano factor of a constitutive gene becomes a quantitative readout of degradation capacity: fitting the distribution to Eq. (5) yields an estimate of K.
  • Reducing catabolic enzyme or RNase levels in a cell should raise the Fano factor monotonically, a direct prediction for perturbation experiments.
  • In gene-regulatory switches, saturating degradation slows the transition from the high state to the low state and can eliminate bistability entirely when K is small.
  • Periodic forcing of a switch becomes less effective as K decreases: the stochastic-resonance peak is suppressed and displaced toward stronger noise.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The distribution in Eq. (5) is normalizable only if x is restricted to x < K; a refined model would impose an absorbing or reflecting boundary at K or replace logistic saturation with a Michaelis–Menten form that stays positive for all x.
  • A testable extension is to compare Eq. (5) against alternative super-Poissonian models (e.g., negative-binomial or Gamma distributions) on the same single-cell data; the cubic-exponent shape is distinctive and would show up in the skewness of the mRNA distribution.
  • If the K effect is confirmed, it would connect two currently separate literatures: noise in gene expression and saturating enzymatic degradation of mRNA, suggesting that resource limitation is not just a metabolic constraint but a sculpting force on expression noise.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

5 major / 5 minor

Summary. The paper proposes a logistic-type modification of the mRNA/protein degradation rate, replacing γx by γx(1−x/K), where K is interpreted as the environmental carrying capacity for degradation. For constitutive gene expression, the authors write a Langevin equation with intrinsic and extrinsic Gaussian white noises, derive the corresponding Fokker-Planck equation, and present a stationary distribution. They claim that for intrinsic noise alone the Fano factor becomes slightly larger than 1 for finite K, and suggest that this K effect may explain the Fano factor deviation from 1 measured in Science 346, 1533 (2014). For a Smolen-type genetic regulatory circuit, they similarly derive a stationary distribution and compute the mean first passage time (MFPT) and signal-to-noise ratio (SNR) for stochastic resonance, reporting that a stronger K effect increases the MFPT and reduces the SNR.

Significance. If the central claim holds, the paper offers a simple, analytically tractable mechanism for the observed Fano factor above 1 in constitutive gene expression, and it extends the analysis to transient properties of a regulatory switch. The Fano factor is computed rather than fitted, and the MFPT/SNR results are additional derived outputs. However, the significance is limited by two issues. First, the logistic degradation form γx(1−x/K) is introduced ad hoc, and for x>K it turns degradation into production, which is physically questionable for mRNA decay; the authors provide no direct biological justification or independent measurement of K. Second, the central Fano-factor claim rests on the stationary density Eq. (5), which is not normalizable on the stated domain [0,∞) for finite K, so without an explicit support restriction the claim is not well-defined. The paper thus has a promising idea but needs substantial revision before the central claim can be evaluated.

major comments (5)
  1. [Sec. II, Eq. (5)] The stationary density in Eq. (5) is not normalizable on [0,∞) for finite K because the factor exp(γx^3/(3KQ)) diverges as x→∞. The paper never states a support restriction. This is not a cosmetic detail: the deterministic drift in Eq. (1) has A(K)=α>0, so at x=K the drift is outward, and for x>K the degradation term γx(1−x/K) becomes negative, i.e., it acts as production rather than degradation. Without an explicit reflecting boundary at x=K (or another cutoff) and a corresponding boundary condition for the Fokker-Planck equation, the stationary distribution, the Fano factor in Fig. 2, and the results in Sec. III are not well-defined. The authors must specify the domain, impose a boundary condition, and recompute the normalization and the Fano factor under that boundary condition.
  2. [Sec. II, Eq. (4)] For the Fokker-Planck equation (3), the stationary solution is P(x) ∝ G^{-2}(x) exp(∫ A(x)/G^2(x) dx), but Eq. (4) uses G^{-1}(x), not G^{-2}(x). This is an error in the stationary solution. In the D=0 limit G=√Q is constant, so the missing factor is absorbed into the normalization constant and Eq. (5) is unaffected, but all D>0 stationary distributions (Fig. 1 upper panel and Fig. 3) and their normalization are affected. The prefactor should be corrected.
  3. [Sec. II, Eq. (1)] The logistic degradation form γx(1−x/K) is introduced without a physical justification for mRNA or protein decay. For x>K this term becomes negative, meaning the degradation rate is actually a production rate. This behavior is unphysical for molecular degradation. If the intended effect is saturation of the degradation rate, a more natural form would be γx/(1+x/K) (Michaelis-Menten-like); the authors should either provide a biological rationale for the sign change or test whether their conclusions survive under a saturating degradation form.
  4. [Sec. II, Fig. 2 and discussion of Science 346, 1533 (2014)] The central claim that the observed Fano factor slightly larger than 1 'probably originates from the K effect' is not supported by a quantitative comparison with the experimental data. The parameter K is a free parameter and is never measured or constrained. Since different values of K produce different Fano factors, the explanation is not falsifiable unless the authors show that the K values required to match the Science 2014 data are biologically plausible. A direct comparison to the experimental Fano-factor data, with a discussion of the inferred K, is needed.
  5. [Sec. II, Eqs. (2) and (3)] The paper does not specify the stochastic interpretation of the multiplicative noise term x(1−x/K)ξ(t) in Eq. (2). The Fokker-Planck equation (3) corresponds to the Ito interpretation; under Stratonovich, an additional drift term D g(x)g'(x) appears. The stationary distribution for D>0 depends on this choice. The authors should state which convention they use and, if relevant, justify it for the biological noise source.
minor comments (5)
  1. [Sec. II, text near Eq. (5)] The notation 'the variance σ^2 is equal to the mean x' uses x for both the random variable and its expectation; please use ⟨x⟩ for the mean.
  2. [Sec. III, Fig. 3 caption] The caption contains the duplicated phrase 'different different K values'; please correct this typo.
  3. [Fig. 1 caption] The caption contains 'K effct' instead of 'K effect', and the axis labels in the figures appear garbled; please ensure all figure text is legible.
  4. [Secs. II and III, Eqs. (4) and (7)] The normalization constants N0 and N are not defined; please state explicitly that they are determined by numerical integration, or provide their analytic expressions.
  5. [Sec. III, Eq. (8)] The mean first passage time formula uses U0''(xu) and U0''(x2), but the potential U0(x) is not defined in the paper; please define it and its relation to the drift A(x).

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the Fano-factor and MFPT/SNR results are computed outputs of the stated model, with self-citations used only for comparison and parameters, not as load-bearing justification.

full rationale

The paper's core derivation chain is self-contained: Eq. (1) posits a logistic degradation term, Eq. (2) adds white noises, Eq. (3) writes the corresponding Fokker-Planck equation, Eq. (4) gives the stationary solution, and Eq. (5) is the D=0 specialization used for the Fano factor. The Fano factor is then computed from Eq. (5) and shown in Fig. 2 as a function of K and mean mRNA copy number. No parameter is fitted to the Science 2014 experimental Fano-factor data that the paper discusses; alpha and gamma are assigned values, K is scanned, and the experimental comparison is qualitative and post hoc. The MFPT and SNR results in Sec. III are additional outputs of the same model and do not reduce to any fitted input. The self-citations [21], [30], and [31] are prior work by the authors on related gene-regulation systems and time-delay effects; they are used for comparison and parameter context, not to justify the central claim. One mathematical concern is that Eq. (5) contains exp(+gamma x^3/(3KQ)), which is not normalizable on [0,infinity) for finite K unless an implicit x in [0,K] domain and boundary conditions are supplied; this is a correctness and well-posedness flaw, but it is not circularity because the output is still genuinely computed from the stated model rather than being equivalent to the input by construction.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

The central claim rests on an ad hoc logistic degradation model whose state space is not specified, a noise interpretation that is not stated, and standard FPE/Kramers/SNR results from the literature. The only free parameter specific to the claimed new effect is K, which is varied rather than measured.

free parameters (3)
  • K = varied (e.g., K=50, 200, 500 in figures)
    Carrying capacity of environment for degradation; central parameter of the model, not measured independently.
  • Q = set to alpha=20 for constitutive expression
    Intrinsic noise amplitude calibrated so that the K->infinity limit yields Fano factor 1; no microscopic derivation.
  • D = varied (e.g., D=0.01, 0.02)
    External noise intensity; varied to explore distribution shapes and stochastic resonance.
assumptions (5)
  • ad hoc to paper Degradation rate follows the logistic form gamma*x*(1 - x/K).
    Eq. (1) replaces linear degradation with a logistic term without derivation from enzyme kinetics; for x>K the term becomes negative, acting as production.
  • ad hoc to paper The state space is restricted so that x<K (or some cutoff) to make distributions normalizable.
    Eq. (5) is not normalizable on [0, infinity) for finite K; the paper never states the implied domain restriction.
  • domain assumption The Langevin equation with multiplicative noise is interpreted in the Ito sense.
    The Fokker-Planck equation Eq. (3) has drift A(x) and diffusion G^2(x), which corresponds to Ito convention; the paper does not specify this.
  • domain assumption Internal and external noises are independent Gaussian white noises.
    Eq. (2) assumes delta-correlated Gaussian noises with zero mean and independence; standard in the gene expression noise literature.
  • standard math Kramers MFPT formula and McNamara-Wiesenfeld SNR formula apply.
    Eqs. (8)-(12) are taken from Refs. [32,33,43] and assume a bistable potential and adiabatic limit.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Effects of restrained degradation on gene expression and regulation." pith.science (2026). https://pith.science/paper/5UCP3TQX

@misc{pith2026190810511,
  author       = {Pith},
  title        = {Pith review of: Effects of restrained degradation on gene expression and regulation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5UCP3TQX}},
  note         = {Machine review of arXiv:1908.10511}
}
abstract

The effects of carrying capacity of environment $K$ for degradation (the $K$ effect for short) on the constitutive gene expression and a simple genetic regulation system, are investigated by employing a stochastic Langevin equation combined with the corresponding Fokker-Planck equation for the two stochastic systems subjected to internal and external noises. This $K$ effect characterizes the limited degradation ability of the environment for RNA or proteins, such as insufficient catabolic enzymes. The $K$ effect could significantly change the distribution of mRNA copy-number in constitutive gene expression, and interestingly, it leads to the Fano factor slightly larger than 1 if only the internal noise exists. Therefore, that the recent experimental measurements suggests the Fano factor deviates from 1 slightly (Science {\bf 346} (2014) 1533), probably originates from the $K$ effect. The $K$ effects on the steady and transient properties of genetic regulation system, have been investigated in detail. It could enhance the mean first passage time significantly especially when the noises are weak and reduce the signal-to-noise ratio in stochastic resonance substantially.

Figures

Figures reproduced from arXiv: 1908.10511 by the authors.

Figure 1
Figure 1. FIG. 1: The stationary probability distribution [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: The Fano factor versus the (a) carrying capacities of [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: The stationary probability distribution [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: FIG. 4: The mean first passage time versus the noise strength [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5: The signal-to-noise ratio (SNR) versus the external [PITH_FULL_IMAGE:figures/full_fig_p007_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

43 extracted references · 42 canonical work pages

  1. [1]

    Becskei, L

    A. Becskei, L. Serrano, Nature 405, 590 (2000)

    work page 2000

  2. [2]

    Ahmad and S

    K. Ahmad and S. Heniko ff, Cell 104, 839 (2001)

    work page 2001

  3. [3]

    Elowitz, A

    M. Elowitz, A. Levine, E. Siggia, and P . Swain, Science 297, 1183 (2002)

    work page 2002

  4. [4]

    J. M. Pedraza and A. van Oudenaarden, Science 307, 1965 (2005)

    work page 2005

  5. [5]

    Bratsun, D

    D. Bratsun, D. V olfson, L. S. Tsimring, and J. Hasty, PNAS 102, 14593 (2005)

    work page 2005

  6. [6]

    Sanchez and J

    A. Sanchez and J. Kondev, PNAS 105, 5081 (2008)

    work page 2008

  7. [7]

    L.-h. So, A. Ghosh, C. Zong, L. A Sepulveda, R. Segev, I. Go lding, Nat. Gene. 43, 554 (2011)

    work page 2011

  8. [8]

    Munsky, G

    B. Munsky, G. Neuert, and A. van Oudenaarden, Science 336, 183 (2012)

    work page 2012

Show all 43 references
  1. [9]

    L. S. Tsimring, Rep. Prog. Phys. 77, 026601 (2014)

    2014

  2. [10]

    Anandamohan, J

    G. Anandamohan, J. Biol. Phys. 41, 49 (2015)

    2015

  3. [11]

    D. T. Gillespie, J. Phys. Chem. 81, 2340 (1977)

    1977

  4. [12]

    M. A. Gibson and J. Bruck, J. Phys. Chem. 104, 1876 (2000)

    2000

  5. [13]

    D. T. Gillespie, J. Chem. Phys. 115, 1716 (2001)

    2001

  6. [14]

    Kepler and T

    T. Kepler and T. Elston, Biophys. J. 81, 3116 (2001)

    2001

  7. [15]

    Tian, et al., J

    T. Tian, et al., J. Comput. Appl. Math. 205, 696 (2007)

    2007

  8. [16]

    Bretta and T

    T. Bretta and T. Gallab, J. Chem. Phys. 140, 124112 (2014)

    2014

  9. [17]

    W. R. Zhong, Y . Z. Shao, Z. H. He, Phys. Rev. E 73, 060902(R) (2006)

    2006

  10. [18]

    C. J. Wang, D. Li, D. C. Mei, Commun. Theor. Phys. 52, 463 (2009)

    2009

  11. [19]

    L. C. Du, D. C. Mei, Phys. Lett. A 374, 3275 (2010)

    2010

  12. [20]

    Smolen, D

    P . Smolen, D. A. Baxter and J. H. Byrne, Am. J. Physiol 274, 531 (1998)

    1998

  13. [21]

    Feng, J.-M

    Y .-L. Feng, J.-M. Dong, D. Wang, and X.-L. Tang, Commun. Theor. Phys. 68, 357 (2017)

    2017

  14. [22]

    D. L. Jones, R. C. Brewster, R. Phillips, Science 346, 1533 (2014)

    2014

  15. [23]

    Liu and Y

    Q. Liu and Y . Jia, Phys. Rev. E 70, 041907 (2004)

    2004

  16. [24]

    C. H. Zeng and X. C. Wei, Chin. Phys. Lett. 25, 1587 (2008)

    2008

  17. [25]

    C. H. Zeng and C. W. Xie, Phys. Scr. 78, 035801 (2008)

    2008

  18. [26]

    C. J. Wang, Chin. Phys. B 19, 030503 (2010)

    2010

  19. [27]

    C. Y . Bai, Y . Yan, and D. C. Mei, Chin. Phys. B 19, 060503 (2010). 10

    2010

  20. [28]

    C. J. Wang, Acta Phys. Sin. 61, 010503 (2012)

    2012

  21. [29]

    Yang, et al., J

    T. Yang, et al., J. Stat. Mech. 12, 12015 (2014)

    2014

  22. [30]

    Feng, J.-M

    Y .-L. Feng, J.-M. Dong, and X.-L. Tang, Chin. Phys. Lett . 33, 108701 (2016)

    2016

  23. [31]

    Y . L. Feng, J. Zhu, M. Zhang, L. L. Gao, Y . F. Liu, and J. M. D ong, Int. J. Mod. Phys. B 30, 1650067 (2016)

    2016

  24. [32]

    Jia and J

    Y . Jia and J. R. Li, Phys. Rev. E 53, 5764 (1996)

    1996

  25. [33]

    D. C. Mei, G. Z. Xie, L. Cao, and D. J. Wu, Phys. Rev. E 59, 3880 (1999)

    1999

  26. [34]

    Benzi, A

    R. Benzi, A. Sutera, and A. Vulpiani, J. Phys. A 14, 453 (1981)

    1981

  27. [35]

    Nicolis and G

    C. Nicolis and G. Nicolis, Tellus 33, 225 (1981)

    1981

  28. [36]

    Y . R. Zhou, Chin. Phys. B 20, 010501 (2011)

    2011

  29. [37]

    K. K. Wang and X. B. Liu, Chin. Phys. Lett. 30, 70504 (2013)

    2013

  30. [38]

    L. F. Lin, Y . Tian, and H. Ma, Chin. Phys. B 23, 080503 (2014)

    2014

  31. [39]

    K. K. Wang and X. B. Liu, Chin. Phys. B 23, 010502 (2014)

    2014

  32. [40]

    J. H. Li, Chin. Phys. Lett. 31, 030502 (2014)

    2014

  33. [41]

    D. X. Yang, F. S. Gu, G. J. Feng, Y . M. Yang, and B. Andrew, C hin. Phys. B 24, 110502 (2015)

    2015

  34. [42]

    Y . X. Wang, J. Q. Zhai, W. W. Xu, G. Z. Sun, and P . H. Wu, Chin. Phys. Lett. 32, 097401 (2015)

    2015

  35. [43]

    McNamara and K

    B. McNamara and K. Wiesenfeld, Phys. Rev. A 39, 4854 (1989)

    1989

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.