{"id":"045faf5a-3876-4d6e-9c89-9b0e04099e9f","arxiv_id":"2509.01817","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Colonial taxis based on individual cell steering produces a counterproductive down-gradient drift that does not shrink with colony size, while kinesis yields up-gradient colony drift with speed independent of colony size.","lead":"Protozoan colonies were modeled as rigid discs of cells whose flagella either steer or jitter in response to a chemical gradient. The analysis shows steering-based taxis produces counterproductive drift that does not shrink with colony size, while noise-based kinesis keeps the colony moving up-gradient at a speed independent of colony size.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Taxis-failure conclusion is model-specific: it rests on the linear, unsaturated, orientation-only steering response in Eq. (8), and the authors' own caveat that saturation or force modulation removes the down-gradient drift.","rationale":"The paper is a careful multiscale derivation backed by Monte Carlo simulations, with explicit self-consistency checks and honest caveats. The taxis down-gradient drift is a real consequence of the stated model: for perfectly symmetric flagellar placement, all steering torques cancel while the deflected force directions produce a net force opposite the gradient. The kinesis drift for symmetric colonies is also correctly derived from the multiplicative-noise modulation of flagellar orientation variance. Both results are internally consistent. The central weakness is external validity: the headline conclusion that colonial taxis fails relies on a specific, unsaturated, orientation-only steering response, and the authors themselves identify the exact variations that would remove the counterproductive drift. Because the parameters k_T and k_K are fit to the same experimental dataset used for qualitative support, the model-to-organism step is not independently constrained. The reader's CONDITIONAL verdict already captures this limitation; our stress-test confirms it without finding an internal mathematical error. No single assumption appears fatally wrong within the stated scope, but the biological generalizability of the taxis-failure claim is contingent on the response structure. Hence the verdict should remain CONDITIONAL, with no change from the reader's assessment.","tokens_in":40421,"tokens_out":16733,"duration_ms":167958,"concrete_test":"Test the model-dependence of the taxis failure by Monte Carlo simulation of the taxis model (8) with a saturating steering response: replace the linear bias sin(Θ(c)+α_j−θ_g) by sinh(k_T g sin(...))/k_T g, or equivalently cap the mean flagellar deflection at a few σΘ, and repeat the Figure 6 protocol for N=10 over m_T ∈ {0.5, 1, 2, 3}. If the mean drift projection along ˆe_g no longer turns negative or plateaus at a positive value, then the counterproductive term is an artifact of the unsaturated linear-response assumption. If negative drift persists under saturation, the central taxis-failure claim is robust to this model variation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The biological message — steering-based taxis fails in colonies while kinesis succeeds — is load-bearing for the paper. The taxis half depends entirely on Eq. (8): each flagellum exerts a constant-magnitude force whose orientation is biased linearly and without saturation. The counterproductive down-gradient drift in Eq. (16) is a first-order-in-k_T g effect that survives only because the response reorients flagella while leaving force magnitude fixed and does not saturate. The authors explicitly concede in Section 3.2 that this drift “would drop out under various model variations, such as having a quick saturation to the steering response... or modulating the flagellar force rather than orientation.” Thus the paper does not establish a general failure of colonial taxis; it establishes a failure of one particular response architecture. Whether that architecture describes choanoflagellates is not independently supported: k_T and k_K are fitted to the same S. rosetta measurements used for qualitative comparison (Section 6 preamble), and no direct evidence rules out force modulation or saturation. The kinesis half, while internally consistent, similarly inherits the 2D geometry: the claimed independence of colony size in Eq. (19) follows from γ_t ∝ N with a ∝ N in the disc model; for a 3D spherical colony (N ∝ a^2, γ_t ∝ a) the shape-independent term would scale as N^{1/2}, changing the quantitative claim. The abstract presents the N-independence as a general result despite the 2D restriction.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a two-dimensional stochastic model of a rigid circular colony of flagellated cells, where each flagellum exerts a constant-magnitude force with orientation governed by either a taxis drift-diffusion process (Eq. 8) or a kinesis noise-modulation process (Eq. 9). Through nondimensionalization and homogenization with formal small parameter ε = σΘ ζ/N^2, the authors derive effective colony-level SDEs and long-time drift formulas. The taxis analysis yields an up-gradient drift that decays with colony size plus a counterproductive down-gradient drift; for perfectly symmetric flagellar placement the colony drifts down-gradient. The kinesis analysis yields an up-gradient drift whose leading term is independent of colony size in the two-dimensional disc geometry. Monte Carlo simulations are compared to the asymptotic formulas, with good agreement for moderate parameters.","tokens_in":40723,"tokens_out":7514,"duration_ms":86028,"significance":"If correct, the paper gives a concrete mechanical explanation for why colonial protozoa such as S. rosetta may rely on kinesis rather than steering-based taxis: the steering response of individual cells can produce a colony-level drift in the wrong direction, whereas noise modulation produces a robust up-gradient drift. The derivation is a strength: the effective drift formulas are obtained from the stated stochastic models rather than fitted, with explicit self-consistency conditions (42) and (53) and Monte Carlo validation. The limitations are equally clear: the taxis conclusion depends on the specific linear, unsaturated, orientation-only response in Eq. (8), and the kinesis size-independence is a consequence of the two-dimensional disc scaling. These scope restrictions must be reflected in the abstract and conclusions.","major_comments":[{"comment":"The counterproductive down-gradient drift for χ=0 (or as the first term in the expansion) is an artifact of the linear, unsaturated, orientation-only steering response in Eq. (8). The authors explicitly state that this drift \"would drop out under various model variations, such as having a quick saturation to the steering response... or modulating the flagellar force rather than orientation.\" Because Eq. (8) is not independently constrained by data, the paper cannot claim a general failure of colonial taxis; it establishes a failure of one response architecture. Please either analyze robustness (e.g., saturating response or force modulation) or qualify the abstract/conclusions to \"within the linear orientation-steering model.\"","section":"Section 3.2, Eq. (16)"},{"comment":"The claimed colony-size independence of kinesis drift follows from the two-dimensional disc scaling a ∝ N and γ_t ∝ a. For a three-dimensional spherical colony with cells on the surface, N ∝ a^2 while γ_t ∝ a, so the factor N/γ_t in Eq. (19) scales as a ∝ N^{1/2}, not as N^0. The discussion in Section 7 only says \"general aspects\" carry over and does not retract the abstract's unconditional statement. Please derive or at least state the 3D scaling and qualify all N-independence claims as 2D disc-geometry results.","section":"Section 3.3, Eq. (19) and Section 7"},{"comment":"For the large σΘ^2 values used in the kinesis Monte Carlo comparison, the formal self-consistency condition N ≫ (ζδ^2 + σΘ^2 ζ^2)^{1/3} is violated at N=7, and the asymptotic theory overpredicts the simulation drift (as the authors note). This does not invalidate the physical parameter regime (σΘ^2 ≈ 0.002), but it overstates the claim of \"quantitative agreement.\" Please add simulations at physically relevant σΘ^2 and/or a higher-order correction, and clearly separate validated from exploratory parameter regimes.","section":"Section 6.2, Eqs. (42)/(53), Fig. 9"}],"minor_comments":[{"comment":"The response coefficients k_T and k_K are taken from fits to the same S. rosetta colony data used for qualitative comparison; the later claim of consistency with Kirkegaard et al. is therefore not an independent test. Please clarify that this is a parameter-matching exercise, not validation.","section":"Section 6 preamble"},{"comment":"The second component uses Θ(c)(\\tilde t) with a tilde, likely a typo for Θ(c)(t); please correct.","section":"Eq. (4b)"},{"comment":"The horizontal axis label appears truncated: the explicit variable σΘ^2 is missing. This makes the figure hard to read without referring to the text.","section":"Figure 9"},{"comment":"The abstract states the kinesis drift is \"independent of colony size\" without the two-dimensional qualifier that the model and Section 7 rely on. A brief qualifier such as \"in the two-dimensional disc model\" would align the abstract with the actual results.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"This is a solid mathematical modeling contribution with clear derivations and Monte Carlo checks. The main risk is overgeneralization: the taxis-failure conclusion is tied to a particular steering architecture, and the kinesis size-independence is a 2D scaling result. These are fixable within the manuscript's scope by tightening the claims and adding robustness/scaling discussion; I would not reject."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should read this for the homogenization result. The paper takes cell-level stochastic models of flagellar steering (taxis) and noise modulation (kinesis), coarse-grains them to colony-level dynamics, and derives closed-form drift formulas. The counterproductive down-gradient drift under steering taxis, which survives for large symmetric colonies, and the kinesis drift that does not shrink with colony size in the 2D disc model are new and non-obvious. The derivation is careful: explicit nondimensionalization, a defined small parameter, self-consistency conditions, and Monte Carlo validation across colony sizes and response strengths. The taxis predictions track simulations including the drift reversal at larger mT; the kinesis theory works where its validity condition holds. The authors also state their main caveats plainly, which I respect.\n\nThe soft spots are real but mostly scoped by the authors themselves. The taxis-failure conclusion is tied to Eq. (8): reorientation of a constant-magnitude flagellar force with linear, unsaturated response. In Section 3.2 they concede that the down-gradient drift would drop out under quick saturation or force modulation. Since they don't independently show choanoflagellates use this architecture — the response parameters are fitted to the same S. rosetta data used for qualitative support — the biology is plausible but not established. That's a limit, not a flaw in the math. The kinesis N-independence is a 2D result; a ∝ N in the disc gives N/γt constant, while a 3D colony would scale differently. They note the 2D restriction in the introduction and discussion but the abstract drops it. Minor. The self-consistency condition (53) is violated for the larger σΘ explored, and the theory overpredicts at N=7; they acknowledge this too. No code or error bars for the simulations — minor given the comparisons are convincing.\n\nShould this go to peer review? Yes. It is a serious, technically detailed paper with genuine new results. A referee should push for clearer scoping of the biological claim and explicit acknowledgment of the 2D/3D issue, but the core derivation deserves a venue. I'd bring it to a reading group on active matter or multiscale stochastic dynamics. Worth citing for the homogenization formulas.","headline":"Careful multiscale asymptotic derivation of colony-level taxis/kinesis drift; new and honest about its limits, but the headline biological claim is architecture-specific.","tokens_in":41268,"tokens_out":2691,"would_cite":true,"duration_ms":31960,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["92C17","60H10","92C05"],"pacs":[],"model":"deepseek-v4-flash","headline":"Steering-based taxis fails in protozoan colonies: up-gradient drift shrinks with colony size while a size-independent down-gradient drift dominates; noise-based kinesis keeps working at any size.","keywords":["colonial protozoa","taxis","kinesis","stochastic differential equations","asymptotic homogenization","flagellar asymmetry","choanoflagellates","chemotaxis"],"falsifier":"Track colonies of S. rosetta of different cell counts in a steady oxygen gradient and measure net drift: the taxis prediction is up-gradient speed that decreases with N and turns negative for large or unusually symmetric colonies, while kinesis predicts positive, approximately size-independent drift. A clean sign reversal with colony size would support the taxis-failure mechanism; a more direct test varies flagellar placement disorder u and checks whether taxis effectiveness scales as (u/l)²/N as Eq. (17) predicts.","tokens_in":40196,"feed_emoji":"🦠","tokens_out":9723,"duration_ms":98566,"temperature":0.7,"pith_summary":"The paper asks whether a protozoan colony can still climb an environmental gradient when its cells respond independently, each sensing the local gradient and adjusting only its own flagellum. Building a stochastic model of a rigid circular colony of N cells, the authors show that steering-based taxis degrades with colony size: the useful up-gradient drift shrinks as the flagellar-placement asymmetry becomes relatively smaller, while a counterproductive down-gradient drift stays roughly constant, so large or symmetric colonies swim the wrong way. They further show that a kinesis response — modulating the noisiness of flagellar motion rather than steering — produces up-gradient drift at a speed independent of colony size, even with perfectly symmetric flagella. If correct, this explains why colonial choanoflagellates such as S. rosetta rely on kinesis rather than steering-based taxis, and suggests that noise-based responses are the ones that survive the transition from single cells to colonies.","feed_headline":"Steering-based taxis fails in protozoan colonies","feed_subtitle":"Flagellar steering drifts down-gradient at colony scale; noise-modulated kinesis keeps working at any colony size.","key_machinery":"Central object: a rigid circular colony of N cells, each exerting a flagellar force of fixed magnitude at a stochastic angle to the local surface normal, with flagellar attachment points randomly displaced and force orientations fluctuating in time. The argument rests on a time-scale separation: fast flagellar orientation dynamics (rate γΘ) are averaged out to give slow colony rotation and translation, with ε = σΘζ/N² the small parameter. Results are organized by two asymmetry measures χ and χ₂ — the vector sums of flagellar positions and their second harmonics — since uniform placement cancels all gradient response; drift scales with χ/√N (taxis) or contains a shape-independent part (kinesi","core_discovery":"Under the taxis model (Eq. 8), where each cell only reorients a constant-magnitude flagellar force to steer itself up the gradient, colony drift has a constructive part proportional to flagellar-placement asymmetry χ that shrinks with colony size N and a counterproductive down-gradient part that does not; at χ = 0 the colony drifts purely down-gradient (Eq. 16). Under the kinesis model (Eq. 9), where rotational noise grows when a flagellum faces away from the gradient, the colony drifts up-gradient at a speed independent of N even with symmetric flagella, because noisier wrongly-oriented flagella lose more directed force (Eq. 19). Both effective colony response coefficients scale as χ N^{-1/","pith_inferences":["If real choanoflagellate flagella modulate force magnitude rather than only orientation — or saturate their steering response quickly — the paper's own analysis implies the down-gradient taxis drift disappears; measuring beat asymmetry versus force magnitude in a gradient would discriminate the two response classes.","The size-independence of kinesis drift may be a general principle for the transition to multicellularity: a noise-modulated response keeps functioning as cells are added and geometry becomes disordered, whereas a steering response demands sustained symmetry, so kinesis should be favored in lineages where colonies grow by adding independently beating cells.","The formula predicting taxis effectiveness scaling as (u/l)²/N suggests a quantitative experiment: colonies with deliberately controlled flagellar placement disorder u should show net up-gradient drift only when N(u/l)² exceeds a threshold.","The 2D model likely underestimates the minimal colony size needed for the asymptotic analysis in 3D (colony radius grows as √N rather than N), so the taxis-failure and kinesis-success predictions should be most robust for larger 3D colonies."],"forward_implications":["Colonies that steer by reorienting constant-strength flagellar forces lose the ability to navigate up gradients as cell number grows; at large N the down-gradient drift dominates and the colony moves toward lower stimulus.","Perfectly symmetric flagellar placement makes steering-based taxis actively harmful: with torques cancelled, all the steering deflections rotate the flagellar forces to push down the gradient.","Kinesis delivers up-gradient motion independent of colony size, so a colony can keep responding to gradients without coordination or geometric regularity, consistent with the observed dominance of kinesis in S. rosetta.","The asymmetry-dependent part of the kinesis drift decays roughly as N^{-3} for the parameter range studied, so small colonies get a geometric boost but large ones rely on the shape-independent mechanism.","The derived drift formulas give a direct, parameterized route to Keller-Segel-type continuum models for suspensions of colonies."],"supporting_citations":[{"why":"Supplies the phenomenological taxis and kinesis models the paper adapts to flagellar force orientation, plus the S. rosetta aerotaxis observations showing kinesis outperforms taxis and the fitted response strengths used in simulations.","marker":"Kirkegaard et al. [2016a]"},{"why":"Provides the measured statistics of flagellar beating (correlation time, orientation variance, uncoordinated beating) that set the parameter values and motivate the stochastic flagellar model.","marker":"Kirkegaard et al. [2016b]"},{"why":"Establishes the theoretical approach of quantifying colonial swimming with independent flagella and supplies flagellar force magnitude estimates.","marker":"Roper et al. [2013]"},{"why":"The companion paper's colony dynamics (drag coefficients, flagellar force kinematics, asymmetry measures) that this work augments with taxis and kinesis responses.","marker":"Ashenafi and Kramer [2024]"},{"why":"Defines the chemotactic index used to measure navigation success and provides the emergent-chemotaxis cluster model whose force-cancellation logic the paper contrasts with its own.","marker":"Fancher et al. [2017]"},{"why":"A similar autonomous Langevin model for cell-level taxis responses whose structure the paper's per-cell steering model parallels.","marker":"Hopkins and Camley [2019]"},{"why":"Central limit theorem for fast-slow stochastic systems used to derive the effective noise on colony orientation in the taxis reduction.","marker":"Bouchet et al. [2016]"},{"why":"Averaging and homogenization theorems applied to coarse-grain the kinesis flagellar dynamics into colony-scale equations.","marker":"Stuart and Pavliotis [2008]"}],"fun_headline_variants":["Taxis drift reverses in colonies; kinesis size-proof","Steering-based taxis fails in colonies; kinesis works","Colony taxis down-gradient, kinesis up-gradient","Flagellar steering fails at scale; noise kinesis prevails","Taxis fails for colonies, kinesis independent of size"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The central conclusions follow from the assumption that each cell responds to the gradient only by reorienting a flagellar force of fixed strength, with flagellar dynamics much faster than colony rotation; the authors note the down-gradient taxis drift would disappear if the steering response saturated quickly or modulated force magnitude, and the time-scale separation is only marginally satisfied at the largest forces and fluctuation levels tested.","fun_headline_variants_meta":{"raw":{"variants":["Taxis drift reverses in colonies; kinesis size-proof","Steering-based taxis fails in colonies; kinesis works","Colony taxis down-gradient, kinesis up-gradient","Flagellar steering fails at scale; noise kinesis prevails","Taxis fails for colonies, kinesis independent of size"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000298,"raw_usage":{"total_tokens":1573,"prompt_tokens":767,"completion_tokens":806,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":511,"completion_tokens_details":{"reasoning_tokens":719}},"tokens_in":511,"tokens_out":806,"duration_ms":7863,"temperature":1.0,"reasoning_tokens":719,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T12:09:37.347567+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Track colonies of S. rosetta of different cell counts in a steady oxygen gradient and measure net drift: the taxis prediction is up-gradient speed that decreases with N and turns negative for large or unusually symmetric colonies, while kinesis predicts positive, approximately size-independent drift. A clean sign reversal with colony size would support the taxis-failure mechanism; a more direct test varies flagellar placement disorder u and checks whether taxis effectiveness scales as (u/l)²/N as Eq. (17) predicts.","supporting_citations":[{"cited_title":"Dayel, Rachel E","cited_arxiv_id":null,"evidence_quote":"Establishes the theoretical approach of quantifying colonial swimming with independent flagella and supplies flagellar force magnitude estimates."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The companion paper's colony dynamics (drag coefficients, flagellar force kinematics, asymmetry measures) that this work augments with taxis and kinesis responses."},{"cited_title":"Emergent versus Individual - Based Multicellular Chemotaxis","cited_arxiv_id":null,"evidence_quote":"Defines the chemotactic index used to measure navigation success and provides the emergent-chemotaxis cluster model whose force-cancellation logic the paper contrasts with its own."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"A similar autonomous Langevin model for cell-level taxis responses whose structure the paper's per-cell steering model parallels."},{"cited_title":"Large Deviations in Fast -- Slow Systems","cited_arxiv_id":null,"evidence_quote":"Central limit theorem for fast-slow stochastic systems used to derive the effective noise on colony orientation in the taxis reduction."}],"review_version":1}