REVIEW 4 major objections 5 minor 38 references
A non-reciprocal model for morphogenesis in symbiosis
T0 review · 4 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read A fluid vesicle and hard-sphere symbionts with unequal two-way interaction strengths generate protrusions, invaginations, budding, cyclic blebs, and a ballistically moving polar vesicle.
desk verdict New non-reciprocal membrane-remodeling mechanism with solid simulations, but the abstract's uniqueness claim needs a same-model reciprocal-active control before publication. 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 the symbiotic force $\vec f_S$, defined as the net force on the cluster of membrane beads and symbiont particles within interaction range. In the reciprocal limit $\Delta\varepsilon=0$, this force cancels by action-reaction; when $\Delta\varepsilon\neq 0$, it is nonzero, grows with $|\Delta\varepsilon|$ and with cluster size, and carries a direction set by the sign of $\Delta\varepsilon$. Its dependence on the wrapping geometry closes the feedback loop: as the membrane deforms, the cluster changes and $\vec f_S$ changes with it. The machinery is implemented as a dynamically triangulated membrane with constrained area and volume, Monte Carlo bond swaps that make the membrane fluid, and a Langevin thermostat that absorbs the unbalanced momentum.
What would settle it
Re-run the triangulated-membrane simulations with $\Delta\varepsilon=0$ and apply to the symbiont an external constant traction equal to the steady-state symbiotic force measured at finite $\Delta\varepsilon$; if the same morphology diagram (protrusions, cyclic blebs, ballistic polar vesicle) still appears, the claim that shape-dependent non-reciprocal feedback drives these morphologies is falsified.
Extended reading notes
Core claim
The central claim is that the non-reciprocal pairwise force law $\vec F_{k,i\to l,j}=-\phi\varepsilon_{ij}/r^\gamma\,\hat r_{kl}$ for $r<r_c$, with $\varepsilon_{sh}\neq\varepsilon_{hs}$, turns passive contact into an active shape generator for fluid membranes. For $\Delta\varepsilon>0$ the net force on the host-symbiont cluster points toward the symbiont and, once $\Delta\varepsilon/k_BT\gtrsim 5$, extrudes dynamic protrusions; for $\Delta\varepsilon<0$ the force instead pushes symbionts into wrapped, invaginated, or budded states, and with several symbionts produces cyclic blebs and a polar vesicle that moves ballistically. The key point is that the symbiotic force $\vec f_S$ is not fixed: it depends on the geometry of the cluster, so the local membrane deformation alters the very force that drives it. The paper presents this shape-force feedback as the mechanism behind morphologies that are not accessible to reciprocal membranes or to membranes deformed by self-propelled particles.
Load-bearing premise
The model stands on the assumption that a symbiont can push on a host more strongly than the host pushes back, as an instantaneous pairwise force with the surrounding fluid absorbing the missing momentum; if real host-symbiont coupling is mechanically reciprocal or mediated mainly by chemical fields, the predicted shapes would not arise.
Editorial extensions
If this is right
- A single symbiont with $\Delta\varepsilon/k_BT\gtrsim 5$ extrudes a persistent membrane protrusion even though no agent explicitly pulls the membrane.
- For $\Delta\varepsilon<0$, symbionts become wrapped, invaginated, or budded by the host, and at intermediate symbiont numbers the membrane produces blebs in a cyclic sequence.
- The number of symbionts controls the mode: independent protrusions at low coverage, cooperative branched protrusions at intermediate coverage, and suppression by membrane tension at high coverage.
- Asymmetric symbiont distributions can create a polar vesicle that moves ballistically.
- These regimes are not found in reciprocal systems with constant intrinsic activity, so their observation would indicate a genuinely interaction-derived driving force.
Reading between the lines
- Because $\gamma=2$ maps the same force law onto phoretic interactions between catalytically coated colloids, the model suggests a direct experimental test: two colloid species with independently tunable surface chemistry paired with a giant vesicle should show protrusion for one sign of $\Delta\varepsilon$ and invagination or blebbing for the other; this goes beyond the paper, which does not propo
- The shape-dependent $\vec f_S$ could be coarse-grained into a local active tension or curvature-dependent force in a continuum membrane theory, letting future work predict phase boundaries for cyclic blebbing or polar motion without resolving every particle; the authors do not make this connection.
- If the interaction strengths track metabolic exchange rates, the results imply that a biological cell could switch between protrusive, invaginated, and blebbing morphologies purely by tuning what it gives to and receives from a partner, without invoking cytoskeletal machinery; this is an inference from the authors' discussion, not a claim they test.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. Muñoz-Basagoiti et al. introduce a coarse-grained model in which a host vesicle (dynamically triangulated membrane) and symbiont particles interact via pairwise forces that violate action–reaction symmetry. The force on a particle of species j from species i is given by Eq. (1), with strength ε_ij; the asymmetry is parameterized by Δε = ε_sh − ε_hs. The authors show by molecular dynamics/HMC simulations that a single symbiont produces stable protrusions for Δε > 0 and invaginations, catenoids, or budding for Δε < 0, and that the net 'symbiotic force' fS on the interacting cluster is proportional to |Δε| and depends on the cluster geometry, establishing a feedback between membrane shape and drive. With multiple symbionts, they report cooperative and suppressed protrusions for Δε > 0, and for Δε < 0 cyclic blebbing, ballistic polar vesicles, and 'inside-out' protrusions. Results are checked with a second membrane model (one-particle-thick) and against an equilibrium wrapping limit (γ=7, φ=6/σ̃, Δε=0) and a tether-pulling control.
Significance. If the results hold, the paper offers a minimal mechanism—non-reciprocal interaction alone—for generating a rich set of membrane morphologies relevant to symbiosis, and may inspire synthetic cell experiments. Strengths: explicit non-reciprocal pair style in LAMMPS; two independent membrane models; multiple replicas for the single-symbiont phase diagram; clear definition of fS and its shape feedback; code availability. The central open issue is the abstract's uniqueness claim relative to reciprocal active systems, which is not tested with a same-model baseline.
major comments (4)
- [Abstract and Sec. III.C] The claim that the observed morphologies are 'not reported in reciprocal systems with constant activity' and 'not accessible to systems which rely on intrinsic force production via motile filaments or other self-propelled agents' is not supported by a same-model control. The only quantitative control (SI Fig. S5) uses a reciprocal interaction plus an external constant pulling force on a tethered symbiont; this is not equivalent to a self-propelled (active Brownian) particle. No simulation with ε_hs = ε_sh and a constant self-propulsion force on the symbionts is reported. Please add such simulations (e.g., active Langevin dynamics of symbionts with reciprocal membrane interactions) at overlapping parameters, or soften the claims to 'not observed in our reciprocal control' and 'differ from previously reported active-membrane studies.' Without this, the abstract's novelty statement is not verifiable from the presented data.
- [Sec. III.A and Fig. 2A] The protrusion region is defined by heuristic thresholds (p ≥ 30% protrusion probability, acylindricity a ≤ 1, cutoff d = 5σ) and the phase boundaries are only guides for the eye. The sensitivity of the phase diagram to these thresholds is not characterized. Since the existence of a distinct protrusion regime is a central result, please report the underlying probability distributions or a threshold-sensitivity analysis so the reader can judge whether the boundaries correspond to sharp or gradual crossovers.
- [Sec. III.C, Figs. 4B and 5] The multi-symbiont phase diagrams are based on only three replicas per parameter set, yet they identify several dynamical phases (cyclic blebbing, polar vesicle, inside-out protrusions). For stochastic classification of transient morphologies, three replicas is low, especially for distinguishing cyclic blebbing from stochastic membrane fluctuations. Please provide additional replicas for the dynamic phases or give quantitative order parameters with error bars, and state explicitly how 'cyclic' blebbing is defined and measured.
- [Sec. II.C and Sec. III.B] The model relies on the Langevin thermostat to absorb the momentum imbalance produced by the non-reciprocal forces. The sensitivity of the results to the thermostat damping coefficient and to simulation box size is not reported, and these choices could affect the persistence or magnitude of the symbiotic force fS. Please add a short parameter check (e.g., two damping values and two box sizes) to show that the morphologies and the fS scaling in Fig. 3B are robust to the thermostatting protocol.
minor comments (5)
- [Sec. III.C and Fig. 5 caption] The caption of Fig. 5 states ε_sh/kBT = 2 while Sec. III.C states ε_hs/kBT = 2; please clarify which interaction strength is held fixed in the Δε < 0 multi-symbiont simulations, since the wrapping state and hence the available symbiotic force depend on it.
- [Fig. 2A caption and Sec. III.A] The caption refers to 'white data points' while the text later refers to 'white diamond symbols' for the same region; please make the symbol terminology consistent.
- [SI Fig. S5 caption] The caption calls the control 'reciprocal and non-reciprocal pulling'; the control is actually a reciprocal interaction plus an external force, which is a different driving protocol. Please rephrase to avoid implying that it is an active self-propulsion baseline.
- [Abstract and Introduction] The sentence 'The shape of a cell influences, and it is influenced by its interactions with its neighbours' is grammatically awkward; consider rephrasing for clarity.
- [Sec. II and SI Sec. I] The initial placement of symbionts (distance from the vesicle, equilibration protocol) is not fully specified; please add the exact initialization details used for the single- and multi-symbiont simulations.
Circularity Check
No significant circularity: the model's morphologies are simulation outputs of the stated non-reciprocal interaction, and the one definitional relation (fS proportional to |Δε|) is explicitly acknowledged by the authors.
full rationale
The paper's derivation chain is a simulation study: Eq. (1) defines a non-reciprocal pairwise force, and the 'symbiotic force' fS is then defined as the net force on the host-symbiont cluster. The authors explicitly state that the proportionality of fS to |Δε| follows by definition and that the shape-force coupling is 'by design' (Section IV, Discussion). These are model characterizations, not fitted predictions or disguised inputs; the morphology diagrams in Figs. 2, 4, and 5 are simulation outputs rather than fits to target data. The reciprocal-equilibrium comparison and the one-particle-thick membrane cross-check in the SI provide independent controls. The only caveat is the abstract's uniqueness claim that the morphologies are 'not reported in reciprocal systems with constant activity': this rests on a literature comparison and on SI Fig. S5, which uses an external tether pull rather than a reciprocal-interaction-plus-constant-activity baseline. That is an evidentiary/completeness gap in the comparative claim, but it is not a circular reduction, because the central non-reciprocal phenomenology is not obtained by fitting or by defining the output into the input. Self-citations concerning the simulation method and equilibrium adhesion benchmarks are not load-bearing for the non-reciprocal results.
Assumptions & free parameters
free parameters (8)
- epsilon_hs (host-on-symbiont force strength) =
scanned; e.g., 10 kBT in protrusion runs, 2 kBT for Delta-epsilon < 0 diagrams
- Delta-epsilon = epsilon_sh - epsilon_hs =
varied from -12 to +12 kBT
- gamma (force exponent) =
7 (checked at 2)
- phi (force scale) =
6 / sigma_tilde
- r_c (interaction cutoff) =
1.5 sigma_tilde
- sigma_s (symbiont diameter) =
5 sigma, 10 sigma, 20 sigma
- kappa_A,V (area/volume constraint stiffness) =
2.5e5 kBT
- kappa_B (bending rigidity) =
20 kBT
assumptions (5)
- domain assumption Non-reciprocal pairwise forces can drive particle dynamics when the momentum imbalance is absorbed by an implicit solvent (Langevin thermostat).
- domain assumption The host cell can be modeled as a fluid vesicle with fixed topology; topological changes (budding, rupture) are secondary for the reported regimes.
- domain assumption Symbiont biology is representable by hard spheres with pairwise distance-dependent forces, without chemical fields or internal degrees of freedom.
- ad hoc to paper The parameter mapping gamma = 7 and phi = 6/sigma_tilde connects the non-reciprocal force to Lennard-Jones adhesion, so Delta-epsilon = 0 recovers equilibrium wrapping.
- ad hoc to paper Morphology classification thresholds (30% protrusion probability, d = 5 sigma, acylindricity <= 1, N/Nmax >= 0.9 for budding) correctly separate the reported phases.
Cite this review
Pith. "Pith review of A non-reciprocal model for morphogenesis in symbiosis." pith.science (2026). https://pith.science/paper/GWTNE3A6
@misc{pith2026250613299,
author = {Pith},
title = {Pith review of: A non-reciprocal model for morphogenesis in symbiosis},
year = {2026},
howpublished = {\url{https://pith.science/paper/GWTNE3A6}},
note = {Machine review of arXiv:2506.13299}
}
read the original abstract
The shape of a cell influences, and it is influenced by its interactions with its neighbours. Here, we introduce a coarse-grained computational model of non-reciprocal interactions between single-cell organisms to study emergent morphologies during symbiotic association. We show that the cell membrane can be remodelled into branched protrusions, invaginations, transient blebs and other dynamical morphologies that depend on the number of interacting partners, the asymmetry, and the magnitude of partnership activity. Our model finds a dynamical feedback between the local deformation of the membrane and its driving force, leading to membrane morphologies not reported in reciprocal systems with constant activity.
Figures
Reference graph
Works this paper leans on
-
[1]
P. Nardon and H. Charles, Morphological Aspects of Symbiosis, inSymbiosis: Mechanisms and Model Sys- tems, edited by J. Seckbach (Springer Netherlands, Dordrecht, 2002) pp. 13–44
work page 2002
-
[2]
A. K. Baidya, S. Bhattacharya, G. P. Dubey, G. Mamou, and S. Ben-Yehuda, Bacterial nanotubes: a conduit for intercellular molecular trade, Cur- rent Opinion in Microbiology Cell Regulation,42, 1 (2018)
work page 2018
-
[3]
M. Parniske, Arbuscular mycorrhiza: the mother of plant root endosymbioses, Nature Reviews Microbi- ology 6, 763 (2008), publisher: Nature Publishing Group
work page 2008
-
[4]
J. C. Cross, H. Nakano, D. R. C. Natale, D. G. Sim- mons, and E. D. Watson, Branching morphogenesis during development of placental villi, Differentiation 74, 393 (2006)
work page 2006
-
[5]
B. S. Leander and P. J. Keeling, Symbiotic Innova- tion in the OxymonadStreblomastix strix, Journal of Eukaryotic Microbiology51, 291 (2004)
work page 2004
-
[6]
L. Margulis, Symbiosis and Evolution, Scientific American 225, 48 (1971), publisher: Scientific Amer- ican, a division of Nature America, Inc
work page 1971
-
[7]
D. A. Baum and B. Baum, An inside-out origin for the eukaryotic cell, BMC Biology12, 76 (2014)
work page 2014
-
[8]
M. Fruchart, R. Hanai, P. B. Littlewood, and V. Vitelli, Non-reciprocal phase transitions, Nature 592, 363 (2021), publisher: Nature Publishing Group
work page 2021
Show all 38 references
-
[9]
S. H. L. Klapp, Non-reciprocal interaction for living matter, Nature Nanotechnology 18, 8 (2023), pub- lisher: Nature Publishing Group
2023
-
[10]
R.SotoandR.Golestanian,Self-AssemblyofCatalyt- ically Active Colloidal Molecules: Tailoring Activity Through Surface Chemistry, Physical Review Letters 112, 068301 (2014), 1306.6596
2014 arXiv
-
[11]
C. H. Meredith, P. G. Moerman, J. Groenewold, Y.- J. Chiu, W. K. Kegel, A. van Blaaderen, and L. D. Zarzar, Predator–prey interactions between droplets driven by non-reciprocal oil exchange, Nature Chem- istry 12, 1136 (2020), publisher: Nature Publishing Group
2020
-
[12]
S. Saha, J. Agudo-Canalejo, and R. Golestanian, Scalar Active Mixtures: The Nonreciprocal Cahn- Hilliard Model, Physical Review X10, 041009 (2020), 7 publisher: American Physical Society
2020
-
[13]
Brauns and M
F. Brauns and M. C. Marchetti, Nonreciprocal Pat- tern Formation of Conserved Fields, Physical Review X 14, 021014 (2024), publisher: American Physical Society
2024
-
[14]
Osat and R
S. Osat and R. Golestanian, Non-reciprocal multifar- ious self-organization, Nature Nanotechnology18, 79 (2023), publisher: Nature Publishing Group
2023
-
[15]
Ouazan-Reboul, J
V. Ouazan-Reboul, J. Agudo-Canalejo, and R. Golestanian, Self-organization of primitive metabolic cycles due to non-reciprocal interactions, Nature Communications 14, 4496 (2023), publisher: Nature Publishing Group
2023
-
[16]
S. A. Loos, S. H. Klapp, and T. Martynec, Long- Range Order and Directional Defect Propagation in the Nonreciprocal XY Model with Vision Cone Inter- actions, Physical Review Letters130, 198301 (2023), publisher: American Physical Society
2023
-
[17]
H. R. Vutukuri, M. Hoore, C. Abaurrea-Velasco, L. v. Buren, A. Dutto, T. Auth, D. A. Fedosov, G. Gomp- per, and J. Vermant, Active particles induce large shapedeformationsingiantlipidvesicles,Nature 586, 52 (2020)
2020
-
[18]
P. Iyer, G. Gompper, and D. A. Fedosov, Non- equilibrium shapes and dynamics of active vesicles, Soft Matter18, 6868 (2022)
2022
-
[19]
P. Iyer, G. Gompper, and D. A. Fedosov, Dynamic shapes of floppy vesicles enclosing active Brown- ian particles with membrane adhesion, Soft Matter 10.1039/d3sm00004d (2023)
2023 doi
-
[20]
Siggel, S
M. Siggel, S. Kehl, K. Reuter, J. Köfinger, and G. Hummer, TriMem: A parallelized hybrid Monte Carlo software for efficient simulations of lipid mem- branes, The Journal of Chemical Physics157, 174801 (2022)
2022
-
[21]
Muñoz-Basagoiti, F
M. Muñoz-Basagoiti, F. Frey, B. Meadowcroft, M. Amaral, A. Prada, and A. Šarić, A tutorial for mesoscale computer simulations of lipid membranes: tether pulling, tubulation and fluctuation (2025), arXiv:2502.09798
2025 arXiv
-
[22]
Deserno, Elastic deformation of a fluid membrane upon colloid binding, Physical Review E69, 031903 (2004), publisher: American Physical Society
M. Deserno, Elastic deformation of a fluid membrane upon colloid binding, Physical Review E69, 031903 (2004), publisher: American Physical Society
2004
-
[23]
Šarić and A
A. Šarić and A. Cacciuto, Self-assembly of nanopar- ticles adsorbed on fluid and elastic membranes, Soft Matter 9, 6677 (2013), publisher: Royal Society of Chemistry
2013
-
[24]
Šarić and A
A. Šarić and A. Cacciuto, Fluid Membranes Can Drive Linear Aggregation of Adsorbed Spherical Nanoparticles, Physical Review Letters108, 118101 (2012), 1201.0036
2012 arXiv
-
[25]
S. A. Levin,The Princeton Guide to Ecology (Prince- ton University Press, 2012) 9780691156040
2012
-
[26]
H. Yuan, C. Huang, J. Li, G. Lykotrafitis, and S. Zhang, One-particle-thick, solvent-free, coarse- grained model for biological and biomimetic fluid membranes, Physical Review E82, 011905 (2010)
2010
-
[27]
Derényi, F
I. Derényi, F. Jülicher, and J. Prost, Formation and Interaction of Membrane Tubes, Physical Review Let- ters 88, 238101 (2002), publisher: American Physical Society
2002
-
[28]
Rodrigues-Oliveira, F
T. Rodrigues-Oliveira, F. Wollweber, R. I. Ponce- Toledo, J. Xu, S. K.-M. R. Rittmann, A. Klingl, M. Pilhofer, and C. Schleper, Actin cytoskeleton and complex cell architecture in an Asgard archaeon, Na- ture 613, 332 (2023)
2023
-
[29]
Muller, T
F. Muller, T. Brissac, N. Le Bris, H. Felbeck, and O. Gros, First description of giant Archaea (Thaumarchaeota) associated with putative bacterial ectosymbionts in a sulfidic marine habitat, Envi- ronmental Microbiology 12, 2371 (2010), _eprint: https://onlinelibrary.wiley.com...
2010
-
[30]
A. H. Bahrami and G. Hummer, Formation and Sta- bility of Lipid Membrane Nanotubes, ACS Nano11, 9558 (2017), publisher: American Chemical Society
2017
-
[31]
Hamada, K
M. Hamada, K. Schröder, J. Bathia, U. Kürn, S. Fraune, M. Khalturina, K. Khalturin, C. Shin- zato, N. Satoh, and T. C. Bosch, Metabolic co- dependence drives the evolutionarily ancient Hy- dra–Chlorella symbiosis, eLife7, e35122 (2018), pub- lisher: eLife Sciences Publications, Ltd
2018
-
[32]
V. C. Clarke, P. C. Loughlin, D. A. Day, and P. M. C. Smith, Transport processes of the legume symbiosome membrane, Frontiers in Plant Science5, 10.3389/fpls.2014.00699 (2014), publisher: Frontiers
2014
-
[33]
S. Shim, S. Khodaparast, C.-Y. Lai, J. Yan, J. T. Ault, B. Rallabandi, O. Shardt, and H. A. Stone, CO 2 -Driven diffusiophoresis for maintaining a bacteria-free surface, Soft Matter 17, 2568 (2021), publisher: Royal Society of Chemistry
2021
-
[34]
Imachi, M
H. Imachi, M. K. Nobu, N. Nakahara, Y. Morono, M. Ogawara, Y. Takaki, Y. Takano, K. Uematsu, T. Ikuta, M. Ito, Y. Matsui, M. Miyazaki, K. Mu- rata, Y. Saito, S. Sakai, C. Song, E. Tasumi, Y. Ya- manaka, T. Yamaguchi, Y. Kamagata, H. Tamaki, and K. Takai, Isolation of an archae...
2020
-
[35]
Spang, J
A. Spang, J. H. Saw, S. L. Jørgensen, K. Zaremba- Niedzwiedzka, J. Martijn, A. E. Lind, R. van Eijk, C. Schleper, L. Guy, and T. J. G. Ettema, Complex archaea that bridge the gap between prokaryotes and eukaryotes, Nature 521, 173 (2015), publisher: Na- ture Publishing Group
2015
-
[36]
Y. Ji, T. Chakraborty, and S. V. Wegner, Self- Regulated and Bidirectional Communication in Syn- thetic Cell Communities, ACS Nano17, 8992 (2023)
2023
-
[37]
Wassermair, M
M. Wassermair, M. Muñoz-Basagoiti, M. Amaral, and A. Saric, TriLMP (TriMEM + LAMMPS): Dynami- cally Triangulated Networks in LAMMPS
-
[38]
A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. i. t. Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton, LAMMPS - a flexible simulation tool for particle-based mater...
2022
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.