单细胞轨迹推断与加速匹配：用动态趋势重建分化路径

单细胞转录组数据分析中轨迹推断Trajectory Inference是一类专门回答“细胞从哪里来、往哪里去”的方法。它从成千上万个细胞的静态快照中重建出分化、激活或应答的连续路径帮助研究者判断祖细胞如何进入不同命运分支。加速匹配Acceleration Matching是这类思路中的一个细分视角不仅比较细胞状态的“位置”还比较状态变化的“速度”和“加速度”让轨迹方向由动态趋势而非静态距离来约束。这篇文章从轨迹推断的基本问题讲起给出可以在自己数据上运行的最小实现并解释加速匹配的数学直觉、代码写法、验证方式和常见坑。读完以后你可以把标准轨迹推断流程跑通并用加速度修正的思路改善分支区域的图结构。1. 轨迹推断要解决什么问题加速匹配又带来什么1.1 单细胞快照为什么需要轨迹推断单细胞 RNA 测序scRNA-seq得到的是某一时刻的细胞快照。同一个组织里刚分化的细胞、处于中间状态的细胞和终末细胞可能同时存在但每类细胞都被测序成独立样本。轨迹推断要做的就是把这些无序的细胞按真实的发育进程重新排序输出一个伪时间pseudotime和一个分支结构。它不依赖实际的连续时间采样因此特别适合人类胚胎发育、造血分化、肿瘤异质性这类难以按时间全程取样的场景。在具体工具层面Monocle、SCORPIUS、Wishbone、PAGA 加 Diffusion PseudotimeDPT都是常见选择。它们的基本思路类似先用表达相似性构造细胞图再在图上寻找最短路径或扩散过程。问题是只靠“表达距离近”并不等于“发育路径相邻”尤其在分支点和高曲率区域k 近邻图很容易产生跨分支的短连接。这也是为什么轨迹推断算法不只是一个距离算法而是一个需要仔细构造图结构的问题。1.2 从“位置相似”到“变化趋势相似”标准轨迹推断使用的相似度本质上是表达状态的静态距离。比如 PCA 空间中两个细胞距离很近可能是因为它们都处于某个过渡状态也可能是因为它们共享了高表达的管家基因。静态距离无法区分“我下一步会变成你”和“我们只是长得像”。加速匹配改变了这个判断角度。它把每个细胞看成一个在表达空间中运动的点先估计这个点的局部运动方向速度再估计运动方向本身的变化加速度最后把加速度方向一致的邻居保留为轨迹边。这样选出来的边不只是“位置近”而且是“运动趋势一致”。在分化轨迹上趋势一致意味着两个细胞更容易处于同一条发育路径而非相邻但不同命运的平行轨迹。1.3 加速匹配的直觉用拐弯方向判断路径回想一下开车在一个连续的弯道里车的位置每时每刻都在变方向也在变。如果只看两个时刻车的距离可能会把对面车道的车误认为“下一个状态”。但如果同时比较速度方向是否和加速度方向一致就能自然地排除明显跑偏的候选点。单细胞轨迹推断