跳到主要内容
返回时间线
bioRxiv来源发表:

发育轨迹并不走最小作用量路线:校准掉测量底噪后,线虫谱系仍快出约两倍

核心概要

单细胞图谱只能把每个细胞测一次,因此“发育轨迹”是跨胚胎拼接的群体质心路径,而非任何细胞的真实路径。LADD 从时序单细胞转录组学出一个随机发育向量场,用大偏差理论算出最小作用量与最短时间路径,并用三层校准把质心测量本身的代价与真实路径偏离分开。在线虫谱系图谱中,原始作用量比为 3.79,但模型自身前向模拟的集合就有 2.07,校准后仍留下约两倍的路径特异超出;斑马鱼中观测值完全落在模型底噪上。

Source-provided article image: Calibrating developmental dynamics from single-cell atlases: lineage paths outrun least-action routes
bioRxiv · 第 6 页

深度剖析

发育轨迹与推断动力学所认为的“最自然”路线之间差多少,现在有了带对照的测量:线虫谱系质心路径的作用量约为同端点最小作用量路径的 3.79 倍。 此前把发育读成最小作用量或变分自由能最小化的做法,是假定一条最优性原理,而没有先算出学习到的动力学所蕴含的最优路径、再带对照地与胚胎实际走的路径比较。 在 179 条线虫谱系分支路径上,观测质心路径的作用量比中位数为 3.79,17 种模型配置的总中位数为 3.79,嵌套自助法 95% 置信区间为 [3.56, 4.04],每种配置的中位数都高于 3.1。

这个比值不能直接读成“发育不最优”,因为把一群正在扩散、分叉的细胞取质心,本身就会带来一份与路线无关的作用量代价。 同一管线若不做校准,会在两个物种上都报出“发育比最优慢约 4–9 倍”,而其中大部分甚至全部来自测量方式。 在学习到的场下,从每条分支起始群体前向模拟 2000 个粒子得到的质心比为 2.07(四分位距 1.61–2.39),观测值显著高于这一底噪(逐路径中位倍数 1.96,Wilcoxon p = 4×10⁻³¹),又显著低于打乱谱系结构的随机划分质心 4.77(p = 3×10⁻⁵)。

线虫中校准后剩下的约两倍超出,来自谱系走得比最小作用量路线更快,而不是绕路或逆着流走。 几何与时间的分解把 87% 的对数超出归给几何因子,但集合层面的路线比较无法区分“换了路线”和“同一路线换了速度”,因此需要直接检验速度。 观测路径的遍历速度比同端点最小作用量路径快 43%(Wilcoxon p = 1×10⁻²⁷),经过的场强略弱(中位 |f| 1.80 对 1.92,p = 8×10⁻⁹),沿运动方向的场投影为正且与最小作用量路径无差别(p = 0.67);直线路线反而在 93% 的步上与场反向。

斑马鱼复现给出的是边界而非同一结论:胚后窗口的观测作用量比落在模型自身集合底噪上,因此那里的超出被归给质心几何。 线虫具有不变且高度确定的谱系,斑马鱼图谱则是阶段级时间分辨率、以 URD 分支而非不变谱系定义路径,两者不能直接比原始数值。 斑马鱼胚后窗口 25 条分支的观测比为 7.95,前向模拟集合底噪为 10.03(p = 5×10⁻³),且与随机划分质心不可区分;在统一 1 小时网格上,线虫为 2.30、斑马鱼为 2.41。

启示与展望

已展示的范围是:线虫 150–750 分钟、87,702 个细胞、50 维主成分潜空间中的 179 条谱系分支质心路径,以及斑马鱼 38,731 个细胞、胚后 4.7–12 小时窗口的 25 条 URD 分支路径;作用量是相对于学习到的自主动力学定义的量,不是适应度或能量代价的度量。作者提出的下一步是把同一框架用于蛙和小鼠图谱,前提是采用匹配的时间分辨率并与校准过的底噪比较。

作者明确说明,校准后的结论是一个测量而非机制:它确立了路径特异的超出并把直接相关量定位到速度,但这份超出是被选择还是被容忍、以及它在单细胞分辨率下的形态仍是开放问题。作用量比随潜空间维度单调上升(20 维约 3.2、50 维 3.79、100 维 5.0–5.5),作者称 50 维的选择依据的是与作用量无关的预测与谱系一致性标准。能量解读只是量级校准:与作用量率耦合的转录重塑约占胚胎 ATP 预算的 0.1%,但翻译相关成本可达周转的 10–50%,因此不能排除能量学解释。