{"id":"7aa60c8c-1700-4b21-9e65-15374730da10","arxiv_id":"1908.02532","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Spatial disorder in movement makes invasion fronts wander and has a nonmonotonic effect on genetic loss: weak disorder preserves diversity, strong disorder accelerates its loss.","lead":"This paper studies how random spatial variation in how fast individuals move affects a population wave and the genetic diversity it carries. It reports that weak environmental disorder slows genetic loss while strong disorder speeds it up, and it derives a wandering motion for the invasion front.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The genetic-loss claim rests on a false homogeneous baseline: under Eq. (13), the difference C2−C1 obeys a noiseless linear equation, so the less abundant population cannot go extinct deterministically; Fig. 3's ξ=0 curve is likely a numerical-threshold artifact.","rationale":"The reader's weakest_assumption identifies exactly the load-bearing point. The paper's headline novelty is environmental regulation of genetic loss, and the central comparison in Fig. 3 depends on a homogeneous-reference extinction that the paper's own deterministic PDE, Eq. (13), forbids. Because f(x) is common to both populations, the difference equation for Υ has no stochastic forcing; the only asymmetry is the initial condition, so the less abundant population must persist indefinitely in the continuum model. The apparent extinction in the ξ=0 baseline can only arise from the integer-multiple-of-10^-10 flooring in the finite-difference scheme, which is a numerical threshold effect, not a biological mechanism. If the proposed test confirms the absence of deterministic extinction, the nonmonotonic heterozygosity curves in Fig. 3 are measuring how noise interacts with this numerical threshold, not a genuine suppression or intensification of genetic loss. The front-wandering result (Eq. 11 and Fig. 2) is a separate and more defensible contribution with numerical support, so the paper contains a useful piece of work, but that does not rescue the abstract's central genetic-loss claim. The appropriate verdict for the current manuscript is therefore REJECT rather than CONDITIONAL, since the stated model contradicts the baseline on which the main claim is built; a corrected analysis would require introducing demographic stochasticity or another explicit extinction mechanism, which is a substantial reworking rather than a minor revision.","tokens_in":7107,"tokens_out":9218,"duration_ms":114648,"concrete_test":"Re-run the homogeneous ξ=0 case of Eq. (13) with the same parameters and initial ε used in Fig. 3, but with a standard finite-difference scheme that omits the 'flows are integer multiples of 10^-10' rounding; track min_x C1(x,t) over the same simulated horizon. Independently, integrate the linear Υ equation to confirm Υ remains positive. If C1 never reaches the numerical floor, the extinction baseline is false, and Fig. 3 must be reinterpreted as an artifact of the 10^-10 cutoff rather than a biological genetic-loss process.","verdict_should_be":"REJECT","load_bearing_attack":"The weakest load-bearing assumption is the homogeneous-environment baseline for genetic loss. The text asserts before Fig. 3 that 'for non-zero epsilon in an homogeneous environment we know that C1 will be extinct, due to its smaller initial value.' This contradicts Eq. (13). Since the environmental noise f(x) is the same for C1 and C2, the difference Υ=C2−C1 satisfies ∂_tΥ = ∂_x[D0(1+ξf)∂_xΥ] + RΥ(1−C), a linear equation with no stochastic source and with a nonnegative growth coefficient R(1−C). If Υ(x,0)=2ε>0, Υ remains strictly positive: behind the front, C→1, the growth term vanishes, and Υ only diffuses; ahead of the front, both populations grow with the same per-capita rate R(1−C), preserving the ratio. Thus C1=(C−Υ)/2 stays at approximately 0.5−ε and cannot go extinct. The only way the ξ=0 curve in Fig. 3(a) can show extinction is the numerical discretization rule that rounds flows to integer multiples of 10^-10, an artificial cutoff. If the baseline is an artifact, then the statement that weak disorder 'delays' extinction and the nonmonotonic heterozygosity in Fig. 3(b) lack a meaningful reference. Moreover, because f(x) multiplies both populations' diffusion identically, environmental disorder alone provides no differential force between the two populations; ordinary genetic drift from demographic stochasticity is absent from the model. The central genetic-loss claim is therefore not established by the reported equations and simulations.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a one-dimensional FKPP-type reaction-diffusion equation with a spatially fluctuating diffusion coefficient, D(x) = D0(1 + ξ f(x)), where f(x) is introduced as uniform noise in [-1,1]. The authors derive a Langevin description for the position of the invasion front and obtain an effective diffusion coefficient D_C = (1/8) ξ^2 sqrt(R D0) for front wanderings, which they support with numerical fits reproducing the scaling in ξ, R, and D0. The second part addresses genetic loss by considering two neutral populations with slightly different initial densities and analyzes how spatial disorder in diffusivity affects their coexistence. The paper claims that weak environmental fluctuations (small ξ) delay genetic loss, while strong fluctuations (large ξ) intensify it, and it interprets this as a potential explanation for regional differences in human genetic diversity and tumor heterogeneity.","tokens_in":7517,"tokens_out":3158,"duration_ms":38316,"significance":"If correct, the front-wandering result would be a useful quantitative contribution to the theory of reaction-diffusion fronts in heterogeneous media, and the genetic-loss result would be a striking and biologically relevant prediction. The numerical validation of the D_C scaling is, on its face, well executed: the paper reports scaling exponents and R² values, and the perturbation derivation follows a standard route. The genetic-loss claim, however, is the central novelty of the paper, and it rests on assumptions that contradict the model equations. Because the claimed regulatory role of environmental disorder on genetic loss is not established by the reported analysis, the overall significance of the paper is substantially diminished. The front-wandering part may be salvageable, but the genetic-loss claim as presented is not.","major_comments":[{"comment":"The noise prescription is internally inconsistent: the model defines f as 'uniform white noise in the range [-1,1]', but the ensemble average leading to Eq. (9) uses <f(x)f(y)> = δ(x-y), which is the correlation of Gaussian white noise with unit variance, not of bounded uniform noise. For uniform noise in [-1,1] one would have <f(x)f(y)> = (1/3)δ(x-y), changing the prefactor in Eq. (11). The numerical fits in Fig. 2 compare only scaling exponents and do not verify the prefactor 1/8. The derivation of D_C is therefore not a parameter-free quantitative prediction for the simulated model.","section":"Model and Invasion Front Wanderings, Eq. (2) and Eq. (9)"},{"comment":"The baseline for the genetic-loss claim is incorrect. The text states that for non-zero ε in a homogeneous environment C1 will be extinct 'due to its smaller initial value', but this contradicts Eq. (14b). The difference Υ = C2 - C1 satisfies a linear equation with diffusion and growth term R(1-C)Υ, with R(1-C) ≥ 0, so starting from Υ(x,0) = 2ε > 0 yields Υ(x,t) > 0 for all time. The minority population cannot go extinct in the deterministic continuum model. The ξ=0 extinction curve in Fig. 3(a) is therefore a numerical artifact, most plausibly the discretization rule that rounds all flows to integer multiples of 10^-10. The comparison showing that weak disorder 'delays' extinction has no valid baseline.","section":"Genetic loss, paragraph before Fig. 3"},{"comment":"The analytic derivation of noise-induced fluctuations in κ rests on an uncontrolled approximation. In Eq. (16) the authors neglect the term 2R∫ C Υ² dx by asserting that C≈1 only where Υ≈0 and C≈0 only where Υ≈1. In a traveling front, however, there is a finite interfacial region where C and Υ are both of order unity, and there is no separation of scales that makes the product negligibly small. Since this approximation is the basis for Eqs. (17)-(19) and for the claimed fluctuation <Δκ²>, the analytic support for the genetic-loss result is not established.","section":"Genetic loss, Eqs. (16)-(19)"}],"minor_comments":[{"comment":"The text writes f(x) as spatial noise but Eq. (2) uses f(x,t); the time dependence is not defined and is later dropped. Please clarify the notation.","section":"Model section and Eq. (2)"},{"comment":"The phrase 'uniform white noise in the range [-1,1]' is ambiguous: it could mean uncorrelated uniform spatial noise, but 'white noise' usually implies Gaussian delta-correlated noise. The normalization and correlation function of f should be stated precisely.","section":"Model section"},{"comment":"The caption uses 'adimensional time'; the correct term is 'dimensionless time'.","section":"Fig. 3 caption"},{"comment":"The paper does not provide a clear definition of the ensemble average notation; Fig. 2 uses <(X−Xbar)²> and <(C_F−C_Fbar)²> without specifying how the average over disorder realizations is taken or how many realizations are used. Reporting the number of realizations and error bars would strengthen the numerical claims.","section":"Throughout"}],"recommendation":"reject","confidential_remarks":"The front-wandering scaling may be sound, but the genetic-loss result is the paper's central biological claim and it appears to be an artifact of the numerical cutoff in the homogeneous baseline. The authors should be invited to resubmit if they can demonstrate that their model, without artificial round-off, exhibits genuine extinction and that the nonmonotonic dependence on ξ is robust to discretization and to a proper stochastic formulation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper has one genuinely useful piece—front wandering in a FKPP front with spatially quenched diffusivity noise—and one headline result that appears to rest on a numerical artifact. The claim that small disorder suppresses genetic loss and large disorder intensifies it is not supported by the model as written.\n\nWhat is new and good: the effective front diffusivity D_C ~ xi^2 sqrt(R D0) is a parameter-free scaling law, numerically confirmed over a wide range. The zero-mode projection method is standard, but applying it to quenched diffusivity noise is a fair new application, and the simulation support for the exponents is solid.\n\nThe soft spots: first, the noise prescription is inconsistent. The model defines f as bounded uniform noise on [-1,1], while the analytics use Gaussian delta-correlated white noise with <f(x)f(y)> = delta(x-y). Those are not interchangeable, and the difference affects the constant in D_C. The step from Eq. (10) to (11) is also not shown, so the stated 1/8 prefactor is not derived in any transparent way.\n\nSecond, and more serious: the genetic-loss section claims that in a homogeneous environment C1 goes extinct because of its smaller initial value. But Eq. (14b) for Upsilon = C2 - C1 is linear, with no stochastic source and a nonnegative growth term. If Upsilon starts at 2epsilon, it remains strictly positive; there is no deterministic extinction mechanism. The zeta=0 curve in Fig. 3 likely shows extinction only because the numerics round flows to integer multiples of 10^-10, an artificial cutoff. That makes the baseline for the 'weak disorder delays extinction' story false, and the nonmonotonic heterozygosity curve is measured against the wrong reference. The model also has no demographic stochasticity, so 'genetic drift' in the usual sense is absent; environmental disorder acting identically on both populations cannot break the symmetry between equal-fitness types.\n\nThe front-wandering result may survive additional work, but the paper's central biological claim does not. A proper stochastic branching process or an analytically justified extinction threshold could revive the genetic-loss story, but as submitted it is one good scaling law attached to a broken headline.\n\nRecommendation: send it to peer review if you want the authors to fix the baseline and re-examine the noise prescription. It deserves a serious referee, but I would expect heavy revision.","headline":"Front-wandering scaling may be real, but the genetic-loss headline rests on a false homogeneous-extinction baseline.","tokens_in":7978,"tokens_out":2655,"would_cite":false,"duration_ms":29541,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["92D15","92D25","35K57","60H15"],"pacs":[],"model":"deepseek-v4-flash","headline":"Spatial fluctuations in movement speed regulate genetic loss during range expansion, with weak disorder slowing it and strong disorder accelerating it.","keywords":["range expansion","genetic drift","reaction-diffusion equation","spatial disorder","invasion front","heterozygosity","stochastic diffusivity","FKPP equation"],"falsifier":"A direct test is to integrate the deterministic two-population PDE (13) with $\\xi=0$, equal fitness and no demographic noise, starting from $C_1=0.49$, $C_2=0.51$: if the minority population does not go extinct asymptotically, the paper's reference behavior for 'genetic loss' is invalid. A complementary experiment would track heterozygosity in a microbial or cell range expansion on substrates with controlled diffusivity variance, looking for the predicted non-monotonic response with a minimum of genetic loss near $\\xi\\approx0.2$.","tokens_in":6879,"feed_emoji":"🧬","tokens_out":10582,"duration_ms":102979,"temperature":0.7,"pith_summary":"This paper asks whether randomness in the environment—specifically, spatial fluctuations in how fast individuals move—can change the fate of genetic diversity during a biological invasion. The authors model an invading population with the Fisher-Kolmogorov-Petrovsky-Piskunov (FKPP) equation, replacing the constant diffusion coefficient by $D_0(1+\\xi f(x))$ with white-noise $f(x)$. They derive a Langevin equation for the front position whose effective diffusion coefficient is $D_C = \\frac{1}{8}\\,\\xi^2\\sqrt{RD_0}$, and they show numerically and analytically that the same noise controls genetic loss: weak fluctuations (small $\\xi$) delay the extinction of a minority population and raise heterozygosity at the front, while strong fluctuations accelerate extinction. If true, environmental disorder is a regulator of genetic drift in its own right, independent of fitness differences or population-size fluctuations.","feed_headline":"Weak disorder slows genetic loss; strong disorder speeds it up","feed_subtitle":"How environmental variability alone can shape genetic drift during invasion and in tumors.","key_machinery":"The argument is carried by a perturbative treatment of the stochastic reaction-diffusion equation in a comoving frame. Writing the density as $C(\\zeta,t)\\approx C_0(\\zeta+\\eta(t),t)+\\delta C_1(\\zeta,t)$, with $C_0$ the solution of the unperturbed linearized equation, the authors project the noise term onto the zero eigenfunction $e^{\\sqrt{R/D_0}\\,\\zeta}$ of the adjoint operator. This yields an equation of motion for the front displacement $\\eta(t)$ and, after ensemble averaging $\\langle f(x)f(y)\\rangle=\\delta(x-y)$, the effective diffusion coefficient $D_C$. For genetic loss, the key object is the difference field $\\Upsilon=C_2-C_1$ and its squared integral $\\kappa(t)=\\int \\Upsilon^2\\,dx$; a Green-function calculation gives $\\langle\\Delta\\kappa^2\\rangle$, and direct simulation of the two-population equation supplies the heterozygosity $H=\\langle\\int C_1C_2\\,dx\\rangle$ that exhibits the non-monotonic dependence on $\\xi$.","core_discovery":"The central discovery is that spatially random diffusivity acts non-monotonically on genetic loss during range expansion. In the paper's two-population model, starting from nearly equal densities $C_1=0.5-\\epsilon$ and $C_2=0.5+\\epsilon$, the minority population eventually disappears in a uniform environment; introducing disorder of amplitude $\\xi$ delays that loss for $\\xi\\lesssim0.2$ and speeds it up for $\\xi\\gtrsim0.2$. The front itself performs an unbiased random walk, so the invasion border wanders as $\\langle (X_F-\\bar X_F)^2\\rangle \\sim 2D_C t$ with $D_C=\\frac18\\xi^2\\sqrt{RD_0}$. The authors interpret this as evidence that regional differences in the intensity of genetic drift—such as those observed in human expansions and tumor mutational landscapes—can be produced by environmental heterogeneity alone.","pith_inferences":["Inference: a similar non-monotonic response should appear in two- and three-dimensional expansions, where the disorder variance replaces the one-dimensional $\\xi^2$; the transition point may shift with dimension and with the correlation length of the noise.","Inference: the model predicts an experimentally accessible control variable—substrate heterogeneity in microbial range expansions—so the claim can be tested by growing populations on surfaces with engineered patches of different motility and measuring heterozygosity.","Inference: if this mechanism operates in humans, genetic-loss rates should correlate with measurable environmental variance (terrain ruggedness, resource patchiness) along the expansion path, not just with distance from the origin.","Inference: because the theory only requires a fluctuating diffusion constant, any process that modulates dispersal—seasonality, climate variability, or tissue stiffness—could play the same regulatory role, broadening the relevance beyond the paper's examples."],"forward_implications":["An invasion front in a habitat with random motility performs an unbiased random walk, so the position variance grows linearly in time with coefficient $\\frac18\\xi^2\\sqrt{RD_0}$; this can be tested by tracking range edges in experiments or field data.","Weak environmental disorder (roughly $\\xi<0.2$) preserves genetic diversity at the advancing front by delaying the extinction of less abundant populations.","Strong disorder ($\\xi>0.2$) does the opposite, accelerating genetic loss and making the loss spatially patchy even when fitness is equal.","Observed regional differences in human genetic loss intensity could be explained by differences in environmental variability rather than by distinct demographic histories alone.","In tumors, clonal diversity (and hence mutational heterogeneity) may be shaped by the physical disorder of the surrounding tissue through the same mechanism."],"supporting_citations":[{"why":"Supplies the perturbative method for treating front dynamics in a fluctuating medium and the effective-diffusion formula that leads to Eq. (11).","marker":"[5]"},{"why":"Provides the perturbation ansatz $C=C_0(\\zeta+\\eta(t),t)+\\delta C_1$ used to derive the Langevin equation for the front displacement.","marker":"[20]"},{"why":"Defines the FKPP reaction-diffusion framework that the paper extends to a spatially random diffusion constant.","marker":"[2]"},{"why":"Gives the asymptotic front speed $v=2\\sqrt{RD_0}$ used to move to the comoving frame.","marker":"[21]"},{"why":"Supplies the finite-difference master-equation discretization used in the numerical simulations.","marker":"[22]"},{"why":"Establishes genetic drift during range expansion as the biological process the model addresses.","marker":"[14]"},{"why":"Documents the East Asian versus European differences in human genetic loss that the paper proposes to explain.","marker":"[19]"},{"why":"Provides the human expansion context and the observed loss of genetic diversity during migration out of Africa.","marker":"[18]"}],"fun_headline_variants":["Weak disorder slows genetic loss, strong disorder speeds it up","Environmental noise: a switch for genetic loss in invasions","Spatial disorder tunes genetic drift: small slows, large accelerates","Non-monotonic disorder: genetic loss delayed then hastened"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that, in a perfectly uniform environment, the initially less abundant population is doomed to extinction; if that deterministic baseline is not accurate, the paper's claim that weak disorder delays genetic loss is measured against the wrong reference.","fun_headline_variants_meta":{"raw":{"variants":["Weak disorder slows genetic loss, strong disorder speeds it up","Environmental noise: a switch for genetic loss in invasions","Spatial disorder tunes genetic drift: small slows, large accelerates","Non-monotonic disorder: genetic loss delayed then hastened"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000623,"raw_usage":{"total_tokens":2820,"prompt_tokens":817,"completion_tokens":2003,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":433,"completion_tokens_details":{"reasoning_tokens":1934}},"tokens_in":433,"tokens_out":2003,"duration_ms":14016,"temperature":1.0,"reasoning_tokens":1934,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:42:53.870026+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test is to integrate the deterministic two-population PDE (13) with $\\xi=0$, equal fitness and no demographic noise, starting from $C_1=0.49$, $C_2=0.51$: if the minority population does not go extinct asymptotically, the paper's reference behavior for 'genetic loss' is invalid. A complementary experiment would track heterozygosity in a microbial or cell range expansion on substrates with controlled diffusivity variance, looking for the predicted non-monotonic response with a minimum of genetic loss near $\\xi\\approx0.2$.","supporting_citations":[{"cited_title":"Birzu, O","cited_arxiv_id":null,"evidence_quote":"Supplies the perturbative method for treating front dynamics in a fluctuating medium and the effective-diffusion formula that leads to Eq. (11)."},{"cited_title":"Mikhailov, L","cited_arxiv_id":null,"evidence_quote":"Provides the perturbation ansatz $C=C_0(\\zeta+\\eta(t),t)+\\delta C_1$ used to derive the Langevin equation for the front displacement."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the FKPP reaction-diffusion framework that the paper extends to a spatially random diffusion constant."},{"cited_title":"Brunet and B","cited_arxiv_id":null,"evidence_quote":"Gives the asymptotic front speed $v=2\\sqrt{RD_0}$ used to move to the comoving frame."},{"cited_title":"Sahimi, B","cited_arxiv_id":null,"evidence_quote":"Supplies the finite-difference master-equation discretization used in the numerical simulations."},{"cited_title":"Keinan, J","cited_arxiv_id":null,"evidence_quote":"Documents the East Asian versus European differences in human genetic loss that the paper proposes to explain."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the human expansion context and the observed loss of genetic diversity during migration out of Africa."}],"review_version":1}