6.2/Q2,郑州大学生信力作|NMF+SCENIC 多算法拆解神经免疫轴,挖掘乳腺癌转移专属预后标志物 COL1A2

医药数据库是什么？今天，我们为你揭开「医药数据宝库」的神秘面纱！想知道如何免费使用这个「医药神器」？往下看，3分钟教你用对数据库，健康少走弯路！”

还在发愁乳腺癌生信课题没创新点？郑州大学这篇发在《Breast Cancer Research》高分文章直接给完整可复制模板！作者直接薅 GEO、UCSC Xena 三大免费公共数据集，瞄准当下超级火热的神经免疫热点分子 GABA，一套组合 AI 算法拉满：NMF 细分免疫亚群、Monocle 画分化轨迹、SCENIC 抓转录因子、CellChat 拆细胞通讯，单细胞 + 批量双组学联动深挖乳腺癌转移微环境隐藏机制，纯生信直接拿下权威肿瘤期刊，思路通用性拉满！

中文标题：乳腺癌转移肿瘤内异质性单细胞景观及GABA介导免疫微环境重塑

发表期刊：Breast Cancer Research

发布时间：2026年3月

影响因子：6.2

01

研究背景

Hot

乳腺癌转移是女性肿瘤致死首要诱因，肿瘤神经 - 免疫轴重塑微环境是新兴热点。GABA 作为核心抑制性神经递质，可调控免疫细胞极化，但 GABA 介导转移灶肿瘤免疫微环境（TIME）的单细胞分子机制尚不清晰。现有研究多批量测序，缺少单细胞细分免疫亚群、细胞通讯与临床关联整合分析。本研究整合单细胞、批量转录组，借助 NMF、SCENIC、Monocle 等 AI 生信工具，解析 GABA 调控转移微环境的完整通路，挖掘预后标志物。

02

研究方法

Hot

1.多公共数据库数据获取：单细胞数据集 GSE158399（原发、转移阳性 / 阴性淋巴结）、独立验证单细胞队列 GSE180286；1203 例乳腺癌批量 RNA-seq 数据取自 UCSC Xena；1417 个 GABA 相关基因从 GeneCards 数据库检索。

2.单细胞基础生信流程Seurat 完成质控、标准化、PCA、UMAP 降维，Harmony 去除批次效应；基于 marker 基因注释 9 大类细胞；Monocle DDRTree 算法拟时序分析；CellChat 解析细胞配体 - 受体通讯网络。

3.AI 聚类与调控网络算法：NMF 非负矩阵分解细分 B/T/ 巨噬细胞 / 肌成纤维细胞 GABA 相关亚群；pySCENIC 预测各亚群转录因子活性；BEAM 筛选分化关键分支基因。

4.多组交叉临床验证：单细胞差异基因、GABA 基因集、批量肿瘤差异基因三者取交集，锁定核心基因 COL1A2，关联肿瘤分期、淋巴结转移状态评估临床价值；通路采用 HALLMARK、REACTOME 富集分析。

1

绘制乳腺癌转移单细胞全景

图 1 展示单细胞 UMAP 分群、细胞互作圈图、GABA 基因表达热图、NMF 亚群细分结果。质控后保留 30529 个细胞，注释出 T、B、巨噬细胞、肌成纤维细胞等 9 大细胞类型。基于 1417 个 GABA 基因集结合 NMF 算法，细分四大关键效应亚群：B_C2、CD4+Tcm、M2 巨噬细胞、Myo_C2，四类细胞 GABA 基因平均表达显著高于同谱系其他亚群。原发灶与淋巴结细胞占比差异明显：淋巴结样本 B、T 细胞占比上升，肌成纤维细胞、M2 巨噬细胞占比下降；独立数据集 GSE180286 完全复现该细胞分布规律，证明分型稳定可靠，明确 GABA 主要作用的靶细胞亚群，为后续机制解析划定核心细胞对象。

2

拟时序揭示 GABA 调控免疫细胞分化轨迹，转录因子存在亚群特异性

图 2-5 分别为 B 细胞、T 细胞、巨噬细胞、肌成纤维细胞拟时序轨迹、GABA 基因表达、亚群占比、BEAM 分支差异基因。Monocle 轨迹分析显示：B 细胞分化末端 B_C4 在转移淋巴结占比大幅降低；CD4+Tcm 全程高表达 GABA 基因；M2 巨噬细胞 GABA 相关分子显著高于 M1；Myo_C2 仅存在原发灶。SCENIC 转录因子分析发现各亚群核心调控分子完全区分：B_C4 富集 XBP1，M1 高表达 BACH1、IRF7，M2 富集 NR1H3、STAT1，Myo_C1 激活 WNT、Hedgehog 通路。BEAM 筛选上百个分化关键基因，HALLMARK 富集证实 M2、Myo_C2 富集血管生成、EMT、免疫抑制通路，阐明 GABA 通过调控转录程序重塑免疫细胞分化，诱导促转移表型。

3

GABA 重塑细胞通讯网络，三大核心配体受体轴驱动免疫抑制微环境

图 6 展示全细胞通讯强度、信号模式分类、关键配体受体气泡图、PTN/MDK 通路调控网络。肌成纤维细胞是整个微环境核心信号输出细胞，通过多通路和 T、B、巨噬细胞广泛互作。NMF 将细胞通讯分为 2 大模式，模式 1 以肌成纤维细胞为主，富集 PTN/MDK、CXCL12、APP 信号；模式 2 以免疫细胞为主，补体、SPP1 通路占主导。其中 PTN/MDK-NCL、CXCL12-CXCR4、APP-CD74 是跨细胞最关键互作对；ANXA1 介导肌成纤维细胞与巨噬细胞交流。PTN、MDK 信号主要由 Myo_C2/3 释放，作用于 Myo_C1 放大促癌信号。上述通路协同抑制 CD8+T 功能、促进 M2 极化，完整解释 GABA 如何通过细胞间信号串扰构建转移允许型微环境。

4

交叉多组筛选 COL1A2，可作为乳腺癌转移与分期独立标志物

图 7 为三组差异基因韦恩图、COL1A2 在不同分期 / 淋巴结分组表达箱线。将单细胞 BEAM 分化基因、批量癌正常差异基因、GABA 基因集三者交集，唯一得到核心靶基因 COL1A2。1203 例大样本临床数据验证：Ⅲ 期乳腺癌 COL1A2 表达显著高于 Ⅰ、Ⅱ 期；N2、N3 高转移淋巴结组表达远高于无转移 N0 组。COL1A2 编码胶原，参与 ECM 重塑、血管生成，与前文肌成纤维细胞、M2 巨噬细胞促转移通路高度契合。该基因同时连接 GABA 神经免疫轴与临床肿瘤进展，兼具分子机制与临床分层双重价值，为乳腺癌转移预后检测提供单一候选标志物。

04

文章小结

Hot

本研究依托 GEO、UCSC Xena、GeneCards 多公共数据库，融合单细胞与批量转录组，运用 NMF、SCENIC、Monocle 等 AI 生信算法，首次系统解析 GABA 调控乳腺癌转移免疫微环境机制。锁定 4 类 GABA 高表达关键细胞亚群，阐明 PTN/MDK 等三大促转移通讯轴，筛选出临床标志物 COL1A2。研究仅公共数据集、缺少体内外实验，后续需功能验证，为乳腺癌神经免疫靶向治疗提供全新生信依据。