REVIEW 3 major objections 5 minor 31 references
Chemotaxis of branched cells in complex environments
T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read A model of branched-cell chemotaxis predicts a speed–accuracy tradeoff and places neutrophils in the fast, less-accurate regime.
desk verdict A solid, honest modeling paper whose central biological claim—neutrophils are fast cells—rests on the same FMI data used to calibrate the model; the speed-accuracy tradeoff itself is a genuine and reproducible model prediction. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is a branched cell on a hexagonal network, described by coupled variables for each arm: arm length $x_i$, fraction of engaged slip-bond adhesions $n_i$, and local actin treadmilling velocity $v_i$. The chemoattractant enters by increasing the maximal polymerization speed $\alpha$ at each arm tip according to $\alpha(r)=\alpha_0[1+\epsilon(1+C/c_0)c(r)/(c(r)+C)]$, so a gradient biases the competition between arms without any explicit gradient-sensing module. The parameter $\alpha_0$ is the free control knob: below a critical value cells cannot polarize and migrate, while above a higher critical value they get trapped in slow-mode decision events. This machinery produces the tradeoff because faster dynamics make the cell more likely to lose polarity and choose arms based on noise rather than on the weak signal bias.
What would settle it
In the same zebrafish wound assay, measure per-cell speed and wrong-turn rate before and after the chemokine appears: the model predicts a clear negative correlation, with faster cells making more wrong turns, and predicts that drug-treated cells with reduced actin activity become slower but more accurate; observing instead that the fastest cells are also the most accurate would falsify the tradeoff.
Extended reading notes
Core claim
The paper claims that the chemotactic performance of a branched cell is set by the same directional decision-making machinery that governs its random migration, with the chemokine entering only through a local enhancement of actin-polymerization activity at each arm tip. On a single junction and on a hexagonal network, raising the baseline actin activity $\alpha_0$ makes cells move faster but also makes them choose the wrong arm more often, and at high $\alpha_0$ they can enter slow-mode states with two competing elongated arms. Comparing the model with neutrophils migrating to a laser wound in zebrafish and with cells migrating between pillars, the paper finds that the experimental cells sit in the high-$\alpha_0$ regime. The authors conclude that wild-type neutrophils behave as fast cells that compromise chemotaxis accuracy for arrival speed, a strategy that works because some cells are always close to the wound.
Load-bearing premise
The load-bearing premise is that a chemoattractant gradient influences a branched cell only by locally increasing actin-polymerization activity at each arm tip, with no receptor-binding kinetics, adaptation, or cell-wide gradient comparison, and that the wild-type activity $\alpha_0$ is set by fitting the measured forward-migration-index increase; if real neutrophils sense gradients through additional mechanisms, the prediction that they are fast cells may not transfer.
Editorial extensions
If this is right
- Wild-type neutrophils are predicted to prioritize arrival time over directional accuracy, so their recruitment efficiency depends on being densely distributed rather than on each cell navigating perfectly.
- Near the wound, where the chemokine signal is strong, high actin activity is the optimal strategy because fast arrival outweighs occasional wrong turns.
- Far from the wound, the same fast cells meander and sometimes enter slow-mode traps, so effective long-range recruitment requires signal amplification, for example secondary chemokine secretion by other neutrophils.
- Drug-treated cells with reduced actin polymerization or myosin-II activity are predicted to become slower but more accurate, matching the qualitative changes in forward migration index reported for CK666- and blebbistatin-treated cells.
- On larger networks, high-activity cells are increasingly slowed by slow-mode events, so the benefit of speed may reverse as tissue geometry becomes more complex.
Reading between the lines
- A testable extension is to look for the same speed–accuracy tradeoff within a single, genetically identical cell population: individual faster cells should have lower directional persistence and more wrong turns than slower cells.
- Because the model treats chemokine sensing only as local actin enhancement, real receptor adaptation or cell-wide gradient comparison could weaken the predicted tradeoff; experiments with cells whose adaptation machinery is perturbed would test this directly.
- The mechanism suggests a general design principle for cells moving through crowded geometries: speed and directional accuracy are coupled through decision sampling time, so any intervention that slows a cell may improve its navigation accuracy in complex environments.
- On pillar lattices with different spacings, the model appears to predict that denser or more obstructive geometries will increase slow-mode trapping of fast cells, which could be measured as a drop in arrival success at high pillar density.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript extends a previously developed coarse-grained model of branched-cell migration on hexagonal networks to include chemokine sources. The chemokine acts by locally upregulating the actin polymerization parameter α at arm tips (Eqs. 1 and 3). The model reports a speed-accuracy tradeoff at single junctions and on networks, with an optimal intermediate activity for weak signals. The authors compare the model to in vivo zebrafish neutrophils and in vitro PLB-985 cells, and conclude that neutrophils operate in the high-α0 (fast) regime, sacrificing accuracy for speed.
Significance. If valid, the paper offers a mechanistic explanation for why neutrophils are fast but inaccurate chemotaxers, linking intracellular actin dynamics to tissue-scale recruitment. The speed-accuracy tradeoff is a nontrivial emergent prediction of the model, not imposed by construction. The model reproduces several qualitative features of neutrophil behavior, including the transient slowdown and reorientation after wounding and the occurrence of slow-mode events. The quantitative support for the fast-cell inference, however, rests on a circular calibration and lacks uncertainty quantification.
major comments (3)
- [Comparison of the model with the chemotaxis of neutrophils (Fig. 4E-F; SI S-6, S-7)] The conclusion that wild-type neutrophils behave as fast cells is a restatement of the calibration rather than an independent prediction. The paper states that the experimental FMI increase is used 'to calibrate the value of the parameter ´0' and that this increase 'is captured by the model only when we are in the regime of large ´0' (Fig. 4E-F); the simulations in Fig. 4 use ´0=12. The later conclusion that 'the neutrophils correspond to the high-activity cells of our model' is therefore supported by the same observable that fixed the parameter. To break the circularity, the authors should either constrain ´0 using an independent observable (e.g., absolute cell speed, meandering-path frequency, or slow-mode statistics) and then predict the FMI increase out-of-sample, or provide an identifiability analysis showing that the FMI increase uniquely selects the high-´0 regime over the plausible parameter range.
- [Comparison of the model with the chemotaxis of neutrophils (Fig. 4C-F; SI S-6, S-7)] The quantitative comparison to experiment is underdetermined. The calibration targets a single summary statistic (the relative FMI increase after BM time), while the model has several free parameters (´0, ϵ, C/c0, y0/yend, d, Ã); no simulation error bars, confidence intervals, or parameter sensitivity analysis are reported. This also applies to the central tradeoff curves in Figs. 2B and 3C, where the non-monotonic features (e.g., the minimum of Tarr in Fig. 3C(i) and the sharp rise in P(wrong) near ´0≈8 in Fig. 2B(i)) could be within noise. Furthermore, the slow-mode comparison in Fig. 5C uses a different parameter set (´0=12.6, d=3.7, Ã=0.6) from the calibration (´0=12, d=3, Ã=0.5), so it is unclear whether the qualitative agreement is robust. Please add error bars and a sensitivity analysis, and demonstrate that the fast-cell interpretation is not an artifact of the chosen parameter point.
- [Eqs. (1) and (3); SI S-3] The transfer of the speed-accuracy tradeoff to real neutrophils depends critically on the assumption that chemotaxis is mediated solely by local upregulation of α at each arm tip, with no receptor binding kinetics, adaptation, or cell-wide gradient comparison. While this coarse-graining is a reasonable modeling choice, the paper's claim that the model 'explains the functionality of these immune cells' is stronger than the evidence. I recommend tempering the conclusion or testing whether a different sensing rule (e.g., adding adaptation or a front-back comparison) preserves the predicted tradeoff; such a test would substantially strengthen the biological relevance.
minor comments (5)
- [Fig. 4C] The caption contains a typo: 'simlation' should be 'simulation'.
- [Fig. 5 caption] The panel labeled '(F)' should likely be '(iii)' for consistency with the (i)-(iii) panels in (C).
- [Bullet list on p. 6] In the bullet list, 'Fig. 3D(i)' should probably be 'Fig. 3C(i)' because the arrival-time minimum is shown in Fig. 3C(i), not in the trajectory panel 3D(i).
- [SI S-6] The sentence 'The effect of the chemokine on the FMI of the cells is given in Fig. 5E' should refer to Fig. 4E, since Fig. 4E displays the experimental FMI data.
- [Comparison of the model with the chemotaxis of neutrophils] The mapping of simulation time to experimental time ('choosing a simulation time interval of ∼ 2.5') is stated without derivation; please provide the calculation or a reference for this conversion.
Circularity Check
The speed-accuracy tradeoff is an emergent simulation result, but the conclusion that neutrophils are fast cells is largely a restatement of the calibration of alpha0 to the experimental FMI increase.
-
fitted input called prediction
[Main text, 'Comparison of the model with the chemotaxis of neutrophils', Fig. 4E-F; SI S-6; Conclusion]
"In Fig. 4E,F, we compare the effects of the chemokine presence on the FMI of the neutrophils [4] with our simulation model. We use this experimental data to calibrate the value of the parameter ´0 in our simulations. The FMI of the wild-type cells is observed to significantly increase after the BM time (by ∼ 20% − 30%, Fig. 4E). This behavior is captured by the model only when we are in the regime of large ´0."
The paper's central biological inference, that wild-type neutrophils correspond to the fast, high-alpha0 regime, is not an out-of-sample prediction. The same experimental FMI increase after BM is used both to select alpha0 and as evidence that only large alpha0 reproduces it. The later statement that neutrophils 'behave as fast cells' therefore restates the calibration target. The drug-treated comparison in S-7 likewise fits reduced alpha0 (and k) to the same FMI/length/speed data and then cites the agreement as confirmation.
-
fitted input called prediction
[SI S-7, Fig. S-7C-F]
"In Fig. S-7C,D, we compare the experimental [6] and simulated changes to the FMI for the WT (DMSO) and drug-treated cells. The relative changes are qualitatively captures by the model, especially the increase in FMI due to the chemokine for the WT cells, and the vanishing of this effect upon drug application. This, again, points to the WT cells residing in the high-´0 regime of our model, which is the only regime where we see that the FMI is significantly increased by the chemokine gradient."
The drug parameters (alpha0 and k) are chosen to fit the experimental cell length, speed, and FMI data from the drug-treated cells, and then the qualitative agreement is used as evidence that wild-type cells sit in the high-alpha0 regime. Since the model's FMI increase appears only in that regime, the agreement is a byproduct of the fit rather than an independent test.
full rationale
The speed-accuracy tradeoff itself is an emergent output of the stochastic branched-cell equations and is not written into the model: the single-junction and network simulations show Tesc, P(wrong), Lpath, and FMI varying with alpha0 in a nontrivial way. That part of the paper is self-contained and not circular. The circularity concerns the application to neutrophils. The abstract's closing claim, that neutrophils behave as fast cells compromising chemotaxis accuracy, is reached by calibrating alpha0 to the experimental FMI increase after the BM time and then asserting that this behavior occurs only at large alpha0. The FMI increase is thus both the fitting target and the evidence for the high-alpha0 assignment. The supporting observations (meandering trajectories, slow-mode events) are qualitative and are not quantitatively tied to alpha0 with error bars or parameter-identifiability analysis, so they do not break the circularity. Overall, the model's internal tradeoff is independent, but the headline experimental conclusion is partially forced by the calibration, giving a partial circularity score of 6.
Assumptions & free parameters
free parameters (7)
- alpha0 (baseline actin polymerization activity) =
12.0 for WT comparison; 10.0 for CK-666; 10.0 for blebbistatin; scanned 5.5-10.5
- epsilon (max chemokine enhancement) =
0.2 in parametric study; 0.1 in experimental comparisons
- C/c0 (chemokine saturation ratio) =
0.01 (strong) and 1 (weak)
- y0 (exponential decay length) / yend (linear cutoff) =
y0 = 1.5 d; ysource = -yend = 8 d
- d (hexagonal edge length) =
3, 3.7, 7.5, 9 in different runs
- noise amplitude Gamma =
0.1, 0.5, 0.6, 1.0, 1.2
- adopted model parameters (Table S-1) =
c=3.85, D=3.85, k=0.8, fs=5, r=5, eta=20, beta=250
assumptions (5)
- domain assumption A branched cell can be represented as N arms on a hexagonal network with actin flow, slip-bond adhesion, and an advected polarity inhibitor described by Eqs. S-1 to S-5.
- domain assumption The chemokine signal modulates cell behavior solely by enhancing local actin polymerization speed alpha at each arm tip via Eq. 3, with no receptor dynamics or internal gradient-sensing machinery.
- domain assumption The in-vivo tissue can be approximated by a regular hexagonal lattice with a linear or exponential chemokine profile.
- standard math The numerical integration of the stochastic ODEs faithfully samples model behavior over Tmax = 1000.
- ad hoc to paper The mapping of simulation units to experiment (one cell length in 3 min corresponds to simulation time ~2.5) preserves the dynamics being compared.
Cite this review
Pith. "Pith review of Chemotaxis of branched cells in complex environments." pith.science (2026). https://pith.science/paper/AIX2JVJI
@misc{pith2026250521949,
author = {Pith},
title = {Pith review of: Chemotaxis of branched cells in complex environments},
year = {2026},
howpublished = {\url{https://pith.science/paper/AIX2JVJI}},
note = {Machine review of arXiv:2505.21949}
}
read the original abstract
Cell migration in vivo is often guided by chemical signals. Such chemotaxis, such as performed by immune cells migrating to a wound site, is complicated by the complex geometry inside living tissues. In this study, we extend our theoretical model of branched-cell migration on a network by introducing chemokine sources to explore the cellular response. The model predicts a speed-accuracy tradeoff, whereby slow cells are significantly more accurate and able to follow efficiently a weak chemoattractant signal. We then compare the model's predictions with experimental observations of neutrophils migrating to the site of laser-inflicted wound in a zebrafish larva fin, and migrating in-vitro inside a regular lattice of pillars. We find that the model captures the details of the sub-cellular response to the chemokine gradient, as well as the large-scale migration response. This comparison suggests that the neutrophils behave as fast cells, compromising their chemotaxis accuracy, which explains the functionality of these immune cells.
Figures
Reference graph
Works this paper leans on
-
[1]
D. L. Bodor, W. P¨ onisch, R. G. Endres, and E. K. Paluch, Develo pmental cell 52, 550 (2020)
work page 2020
-
[2]
P. J. S´ aez, L. Barbier, R. Attia, H.-R. Thiam, M. Piel, and P . Vargas, Cell Migration: Methods and Protocols pp. 361–373 (2018)
work page 2018
- [3]
- [5]
-
[6]
J. M. Garcia-Arcos, R. Chabrier, M. Deygas, G. Nader, L. Barbier, P . J. S´ aez, A. Mathur, P. Vargas, and M. Piel, Journal of cell science 132, jcs225565 (2019)
work page 2019
- [7]
-
[8]
K. M. Yamada and M. Sixt, Nature Reviews molecular cell bio logy 20, 738 (2019)
work page 2019
-
[10]
J. E. Ron, M. Crestani, J. M. Kux, J. Liu, N. Al-Dam, P. Monzo , N. C. Gauthier, P. J. S´ aez, and N. S. Gov, Nature Physics pp. 1–11 (2024)
work page 2024
Show all 31 references
-
[12]
Poplimont, A
H. Poplimont, A. Georgantzoglou, M. Boulch, H. A. Walker, C. Coombs, F. Papaleonidopoulou, and M. Sarris, Current Biology 30, 2761 (2020)
2020
-
[13]
J. E. Ron, P. Monzo, N. C. Gauthier, R. Voituriez, and N. S. G ov, Physical Review Research 2, 033237 (2020)
2020
-
[14]
R. H. Insall and L. M. Machesky, Developmental cell 17, 310 (2009)
2009
-
[15]
Ivaska, Nature cell biology 14, 981 (2012)
J. Ivaska, Nature cell biology 14, 981 (2012)
2012
-
[16]
X.-D. Ren, W. B. Kiosses, D. J. Sieg, C. A. Otey, D. D. Schlaepfe r, and M. A. Schwartz, Journal of cell science 113, 3673 (2000)
2000
-
[17]
M. E. Berginski, E. A. Vitriol, K. M. Hahn, and S. M. Gomez, P loS one 6, e22025 (2011)
2011
-
[18]
S. L. Gupton and C. M. Waterman-Storer, Cell 125, 1361 (2006)
2006
-
[19]
P. V. Afonso, M. Janka-Junttila, Y. J. Lee, C. P. McCann, C . M. Oliver, K. A. Aamer, W. Losert, M. T. Cicerone, and C. A. Parent, Developmental cell 22, 1079 (2012)
2012
-
[20]
Re´ ategui, F
E. Re´ ategui, F. Jalali, A. H. Khankhel, E. Wong, H. Cho, J. Lee, C. N. Serhan, J. Dalli, H. Elliott, and D. Irimia, Nature biomedical engineering 1, 0094 (2017)
2017
-
[21]
Z. Song, S. Bhattacharya, R. A. Clemens, and M. C. Dinauer, I science 26 (2023)
2023
-
[22]
Strickland, D
E. Strickland, D. Pan, C. Godfrey, J. S. Kim, A. Hopke, W. Ji, M . Degrange, B. Villavicencio, M. K. Mansour, C. S. Zerbe, et al., Developmental Cell 59, 2659 (2024). Supplementary Information - Chemotaxis of branched cells in complex envi ronments Jiayi Liu 1,2, Jonathan E. R...
2024
-
[24]
the local actin treadmilling flow velocity vi at the leading edge. The dynamic equations for these variables are given by [8] ˙xi = 1 Γ i [vi − k (L − 1)] (S-1) ˙ni = r (1 − ni) − ni exp [ −vi + k (L − 1) fs ni ] (S-2) 4 ˙vi = −¶ (vi − v∗ i ) + ÃÀt (S-3) In Eq.S-1, Γ i is a non...
-
[25]
Westerfield, The Zebrafish Book; A guide for the laboratory use of zebrafish (Danio rerio) (2007)
M. Westerfield, The Zebrafish Book; A guide for the laboratory use of zebrafish (Danio rerio) (2007)
2007
-
[26]
Williantarra, A
I. Williantarra, A. Georgantzoglou, and M. Sarris, Bio-protoc ol 14, e4997 (2024)
2024
-
[27]
Poplimont, A
H. Poplimont, A. Georgantzoglou, M. Boulch, H. A. Walker, C . Coombs, F. Papaleonidopoulou, and M. Sarris, Current Biology 30, 2761 (2020)
2020
-
[28]
T.-W. Chen, T. J. Wardill, Y. Sun, S. R. Pulver, S. L. Renning er, A. Baohan, E. R. Schreiter, R. A. Kerr, M. B. Orger, V. Jayaraman, et al., Nature 499, 295 (2013)
2013
-
[29]
Schindelin, I
J. Schindelin, I. Arganda-Carreras, E. Frise, V. Kaynig, M. Lo ngair, T. Pietzsch, S. Preibisch, C. Rueden, S. Saalfeld, B. Schmid, et al., Nature methods 9, 676 (2012)
2012
-
[30]
Georgantzoglou, H
A. Georgantzoglou, H. Poplimont, H. A. Walker, T. L¨ ammermann, and M. Sarris, Journal of Cell Biology 221, e202103207 (2022)
2022
-
[31]
Renkawitz, A
J. Renkawitz, A. Reversat, A. Leithner, J. Merrin, and M. Sixt , in Methods in cell biology (Elsevier, 2018), vol. 147, pp. 79–91
2018
-
[32]
J. Liu, J. Boix-Campos, J. E. Ron, J. M. Kux, N. S. Gov, and P . J. S´ aez, arXiv preprint arXiv:2404.00118 (2024)
2024 arXiv
-
[33]
Mukherjee, J
A. Mukherjee, J. E. Ron, H. T. Hu, T. Nishimura, K. Hanawa-Sue tsugu, B. Behkam, Y. Mimori-Kiyosue, N. S. Gov, S. Suetsugu, and A. S. Nain, Advanced science 10, 2207368 (2023)
2023
-
[34]
J. E. Ron, M. Crestani, J. M. Kux, J. Liu, N. Al-Dam, P. Monzo , N. C. Gauthier, P. J. S´ aez, and N. S. Gov, Nature Physics pp. 1–11 (2024). 12 Fig. S-9: Illustration of the second type of slow mode during ce ll migration on the large grids. (A) A representative simulation tr...
2024
-
[226]
EP/Y02799X/1
and the European Research Council (ERC) under the Horizon 2020 program and UKR I, Grant agreement No. EP/Y02799X/1. M.S. and I.d.V acknowledge support by the European Research Counc il (grant ERC-SyG 101071793 to M.S), and thank Jack Merrin for support with microfluidic engineer ing
2020
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.