最近在带组里的研究生读文献，看到一篇关于肝癌微血管侵犯（MVI）的多组学研究，我特意把它放在组会上反复拆了两遍。这篇文献最有意思的地方，就是把 bulk 转录组、单细胞转录组、空间转录组三个层面，像拼图一样严丝合缝地拼出了一个完整的机制故事。MVI 这个临床病理指标，很多人都知道它和术后复发强相关，但你要真想搞清楚“为什么有些癌细胞特别容易钻到血管里”，单靠一张 bulk 热图远远不够。
这篇文献之所以值得当作机制挖掘的范本，是因为它没有堆测序技术，也没有把三种组学数据做成“三张独立的图”，而是让每一层数据都承担一个明确的角色，层层递进，最后落到一个可验证的生物学模型。对还在学习阶段的研究生来说，读懂这种“顶刊逻辑”，比单纯跑一套分析流程有价值得多。这篇内容我就把它的骨架拆开，顺带聊聊我们自己在实操里踩过的坑。
1. 三个数据维度，怎么搭出一套闭环逻辑
1.1 为什么研究 MVI 要同时动用三种组学
先说清楚 MVI 是什么。在肝癌病理切片上，如果显微镜下能看到癌细胞成团地出现在门静脉或肝静脉的血管腔内，就定义为微血管侵犯。它不只是一个病理描述，而是和术后早期复发、无瘤生存期缩短高度相关的独立危险因素。临床上很早就知道 MVI 阳性患者的预后明显更差，但背后的细胞和分子机制却一直没有被完整解释。
要研究这种机制，任何单一组学都有它的盲区。如果只用 bulk 转录组，你会看到 MVI 阳性组里一堆基因上调，比如间质标志物、细胞外基质相关基因，但你不知道这些信号来自哪种细胞，是肿瘤细胞自己启动了什么程序，还是周围的成纤维细胞在捣乱。如果只用单细胞转录组，你能把细胞类型分得很细，甚至找到一群可疑的促血管生成巨噬细胞，但数据一旦离体，细胞原来的空间位置就丢掉了，你没法判断这群细胞到底是不是真的待在血管旁边。如果只用空间转录组，虽然保留了组织坐标，但单点分辨率不足，不借助单细胞参考的话，很难把转录信号对应到明确的细胞亚群上。
所以最合理的做法就是让三种数据互相补位：bulk 负责“发现现象和临床关联”，单细胞负责“找出执行这一现象的细胞亚群”，空间转录组负责“证明这些细胞确实在关键位置干坏事”。这个三角逻辑一旦成立，文章的机制故事就立住了。
1.2 Bulk 转录组：临床表型与大规模筛选的“粗筛网”
bulk 转录组在这个框架里承担的是“粗筛”任务。它的优势很实在：样本获取成本低，测序数据稳定，统计分析成熟，而且可以方便地和临床病理参数做关联。在做 MVI 研究时，最常见的做法是把患者分成 MVI 阳性组和阴性组，做差异表达分析，然后对差异基因做功能富集。很多文献里都会看到，差异基因集中在细胞外基质组织、上皮间质转化、血管生成、缺氧应答这些通路，这个结果本身就在提示：MVI 的过程可能不仅仅是肿瘤细胞“跑得快”，而是整个微环境“给了它一条通道”。
但 bulk 的缺陷也必须清醒认识。肿瘤组织本身是混合体，包含肿瘤细胞、免疫细胞、成纤维细胞、内皮细胞，还有大量细胞外基质成分，所谓“差异基因”其实是所有细胞信号的平均值。你可能筛到一个高表达的基因，但不知道它究竟在哪类细胞里高表达。这正是下一步单细胞转录组要回答的问题。
我的经验是，做这类文献阅读时，不要只盯着 bulk 的显著基因列表，更要关注 GSEA 这类排序富集分析的结果，因为 MVI 相关的表型往往是整个程序性变化，而不是单个明星基因突变。通路层面的信号会给你一个更稳的方向，后续做机制模型时也能少走弯路。
1.3 单细胞转录组：细胞“户口普查”与亚群锁定
单细胞转录组把组织拆成了“单个细胞”，每个细胞都能得到一套转录组。通过聚类，你可以把组织里所有细胞分成不同的类型，比如肝细胞、恶性细胞、T 细胞、巨噬细胞、成纤维细胞、内皮细胞等。
在 MVI 研究的单细胞分析中，有几个细节特别值得学习。第一，不同样本间做细胞比例比较，比如 MVI 阳性样本里的巨噬细胞是不是比阴性样本更多，免疫抑制性 T 细胞比例是否升高。第二，把 bulk 里富集到的通路信号投射到单细胞图谱上，看这些信号主要来自哪些亚群。第三，对关键群体做进一步的亚群分析，比如把巨噬细胞分型后，可能会发现某个高表达 SPP1 或 TREM2 的亚群在 MVI 阳性肿瘤中显著富集，而这类细胞常常和促血管生成、免疫抑制相关。
这一步的产出，是把“现象”落到了“细胞”。此时研究人员心里就有了一个候选主角名单：可能是某个巨噬细胞亚群，可能是某个炎性成纤维细胞亚群，也可能是恶性细胞本身处于一个特殊状态。有了名单，才能带着问题去看空间数据的坐标。
单细胞分析最大的坑是注释。很多初学者拿到 UMAP 图就开始凭着 known markers 打标签，但遇到巨噬细胞、树突状细胞这种转录状态高度重叠的群体，很容易翻车。后面我会专门展开讲这一点。
1.4 空间转录组：细胞互作与组织生态位的“原位证据”
空间转录组与前两种数据最大的不同是保留了组织切片上的物理位置。每一个捕获点（spot）既有转录组信息，又有二维坐标。分析时可以通过去卷积算法，把单细胞里定义好的细胞亚群映射回空间表达谱中，从而知道某个细胞亚群在切片上是集中在肿瘤边界，还是弥散在癌巢内部，或者恰好聚集在血管周围。
在 MVI 研究里，空间数据最漂亮的用法，是直接回答“关键生态位”问题。你可能会看到：在 MVI 阳性样本中，某个巨噬细胞亚群和肿瘤细胞形成了特定邻域结构，两者接触的位置恰好靠近血管；在这个空间邻域里，还检测到了配体-受体对的高表达，比如肿瘤细胞高表达某一种生长因子，而巨噬细胞高表达对应受体，这就在空间上支持了旁分泌互作假说。这种证据是任何一种非空间测序都给不了的。
有了空间证据，机制模型就从“谁做了什么”升级到了“谁在什么位置和谁互动，做了什么”。这一步让审稿人有数据可查、有空间可感，文章的可信度和生物学深度都会有明显提升。
2. 从课题设计看 MVI 研究的“顶刊逻辑”
2.1 现象驱动的设问：从临床问题到科学假设
我带着组里同学拆解这篇文献时，反复强调一个观点：好的多组学机制研究，出发点绝对不是“我手头有三套测序数据”，而是“我有一个非解决不可的临床问题”。MVI 的临床重要性毋庸置疑，但它背后的生物学过程仍然模糊，所以就有了一个清晰的研究问句：MVI 阳性肿瘤的微环境究竟有什么特殊性，让癌细胞获得了血管侵袭能力？
从这个问题出发，合理的假设可能包括：特定的细胞亚群在侵袭前沿富集；肿瘤细胞与微环境细胞之间存在异常的配体-受体互作；免疫抑制生态位为肿瘤细胞逃避免疫监视并侵入血管提供了有利条件。假设不一定要全对，但一定要能被后续的多组学数据检验或修正。这套“现象—问题—假设—数据验证”的链路，正是顶刊 MVI 研究最值得模仿的地方。
说实话，不少学生拿到数据后第一反应是让生信同学“跑个差异看看”，然后对着几千个差异基因发呆。原因就是没有在第一步把假设收窄。你一旦有了要验证的机制模型，看结果时就会有方向感，不会眉毛胡子一把抓。
2.2 差异筛选与验证：如何从 Bulk 数据锁定方向
从 bulk 数据锁定方向，技术流程其实是标准化的：先做差异表达分析，然后做 GO/KEGG 富集，再做 GSEA 查看通路级趋势。MVI 阳性组里通常会看到 EMT、血管生成、ECM 受体互作、TGF-β 等信号被激活，而这些通路本身就有非常明确的生物学内涵。
我比较关注的是“表达方向是否自洽”这一点。比如 MVI 阳性组里和间质相关的基因明显上调，同时上皮分化相关的基因下调，这就是一个自洽的 EMT 趋势；如果数据里出现“EMT 基因上调、同时 Notch 通路也全面激活”，那就需要看这个 Notch 信号到底在哪些细胞里，是肿瘤细胞自己的，还是微环境细胞的，这股信号有没有功能上的协同。只看基因列表而不看通路协同，容易把文章做成“基因罗列”。
另外还有一点实操建议：尽量利用公共数据库的 bulk 队列做外部验证。MVI 阳性和阴性分组合并后要留出一部分样本做验证，或者直接用公共肝癌数据集来检验核心基因的表达差异与预后价值。这一步能显著提升差异筛选结果的可信度，也为后面机制模型的“临床意义”提供支持。
2.3 单细胞层面的细化与落点
单细胞分析在这一步的任务，是把 bulk 的粗信号“拆分归因”。具体来说，要回答两个问题：第一，核心通路或核心基因主要由哪类细胞表达；第二，这些细胞在 MVI 阳性样本中是否更富集、是否处于更极端的功能状态。
比较典型的分析路径是：先做全局的细胞注释，然后在各细胞大类里重新聚类，找到细分亚群。以巨噬细胞为例，你会看到经典的单核细胞来源亚群和组织驻留亚群，每个亚群又可以在炎症、血管生成、免疫抑制等程序上有不同倾向。接下来就可以统计每个亚群的细胞比例，比较 MVI 阳性组与阴性组的差异，同时用差异表达分析找出亚群特异性的 marker。
在机制落点上，单细胞分析还能往前推进一步：对亚群做细胞通讯分析。比如发现肿瘤细胞高表达配体 A，而巨噬细胞高表达受体 B，那么就可以初步提出一个假说——肿瘤细胞通过 A-B 互作招募或极化巨噬细胞，使其获得促血管生成表型，最终协助肿瘤侵入血管。这个假说非常具体，后续空间数据和实验验证都是围绕这个核心互作展开的。
这个阶段特别容易犯的错误是“把每一种细胞都描述一遍”。真正高分的文献只围绕和 MVI 机制相关的关键细胞展开，不相关内容的图形和描述果断舍弃，叙事才能紧凑。
2.4 空间维度的验证与机制模型
空间转录组为机制模型提供“定位证据”。怎么实现这一步？首先把单细胞注释得到的细胞亚群作为参考，用去卷积算法映射到空间 spot 上。这时你会得到整张切片上每一个 spot 的细胞类型比例。接下来可以做两类非常出效果的分析。
一类是“细胞邻域分析”。计算每个 spot 周围一定半径内其他 spot 的细胞类型组成，从而发现哪些细胞类型倾向于共定位。这种共定位不是随机的，往往反映了组织微环境中的功能单元或生态位。在 MVI 研究中，你可能会发现某种巨噬细胞亚群与恶性细胞在肿瘤边界高度共定位，且这个邻域靠近血管标记阳性区域。
另一类是“空间配受体共表达分析”。把之前单细胞通讯分析里筛选出的配体-受体对，映射到空间位置上，看它们在肿瘤血管界面是不是同时高表达。如果在同一空间生态位中确实存在配体和受体的共表达，那这个互作机制的证据链就完整了。最后再叠加公共队列验证核心基因预后价值，整个多组学故事就有了“宏观现象—细胞归因—空间验证—临床关联”的闭环。
这样一套流程走下来，文献从头到尾讲的都是同一个故事，而不是三个独立的数据报告。这也是它为什么能发到高分期刊的原因。
3. 实操拆解：跟着文献学多组学分析的重点环节
3.1 样本质量控制：多组学分析最容易翻车的一步
别觉得质控不是新东西。我见过太多研究，生信同学跑出来的图漂亮得像模板，但大家心里都清楚样本质量根本不过关，最后审稿人随便一问就问穿了。
先说 bulk 数据。建议在做差异分析之前，先检查样本的 RNA 完整性、比对率、基因检出数等指标，并对所有样本做主成分分析，看是否有明显的批次聚集效应。如果 MVI 阳性和阴性两组样本分别来自不同医院或不同测序批次，那么第一步就应该考虑用批次校正方法处理，并明确写出批次的解释比例。
单细胞数据的质控稍微复杂一些。低质量细胞往往表现为线粒体基因比例特别高、检出基因数目异常少或异常多，需要做过滤。同时还要做双细胞去除，这一步可以显著降低后期注释的错误率。空间转录组则要额外检查组织形态，切得不完整、边缘破碎都会直接影响 spot 的信号质量。
从我的经验来看，多组学的质控标准和“去批次”策略最好在做实验设计时就确定下来，而不是等全部数据拿到手再临时救火。尤其是单细胞和空间转录组，样本的冷冻切片的时机、运输条件、保存方式都会显著影响数据质量，这些如果有条件，应该提前跑一两个测试样验证流程。
3.2 不同层次的主流分析流程
如果你是为了读文献需要了解流程，不需要会敲所有代码，但至少要知道哪些环节是关键节点。如果是自己要做分析，那下面这套流程是相对稳妥的思路。
Bulk 转录组分析流程可以概括为：原始测序数据质控、比对、基因计数、差异表达分析、富集分析、通路级评分、外部队列验证。常用的开源工具这里不作具体版本推荐，但差异分析这一步，基于负二项分布的算法和基于线性模型的算法都可以，重点是评估数据分布后选择合适的模型。
单细胞转录组的核心流程更长一些：质控过滤、去双细胞、标准化、高变基因筛选、主成分分析、批次整合、聚类、UMAP 可视化、细胞注释、亚群分析、细胞通讯、轨迹分析。细胞注释建议用自动化辅助数据库打初标，再结合 canonical markers 做人工校准，绝不可以把参考数据库的结果直接当最终结论。空间转录组的流程则是：质控、归一化、单细胞参考映射（去卷积）、空间聚类、邻域分析、生态位识别、配体-受体共表达分析。
读这类文献时，你重点要看的是方法部分的几个参数：去卷积用的什么参考矩阵、邻域分析定义的半径是多少、配受体库来自哪个版本。把这几条搞明白，你就能判断作者的分析是不是自洽。尤其是邻域半径的选择，半径太小只看到点状噪声，半径太大又失去空间特异性，不同组织密度下最优参数差别很大。
3.3 如何把三个层面的结论“拧”在一起
实操中，很多人做完三套分析后发现自己手里是三个“平行宇宙”：bulk 筛到的基因和单细胞热点基因对不上，空间结果里看着也不明显。如何把结论拧在一起，其实有规律可循。
核心做法是“交叉检索”。先定主细胞群，再定主配受体对，然后回到 bulk 数据看这些配受体对应的通路是否在 MVI 阳性组激活，接着回到单细胞数据验证表达来源，最后回到空间数据验证共定位，形成一个不断校验的循环。比如单细胞注释出某个 MVI 相关巨噬细胞亚群，就把该亚群的特征基因列表拿去做基因集评分，在 bulk 队列里对每个样本计算评分，看评分能否区分 MVI 阳性和阴性。同时可以把该亚群的特征基因映射到空间转录组上去，看每个 spot 的评分分布是否集中在肿瘤边界或血管周围。这样三套数据之间就形成了横向、纵向的交叉验证。
另一个实用思路是“评分先行”。给每个细胞亚群或每条通路做一个连续变量评分，用于在不同数据层级里关联临床指标。这样做既保留了数据连续性，又提高了跨平台可比性，比单纯做“阳性/阴性”分组要稳健得多。
读文献时你可以留意，高分文章里几乎都有一个“核心特征评分图”，它往往是把单细胞亚群特征或空间生态位特征泛化到更大样本集上的桥梁。这个技巧学好了，自己的课题也能用上。
3.4 图与叙事：顶刊里多组学结果怎么排兵布阵
拆解这篇 MVI 研究的图片结构，你会发现非常明显的节奏感。第一张图通常回答“临床问题是什么”：MVI 阳性患者预后差、病理特征明确，再加上研究的总体验设计和样本信息。第二张图开始进入单细胞景观，展示主要细胞类型和关键亚群。第三张图是空间验证的核心，展示关键细胞的分布、邻域分析、生态位识别。第四张图则汇总机制模型并配以功能验证或预后关联。
这里有一个非常值得学习的细节：空间转录组的结果不会被拆得七零八落，而是被集中放在一起，形成一个“空间证据组”。这样做的好处是让审稿人和读者能够一次性捕捉到空间数据最核心的贡献。如果空间结果被散落地塞在不同章节，叙事的说服力会明显下降。
做机制图时，也不要陷入过度设计的陷阱。一张好的机制模式图只要画清楚三件事：哪些细胞、在哪里、通过什么分子互作。其他细节全部留给正文。组里不少同学画机制图时恨不得把每一条通路都画上去，结果画面越满，信息越模糊。这个毛病一定要改。
4. 常见问题与避坑心得实录
4.1 批效应处理：不同来源的样本怎么合并比较
MVI 研究的临床样本往往来自多个中心，这就引入了批次效应。我看到过有人直接把不同医院来源的样本放在一起做差异分析，结果前几个主成分反映的是医院分组而不是 MVI 分组，后面的结果全都受影响，这是很严重的错误。
处理批效应的正确姿势分两步。第一步是做严格的质控与探索性分析，画出主成分图、样本相关性热图，确认批次是否明显聚集。第二步是选择合适的批次校正策略。对于 bulk 数据，可以用比较经典的 limma 去除批次效应；对单细胞数据，更推荐先做标准流程，再用集成算法进行批次整合，然后检查校正后的细胞类型比例是否合理，避免过度校正把真实的生物学差异也抹掉了。
这里有一个经验之谈：校正后的数据一定要回到生物学层面做 sanity check。比如某一群细胞理论上只应该在某类样本里出现，校正后如果这种分布消失了，说明批次校正力度过大，应该调整参数或改用更保留生物学差异的方法。
4.2 细胞注释的“假阳性”陷阱
细胞注释是整个单细胞分析里最依赖经验、也最容易出错的环节。举个例子，巨噬细胞和树突状细胞其实共享很多标记基因，单靠一个 CD68 或 CD163 就下结论，很容易把 DC 当成巨噬细胞，或者相反。更麻烦的是肿瘤组织里的巨噬细胞还带有高度可塑性，其转录状态和单核细胞非常接近，需要结合亚群 marker、转录因子活性等综合判断。
恶意细胞和上皮细胞也会混淆。如果不对恶性细胞做拷贝数变异推断，直接用上皮 marker 去定义“肿瘤细胞”，很可能把一小群正常上皮混进来。反过来，如果恶性细胞经历过 EMT，上皮 marker 表达下调，你又可能会把它漏掉。所以注释一定要遵循“多标记投票+功能状态辅助+CNV 佐证”的流程。
我们组里之前就出过一次这种问题：一个实习生把一群高表达内皮 marker 的细胞注释成了内皮细胞，后来仔细做 CNV 推断，才发现那其实是某种特殊状态的肿瘤细胞亚群。整个下游通讯分析和空间映射都因此重跑了一遍。这个教训很昂贵：注释错一步，后面三步全白干。
4.3 空间转录组分辨率受限的补偿办法
空间转录组的每个 spot 往往包含多个细胞，所以单一 spot 的转录组是周围几种细胞的混合信号。这种情况下，如果你直接拿一个 marker 看起来空间上“没有表达”，就断定该细胞不在那里，很容易误判。我们需要通过单细胞参考去卷积，获得每个 spot 上各种细胞类型的比例，再基于比例做空间分布分析。
另一个常用的补偿办法是把空间数据和病理形态学结合起来。比如在空间切片上同时拿到 HE 染色图像，利用病理医生的标注框出血管区域、肿瘤边界、坏死区等结构，再与转录组空间特征叠加。这种方法能有效弥补 spot 分辨率不足导致的解剖位置判断困难，也便于把结果翻译成临床病理语言。
如果空间 spot 比较稀疏，建议减少使用单个 spot 层面的统计，转而对空间域（多个 spot 的聚类结果）进行统计。空间域比单个 spot 更稳定，也不容易被背景信号干扰。很多初学者不懂这个差别，做了单 spot 层面的检验，最后 q 值全都很大，还以为自己的数据不行，其实只是分析单元选错了。
4.4 审稿人视角：多组学文章的“验证闭环”
多组学探索性文章如果没有任何湿实验验证，很容易被审稿人问住。单纯靠公共数据库做生存分析也不够。顶刊里的 MVI 多组学研究，通常会搭配多重免疫荧光染色、体外细胞共培养实验或类器官实验，验证关键配体-受体互作的功能方向。
这种验证不一定要非常深入。比如，可以用多重免疫荧光验证关键细胞亚群的空间共定位，证明单细胞和空间转录组的结论在蛋白水平上也成立；用共培养实验证明某细胞群的条件培养基能够促进内皮细胞的血管形成能力；再用通路抑制剂或者中和抗体回补，验证机制因果性。这样一套相对轻量的功能实验，就能把机制研究的证据链补完整。
读文献的时候，我会格外留意关键配体受体对的验证思路。看作者是在空间上证明了共表达，还是在功能实验里用阻断抗体直接证明了因果关系。前者给出“位置线索”，后者给出“因果证据”，两者结合才能把机制讲到及格线以上。
我自己的体会是，读完一篇好的多组学研究，收获最大的并不是某个具体的分子靶点，而是那条层层递进的逻辑链。Bulk、单细胞、空间转录组的组合不是数据军备竞赛，而是对同一个生物学问题的多角度回答，每一层都要在前一层的地基上继续推进。这篇文章之所以能成为范本，正在于它把每个数据维度都用在刀刃上，没有脱离临床问题去做单细胞注释大杂烩。
如果你也正在设计类似的多组学课题，我的建议很简单：先把要回答的科学问题写成一句话，再拆成三个小问题，分别安排给 bulk、单细胞和空间数据，让三套数据各司其职。文献里那种环环相扣的感觉，本质上就是提前设计出来的。后面等我整理完空间组学的分析复盘，再来跟你们细聊具体参数怎么调。
