{"id":"75d0ec74-7097-47f2-926f-4df04a722df5","arxiv_id":"1908.11144","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In the bounded-growth regime, microtubule amplification via templated severing requires a rescue-after-severing probability above a critical threshold, predicted analytically and checked by simulation.","lead":"Plant cells reorient their internal microtubule scaffolding when exposed to light. This paper builds a mathematical model showing that, when tubulin is limiting, the rescue protein must exceed a critical probability for the reorientation to spread.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The threshold prediction rests on a branching criterion applied to an offspring count that includes lagging-daughter severings as direct daughters; this can shift the predicted p+_crit independently of the FPTD approximation.","rationale":"The reader's conditional verdict is supported by the FPTD and two-crossovers approximation issues, but the most load-bearing weakness is more fundamental: the analytical threshold condition in Eq. (5) is applied to a quantity m_i that is not the mean number of direct daughters in the branching process. This matters because the central quantitative claim is that the two-crossovers theory predicts p+_crit with small relative error. If the branching criterion is misapplied, the agreement with simulations for most parameter sets could be numerically reasonable but theoretically accidental, and the prediction cannot be trusted outside the tested sets. The simulations establishing the existence of a threshold and its approximate value near 0.36 remain valuable evidence, so the paper should not be rejected outright. However, the analytical derivation needs correction or explicit justification for why counting lagging-daughter severings as direct daughters is legitimate for the supercriticality condition. A targeted simulation measurement of E[s_i] at the threshold would settle whether this concern is real. I therefore recommend keeping the conditional verdict rather than accepting the analytical prediction as presently justified.","tokens_in":11,"tokens_out":15694,"duration_ms":354087,"concrete_test":"In the full-grid simulations at the reference bounded-growth parameters and p+ = 0.36 (the simulated threshold), measure the mean number of direct severing events per microtubule lifetime, E[s_i], counting only severing events experienced by each individual microtubule body before its death. If the branching criterion M>1 is correctly applied, this quantity should be approximately 1 at threshold. If E[s_i] is substantially different from 1, then the m_i-based criterion in Eq. (8) is not the true reproduction number, and the analytical prediction of p+_crit should be re-derived using E[s_i] as the offspring mean.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"In Sec. III.B.2 the paper states that amplification occurs when M>1, where M is the number of severing events that a newly-created microtubule undergoes. The quantity actually computed in Eq. (8), however, is m_i = s_i + Σ_j c_ji, where c_ji counts severing events at pre-existing crossovers on the lagging daughter microtubules, and the text explicitly says these are treated as 'direct daughter microtubules of the mother'. In a branching process, the supercriticality condition is that the expected number of direct daughters per individual exceeds 1; adding the reproduction of lagging daughters to the mother's count inflates M. The two-crossovers correction in Eq. (11) subtracts only the fraction pcr of such crossovers that are resolved by shrinkage; it does not remove the fundamental misattribution of lagging-daughter reproduction to the mother. Thus p+_crit,(2) is not derived from a valid reproduction number, and the reported agreement with simulated thresholds for 9 of 11 parameter sets is not theoretically explained. The simulated threshold itself may be robust, but the paper's quantitative claim that the analytical theory predicts the threshold lacks a self-consistent foundation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper investigates a stochastic one-dimensional model of cortical microtubule amplification via templated severing at crossovers with a stable transverse array. The model couples dynamic instability of longitudinal microtubules with a Gamma-distributed severing waiting time and a probability p+ of rescue-after-severing of the newly created lagging plus end. Stochastic simulations show that in the unbounded-growth regime amplification occurs even at p+ = 0, while in the bounded-growth regime amplification only occurs above a critical threshold p+_crit ≈ 0.36 for the reference parameters. The authors construct an approximate analytical theory based on truncating the number of crossovers considered (one- and two-crossovers approximations) and a first-passage-time distribution, and report that the two-crossovers theory predicts the simulated threshold with relative errors below 3% for most of the eleven parameter sets in Table II.","tokens_in":17348,"tokens_out":12477,"duration_ms":123036,"significance":"If valid, the result provides a quantitative, falsifiable prediction for when rescue-after-severing is required for light-induced microtubule array reorientation, connecting the CLASP-mediated rescue effect to a sharp threshold in the bounded-growth regime. The paper combines large-scale stochastic simulations (N = 10^5 trials) with an approximate analytical closure, and the central threshold phenomenon is robustly demonstrated in the simulations. The comparison across eleven parameter sets is a useful test of the analytical approximation. However, the analytical derivation of the critical point relies on a branching-process criterion whose formulation is not self-consistent, which undermines the theoretical explanation of the observed threshold even if the simulation result stands.","major_comments":[{"comment":"The offspring count m_i is not a standard reproduction number. In Eq. (8), m_i = s_i + Σ_j c_ji, where c_ji counts the crossovers behind a severing event that may be resolved by severing on the lagging daughter microtubule. The text explicitly states that these lagging-daughter severings are treated as 'direct daughter microtubules of the mother.' In a branching process, the criticality condition is that the expected number of direct daughters per individual exceeds one; each severing event on a microtubule creates exactly one direct daughter (the lagging microtubule). Severing events on that lagging daughter are granddaughters and must be counted in the daughter's reproductive output, not the mother's. Adding them to m_i inflates M and makes the condition M > 1 (Eq. (5)) neither necessary nor sufficient for supercriticality. The two-crossovers correction in Eq. (11) subtracts only the fraction pcr of such crossovers resolved by shrinkage; it does not remove the fundamental misattribution. Therefore p+_crit,(2) is not derived from a valid reproduction number, and the agreement with simulated thresholds in Table II is not theoretically explained. The authors should either derive the threshold using the correct per-microtubule direct-daughter reproduction number or explicitly justify why their two-generation count acts as a valid criticality condition in this model.","section":"Sec. III.B.2, Eq. (8)"},{"comment":"The simulated critical values p+_crit are reported without error bars or confidence intervals. The relative errors of the two-crossovers theory range from 0% to 22.2% (row with v+=0.10, v-=0.250, rc=0.010, rr=0.020), with a second error of 15.7% (v+=0.10, v-=0.250, rc=0.015, rr=0.030). Without uncertainty estimates on the simulated thresholds, it is impossible to judge whether these deviations are statistically significant or whether the 'reasonable accuracy' claim holds uniformly across parameter sets. Please report standard errors computed from the simulation ensemble, or confidence intervals from a bootstrap analysis.","section":"Table II"},{"comment":"The FPTD approximation assumes at most one catastrophe-rescue event before reaching the target (Eq. (14)) and is shown to fail for distant targets (d = 30 and 60 μm in Fig. 4E,F). The authors argue that this failure is irrelevant because the arrival probability decays as exp(-d/l) for d ≫ l. However, this argument is not quantified for the parameter sets used in Table II, and the two-crossovers theory relies on FPTDs to both the first and second crossover. For parameter sets where the second crossover at 2d may be in a regime where the approximation degrades, the predicted p+_crit could inherit additional error beyond the branching-process issue. Please state the regime of validity of Eq. (14) for the parameters of Table II and discuss whether any of the reported predictions are affected by the FPTD failure.","section":"Sec. III.C.1, Eq. (14) and Fig. 4E-F"}],"minor_comments":[{"comment":"The abstract claims the analytical theory predicts the critical threshold with 'reasonable accuracy,' but the relative errors in Table II include values of 15.7% and 22.2%. Reporting the range of errors in the abstract and discussion would give a more accurate picture.","section":"Abstract / Discussion"},{"comment":"The definition of M as 'the number of severing events that a newly-created microtubule undergoes' is ambiguous because the preceding text and Eq. (8) also count severing events on lagging daughters. A precise definition stating that each severing event creates exactly one lagging daughter, and clarifying what is included in M, would help the reader.","section":"Sec. III.B.2, Eq. (5)"},{"comment":"The Heaviside theta condition d(2/v+ + 1/v-) in Eq. (14) is not explained. A short derivation or reference would improve readability.","section":"Eq. (14)"},{"comment":"The correlation between s_i and ⟨1-δ_{c_i,0}⟩ is stated to be negligible without numerical support. Given that this factorization is used in Eq. (36), a brief test of its validity would strengthen the derivation.","section":"Appendix B"}],"recommendation":"major_revision","confidential_remarks":"The branching-process misattribution in Eq. (8) is the most serious issue. It is not a matter of approximation accuracy but of the correctness of the criticality criterion. The simulation-based threshold claim is likely sound, but the analytical theory as presented does not provide a self-consistent derivation. I would want to see either a corrected reproduction-number calculation or a clear statement that the analytical theory is a heuristic mean-field condition whose numerical agreement with simulations is empirical, not derived. The lack of error bars in Table II also needs attention. Overall, the manuscript has interesting ideas and a novel FPTD technique, but the theoretical framework requires substantial revision before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper you asked about is worth a serious look. The main new result is that in the bounded-growth regime, templated severing amplifies microtubules only if the rescue-after-severing probability p+ exceeds a threshold, and the authors give a two-crossovers approximation that predicts this threshold with good accuracy for most parameter choices. The unbounded-regime result is a useful contrast to the CLASP-focused experiments, and the model is simple but the simulations are systematic. The analytical approximation is not curve-fitting: M+ and S are measured from independent simulations of the same model, not fitted to the target.\n\nOne concern from the stress-test was that the branching criterion M>1 is invalid because Eq. (8) counts severing events at pre-existing crossovers on lagging daughters as direct offspring of the mother. On reading the text, I think that concern misfires. The 'mother' whose offspring is being counted is the newly-created lagging microtubule itself, and the crossovers at (n-1)d, ... are on its own filament. They are legitimate sites for direct severing events that produce daughters. So the supercriticality condition is conceptually fine.\n\nThe real soft spots are more modest. The simulated critical values in Table II have no error bars; with 10^5 trials you could estimate them cheaply and should. The analytical theory uses simulation-measured M+ and S, so it is not fully first-principles, though the paper is honest about that. The two-crossovers approximation loses accuracy for two of eleven parameter sets (relative errors 16% and 22%), and the FPTD approximation is shown to fail for distant targets (Fig. 4E,F). Those are not fatal, but they mean the quantitative agreement should be read as 'good for most of the tested range', not universal.\n\nWho is this for? Cytoskeleton biophysicists and anyone modeling stochastic amplification in biological systems will get value from the threshold concept and the approximate FPTD technique. I would not desk-reject. The paper deserves a proper referee; an editor should send it out, and the referee should ask for simulation error bars, a discussion of the two failure cases, and ideally a statement of what parameters are accessible to experiment. I would engage with it.","headline":"Worth reading and citing: the bounded-growth threshold for templated severing is well supported by simulations and the two-crossovers theory lands for most parameters; the stress-test's offspring-count objection collapses on close reading.","tokens_in":17771,"tokens_out":12363,"would_cite":true,"duration_ms":116518,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["87.16.Ka","87.10.Mn"],"model":"deepseek-v4-flash","headline":"Microtubule copying flips on above a rescue threshold of ~0.36.","keywords":["microtubule amplification","templated severing","rescue-after-severing","critical threshold","dynamic instability","bounded-growth regime","unbounded-growth regime","first-passage time distribution"],"falsifier":"In a bounded-growth reconstituted assay with transverse spacing fixed at 1.5 µm, tune $p_+$ (e.g., by varying CLASP concentration) and track individual severing lineages; the two-crossovers theory predicts extinction below $p_+^*=0.361$ and exponential amplification above it. Observing any lineage that survives repeated severing at $p_+=0.25$, or an extinction-to-amplification transition outside roughly 0.34–0.38, would falsify the predicted threshold.","tokens_in":16849,"feed_emoji":"🌱","tokens_out":14336,"duration_ms":133024,"temperature":0.7,"pith_summary":"This paper asks when the severing of one microtubule at a crossover with another can multiply the severed population fast enough to reorient the cortical microtubule array of a plant cell. The authors build a one-dimensional stochastic model of a longitudinal microtubule undergoing dynamic instability—stochastically switching between growth and shrinkage—on a grid of stable transverse filaments, with a tunable probability $p_+$ that the lagging fragment is rescued immediately after severing. They find a sharp dichotomy: in the bounded-growth regime, where filaments have finite lifetime and would die out, amplification occurs only if $p_+$ exceeds a critical value (about 0.36 for the reference parameters), while in the unbounded-growth regime, where filaments on average grow forever, amplification occurs even at $p_+=0$ and $p_+$ only tunes its speed and reliability. The paper then derives an approximate analytical formula for the critical $p_+$ by replacing the full transverse grid with two filaments and counting how many severing descendants a newborn filament produces on average, showing the predicted threshold matches simulations within a few percent for most parameter sets.","feed_headline":"Microtubule copying flips on above a rescue threshold of ~0.36","feed_subtitle":"In growth-limited cells, severed filaments multiply only above a rescue-after-severing threshold that theory can predict.","key_machinery":"The machine is a one-dimensional dynamic-instability model of a longitudinal microtubule in a periodic grid of stable transverse filaments, using the classical two-state growth-shrinkage kinetics with catastrophe and rescue rates. The load-bearing piece is the offspring-size criterion $M = [p_+ + (1-p_+) R^-_d(d)] M_+ > 1$, where $M_+$ is the mean number of severing descendants of a filament born in the growing state and $R^-_d(d)$ is the probability that a shrinking filament recovers the birth length $d$; from this condition the authors obtain an explicit formula for the critical $p_+^*$. To evaluate $M_+$ they develop a two-crossovers approximation in which only the first two transverse filaments matter, supplemented by a first-passage-time distribution for reaching a nearby crossover built from legal paths with at most one catastrophe-rescue cycle. The approximation's key output is the probability $p_{cr}(p_+)$ that a crossover at the first filament is erased by shrinkage after a severing at the second, which injects the $p_+$ dependence into the threshold.","core_discovery":"The central claim is that templated severing is self-amplifying in the unbounded-growth regime but requires a threshold probability of rescue-after-severing in the bounded-growth regime. Writing the mean number of severing descendants of a newborn filament as $M = [p_+ + (1-p_+) R^-_d(d)] M_+$, where $R^-_d(d)$ is the splitting probability that a shrinking filament recovers the birth length $d$ and $M_+$ is the offspring count for a filament born growing, the authors reduce the amplification condition to $M>1$ and solve it for a closed-form critical value $p_+^*$. The evaluation of $M_+$ relies on a two-crossovers approximation in which only the first two transverse filaments are retained and at most one catastrophe-rescue event is allowed before a crossover; this yields $p_+^*=0.361$ against a simulated $0.360$ for the reference parameters, with relative errors under 3% for most of the other tested parameter sets. In the unbounded-growth regime the average lifetime is infinite, so a filament can be severed without bound and descendants accumulate even when every lagging plus end shrinks, which removes the need for any threshold.","pith_inferences":["A testable corollary not drawn in the paper: the threshold should depend mainly on the probability that a severed lagging filament reaches the next transverse filament, so a single measurement of that crossing probability in vivo could predict whether reorientation will proceed.","The one-catastrophe first-passage-time approximation is shown to fail for distant targets; this suggests the predicted threshold would drift if the transverse array is sparse (large spacing), a regime the paper does not simulate.","The model freezes the transverse array as an inert template, but real transverse filaments are themselves dynamic and consume tubulin; extending the model to let the template erode during amplification would likely raise the effective threshold over time and could explain why reorientation slows before completion."],"forward_implications":["In the bounded-growth regime, a newly severed filament and all its descendants go extinct with probability one when $p_+$ lies below the critical value, so the threshold is a true sharp phase boundary in the model.","Above the threshold, the number of longitudinal microtubules grows exponentially; the growth rate and the probability that amplification succeeds both rise monotonically with $p_+$.","In the unbounded-growth regime, amplification occurs even at $p_+=0$; rescue-after-severing only speeds up the process and lowers the extinction probability.","Since tubulin depletion slows growth, a system that begins in the unbounded-growth regime can be pushed into the bounded-growth regime, where the threshold becomes decisive for whether reorientation completes."],"supporting_citations":[{"why":"Supplies the classical dynamic-instability model and the definition of bounded vs unbounded growth that the paper's regimes rest on.","marker":"[8]"},{"why":"Provides the experimentally measured parameter set for the unbounded-growth regime and the Gamma severing-waiting-time distribution, plus the observation that rescue-after-severing speed controls amplification.","marker":"[9]"},{"why":"Establishes that katanin-mediated severing at crossovers drives cortical-array reorientation, the biological phenomenon being modelled.","marker":"[5]"},{"why":"Shows the in vivo role of rescue of the severed lagging plus end, motivating the central parameter $p_+$.","marker":"[6]"},{"why":"Provides the splitting probabilities for a microtubule in an interstitial strip used in the critical-threshold formula.","marker":"[13]"},{"why":"Supplies the lifetime density functions and ultimate survival probabilities used in the first-passage-time and crossover-resolution calculations.","marker":"[14]"},{"why":"Provides the dynamic parameters for the bounded-growth regime used in the simulations that locate the critical threshold.","marker":"[12]"},{"why":"Supports the steady-state length distribution of microtubules confined to a finite strip, used to treat the interstitial dynamics.","marker":"[11]"}],"fun_headline_variants":["Microtubule amplification needs rescue probability above threshold","Templated severing amplifies only above critical rescue threshold","Bounded growth: microtubule amplification requires critical rescue","Predicting the critical threshold for microtubule amplification","Rescue-after-severing threshold controls microtubule amplification"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analytical threshold assumes the full grid of transverse microtubules can be replaced by just two filaments, and that at most one catastrophe-rescue event occurs before a crossover is reached; the paper's own comparisons show the resulting first-passage-time approximation fails for distant targets, so the closure is only locally valid.","fun_headline_variants_meta":{"raw":{"variants":["Microtubule amplification needs rescue probability above threshold","Templated severing amplifies only above critical rescue threshold","Bounded growth: microtubule amplification requires critical rescue","Predicting the critical threshold for microtubule amplification","Rescue-after-severing threshold controls microtubule amplification"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00023,"raw_usage":{"total_tokens":1515,"prompt_tokens":1009,"completion_tokens":506,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":625,"completion_tokens_details":{"reasoning_tokens":430}},"tokens_in":625,"tokens_out":506,"duration_ms":5334,"temperature":1.0,"reasoning_tokens":430,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:23:22.158042+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"In a bounded-growth reconstituted assay with transverse spacing fixed at 1.5 µm, tune $p_+$ (e.g., by varying CLASP concentration) and track individual severing lineages; the two-crossovers theory predicts extinction below $p_+^*=0.361$ and exponential amplification above it. Observing any lineage that survives repeated severing at $p_+=0.25$, or an extinction-to-amplification transition outside roughly 0.34–0.38, would falsify the predicted threshold.","supporting_citations":[{"cited_title":"MAPs in plant cells: delineating microtubule growth dynamics and organization","cited_arxiv_id":null,"evidence_quote":"Supplies the classical dynamic-instability model and the definition of bounded vs unbounded growth that the paper's regimes rest on."},{"cited_title":"The cortical microtubule array: from dynamics to organization","cited_arxiv_id":null,"evidence_quote":"Provides the experimentally measured parameter set for the unbounded-growth regime and the Gamma severing-waiting-time distribution, plus the observation that rescue-after-severing speed controls amplification."},{"cited_title":"severing","cited_arxiv_id":null,"evidence_quote":"Establishes that katanin-mediated severing at crossovers drives cortical-array reorientation, the biological phenomenon being modelled."},{"cited_title":"Molecular Biology of the Cell","cited_arxiv_id":null,"evidence_quote":"Shows the in vivo role of rescue of the severed lagging plus end, motivating the central parameter $p_+$."},{"cited_title":"Physical aspects of the growth and regulation of microtubule structures","cited_arxiv_id":null,"evidence_quote":"Provides the splitting probabilities for a microtubule in an interstitial strip used in the critical-threshold formula."},{"cited_title":"CLASP stabilization of plus ends created by severing promotes microtubule creation and reorientation","cited_arxiv_id":null,"evidence_quote":"Supplies the lifetime density functions and ultimate survival probabilities used in the first-passage-time and crossover-resolution calculations."},{"cited_title":"Tran, R.A","cited_arxiv_id":null,"evidence_quote":"Provides the dynamic parameters for the bounded-growth regime used in the simulations that locate the critical threshold."},{"cited_title":"SPR2 protects minus ends to promote severing and reorientation of plant cortical microtubule arrays","cited_arxiv_id":null,"evidence_quote":"Supports the steady-state length distribution of microtubules confined to a finite strip, used to treat the interstitial dynamics."}],"review_version":1}