《Journal of Cellular and Molecular Medicine》：Integrated Multi-Omics Analyses Reveal Causal DNA Methylation-Breast Cancer Associations Mediated by Plasma Proteins
编辑推荐：
尽管异常DNA甲基化（DNA methylation, DNAm）与乳腺癌（breast cancer, BC）发病相关，但因果证据不足且下游介质不明。研究人员通过全基因组血液甲基化数据鉴定差异甲基化位点（differentially methylated p
尽管异常DNA甲基化（DNA methylation, DNAm）与乳腺癌（breast cancer, BC）发病相关，但因果证据不足且下游介质不明。研究人员通过全基因组血液甲基化数据鉴定差异甲基化位点（differentially methylated positions, DMPs），并以甲基化数量性状位点（methylation quantitative trait locus, meQTL）为工具变量，结合BCAC与FinnGen两个独立联盟的BC遗传数据开展两样本孟德尔随机化（two-sample Mendelian randomization, TSMR），辅以广义汇总数据孟德尔随机化（Generalized Summary-data-based Mendelian Randomization, GSMR）、贝叶斯共定位（Bayesian colocalization）与表型组宽孟德尔随机化（Phenome-wide Mendelian randomization, PheWAS）验证。进一步通过蛋白质组范围中介分析与单细胞RNA测序（single-cell RNA sequencing, scRNA-seq）解析下游通路。研究人员共识别19,169个DMPs，MR支持95个CpG位点与BC一致相关；荟萃（meta）分析确认三个高置信位点（cg15146122，OR=0.91；cg01529201，OR=1.14；cg24704912，OR=1.11），其中两个获共定位支持。中介分析鉴定203个BC相关蛋白，三个CpG分别调控41至314个蛋白，涉及免疫调节、凋亡与溶酶体通路，scRNA-seq将介质蛋白定位于肿瘤微环境特定细胞群。结果表明，血液DNAm经特定血浆蛋白以MR支持的方式关联BC风险，为表观遗传标志物提供依据。
研究背景与问题提出
乳腺癌（breast cancer, BC）是全球女性高发恶性肿瘤，现有筛查手段在灵敏度与假阳性控制上仍有限。DNA甲基化（DNA methylation, DNAm）作为连接遗传易感与环境暴露的表观遗传机制，其异常——尤其启动子高甲基化导致的抑癌基因沉默——参与肿瘤起始与进展。然而既往证据多来自观察性研究，难以区分因果驱动与疾病继发改变，且缺乏系统性的全基因组因果定位，以及从甲基化到疾病的可解释下游通路。
研究人员开展的工作与总体结论
研究人员开展一项多组学因果推断研究：从外周血甲基化谱鉴定差异甲基化位点（differentially methylated positions, DMPs），以甲基化数量性状位点（methylation quantitative trait locus, meQTL）为工具变量（instrumental variable, IV），整合BCAC与FinnGen的汇总级遗传数据，经两样本孟德尔随机化（two-sample Mendelian randomization, TSMR）、广义汇总数据孟德尔随机化（Generalized Summary-data-based Mendelian Randomization, GSMR）、贝叶斯共定位（Bayesian colocalization）、表型组宽关联（Phenome-wide association study, PheWAS）识别高置信因果CpG；继而用蛋白质数量性状位点（protein quantitative trait locus, pQTL）数据做蛋白质组范围中介，并用单细胞RNA测序（single-cell RNA sequencing, scRNA-seq）定位介质蛋白的细胞来源。研究结论为：血液DNAm经特定血浆蛋白介导与BC风险存在MR支持的因果关系，三个高置信CpG可作为候选表观遗传标志物。
主要关键技术方法
研究基于公开汇总数据与非独立样本队列：甲基化阵列样本来自GEO的GSE104942（87例BC与123例对照外周血）；meQTL来自GoDMC（27,750人）；BC结局来自BCAC与FinnGen R11；pQTL来自UKB-PPP（33,628人，2,940个血浆蛋白）；scRNA-seq来自GEO的GSE161529（52例乳腺肿瘤与癌旁）。方法上依次采用DMP差异分析、TSMR（IVW/Wald ratio）、GSMR与HEIDI、双队列meta分析、共定位（PP.H4>0.8）、PheWAS、两步MR中介比例计算，以及Seurat标准化、Harmony批次校正与DoubletFinder双细胞剔除的单细胞注释。
研究结果
3.1 全基因组差异甲基化分析
经预处理与线性模型检验，研究人员在BC与对照间鉴定19,169个DMPs，其中14,118个高甲基化、5,051个低甲基化。
3.2 两样本MR交叉验证因果CpG
以GoDMC meQTL为IV，分别在FinnGen与BCAC中运行TSMR，得到705与872个显著CpG，交集为95个跨队列一致因果CpG，敏感性分析未见明显多效性或异质性。
3.3 GSMR验证
对95个CpG按IV数量质控后，25个进入GSMR，其中13个被复制，8个在两队列方向一致，证明估计非方法依赖假象。
3.4 Meta分析整合
对8个高置信位点合并两队列估计，三个位点显著：cg15146122（OR=0.91，保护效应）、cg01529201（OR=1.14）、cg24704912（OR=1.11，风险效应）。
3.5 遗传共定位
cg15146122（rs440155，PP.H4=0.99）与cg24704912（rs7465455，PP.H4=0.99）存在共享因果变异；cg01529201未达阈值。
3.6 PheWAS多效性筛查
rs440155无显著脱靶；rs595224关联前列腺炎与药物不良事件；rs7465455关联卵巢功能障碍与神经痛，提示后两者解释需谨慎。
3.7 蛋白介导通路
三个CpG分别影响41、314、126个血浆蛋白；顶部介质为CD207、RPE、OCLN，功能富集于免疫应答、C型凝集素受体信号、凋亡、焦点黏附与溶酶体通路。
3.8 中介网络
两步MR确认15对显著间接效应：cg01529201经ANXA10，cg15146122经ALPI、ANXA10、FGF2、OLR1等13个蛋白，cg24704912经ANGPTL2介导BC风险。
3.9 染色体定位
PhenoGram显示CpG与介质基因多为远端trans调控，如cg01529201（chr22）调控ANXA10（chr4）。
3.10 差异甲基化验证
三个CpG在BC患者中均呈显著低甲基化，支持其作为检测标志物的可检测性。
3.11 单细胞微环境定位
UMAP分七类细胞；13个介质基因（除ALPI）呈细胞类型特异表达，FGF2与OLR1在多种细胞中差异分布，定位于上皮、免疫与基质细胞群体。
讨论与结论翻译
讨论指出，本研究首次将血液甲基化、血浆蛋白质组与BC风险纳入统一因果框架，发现三个新颖高置信CpG及下游蛋白介质，弥补了既往仅关注吸烟/饮酒相关位点的不足。优势在于大样本meQTL、双结局队列、多敏感性与单细胞验证；局限为欧洲血统偏倚、残留多效性、线性中介假设及缺实验功能验证。
结论原文意译：
综上，本多组学研究提供稳健且新颖的证据，支持血液DNAm在BC发病中的因果作用。TSMR经BCAC与FinnGen复制，并结合共定位与敏感分析，确认三个高置信CpG——cg15146122（保护）、cg01529201（风险）与cg24704912（风险）——为关键表观遗传决定因素；其效应经14个血浆蛋白、15对CpG-蛋白中介对实现，涉及免疫调节、凋亡与代谢通路。scRNA-seq将介质精确定位于肿瘤微环境细胞区室。未来需在多样人群中验证、用实验模型阐明机制，并推动表观遗传风险预测与靶向策略进入精准肿瘤学。
这篇论文解读已完成，接下来可以选择一个方向继续深化：
评估该研究的方法学局限与改进路径 设计三站点标志物的临床验证方案 评估该研究的方法学局限与改进路径 设计三站点标志物的临床验证方案
