空间转录组三大核心分析模式：细胞映射、结构域识别与通讯推断-北京尧图网络

空间转录组这几年在科研圈的热度不用我多说几乎每个做组织微环境、发育生物学、肿瘤免疫的课题组都在往这个方向上靠。但很多人拿到数据后第一反应是懵——这玩意到底怎么分析它跟单细胞转录组到底有什么区别我到底该往哪个方向挖掘去年我前后跑了将近二十套空间转录组数据从一开始对着Seurat的教程硬啃到后来能够针对不同的科学问题快速搭出分析方案中间踩了不少坑也总结出一些比较稳定的思路。今天就把我对空间转录组三类核心分析模式的理解整理出来希望对准备入门或者已经被数据折磨得头秃的朋友有点帮助。先说一句题外话空间转录组本质上是在单细胞转录组的基础上多给了你一个坐标。有了坐标意味着什么意味着你可以回答过去单细胞数据回答不了的问题——那些细胞在组织里到底待在哪、旁边是谁、它们之间有没有物理接触、不同区域在转录水平上有多大的差异。说得直白点单细胞数据是告诉你“有哪些演员、各自说了什么台词”空间转录组是额外告诉你“舞台上每个人站在哪个位置、谁跟谁挨着”。但位置带来便利的同时也带来了分析上的复杂度这也是为什么会出现不同的分析模式。每一种模式的出发点、算法逻辑、适用场景都不一样没有一套代码能同时解决所有问题。你要做的第一件事就是想清楚你的生物学问题到底属于哪一类。1. 空间转录组数据的基本盘数据结构和预处理在拆三种分析模式之前我必须先把数据预处理这部分讲清楚因为不管你后续走哪条分析路线这一步都是地基。很多时候分析结果不理想问题不是出在后面的算法选择上而是前面的数据过滤跟标准化就没做对。空间转录组的原始数据长什么样以最主流的10x Visium平台为例一张切片上排列着大约5000个spot斑点每个spot的直径是55微米spot之间的中心距离是100微米。一个spot里捕获的是周围多个细胞的混合转录本理论上包含1到10个细胞不等具体取决于组织类型和细胞密度。这意味着什么空间转录组的每一个spot本质上是一个“小规模的bulk转录组”它的分辨率是“区域级”而不是“单细胞级”的。这一点贯穿所有分析模式无论你后面做什么都不能忘了这个前提。数据预处理的核心步骤和单细胞基本一致包括质量控制、标准化、特征选择、降维和聚类。但在实际操作中有几个空间数据特有的坑需要特别注意。第一个坑是线粒体基因比例的问题。在单细胞数据里线粒体基因比例高通常意味着细胞状态差或者正在凋亡所以我们会设定一个阈值比如20%把这种细胞过滤掉。但在空间转录组里情况不完全一样——如果一个spot正好落在组织边缘或者切片处理过程中细胞受损线粒体比例也会很高。我的做法是先看整体分布如果大部分spot的线粒体比例集中在10%以下可以把阈值放宽到25%但如果数据质量本身不错还是建议卡在15%-20%之间。关键是不要无脑过滤要看具体切片的QC指标。第二个坑是spot的过滤。有些spot可能落在组织切片以外或者几乎没有捕获到任何转录本。Seurat里会建议过滤掉count小于某个阈值的spot但我会额外检查这些spot在组织图像上的位置。如果你用的平台带有HE染色图像建议把filtered的spot直接可视化到图像上看看它们是不是正好对应组织边缘或者切片有缺损的区域。这样做的好处是你能直观地判断过滤是否合理也能在后续分析里避开通孔区域或边缘区域的假信号。第三个坑是标准化方法的选择。这一点我必须强调SCTransform在空间转录组数据分析里的表现比传统的LogNormalize要好不少。原因有两个一是空间数据本身的捕获效率在不同spot之间存在比较明显的差异SCTransform的正则化框架能更好地处理这种技术噪声二是SCTransform可以同时回写到表示系统误差的变量比如UMI总数避免我们把技术变量当成生物学信号来聚类。我个人的建议是无论你后续用Seurat还是Scanpy预处理环节直接用SCTransform不要再用老一套的NormalizeData加FindVariableFeatures。library(Seurat) library(SeuratData) # 读取空间数据 sp_data - Load10X_Spatial(data.dir ./filtered_feature_bc_matrix/, filename filtered_feature_bc_matrix.h5, assay Spatial) # QC指标计算 sp_data[[percent.mt]] - PercentageFeatureSet(sp_data, pattern ^MT-) # 过滤 sp_data - subset(sp_data, subset nFeature_Spatial 200 nCount_Spatial 500 percent.mt 20) # 标准化 sp_data - SCTransform(sp_data, assay Spatial, verbose FALSE)预处理跑通之后你会得到一套可以用于下游分析的归一化表达矩阵和降维结果。这时候就可以正式进入下一步——根据你的生物学问题选择分析模式了。2. 为什么会有“三种分析模式”科学问题决定分析路线很多初学者最容易犯的错误是拿到数据后直接套用单细胞的分析流程做完聚类、UMAP画图之后就不知道下一步该干什么了。这背后的深层原因是没有理解空间转录组分析的“套路边界”——空间数据能回答的科学问题是有限的而且不同问题的分析路径完全不同。根据我自己的实战体会空间转录组数据在绝大多数课题里回答的是以下三类问题细胞类型在组织切片上如何分布模式一、组织微环境由哪些空间区域构成以及各区域的分子特征模式二、不同细胞或不同区域之间通过配体-受体等机制如何相互交流模式三。这三类问题对应三种分析模式核心逻辑分别是“映射已知”、“从头划分”和“关系推断”。模式一的核心思路是我手里已经有一套单细胞转录组数据可能是同一个样本的也可能来自公共数据库我已经通过无监督聚类定义了细胞类型现在要做的就是把这些细胞类型的标签“搬”到空间转录组的spot上看看每种细胞在组织里长什么样。这种分析模式的关键词是“注释”和“映射”。模式二的核心思路是我暂时不管具体的细胞类型而是把组织切片按照转录组表达模式的相似性划分成若干个区域domain或niche每个区域代表一个相对均质的状态。这种分析模式比较适合探索性研究比如你想看看肿瘤组织里有没有免疫排斥型的区域结构或者发育中的组织有哪些未被定义过的功能区域。它的关键词是“区域化”或“结构域”。模式三的核心思路是我关心的是“谁跟谁之间存在相互作用”。空间数据给了我们坐标信息两个spot如果靠得很近并且一个表达配体另一个表达受体我们就有理由推断它们之间存在潜在的通讯关系。这种分析模式的关键词是“邻域”和“通讯”。要注意的是这三种模式不是互斥的一个完整的课题通常会组合使用两到三种。比如先通过模式二识别出肿瘤边界区域再通过模式一确认边界区域的细胞构成最后用模式三分析边界区域中肿瘤细胞与免疫细胞的通讯关系。但如果你一开始没有想清楚自己的问题归属很容易在中间步骤迷失方向做出一些看起来很炫但实际解释不了生物学问题的图。在后面三个章节里我分别针对三种模式展开讲。3. 模式一基于参考注释的细胞类型映射这种模式说白了就是做标签迁移。你有一组带注释的单细胞数据作为参考要把这些标签映射到空间转录组的spot上。最常见的做法是使用Seurat的TransferData函数整体思路是把两个数据集的共享基因锚定起来然后通过典型相关分析找相似性。先说说什么时候应该用模式一。如果你的课题目标是回答“某种细胞在组织里的空间分布是怎样的”例如“CD8阳性T细胞在肿瘤组织里分布在边缘还是浸润到内部”那你基本就是走模式一。它是三种模式里最成熟、也最容易出结果的一种。但有一个前提条件容易踩坑你的参考单细胞数据最好来自于同一个样本或者至少是同一个组织的相似样本。如果参考数据和空间数据来自完全不同的组织类型比如拿胎肝的单细胞数据去注释成人肝组织的空间切片迁移效果会大打折扣甚至出现系统性的错误注释。我遇到过有同学用公共数据库里的数据直接做映射结果把空间切片上所有上皮细胞都注释成了T细胞原因就是参考数据里免疫细胞比例偏高anchoring的时候产生了偏差。实际操作时我通常会把映射结果做一个“信心度”评估不要只拿一个标签就往图上画。Seurat的TransferData会返回预测分数建议把预测分数比较低的spot标记为“未定义”或者“低可信度”而不是强行归到某一类。这一步虽然会增加一些后期处理的功夫但对结果的可靠性提升非常明显。# 参考单细胞数据 s_cell_ref已包含细胞类型注释列 cell_type DefaultAssay(sp_data) - SCT anchors - FindTransferAnchors( reference s_cell_ref, query sp_data, normalization.method SCT, reference.reduction pcaproject, recompute.residuals TRUE ) # 预测细胞类型 predictions - TransferData( anchorset anchors, refdata s_cell_ref$cell_type, weight.reduction sp_data[[pcaproject]] ) # 合并预测结果 sp_data - AddMetaData(sp_data, metadata predictions)映射完成后一定要做两件验证的事情。第一把注释结果画到组织切片上肉眼检查不同细胞类型的空间分布是否符合组织学预期。比如上皮细胞应该主要集中在组织实质区巨噬细胞可能更多分布在间质区或坏死边缘如果你看到“心肌细胞”分布在大脑皮层切片的正中间那大概率是注释出了问题。第二抽几个核心marker基因做空间可视化验证marker表达最强的区域是不是跟注释结果一致。比如你注释出了T细胞那CD3D、CD2这些基因的表达位置应该和你标注的T细胞spot区域高度重叠。还有一个比较进阶但很实用的思路如果手头没有匹配的单细胞数据可以考虑用公共数据做参考但必须先用Harmony或者CCA做跨数据集的批次校正。我实测过用不同平台10x Chromium和Smart-seq2的参考数据去映射spatial数据Harmony校正后再做TransferData比直接TransferData的准确率要高不少。不要嫌麻烦这十几行代码有时候决定了你整个分析的成败。4. 模式二空间结构域识别与微环境分型如果你拿到数据后脑子里想的是“这个组织里有哪些不同的功能区域”“肿瘤浸润边界和肿瘤核心在分子层面有什么差异”那你需要的不是简单的标签迁移而是从头对spot进行区域划分。这就是模式二要解决的问题。空间结构域识别的思路本质上是一种带空间约束的聚类。传统的非空间聚类比如Louveain或Leiden聚类只考虑表达谱的相似性不考虑位置信息这在空间数据上会产生一个问题——同一个生物学区域可能被切碎成好几块分散在切片的不同位置。为了解决这个问题空间聚类算法会在聚类过程中加入一个空间约束项让相邻的spot更倾向于被分为同一类。目前最常见的两种工具是Seurat的FindSpatialClusters实际上是对SCT变换后的表达数据进行加权的空间聚类和BayesSpace。Seurat的实现相对快适合跑通流程BayesSpace则是基于贝叶斯统计的方法引入了马尔可夫随机场来平滑空间上的聚类标签效果通常更细更好但速度慢一些对计算资源有一定要求。我自己在实践里的做法是先用Seurat跑一遍快速的空间聚类看一下大概的区域结构如果有比较强的区域边界信号再上BayesSpace做细化。两个结果对比着看选一个跟组织形态更吻合的结果作为最终的区域定义。# 使用 BayesSpace 识别空间结构域 library(BayesSpace) # 转换 Seurat 对象为 SingleCellExperiment sce - as.SingleCellExperiment(sp_data) # 设置空间邻接信息spatialLocation 是 spot 坐标 set.seed(100) sce - spatialPreprocess(sce, platform Visium, n.PCs 15) # 指定聚类数 k通常先用较小值探索 sce - spatialCluster(sce, q 7, platform Visium, nrep 5000, burnin 1000) # 拿到聚类结果 spatial_cluster - colData(sce)$spatial.clusterBayesSpace里有一个比较关键的超参数——聚类数q。这个数字怎么定我的经验是分两步先跑非空间聚类的Leiden看轮廓系数和分群稳定性得到一个大致的簇数范围然后用这个范围附近的几个q值分别跑BayesSpace比较不同q下的聚类结果在组织切片上的分布是否合理以及各个结构域之间的差异表达基因是否有生物学意义。不要盲目追求大的q值结构域越多单个结构域的样本量越小后续差异分析的统计功效会变得很弱。结构域识别完之后下一步通常是做区域间的差异表达分析找每个结构域的标记基因。这里有个建议用FindAllMarkers默认参数跑出来的结果有时会包含一些核糖体蛋白基因或热休克蛋白基因这些往往反映的是细胞活性状态的差异而不是区域身份特征。你可以通过设定features排除这些基因或者在后续分析里特别关注那些已知跟组织功能相关的基因。否则你画出来的图虽然很热闹但实际上解释不了任何生物学问题。区域分型的下游还有两个比较常见的延伸方向一是做“结构域之间的轨迹分析”比如从正常组织区域到肿瘤区域的过渡区域看看基因表达有没有渐进式的变化趋势二是结合临床信息做“结构域组成比例与预后”的分析比如某个免疫结构域的丰度是否跟患者生存相关。这些方向做出来都是不错的课题亮点。模式二最忌讳的是把聚类结果当成终点。识别出十个结构域这只是起点。关键是你能否给每个结构域一个生物学身份的命名比如“富血供区域”“免疫排斥区域”“增殖活跃区域”。这个命名的过程需要结合marker基因、通路富集分析、组织形态学验证三管齐下才能让人信服。5. 模式三配体-受体互作与空间通讯推断第三种分析模式是最能体现空间转录组“空间价值”的部分因为如果你只有单细胞数据做配体-受体互作分析时是缺乏空间约束的——两个细胞群虽然在转录组层面表达了配体和受体但在组织里可能根本相距十万八千里。空间数据天然解决了这个问题spot之间的物理距离就是最直接的约束条件。核心思路并不复杂先定义“细胞类型A和细胞类型B在空间上是相邻的”然后检查A是否表达配体L、B是否表达受体R最后通过统计检验判断这种配受体共表达是否显著地富集在相邻的spot对之间。听起来不难但实际操作中有几个关键问题。第一个问题是“什么是相邻”。10x Visium的spot之间间距是固定的spot中心距100微米一般定义相邻spot为“切比雪夫距离为1”的六个邻近spot。如果你的分析用的是更高分辨率的平台比如Stereo-seq邻域的定义可能需要调整。最稳妥的方式是先对你的数据跑一下FindSpatialNeighbors把邻接矩阵建好再结合你的平台特点设定邻域范围。第二个问题是配体受体数据库的选择。目前主流的选择有CellPhoneDB基于已知的配受体互作关系和CellChat有内置的文献支撑的通讯网络数据库。空间转录组分析里我更常用CellChat因为它提供了将空间位置信息纳入通讯概率计算的函数可以直接接受Seurat对象作为输入使用起来非常顺手。library(CellChat) # 准备 CellChat 对象meta 包含细胞类型/结构域注释列 cellchat - createCellChat(object sp_data, group.by cell_type, assay Spatial) cellchatDB - CellChatDB.human # 识别过表达配体受体基因 cellchat - subsetData(cellchat) cellchat - identifyOverExpressedGenes(cellchat) cellchat - identifyOverExpressedInteractions(cellchat) # 推断空间近邻关系中的通讯概率 cellchat - computeCommunProb(cellchat, type triMean) # 聚合通讯网络 cellchat - aggregateNet(cellchat)关于结果解读一个特别容易犯的错误是只看配体受体表达量就下结论。即使两个spot相邻一个spot表达了配体另一个spot表达了受体也不能直接说两者“一定”存在通讯。空间转录组的分辨率决定了配体受体可能来自同一个spot里的两个不同细胞而不是来自两个不同spot。换句话说共定位是通讯的前提但不等于通讯本身。要想把结论做扎实至少需要补充蛋白层面的证据比如免疫荧光共染色或者通路下游靶基因的激活证据。另一个很实用的思路是把模式三和模式二组合起来用而不是直接跳到细胞类型水平做通讯分析。我的习惯是先用BayesSpace识别出结构域然后把结构域之间的通讯关系画成一个网络图。这样做的优势是结构域级别的通讯分析稳定性远高于spot级别而且结果更便于跟组织形态学对照。你可以在图上看到“肿瘤核心区域”如何通过某些配受体通路“招募”免疫细胞到周边区域这种宏观视角恰恰是审稿人比较认可的呈现方式。如果你做的是肿瘤相关课题我强烈推荐关注“巨噬细胞-肿瘤细胞”互作中的MIF信号通路。很多空间数据里这条通路的活性都很高但CellChat默认有时候不会提示它需要手动查看。另外SEMA3家族和PTN通路也可能在空间通讯分析中反复出现如果看到不要惊讶这是很多实体瘤组织的共性信号。6. 三种模式的整合流程与常见坑位速查讲到这里三类模式的逻辑和实操重点都过了一遍。但现实中的课题很少只用到单一种模式更常见的是一个完整的数据分析流程里三种模式各司其职、环环相扣。这里我给出一个我自己比较常用的整合分析流程以及这一路跑下来总结出的常见问题和避坑心得。6.1 一个完整的组合分析流程参考第一步先跑模式二对空间转录组切片做结构域识别。这个过程中你对该组织有哪些大块区域、边界在哪、区域之间的表达差异有多大会有一个全局的印象。第二步如果手头有匹配的单细胞数据跑模式一把细胞类型标签映射到每个spot上。做完之后把细胞类型比例叠加到结构域上你就能看到每个结构域主要由哪些细胞构成。这一步相当于把“区域”和“细胞”两层信息接上了。第三步综合前两步的结果筛选出你关心的结构域或细胞类型组合做模式三的通讯分析。比如发现一个免疫排斥型结构域那就分析这个结构域里的肿瘤细胞通过哪些配受体通路抑制了T细胞浸润。最后所有的分析结果都要回到切片图像上做可视化。空间转录组的一个显著优势就是结果可以做漂亮的可视化既然有这个优势不要浪费。SpatialDimPlot、SpatialFeaturePlot这些函数画出来的图只要调好配色直接能放论文。6.2 常见问题与排查思路问题现象可能原因排查与解决办法聚类结果在空间上呈现“椒盐状”散点分布聚类时未加入空间约束普通Leiden导致的换成BayesSpace等空间聚类算法检查邻接矩阵是否构建正确标注的细胞类型在空间分布上没有逻辑参考数据与空间数据组织类型不匹配核对参考数据来源做Harmony批次校正检查预测分数低的比例高不高配体受体分析结果信号过于杂乱到处都是未设定合理的邻域范围调整邻接矩阵的构建参数尝试减少邻域范围只在特定结构域内部做通讯分析肿瘤组织切片里spot的UMI数普遍偏低组织透化时间的优化出了问题或者切片选得不好回到湿实验环节重新优化透化时间UMI低的区域在后续分析里标记为低质量区域用公共单细胞数据做映射时出现系统性错误参考数据里的细胞类型构成跟空间数据差距过大增加参考数据的样本量或者改用多数据集整合后再做映射6.3 关于平台选择和计算资源的一点经验文章开头我主要以10x Visium为例但其实现在可选的空间转录组平台越来越多包括华大的Stereo-seq、Nanostring的CosMx、Vizgen的MERSCOPE等等。这些平台在分辨率上远超Visium但分析思路的底层逻辑是相通的。如果你用的是超高分辨率平台模式一的标签迁移基本可以做到单细胞级别模式二的结构域识别会更精细模式三的邻域定义则需要参考实际的细胞密度做调整。还有一个经常被忽视的点是计算资源。空间转录组数据比单细胞数据要重得多一张切片往往就要占用十几个G的内存。跑BayesSpace或者大范围的配体受体分析时建议至少准备32G内存的服务器否则很容易跑一半卡死。如果只有个人电脑我的建议是优先跑Seurat流程BayesSpace可以用但把nrep参数调低一点比如2000。不是每个分析都值得全量跑很多时候一个切片的分析够了不用硬上多个生物学重复的全量数据——先在一个切片上把你的分析流程打磨好再批量跑其他样本。7. 写在最后从“跑通脚本”到“问对问题”空间转录组分析工具更新得很快今天流行的函数可能半年后就半废弃了但三种分析模式背后的逻辑是非常稳定的。细胞类型映射、结构域识别、空间通讯推断无论工具怎么迭代你最终要做的都是在回答这三个层面的问题。我自己踩过最大的一个坑是早期拿到数据后迫不及待地想用“最新的算法”跑一种看起来很高端的分析结果输出的图表面很美但跟生物学问题对不上浪费了一个多星期的时间才回头重新梳理思路。所以我在带新人的时候都会先让他们把要回答的科学问题写在一张纸上然后用自己的话讲清楚属于哪种分析模式最后才开始写代码。这个习惯救了不少项目的命。如果你手头正好有一批空间转录组数据我建议你也试试这个方法先定义问题再选模式最后才是抠代码。等三种模式都各跑过一遍你大概就会有一种感受——空间转录组最好玩的地方在于它逼着你去思考组织的空间结构如何决定功能这种“空间思维”一旦建立起来对你的科研思路绝对是有长远帮助的。

空间转录组三大核心分析模式：细胞映射、结构域识别与通讯推断

相关推荐

列式存储为何能实现秒级聚合？从原理到ClickHouse选型实践

Java大厂面试实战：JVM异常、并发锁与Redis踩坑全解析

2026洛阳化工产品成分分析检测排名 TOP5 CMA 资质提供含量检测、纯度检测、元素分析联系方式推荐

更多资讯

Matlab弹道仿真：从理想抛物线到空气阻力建模与数值求解

嵌入式代码重构实战：从超级大循环到事件驱动

SEO优化最常见的5个致命错误：从关键词到外链的完整避坑指南

校园卡系统：SSM框架练手项目的经典之选

SEO优化五大常见错误：从关键词堆砌到技术隐患的避坑指南

React入门实战：从概念到组件化与状态管理

今日更新

在S3C2410上移植Nucleus RTOS：启动、驱动与实时性调优实战

2026年Turnitin降重降AI二合一红黑榜：能过和别碰

2026期刊投稿AI写作工具实测：6款打分谁更值得用

本周精选

超人会飞不算本事：系统稳定依赖清晰规则与边界设计

超人VS蜘蛛侠：拆解超级IP的影响力与传播方法论

基于CNN的调制信号识别：MATLAB实现时频图分类实战

本月热门

基于YOLO的焊缝缺陷检测：从数据标注到工业部署全流程实践

分布式服务的决策记录

创业产品的首版范围

读完文章，想聊聊您的网站？