慢性阻塞性肺疾病与骨质疏松症遗传因果关联的孟德尔随机化与共定位整合分析

《Biomedicines》：A Synthetic/Fibrotic Vascular Remodeling Program Associated with Arteriovenous Fistula Maturation Failure

　　

摘要 背景：慢性阻塞性肺疾病（COPD）与骨质疏松症高度共病，但其遗传因果方向及共享机制尚未明确。方法：本研究整合了全球生物库荟萃分析（GBMI）的COPD汇总统计数据（58,559例欧洲血统病例/937,358例对照）、ReGGen东亚骨密度（BMD）荟萃分析（365,403例个体）及FinnGen R12数据（骨质疏松症10,461例病例；病理性骨折2,429例病例）。工具变量选择阈值为p < 5 × 10?8，并通过连锁不平衡（LD）剪枝（基于1000 Genomes参考面板，r2 < 0.001，窗口10 Mb）；预先设定了每工具变量的F统计量、校正肝脏疾病遗传倾向的多变量孟德尔随机化（MVMR）以及排除PIZ变异（rs28929474）的敏感性孟德尔随机化分析。共定位分析采用贝叶斯近似贝叶斯因子（ABF）方法的病例-对照形式，并对先验概率p12进行敏感性分析（取值集合{1 × 10?6, 1 × 10?5, 1 × 10?4}）。结果：遗传预测的COPD显示出方向一致的阳性关联：固定效应逆方差加权法（IVW）比值比（OR）= 1.226（95%置信区间CI 1.095–1.373，p* = 4.2 × 10?4）；加权中位数法p = 0.016；校正肝脏疾病的多变量IVW法p = 0.040；排除PIZ后OR = 1.160，p* = 0.013。在乘法随机效应模型下，估计值减弱至边缘不显著（OR = 1.169，95% CI 0.997–1.371，p = 0.055；Cochran Q = 46.6，p = 0.008）；骨折关联呈名义显著性，未通过针对五项主要检验的Bonferroni校正（p = 0.042，alpha = 0.01）（OR = 1.390，p = 0.042）。反向因果关系检验结果为阴性（BMD对COPD的OR = 0.930，p = 0.127；FinnGen内骨质疏松症对COPD的OR = 1.031，p = 0.503）。共定位分析不支持任何区域在任何先验设定下达到后验概率H4（PP.H4）> 0.5；在SERPINA1基因座，后验概率倾向于H3（PP.H3 = 0.72），即存在不同的因果变异而非单一共享变异，因此PIZ被提出作为候选的方向一致基因座，而非经证实的共享因果变异。与FEV1%预测值相关的肺部共表达模块（32个模块中有21个经Benjamini-Hochberg校正显著；最强模块beta = ?11.06每标准差，p* = 5.4 × 10?21），而骨髓间充质干细胞差异表达及外周血单细胞分析（补充材料，探索性）仅为探索性和假设生成性质。结论：本分析提供了微弱、方向一致且受遗传学支持的证据，将COPD遗传易感性与骨质疏松症风险联系起来，但该结果对工具变量间的异质性敏感，且未发现统计学上显著的共享因果基因座。临床转化应视为假设生成性质；在COPD人群中进行骨骼健康监测以及对SERPINA1携带者队列进行随访是未来的研究方向，而非既定推荐。稳健的敏感性估计量（加权中位数、加权众数、MR-Egger、基于1000次自助法分布的MR-PRESSO、MR-RAPS、污染混合模型）均被预先设定并在文中完整报告。在SERPINA1基因座，一项靶向条件分析利用来自1000 Genomes高覆盖度欧洲面板（n = 503）的LD，对骨骼侧（M13_OST）统计量基于rs28929474进行了条件化处理，并在条件化后的统计量上重复了共定位分析。五项预先设定的主要检验在alpha = 0.05/5 = 0.01水平（Bonferroni）进行校正：F1检验仍显著（p = 4.2 × 10?4），而骨折关联（F2，p = 0.042；Cochran Q p = 0.997，无异质性）未通过校正，故被解释为名义显著性。所有针对F1检验的离群值鲁棒估计量均保持方向一致：MR-RAPS OR = 1.211（p = 0.0017），污染混合模型OR = 1.294（95% CI 1.039–1.565，p = 0.040），加权中位数OR = 1.297（p = 0.0034），以及MR-PRESSO离群值校正估计OR = 1.190（95% CI 1.060–1.337，p = 0.0064；全局异质性p = 0.006）。

论文解读：慢性阻塞性肺疾病与骨质疏松症遗传因果关联的孟德尔随机化与共定位整合分析

研究背景与问题

慢性阻塞性肺疾病（COPD）是全球第四大死因，以持续性气流受限和慢性炎症为特征；骨质疏松症是一种以低骨量和骨组织微结构退化为特征的全身性骨骼疾病，可导致骨折风险增加。流行病学研究报道，COPD人群中骨质疏松症及椎体/髋部骨折的风险升高2至3倍，糖皮质激素暴露、体力活动减少、维生素D缺乏、吸烟及慢性炎症被认为是可能的候选机制。然而，观察性研究无法解决COPD与骨质疏松症之间是否存在因果关系及其作用方向。先前的孟德尔随机化（MR）研究聚焦于肺功能与骨密度（BMD），但信号混杂，提示需要针对特定性状和结局进行设计，并辅以明确的敏感性分析。此外，α-1抗胰蛋白酶（SERPINA1）基因轴为肺-骨耦合提供了一个具体的生物学假说：由PI*Z变异（rs28929474 Glu342Lys）驱动的蛋白酶-抗蛋白酶失衡导致早发性肺气肿，且α-1抗胰蛋白酶缺乏症患者据报道伴有包括骨质疏松症和骨折在内的额外肌肉骨骼疾病负担。然而，该基因座是否在群体水平上贡献于肺-骨共病尚不清楚。因此，本研究旨在通过改进的分析框架，利用大规模遗传数据进行双向MR和共定位分析，以评估COPD与骨质疏松症之间的潜在因果联系。

研究方法概述

研究人员采用了多项关键技术方法。首先，利用孟德尔随机化（MR） 框架，从多个公共数据库中获取汇总统计数据，包括：全球生物库荟萃分析（GBMI）的COPD数据（58,559例欧洲血统病例/937,358例对照）、ReGGen东亚BMD荟萃分析（365,403例东亚个体）以及FinnGen R12数据库（骨质疏松症10,461例病例；病理性骨折2,429例病例；COPD 13,535例病例）。工具变量选自全基因组显著性阈值（p < 5 × 10^?8）的变异，并使用1000 Genomes参考面板进行连锁不平衡（LD）剪枝。其次，应用了多种敏感性分析方法，包括加权中位数、MR-Egger、MR-PRESSO、MR-RAPS和污染混合模型，并进行排除PIZ变异（rs28929474）的敏感性分析和多变量MR（MVMR），以调整肝脏疾病遗传倾向。第三，采用贝叶斯近似贝叶斯因子（ABF）共定位分析*的病例-对照形式，在先验概率p12的不同设定下评估共享因果变异的可能性。此外，还利用了肺组织转录组数据集（GSE76925，151个样本）进行共表达聚类分析，以及骨髓间充质干细胞（MSC）转录组数据集（GSE35956，5例骨质疏松症 vs. 5例对照MSC）进行差异表达和跨组织富集分析。

研究结果

3.1 研究设计、工具变量构建与损耗

在GBMI数据流中，1216个变异达到显著性阈值（584个经聚类后的先导变异），经LD剪枝保留30个独立候选工具变量，其中27个与FinnGen骨质疏松症结局成功协调。对于东亚测试，120,100个ReGGen先导变异进入流程，最终162个独立工具变量与COPD结局协调。反向测试保留了13个骨质疏松症先导变异，其中10个成功协调。

3.2 正向与双向MR分析，含PIZ排除与MVMR*

测试F1（GBMI COPD对骨质疏松症）显示，固定效应IVW的OR为1.226（95% CI 1.095–1.373，p = 4.2 × 10^?4），但存在显著的异质性（Cochran Q p = 0.008），乘法随机效应估计值变为边缘不显著（OR = 1.169，p = 0.055）。加权中位数估计方向一致（p = 0.016），MR-Egger截距无显著定向多效性（p = 0.164）。排除PIZ变异后关联减弱但未消失（OR = 1.160，p* = 0.013）。在MVMR中，调整肝脏疾病遗传倾向后，COPD对骨质疏松症的效应仍显著（p = 0.040）。测试F2（COPD对病理性骨折）呈名义显著性（OR = 1.390，p = 0.042），但未通过多重检验校正。东亚测试F3（BMD对COPD）为阴性（OR = 0.930，p = 0.127）。反向测试R1（骨质疏松症对COPD，FinnGen内）为阴性（OR = 1.031，p = 0.503）。稳健估计量堆栈在所有测试中均保持方向一致。

3.3 共定位分析：无共享因果支持；SERPINA1为H3主导

在标准先验设定下，98个区域中无一达到PP.H4 > 0.5。在SERPINA1基因座（先导变异rs28929474），后验概率倾向于H3（不同因果变异）= 0.72，而非H4（共享因果变异）。该区域包含强COPD侧信号和较弱的、方向一致的骨质疏松症侧关联（beta_OUT = 0.221，p = 4.6 × 10^?6），但共定位模型将其归因于不同的因果变异。先验敏感性分析未改变结论。因此，研究人员撤回此前关于SERPINA1存在共享因果变异的声明，将PI*Z报告为方向一致的候选基因座。

3.4 肺部共表达模块与肺功能

在32个肺部共表达模块中，21个与FEV1%预测值显著关联（经Benjamini-Hochberg校正）。模块M133（34个探针）显示出最强的负相关性（beta = ?11.06 FEV1%每模块特征基因标准差，p = 5.4 × 10^?21）。无模块与二元COPD状态显著关联。

3.5 骨骼侧差异表达与跨组织富集（名义显著性）

骨髓MSC对比产生了354个名义显著的探针（p < 0.001），但无探针达到预设的FDR < 0.1阈值（最小FDR为0.108）。在87项跨组织富集测试中，最强的是下调基因集在模块M29中的富集（OR = 5.13，p = 0.0016），但经校正后仅为名义显著性（最小BH-q = 0.076）。先前重点关注的六基因集（CASP6、DIAPH3、DEK、SPHAR、SUMO1、SNAP23）的功能富集分析未产生任何经校正后显著的条目。

讨论与结论总结

讨论：本修订版分析报告了微弱、方向一致的遗传证据，表明COPD易感性可能与骨质疏松症风险相关。该关联在固定效应IVW、加权中位数、MVMR及排除PIZ等分析中呈名义上的稳健性，但对工具变量间的异质性敏感（随机效应p* = 0.055）。反向因果关系检验结果为阴性或不可估计。研究人员认为，此模式与一个适度的肺到骨效应成分相符，但在当前样本量下无法将其与吸烟相关多效性的残余混杂区分开来。未测量的混杂因素，特别是吸烟表型和系统性糖皮质激素暴露，仍是此类MR设计的主要威胁。共定位再分析是本次修订的核心，结果显示无区域支持共享因果变异，SERPINA1基因座的后验概率倾向于存在不同因果变异（H3）。此前声称的共享因果变异源于实施错误，研究人员对此表示感谢审稿过程的监督。PI*Z本身因其方向一致的骨质疏松症侧关联及在α-1抗胰蛋白酶缺乏症中观察到的额外肌肉骨骼疾病负担而仍具生物学趣味，但群体水平的共享因果解释不受支持。转录组学分析层被明确视为假设生成性质。

结论：全部证据与一个微弱、方向一致的COPD到骨骼的遗传效应相一致，在当前样本量下无法将其与吸烟相关多效性的残余混杂区分开来；这些发现属于假设生成性质，需要在无重叠的生物库中进行独立复制。在COPD人群中前瞻性评估骨骼健康监测以及对SERPINA1携带者队列进行随访是未来的研究方向，而非既定的临床推荐。