空间转录组学简介

各位同学好,我是老周。在生物信息这行摸爬滚打了快十年,前五年做单细胞,后五年扎进空间转录组。今天咱们聊第一章——空间转录组学的来龙去脉。

说白了,空间转录组就是在 组织切片 的原始位置上,把每个细胞的基因表达量给测出来。你想想看,传统转录组把组织磨碎了测,什么位置信息都丢了。单细胞转录组虽然能分辨细胞类型,但细胞原本在组织里挨着谁、排在哪,你完全不知道。空间转录组补的就是这个短板。

我记得2019年第一次拿到10x Visium数据时,对着组织切片上那些斑点(spots)看了半天——原来肿瘤边缘的T细胞真是这样一圈一圈围着癌细胞的。那种空间上的直观感,是任何热图和聚类图都替代不了的。

从传统转录组到空间转录组

技术发展大致走了这么几步:

  • 2000年代:芯片和RNA-seq让全转录组定量成为现实,但样品得匀浆。
  • 2010年代初:单细胞RNA-seq出现,细胞类型分辨率上去了,空间位置还是盲区。
  • 2015年左右:基于原位杂交(FISH)的技术开始增多,比如MERFISH、seqFISH,能测几十到几百个基因,但通量有限。
  • 2019年:10x Genomics收购了Visium(前身是Spatial Transcriptomics),推出了第一代商业化的空间转录组解决方案。
  • 2020年后:Slide-seq、HDST、Stereo-seq等技术纷纷登场,分辨率从100μm推向亚细胞级别。
我个人的看法:现在做研究,不丢空间位置信息已经快成了刚需。你想发好点的文章,审稿人大概率会问一句“细胞在组织里的空间分布是什么样的?”。

10x Visium 原理——一张图讲清楚

Visium的核心是那张 捕获玻片,上面印满了带空间条形码(spatial barcode)的探针。流程是这样的:

  1. 新鲜冻存或FFPE组织切成10μm薄片,贴在捕获区上。
  2. 组织被透化处理,mRNA释放出来,被下方的探针捕获。每个探针带一个唯一的空间条形码,对应玻片上的某个坐标。
  3. 反转录成cDNA,建库测序。
  4. 根据条形码还原出每个spot的基因表达矩阵。

每个spot直径55μm,覆盖1~10个细胞(取决于组织密度)。分辨率不算极致,但通量和稳定性让它在应用上占了大头。

避坑指南:透化时间一定要优化好。我曾遇到过一批样本,透化过头导致mRNA扩散到相邻spot,空间信息全乱了。建议每个新组织类型都先跑透化时间梯度实验。

下面这张图总结了空间转录组技术的主要发展脉络:

空间转录组技术发展历程 传统RNA-seq 无空间信息 scRNA-seq 单细胞分辨率 MERFISH/seqFISH 原位捕获,通量有限 10x Visium 商业化空间全转录组 Slide-seq / HDST 高分辨率,~10μm Stereo-seq 亚细胞分辨率 共识:未来趋势是 更高分辨率 + 更大视野 10x Visium是目前 最广泛采用的 商业方案

核心技术参数对比

技术方法原理分辨率基因通量商业化程度
10x Visium空间条形码捕获55μm(spot)全转录组★★★★★
Slide-seqV2条形码微珠~10μm全转录组★★★☆☆
MERFISH多重FISH~200nm100~500基因★★☆☆☆
Stereo-seqDNA纳米球~0.5μm全转录组★★☆☆☆
HDST条形码微阵列~2μm全转录组★★☆☆☆
选型建议:如果组织样本充足、预算允许,我建议先用Visium跑一轮,找到感兴趣的细胞互作区域。后续再用MERFISH或Stereo-seq做高分辨率验证。Visium最大的优势是 出数据快、分析工具成熟,适合快速建立空间图谱。

应用场景,远比你想象的广

这几年我参与过的Visium项目覆盖了好几个领域,挑几个典型的说说:

  • 肿瘤微环境:最常见的应用。比如三阴性乳腺癌,看T细胞是富集在肿瘤内部还是外周基质。我遇到过一项研究用Visium发现了一个之前被忽视的“冷肿瘤-热肿瘤过渡区”,那里免疫抑制信号特别强。
  • 神经科学:大脑皮层分层、海马区空间表达模式。Visium目前分辨率在脑组织上勉强够用,但结合scRNA-seq解卷积后效果很好。
  • 发育生物学:胚胎器官发生的空间梯度。我曾帮一位合作者分析小鼠心脏发育数据,找到了心外膜到心内膜的渐进式基因表达波。
  • 疾病机制:比如肝纤维化、肾纤维化、阿尔茨海默症。空间数据能告诉你哪些区域先发生基质重塑,哪些区域代偿。

你看,有了空间位置,很多传统RNA-seq解决不了的“哪里表达”问题就迎刃而解了。

说说实战中的数据产出

用10x的Space Ranger流程处理一个标准Visium样本,最终产出几个关键文件:

sample/
├── filtered_feature_bc_matrix.h5      # 基因表达矩阵(HDF5)
├── filtered_feature_bc_matrix/         # 同上,MEX格式
├── spatial/
│   ├── tissue_hires_image.png          # 组织切片高清图
│   ├── tissue_positions_list.csv      # 每个spot的坐标与组织下
│   ├── scalefactors_json.json          # 图像缩放参数
│   └── aligned_fiducials.jpg          # fiducial标记对齐图
├── metrics_summary.csv                 # 数据质量摘要
└── analysis/                           # Space Ranger内置分析结果

嗯,这里要强调一下:拿到这些文件,第一步不是急着跑聚类,而是先检查 metric_summary.csv 里的几个关键指标。

指标正常范围危险信号
Number of Spots Under Tissue视组织大小而定低于预期(可能组织贴歪了)
Mean Reads per Spot10,000~50,000< 5,000(测序深度不足)
Median Genes per Spot1,000~3,000< 500(组织质量或透化问题)
Fraction Reads in Spots Under Tissue> 80%< 60%(污染或边缘效应严重)
我曾经踩过的坑:有一次一个样本 median genes 只有400多,大部分spots都是空转。我反复查了Space Ranger日志,最后发现是组织切片在捕获玻片上贴的时候有气泡,导致大片区域没有接触。从那以后,每次贴片我都要亲自在显微镜下确认组织平整度。

好了,这一章我们大致建立了空间转录组的全景图。下面对你来说最重要的事:记住Visium的核心价值是 在保留组织形态的前提下获得基因表达空间坐标。后面的所有分析技巧都是围绕如何用好这层空间信息展开的。