第1章 空间转录组学概览:从单细胞到空间,技术平台简介

1.1 为什么要做空间转录组?

大家好,我是老凌。做生物信息这些年,碰到的第一个灵魂拷问就是:单细胞RNA测序都这么强了,为啥还要搞空间转录组?

我给大家举个例子。你说你拿到了一套单细胞数据,找到了一个稀有细胞类型,看起来它和肿瘤免疫逃逸有关。然后呢?你能确定它在肿瘤组织里到底待在哪吗?是待在肿瘤核心区,还是躲在血管旁边?单细胞数据告诉你的是“谁在场”,但没法告诉你“谁坐在谁旁边”。

空间转录组学要解决的,就是这个核心问题:

  • 保持细胞的空间位置信息——不破坏组织切片,原位捕获RNA
  • 实现分子表达与组织形态的关联——看见基因在哪表达,才能理解它的功能

说白了,空间转录组就是把“细胞地图”和“基因表达谱”合二为一。我2019年第一次拿到10x Visium的数据时,看到基因表达信号叠加在H&E染色的组织图像上,真的觉得挺震撼的。不过,嗯,随后就掉进了数据处理的坑里……

1.2 从单细胞到空间的思路转变

这里有个概念陷阱,我刚开始做的时候也栽过。单细胞测序的输入是悬液——细胞被打散了,每个细胞独立测。而空间转录组的数据单元是捕获点(spot)或像素(pixel),一个点里可能混了几个到几十个细胞。

核心差异:单细胞 ≈ 单个细胞的快照;空间转录组 ≈ 局部细胞团的混合信号 + 位置坐标。

为什么会这样?因为目前的捕获技术,分辨率还没法做到每个点只抓一个细胞。比如10x Visium,每个点的直径是55微米,点间距100微米。一个点下覆盖的细胞数取决于组织密度,3到30个不等。

我遇到过一个小白问:那是不是可以直接拿空间数据的每个点当“伪单细胞”分析?千万别。伪单细胞不如不做,因为你根本不知道一个点里的细胞组成,做出来的差异表达很可能全是假阳性。

所以,实际分析中,我们通常会把空间数据和单细胞数据结合起来:

  1. 用单细胞数据做精细的细胞类型注释
  2. 用空间数据做细胞类型的空间分布推断
  3. 用反卷积算法(比如我们这门课的SPOTlight)把混合信号拆解开

1.3 常见空间转录组平台

做空间分析前,先得知道你手里的数据长啥样。不同平台的数据格式、分辨率、通量差别挺大的。我整理了一个对照表,方便大家快速理解:

平台 捕获方式 点直径 / 分辨率 通量(基因数) 是否有组织图像
10x Visium poly(dT) 捕获mRNA 55 μm ~2000 genes/spot H&E 或 IF
Slide-seq V2 DNA条形码微珠 10 μm ~500 genes/bead 无(需匹配)
MERFISH 单分子荧光原位杂交 单细胞级 ~100-10000 genes DAPI 或 IF
STARmap 原位测序 单细胞 / 亚细胞 ~160-1020 genes DAPI

从表里能看出来:

  • 高分辨率(MERFISH、STARmap)能做到单细胞甚至亚细胞,但通量相对受限,适合研究特定基因集
  • 高通量(10x Visium、Slide-seq)能覆盖全转录组,但分辨率较粗,必须依赖反卷积

我个人习惯的做法:如果项目目标是发现新细胞状态或新的表达程序,我会选10x Visium,因为可以测全基因组;如果是对已知的几十个基因做精细空间定位,MERFISH更合适。有一次在脑组织项目上,我坚持用Visium做全转录组探索,结果在皮层发现了一个新的边界层标记基因——如果用靶向方法肯定就错过了。

1.4 核心分析流程速览

先搭个框架,让大家对后续章节要讲什么有个心理准备。空间数据的标准处理流程大致是这样:

  1. 读入与质量控制——过滤掉低质量的点(比如空点、坏死区域)
  2. 标准化与降维——类似单细胞,用SCTransform或对数标准化
  3. 空间变量基因检测——找哪些基因有明确的空间表达模式
  4. 细胞类型反卷积——用参考单细胞数据拆解每个点的细胞组成
  5. 空间邻域分析——看细胞类型之间的空间共定位关系

我们这门课的核心——SPOTlight——就卡在第4步。它的任务很简单:给定一个空间点和一个单细胞参考矩阵,估算这个点里每种细胞类型占了多少比例。

很多人觉得反卷积就是把信号拆开,听起来像解方程。但实际操作中,参考数据的选择和预处理对结果影响巨大。一个好的参考谱,比一个好的算法更重要。我踩过这个坑,后面章节会跟大家细说。

1.5 本章知识体系结构

为了让大家更直观地理解这一章的知识脉络,我画了一张概览图。它把我们从单细胞到空间的思路转变、主流平台、以及分析流程串在了一起:

空间转录组学概览:从单细胞到空间 空间转录组学:从单细胞到空间的演化路径 单细胞 scRNA-seq 空间转录组 单细胞 → 悬液,丢失位置 空间 → 保留组织原位坐标 输出:基因×细胞矩阵 分辨率:单细胞 输出:基因×spot矩阵 + 坐标 分辨率:spot(多细胞混合) 平台:10x Visium 平台:Slide-seq 平台:MERFISH 核心分析:反卷积(SPOTlight) 用单细胞参考谱拆解空间混合信号 💡 关键点: 空间数据的高通量(全转录组) vs 高分辨率(单细胞级) 不可兼得 SPOTlight 就是为了弥补这条“分辨率鸿沟”而设计的

1.6 避坑指南 & 我的经验

最后跟大家分享几个教训,都是真金白银换来的:

我曾经踩过的一个坑:在分析一个新的肾脏空间数据时,直接用公共数据库的整个健康肾单细胞数据做参考。结果反卷积出来的细胞比例严重偏离实际情况。后来发现,我的空间样本是纤维化肾,很多成纤维细胞状态和正常肾完全不一样。

这件事给我一个教训:参考数据一定要和空间数据来自相似的组织状态和个体背景。你拿正常肝的单细胞参考去解纤维化肝的空间数据,大概率是错的。

另一个经验是:不要只看空间数据里的热点图就下结论。很多时候,一个基因看起来在某个区域高表达,可能不是因为该区域富集了高表达的细胞,而是因为那个区域细胞密度高。SPOTlight做反卷积的同时,也要结合细胞密度估计来校正。

好了,本章的概览就到这里。接下来我们会一步步深入SPOTlight的算法细节、数学原理和参数调优。大家先把这张“从单细胞到空间”的图谱记在心里,后面每讲一个环节,你都能找到它的位置。