单细胞空间转录组联合分析硬核出击--CellTrek

自单细胞测序技术横空出世以来，由于其可以将测序缩小到单个细胞的单位内进行，对于更加全面精确地了解生物基因表达情况具有重要意义。但是，由于单细胞测序本身的技术特点，尽管我们可以得到单个细胞的表达信息，但是细胞在原本组织内的空间位置信息是完全丢失的，这些信息在理解细胞微环境或细胞间互作时同样关键。为了解决这个问题，空间转录组技术出现了，这项技术可以在组织原位上以接近单个细胞大小的尺度来进行转录组测序，但因为技术限制以及生物组织本身的复杂性，空间转录组还不能精确提供单个细胞尺度上的转录组信息。

单细胞测序技术具有单个细胞的精度，空间转录组提供了组织原位的信息，因此将两者结合起来就能近似得到完整的信息。目前已经有许多分析工具试图将空间信息与单细胞信息整合起来，并且取得了一定的进展。本文将介绍的CellTrek就是其中之一。

目前常用的整合算法有许多，例如去卷积，通过推断spot内的细胞组成来将两者整合在一起。而CellTrek则另辟蹊径，通过直接推断单细胞的空间坐标来实现信息的补全。Cell Trek的计算原理可见下图。简单来说，该方法首先将单细胞数据与空间转录组数据整合，共嵌入到（coembed）共同的隐空间内，然后从中抽出空间转录组的数据作为训练数据集，使用随机森林模型（Random Forests，RF）从spot的基因表达情况预测其空间坐标。这里所依据的基本思路就是，细胞在组织中的分布并不是随机混乱的，同一类细胞大多都是集中于同一片区域。空间转录组的数据还会同时进行非线性插值增强其空间分辨率。然后，将整个共嵌入矩阵输入训练好的模型，生成一个衡量单细胞与空间spot之间表达相似性的RF距离矩阵，基于这个矩阵使用相互最近邻域算法（mutual nearest neighbors，MNN）生成一个spot-细胞稀疏图，最后根据细胞在稀疏图上的相邻spot推断出细胞的空间坐标。因为相同类型的细胞拥有较为相似的基因表达，所以基于空间转录组数据训练的模型可以用来预测细胞的位置坐标。