单细胞个性化分析之轨迹分析篇

作者，追风少年i

随着单细胞测序技术的飞速发展和相应分析工具的大量出现，单细胞测序的应用逐渐走向成熟，能解决的生物学问题也越来越广泛。而其中，轨迹分析（拟时分析）是单细胞分析中重要的一环，细胞轨迹分析可以通过构建细胞间的变化轨迹来重塑细胞随着时间的变化过程，帮助研究者从单细胞水平推断细胞之间的演化及分化过程，尤其在组织发育、疾病免疫细胞演化方面的研究起到了至关重要的作用。

轨迹分析的分析基础

拟时序分析，即根据不同细胞亚群基因表达量随时间的变化情况，构建细胞谱系发育，但这里的时间并不是真时间，而是一个虚拟的时间，是指的细胞与细胞之间的转化和演替的顺序和轨迹。机体为响应各种应激，其细胞会从一种功能“状态”转变为另一种功能“状态”；当细胞在不同状态之间转变时，往往会经历转录重组，导致一些基因被沉默，一些基因被重新激活，但纯化这些瞬态细胞进行研究是很困难或不可能的，而单细胞数据完美解决了这个问题。理论上当研究若干种细胞亚型，其潜在的细胞分化轨迹路径有7种（如下图）。包括环形（cycle）、线性（linear），分叉（bifurcation）、多分叉（Multifurcation）,树型（Tree），以及更复杂的连接图（相当于树结构内部又形成了回路）、分离图（图形连接中存在断点）。

图1 轨迹分化类型示意图

其中通常做的单细胞或多样本的单细胞数据，轨迹分析大部分符合上图第七种分化轨迹，即分离图，但是轨迹分析的方法很多并不能识别这种具有多种分化轨迹的分化结果，典型如monocle[1]，构建的树形图会把免疫细胞和组织细胞混合在一起，这样的结果明显是错误的。因而其实做轨迹分析最大的基础，就是首先通过细胞注释识别具有分化关系的细胞类型。

轨迹分析的主流分析思路

图2 轨迹分析的示意图

轨迹分析主要基于以下3个步骤：
（1）基因筛选：寻找以“拟时”(即不只是嘈杂)方式变化的基因，并利用这些基因来构造数据。
（2）降低维度：一旦选择了用于细胞排序的基因，就会对数据进行降维处理。
（3）pseudotime对细胞排序：通过将表达数据投影到较低维空间，构建细胞间的分化轨迹。

基因筛选

构建的轨迹分析首先是要选择用于构建轨迹的基因，当然，选择轨迹分析时用到的基因有很多方法。
（1）离散度高的基因（monocle自带的方法，默认前1000）:缺点是a、基因是否与发育相关不清楚；b、不同细胞类型的发育选取的基因数量不可能一致；c、基因断层（即基因并不是连续变化）；d、软件并不能依据轨迹基因来判断细胞是否具有多种分化路径。
（2）Seurat[2]本身挑选高变基因的三种方法（vst、mean.var.plot、dispersion）:因为Seurat降维聚类的关系，Seurat选择的高变基因也可以用于做轨迹分析，但缺点也很明显a、基因是否与发育相关不清楚；b、不同细胞类型的发育选取的基因数量不可能一致；c、基因断层；d、Seurat本身挑选的高变基因基于样本整体，没有分化关系的也纳入了分析。
（3）如果背景很强，最好的解决方式是根据生物学背景选取发育的相关基因（例如采取多样本、多时间点的策略推断发育基因，a、对比不同时间相同细胞类型的基因变化关系。b、挑选表征分化关系的基因进行轨迹分析），缺点很明显，难度特别大。
（4）寻找细胞类型之间具有连续变化的基因，理论上这是最优的选择。下图