{"id":"b72175a3-9452-4a5d-a154-a45cb750b1f6","arxiv_id":"2506.04219","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A Hopfield-inspired model maps gene expression to cell fate coordinates and identifies straight, curved, or clustered trajectories as signatures of three decision landscape classes.","lead":"Cells commit to a fate by rolling toward stable states in an invisible landscape. This paper builds a computer model that projects single-cell genetic readouts onto cell fate coordinates, simulates three classic decision landscapes, and matches their signatures to mouse blood and lung development data. The approach could let researchers classify cell fate decisions across organs using existing single-cell datasets.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Class-specific trajectory signatures are not established as diagnostic: the same observable features can be produced by multiple decision-making classes once topographies and signal timing are varied.","rationale":"The reader's weakest assumption concerns the completeness of the reference basis used to define cell fate coordinates. That is a valid modeling assumption, but the more load-bearing issue for the central claim is that, even granting the coordinates, the qualitative matching of trajectories to simulated archetypes does not identify a landscape class unless the signatures are shown to be class-specific. The paper's own SI C concedes that trajectories are highly dependent on potential choices and that different landscapes can produce the same trajectories. This makes the experimental conclusions underdetermined. The concern is not an internal inconsistency or a disagreement with consensus; it is a missing quantitative link between model realizations and data. The proposed test—a systematic sweep over within-class topographies and signal timing, plus a quantitative model comparison against the experimental trajectories—would settle whether the signatures are diagnostic. If the test shows overlap across classes, the paper's assigned class labels are unsupported; if it shows separation, the central claim is substantially strengthened. Either way, the appropriate verdict remains conditional, as the reader concluded, pending these additional analyses.","tokens_in":27167,"tokens_out":15854,"duration_ms":139916,"concrete_test":"Systematically perturb within-class topographies and signal timing, then test whether the three signatures separate classes. For each landscape class, sample the polynomial coefficients of V_double cusp, V_triple cusp, and V_heteroclinic flip over a plausible range (e.g., +/-30% of the quoted values) and vary the signaling schedules (gap lengths, ramp durations, parameter ranges) as in Methods B. For each simulation, compute three statistics: trajectory straightness (displacement over path length), maximum curvature, and dwell time in the central multilineage region. If the distributions of these statistics overlap across classes—for instance, if a simple threshold classifier achieves near-chance accuracy—the signatures are not class-specific.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim depends on matching observed cell fate trajectories (straight basophil-neutrophil paths, curved monocyte-neutrophil paths, AT1/AT2 intermediate cluster) to simulated archetypes of the double cusp, heteroclinic flip, and triple cusp. This inference is sound only if the mapping from landscape class to trajectory shape is sufficiently discriminative—that is, if each signature appears for essentially all topographies of one class and essentially none of the others. The paper does not establish this. Instead, it simulates one hand-picked topography per class with hand-picked signaling schedules (Methods B, SI B1). The SI C acknowledges that 'the paths of transition predicted by our model are highly dependent on the choice of potential, control parameters, and other mathematical details' and that 'different landscapes can produce the exact same trajectories in certain regimes.' The triple-cusp intermediate cluster, for example, is produced by inserting a 200-step gap between the two signals; shortening that gap makes the intermediate state undetectable. Conversely, a double-cusp landscape with appropriately placed saddles, or a heteroclinic flip with slow passage through the saddle region, can generate curved paths or transient central densities. Since no quantitative comparison is made between the experimental m^mu trajectories and the model output, the observed features do not select a class. The conclusion 'consistent with' is too weak to support the assigned class labels, and the claimed signatures may be artifacts of the particular potentials and signal schedules chosen.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a phenomenological framework that links high-dimensional gene expression dynamics to low-dimensional cell fate landscapes. It generalizes modern Hopfield networks by adding signal-dependent potentials built from normal forms of elementary bifurcations, and defines scTOP generalized order parameters as cell fate coordinates. The authors simulate three decision-making classes—double cusp, triple cusp, and heteroclinic flip—and propose qualitative trajectory signatures for each. They then apply the framework to two published scRNA-seq time-series datasets (hematopoietic lineage tracing and developing mouse lung) and to a spatial patterning model of Notch-mediated airway differentiation, concluding that the measured cell fate dynamics are consistent with landscapes containing intermediate progenitors and saddle points.","tokens_in":1382,"tokens_out":1730,"duration_ms":50457,"significance":"If the proposed mapping from landscape class to trajectory shape were quantitatively validated, the framework would provide a useful bridge between bifurcation-theoretic decision classes and transcriptome-wide single-cell data, potentially enabling universal comparisons across cell fate transitions. The mathematical derivation in Eqs. 1–12 is internally consistent, the model is self-contained and uses externally defined scTOP coordinates without fitting model parameters to the experimental data, and the code and data availability statements are concrete. The spatial patterning prediction in Section II D is a falsifiable extension. However, the significance is currently limited because the central experimental inferences rest on visual pattern matching rather than quantitative tests, and the manuscript itself acknowledges that different landscapes can produce identical trajectories in certain regimes.","major_comments":[{"comment":"The central inference from experimental cell fate trajectories to specific decision-making classes is not quantitatively established. The assignments in Section II C (basophil-neutrophil as double cusp, monocyte-neutrophil as heteroclinic flip, AT1/AT2 as triple cusp) are based on visual inspection of scatter plots without error bars, null models, or alternative-class simulations. Since SI C states that 'different landscapes can produce the exact same trajectories in certain regimes' and that it is easier to eliminate classes than to identify them, the qualitative signatures described in Figure 3 do not by themselves select a class. A quantitative comparison—for example, computing summary statistics such as path curvature, intermediate-region density, or time spent in intermediate states for simulated ensembles of each class and comparing them with the experimental m^mu trajectories—is needed to support the class assignments.","section":"II C and SI C"},{"comment":"The triple-cusp interpretation of the developing lung data relies on the presence of an intermediate AT1/AT2 cluster, but the authors themselves note that this cluster could be a saddle point rather than a transient attractor. The claim that it is a transiently stable mixed-state progenitor is supported only by the hand-picked simulation schedules: SI B1 states that the prominence of the intermediate cluster depends on inserting a 200-step gap between the two signals and that shortening the gap makes the intermediate state difficult to detect. No statistical test is provided to show the cluster is significantly more populated than expected under a double-cusp or heteroclinic-flip model with a slowly passing trajectory. The 'transient AT1/AT2 mixed-state progenitor' is therefore an invented entity whose existence is not demonstrated by the data presented.","section":"II C 2 and SI B 1"},{"comment":"The model's projection onto the p-dimensional cell fate subspace is load-bearing for all three experimental conclusions, yet the manuscript does not validate that the perpendicular component x_perp is dynamically irrelevant. Equations (5)–(7) show that x_perp is static under the model dynamics, but this is a modeling assumption, not an empirical fact. If genes or regulatory programs outside the scTOP reference basis contribute to fate transitions, the observed m^mu trajectories would not faithfully reflect the true landscape. The paper should provide evidence that conclusions are robust to the choice of reference basis, for example by repeating the analysis with shuffled or augmented reference profiles and showing that the qualitative signatures persist.","section":"II A, Eq. (4)"}],"minor_comments":[{"comment":"Several cross-references are unresolved placeholders, including 'SI section??' in Section II B, 'Figure??' in the caption of Figure 4, and 'section??' in Section II A. These should be replaced with actual references before publication.","section":"Throughout"},{"comment":"The notation for the potential is confusing: V in Eq. (8) is the inverted parabola, while V with a tilde in the following paragraph is the signal-dependent term and the full potential in Eqs. (10)–(12) appears to include both. Defining the full potential explicitly and distinguishing it from the bare Hopfield potential would improve readability.","section":"II A, Eqs. (8)–(12)"},{"comment":"The analysis of hematopoietic data is restricted to clonal families with exactly two final fates, but the text does not state how many clones met this criterion or how the 'effectively bipotent' families were selected from the lineage-tracing data. Reporting the number of clones and the selection criteria would aid reproducibility.","section":"II C 1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is transparent about its limitations, which is commendable, but the abstract and discussion overstate the strength of the experimental support. As a theory-and-phenomenology paper introducing a framework, the modeling contribution is reasonable; however, the load-bearing inference from data to landscape class needs quantitative backing or a clearly weakened claim. If the authors add a quantitative discriminability analysis or reframe the conclusions as illustrative rather than evidential, the paper could be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The construction is the reason to read this paper: scTOP order parameters meet Rand–Sáez normal-form landscapes in a modern Hopfield update, giving a principled map from high-dimensional gene expression to decision-class dynamics. The derivation of Eqs. 1–12 is clean, the simulations are clearly described, and the code and data are public. That part is solid and reproducible. There is also a novel spatial-patterning extension where lateral inhibition tilts each cell's landscape; the two-class comparison (double vs triple cusp) is a nice concrete prediction.\n\nThe soft spot is the experimental inference. The signatures are proposed but not proven diagnostic. The authors simulate one topography per class with hand-picked signal schedules, and the SI explicitly concedes that paths depend heavily on the potential and that different landscapes can produce the same trajectories. The triple-cusp intermediate cluster, for instance, appears only when a 200-step gap separates the two signals; shortening the gap erases it. No quantitative comparison to the experimental m^mu trajectories is attempted, and no alternative-model baselines are fit. The straight basophil–neutrophil paths, curved monocyte–neutrophil paths, and central AT1/AT2 cluster are therefore plausible readings, not demonstrated class signatures. The phrase 'consistent with' is accurate but weak.\n\nThere is also a small mathematical slip: the text says x_perp is static, but the update rule (Eq. 5) gives dx_perp/dt = -x_perp/tau, so it decays. The m-dynamics are unaffected, so this does not change the conclusions, but it should be corrected.\n\nDespite the weak experimental leg, the framework is new, honestly discussed, and reproducible. The authors acknowledge the identification problem in the SI. What is missing is a quantitative protocol to go from trajectories to class assignments. That is exactly what a serious referee should ask for: statistical comparison of model output to data, robustness checks over topographies and signal schedules, and analysis of the full clonal families rather than only bipotent ones.\n\nI would send this to peer review. It deserves referee time, with the expectation of substantial revision. The theoretical core can stand; the empirical claims need to be either sharpened or explicitly downgraded to illustrative.","headline":"A genuinely new Hopfield–landscape construction whose experimental class signatures are interesting but not yet diagnostic.","tokens_in":28001,"tokens_out":4468,"would_cite":true,"duration_ms":42903,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["37N25","92C42"],"pacs":[],"model":"deepseek-v4-flash","headline":"Cell fate trajectories in single-cell data carry signatures of universal decision landscapes, with straight paths flagging double cusps, curved paths flagging heteroclinic flips, and mixed clusters flagging triple cusps.","keywords":["cell fate landscapes","Hopfield networks","scTOP order parameters","single-cell RNA-seq","decision-making classes","bifurcation theory","Waddington landscape","hematopoiesis"],"falsifier":"Re-analyze the lineage-traced hematopoiesis and lung data while tracking the perpendicular component $x_i^\\perp(t) = x_i(t) - \\sum_\\mu m^\\mu(t)\\,\\xi_{\\mu i}$; if this component changes systematically during any of the three transitions, or if adding a few reference profiles orthogonal to the current basis turns a straight trajectory curved, the projection assumption fails and the class assignments are not trustworthy.","tokens_in":26955,"feed_emoji":"🧬","tokens_out":5862,"duration_ms":55886,"temperature":0.7,"pith_summary":"Single-cell RNA-sequencing time series can be read as trajectories in a low-dimensional cell fate space, and this paper tries to show that the shapes of those trajectories reveal which universal decision-making class a cell fate transition belongs to. Using generalized Hopfield order parameters as coordinates, the authors simulate archetypal trajectories for three landscape classes and then match them to experimental data: straight basophil-to-neutrophil paths, curved monocyte-to-neutrophil paths, and a mixed AT1/AT2 cluster in the developing lung are consistent with a double cusp, a heteroclinic flip, and a triple cusp, respectively. If the correspondences hold, cell fate decisions can be classified directly from scRNA-seq time series, and the same landscape classes can explain spatial patterning through lateral inhibition signals.","feed_headline":"Cell fate turns read out as cusp and flip landscapes","feed_subtitle":"Straight, curved, and clustered trajectories in mouse blood and lung cells match three predicted classes.","key_machinery":"The load-bearing object is the generalized Hopfield order parameter $m^\\mu = \\sum_\\nu (A^{-1})^{\\mu\\nu} m_\\nu$, with $A_{\\mu\\nu} = \\sum_i \\xi_{\\mu i} \\xi_{\\nu i}$, which projects a cell's gene expression vector onto the subspace spanned by reference cell type profiles and defines cell fate coordinates. The dynamics are $\\tau \\, dm^\\mu/dt = \\sigma^\\mu(-\\beta \\, \\partial V/\\partial m^\\mu) - m^\\mu$, where $\\sigma^\\mu$ is a softmax nonlinearity and the potential $V$ combines an inverted parabola with a signal-dependent landscape $\\tilde V$ built from normal forms of elementary bifurcations. The signaling parameters $f, k$ tilt the landscape and destabilize attractors, producing the bifurcations that move cells from one fate to another.","core_discovery":"The paper's central claim is that the measured cell fate dynamics are consistent with developmental landscapes containing intermediate progenitors and saddle points, and that the specific trajectory geometry encodes the class. In cell fate coordinates, straight trajectories to final fates indicate a double cusp with no mixed state; curved trajectories through multilineage states indicate a heteroclinic flip with an unstable manifold; and a dense intermediate cluster between fates indicates a triple cusp with a transiently stable progenitor. Applied to lineage-traced hematopoiesis and developing lung alveolar cells, the data show all three signatures. The same landscape logic is used to distinguish two Notch-dependent models of airway patterning, where a single signal corresponds to a double cusp and two sequential signals to a triple cusp.","pith_inferences":["An editor's extension: the straight-versus-curved distinction may serve as a practical lineage-relationship diagnostic, with straight paths suggesting distantly related fates and curved paths suggesting adjacent fates connected by a saddle or progenitor.","The projection-based coordinates depend on the reference basis chosen; systematically perturbing the basis, for example by adding or removing closely related cell types, would quantify how robust each class assignment is.","Because the same potential can be tilted by different signaling schedules, controlled in vitro differentiation with measured time courses could be used to decide between competing landscape classes for the same pair of fates.","The framework could be extended to multi-fate decisions by building higher-dimensional normal forms, though the paper restricts itself to three-attractor classes."],"forward_implications":["If the correspondences hold, scRNA-seq time series can be classified by trajectory shape alone, without fitting a full landscape or choosing marker genes by hand.","The triple-cusp assignment for alveolar maturation implies a transient AT1/AT2 progenitor that is stable before birth and destabilized when air breathing begins, a testable in vitro prediction.","In airway injury, the two Notch-signaling models predict different commitment timing: double cusp commits cells early, triple cusp keeps them in a mixed state until a second signal.","The method is designed to scale to atlas-level data, so the same signatures could be used to survey developmental transitions across organs and species.","Different classes predict different signal sensitivities: heteroclinic flips are the most sensitive to fate-biasing signals, whereas double cusps resist multilineage expression."],"supporting_citations":[{"why":"Supplies the mathematical decision-making classes built from fold and heteroclinic flip bifurcations used to construct the landscapes.","marker":"[23]"},{"why":"Grounds the claim that fold bifurcations and heteroclinic flips suffice to connect Morse-Smale systems, justifying the landscape catalog.","marker":"[24]"},{"why":"Defines the scTOP generalized Hopfield order parameters used as cell fate coordinates and the reference-basis procedure.","marker":"[29]"},{"why":"Provides the lineage-traced hematopoietic scRNA-seq data whose straight and curved trajectories are matched to double cusp and heteroclinic flip classes.","marker":"[53]"},{"why":"Provides the embryonic-to-postnatal mouse lung scRNA-seq time series exhibiting the AT1/AT2 mixed cluster.","marker":"[56]"},{"why":"Establishes statistically derived geometrical landscapes for decision-making classes used as the comparison baseline for trajectory signatures.","marker":"[21]"},{"why":"Provides an algorithm for inferring landscapes from continuous coordinates, which the paper's cell fate coordinates could plug into.","marker":"[26]"},{"why":"Shows that landscapes do not fully specify dynamics and provides inference methods for fitting class probabilities.","marker":"[27]"},{"why":"Supplies the two Notch-signaling models, single versus two signals, that the spatial patterning simulations map onto double cusp versus triple cusp.","marker":"[58]"}],"fun_headline_variants":["Cell fate paths expose hidden decision landscapes","Trajectory shapes reveal cell fate decision class","Single-cell paths map gene decision landscapes","Cusp and flip signatures in cell fate data","Cell fate dynamics decode landscape classes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis assumes that the only gene expression relevant to a fate decision is the part that lies in the subspace spanned by the chosen reference cell types; the perpendicular component is treated as static and ignored, so if genes or states outside that reference basis drive a transition, the observed trajectories would not reflect the true landscape.","fun_headline_variants_meta":{"raw":{"variants":["Cell fate paths expose hidden decision landscapes","Trajectory shapes reveal cell fate decision class","Single-cell paths map gene decision landscapes","Cusp and flip signatures in cell fate data","Cell fate dynamics decode landscape classes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000133,"raw_usage":{"total_tokens":1104,"prompt_tokens":880,"completion_tokens":224,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":496,"completion_tokens_details":{"reasoning_tokens":160}},"tokens_in":496,"tokens_out":224,"duration_ms":2908,"temperature":1.0,"reasoning_tokens":160,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T10:45:54.151690+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-analyze the lineage-traced hematopoiesis and lung data while tracking the perpendicular component $x_i^\\perp(t) = x_i(t) - \\sum_\\mu m^\\mu(t)\\,\\xi_{\\mu i}$; if this component changes systematically during any of the three transitions, or if adding a few reference profiles orthogonal to the current basis turns a straight trajectory curved, the projection assumption fails and the class assignments are not trustworthy.","supporting_citations":[{"cited_title":"S´ aez, R","cited_arxiv_id":null,"evidence_quote":"Supplies the mathematical decision-making classes built from fold and heteroclinic flip bifurcations used to construct the landscapes."},{"cited_title":"Wang, Journal of Biological Physics48, 1 (2022)","cited_arxiv_id":null,"evidence_quote":"Grounds the claim that fold bifurcations and heteroclinic flips suffice to connect Morse-Smale systems, justifying the landscape catalog."},{"cited_title":"Mochulska and P","cited_arxiv_id":null,"evidence_quote":"Defines the scTOP generalized Hopfield order parameters used as cell fate coordinates and the reference-basis procedure."},{"cited_title":"Bialek,Biophysics: searching for principles(Princeton University Press, 2012)","cited_arxiv_id":null,"evidence_quote":"Provides the lineage-traced hematopoietic scRNA-seq data whose straight and curved trajectories are matched to double cusp and heteroclinic flip classes."},{"cited_title":"Alysandratos, M","cited_arxiv_id":null,"evidence_quote":"Provides the embryonic-to-postnatal mouse lung scRNA-seq time series exhibiting the AT1/AT2 mixed cluster."},{"cited_title":"Camacho-Aguilar, A","cited_arxiv_id":null,"evidence_quote":"Shows that landscapes do not fully specify dynamics and provides inference methods for fitting class probabilities."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the two Notch-signaling models, single versus two signals, that the spatial patterning simulations map onto double cusp versus triple cusp."}],"review_version":1}