pith. sign in

arxiv: 2606.09040 · v1 · pith:4L7D2BBYnew · submitted 2026-06-08 · 🧬 q-bio.PE · cond-mat.stat-mech· physics.bio-ph

Natural Selection in the Wake of Catastrophe

Pith reviewed 2026-06-27 14:20 UTC · model grok-4.3

classification 🧬 q-bio.PE cond-mat.stat-mechphysics.bio-ph
keywords natural selectioncatastrophe recoveryfitness landscapesLevenberg-Marquardt optimizationE. coli adaptationmean fitnessevolutionary trajectories
0
0 comments X

The pith

Natural selection after catastrophe drives mean fitness to relax inversely with time, scaled by the number of coupled traits.

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

The paper develops a theory of how populations adapt via natural selection right after a catastrophe removes prior niche adaptations. It derives a simple law in which mean fitness relaxes inversely with time, and the prefactor scales with the number of traits that couple to the new environment. The theory is tested on experimental fitness landscapes from E. coli recovering after antibiotic exposure. Trait adaptation on these landscapes follows Levenberg-Marquardt optimization instead of gradient ascent, producing trajectories that avoid greedy steps near fitness peaks.

Core claim

In the wake of catastrophe, natural selection causes the mean fitness of the population to relax inversely with time, with a prefactor proportional to the number of traits coupled to the post-catastrophe environment. The mean trait adaptation follows Levenberg-Marquardt optimization rather than gradient ascent on the fitness landscape. Near fitness peaks, evolutionary trajectories are biased against greediness, making post-catastrophic selection optimistic from an optimization perspective.

What carries the argument

The inverse-time relaxation law for mean fitness together with the Levenberg-Marquardt optimization rule that governs how mean traits change.

Load-bearing premise

That fitness landscapes measured in E. coli after antibiotic exposure represent general post-catastrophe recovery and that the trajectories match Levenberg-Marquardt optimization rather than some other process.

What would settle it

A post-catastrophe experiment in which mean fitness is tracked over time and fails to show inverse scaling with time, or in which fitted trajectories on comparable landscapes match gradient ascent instead of Levenberg-Marquardt steps.

Figures

Figures reproduced from arXiv: 2606.09040 by David Pincus, Jesse Young Lin, Joshua Sodicoff, Omer Granek, Seppe Kuehn, Vincenzo Vitelli.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
read the original abstract

Living organisms, from bacteria to humans, are more likely to survive if their traits enhance fitness. In populations well adapted to their environmental niches, natural selection proceeds via rarely beneficial mutations. But when a catastrophe wipes out niche diversity, sudden adaptation often follows. Here, we present a data-validated theory of natural selection in the wake of catastrophe and unveil a simple law that emerges during recovery: the mean fitness relaxes inversely with time, with a prefactor proportional to the number of traits coupled to the post-catastrophe environment. We put our approach to test using experimental fitness landscapes measured following antibiotic administration to E. coli. The resulting mean trait adaptation is not described by gradient ascent on a fitness landscape, instead it follows an algorithm known as Levenberg-Marquardt optimization. Near fitness peaks, evolutionary trajectories are biased against greediness - from an optimization perspective, post-catastrophic selection is optimistic.

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

3 major / 0 minor

Summary. The paper presents a theory of natural selection following catastrophe, deriving that mean fitness relaxes inversely with time (prefactor proportional to the number of traits coupled to the post-catastrophe environment). It tests this on experimental fitness landscapes from antibiotic-treated E. coli and claims that mean trait adaptation follows Levenberg-Marquardt optimization rather than gradient ascent, with trajectories biased against greediness near peaks.

Significance. If the inverse-time relaxation law and the specific mapping to Levenberg-Marquardt optimization are substantiated, the work would supply a simple, potentially parameter-free prediction for post-catastrophe recovery dynamics and a concrete link between evolutionary trajectories and a damped least-squares algorithm, which could be tested in other systems and inform optimization theory.

major comments (3)
  1. [Abstract] Abstract: the claim of a 'data-validated theory' is not accompanied by any derivation steps for the inverse-time law, data details, exclusion criteria, or error analysis, preventing evaluation of the central claims from the given information.
  2. [Abstract] Abstract: the assertion that trajectories follow Levenberg-Marquardt optimization (rather than gradient ascent or other damped methods) requires a quantitative metric or error comparison on the E. coli data demonstrating superiority; qualitative similarity alone does not secure the distinction.
  3. [Abstract] Abstract: the proportionality of the prefactor to the number of coupled traits is presented as emerging from the model, but without the explicit assumptions or derivation shown, it is unclear whether this is a genuine prediction or inherits post-hoc choices.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments on our manuscript. We respond point-by-point to the major comments below. Revisions have been made to the abstract and supporting analyses to improve clarity and provide the requested details.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the claim of a 'data-validated theory' is not accompanied by any derivation steps for the inverse-time law, data details, exclusion criteria, or error analysis, preventing evaluation of the central claims from the given information.

    Authors: The abstract is a high-level summary; the full derivation of the inverse-time relaxation (from the multi-trait dynamical equations) appears in Section 2 of the main text. The E. coli fitness landscape data source, trajectory exclusion criteria, and error analysis are described in the Methods and Supplementary Information. To address the concern, we have revised the abstract to briefly note the key modeling assumptions and data origin. revision: yes

  2. Referee: [Abstract] Abstract: the assertion that trajectories follow Levenberg-Marquardt optimization (rather than gradient ascent or other damped methods) requires a quantitative metric or error comparison on the E. coli data demonstrating superiority; qualitative similarity alone does not secure the distinction.

    Authors: The original manuscript presents trajectory comparisons but relies primarily on visual and qualitative agreement. We agree a quantitative metric strengthens the distinction. The revised manuscript adds a root-mean-square deviation metric between observed mean trait paths and those predicted by Levenberg-Marquardt versus gradient ascent (and other damped methods), with statistical tests confirming superior fit for Levenberg-Marquardt on the experimental data. revision: yes

  3. Referee: [Abstract] Abstract: the proportionality of the prefactor to the number of coupled traits is presented as emerging from the model, but without the explicit assumptions or derivation shown, it is unclear whether this is a genuine prediction or inherits post-hoc choices.

    Authors: The proportionality follows directly from integrating the multi-trait response equations under the assumption of linear trait-environment coupling (see Eq. 3 and derivation in Section 2). It is a model prediction, not a post-hoc fit. We have revised the abstract to state this assumption explicitly and reference the derivation. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation self-contained against external data

full rationale

The paper derives the inverse-time mean-fitness relaxation from explicit model assumptions on trait coupling to the post-catastrophe environment and then compares resulting trajectories to experimental E. coli fitness landscapes. No quoted equation or step reduces a claimed prediction to a fitted parameter or self-citation by construction; the LM-optimization identification is presented as an empirical match rather than a definitional renaming. The central claims remain independently falsifiable on the external dataset.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract only; no explicit free parameters, axioms, or invented entities can be extracted from the provided text.

pith-pipeline@v0.9.1-grok · 5702 in / 1055 out tokens · 22573 ms · 2026-06-27T14:20:14.161112+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

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

  1. [1]

    Cox and D

    D. Cox and D. Oakes,Analysis of Survival Data(Chap- man and Hall/CRC, 1984)

  2. [2]

    P. K. Andersen, Ø. Borgan, R. D. Gill, and N. Keiding, Statistical Models Based on Counting Processes, Springer Series in Statistics (Springer US, New York, NY, 1993)

  3. [3]

    Narteau, P

    C. Narteau, P. Shebalin, and M. Holschneider, Journal of Geophysical Research: Solid Earth107, 2359 (2002)

  4. [4]

    J. W. Vaupel, K. G. Manton, and E. Stallard, Demogra- phy16, 439 (1979)

  5. [5]

    O. O. Aalen, Statistics in Medicine7, 1121 (1988)

  6. [6]

    Hougaard, Lifetime Data Analysis1, 255 (1995)

    P. Hougaard, Lifetime Data Analysis1, 255 (1995)

  7. [7]

    T. I. Missov and M. Finkelstein, Theoretical Population Biology80, 64 (2011)

  8. [8]

    T.A.BalanandH.Putter,StatisticalMethodsinMedical Research29, 3424 (2020)

  9. [9]

    Holmlund, The Scandinavian Journal of Economics 100, 113 (1998)

    B. Holmlund, The Scandinavian Journal of Economics 100, 113 (1998)

  10. [10]

    N. M. Kiefer, Journal of Economic Literature26, 646 (1988)

  11. [11]

    Rogerson, R

    R. Rogerson, R. Shimer, and R. Wright, Journal of Eco- nomic Literature43, 959 (2005)

  12. [12]

    Lancaster,The Econometric Analysis of Transition Data(Cambridge University Press, Cambridge, 1990)

    T. Lancaster,The Econometric Analysis of Transition Data(Cambridge University Press, Cambridge, 1990)

  13. [13]

    J. W. Vaupel and A. I. Yashin, The American Statistician 39, 176 (1985)

  14. [14]

    T. A. Balan and H. Putter, Statistics in Medicine38, 3405 (2019)

  15. [15]

    B. H. Good and O. Hallatschek, Current Opinion in Microbiology45, 203 (2018)

  16. [16]

    M. T. Wortel, D. Agashe, S. F. Bailey, C. Bank, K. Biss- chop, T. Blankers, J. Cairns, E. S. Colizzi, D. Cusseddu, 10 M. M. Desai, B. van Dijk, M. Egas, J. Ellers, A. T. Groot, D. G. Heckel, M. L. Johnson, K. Kraaijeveld, J. Krug, L. Laan, M. Lässig, P. A. Lind, J. Meijer, L. M. Noble, S. Okasha, P. B. Rainey, D. E. Rozen, S. Shitut, S. J. Tans, O. Tenail...

  17. [17]

    J.A.AscensaoandM.M.Desai,NatureReviewsGenetics 27, 81 (2026)

  18. [18]

    D. J. Hernandez, A. S. David, E. S. Menges, C. A. Searcy, and M. E. Afkhami, The ISME Journal15, 1722 (2021)

  19. [19]

    O. F. A. Larsen and L. H. M. van de Burgwal, Frontiers in Microbiology12, 784797 (2021)

  20. [20]

    O. F. A. Larsen and R. J. M. Brummer, Frontiers in Microbiology15, 1449133 (2024)

  21. [21]

    Y. Hong, H. Li, L. Chen, H. Su, B. Zhang, Y. Luo, C. Li, Z. Zhao, Y. Shao, and L. Guo, Microbiome12, 80 (2024)

  22. [22]

    Patnaik and G

    S. Patnaik and G. Baruah, Theoretical Ecology18, 10 (2025)

  23. [23]

    C. C. Xu, V. Fugère, N. Barbosa da Costa, B. E. Beisner, G. Bell, M. E. Cristescu, G. F. Fussmann, A. Gonzalez, B. J. Shapiro, and R. D. Barrett, Current Biology35, 1061 (2025)

  24. [24]

    Hildebrand, L

    F. Hildebrand, L. Moitinho-Silva, S. Blasche, M. T. Jahn, T. I. Gossmann, J. Huerta-Cepas, R. Hercog, M. Luetge, M. Bahram, A. Pryszlak, R. J. Alves, S. M. Waszak, A. Zhu, L. Ye, P. I. Costea, S. Aalvink, C. Belzer, S. K. Forslund, S. Sunagawa, U. Hentschel, C. Merten, K. R. Patil, V. Benes, and P. Bork, Gut68, 1781 (2019)

  25. [25]

    M. P. Thakur, W. H. van der Putten, F. Apon, E. An- gelini, B. Vreš, and S. Geisen, Functional Ecology35, 216 (2021)

  26. [26]

    L. Dai, D. Vorselen, K. S. Korolev, and J. Gore, Science 336, 1175 (2012)

  27. [27]

    Dakos, B

    V. Dakos, B. Matthews, A. P. Hendry, J. Levine, N. Loeuille, J. Norberg, P. Nosil, M. Scheffer, and L. De Meester, Nature Ecology & Evolution3, 355 (2019)

  28. [28]

    Alcantara-Diaz, Mutagenesis19, 349 (2004)

    D. Alcantara-Diaz, Mutagenesis19, 349 (2004)

  29. [29]

    Javed, M

    S. Javed, M. Asgher, M. A. Sheikh, and H. Nawaz, Food Biotechnology24, 165 (2010)

  30. [30]

    Shibai, Y

    A. Shibai, Y. Takahashi, Y. Ishizawa, D. Motooka, S. Nakamura, B.-W. Ying, and S. Tsuru, Scientific Re- ports7, 14531 (2017)

  31. [31]

    F. K. Richtmyer and E. H. Kennard,Introduction to Modern Physics, 4th ed., International Series in Pure and Applied Physics (McGraw-Hill Book Company, New York, 1947)

  32. [32]

    M. A. Nowak,Evolutionary Dynamics: Exploring the Equations of Life(Harvard University Press, 2006)

  33. [33]

    Without loss of generality, we considern non-neutral phenotypic axes in the post-catastrophe environment

  34. [34]

    W. H. Fleming and M. Viot, Indiana University Mathe- matics Journal28, 817 (1979)

  35. [35]

    A. S. Bhat, Theoretical Population Biology161, 1 (2025)

  36. [36]

    Kimura, Proceedings of the National Academy of Sciences54, 731 (1965)

    M. Kimura, Proceedings of the National Academy of Sciences54, 731 (1965)

  37. [37]

    L. S. Tsimring, H. Levine, and D. A. Kessler, Physical Review Letters76, 4440 (1996)

  38. [38]

    P. J. Gerrish, A. Colato, A. S. Perelson, and P. D. Sniegowski, Proceedings of the National Academy of Sciences104, 6266 (2007)

  39. [39]

    P. D. Sniegowski and P. J. Gerrish, Philosophical Trans- actions of the Royal Society B: Biological Sciences365, 1255 (2010)

  40. [40]

    M. M. Desai and D. S. Fisher, Genetics188, 997 (2011)

  41. [41]

    Alfaro and R

    M. Alfaro and R. Carles, Proceedings of the American Mathematical Society145, 5315 (2017)

  42. [42]

    Hamel, F

    F. Hamel, F. Lavigne, G. Martin, and L. Roques, Non- linear Analysis: Real World Applications54, 103107 (2020)

  43. [43]

    Mahdisoltani, P

    S. Mahdisoltani, P. Murugan, A. K. Chakraborty, and M. Kardar, Physical Review E110, 064137 (2024)

  44. [44]

    Tenaillon, Annual Review of Ecology, Evolution, and Systematics45, 179 (2014)

    O. Tenaillon, Annual Review of Ecology, Evolution, and Systematics45, 179 (2014)

  45. [45]

    Huang,Statistical Mechanics, 2nd ed

    K. Huang,Statistical Mechanics, 2nd ed. (J. Wiley & Sons, New York, 1987)

  46. [46]

    G. I. Barenblatt,Scaling, self-similarity, and intermediate asymptotics, Cambridge Texts in Applied Mathematics No. 14 (Cambridge University Press, Cambridge, 1996)

  47. [47]

    Lande and S

    R. Lande and S. J. Arnold, Evolution37, 1210 (1983)

  48. [48]

    R. M. Otto, A. Turska-Nowak, P. M. Brown, and K. A. Reynolds, Cell Systems15, 134 (2024)

  49. [49]

    Papkou, L

    A. Papkou, L. Garcia-Pastor, J. A. Escudero, and A. Wag- ner, Science382, eadh3860 (2023)

  50. [50]

    Clauset, C

    A. Clauset, C. R. Shalizi, and M. E. J. Newman, SIAM Review51, 661 (2009)

  51. [51]

    J.-P. Côté, S. French, S. S. Gehrke, C. R. MacNair, C. S. Mangat, A. Bharat, and E. D. Brown, mBio7, e01714-16 (2016)

  52. [52]

    Kacser and J

    H. Kacser and J. A. Burns, Genetics97, 639 (1981)

  53. [53]

    A. M. Dean, D. E. Dykhuizen, and D. L. Hartl, Genetical Research48, 1 (1986)

  54. [54]

    D. E. Dykhuizen and A. M. Dean, Trends in Ecology & Evolution5, 257 (1990)

  55. [55]

    M. S. Newton, V. L. Arcus, M. L. Gerth, and W. M. Patrick, Current Opinion in Structural Biology48, 110 (2018)

  56. [56]

    Bar-Even, E

    A. Bar-Even, E. Noor, Y. Savir, W. Liebermeister, D. Da- vidi, D. S. Tawfik, and R. Milo, Biochemistry50, 4402 (2011)

  57. [57]

    Y. T. Tamer, I. K. Gaszek, H. Abdizadeh, T. A. Batur, K. A. Reynolds, A. R. Atilgan, C. Atilgan, and E. Toprak, Molecular Biology and Evolution36, 1533 (2019)

  58. [58]

    Krucinska, M

    J. Krucinska, M. N. Lombardo, H. Erlandsen, A. Estrada, D. Si, K. Viswanathan, and D. L. Wright, Communica- tions Biology5, 459 (2022)

  59. [59]

    J. V. Rodrigues, S. Bershtein, A. Li, E. R. Lozovsky, D. L. Hartl, and E. I. Shakhnovich, Proceedings of the National Academy of Sciences113, E1470 (2016)

  60. [60]

    T. N. Nguyen, C. Ingle, S. Thompson, and K. A. Reynolds, Nature Communications15, 3351 (2024)

  61. [61]

    S. G. Das, S. O. Direito, B. Waclaw, R. J. Allen, and J. Krug, eLife9, e55155 (2020)

  62. [62]

    S. G. Das, J. Krug, and M. Mungan, Physical Review X 12, 031040 (2022)

  63. [63]

    S. G. Das, M. Mungan, and J. Krug, Proceedings of the National Academy of Sciences122, e2422520122 (2025)

  64. [64]

    R. R. Regoes, C. Wiuff, R. M. Zappala, K. N. Garner, F. Baquero, and B. R. Levin, Antimicrobial Agents and Chemotherapy48, 3670 (2004)

  65. [65]

    Ziegler, J

    C. Ziegler, J. Martin, C. Sinner, and F. Morcos, Nature Communications14, 2222 (2023)

  66. [66]

    I. K. Gaszek, M. S. Yildiz, D. Meng, J. A. de la Paz, S. M. Alvarez, F. Morcos, M. Lin, and E. Toprak, (2025), bioRxiv:2025.07.08.663783

  67. [67]

    C. L. David, E. E. Howell, M. F. Farnum, J. E. Vil- lafranca, S. J. Oatley, and J. Kraut, Biochemistry31, 9813 (1992). 11

  68. [68]

    Le Nagard and O

    H. Le Nagard and O. Tenaillon, inAdvances in Network Complexity, Vol. 4 (Wiley, 2013) pp. 41–61

  69. [69]

    Schneemann, B

    H. Schneemann, B. De Sanctis, and J. J. Welch, Cold Spring Harbor Perspectives in Biology16, a041442 (2024)

  70. [70]

    Blanquart, G

    F. Blanquart, G. Achaz, T. Bataillon, and O. Tenaillon, Evolution68, 3537 (2014)

  71. [71]

    Blanquart and T

    F. Blanquart and T. Bataillon, Genetics203, 847 (2016)

  72. [72]

    Hwang, S.-C

    S. Hwang, S.-C. Park, and J. Krug, Genetics206, 1049 (2017)

  73. [73]

    T. D. Schneider, Nucleic Acids Research28, 2794 (2000)

  74. [74]

    F. W. Stearns, Genetics186, 767 (2010)

  75. [75]

    Levenberg, Quarterly of Applied Mathematics2, 164 (1944)

    K. Levenberg, Quarterly of Applied Mathematics2, 164 (1944)

  76. [76]

    D. W. Marquardt, Journal of the Society for Industrial and Applied Mathematics11, 431 (1963)

  77. [77]

    Attouch and B

    H. Attouch and B. F. Svaiter, SIAM Journal on Control and Optimization49, 574 (2011)

  78. [78]

    M. K. Transtrum and J. P. Sethna, (2012), arXiv:1201.5885 [physics.data-an]

  79. [79]

    Nesterov,Lectures on Convex Optimization, Springer Optimization and Its Applications, Vol

    Y. Nesterov,Lectures on Convex Optimization, Springer Optimization and Its Applications, Vol. 137 (Springer International Publishing, Cham, Switzerland, 2018)

  80. [80]

    D. P. Bertsekas,Reinforcement Learning and Opti- mal Control, 2nd ed. (Athena Scientific, Belmont, Mas- sachusetts, 2019)

Showing first 80 references.