空间转录组+scRNA+bulk RNA-seq：垂体瘤免疫微环境多组学解析

我最近跟某课题组一起梳理了一个垂体瘤多组学项目核心就是标题里这句空间转录组 scRNA bulk RNA-Seq组合起来看免疫细胞怎么影响肿瘤生长。做完这轮分析我最大的感受是垂体瘤这病太容易被良性两个字带偏了。很多病例在病理上确实长得慢、边界清但总有那么一部分就是往海绵窦里钻、手术切不干净、三五年后又冒出来。单靠激素分型和增殖指数Ki-67其实说明不了太多。而免疫微环境这个角度恰好能补上这块解释力——只是单独用任何一种组学都讲不完整。这篇文章主要写给三类人一是准备设计肿瘤微环境多组学项目的同学二是做垂体瘤方向但还没碰过空间转录组的研究者三是已经跑通了scRNA流程、想往上叠加空间和bulk验证的生信新人。我会把这几种技术各自的价值、怎么串联成一条证据链、具体跑分析时容易踩的坑都拆开讲。1. 为什么偏偏是这三种组学组合1.1 单看bulk只得到一个平均口味的结果bulk RNA-Seq干的事情是把一整块垂体瘤组织磨碎、提RNA、测序最后得到的是整个组织里所有细胞表达信号的混合平均。这个平均会掩盖掉非常关键的结构信息。你可以理解成把一锅汤搅匀之后尝味道能知道大概咸了还是淡了、放没放姜但你是分不清里面哪块肉是哪个部位的也看不出哪块肉附近围着哪些配料。在垂体瘤里尤其明显。肿瘤组织不是均一的团块里面有增殖状态的肿瘤细胞、散在的基质细胞、浸润的T淋巴细胞、巨噬细胞还有不少血管内皮和周细胞。如果某类免疫细胞只在组织的某个角落富集或者只占整体细胞数的很小比例bulk信号很容易把它淹没在主要细胞类型的表达背景里。你可能会看到巨噬细胞标志物整体上调但根本不知道是哪个亚群在贡献这个信号也不知道它们待在肿瘤的什么位置。不过bulk也有不可替代的优点样本处理简单、便宜、可重复性高适合做大样本队列拿来做临床关联——比如某个基因高表达的患者是不是更容易复发。所以我把bulk定位成横向铺开的工具给单细胞和空间里发现的候选分子做群体层面的验证。1.2 scRNA-seq看清了细胞类型却弄丢了空间单细胞RNA测序解决的是异质性问题。组织被解离成单细胞悬液后逐个细胞捕获转录组你就能非常直白地看到垂体瘤里到底有哪些细胞类型GH阳性的肿瘤细胞、PRL阳性的肿瘤细胞、CD3阳性的T细胞亚群、CD68阳性的巨噬细胞、浆样树突细胞以及一小部分B细胞和NK细胞。更进一步同一大类里还能分出功能状态——比如巨噬细胞可以按SPP1、CCL18、C1Q等标志基因分出不同的亚群T细胞里可以分出初始态、效应态和耗竭态。但scRNA有一个天生缺陷解离动作把组织彻底拆散了细胞之间的空间邻域关系全部丢失。你拿到一张聚类图知道某个巨噬细胞亚群和某种耗竭T细胞都存在可它们到底是不是挨在一起的是不是真的有直接接触——这个问题在单细胞数据里是回答不了的。而恰恰这个距离非常关键因为免疫细胞对肿瘤的作用很多是通过短距离接触或旁分泌实现的。我在这个项目里很早就意识到单细胞数据适合用来提名候选人找出差异显著的免疫细胞亚群、锁定亚群特异性的标志基因和配体受体对但要把这些候选人和肿瘤生长的关系坐实必须有空间信息。1.3 空间转录组把表达信号贴回组织切片空间转录组是我觉得这整个项目里最有意思的一块。它的原理是在组织切片上直接捕获mRNA然后测序让每个表达值都带一个切片上的坐标位置。Visium这一类捕获平台spot直径大概是55微米量级每个spot里面通常会有几个到十几个细胞虽然不是真正单细胞分辨率但已经足够把细胞类型的大致分布还原出来。也有更高分辨率的平台能做到近似单细胞级别但成本和处理门槛都会翻倍。对垂体瘤而言空间转录组能回答最核心的问题某类免疫细胞是不是富集在肿瘤侵袭前沿是不是跟肿瘤细胞形成了聚集生态位还是只在肿瘤间质区域里均匀散落这些信息一旦拿到bulk和scRNA里那些悬浮的发现才能落到组织学证据上。所以我的结论是这三种组学不是简单的重复验证而是互补的三个维度bulk负责广度和临床相关性scRNA负责细胞种类和状态的分辨空间转录组负责把细胞放回真实的组织坐标系里看关系。三者组合就能形成发现亚群—验证关联—定位机制的完整证据链。2. 垂体瘤研究设计思路2.1 样本收集新鲜组织是一切的前提这类项目里样本是整个流程里最容易翻车也最没得补救的环节。垂体瘤手术标本拿到手之后一条非常紧的时间线就开始了。scRNA需要新鲜组织在离体后尽量短的时间内完成解离通常我会控制在30分钟甚至更短时间拖得越久活细胞比例越低神经元和内分泌细胞这类脆弱的细胞最容易先死掉。空间转录组则需要把组织用OCT包埋后快速冷冻最好切成10微米左右的连续切片一部分做RNA捕获一部分做HE染色用于病理注释。bulk的部分可以最简单从同一块标本上取一小块组织提RNA就行但要注意如果标本不够大优先保scRNA和空间因为这两者对样本新鲜度最敏感。分组设计上不能拍脑袋。垂体瘤按侵袭行为分侵袭性和非侵袭性按复发情况可以分复发组和初发组按体积可以分微腺瘤和大腺瘤按激素分泌可以分GH瘤、PRL瘤、ACTH瘤、无功能瘤。关键是临床指标必须在收集样本的时候就定好标准。比如侵袭性这件事到底以影像上海绵窦侵犯为准还是以术中所见为准还是以Knosp分级为准这个不统一后面所有组间比较都会出现不一致的结果。我建议用影像加术中所见综合评估至少每组的样本量不要低于3例条件允许做到5到8例一组这样后续差异分析才有基本的统计效力。2.2 免疫细胞注释的底表准备注释是单细胞分析里最考验经验的一步尤其垂体瘤有它自己的特殊性。垂体前叶本身有大量内分泌细胞GH、PRL、ACTH、TSH、LH/FSH阳性的细胞这些细胞的转录组特征和肿瘤细胞非常接近但它们是正常组织成分要跟真正的肿瘤上皮仔细区分。常用的做法是先看PIT1、TPIT、SF1这类谱系转录因子的表达模式再结合激素基因GH1、PRL、POMC等把内分泌谱系和免疫谱系分开。免疫细胞的marker组合我习惯用一个精简的底表T细胞看CD3D、CD3E、CD2CD8阳性T细胞看CD8A、CD8BCD4阳性T细胞看IL7R、FOXP3Treg的核内标记RNA水平也靠谱NK细胞看NKG7、KLRD1、FCGR3AB细胞看MS4A1、CD79A巨噬细胞看LYZ、CD68、C1QA/C1QB树突细胞看CLEC9A、LILRA4中性粒细胞看S100A8、S100A9。但注意中性粒细胞在解离过程中很容易丢失scRNA里看到的中性粒细胞比例通常偏低不要因此就下结论说垂体瘤里没有中性粒细胞浸润。我踩过一个坑最早做注释时只看单个marker比如看到CD68阳性就把整个簇标成巨噬细胞结果那个簇里其实混了一部分内皮细胞导致后面差异基因全乱了。后来改成多marker交叉验证同时看了CD68、C1QA、LYZ、CD163、MRC1这些才算把巨噬和内皮彻底分开。这一步千万别图快。2.3 三种数据怎么串成一条逻辑链我在这类多组学项目里用的分析策略是分三步走。第一步scRNA先上把所有样本的免疫细胞单独提出来重新聚类找到侵袭组和非侵袭组之间比例有明显差异的亚群。第二步把scRNA里找到的关键亚群marker和关键配受体放到更大的bulk队列里验证——看这些基因的高表达是否确实与肿瘤体积、复发等临床指标关联。第三步利用空间转录组把前两步得到的明星分子和细胞类型定位到具体组织区域看它们的空间富集位置是否与侵袭前沿等结构吻合。这样做的好处是每一步的产出都能指导下一步不会出现三种组学结果互相打架、最后不知道信谁的情况。还有一个更省力的地方空间反卷积需要一套高质量的细胞类型参考表达谱而scRNA正好可以提供这个参考不需要额外去公共数据库里找或重新注释。技术维度bulk RNA-SeqscRNA-seq空间转录组分辨率组织块整体单细胞空间spot/单细胞空间位置无无有核心产出差异表达、预后关联细胞亚群、状态、通讯细胞分布、邻域富集、生态位优势样本宽容、成本低、适合队列异质性精细、亚群发现能力强还原组织结构、定位证据局限掩盖异质性丢空间、有解离偏好深度有限、免疫信号偏弱3. 核心分析流程从原始数据到空间证据3.1 scRNA-seq质控、聚类与亚群锁定单细胞分析的第一个硬门槛是质控。垂体瘤组织的解离产物里经常夹杂细胞碎片和双细胞我一般设置基因数在500到2500左右、UMI数大于1000、线粒体基因比例低于20%。有些垂体瘤样本如果解离时间太长线粒体比例会飙到30%以上这类样本建议直接剔除或重做不要强行保留。质控之后做标准化和降维用SCTransform或者常规的LogNormalize都行重点是要检查前30到50个主成分是否能覆盖足够多的方差。聚类这一步我先分大类再分小类。第一轮聚类分辨出上皮/内分泌细胞、T/NK、髓系、B/浆细胞、内皮和基质这几个大类第二轮把免疫细胞子集抽出来重新聚类。巨噬细胞亚群一般用0.5到1.0的分辨率来切切完之后逐个看亚群之间的差异基因确认每个簇都有独立的生物学解释而不是同一个细胞类型被随机切开。亚群锁定之后还要做一步嫌疑排除看目标亚群的平均增殖评分是不是偏高。如果某个巨噬细胞亚群的MKI67、PCNA高表达说明它可能是正在分裂的细胞不一定是功能表型稳定的亚群这种情况需要小心处理往往要把它单独拎出来看是不是由细胞周期主导了聚类。3.2 空间转录组从表达矩阵到空间证据空间部分我用的流程是比对得到每个spot的基因表达矩阵后首先做基础质控。每个spot的UMI中位数要在合理范围内线粒体基因比例不能过高组织区域的HE染色和捕获区域能对应上。下一步是把病理医生的人工注释和转录组数据对齐比如标出肿瘤主体区、侵袭边缘、间质区、坏死区后续分析都在这几张注释上做。细胞类型定位这一步我优先用空间反卷积工具比如Cell2location或SpatialDecon。输入是scRNA得到的细胞类型表达均值输出是每个spot上各细胞类型的估计比例。但这里要特别提醒空间转录组的每个spot深度有限免疫细胞相关的信号往往较弱反卷积结果只能看趋势和相对富集不能当成绝对定量。我在项目里对巨噬细胞亚群的重建效果一般对T细胞这类表达丰度偏低的分辨会更差一些。邻域富集分析是空间数据里我最看重的产出。先给每个spot计算目标细胞类型的富集分数然后看这些高分区是不是在侵袭前沿附近成簇出现。具体做法可以用置换检验把细胞类型标签打乱1000次看实际观察到的空间共定位分数是不是显著高于随机分布。如果某种巨噬细胞亚群在肿瘤侵袭边缘spot中的富集显著高于内部区域这个证据强度就上来了。3.3 bulk RNA-Seq差异表达与免疫浸润估算bulk这边如果样本量允许我建议自己做测序如果资源紧张也可以用公共数据库里已有的垂体瘤转录组队列做补充验证。差异表达分析我用DESeq2或limma这套流程取FDR小于0.05、|log2FC|大于1的基因作为显著差异基因然后跑GO、KEGG和HALLMARK富集重点关注炎症反应、补体系统、TNF-α/NF-κB、上皮间质转化这几类通路的富集情况。免疫浸润估算这块我习惯同时跑几个算法对比比如ssGSEA、CIBERSORTx、MCPcounter、xCell。每个算法的原理和参考矩阵不一样结果会存在差异所以不能只看某一个工具的输出。我的经验是如果多个算法都指向同一类免疫细胞在侵袭组中显著富集这个结论才比较稳如果出现矛盾先看是不是参考矩阵里没有垂体相关细胞类型导致的偏差。另外把scRNA里锁定亚群的marker基因做成一个基因集在bulk样本里用ssGSEA打分这个操作很实用。它能直接回答一个问题侵袭性垂体瘤组织里是不是整体上更偏向那个巨噬细胞亚群的表型如果分数差异显著就算衔接上了单细胞和群体转录组的证据。3.4 细胞通讯配体受体的空间证据细胞通讯分析是这个项目的临门一脚。我用CellChat在scRNA数据上做免疫细胞和肿瘤细胞之间的配体-受体推断再结合空间共定位结果来判断这些通讯可能性有没有空间支撑。关键在两点。第一点配体-受体对的选择要有方向性。比如某个巨噬亚群高表达CCL18和SPP1肿瘤细胞高表达CCR8和CD44这才构成完整的信号通路如果受体在肿瘤细胞里根本不表达那这个通路的假设就站不住。第二点通讯结果一定要回到空间里去看。如果CellChat提示某个配受体对很强但空间数据里这个巨噬亚群和肿瘤细胞各在各的区域、根本不聚在一起那最有可能是远端旁分泌或者干脆是数据噪音。反过来如果空间上两者确实紧挨着这个证据链就相当完整了。我不推荐只跑一遍CellChat输出一堆气泡图就收工。实际操作里要反复交叉先缩小到几个最值得关注的配受体对然后回到空间数据里为每一对做共富集检验最后在组织切片图上画出来。这样出来的结果评审和合作病理医生都更容易信服。4. 实操中常见的坑和排查思路4.1 免疫细胞比例过低如何在淹没之前抓住信号垂体瘤组织里肿瘤上皮细胞往往占绝对优势免疫细胞比例按经验看通常只有个位数到百分之十几。如果直接全组织解离做scRNA很可能测了几千个细胞真正的免疫细胞只有几百个亚群分析根本没法做。我的解决办法有两个方向。一是如果富集条件允许解离后做CD45阳性分选把免疫细胞富集出来再做单细胞测序这样免疫细胞亚群的分辨率会大幅提升。二是如果不想加实验步骤就把每个样本的测序细胞数提到尽量高同时在数据过滤时对免疫细胞相关marker更宽容一些。就算这样稀有细胞类型也还是容易在聚类时被并入其他大簇这时候需要把表达目标marker的细胞单独抽出来做一个冷门亚群存档宁可单独聚类也不要随便牺牲掉。4.2 批次效应整合不是万能的多患者样本合并后批次效应几乎是必然出现的问题。不同患者、不同测序批次、不同样本处理时间都会在UMAP上形成明显的患者在各自小岛的结构。这时候用Harmony或者Seurat的整合流程确实能抹掉一部分批次差异但我要提醒一句整合算法也会把真实生物学差异一并抹平。我踩过这个坑。有一轮整合之后侵袭组和非侵袭组之间的巨噬细胞亚群差异变得非常不显著我当时还以为是设计出了问题后来回到单个样本里单独看发现差异其实是存在的只不过整合过猛给平均掉了。从那以后我养成了一个习惯整合前先按样本各自聚类看一遍结构整合后再分别跑两遍差异分析对比整合前后差异基因的方向和显著性。如果整合后某些已知的临床相关差异消失了宁可手工校正批次也不要一味依赖算法。4.3 多组学结论不一致多数时候不是bug项目做到中后段最容易让人自闭的是bulk、scRNA、空间三边的结果对不上。比如某个基因在scRNA的特定亚群里明显高表达但在bulk的整体差异分析里不显著或者某个免疫细胞在空间里明明富集在侵袭前沿但bulk的CIBERSORTx结果里该细胞类型组间差异却不明显。这通常不是分析错误而是三组学数据结构不同造成的理所当然的不一致scRNA看的是亚群细分bulk把亚群信号摊到了整个组织里空间只反映局部区域。处理方式是重新定义验证粒度——bulk里不显著就换成看这个亚群的marker基因集整体打分CIBERSORTx不显著就回到反卷积结果里看这个细胞类型的平均比例趋势是否与空间富集一致。只要方向大体一致就不要轻易否定数据。4.4 空间切片质检别让数据死在第一步空间转录组实验环节里组织切片质量比后面任何分析参数都重要。我见过最可惜的情况是样品解离和scRNA都做得很好但OCT包埋的时候留了气泡切片一刮就碎或者组织在运输过程中部分解冻RNA降解严重捕获后总UMI数低得没法用。切片做完第一件事是看镜下形态组织有没有完整细胞核结构、有没有大片坏死、有没有折叠或者刮痕。第二步看RIN值空间平台通常要求RNA完整性在7以上。第三步在正式上机前先跑一次小规模文库或至少测试染色确认组织形态和RNA信号能对应上。钙化或出血区域多的垂体瘤样本建议在注释时显式标记这些区域下游分析里排除否则会出现假的低表达区域。5. 结果解读与机制假设的推进5.1 从亚群表型到空间排布形成故事分析做完之后怎么把散落的图表串成一个有说服力的故事这步往往决定项目最终能发到哪个档次。我以模拟项目X为例完整走一遍推导链条。假设scRNA里发现在侵袭性垂体瘤样本中富集了一类SPP1和CCL18高表达的巨噬细胞亚群同时T细胞亚群中出现耗竭特征基因如TIGIT、LAG3上调。这个结果本身只能说明有这种巨噬亚群和T细胞耗竭两者之间没有直接联系。接下来用空间转录组看这巨噬亚群在侵袭前沿的富集分数显著高于肿瘤内部而在侵袭前沿周围T细胞的密度并没有同步上升反而在空间上有负相关趋势。这时再回到bulkSPP1和CCL18高表达的患者组T细胞活化相关基因集评分偏低肿瘤增殖评分偏高。三条证据合在一起就能形成一个值得进一步验证的机制假说这类巨噬细胞可能通过分泌趋化因子和抑制性信号在肿瘤侵袭前沿构建了一个局部免疫抑制生态位。我不建议在空间数据里强行叠加一个因果帽子。空间共定位只能说明相关真正机制要靠实验去闭环。5.2 阴性、矛盾和看似不预期的结果怎么处理项目里一定会出现阴性结果。比如某个常规的好预后免疫指标——细胞毒性T细胞浸润评分——在垂体瘤里跟复发没有显著关联。这种阴性结果其实有它的价值。如果同时观察到T细胞是耗竭状态那说明问题不在于T细胞有没有来而是来了但被抑制了。这种情况下阴性比阳性更能支撑微环境抑制假说。再比如B细胞和浆细胞在某些样本里数量很低低到没法做亚群分析那就不要硬做直接在这类样本里删掉B细胞分群或者只在有足够数量的亚组里进行分析。强行分析稀有类型只会得到一堆不可靠的小簇。5.3 实验验证方案推荐多组学分析给出的永远是需要实验验证的假说而不是结论。我按成本和可操作性排序推荐几个验证手段。第一个是多重免疫荧光或多色免疫组化直接在垂体瘤组织切片上共染目标巨噬细胞marker、肿瘤marker和T细胞marker验证侵袭前沿是否真的存在这些细胞的空间聚集。这个是性价比最高的方案样本可以复用。第二个是条件培养基或共培养实验。用巨噬细胞的条件培养基处理垂体瘤细胞系或者原代肿瘤细胞检测增殖、迁移、侵袭能力的变化还可以加中和抗体来验证某个关键细胞因子的作用。第三个是用已有的垂体瘤bulk队列做临床关联。如果样本量够直接对关键巨噬亚群marker高表达组和低表达组做无复发生存曲线这一步往往能显著提升说服力。如果空间结果提示特定的配受体关系能加做一个小干扰RNA或者抑制剂实验验证受体在垂体瘤细胞上的功能那就更完整了。6. 几点实操心得6.1 项目流程管理组学项目最容易被忽略的坑这种多组学项目样本信息复杂、测序批次多、分析环节长我在做的时候最深刻的体会是流程管理决定了项目能不能说得清。每一个样本从离体时间、OCT包埋编号、切片的HE图、RNA提取浓度到测序文件名都要贯穿着同一个样本编号体系。如果前期记录不规范后期碰到一个异常数据你根本没办法倒推是实验问题还是分析问题。我建议从一开始就按批次建好样本清单把临床信息脱敏后单独存一份映射表分析代码里统一用编号引用样本不硬编码路径。所有软件的版本、参数和运行时间都记在日志文件里。这件事看起来很简单但在项目收尾复盘的时候能省出整整一周的返工时间。6.2 后续可以怎么扩展这个项目框架本身有很强的可扩展性。如果样本条件允许可以在单细胞层面加做TCR/BCR组库分析看肿瘤浸润T细胞的克隆扩增情况能更好地区分旁观者T细胞和真正靶向肿瘤抗原的T细胞配合空间数据能进一步确认克隆扩增细胞的位置。空间维度上如果条件允许可以升级到真正单细胞分辨率的空间平台这样对侵袭前沿微环境的解析会更细或者用多重免疫荧光做高内涵定量把空间转录组里的spot结果在蛋白层面做大范围验证。如果往后想往转化方向走还可以把找到的候选配受体对整理成药物靶点谱结合公共药物数据库做优先级排序为后续的功能实验选择提供参考。当然这只是一种可能的延伸前提是前面的分析结论足够扎实。这个项目整个过程做下来我个人最大的体会是多组学不是越多越好而是每一组都要完成别的组替代不了的任务。bulk做广度scRNA做深度空间做位置三者各司其职才能把免疫细胞影响垂体瘤生长这顶帽子稳稳地戴在数据头上。如果让我给一个建议那就是无论分析做了多少层永远留一部分精力回到组织切片上把数字翻译成组织学里看得见、摸得着的事实。