{"id":"d67a1ca3-d81e-498b-a6f1-6351cbe38089","arxiv_id":"2411.15293","paper_version":3,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"A rigorous linear stability analysis shows the mean-field survival model for zebrafish stripes exhibits Turing instability exactly when birth rate is below death rate, interaction range h is at least 3, and cell count N is sufficiently large, with stripe widths between 2h and 4h.","lead":"This paper analyzes a discrete ring model of zebrafish stripe formation and proves, with explicit parameter conditions, that the model undergoes a Turing instability: stripes emerge when melanophore death outpaces birth, the interaction range is at least three cells, and the ring is large enough. It matters because it shows that analyzing the discrete equations directly yields different predictions from the usual continuum PDE limit.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified; the linear-stability argument is self-consistent and the main theorems are supported.","rationale":"The reader's weakest assumption is the mean-field closure: the ODE system assumes no spatial correlations in the stochastic transitions. This is a legitimate limitation for the paper's biological framing, but the central claim is a theorem about the ODE system itself, and that theorem is proven rigorously. I checked the key algebraic steps: the coexistence equilibrium, the Fourier-mode Jacobian, the trace and determinant formulas, and the equivalence between F<0 and instability. The proofs of Lemma 1 and Theorems 1-2 are internally consistent, and the numerical simulations support the linear predictions. The only substantive caveat is that the mean-field closure is not validated against the stochastic model, which would matter if the title is read as a claim about the original survival process. That caveat was already noted by the reader and does not change the verdict on the mathematical results.","tokens_in":31474,"tokens_out":23203,"duration_ms":221210,"concrete_test":"Run the original stochastic survival model from Konow et al. (2021) at parameters inside P, e.g., b=4, d=16, h=4, N=23, and compare the emergence of spatial modes with the ODE predictions; if no pattern forms, the biological interpretation of the title would need qualification, though the ODE theorem stands.","verdict_should_be":"UNCHANGED","load_bearing_attack":"No load-bearing concern with the central mathematical claim. I re-derived the Jacobian, trace, and determinant in Section 2.1 and confirmed that linear instability is controlled by the sign of F in Eq. (5); the F<0 criterion and the containment of R(b,d) in Lemma 1 are internally consistent, and the density arguments in Theorems 1 and 2 are valid. The mean-field closure noted in Section 1.1 is a real limitation for the biological extrapolation, but it is explicitly disclosed and does not undermine the ODE theorem.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper performs a linear stability analysis of the mean-field survival model (1), a 2N-dimensional coupled ODE system on a ring describing xanthophore and melanophore probabilities, with nearest-neighbor inhibitory coupling and distance-h excitatory coupling. The main results characterize when the nontrivial homogeneous equilibrium is unstable to spatially periodic perturbations: Theorem 1 states that for (b,d) in P={0<b<d}, instability occurs for sufficiently large h and N, with the most unstable wavenumbers lying between ceil(N/4h) and ceil(N/2h), and that the equilibrium is linearly stable outside P; Theorem 2 shows that h>=3 is necessary and sufficient for instability to be possible for large N; Theorem 3 and the corollaries give vertical and slant asymptotes of the neutral stability curves, including the N->infinity limit. The paper also compares the ODE predictions with the earlier PDE limit, reports numerical simulations consistent with the linear analysis, and draws qualitative comparisons with zebrafish stripe widths and melanophore projection lengths.","tokens_in":31507,"tokens_out":19037,"duration_ms":176219,"significance":"If the main theorems are correct, this is a valuable contribution: it provides a complete, essentially parameter-free linear-stability characterization of a coupled ODE model that is not obtained by discretizing a PDE, and it shows that the N->infinity ODE limit differs qualitatively from the previously studied PDE limit. The derivations are self-contained and the linear algebra is checked in detail: the sign of the determinant is controlled by the explicit function F in Eq. (5), Lemma 1 characterizes its negativity region, and the density arguments in Theorems 1 and 2 are valid. The numerical simulations corroborate the predicted unstable modes. The biological discussion is appropriately cautious, although the mean-field closure assumption, disclosed in Section 1.1, is not tested against stochastic simulations of the underlying survival model; this limits the strength of the biological extrapolation but does not affect the ODE theorem.","major_comments":[{"comment":"","section":"Theorem 3, Eq. (10c); Corollary 1, Eq. (12c); Corollary 2, Eq. (16d)"},{"comment":"","section":"Corollary 2, Eq. (16b)"}],"minor_comments":[{"comment":"","section":"Proof of Theorem 2"},{"comment":"","section":"Proof of Theorem 3"},{"comment":"","section":"Figure 4 and Figure 5 captions"},{"comment":"","section":"Section 1, affiliation"}],"recommendation":"major_revision","confidential_remarks":"The mathematical core of the paper is sound, and I verified the trace/determinant computations and the sign criterion F<0. The main issue is that two displayed asymptote formulas in the main theorem statements (Theorem 3 and Corollaries 1-2) are inconsistent with the proofs and with the numerical values quoted in the text. These are straightforward sign/transposition typos, but because they appear in the central statements of the paper, they need to be fixed before publication. I do not see any load-bearing error in the instability criteria or the wavenumber bounds, so after correction I would support acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a rigorous and self-contained linear stability analysis of a discrete ring ODE model for zebrafish stripe formation, and the main theorems hold up. It deserves a serious referee.\n\nThe genuinely new part is the direct LSA of the mean-field survival model, rather than its PDE limit. The paper gives necessary and sufficient conditions for Turing instability in terms of b, d, h, N (Theorems 1 and 2), explicit bounds on the unstable wavenumbers, and the asymptotics of the neutral stability curves (Theorem 3). I re-derived the Jacobian, trace, and determinant for the coexistence equilibrium and found them correct; the F<0 instability criterion and Lemma 1's characterization of the instability region are internally consistent. The numerical simulations agree with the linear predictions, and the wavelength bounds match observed stripe widths in a qualitative but non-circular way—the biological comparisons are derived bounds, not fitted parameters.\n\nThe soft spots are real but not load-bearing. The mean-field closure is stated in Section 1.1—no spatial correlations in the stochastic transition events—but the paper never tests it against stochastic simulation. If correlations matter at biological scales, the predicted instability may not appear in the actual process the model is meant to describe. This caveat is disclosed, but a referee should ask how robust the instability is to correlations. Second, there is no shared code or data, which makes the numerical part harder to reproduce, though the analytic results are checkable by hand. Third, the 'PDE limit misleads' narrative is tied to the specific second-order Taylor expansion used in Konow et al. 2021; the paper acknowledges this, but the case-study conclusion is only as strong as that particular continuum approximation. That is a minor point, not a flaw in the ODE analysis.\n\nBottom line: the mathematics is solid, the novelty is real, and the limitations are honestly disclosed. This is a paper for mathematical biologists working on discrete Turing models and for zebrafish researchers who want a rigorous account of this particular mean-field model. I would cite it and would bring it to a reading group. Send it to peer review; ask the referees to focus on the mean-field closure and on requesting simulation or code availability.","headline":"Rigorous, self-contained Turing analysis of a discrete cell model; the main theorems check out and the biological claims are honest, though the mean-field closure deserves scrutiny.","tokens_in":32025,"tokens_out":2396,"would_cite":true,"duration_ms":23494,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["34C23","92C15"],"pacs":[],"model":"deepseek-v4-flash","headline":"A ring of coupled cell equations provably reproduces zebrafish stripe formation through a Turing instability, with exact parameter bounds.","keywords":["Turing instability","mean-field survival model","coupled ODE system","zebrafish stripe formation","linear stability analysis","discrete Fourier transform","Chebyshev polynomials","pattern formation"],"falsifier":"Run the stochastic survival model (system 2) at parameters inside $P$ with $h\\ge 3$ and $N$ large enough for Theorem 1, starting from a uniform random state; if the homogeneous state persists and no periodic stripe mode grows, the mean-field closure—and with it the paper's central prediction—would be falsified.","tokens_in":31244,"feed_emoji":"🐟","tokens_out":7073,"duration_ms":65305,"temperature":0.7,"pith_summary":"This paper argues that the black and yellow stripes of zebrafish can be generated by a Turing instability—the classic mechanism in which a uniform state is destabilized by spatially periodic perturbations—in a ring of coupled ordinary differential equations. The model tracks the expected presence of melanophores and xanthophores at each cell location, with short-range inhibition and long-range promotion at distance h. The main theorems give a complete parameter classification: if the melanophore death rate d exceeds the melanophore birth parameter b (the set $P=\\{0<d,\\; 0<b<d\\}$), then for melanophore projections at least three cells long and for sufficiently many cells N, the uniform coexistence state is unstable to periodic patterns, and the most unstable stripe wavenumber lies between $N/4h$ and $N/2h$. Outside $P$, or when $h<3$, the uniform state is linearly stable. The analysis is performed directly on the finite-N ODE system, and the paper shows that the continuum PDE limit proposed earlier predicts instability in parameter regions where the finite model is stable.","feed_headline":"Zebrafish stripes emerge from a small ring of coupled cells","feed_subtitle":"New proof fixes the exact parameter window: death rate above birth rate and projections at least three cells long.","key_machinery":"The load-bearing object is the linearization of the ring ODE system in discrete Fourier space. Because the ring is translation invariant, the $2N$-dimensional Jacobian splits into $N$ independent $2\\times 2$ blocks $\\hat{L}_k$ labeled by Fourier mode $k$, whose trace is always negative and whose determinant is $F(x_k, T_h(x_k), b, d)/((b+d)(b+2d))$. The sign of $F$ therefore decides stability: $F<0$ means an eigenvalue with positive real part and a Turing bifurcation. Lemma 1 establishes the geometry of the set where $F<0$ in the $(x,y)$-plane, including that it is nonempty exactly for $(b,d)\\in P$ and that its left edge has infimum $1/3$, which is what forces the threshold $h\\ge 3$. The proof then uses properties of Chebyshev polynomials—their rightmost local minimum at $x=\\cos(\\pi/h)$ approaching $x=1$ for large $h$, and the density of the points $x_k$ for large $N$—to turn the geometry into the parameter conditions of Theorems 1 and 2.","core_discovery":"On the paper's own terms, the central discovery is that the mean-field survival model has a Turing bifurcation whose onset is exactly encoded by the sign of a single polynomial, $F(x_k, T_h(x_k), b, d)$, where $x_k=\\cos(2\\pi k/N)$ and $T_h$ is the $h$-order Chebyshev polynomial. The coexistence equilibrium $(d/(b+d), b/(b+2d))$ is linearly unstable exactly when some Fourier mode $k\\neq 0$ makes $F<0$. Theorem 1 proves that for $(b,d)\\in P$, sufficiently large $h$ and $N$ guarantee such a mode exists, and the fastest-growing mode satisfies $\\lceil N/4h\\rceil \\le k \\le \\lceil N/2h\\rceil$; for $(b,d)\\notin P$ the equilibrium is linearly stable. Theorem 2 proves that $h\\ge 3$ is necessary and sufficient for instability to be possible for large $N$, and Theorem 3 gives explicit vertical and slant asymptotes for the neutral stability curves. The smallest ring that can produce stripes is $N=6$ with $h=3$. These results imply that melanophore and xanthophore stripes always form out of phase with each other, and that stripe width should be between two and four projection lengths, matching the observed values in zebrafish.","pith_inferences":["Editorial inference: the same discrete-Fourier-plus-Chebyshev method should transfer to any ring or network model whose coupling is a translation-invariant average, with the curve $y=T_h(x)$ replaced by the ratio of the Fourier symbols of the short-range and long-range coupling operators.","Editorial inference: the discrepancy with the PDE limit suggests a testable mathematical question—whether a higher-order Taylor expansion or a nonlocal integral kernel converges to the ODE instability region as $N\\to\\infty$; the paper states no such general result exists.","Editorial inference: the mean-field assumption of no spatial correlation could be checked by simulating the stochastic survival model directly; if correlations suppress the predicted modes, the deterministic predictions would fail precisely at biologically relevant cell densities.","Editorial inference: the model's stripe-count growth prediction could be compared with time-lapse observations of regenerating or growing fin stripes, where the wavenumber should increase in stepwise fashion as circumference crosses multiples of the preferred wavelength."],"forward_implications":["If the theorem is right, any zebrafish large enough to have at least six cells around its circumference and projections at least three cells long lies in a parameter regime where the uniform pigment state is unstable, so stripes can form without cell movement or iridophores.","The bound $\\lceil N/4h\\rceil \\le k \\le \\lceil N/2h\\rceil$ fixes stripe width between two and four projection lengths; combined with observed stripe widths of 7–12 cell diameters, it selects the realistic range $4\\le h\\le 13$, $112\\le N\\le 192$.","Because $h=1,2$ never yield instability for any birth and death rates, the model makes a sharp biological prediction: melanophore projections shorter than three cell diameters cannot produce stripes by this mechanism.","The finite-$N$ ODE analysis and its $N\\to\\infty$ limit differ qualitatively from the earlier PDE continuum model, which predicts instability for arbitrarily large birth parameter $b$; the ODE analysis says the unstable region is confined to $0<b<d$.","For a growing fish, the bounds imply stripe width stays near a fixed chemical wavelength while the number of stripes increases roughly in proportion to circumference."],"supporting_citations":[{"why":"Supplies the survival model, its mean-field ODE and PDE reductions, and the numerical analysis that the present work extends with a direct linear stability analysis.","marker":"[26]"},{"why":"Establishes the Turing-instability framework and the ring geometry used throughout.","marker":"[1]"},{"why":"Provides the ablation observations motivating the death and growth rates and the stripe-width measurements used for biological comparison.","marker":"[22]"},{"why":"Supplies the measured melanophore projection lengths that set the thresholds $h\\ge 3$ and stripe-width bounds.","marker":"[23]"},{"why":"The promotion model whose parameter simplifications are inherited by the mean-field survival model.","marker":"[25]"},{"why":"Agent-based model used as a biological comparison for stripe width and cell diameter.","marker":"[27]"},{"why":"Provides staging data on stripe count versus fish length used in the growth prediction.","marker":"[30]"}],"fun_headline_variants":["Turing instability proven in mean-field zebrafish stripe model","Exact conditions found for stripes in coupled-cell model","Zebrafish stripes: math pins down parameter window","Mean-field model shows Turing pattern at ring size six","Proof: stripe formation needs death rate above birth rate"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the expected values at neighboring cell locations are uncorrelated, so the stochastic survival model closes exactly to the ODE system; if spatial correlations are significant, the predicted Turing stripes may not occur in the stochastic process.","fun_headline_variants_meta":{"raw":{"variants":["Turing instability proven in mean-field zebrafish stripe model","Exact conditions found for stripes in coupled-cell model","Zebrafish stripes: math pins down parameter window","Mean-field model shows Turing pattern at ring size six","Proof: stripe formation needs death rate above birth rate"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000465,"raw_usage":{"total_tokens":2341,"prompt_tokens":984,"completion_tokens":1357,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":600,"completion_tokens_details":{"reasoning_tokens":1281}},"tokens_in":600,"tokens_out":1357,"duration_ms":10232,"temperature":1.0,"reasoning_tokens":1281,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:31:15.193403+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the stochastic survival model (system 2) at parameters inside $P$ with $h\\ge 3$ and $N$ large enough for Theorem 1, starting from a uniform random state; if the homogeneous state persists and no periodic stripe mode grows, the mean-field closure—and with it the paper's central prediction—would be falsified.","supporting_citations":[{"cited_title":"Konow, Z","cited_arxiv_id":null,"evidence_quote":"Supplies the survival model, its mean-field ODE and PDE reductions, and the numerical analysis that the present work extends with a direct linear stability analysis."},{"cited_title":"Hamada, M","cited_arxiv_id":null,"evidence_quote":"Supplies the measured melanophore projection lengths that set the thresholds $h\\ge 3$ and stripe-width bounds."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides staging data on stripe count versus fish length used in the growth prediction."}],"review_version":1}