Kinetic theories of state- and generation-dependent cell populations
Pith reviewed 2026-05-24 03:26 UTC · model grok-4.3
The pith
A kinetic theory describes internal cell states and generation numbers in stochastically dividing populations.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We formulate a general, high-dimensional kinetic theory describing the internal state of cells in a stochastically evolving population. The resolution of our kinetic theory also allows one to track subpopulations associated with each generation. Both intrinsic noise of the cell's internal attribute and randomness in a cell's division times are fundamental to the development of our model. Based on this general framework, we are able to marginalize the high-dimensional kinetic PDEs in a number of different ways to derive equations that describe the dynamics of marginalized or macroscopic quantities such as structured population densities, moments of generation-dependent cellular states, andmom
What carries the argument
High-dimensional kinetic partial differential equations that couple stochastic evolution of cell internal states with birth-death processes indexed by generation.
If this is right
- Structured population densities can be obtained by marginalizing the high-dimensional kinetic PDEs.
- Moments of generation-dependent cellular states follow from the same marginalization procedure.
- Moments of the total population size are likewise derivable.
- Nonlinear interaction terms appear in lower-dimensional integrodifferential equations when cell rates depend on variables of other cells.
- The framework resolves the coevolution of cell populations and cell states for applications such as gene expression in developing tissues.
Where Pith is reading between the lines
- The same marginalization steps could be applied to models that add spatial position as an additional cell attribute.
- Lineage-specific accumulation of mutations might be tracked by extending the generation index to include state-dependent mutation rates.
- Comparison against individual-based simulations would test whether the reduced equations preserve fluctuation statistics for small populations.
- Feedback from total population density back to individual rates could be inserted at the high-dimensional level before marginalization.
Load-bearing premise
The high-dimensional kinetic PDEs can be marginalized in multiple ways to obtain closed or useful equations for structured population densities, moments, and total population without introducing uncontrolled approximations or losing essential dynamical features.
What would settle it
Direct stochastic simulation of a cell population with state-dependent division rates where the derived marginalized moment equations deviate measurably from the full-model statistics would show that the marginalization step loses essential features.
Figures
read the original abstract
We formulate a general, high-dimensional kinetic theory describing the internal state (such as gene expression or protein levels) of cells in a stochastically evolving population. The resolution of our kinetic theory also allows one to track subpopulations associated with each generation. Both intrinsic noise of the cell's internal attribute and randomness in a cell's division times (demographic stochasticity) are fundamental to the development of our model. Based on this general framework, we are able to marginalize the high-dimensional kinetic PDEs in a number of different ways to derive equations that describe the dynamics of marginalized or "macroscopic" quantities such as structured population densities, moments of generation-dependent cellular states, and moments of the total population. We also show how nonlinear "interaction" terms in lower-dimensional integrodifferential equations can arise from high-dimensional linear kinetic models that contain rate parameters of a cell (birth and death rates) that depend on variables associated with other cells, generating couplings in the dynamics. Our analysis provides a general, more complete mathematical framework that resolves the coevolution of cell populations and cell states. The approach may be tailored for studying, e.g., gene expression in developing tissues, or other more general particle systems which exhibit Brownian noise in individual attributes and population-level demographic noise.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript formulates a high-dimensional linear kinetic PDE framework for stochastically evolving cell populations that tracks both internal cell states (e.g., gene expression levels) and generation-specific subpopulations. It incorporates intrinsic Brownian noise in cell attributes and demographic stochasticity in division times, then derives marginalized equations for structured population densities, generation-dependent moments, and total population size. The work also shows how nonlinear interaction terms can emerge in lower-dimensional integrodifferential equations from the underlying linear high-dimensional model when birth/death rates depend on variables associated with other cells.
Significance. If the marginalizations are rigorously shown to be exact (or controlled) without loss of essential cross-generation correlations or introduction of unclosed hierarchies, the framework would supply a general, parameter-free route from microscopic stochastic assumptions to macroscopic structured-population equations. This is potentially useful for modeling coevolution of cell states and demographics in tissues or other particle systems with both intrinsic and demographic noise; the explicit separation of linear high-dimensional dynamics from emergent nonlinearities is a conceptual strength.
major comments (2)
- [Abstract and marginalization sections] The central claim that multiple marginalizations over internal states and generations produce closed or useful macroscopic equations without uncontrolled approximations is load-bearing but not yet verified in the supplied abstract; the skeptic concern that marginalization over state-dependent division times typically yields integro-differential or unclosed moment equations must be addressed by explicit derivation steps (e.g., in the section presenting the marginalization procedure).
- [Section deriving nonlinear terms] It is unclear whether the emergence of nonlinear interaction terms preserves all cross-generation correlations or requires additional closure assumptions when rates depend on other cells' states; this needs to be shown with a concrete example (e.g., a two-cell or low-dimensional case) that compares the marginalized equation to the full high-dimensional solution.
minor comments (2)
- [Model formulation] Notation for the high-dimensional state variable and generation index should be introduced with a clear table or list of symbols to aid readability.
- [Abstract] The abstract states that the approach 'resolves the coevolution' but does not specify the precise sense in which the marginal equations are closed or exact; a short clarifying sentence would help.
Simulated Author's Rebuttal
We thank the referee for their constructive and detailed comments, which help clarify the presentation of our kinetic framework. We address each major comment below and will revise the manuscript to strengthen the exposition of the marginalization procedures and the emergence of nonlinear terms.
read point-by-point responses
-
Referee: [Abstract and marginalization sections] The central claim that multiple marginalizations over internal states and generations produce closed or useful macroscopic equations without uncontrolled approximations is load-bearing but not yet verified in the supplied abstract; the skeptic concern that marginalization over state-dependent division times typically yields integro-differential or unclosed moment equations must be addressed by explicit derivation steps (e.g., in the section presenting the marginalization procedure).
Authors: The full manuscript contains explicit step-by-step derivations of the marginalizations (in the sections developing the high-dimensional kinetic PDE and the subsequent marginalization procedures) that demonstrate closure for the structured densities, generation-dependent moments, and total population size. These derivations rely on the linearity of the underlying model and the specific stochastic assumptions, which permit exact integration even when division times depend on internal states, avoiding unclosed hierarchies. The abstract summarizes the outcomes rather than the derivations. We will revise the abstract to note the rigorous, approximation-free nature of the marginalizations and expand the marginalization section with additional explicit derivation steps to address the skeptic concern directly. revision: yes
-
Referee: [Section deriving nonlinear terms] It is unclear whether the emergence of nonlinear interaction terms preserves all cross-generation correlations or requires additional closure assumptions when rates depend on other cells' states; this needs to be shown with a concrete example (e.g., a two-cell or low-dimensional case) that compares the marginalized equation to the full high-dimensional solution.
Authors: The nonlinear interaction terms arise directly from embedding state-dependent rates (coupled across cells) into the otherwise linear high-dimensional kinetic model; the marginalization remains exact and preserves cross-generation correlations because no additional closures are imposed beyond the original linear stochastic dynamics. The lower-dimensional equations inherit the correlations from the full system. To make this transparent, we will add a concrete low-dimensional example (e.g., a two-generation case with explicit state-dependent birth rates) to the section on nonlinear terms, explicitly comparing the marginalized integrodifferential equation against the solution of the corresponding high-dimensional linear system. revision: yes
Circularity Check
No significant circularity; derivation is forward from stochastic assumptions
full rationale
The paper starts from explicit stochastic assumptions on intrinsic cell-state noise and demographic division-time randomness, constructs a high-dimensional linear kinetic PDE, and then performs marginalizations over states or generations to obtain structured densities, moments, and total-population equations. No quoted step reduces a derived macroscopic equation to a fitted parameter or to a self-citation whose content is itself the target result; the abstract and provided text present the marginalizations as exact consequences of the linear high-dimensional model under the stated assumptions. Self-citation load-bearing, ansatz smuggling, or renaming of known results are absent from the given material. The central claim therefore remains a self-contained forward derivation rather than a tautological re-expression of its inputs.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Cell internal states evolve according to stochastic processes and division times are random (demographic stochasticity).
- domain assumption High-dimensional kinetic PDEs admit useful marginalizations to lower-dimensional equations for densities and moments.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We formulate a general, high-dimensional kinetic theory describing the internal state ... marginalize the high-dimensional kinetic PDEs in a number of different ways to derive equations that describe ... structured population densities, moments ... nonlinear 'interaction' terms ... from high-dimensional linear kinetic models
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Both intrinsic noise of the cell's internal attribute and randomness in a cell's division times (demographic stochasticity) are fundamental ... SDE-jump-process hybrid model
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
or protein [14, 15] levels from experimental data, mathematica l models of how expression levels or cell states evolve is often couched in terms of transport along Waddington or fitness landscapes [16, 17]. The value of the landscape may represent an “energy” function that is shaped by different ge nes, or a proliferation rate that is a function different ge...
-
[2]
Some remarks on changing population s
Heinz von Foerster. Some remarks on changing population s. In Jr. F. Stohlman, editor, The Kinetics of Cellular Prolifer- ation, pages 382–407. Grune and Stratton, New York, 1959
work page 1959
-
[3]
Yue Wang, Renaud Dessalles, and Tom Chou. Modelling the i mpact of birth control policies on China’s population and age: effects of delayed births and minimum birth age constrai nts. Royal Society Open Science , 9:211619, 2022
work page 2022
-
[4]
S Taheri-Araghi, S Bradde, J. T. Sauls, N. S. Hill, P. A. Le vin, J Paulsson, M Vergassola, and S Jun. Cell-size control a nd homeostasis in bacteria. Current Biology, 25(3):385–391, 2015
work page 2015
-
[5]
Effective Potential for Cellular Size Control
Stanislav Burov and David Kessler. Effective potential f or cellular size control. arXIv:1701.01725, 2017
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[6]
Division in Escherichia coli is triggered by a size-sensing rather than a timing mechanis m
Lydia Robert, Marc Hoffmann, Nathalie Krell, St´ ephane A ymerich, J´ erˆ ome Robert, and Marie Doumic. Division in Escherichia coli is triggered by a size-sensing rather than a timing mechanis m. BMC Biology , 12:17, 2014
work page 2014
-
[7]
Introduction to Structured Equations in Biology
Benoit Perthame. Introduction to Structured Equations in Biology . CNA Summer School Lecture Notes, Paris, France, 2023
work page 2023
-
[8]
J. A. J. Metz and O. Diekmann. The Dynamics of Physiologically Structured Populations . Springer Berlin, Heidelberg, 1986
work page 1986
-
[9]
Chris D. Greenman and Tom Chou. Kinetic theory of age-str uctured stochastic birth-death processes. Physical Review E , 93(1):012112, 2016
work page 2016
-
[10]
A hierarchical kinetic the ory of birth, death and fission in age-structured interactin g populations
Tom Chou and Chris D Greenman. A hierarchical kinetic the ory of birth, death and fission in age-structured interactin g populations. Journal of Statistical Physics , 164(1):49–76, 2016
work page 2016
-
[11]
Mingtao Xia and Tom Chou. Kinetic theory for structured populations: application to stochastic sizer-timer model s of cell proliferation. Journal of Physics A: Mathematical and Theoretical , 54(38):385601, 2021
work page 2021
-
[12]
Kynaston, Chris Guiver, and Christian A
Joshua C. Kynaston, Chris Guiver, and Christian A. Yate s. Equivalence framework for an age-structured multistage representation of the cell cycle. Physical Review E , 105:064411, 2022
work page 2022
-
[13]
Po-Yi Ho, Jie Lin, and Ariel Amir. Modeling cell size reg ulation: From single-cell-level statistics to molecular m echanisms and population-level effects. Annual Review of Biophysics , 47(1):251–271, 2018
work page 2018
-
[14]
Gioele La Manno, Ruslan Soldatov, Amit Zeisel, Emelie B raun, Hannah Hochgerner, Viktor Petukhov, Katja Lidschrei ber, Maria E Kastriti, Peter L¨ onnerberg, Alessandro Furlan, et al. RNA velocity of single cells. Nature, 560(7719):494–498, 2018
work page 2018
-
[15]
Mapping transcrip tomic vector fields of single cells
Xiaojie Qiu, Yan Zhang, Jorge D Martin-Rufino, Chen Weng , Shayan Hosseinzadeh, Dian Yang, Angela N Pogson, Marco Y Hein, Kyung Hoi Joseph Min, Li Wang, et al. Mapping transcrip tomic vector fields of single cells. Cell, 185(4):690–711, 2022
work page 2022
-
[16]
P rotein velocity and acceleration from single-cell multiom ics 13 experiments
Gennady Gorin, Valentine Svensson, and Lior Pachter. P rotein velocity and acceleration from single-cell multiom ics 13 experiments. Genome Biology , 21(1):1–6, 2020
work page 2020
-
[17]
A deterministic map of Waddington’s epigenetic landscape for cell fate specification
Sudin Bhattacharya, Qiang Zhang, and Melvin E Andersen . A deterministic map of Waddington’s epigenetic landscape for cell fate specification. BMC Systems Biology , 5(1):1–12, 2011
work page 2011
-
[18]
Quantifying the Waddington landscape and biological paths for devel- opment and differentiation
Jin Wang, Kun Zhang, Li Xu, and Erkang Wang. Quantifying the Waddington landscape and biological paths for devel- opment and differentiation. Proceedings of the National Academy of Sciences , 108(20):8257–8262, 2011
work page 2011
-
[19]
Regulated cell death (RCD) in cancer: key pathways and targeted therapies
Fu Peng, Minru Liao, Rui Qin, Shiou Zhu, Cheng Peng, Leil ei Fu, Yi Chen, and Bo Han. Regulated cell death (RCD) in cancer: key pathways and targeted therapies. Signal Transduction and Targeted Therapy , 7(1):286, 2022
work page 2022
-
[20]
William F Marzluff and Robert J Duronio. Histone mRNA exp ression: multiple levels of cell cycle regulation and impor tant developmental consequences. Current Opinion in Cell Biology , 14(6):692–699, 2002
work page 2002
-
[21]
Nathaniel Heintz, Hazel L Sive, and Robert G Roeder. Reg ulation of human histone gene expression: kinetics of accum u- lation and changes in the rate of synthesis and in the half-li ves of individual histone mRNAs during the HeLa cell cycle. Molecular and Cellular Biology , 3(4):539–550, 1983
work page 1983
-
[22]
Mingtao Xia, Chris D. Greenman, and Tom Chou. PDE models of adder mechanisms in cellular proliferation. SIAM Journal on Applied Mathematics , 80(3):1307–1335, 2020
work page 2020
-
[23]
mRNAs, p roteins and the emerging principles of gene expression cont rol
Christopher Buccitelli and Matthias Selbach. mRNAs, p roteins and the emerging principles of gene expression cont rol. Nature Reviews Genetics , 21(10):630–644, 2020
work page 2020
-
[24]
Chenglong Xia, Jean Fan, George Emanuel, Junjie Hao, an d Xiaowei Zhuang. Spatial transcriptome profiling by MERFIS H reveals subcellular RNA compartmentalization and cell cyc le-dependent gene expression. Proceedings of the National Academy of Sciences , 116(39):19490–19499, 2019
work page 2019
- [25]
-
[26]
Noise dis torts the epigenetic landscape and shapes cell-fate decisi ons
Megan Coomer, Lucy Ham, and Michael PH Stumpf. Noise dis torts the epigenetic landscape and shapes cell-fate decisi ons. Cell Systems , 13:83–102, 2022
work page 2022
-
[27]
Landscape and flux theory of non-equilibrium d ynamical systems with application to biology
Jin Wang. Landscape and flux theory of non-equilibrium d ynamical systems with application to biology. Advances in Physics, 64(1):1–137, 2015
work page 2015
-
[28]
P robabilistic representation of a class of non conservative nonlinear Partial Differential Equations
Anthony Le Cavil, Nadia Oudjane, and Francesco Russo. P robabilistic representation of a class of non conservative nonlinear Partial Differential Equations. Latin American Journal of Probability and Mathematical Sta tistics, 13:1189–1233, 2016
work page 2016
-
[29]
Noise di storts the epigenetic landscape and shapes cell-fate decisions
Megan A Coomer, Lucy Ham, and Michael PH Stumpf. Noise di storts the epigenetic landscape and shapes cell-fate decisions. Cell Systems , 13(1):83–102, 2022
work page 2022
-
[30]
Chen Jia, Abhyudai Singh, and Ramon Grima. Concentrati on fluctuations in growing and dividing cells: Insights into the emergence of concentration homeostasis. PLOS Computational Biology , 18(10):1–34, 10 2022
work page 2022
-
[31]
Renaud Dessalles, Maria D’Orsogna, and Tom Chou. Exact steady-state distributions of multispecies birth-death- immigration processes: effects of mutations and carrying ca pacity on diversity. Journal of Statistical Physics , 173:182–221, 2018
work page 2018
-
[32]
Mingtao Xia, Xiangting Li, and Tom Chou. Population ove rcompensation, transients, and oscillations in age-struc tured Lotka-Volterra models. arXiv preprint arXiv:2303.00864 , 2023
work page internal anchor Pith review Pith/arXiv arXiv 2023
-
[33]
Efficient scaling and moving techniques for spectral methods in unbounded domains
Mingtao Xia, Sihong Shao, and Tom Chou. Efficient scaling and moving techniques for spectral methods in unbounded domains. SIAM Journal on Scientific Computing , 43(5):A3244–A3268, 2021
work page 2021
-
[34]
A frequency-dep endent p-adaptive technique for spectral methods
Mingtao Xia, Sihong Shao, and Tom Chou. A frequency-dep endent p-adaptive technique for spectral methods. Journal of Computational Physics , 446:110627, 2021
work page 2021
-
[35]
T. Chou, S. Shao, and M. Xia. Adaptive Hermite spectral m ethods in unbounded domains. Applied Numerical Mathematics, 183:201–220, 2023
work page 2023
-
[36]
M. Xia, L. B¨ ottcher, and T. Chou. Spectrally adapted ph ysics-informed neural networks for solving unbounded doma in problems. Machine Learning: Science and Technology , 4(2):025024, 2023
work page 2023
-
[37]
Bruce P. Ayati and Isaac Klapper. A multiscale model of b iofilm as a senescence-structured fluid. Multiscale Modeling & Simulation, 6:374–365, 2007
work page 2007
-
[38]
Aggregati on of variables and applications to population dynamics
Pierre Auger, Pierre Magal, and Shigui Ruan. Aggregati on of variables and applications to population dynamics. In Pierre Magal and Shigui Ruan, editors, Structured population models in biology and epidemiology , pages 209–263. Springer Berlin, Heidelberg, 2008
work page 2008
-
[39]
Spati al structure, cooperation and competition in biofilms
Carey D Nadell, Knut Drescher, and Kevin R Foster. Spati al structure, cooperation and competition in biofilms. Nature Reviews Microbiology, 14(9):589–600, 2016. 14 Appendix A: Derivation of the differential equation satisfie d by the cell population probability density function To show pn(Xn, t |X(0)n(0), 0) defined in Eq. (4) satisfies Eq. (6), we require th...
work page 2016
-
[40]
It is straightforward to verify that F m(0; X(0)n(0), 0) = 1; therefore, we have F m(t; X(0)n(0), 0) ≡ 1, ∀t ≥ 0, which indicates that ∑ n ∫ pm n (Xn, t |X(0)n(0), 0)dXn ≤ 1. (A9) By induction, Eq. (A7) holds true for all m ∈ N+. Finally, it is easy to show that pm n (X(0)n, t |X(0)n(0), 0) ≥ 0, so by the monotone convergence theorem, ∑ n ∫ p∗ n(Xn, t |X(...
-
[41]
We shall apply Theorem 6.2 in [27]
Proof of Proposition 1 Here, we prove Proposition 1 and provide the needed technical ass umptions. We shall apply Theorem 6.2 in [27]. If n ̸= n(0), then by definition ˆ pn = 0, which solves Eq. (A1). If n(s) ≡ n(0), s ∈ [0, t ], for any smooth function φ ∈ C∞ (R|n|1), |n|1 := ∑k i=1 ni, we define the measure γ m(φ, t ) := ∫ C|n|1 φ(Xn(t; ω ))S ( t; X(t, ω ...
-
[42]
Proof of Proposition 2 We prove Proposition 1 by induction on m. Clearly, when m = 0, 1, p0 and p1 solve Eq. (A6) by using Proposition 1. If the conclusion holds for m ≥ 1, then when n ̸= n(0), we have ∂pm+1 n ∂t =E [ S ( t; X(t)n(t) ) J m(t, t; Xn, n(0)) ⏐ ⏐ ⏐Xn(0)(0), 0; n(0 < s < t ) = n(0) ] + ∫ t 0 E [ S( τ ; Xn ) ∂tJ m( t, τ; Xn, n(0)) ⏐ ⏐ ⏐X(0)n(0)...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.