科学问题
组织功能由哪些基因、哪些剪接异构体在什么位置表达决定,但现有手段分成了两派。测序类空间方法覆盖全转录组,捕获效率却偏低,分辨率常常到不了单细胞,读到的也多是转录本两端。成像类方法定位精度高,但每个转录本通常要几十条杂交探针,既难覆盖全转录组,也难分辨只差几十个碱基的剪接片段。本文要回答的是:能否在同一张完整组织切片上,同时拿到接近测序的覆盖度、成像的定位精度以及可变剪接信息。
背景与领域脉络
单细胞与空间组学重塑了人们对组织微环境的认识,但「看到多少基因」和「看得多准」长期不可兼得。荧光原位杂交(FISH)要靠多条探针叠加信号,探针越多信号越强,成本也越高,短转录本与相近异构体因此容易漏检。原位测序(in situ sequencing)每轮只让四分之一左右的基因发光,分子拥挤会压低检出率。已有的分支DNA(branched DNA)信号放大方案能把成像基因数推到约1.9万至2.3万个,但始终没有做出异构体分辨率。剪接异构体又恰恰决定细胞身份、信号通路与疾病表型,空间上的异构体使用图谱此前几乎是空白。
方法与技术路线
作者先把每个原始RNA分子变成一段可无限复制的模板。步骤是原位反转录时引入带T7启动子的模板转换寡核苷酸(TSO),再由T7 RNA聚合酶在靶点原地转录出高拷贝RNA扩增子。扩增子随后用多重抗错荧光原位杂交(MERFISH)读取,编码采用汉明距离4、权重4的纠错码,每轮杂交后用三(2-羧乙基)膦(TCEP)切断荧光染料,共完成60轮三色成像。基因成像使用约2.3万个基因的编码探针库,分两组测量;异构体成像另建了覆盖约4000个基因、约1万个异构体的探针库,每条异构体配4至6条编码探针。组织为成年小鼠脑冠状切片,前后脑区各两张,基因成像共约15万个细胞,异构体成像约13万个细胞。
主要结果
一条探针即可实现单分子检测
扩增把信号放大约一个数量级,单条编码探针的亮度已经超过未扩增smFISH的30条探针。在抑制性神经元标记基因Sst上,用1条或5条探针检出的Sst阳性细胞比例与未扩增、30条探针的结果相当,而非靶探针对照组几乎没有信号。
全转录组空间成像达到测序级一致性
约2.3万个基因在小鼠脑约15万个细胞中完成原位解析,覆盖脑内约85%的蛋白编码与长非编码RNA。重复切片之间基因计数高度一致,与未扩增bulk RNA-seq的表达量拟合斜率接近1,仅在极高丰度区间略呈次线性。与艾伦脑图谱(Allen Brain Atlas)的ISH能量数据比较,共有高表达基因中约80%的空间皮尔逊相关系数超过0.5。
空间表达模块刻画脑区与细胞类型
热点分析从前脑切片筛出约1万个具非随机空间分布的基因,先分出7个粗粒度解剖模块,再细分为上百个微小模块。例如含34个基因的灰被(IG)模块富集Scgn与Cabp7,含207个基因的外侧隔核吻侧部(LSr)模块除Zic家族外还富集胰高血糖素样肽-1受体基因Glp1r与钙通道基因Trpc4。作者同时用配体-受体分析重建了空间上的细胞通讯网络。
异构体使用具有明显的区域与细胞类型特异性
约1万个异构体在约13万个细胞中成像。同一基因的不同异构体常呈现不同的空间模式,部分基因被分配到多个空间模块,说明异构体而非基因本身才是空间功能的单元。Abl1的207、209、210三种异构体在小脑与皮层区域各有分布,Gfap-201与Sox10等胶质细胞异构体的主导地位也与SMART-seq2及bulk测序结果一致。
异构体特异性在不同脑结构中强弱分明
异构体表达的归一化香农熵分布显示,多数基因同时使用多个异构体,但比值差异明显。作者发现少数脑结构与细胞类型对异构体的选择性特别强,例如L5 IT、L6 IT、L5/6 NP、L6 CT与L6b等深层兴奋性神经元之间存在由浅层到深层的异构体切换。
探针数量与覆盖度之间存在可调空间
六条探针的异构体与测序结果相关性最高,四条或五条探针的异构体相关性偏低。但作者指出,由于扩增后单探针即可检测,用四探针甚至单探针策略可以把异构体库进一步扩展到更大的覆盖范围,这是后续扩容的现实路径。
创新点
核心创新是把模板转换反转录与T7原位转录放大无缝嵌进MERFISH流程,形成RT&T-AMP-MERFISH。与常规MERFISH每转录本需要数十条探针不同,该方案把单分子检测的探针需求压到1至6条,既降低探针合成成本,也让短序列与相近异构体能够在原位被区分。相对分支DNA方案只能做到基因级覆盖,本研究首次在成像平台上实现约3.3万个靶标(约2.3万基因加约1万异构体)的真全转录组级空间测量。
意义与影响
这项技术把空间组学长期存在的「靶向成像高精度」与「无偏测序全覆盖」二选一,变成可以在同一张切片上兼得。研究者因此可以在原位解析选择性剪接、可变转录起始位点与选择性加尾,而不必先把组织打碎。论文释放的小鼠脑单细胞空间全转录组与异构体图谱,为空间共表达网络算法、亚细胞RNA动力学建模以及跨模态组织表征基础模型提供了高质量基准数据集。
局限与待验证
第一,异构体成像目前只覆盖约4000个基因的约1万个异构体,相对全转录组仍是抽样,四条、五条探针的异构体相关性也明显低于六条探针。第二,与未扩增bulk RNA-seq相比,极高表达丰度处拟合斜率略低于1,提示高分子密度下仍有分子拥挤。第三,60轮三色成像对仪器连续运行时间与解码算力要求较高。第四,实验只用小鼠脑冠状切片验证,甲醛固定石蜡包埋(FFPE)等降解样本的适用性尚未证明。
对我们的相关性
对做时空组学的团队,这篇文章给出一个明确的对照参照:成像路线已经把异构体分辨率做出来,测序路线的优势仍在于无偏覆盖与组织通量,两者的互补而非替代关系更清楚了。对华大的实际启发有三点。一是原位扩增把探针成本压到十分之一量级,低成本高覆盖的成像方案值得评估。二是文中按异构体而非基因定义空间模块,提示现有的Stereo-seq数据分析流程可以补上异构体层面的空间模块划分,这与我们「剪接分析必须拿到全长」的口径一致。三是约15万个细胞、约1万个空间模块的数据组织形式,是训练组织空间基础模型时值得对齐的公开基准。
术语速查
全文 PDF 已归档 Google Drive
04-学术文献/Paper/;原文版权归 Cell 及作者所有