400-0755-639
联系微信客服
学术文献
三叶兴科为医学诊断、治疗以及医药行业的发展注入了新的活力,我们是您研发工作中值得信赖的伙伴!
首页 学术文献 Nature Methods | 从空间单组学到空间多组学:NicheTrans如何补齐缺失的分子层

Nature Methods | 从空间单组学到空间多组学:NicheTrans如何补齐缺失的分子层

来源:生物探索 发表时间:2026-07-15 阅读:0

同样的转录本,为什么未必对应同样的蛋白?

传统跨组学预测常把细胞视为独立样本:输入基因表达,输出蛋白质或代谢物。但组织中的细胞并不孤立。即使两个细胞的转录表达完全相同,只要位置、邻近细胞以及与血管、病灶或斑块的距离不同,它们的蛋白和代谢状态就可能不同。

研究团队在EEL FISH、卵巢癌MERSCOPE和乳腺癌Xenium三类数据中,共分析了547,097个细胞,发现“表达谱完全相同”的细胞可以分布在组织的不同位置。若模型只读取细胞自身的转录本,就会给它们输出相同答案;但在组织中,“细胞是谁”还取决于“它在哪里”和“周围是谁”。

另一种办法是测量相邻切片,再将两张切片对齐。然而,研究展示的乳腺癌连续切片间隔仅5微米,宏观结构仍然相似,微观层面却已经出现细胞数量、密度和位置差异。邻片对齐能恢复大体结构,却难以保证细胞一一对应。

这意味着,空间跨组学预测不能只把坐标当作一个附加变量。它需要理解细胞所处的局部环境。

NicheTrans的关键,不是加坐标,而是重建局部生态位

NicheTrans中的“niche”指细胞生态位。模型把中心细胞及其邻近细胞组成一个局部单元:低分辨率spot数据使用8个邻居,单细胞分辨率数据使用12个邻居,近邻与次近邻被编码为不同的空间层级。

模型随后整合基因表达、空间关系、可选的细胞类型信息,以及可选的苏木精-伊红染色(hematoxylin and eosin,H&E)图像特征,通过Transformer学习中心细胞与周围细胞之间的依赖关系,最后输出蛋白质或代谢物。

它回答的不再是“这个转录组对应什么蛋白”,而是“这个细胞处在这样的局部环境中,更可能呈现什么分子状态”。

这种变化看似只是多看了几个邻居,实际上却改变了跨组学预测的基本单位:从单个细胞,变成了细胞与其微环境构成的局部系统。

帕金森病脑切片:基因表达能否“翻译”出代谢地图?

研究首先分析了帕金森病(Parkinson’s disease,PD)小鼠脑组织。三张连续切片同时获得空间转录组、空间代谢组和H&E图像,其中两张用于训练,一张用于测试。测试切片包含1,242个完整侧spot1,433个病损侧spot。研究将NicheTrans与六种既有跨组学翻译方法进行了比较。

31种目标代谢物的整体评估中,NicheTrans以及加入组织学图像的NicheTrans*,在Pearson相关系数和Spearman相关系数上均处于领先位置。对于三种空间分布不同的代谢物,其他方法常出现明显偏移,而NicheTrans得到的空间分布更接近实测结果;加入H&E图像后,不同脑区之间的代谢梯度进一步接近实测图谱。

更关键的问题是:这些预测能否改善组织分区?

仅使用转录组时,模型难以区分完整侧和病损侧尾状壳核;仅使用代谢组时,又难以恢复多个功能脑区。将NicheTrans*生成的代谢信息与实测转录组联合后,空间分区的调整兰德指数(adjusted Rand index,ARI)达到0.73,标准化互信息(normalized mutual information,NMI)达到0.75,明显优于其他翻译方法。

因此,跨组学翻译的价值不能只看某个代谢物预测得“像不像”,还要看补出的分子层能否改变对组织结构和病理边界的认识。

模型给出的不只是预测,还有可检验的机制线索

研究人员利用积分梯度(integrated gradients)追踪哪些输入基因对目标代谢物的预测贡献最大。

针对与PD相关的代谢物m/z 674.28592,模型识别出Tac1、Pcp4、Penk、Gpr88、Pdyn和Hpca等高贡献基因。完整侧与病损侧之间,Tac1的表达差异P值为7.2×10⁻²⁰,Pcp4为1.5×10⁻¹⁸,提示这些差异并非轻微波动。

这些基因进一步富集于信号转导、阿片信号、G蛋白偶联受体(G protein-coupled receptor,GPCR)信号、5-羟色胺递质释放和多巴胺递质释放等通路。

这使模型的作用不再局限于复原一张代谢地图。它还可以反向提出:哪些基因程序可能与特定代谢变化相关?

不过,模型归因只表示预测过程依赖这些基因,并不等同于证明这些基因直接导致了目标代谢物变化。它提供的是值得实验验证的机制入口,而不是因果结论。

阿尔茨海默病:最难的不是识别阴性,而是找出稀少阳性

第二个关键场景是阿尔茨海默病(Alzheimer’s disease,AD)小鼠脑组织。

研究使用两张AD模型切片和一张同龄野生型切片,从空间转录组预测淀粉样β蛋白(amyloid-β,Aβ)和磷酸化tau蛋白(phosphorylated tau,p-tau)。这两类蛋白是AD研究中的重要病理标志,但其阳性信号远少于阴性信号,构成典型的类别不平衡问题。

在这种数据中,一个模型即使把绝大多数位置都判断为阴性,也可能获得很高的准确率和特异性,却漏掉大量具有生物学意义的病理区域。因此,F1分数和灵敏度比表面上的总体准确率更值得关注。

NicheTrans对Aβ的受试者工作特征曲线下面积(area under the ROC curve,AUROC)为0.931,对p-tau为0.841。与表现最好的单细胞翻译方法相比,Aβ和p-tau的F1分数分别提高0.261和0.206,即26.1和20.6个百分点

但疾病切片上的高灵敏度还不够。一个模型也可能只是倾向于“多报阳性”。

野生型切片因此构成了更严格的反向检查。该切片共有8,035个细胞,实测仅有2个Aβ阳性和61个p-tau阳性信号。NicheTrans*预测为2个Aβ阳性和52个p-tau阳性,接近实测结果。

相比之下,sciPENN预测了13个Aβ阳性和421个p-tau阳性,cTP-net虽然没有预测出Aβ阳性,却预测了197个p-tau阳性

一个模型在疾病样本中“看见更多病理”并不一定更好。能否在正常样本中不制造大量假阳性,是稳健性的另一半。

当预测Aβ成为空间地标,病理微环境也被重新拼出

研究人员随后把预测出的Aβ斑块作为空间地标,分析其周围的细胞组成。

置换检验显示,Micro3小胶质细胞亚型、Astro3星形胶质细胞亚型和Oligo1少突胶质细胞亚型均显著富集于Aβ周围,而且使用预测Aβ与使用实测Aβ得到的结论一致。

这一分析重建了Aβ斑块周围的“核心—外壳”结构:Micro3更靠近斑块,Astro3和Oligo1更多分布在外层。

模型还识别出Cst7、Ctsb、Trem2和Apoe等Micro3标志基因。高贡献基因富集于先天免疫系统、免疫系统和Toll样受体级联通路;在疾病富集分析排名前七的条目中,有五个直接涉及AD

这一步很重要。因为计算翻译不再只是补齐一个缺失的分子层,而是开始帮助重建细胞类型、病理地标和分子程序之间的空间关系。

如果预测出的蛋白分布只能在图像上“看起来相似”,它的研究价值仍然有限;如果这些预测还能复原已知的病理微环境,并提出新的空间关联假设,它才可能成为后续机制研究的工具。

它没有让空间多组学实验变得多余

NicheTrans最容易被误读为:“以后只测转录组,就能得到蛋白组和代谢组。”

研究并不支持这样的外推。

首先,模型仍然需要配对的空间多组学数据进行训练。没有高质量的实测多组学数据,就无法建立可靠的翻译关系。

其次,PD分析主要基于同一只小鼠脑的三张连续切片;AD部分也只有两张AD模型切片和一张野生型切片。这些结果证明了方法的可行性,却不能自动保证模型在不同物种、不同疾病、不同实验室和不同检测平台之间都能保持同等性能。

再次,Aβ和p-tau本身位于细胞内外不同空间。研究需要把这些蛋白信号分配给邻近细胞:p-tau按15像素范围内最近的3个细胞标记Aβ按50像素范围内最近的3个细胞标记。这些预处理规则会直接影响训练标签。

因此,NicheTrans生成的是在训练数据和模型假设条件下得到的推断组学不是新的直接测量。它适合用于候选发现、空间分区、机制线索筛选和实验优先级排序,但关键结论仍需独立实验验证。

更值得关注的变化:研究单位正在从“细胞”转向“细胞加环境”

这项研究更重要的启示,是重新定义了跨组学预测的输入边界。

组织中的分子状态并不是只由细胞内部转录本决定。它还受到邻近细胞、组织结构、局部病理状态和空间位置的共同影响。

当空间信息被写入模型后,相同的转录表达不再必然导向相同的蛋白或代谢预测。看似只是增加了“邻居”,实质上却把一对一映射改成了情境化推断。

未来更值得关注的,也许不是“模型能否替代一次昂贵的多组学测量”,而是:在有限预算下,哪些分子层必须直接检测?哪些可以通过计算推断?哪些预测最值得优先进行实验验证?

NicheTrans提供了一种新的研究分工:用少量高质量空间多组学数据建立参照,再让更大规模的空间单组学样本获得可检验的多层分子假设。

它没有消除实验,而是在重新安排实验应该把资源用在哪里。

 
 

参考文献

 
Wang Z, Zou Q, Lin S, Li S, Cui Y, Zhang D, Han C, Li Y, Li J, Zhao Y, Gao R, Song J, Yuan Z. NicheTrans: spatial-aware cross-omics translation. Nat Methods. 2026 Jul 9. doi: 10.1038/s41592-026-03153-3. Epub ahead of print. PMID: 42426404.