Phase transitions in microbial lineage trees
Pith reviewed 2026-05-10 07:12 UTC · model grok-4.3
The pith
A first-order phase transition in a bacterial plasmid model imposes a strict lower bound on the number of plasmids that populations can stably maintain.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Statistical physics predicts that microbial populations can undergo phase transitions in which behavior changes discontinuously with small perturbations. These transitions are visible in the genealogies recorded by lineage trees. In a model that incorporates plasmid replication, segregation, and selection, we rigorously establish the existence of a first-order phase transition and obtain a strict lower bound on the plasmid number that remains stable in the population.
What carries the argument
The model of plasmid replication, segregation, and selection operating inside microbial lineage trees, which produces the discontinuous switch in stable plasmid number.
If this is right
- Plasmid numbers below the calculated bound cannot be maintained at steady state regardless of initial conditions.
- Small changes in replication or segregation rates can trigger an abrupt loss of plasmids across the population.
- The genealogical structure of the population directly determines where the phase boundary lies.
- Synthetic constructs that attempt to introduce too few plasmids will fail to persist under the model's selection regime.
Where Pith is reading between the lines
- The bound may constrain the design space for multi-plasmid synthetic circuits in biotechnology.
- Similar transitions could appear in other mobile genetic elements such as transposons or prophages.
- Lineage-tree statistics might serve as an observable proxy for detecting phase transitions without direct fitness measurements.
Load-bearing premise
The model's rules for how plasmids replicate, are divided between daughter cells, and confer fitness accurately represent the dominant processes operating in real bacterial populations.
What would settle it
Measure plasmid copy numbers across many generations in a controlled bacterial population and check whether the fraction of cells retaining plasmids jumps discontinuously at the parameter value predicted by the model, or whether populations can maintain any plasmid number below the derived bound.
Figures
read the original abstract
Statistical physics can describe the behavior of microbial populations consisting of many heterogeneous individuals. A direct consequence is the existence of phase transitions, where the behavior of a population changes discontinuously upon a small perturbation. While such phase transitions have often been proposed in biology, connecting observed behavior to the underlying physics has remained challenging. We show how phase transitions naturally arise in microbial population dynamics and highlight their connection with genealogies. We rigorously demonstrate the existence of a first-order phase transition in a model of bacterial plasmid engineering and find a strict lower bound on the number of plasmids that can be stably maintained in a population.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses statistical physics to analyze microbial population dynamics and lineage trees, claiming that phase transitions arise naturally in such systems. It rigorously demonstrates the existence of a first-order phase transition within a defined model of bacterial plasmid replication, segregation, and selection, and derives a strict lower bound on the number of plasmids that can be stably maintained in a population.
Significance. If the central derivation holds, the work provides a self-contained mathematical link between statistical mechanics and microbial genealogy, offering a parameter-free explanation for discontinuous population-level changes with relevance to plasmid engineering and evolutionary dynamics. The absence of free parameters and the internal consistency of the model (as confirmed by the axiom ledger) are notable strengths that make the bound falsifiable within the model's scope.
minor comments (3)
- [Abstract] The abstract asserts a 'rigorous demonstration' and 'strict lower bound' but does not preview the key model equations or order parameter; adding a brief reference to the master equation or free-energy functional would improve clarity for readers.
- [Model definition] Notation for the plasmid copy number and segregation bias should be defined consistently at first use to avoid ambiguity when transitioning from the general lineage-tree framework to the specific plasmid model.
- [Discussion] The discussion of experimental testability is brief; expanding it with a concrete prediction (e.g., a threshold plasmid number observable in flow-cytometry time series) would strengthen the connection to biology without altering the mathematical claims.
Simulated Author's Rebuttal
We thank the referee for their positive assessment of our manuscript and for recommending minor revision. The referee's summary accurately captures the central result: a first-order phase transition in the plasmid dynamics model that yields a strict, parameter-free lower bound on stably maintained plasmids. No major comments were raised in the report.
Circularity Check
No significant circularity; derivation self-contained in model
full rationale
The paper defines a mathematical model of plasmid replication, segregation, and selection, then derives the existence of a first-order phase transition and a strict lower bound directly from the model's internal dynamics and equations. No load-bearing step reduces to a fitted parameter, self-definition, or self-citation chain; the claims are presented as consequences of the specified processes. The abstract and skeptic analysis confirm the demonstration is internal to the model without external reduction or renaming of known results.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Plasmid replication and segregation follow a statistical process whose moments can be coarse-grained into an effective population-level dynamics.
Reference graph
Works this paper leans on
-
[1]
P. Piho and P. Thomas, Feedback between stochastic gene networks and population dynamics enables cellular decision- making, Sci. Adv.10(2024)
work page 2024
- [2]
- [3]
- [4]
- [5]
-
[6]
P. W. Anderson, More is different, Science177(1972)
work page 1972
-
[7]
R. V. Sol´ e, S. C. Manrubia, B. Luque, J. Delgado, and J. Bascompte, Phase transitions and complex systems: simple, nonlinear models capture complex systems at the edge of chaos, Complexity1(1996)
work page 1996
-
[8]
M. A. Mu˜ noz, Colloquium: Criticality and dynamical scaling in living systems, Rev. Mod. Phys.90(2018)
work page 2018
-
[9]
S. W. Ordway, D. M. King, D. Friend, C. Noto, S. Phu, H. Huelskamp, R. F. Inglis, W. Olivas, and S. Bahar, Phase transition behaviour in yeast and bacterial populations under stress, R. Soc. Open Sci.7(2020)
work page 2020
-
[10]
E. F. W. Heffern, H. Huelskamp, S. Bahar, and R. F. Inglis, Phase transitions in biology: from bird flocks to population dynamics, Proc. R. Soc. B288(2021)
work page 2021
-
[11]
S. Bhattacharyya, Emergence of first-order and second-order phase transitions in a cyclic ecosystem exposed to environ- mental impact, Phys. Scr.99(2024)
work page 2024
-
[12]
D. Bagchi and P. K. Mohanty, Phase transition in an exactly solvable extinction model, Phys. Rev. E84(2011)
work page 2011
-
[13]
A. Skanata and E. Kussell, Evolutionary phase transitions in random environments, Phys. Rev. Lett.117(2016)
work page 2016
- [14]
-
[15]
J. Lin, J. Min, and A. Amir, Optimal segregation of proteins: phase transitions and symmetry breaking, Phys. Rev. Lett. 122(2019)
work page 2019
-
[16]
Y. Wakamoto, A. Y. Grosberg, and E. Kussell, Optimal lineage principle for age-structured populations, Evolution66 (2012)
work page 2012
-
[17]
R. Garc´ ıa-Garc´ ıa, A. Genthon, and D. Lacoste, Linking lineage and population observables in biological branching processes, Phys. Rev. E99(2019)
work page 2019
-
[18]
T. J. Kobayashi and Y. Sughiyama, Fluctuation relations of fitness and information in population dynamics, Phys. Rev. Lett.115(2015)
work page 2015
-
[19]
Y. Sughiyama, T. J. Kobayashi, K. Tsumura, and K. Aihara, Pathwise thermodynamic structure in population dynamics, Phys. Rev. E91(2015)
work page 2015
-
[20]
Pigolotti, Generalized Euler-Lotka equation for correlated cell divisions, Phys
S. Pigolotti, Generalized Euler-Lotka equation for correlated cell divisions, Phys. Rev. E103(2021)
work page 2021
-
[21]
The two-clock problem in population dynamics
K. ¨Ocal and M. P. H. Stumpf, The two-clock problem in population dynamics (2025), arXiv:2504.20388
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[22]
P. Rugbjerg, K. Sarup-Lytzen, M. Nagy, and M. O. A. Sommer, Synthetic addiction extends the productive life time of engineeredEscherichia colipopulations, PNAS115(2018)
work page 2018
-
[23]
J. C. R. Hernandez-Beltran, V. Mir´ o Pina, A. Siri-J´ egousse, S. Palau, R. Pe˜ na-Miller, and A. Gonz´ alez Casanova, Segre- gational instability of multicopy plasmids: a population genetics approach, Ecol. Evol.12(2022)
work page 2022
-
[24]
T. Stadler, O. G. Pybus, and M. P. H. Stumpf, Phylodynamics for cell biologists, Science371(2021)
work page 2021
-
[25]
A. Genthon, T. Nozoe, L. Peliti, and D. Lacoste, Cell lineage statistics with incomplete population trees, PRX Life1 (2023)
work page 2023
-
[26]
K. B. Athreya and P. E. Ney,Branching Processes(Springer, Berlin, Heidelberg, 1972)
work page 1972
-
[27]
Thomas, Making sense of snapshot data: ergodic principle for clonal cell populations, J
P. Thomas, Making sense of snapshot data: ergodic principle for clonal cell populations, J. R. Soc. Interface14(2017)
work page 2017
-
[28]
J. H. H. Gillespie,Population Genetics: A Concise Guide(Johns Hopkins University Press, Baltimore, 2004)
work page 2004
-
[29]
Chevin, On measuring selection in experimental evolution, Biol
L.-M. Chevin, On measuring selection in experimental evolution, Biol. Lett.7(2011)
work page 2011
-
[30]
F. Jafarpour, E. Levien, and A. Amir, Evolutionary dynamics in non-Markovian models of microbial populations, Phys. Rev. E108(2023)
work page 2023
-
[31]
Kimura, On the probability of fixation of mutant genes in a population, Genetics47(1962)
M. Kimura, On the probability of fixation of mutant genes in a population, Genetics47(1962). 11
work page 1962
-
[32]
B. R. Glick, Metabolic load and heterologous gene expression, Biotechn. Adv.13(1995)
work page 1995
-
[33]
E. Ilker and M. Hinczewski, Modeling the growth of organisms validates a general relation between metabolic costs and natural selection, Phys. Rev. Lett.122(2019)
work page 2019
-
[34]
M. V. Rouches, Y. Xu, L. B. G. Cortes, and G. Lambert, A plasmid system with tunable copy number, Nat. Commun.13 (2022)
work page 2022
-
[35]
J. Rodr´ ıguez-Beltr´ an, J. DelaFuente, R. Le´ on-Sampedro, R. C. MacLean, and´A. San Mill´ an, Beyond horizontal gene transfer: the role of plasmids in bacterial evolution, Nat. Rev. Microbiol.19(2021)
work page 2021
-
[36]
Touchette, The large deviation approach to statistical mechanics, Phys
H. Touchette, The large deviation approach to statistical mechanics, Phys. Rep.478(2009)
work page 2009
-
[37]
A. Genthon and D. Lacoste, Fluctuation relations and fitness landscapes of growing cell populations, Sci. Rep.10(2020)
work page 2020
-
[38]
A. D. Scott, D. M. King, S. W. Ordway, and S. Bahar, Phase transitions in evolutionary dynamics, Chaos32(2022)
work page 2022
-
[39]
R. B. Bapat and T. E. S. Raghavan,Nonnegative Matrices and Applications(Cambridge University Press, 1997)
work page 1997
-
[40]
Kato,Perturbation Theory for Linear Operators(Springer, Berlin, Heidelberg, 1995)
T. Kato,Perturbation Theory for Linear Operators(Springer, Berlin, Heidelberg, 1995)
work page 1995
-
[41]
C. Giardin` a, J. Kurchan, and L. Peliti, Direct evaluation of large-deviation functions, Phys. Rev. Lett.96(2006)
work page 2006
-
[42]
V. Lecomte and J. Tailleur, A numerical approach to large deviations in continuous time, J. Stat. Mech.2007(2007)
work page 2007
- [43]
-
[44]
A. Dembo and O. Zeitouni,Large Deviations Techniques and Applications(Springer, Berlin, Heidelberg, 2010)
work page 2010
-
[45]
P. W. Glynn and W. Whitt, Large deviations behavior of counting processes and their inverses, Queueing Syst.17(1994)
work page 1994
-
[46]
K. Duffy and A. P. Metcalfe, How to estimate the rate function of a cumulative process, J. Appl. Probab.42(2005)
work page 2005
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.