多组学揭示乳腺癌线粒体功能障碍分子网络

原文：Multi-omics studies identified the molecular network of mitochondrial dysfunction in the pathogenesis of breast cancer.
期刊：Genes & diseases（影响因子 14.6） |  PubMed | PMC 全文（含 PDF）

乳腺癌是全球女性最常见的恶性肿瘤之一，其发生发展涉及遗传、表观遗传、代谢和肿瘤微环境等多层面异常。近年来，线粒体不再被仅仅视为“细胞的能量工厂”，其在肿瘤代谢重编程、氧化应激、凋亡抵抗以及免疫调控中的作用日益受到关注。已有研究提示，线粒体功能障碍与乳腺癌的发生有关，但具体是哪些线粒体相关基因、通过怎样的调控层次影响乳腺癌风险，仍缺乏系统性的因果证据。

近日，一项发表于 Genes & Diseases（影响因子 14.6）的研究通过整合多组学数据、孟德尔随机化与共定位分析，系统描绘了线粒体功能障碍参与乳腺癌发病的分子网络，并进一步锁定了一个具有治疗潜力的关键基因——UQCRC1。该研究由来自中国医学科学院肿瘤医院等机构的研究团队完成，为理解乳腺癌中线粒体相关调控机制提供了新视角。

研究背景：从“相关性”走向“因果性”

线粒体拥有自身的基因组，并依赖大量核编码蛋白维持其结构与功能。既往观察性研究发现，线粒体DNA突变、拷贝数变化、呼吸链功能异常以及线粒体相关核基因表达改变，都可能出现在乳腺癌组织中。然而，传统流行病学和转录组关联分析往往难以区分“因果驱动”与“疾病继发改变”。

为解决这一问题，研究团队采用了基于汇总数据的孟德尔随机化（Summary-based Mendelian Randomization, SMR）结合共定位分析（colocalization）的策略。该方法利用遗传变异作为工具变量，从DNA甲基化、基因表达和蛋白质数量性状位点（pQTL）三个层面，推断线粒体相关分子特征与乳腺癌风险之间的因果关系，而非简单相关。

核心发现一：甲基化与表达层面的因果关联

研究首先系统分析了线粒体相关基因的DNA甲基化、表达和蛋白调控与乳腺癌的因果联系。在甲基化层面，经过HEIDI检验和多重校正后，共有33个CpG位点（邻近16个独特基因）达到边际显著阈值。进一步共定位分析显示，其中4个基因具有强共定位证据（PPH4 > 0.70），包括：

• BCL2L11（cg27608154）
• BIK（cg01341487、cg04042468、cg26453588）
• MTERF1（cg17054360）
• TRMT1（cg12014333、cg22237401）

这些基因涉及线粒体凋亡调控、线粒体转录终止以及RNA修饰等过程，提示表观遗传层面的线粒体基因调控可能影响乳腺癌易感性。

在基因表达层面，SMR分析经FDR校正后鉴定出13个与乳腺癌风险相关的基因，包括 MTX1、MRPS18C、FDPS、HSCB、MPC2、MUTYH、UQCRC1、SURF1、MRPS16、NSUN2、NIT2、NUDT13 和 COX5A。其中，共定位分析提供了强证据支持三个基因：UQCRC1、MPC2 和 MRPS18C。这三个基因分别参与线粒体呼吸链复合体III组装、丙酮酸转运以及线粒体核糖体蛋白合成，说明线粒体能量代谢与翻译机器在乳腺癌中具有因果作用。

值得注意的是，在蛋白质层面，虽然初步筛选出 KYAT3、FH、NIT2、PPA2、DNAJC30、NUDT2、ACADVL、DCXR、SARS2 等9个线粒体相关蛋白达到名义显著，但在多重检验校正后，没有蛋白保留显著关联。因此，研究最终将治疗靶点聚焦于转录组层面的 UQCRC1、MPC2 和 MRPS18C。

核心发现二：UQCRC1 是保护性因子与潜在治疗靶点

在三个候选基因中，UQCRC1（泛醇-细胞色素C还原酶核心蛋白1）被进一步作为重点研究对象。临床样本分析显示，UQCRC1 在乳腺癌组织中表达降低，且低表达与不良预后相关；免疫组化结果也验证了这一表达模式。在不同乳腺癌亚型中，UQCRC1 表达水平存在差异，并且随着其表达升高，患者肿瘤分期和分级反而下降，提示 UQCRC1 可能发挥“保护因子”的作用。

功能分析表明，UQCRC1 主要参与能量代谢、RNA调控和线粒体活性。进一步的免疫微环境分析发现，UQCRC1 表达与 CD4⁺ T细胞、肥大细胞、CD8⁺ T细胞及中性粒细胞浸润相关；同时，CTLA4、CD276、LAG3 等免疫检查点分子也与其表达密切相关。基于 UQCRC1 构建的临床模型在多个队列中表现出较好的免疫治疗预测效能。

研究还发现，UQCRC1 与靶向药物 capivasertib 在多个队列中可能存在协同关系，或有助于克服 capivasertib 耐药；此外，PIK3CA 突变与 UQCRC1 密切相关，进一步支持 PI3K 通路抑制剂与 UQCRC1 之间的协调效应。

机制解析：UQCRC1 通过抑制 PI3K/AKT 信号抑制肿瘤

为验证 UQCRC1 的生物学功能，研究团队在 MCF7 和 T47D 乳腺癌细胞系中进行了敲低实验。结果显示，敲低 UQCRC1 后，细胞增殖（MTS实验）、迁移（划痕实验）和侵袭（Transwell实验）能力均显著增强，说明 UQCRC1 能够抑制乳腺癌细胞的恶性表型。

机制上，对转移与非转移乳腺癌组织进行RNA测序，并根据 UQCRC1 中位表达分组进行差异分析，发现差异基因显著富集于代谢通路、PI3K/AKT 和 mTOR 信号。Western blot 进一步证实，在 UQCRC1 敲低细胞中，磷酸化 PI3K（P-PI3K）和磷酸化 AKT（P-AKT）水平升高，而总 PI3K 和 AKT 蛋白量不变。这表明，UQCRC1 通过抑制 PI3K/AKT 信号通路的活化，进而抑制乳腺癌细胞增殖和迁移。

方法亮点：多组学整合与药物靶点预测

该研究的方法学框架具有明显创新：一方面，它同时利用甲基化QTL、表达QTL和蛋白QTL公共数据库，并结合 IEU GWAS 项目中乳腺癌汇总统计数据，实现了从“表观—转录—蛋白”多层次因果推断；另一方面，通过 HEIDI 检验排除连锁不平衡导致的假阳性，并以共定位分析（PPH4 > 0.70）强化因果证据。

在潜在转化方面，研究利用 DSigDB 进行药物富集分析，并借助 Cytoscape 构建蛋白—药物网络。结果显示，针对 UQCRC1 的化合物包括 1-Methyl-4-phenyl-2,3-dihydropyridinium、cube root extract 以及化合物 635-65-4；针对 MRPS18C 的化合物为 4-苯基丁酸（4-phenylbutyric acid）。分子对接进一步确认了配体—受体之间具有较低结合能、较高亲和力的相互作用，为后续药物重定位提供了线索。

应用意义：从生物标志物到联合治疗策略

这项研究的价值不仅在于鉴定了新的乳腺癌风险相关线粒体基因，更在于把“线粒体功能障碍”这一宽泛概念落实到了具体分子和通路上。UQCRC1 同时具备以下特征：在肿瘤中低表达、与良好病理特征相关、受甲基化/表达层面因果支持、可调控 PI3K/AKT 通路、关联免疫微环境及免疫治疗反应，并与现有靶向药存在协同潜力。

因此，UQCRC1 有望成为乳腺癌预后评估的生物标志物，以及克服 PI3K/AKT 通路抑制剂耐药的新靶点。当然，作者也指出，蛋白层面未发现显著因果靶点，且部分候选化合物仍需实验验证；未来仍需要更大样本和功能研究来确认这些多组学发现。

总体而言，该研究通过多组学因果推断框架，绘制了乳腺癌中线粒体相关分子网络，并确立 UQCRC1 作为连接线粒体功能、PI3K/AKT 信号与肿瘤进展的关键节点，为乳腺癌的机制研究和靶向干预提供了新方向。

本文图片转载自 PMC 开放获取全文，仅供学习参考，请以原文为准。