三阴性乳腺癌的"坏细胞"与软肋:单细胞+空间转录组+机器学习锁定 RPN1
一句话看懂:三阴性乳腺癌(TNBC)异质性极强、又缺乏靶向药。这篇研究把单细胞转录组 + 空间转录组 + 机器学习三套武器组合起来,在 TNBC 里揪出一个高度恶性的上皮细胞亚群 Luminal_inflam,并顺藤摸瓜锁定了一个关键基因 RPN1。 更关键的发现:基于 Luminal_inflam 相关基因模块,研究团队训练出一个 12 基因风险模型,能有效区分患者生存;而 RPN1 不仅预后价值突出,更是癌细胞的"命门"——敲低它会让癌细胞陷入内质网应激(ER stress)、活性氧与钙离子失衡,最终走向凋亡。 对研究者的启示:单细胞负责"找到坏细胞",空间转录组负责"还原它在组织里的位置",机器学习负责"把它变成可用的临床标志物"——这套组合拳,正在成为肿瘤机制研究的新范式。
一、研究背景:三阴性乳腺癌为什么这么难缠
三阴性乳腺癌(triple-negative breast cancer,TNBC)是乳腺癌里最凶险的亚型之一——雌激素受体(ER)、孕激素受体(PR)和人表皮生长因子受体 2(HER2)均为阴性,意味着内分泌治疗和抗 HER2 靶向药都无从下手,患者预后差、复发风险高。
更棘手的是它的高度细胞异质性:同一块肿瘤里,不同癌细胞的行为天差地别,有的安分,有的极具侵袭性。想找到真正"作恶"的那群细胞,并针对它们开发治疗策略,就必须先把肿瘤拆到单细胞分辨率去观察。
二、研究设计:单细胞 + 空间转录组 + 机器学习的"组合拳"
研究团队从公共数据库获取了 TNBC 的单细胞转录组数据(GSE161529),共 21 个样本(13 例正常乳腺组织 + 8 例 TNBC 组织),从中解析出上皮细胞、T 细胞、B 细胞、浆细胞、髓系细胞、内皮细胞、成纤维细胞、周细胞等主要细胞类群。
在此基础上,研究分三步走:
单细胞
——在每个细胞里找出恶性上皮亚群,推断拷贝数变异(CNV)、分化潜能与基因调控网络; 空间转录组
——用 Visium 数据把关键细胞和基因"放回"组织原位,验证它们的空间共定位; 机器学习
——把恶性亚群相关的基因模块与生存队列结合,构建并验证预后风险模型。
三、第一步:画出 TNBC 的"细胞地图"

▲ 图1:TNBC 单细胞图谱——UMAP 清晰区分主要细胞类型,肿瘤微环境生态位重塑。
单细胞分析首先给出的是 TNBC 的整体细胞图谱。UMAP 降维把细胞按类型清晰分开,各细胞类型的经典标志基因表达符合预期,说明数据质量可靠。
对比正常组织和 TNBC 的细胞组成差异,可以看到肿瘤内部的生态位发生了明显重塑——某些细胞类型的比例此消彼长。在上皮细胞这个"主角"区室里,TNBC 与正常组织之间还存在大量差异表达基因(如 FGFBP2、NPTX1 上调,CYP24A1 等下调),提示恶性上皮细胞在转录层面早已"改头换面"。
四、锁定"坏细胞":Luminal_inflam 恶性上皮亚群

▲ 图2:锁定“坏细胞” Luminal_inflam——高 CNV、终末分化、细胞周期活跃的高度恶性上皮亚群。
把上皮细胞进一步细分后,研究识别出一个特殊的亚群——Luminal_inflam。它有几个鲜明的"恶性特征":
- 拷贝数变异(CNV)水平高:CNV 是基因组不稳定的标志,Luminal_inflam 的 CNV 得分显著高于其他上皮亚群,说明它的基因组更"乱"、更不稳定; - 处于终末分化状态:伪时间(pseudotime)分析显示它位于分化轨迹的末端,恶性的"定型"程度高; - 细胞周期相关程序富集:增殖相关转录程序活跃,暗示它增殖能力强、侵袭性高。
随后,研究者用 hdWGCNA(高维加权基因共表达网络分析)为 Luminal_inflam 定义了与之相关的基因模块(如 Epi-M1、Epi-M3),为后面的临床转化埋下伏笔。
五、肿瘤微环境的"悄悄话":S100A4–EGFR 轴

▲ 图3:微环境“悄悄话”——成纤维细胞 S100A4–EGFR 轴可能给恶性细胞“递刀子”。
癌细胞不是孤岛。细胞间通讯(cell–cell communication)分析发现,TNBC 的信号传递网络发生了系统性重塑——谁在"发信号"、谁在"收信号"、整体连接强度,都与正常组织不同。
其中最受关注的,是一条成纤维细胞来源的 S100A4–EGFR 轴:成纤维细胞分泌 S100A4,作用于上皮细胞的 EGFR 受体,可能反过来调控 Luminal_inflam 相关基因的表达。换句话说,周围的"邻居细胞"可能正在给恶性细胞"递刀子"。这条轴被研究者列为 TNBC 微环境重编程的重点候选机制。
六、从细胞到病床:12 基因风险模型

▲ 图4:12 基因风险模型——在 TCGA 等多个独立队列有效分层患者总生存期。
找到恶性亚群,下一步是把它变成临床上可用的工具。研究者把 Luminal_inflam 的基因模块与生存数据(TCGA、GSE58812 队列)交叉,先筛出与生存相关的基因,再调用多种机器学习分类器(SVM、LASSO 等)层层筛选。
最终得到一个 12 基因风险模型(12-gene risk signature)。在多个独立队列中,该模型都能有效把患者按总生存期分层——高风险组预后显著更差。此外,高低风险组在免疫特征和潜在药物敏感性上也存在差异,为后续免疫治疗与用药选择提供了线索。
七、锁定主角:RPN1 是个"高危指标"

▲ 图5:主角 RPN1——单细胞、空间、蛋白与生存数据一致显示其在 TNBC 高表达并预示不良预后。
在 12 个基因里,RPN1 因为风险比(hazard ratio)最高而格外抢眼。它到底是什么?
RPN1 编码的是一种内质网 I 型跨膜蛋白,是寡糖基转移酶(OST)复合体的核心亚基,负责新生蛋白的 N-连接糖基化——可以说是细胞"蛋白加工厂"里的关键工位。既往研究已发现 RPN1 在多种恶性肿瘤中高表达、与不良预后相关;在 TNBC 中,它还参与 PD-L1 的糖基化与稳定性,与免疫逃逸有关。
本研究的多层次证据(单细胞、空间转录组、蛋白表达、KMplotter 生存数据)一致指向:RPN1 在 TNBC 中高表达,且高表达预示不良预后。
八、致命软肋:敲低 RPN1 触发内质网应激

▲ 图6:RPN1 的软肋——敲低触发 PERK–ATF4–CHOP 内质网应激,癌细胞陷入凋亡。
最精彩的部分在"功能验证"。研究团队在 MDA-MB-231 细胞(一种经典 TNBC 细胞系)里敲低 RPN1,结果癌细胞的表现急转直下:
- 活力下降(CCK-8 实验),且可被 4-PBA(一种内质网应激抑制剂/化学伴侣)部分挽救; - 活性氧(ROS)与胞内钙离子(Ca²⁺)升高,提示细胞稳态失衡; - 细胞周期分布改变,增殖受阻; - 凋亡显著增加(Annexin V/7-AAD 检测)。
机制上,敲低 RPN1 激活了 PERK–ATF4–CHOP 通路——这是内质网应激(ER stress)/未折叠蛋白反应(UPR)的经典信号轴。也就是说,RPN1 一旦缺失,癌细胞的"蛋白加工厂"就乱了套,最终把自己"逼上绝路"。这正是标题所说的 RPN1–PERK UPR 脆弱性(vulnerability)。
九、给研究者的启示
单细胞 + 空间 + 机器学习的"三件套"正在成为肿瘤研究的标配
。单细胞负责发现稀有恶性亚群,空间转录组负责提供组织学语境,机器学习负责把发现转化为可用的预后/预测模型——三者互补,缺一不可。 CNV + 伪时间 + 共表达模块,是"定性"恶性亚群的一套组合拳
。光看差异表达还不够,把基因组不稳定性、分化状态和调控网络叠加起来,才能把"坏细胞"钉死。 从"标志物"到"可成药靶点",功能验证是关键一步
。12 基因模型解决了"预后分层",但真正让论文立住的,是 RPN1 敲低实验——它把一个统计相关的基因,变成了一个有明确表型和机制的脆弱点。对客户项目而言,多组学发现 + 湿实验闭环,才是把故事讲深的标准动作。
结语
这项研究用单细胞 + 空间转录组 + 机器学习的组合拳,在三阴性乳腺癌中画出了一条清晰的逻辑链:从"找到恶性细胞 Luminal_inflam",到"构建 12 基因风险模型",再到"锁定并验证 RPN1–PERK 内质网应激这一致命软肋"。
对 TNBC 这样缺乏靶点的"硬骨头"来说,RPN1 所代表的蛋白质稳态(proteostasis)脆弱性,或许正是下一个值得深挖的治疗方向。
▲原文标题:Integrative single-cell and spatial transcriptomics with machine learning identify a Luminal-inflam malignant program and reveal an RPN1–PERK UPR vulnerability in triple-negative breast cancer。



English

