Nat Commun 高分中介孟德尔随机化全 R 管线拆解附核心图核心图一览Fig 1–6含 Fig 3 与 Fig 6下图从上到下依次对应分析管线的四个关键环节研究设计流程图Fig 1→ 体型/初潮年龄→乳腺密度 UVMRMVMRFig 2→ 乳腺密度→乳腺癌 UVMRMVMRFig 3→ 中介效应分解Fig 6。分析管线拆解一、数据来源与样本设计结局主数据BCAC 2017 乳腺癌汇总统计OpenGWAS IDieu-a-1126 总样本 / ieu-a-1127 ER / ieu-a-1128 ER−总样本 N 228,951亚型扩展数据BCAC 2020 分子亚型Zhang et al. 发布需向 BCAC 申请暴露数据儿童期体型ieu-b-510、成年体型、初潮年龄ukb-b-3768均为 OpenGWAS 公开 ID中介数据乳腺密度 GWAS致密面积 DA / 非致密面积 NDA / 百分密度 PD仓库内提供预处理与 BMI 校正版本的使用逻辑设计亮点MD GWAS不校正成人 BMI的版本作为主分析避免 MVMR 中双重校正的共线陷阱二、数据处理与特征计算工具变量IV筛选沿用各 GWAS 原文阈值最终 IV 数量儿童期体型 115、成年体型 173、初潮年龄 190、致密面积 21、非致密面积 8、百分密度 11暴露与结局 ID 对齐、等位基因方向核对由仓库内 lookup 表 函数脚本统一完成中介分析所需的总效应/直接效应估计在同一套 IV 定义下产出保证 Product 法可乘三、统计模型总览主模型两样本 IVWUVMR与 MVMR多变量 MR控制共暴露双向 MR排除结局→暴露的反向因果敏感性MR-Egger、weighted median、MR-PRESSO、RadialMR、mrclust中介Product 法分解总效应为经致密面积的间接效应 直接效应报告中介占比及 95% CI四、主要分析目标Aim 1–3Aim 1儿童期/成年体型及初潮年龄对三个 MD 表型的因果效应UVMR 互相控制的 MVMR——Fig 2Aim 2三个 MD 表型对乳腺癌总样本及亚型的因果效应含控制儿童期体型后的直接效应——Fig 3Aim 3致密面积是否中介儿童期体型→乳腺癌的保护效应中介 MR 异质性校正后的稳健估计——Fig 6五、结果呈现逻辑儿童期体型降低致密面积成年体型升高非致密面积——两条通路方向清晰分开致密面积升高因果性增加乳腺癌风险总样本与 ER 亚型均有证据中介 MR 显示儿童期体型对乳腺癌的保护效应中56%95% CI 32%–79%经由致密面积传递敏感性分析发现致密面积 IV 存在异质性MR-PRESSO 校正后中介比例结论方向不变PheWAS mrclust 对致密面积 IV 做了表型组扫描排除明显多效性混杂六、一句话总结OpenGWAS 暴露/结局 联盟 MD GWAS → IVW 锁定双向因果 → MVMR 控共暴露 → MR-PRESSO/mrclust 洗多效性 → Product 法算出 56% 中介占比。写在最后框架已经全部摊开但每一步的参数细节、避坑点、结果核对清单没法在一篇短文里穷尽。如果你正在做同类暴露→中介→结局中介 MR 找不到完整可跑管线毕业论文需复现并迁移到自己的 GWAS 组合系统学 IVW/MVMR/中介分解这一套统计方法孟德尔随机化代做、辅导都在淘宝店「万汇泛科技」拍前可先私信你的数据类型帮你判断能不能套这套管线。后续还会拆更多「开源可复现」的高分文章关注不迷路。
