pith. sign in

arxiv: 2506.16770 · v4 · submitted 2025-06-20 · 🌊 nlin.AO · q-bio.PE

Self-Balancing of Cell Populations via Martingale Turnover with Amplification

Tomoyuki Yamaguchi This is my paper

Pith reviewed 2026-05-19 08:26 UTC · model grok-4.3

classification 🌊 nlin.AO q-bio.PE
keywords martingale turnovercell population balancestochastic dynamicsLangevin equationrandom walkbiological adaptationcompetition decay model
0
0 comments X

The pith

A stochastic martingale turnover process allows cell populations to self-balance through competition and decay without type-specific regulation.

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

The paper proposes that biological systems achieve balanced cell compositions through a stochastic framework of martingale turnover with amplification. Cells proliferate solely via mutual competition and decay, with no cell-type-specific regulatory mechanisms or external signals required. Stochastic simulations and mathematical analysis demonstrate that this process generates population compositions associated with low decay probabilities. The dynamics appear as a random walk with step lengths that shrink in low-decay regions, leading to larger total populations and more compatible microscopic states. Overall, the system follows a modified Langevin equation where effective mass depends on fitness and total population size, allowing regulation of resistance to change rather than direction alone.

Core claim

Through stochastic simulations and mathematical analysis, the process of martingale turnover with amplification autonomously generates balanced population compositions associated with low decay probabilities. The compositional dynamics can be described as a random walk whose step lengths decrease in low-decay regions. Reduced decay leads to larger total population sizes and an increase in the number of compatible microscopic states, which in turn shapes the distribution of compositions under fluctuating conditions. More generally, the dynamics follow a modified Langevin equation, in which constant mass is replaced by a fitness-dependent effective mass proportional to the total population.

What carries the argument

Martingale turnover with amplification: a stochastic process of cell proliferation through mutual competition and decay that produces balanced compositions with low decay probabilities.

If this is right

  • Balanced population compositions emerge autonomously from the competition-decay rules alone.
  • Low-decay regions produce larger total population sizes and greater numbers of compatible microscopic states.
  • Compositional distributions under fluctuating conditions are shaped by the increased number of available states.
  • The system regulates resistance to change through fitness-dependent effective mass in its modified Langevin dynamics.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This mechanism could apply to other biological or ecological systems where balance arises from simple competition without complex signaling networks.
  • Disrupting only the competition rule while preserving decay might shift populations away from low-decay balanced states, providing a targeted experimental probe.
  • The random-walk description with shrinking steps suggests testable predictions for how population variance changes as total size grows.

Load-bearing premise

Proliferation occurs solely through mutual competition and decay without any cell-type-specific regulatory mechanisms or external signals.

What would settle it

A direct test would be whether stochastic simulations or biological experiments with only competition-based proliferation and decay produce stable balanced compositions and low overall decay rates, or whether balance collapses without added regulatory mechanisms.

read the original abstract

Adaptive control in biological systems, such as intestinal immunity, remains poorly understood despite detailed knowledge of underlying regulatory networks. We propose an alternative framework based on stochastic martingale turnover, in which cells proliferate through mutual competition and decay without cell-type-specific regulation. Through stochastic simulations and mathematical analysis, we show that this process autonomously generates balanced population compositions associated with low decay probabilities. The compositional dynamics can be described as a random walk whose step lengths decrease in low-decay regions. Reduced decay leads to larger total population sizes and an increase in the number of compatible microscopic states, which in turn shapes the distribution of compositions under fluctuating conditions. More generally, the dynamics follow a modified Langevin equation, in which constant mass is replaced by a fitness-dependent effective mass proportional to the total population size. Thus, biological systems regulate resistance to change, not merely direction, in shaping their macroscopic behavior.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes a framework for cell population dynamics based on stochastic martingale turnover with amplification. Cells proliferate via mutual competition and decay without cell-type-specific regulatory mechanisms. Stochastic simulations and mathematical analysis are used to argue that this process autonomously generates balanced population compositions associated with low decay probabilities. Compositional dynamics are described as a random walk with step lengths decreasing in low-decay regions, leading to larger total population sizes and more compatible states; more generally, the dynamics follow a modified Langevin equation in which constant mass is replaced by a fitness-dependent effective mass proportional to total population size.

Significance. If the central claims hold, the work offers a mechanism by which balanced compositions can emerge in biological systems (e.g., intestinal immunity) purely from stochastic competition and decay, without invoking explicit regulatory networks. The effective-mass construction in the Langevin description provides a conceptual way to think about regulation of resistance to change rather than direction, which could be of interest to models of self-organization in nonlinear dynamics and adaptation.

major comments (2)
  1. [Mathematical analysis / derivation of the modified Langevin equation] Derivation of modified Langevin equation: the mapping from the variable-step compositional random walk to the continuous modified Langevin form with fitness-dependent effective mass (proportional to total population size) is presented as the general explanation for autonomous balance, yet the validity limits of the approximation (large-population or small-fluctuation regimes) are not explicitly justified or stated. This derivation is load-bearing for the claim that the process regulates resistance to change.
  2. [Definition of effective mass in the Langevin equation] Effective-mass definition: the effective mass is defined proportionally to total population size, which is itself generated by the same turnover process. Without clarification on whether this is an emergent relation or an imposed modeling choice, it is unclear whether the low-decay stability is predicted autonomously or follows by construction.
minor comments (2)
  1. [Simulation results] Simulations are stated to support the claims, but no error bars, details on the choice of decay probabilities per cell type, or sensitivity to fitness-function choices are provided in the abstract or summary description.
  2. [Mathematical analysis] The transition from discrete random-walk description to continuous Langevin dynamics would benefit from an explicit listing of the approximation steps or scaling assumptions used.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive feedback, which helps clarify the mathematical foundations of our framework. We address the two major comments point by point below, agreeing that additional explicit statements on approximation limits and the emergent character of the effective mass will improve the manuscript.

read point-by-point responses

  1. Referee: Derivation of modified Langevin equation: the mapping from the variable-step compositional random walk to the continuous modified Langevin form with fitness-dependent effective mass (proportional to total population size) is presented as the general explanation for autonomous balance, yet the validity limits of the approximation (large-population or small-fluctuation regimes) are not explicitly justified or stated. This derivation is load-bearing for the claim that the process regulates resistance to change.

    Authors: We agree that the validity regime of the continuous approximation should be stated explicitly. The mapping from the discrete variable-step random walk to the modified Langevin equation is derived in the large-population limit where relative fluctuations are small (N ≫ 1) and individual step sizes remain small compared with the scale of compositional variation. Under these conditions the Fokker-Planck description becomes accurate and the effective-mass term appears naturally from the 1/√N scaling of step lengths. We will add a short paragraph in the Methods section (or a new subsection) that states these assumptions and notes that all reported simulations operate inside this regime. This revision will make the load-bearing role of the derivation fully transparent. revision: yes

  2. Referee: Effective-mass definition: the effective mass is defined proportionally to total population size, which is itself generated by the same turnover process. Without clarification on whether this is an emergent relation or an imposed modeling choice, it is unclear whether the low-decay stability is predicted autonomously or follows by construction.

    Authors: The effective mass is strictly emergent. Total population size N is not prescribed; it is the steady-state outcome of the same stochastic proliferation-decay rules that govern composition. In the continuum limit the diffusion coefficient of the compositional walk scales as 1/N, which is mathematically equivalent to an inertia (effective mass) proportional to N. Because N itself grows when decay probabilities are low, the resistance to compositional change increases automatically. We will revise the text to include a concise derivation sketch showing how N enters the Langevin coefficients directly from the microscopic transition rates, thereby demonstrating that the stability is a prediction of the autonomous process rather than an imposed feature. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The abstract describes compositional dynamics as a random walk with step lengths decreasing in low-decay regions, leading to larger total population sizes and more compatible states that shape the distribution. The modified Langevin equation is presented as a general description in which effective mass is proportional to total population size. This modeling choice links fitness (low decay) to population size and thus to reduced fluctuations, but does not reduce the central claim of autonomous balance to a definition or fitted input by construction. No equations are shown that make the stability prediction equivalent to the input assumptions, no self-citations are load-bearing, and no uniqueness theorems or ansatzes are invoked. The derivation from discrete random walk to continuous form is stated without evidence of tautological reduction in the provided text, making the analysis self-contained.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the assumption that decay probabilities can be treated as fixed parameters that are not themselves regulated by cell-type identity, plus the modeling choice that total population size directly scales the effective mass in the Langevin description.

free parameters (1)
  • decay probability per cell type
    Treated as an input that is varied to demonstrate balance at low values; no independent measurement or derivation is supplied in the abstract.
axioms (1)
  • domain assumption Proliferation occurs through mutual competition and decay without cell-type-specific regulation
    Invoked at the outset of the proposed framework in the abstract.

pith-pipeline@v0.9.0 · 5675 in / 1351 out tokens · 28655 ms · 2026-05-19T08:26:04.890817+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

26 extracted references · 26 canonical work pages · 1 internal anchor

  1. [1]

    Meyer and J

    P. Meyer and J. Saez -Rodriguez, Advances in systems biology modeling: 10 years of crowdsourcing DREAM challenges, Cell Syst 12, 636 (2021)

    work page 2021

  2. [2]

    Z. R. Thornburg et al., Fundamental behaviors emerge from simulations of a living minimal cell, Cell 185, 345 (2022)

    work page 2022

  3. [3]

    Eling, M

    N. Eling, M. D. Morgan, and J. C. Marioni, Challenges in measuring and understanding biological noise, Nature Reviews Genetics 20, 536 (2019)

    work page 2019

  4. [4]

    Kedmi and D

    R. Kedmi and D. R. Littman, Antigen-presenting cells as specialized drivers of intestinal T cell functions, Immunity 57, 2269 (2024)

    work page 2024

  5. [5]

    Akagbosu et al

    B. Akagbosu et al. , Novel antigen -presenting cell imparts T reg-dependent tolerance to gut microbiota, Nature 610, 752 (2022)

    work page 2022

  6. [6]

    Lyu et al., ILC3s select microbiota -specific regulatory T cells to establish tolerance in the gut, Nature 610, 744 (2022)

    M. Lyu et al., ILC3s select microbiota -specific regulatory T cells to establish tolerance in the gut, Nature 610, 744 (2022)

    work page 2022

  7. [7]

    Kedmi et al

    R. Kedmi et al. , A ROR γt+ cell instructs gut microbiota-specific T reg cell differentiation, Nature 610, 737 (2022)

    work page 2022

  8. [8]

    P. F. Rodrigues et al. , Rorγt-positive dendritic cells are required for the induction of peripheral regulatory T cells in response to oral antigens, Cell 188, 2720 (2025)

    work page 2025

  9. [9]

    L. Fu, R. Upadhyay, M. Pokrovskii, F. M. Chen, G. Romero -Meza, A. Griesemer, and D. R. Littman, Prdm16-dependent antigen -presenting cells induce tolerance to gut antigens, Nature (2025)

    work page 2025

  10. [10]

    M. C. Campos Canesso et al., Identification of antigen-presenting cell –T cell interactions driving immune responses to food, Science 387, eado5088 (2025)

    work page 2025

  11. [11]

    Silver et al., Mastering the game of Go with deep neural networks and tree search, Nature 529, 484 (2016)

    D. Silver et al., Mastering the game of Go with deep neural networks and tree search, Nature 529, 484 (2016)

    work page 2016

  12. [12]

    DeepSeek-AI et al., DeepSeek-R1: Incentivizing Reasoning Capability in LLMs via Reinforcement Learning, arXiv, 2501.12948 (2025)

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  13. [13]

    Pezzulo and M

    G. Pezzulo and M. Levin, Top-down models in biology: explanation and control of complex living systems above the molecular level, J R Soc Interface 13 (2016)

    work page 2016

  14. [14]

    E. O. Neftci and B. B. Averbeck, Reinforcement learning in artificial and biological systems, Nature Machine Intelligence 1, 133 (2019)

    work page 2019

  15. [15]

    Vanchurin, Y

    V . Vanchurin, Y . I. Wolf, M. I. Katsnelson, and E. V . Koonin, Toward a theory of evolution as multilevel learning, Proc Natl Acad Sci U S A 119 (2022)

    work page 2022

  16. [16]

    Yamaguchi, Learning processes in hierarchical pairs regulate entire gene expression in cells, Sci Rep 12, 7549 (2022)

    T. Yamaguchi, Learning processes in hierarchical pairs regulate entire gene expression in cells, Sci Rep 12, 7549 (2022)

    work page 2022

  17. [17]

    K. B. Athreya and S. Karlin, Embedding of Urn Schemes into Continuous Time Markov Branching Processes and Related Limit Theorems, The Annals of Mathematical Statistics 39, 1801 (1968)

    work page 1968

  18. [18]

    J. L. Doob, Stochastic Processes (Wiley, 1953). 6 *Contact author: t.yamaguchi@tokushukai.jp

    work page 1953

  19. [19]

    See Supplemental Material [url] for derivation of SDE, 1-WGM model, and random walk

    work page

  20. [20]

    D. T. Gillespie, The chemical Langevin equation, The Journal of Chemical Physics 113, 297 (2000)

    work page 2000

  21. [21]

    Yamaguchi, S

    T. Yamaguchi, S. Teraguchi, C. Furusawa, H. Machiyama, T. M. Watanabe, H. Fujita, S. Sakaguchi, and T. Yanagida, Theoretical modeling reveals that regulatory T cells increase T -cell interaction with antigen-presenting cells for stable immune tolerance, Int Immunol 31, 743 (2019)

    work page 2019

  22. [22]

    P. W. Reddien, The purpose and ubiquity of turnover, Cell 187, 2657 (2024)

    work page 2024

  23. [23]

    Black and M

    F. Black and M. Scholes, The Pricing of Options and Corporate Liabilities, Journal of Political Economy 81, 637 (1973)

    work page 1973

  24. [24]

    Katori, in Bessel Processes, Schramm – Loewner Evolution, and the Dyson Model (Springer Singapore, Singapore, 2016), pp

    M. Katori, in Bessel Processes, Schramm – Loewner Evolution, and the Dyson Model (Springer Singapore, Singapore, 2016), pp. 1

    work page 2016

  25. [25]

    Eberl, S

    G. Eberl, S. Marmon, M. J. Sunshine, P. D. Rennert, Y . Choi, and D. R. Littman, An essential function for the nuclear receptor RORgamma(t) in the generation of fetal lymphoid tissue inducer cells, Nat Immunol 5, 64 (2004)

    work page 2004

  26. [26]

    demand per cell

    A. Jarade, J. P. Di Santo, and N. Serafini, Group 3 innate lymphoid cells mediate host defense against attaching and effacing pathogens, Curr Opin Microbiol 63, 83 (2021). END MATTER Appendix A: Simulation Methods — In the model describing the dynamics of cell numbers NⅠ(t) and NⅡ(t) under MTA and CSC conditions, a series of three stochastic process es — ...

    work page 2021