{"id":"c0c908ed-a7eb-4e36-b7a6-783966c9cb66","arxiv_id":"1908.02670","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A rare dormant phenotype like bacterial persisters stabilizes competitive microbial communities, while abundant phenotypes destabilize them, going beyond the classical May stability paradigm.","lead":"This paper shows that ignoring the different states, or phenotypes, that microbial species can switch between leads to wrong stability predictions in ecological models. It finds that a rare, dormant state such as bacterial persisters can make a whole community more stable, which matters for understanding infections and the human microbiome.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The exponential stabilization in Eq. (9) hinges entirely on the unproven inequality q'>q, whose proof is deferred to the missing Supplemental Material; the main-text justification after Eq. (6) is ambiguous and unverified.","rationale":"The reader's weakest_assumption identifies q'>q as the key unproven inequality, and my independent reading of the main text agrees. The paper's derivation of Eq. (9) is a mean-field calculation whose exponent is directly proportional to q'−q; if that difference is not positive, the central claim reverses. The main text offers only a terse, ambiguous justification and defers the proof to Supplemental Material [19], which is not available in the arXiv version. This is not an accusation of error; the mechanism is plausible, and a standard perturbative calculation (using the block-triangular J0) likely delivers q'>q. But because the entire surprising conclusion rests on this point, and because the proof is absent from the submitted text, the appropriate verdict remains CONDITIONAL: acceptance should require the Supplemental Material proof and, ideally, a numerical confirmation of q'>q over the sampling distribution used in Fig. 2. I therefore see no reason to change the reader's verdict, only to emphasize that this specific inequality is the gatekeeper for the paper's main result.","tokens_in":8189,"tokens_out":14331,"duration_ms":159259,"concrete_test":"Compute the first-order eigenvalue correction for J* = J0 + εJ1 using the block-triangular structure of J0 in Eq. (6): for an eigenvalue ν0 of K0 with right eigenvector (u,0) and left eigenvector (v, λ/(λ+ν0)v), the correction is δν = −κ ν0/(λ+ν0). Then, drawing many random K0 and J1 from the sampling distribution in [19], collect eigenvalues with |Re(ν0)| ≤ ε and compare the fractions whose perturbed real part becomes negative for J* versus K* = K0 + εK1 with K1 = B0Eκ/λ from Eq. (7). If the empirical q'−q is not positive, Eq. (9) is false; if it is positive, the missing proof in the Supplemental Material should be supplied and checked.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative claim, p'/p ~ exp[2c(q'−q)ε√N] in Eq. (9), rests on the assertion q'>q, made just after Eq. (6). This is the only step that distinguishes the stability of the full phenotype-resolving model from that of the averaged model in the asymptotic regime. In the arXiv v2 text, the proof is not provided; the sentence 'Eq. (6) shows that J1 acts on K0 as the negative definite matrix −κI, so q′ > q [19]' defers all details to the Supplemental Material, which is absent from the arXiv posting. The assertion is not immediate from Eq. (6): J1 is a 2N×2N block matrix containing off-diagonal blocks (−B0F and κI) and a diagonal block (γ−G·B0)I, so it does not literally act on K0. A first-order perturbation calculation around the block-upper-triangular J0 is needed to show that the relevant eigenvalue shifts are negative and that q′ exceeds q. If q′ were ≤ q, the ratio in Eq. (9) would fail to grow, and the claimed stabilizing effect of a rare phenotype would not follow for the stated model. The numerics in Fig. 2 are consistent with q'>q but do not prove it, and the heuristic spectral argument in Eq. (8) also assumes independent eigenvalue stabilization events. The deferred proof of q'>q is therefore the single most load-bearing unverified step in the paper.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the linear stability of Lotka–Volterra models of microbial communities in which each species has multiple phenotypes. It first shows numerically that an averaged model without subpopulations can disagree with the full phenotype-resolving model about stability. It then analyzes a persister-type model with a small parameter ε, expands the Jacobians of the full and averaged models, and uses a random-matrix argument to derive the asymptotic probabilities that each model is stable. The central quantitative claim is Eq. (9): when switching to a rare phenotype, the full phenotype-resolving model is exponentially more likely to be stable than the corresponding averaged model, p′/p ~ exp[2c(q′−q)ε√N]. The paper also presents numerical simulations of large perturbations, showing that averaged and full models can predict different outcomes.","tokens_in":8479,"tokens_out":9186,"duration_ms":104807,"significance":"If the main result holds, it is a substantive contribution to ecological stability theory: it identifies a concrete mechanism by which subpopulation structure, specifically a rare persister phenotype, can stabilize a competitive community, partially overturning May's instability paradigm. A strength of the paper is that it works with explicit dynamical equations and a model-specific Jacobian structure rather than with a purely abstract random Jacobian, and the numerical exploration is extensive, including up to 10^8 random systems. The falsifiable prediction in Eq. (9) is clearly stated. However, the quantitative claim rests on a heuristic spectral calculation and, most importantly, on an inequality q′ > q whose proof is deferred to a Supplemental Material not included in the arXiv version. The paper's main conclusion is therefore not fully verifiable from the submitted text.","major_comments":[{"comment":"The central claim p′/p ~ exp[2c(q′−q)ε√N] depends entirely on the inequality q′ > q, but the proof is not presented in the main text and the cited Supplemental Material [19] is not included in the arXiv posting. The sentence 'Eq. (6) shows that J1 acts on K0 as the negative definite matrix −κI, so q′ > q [19]' is not a proof, because J1 is a 2N×2N block matrix and does not literally operate on K0. The first-order shift of an eigenvalue of K0 under J0 + εJ1 involves the off-diagonal blocks J1_12 = −B0F and J1_21 = κI as well as the diagonal block −κI, and a standard perturbation calculation is needed to determine the sign of the correction. If q′ ≤ q, the ratio in Eq. (9) would not grow and the claimed stabilizing effect of the rare phenotype would not follow for the stated model. The supplemental derivation should be made available and its main steps summarized in an accessible way.","section":"After Eq. (6), Eq. (9)"},{"comment":"The binomial expression in Eq. (8) treats the N−1 non-outlier eigenvalues as independent trials with common probabilities ϖ and q. This independence and identical-distribution assumption is not established: the eigenvalues of the random matrix K0 are correlated, and the row scaling by B0 in B0E means the matrix does not have i.i.d. entries in the standard circular-law setting. The argument is explicitly heuristic, but the paper should state more clearly that Eq. (8)–(9) are heuristic asymptotic estimates rather than proven asymptotics, and either justify the independence step or provide a numerical test of the predicted exponential form beyond Fig. 2.","section":"Eq. (8)"},{"comment":"The numerical support for the quantitative form of Eq. (9) is weak. The simulations use ε = 0.01 and N ≤ 10, so ε√N ≤ 0.032; the predicted log-ratio is 2c(q′−q)ε√N, which is only of order 0.1 if c(q′−q) is O(1). The trend in Fig. 2 is consistent with the direction of Eq. (9), but the figure cannot distinguish an exponential growth from other monotone trends at such small values of the exponent. This is not a fatal flaw, but the text should avoid claiming that Fig. 2 confirms the exponential asymptotic.","section":"Fig. 2, Eq. (9)"}],"minor_comments":[{"comment":"The Supplemental Material is cited as '[url to be inserted]' and is absent from the arXiv version; this needs to be finalized and posted with the manuscript.","section":"Reference [19]"},{"comment":"The notation 'B0E' should be defined explicitly, since B0 is a vector and the product must be understood as diag(B0)E (or equivalently elementwise multiplication); otherwise Eq. (6) and Eq. (7) are ambiguous.","section":"Notation after Eq. (6)"},{"comment":"The closing sentence 'our conclusions apply not only to competitive, but also to more general interactions' is not demonstrated by the analysis, which relies on the specific block structure of the competitive Lotka–Volterra Jacobian; this sentence should be phrased as a conjecture or conjecture plus supporting heuristic.","section":"Final paragraph"},{"comment":"The caption lists 'Eq. (4)', 'Eq. (5), asympt.', and 'Eq. (5), exact' only in the legend text; for reader clarity, the caption should also state which markers correspond to which model, since the main text refers to large and small markers.","section":"Fig. 2 caption"},{"comment":"In Fig. 4(c) the category '#' combines numerical failures and cases where different species die out; the text says these are a small number of cases, but the figure does not separate them, making the result harder to evaluate.","section":"Fig. 4(c)"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely acceptable once the proof of q′ > q is available and the spectral-heuristic caveats are stated more carefully. Since the arXiv version lacks the Supplemental Material, the central inequality cannot currently be verified; I would ask the editor to ensure the supplemental calculation is included in the review version and, ideally, summarized in the main text."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, this paper makes a genuinely new claim: for a two-phenotype Lotka–Volterra model with switching and separated time scales, a rare persister phenotype can linearly stabilize a community in a way the corresponding averaged model does not capture, and the ratio of stability probabilities grows as exp[2c(q'−q)ε√N] (Eq. 9). That asymptotic formula, and the warning that implicit averaging gives wrong stability predictions, are not in the previous literature; Maynard et al. 2019 uses a different model without switching. Second, the central step of the derivation — the inequality q' > q that makes the exponent positive — is asserted in the text and its proof is deferred to the Supplemental Material, which is absent from the arXiv posting. That is the load-bearing wall, and I cannot verify the claim from the main text alone.\n\nWhat the paper does well: the consistency conditions, Eqs. (3), give a clean, explicit way to associate an averaged model to the full system, and the numerics in Fig. 2 track the asymptotic prediction well. The asymptotic expansion of the Jacobian (Eqs. 6–7) is standard, and the spectral heuristic (circular core plus outlier, eigenvalues shifted by O(ε)) is plausible. The authors are honest that the argument is heuristic — c and q are O(1) constants left uncomputed — and they don't overclaim biological scope.\n\nSoft spots, in order of seriousness. (1) The proof of q' > q is not in the main text. The sentence after Eq. (6) says J1 acts on K0 as −κI, but J1 has off-diagonal blocks (−B0F and κI); it does not literally act like a negative multiple of the identity on the interesting subspace. A perturbative argument is needed to show the near-neutral eigenvalues receive negative shifts. If q' ≤ q, Eq. (9) reverses and the claimed stabilization disappears. The numerics are consistent with q' > q but don't prove it. (2) The spectral calculation treats the near-neutral eigenvalues as independent and ignores correlations between J1 and eigenvectors of K0; that's a heuristic, not a theorem. (3) The averaged model is not unique — Eqs. (3) fix it uniquely given the equilibrium, but that is still a modeling choice, and the paper should discuss robustness to that choice. (4) No code or data archive is listed.\n\nMy overall judgment: the core derivation and numerics make a believable case that the rare-phenotype stabilization is real for this model class, but the paper cannot be fully accepted without the missing SM and a referee's check of the q' > q step. That step is exactly the kind that could quietly fail.\n\nThis is for theoretical ecologists and mathematical biologists working on random-matrix stability and microbial persistence. It deserves a serious referee — I would send it out rather than desk reject — with instructions that the SM must be supplied and that the q' > q proof gets a careful, line-by-line check. If it passes, it's a nice short contribution to the stability literature.","headline":"A novel asymptotic result showing rare phenotypes can stabilize competitive communities, but the key inequality q' > q is deferred to missing supplemental material and needs referee verification.","tokens_in":8990,"tokens_out":3048,"would_cite":true,"duration_ms":31610,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["92D25","92D40","15B52"],"pacs":[],"model":"deepseek-v4-flash","headline":"A community with a rare persister phenotype is much more likely to be stable than the averaged model predicts.","keywords":["phenotypic switching","subpopulation structure","ecological stability","Lotka–Volterra","bacterial persisters","random matrix theory","microbiome","linear stability"],"falsifier":"Compute the stabilization probabilities $q$ and $q'$ directly from Eq. (6) for many random parameter draws and check whether $q' > q$ uniformly; alternatively, measure the ratio $P(J_* \\text{ stable})/P(K_* \\text{ stable})$ in simulations of Eqs. (5) at fixed $\\varepsilon$ and increasing $N$ and test whether it grows like $\\exp[2c(q'-q)\\varepsilon\\sqrt{N}]$ rather than falling.","tokens_in":7970,"feed_emoji":"🦠","tokens_out":9301,"duration_ms":88929,"temperature":0.7,"pith_summary":"This paper claims that the standard practice of averaging over the phenotypes of a microbial species before analyzing ecological stability gives the wrong answers. Using a competitive Lotka–Volterra model in which each species has a normal phenotype and a rare persister phenotype, the authors derive a quantitative asymptotic comparison between the full phenotype-resolving model and the corresponding averaged model. The central quantitative result is that the probability that the full model is linearly stable divided by the probability that the averaged model is linearly stable grows like $\\exp[2c(q'-q)\\varepsilon\\sqrt{N}]$ as the number of species $N$ grows, so a rare, slow-growing phenotype has a stabilizing effect. Abundant phenotypes are instead destabilizing. The paper further shows that for large perturbations, such as antibiotic treatment, the full and averaged models frequently settle into different coexistence states.","feed_headline":"Rare persister phenotypes stabilize microbial communities","feed_subtitle":"Full phenotype-resolving models are exponentially more stable than averaged models as species count grows.","key_machinery":"The load-bearing machinery is the matched pair of Lotka–Volterra systems: the phenotype-resolving model (5), where each species has normal cells $B$ and persister cells $C$ with switching rates $\\kappa,\\lambda$ and a small parameter $\\varepsilon$ separating growth and competition, and the averaged model (1) defined from it by the consistency conditions (3). The argument runs entirely through the asymptotic expansion of the two Jacobians, $J_* = J_0 + \\varepsilon J_1 + O(\\varepsilon^2)$ and $K_* = K_0 + \\varepsilon K_1 + O(\\varepsilon^2)$. Because $J_0$ is block-upper-triangular, the full and averaged models share leading-order eigenvalues; the key object is the first-order correction $J_1$, which acts on $K_0$ as the negative definite matrix $-\\kappa I$ and therefore stabilizes near-neutral eigenvalues more often than the averaged correction $K_1 = B_0 E \\kappa/\\lambda$. The binomial identity (8) turns that per-eigenvalue stabilization into the exponential ratio (9).","core_discovery":"The paper establishes that linear stability of microbial communities cannot be reliably inferred from models that lump subpopulations into a single averaged state. Concretely, for the phenotype-resolving model (5) with normal cells $B$ and persisters $C$ and small switching parameters, the Jacobian splits as $J_* = J_0 + \\varepsilon J_1 + O(\\varepsilon^2)$, while the averaged model has $K_* = K_0 + \\varepsilon K_1 + O(\\varepsilon^2)$; the leading-order spectra agree, but the first-order correction $J_1$ contains a negative definite term $-\\kappa I$ that makes near-neutral eigenvalues of the full model more likely to be stable than those of the averaged model. Summing this per-eigenvalue stabilization over the $N-1$ non-outlier eigenvalues yields the exponential ratio in Eq. (9): the full model becomes overwhelmingly more likely to be linearly stable as $N$ grows. This contradicts the expectation, inherited from random-matrix stability arguments, that resolving more degrees of freedom should make a community less stable. For large perturbations such as antibiotic exposure, the full and averaged models also differ in which species survive, even when both are stable at the coexistence equilibrium.","pith_inferences":["Beyond the paper: the exponential factor $\\exp[2c(q'-q)\\varepsilon\\sqrt{N}]$ implies that even very small persister fractions can have a disproportionately large stabilizing effect in large communities; this could be tested by varying $\\varepsilon$ and $N$ in synthetic microbial communities.","Beyond the paper: the same stabilization mechanism may operate for any rare, slow-growing, stress-tolerant subpopulation, not only antibiotic persisters; dormant spores or slow-growing variants would be natural candidates.","Beyond the paper: since the proof of $q'>q$ is deferred to the Supplemental Material, a direct numerical check of the inequality over a wide parameter range would be the fastest way to confirm the mechanism."],"forward_implications":["Stability conclusions drawn from averaged species models should be rechecked when subpopulation structure is present, since resolved and averaged models disagree systematically.","Communities with switching to a rare phenotype such as persisters become more likely to be linearly stable as the number of species grows, relative to what the averaged model predicts.","Adding abundant phenotypes is destabilizing because it effectively increases the number of interacting species.","Large perturbations, like antibiotic treatments, can drive averaged and phenotype-resolving models to different coexistence states, so predicting microbial community resilience requires resolving subpopulations.","The stabilizing effect only needs generic properties of the spectral distribution, so it may apply beyond competitive interactions, although the explicit proofs in the paper are for the competitive case."],"supporting_citations":[{"why":"Supplies the classical random-matrix stability criterion for large ecosystems that the paper's stability results overturn.","marker":"[1]"},{"why":"Provides the stability criteria and eigenvalue-distribution assumptions used to argue about the circular core of the spectrum.","marker":"[3]"},{"why":"Introduces the resource-explicit microbial community model whose phenotype-resolving variant is used to test that abundant phenotypes destabilize.","marker":"[6]"},{"why":"Documents competitive interactions as stabilizers in the human microbiome, the applied context the results speak to.","marker":"[10]"},{"why":"Establishes the biological persistence phenotype and its role in surviving antibiotic stress, motivating the rare-phenotype model.","marker":"[24]"},{"why":"Gives the single-species phenotypic switching model that the paper generalizes to $N$ competing species.","marker":"[25]"},{"why":"Supplies the empirically measured switching-rate separation ($\\varepsilon$ about $10^{-5}$) that defines the asymptotic regime.","marker":"[29]"},{"why":"Provides the circular law for random-matrix eigenvalue distributions used to estimate how often eigenvalues lie within order $\\varepsilon$ of the imaginary axis.","marker":"[30]"},{"why":"Recent related claim of stabilizing phenotypic variability without switching, which the paper contrasts as a different mechanism.","marker":"[32]"}],"fun_headline_variants":["Persister rarity flips microbial stability predictions","Phenotype resolution reveals rare cells stabilize communities","Subpopulation-aware models show persisters are stabilizers","Microbial stability emerges from phenotype detail","Rare phenotypes: hidden stabilizers of microbial ecosystems"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The key inequality $q' > q$, asserted after Eq. (6), is the load-bearing premise; its proof is deferred to the Supplemental Material, which is not part of the arXiv version, and if it failed the exponential stabilization in Eq. (9) would be reversed.","fun_headline_variants_meta":{"raw":{"variants":["Persister rarity flips microbial stability predictions","Phenotype resolution reveals rare cells stabilize communities","Subpopulation-aware models show persisters are stabilizers","Microbial stability emerges from phenotype detail","Rare phenotypes: hidden stabilizers of microbial ecosystems"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000702,"raw_usage":{"total_tokens":3129,"prompt_tokens":869,"completion_tokens":2260,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":485,"completion_tokens_details":{"reasoning_tokens":2190}},"tokens_in":485,"tokens_out":2260,"duration_ms":20345,"temperature":1.0,"reasoning_tokens":2190,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:38:14.815418+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the stabilization probabilities $q$ and $q'$ directly from Eq. (6) for many random parameter draws and check whether $q' > q$ uniformly; alternatively, measure the ratio $P(J_* \\text{ stable})/P(K_* \\text{ stable})$ in simulations of Eqs. (5) at fixed $\\varepsilon$ and increasing $N$ and test whether it grows like $\\exp[2c(q'-q)\\varepsilon\\sqrt{N}]$ rather than falling.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the classical random-matrix stability criterion for large ecosystems that the paper's stability results overturn."},{"cited_title":"Allesina and S","cited_arxiv_id":null,"evidence_quote":"Provides the stability criteria and eigenvalue-distribution assumptions used to argue about the circular core of the spectrum."},{"cited_title":"Butler and J","cited_arxiv_id":null,"evidence_quote":"Introduces the resource-explicit microbial community model whose phenotype-resolving variant is used to test that abundant phenotypes destabilize."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents competitive interactions as stabilizers in the human microbiome, the applied context the results speak to."},{"cited_title":"Harms, E","cited_arxiv_id":null,"evidence_quote":"Establishes the biological persistence phenotype and its role in surviving antibiotic stress, motivating the rare-phenotype model."},{"cited_title":"Kussell, R","cited_arxiv_id":null,"evidence_quote":"Gives the single-species phenotypic switching model that the paper generalizes to $N$ competing species."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the empirically measured switching-rate separation ($\\varepsilon$ about $10^{-5}$) that defines the asymptotic regime."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the circular law for random-matrix eigenvalue distributions used to estimate how often eigenvalues lie within order $\\varepsilon$ of the imaginary axis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Recent related claim of stabilizing phenotypic variability without switching, which the paper contrasts as a different mechanism."}],"review_version":1}