Omega-3脂肪酸与跟腱炎（Achilles Tendinitis）风险的因果关系：一项欧洲人群中的两样本孟德尔随机化（Two-sample Mendelian Randomization, 2SMR）研究

《Med Mat》：The causal relationship between Omega-3 fatty acids and Achilles Tendinitis risk: A two-sample Mendelian randomization study in European populations

　　

摘要：Omega-3脂肪酸（Omega-3 fatty acids）存在于鱼油、亚麻籽油和核桃中，因其抗炎特性被广泛用于多种炎症性疾病的治疗。尽管其在心血管疾病和关节炎中的获益已有充分记载，但其在肌肉骨骼疾病特别是跟腱炎（Achilles Tendinitis）中的作用尚不清楚。研究人员采用两样本孟德尔随机化（Two-sample Mendelian Randomization, 2SMR）方法探讨Omega-3脂肪酸与跟腱炎风险的因果关系。从公开全基因组关联研究（Genome-Wide Association Study, GWAS）数据中筛选与Omega-3水平相关的遗传变异，使用19个独立单核苷酸多态性（Single Nucleotide Polymorphism, SNP）作为工具变量（Instrumental Variable, IV）预测Omega-3水平。应用逆方差加权（Inverse Variance Weighting, IVW）、MR-Egger回归及加权中位数（Weighted Median）法评估基因预测Omega-3水平与跟腱炎风险的因果关系，并行留一法（Leave-one-out）、漏斗图等敏感性分析评估稳健性。主要分析显示，基因预测的Omega-3水平每升高1个标准差（standard deviation, SD），跟腱炎风险增加18%（比值比 Odds Ratio, OR＝1.18, P＝.020）。多种稳健方法得到一致结果，但仅IVW法达到名义显著性。敏感性分析证实剔除任意单个SNP结果仍稳定，无显著水平多效性（horizontal pleiotropy）或异质性（heterogeneity）。本研究提示Omega-3脂肪酸与跟腱炎风险间可能存在因果关联，尽管效应量较小且多重检验校正后无统计学意义，仍需大样本及不同Omega-3亚型分析进一步验证，以明确其在肌腱疾病中的潜在临床意义。

论文解读：Omega-3脂肪酸与跟腱炎风险的因果关系——欧洲人群两样本孟德尔随机化研究

研究背景与目的

Omega-3脂肪酸（主要为二十碳五烯酸 EPA 和二十二碳六烯酸 DHA）因抗炎特性被广泛研究于心血管病、糖尿病、关节炎等慢性炎症性疾病，既往部分证据提示其可降低关节组织炎性因子并改善软骨与肌腱健康。跟腱炎（Achilles Tendinitis / Achilles tendinopathy）是常见的运动相关过度使用损伤，慢性炎症、肿胀和疼痛严重影响患者生活质量，现有非甾体抗炎药、物理治疗及手术疗效有限且可能伴副作用。然而目前关于Omega-3与跟腱炎的关系多为观察性研究，存在反向因果及混杂偏倚，缺乏因果推断证据。为此，研究人员采用两样本孟德尔随机化（Two-sample Mendelian Randomization, 2SMR）——利用遗传变异作为工具变量（Instrumental Variable, IV）模拟随机对照试验以减少混杂与反向因果——探究基因预测的循环Omega-3脂肪酸水平与跟腱炎发病风险的因果关联，以填补该领域空白，为肌腱疾病的预防与干预提供理论依据。本文发表于《Med Mat》。

主要关键技术方法概述

研究采用2SMR设计：暴露（Exposure）数据为IEU OpenGWAS中met-c-855即UK Biobank/MetaboChip来源的欧洲人群血清/血浆Omega-3脂肪酸水平GWAS汇总数据；结局（Outcome）数据为FinnGen Release 11中跟腱炎（finngen_R11_M13_ACHILLESTEND）表型汇总数据（n＝326,188），均为欧洲裔且无样本重叠。初筛P＜5×10^?6关联的SNP，经连锁不平衡（Linkage Disequilibrium, LD）聚类（r²＜0.001，窗口10,000 kb，参照1KGP欧洲人群）保留19个独立SNP作为IV，计算F统计量确保强工具（F＞10），累计解释方差R²＝5.78%。对齐暴露与结局效应等位基因方向，缺失SNP用高LD（r²＞0.8）代理补全。主要因果效应估计采用乘法随机效应逆方差加权（Inverse Variance Weighting, IVW）法，辅以MR-Egger、加权中位数（Weighted Median）、加权模式（Weighted Mode）及简单模式（Simple Mode）作稳健性检验。异质性用Cochran Q检验，水平多效性用MR-Egger截距检验，留一法（Leave-one-out）逐一剔除SNP重估，漏斗图定性评估小研究偏倚，多重比较用Benjamini–Hochberg（BH）错误发现率（False Discovery Rate, FDR）校正，分析基于R语言two sample MR包。

研究结果

3.1 主要效应分析（Main Effect Analysis）

3.1.1 工具变量筛选与模型设置（Instrumental variable selection and model setup）：经严格质控与LD聚类获19个独立强IV（F_min＝22.50，F_mean＝92.74，R^{2＝5.78%），rs174546位于FADS基因簇解释暴露方差最大（R2＝1.14%，F＝155.93）但不主导整体估计。完成等位基因协调与代理SNP填补，预设IVW为主模型。}

3.1.2 主要分析结果——IVW法（Main analysis results (IVW Method)）：IVW分析显示基因预测Omega-3每升高1 SD，跟腱炎风险OR＝1.182（95% CI：1.026～1.362），β＝0.1676，SE＝0.0722，P＝.020，提示基因高水平Omega-3与跟腱炎风险升高约18%相关（名义显著）。BH FDR校正后Q值＝0.458未达显著，作者认为IVW名义显著及各模型方向一致仍为主要推断依据。

3.1.3 多模型一致性验证（Consistency validation across multiple models）：MR-Egger（OR＝1.366，95% CI跨1，P＝.101）、加权中位数、加权模式及简单模式均呈正向效应（β＞0），与IVW方向一致但未达名义显著，反映非IVW模型效率较低，整体正向趋势非单一模型人为产物。

3.2 可视化分析（Visualization analysis）

3.2.1 散点图（Scatter plot）：各SNP对Omega-3的效应（x轴）与对跟腱炎的效应（y轴）呈近似线性正相关，IVW与MR-Egger拟合线斜率均为正，IVW较保守，无典型扇形离散或单向偏移，少数高杠杆点（如FADS簇rs174546）未改变总体正向趋势，与低异质性及无水平多效性诊断相符。

3.3 敏感性分析（Sensitivity analysis）

3.3.1 留一法分析（Leave-one-out analysis）：依次剔除任一个SNP后重组IVW估计，合并β范围0.1277～0.2221（OR 1.14～1.25），效应方向始终为正且不反转，表明整体效应由多位点同向累积贡献而非单个离群工具驱动。

3.3.2 漏斗图（Funnel plot）：点对称分布于IVW合并效应线两侧，底部低精度区呈近似对称三角扩展，无明显不对称或空缺，MR-Egger参考线略右移但近对称，与Egger截距不显著相印证，提示方向性水平多效性不显著。

3.3.3 森林图（Forest plot）：多数SNP点估计为正，个体CI宽且跨1，FADS簇等高信息位点贡献较大但无主导离群拉动，IVW合并OR＞1且95% CI不含1，MR-Egger合并OR＞1但CI跨1，与前述一致。

3.3.4 MR-Egger截距检验（MR-Egger intercept test）：截距＝?0.0166（SE＝0.0190，P＝.3928），95% CI含0，未检测到显著方向性水平多效性。

3.3.5 Cochran Q异质性检验（Cochran Q heterogeneity test）：IVW Q＝19.754（df＝18，P＝.3468），MR-Egger Q＝18.900（df＝17，P＝.3343），均无显著异质性，支持工具变量效应一致性。

综合上述，19个强IV、低异质性、无显著水平多效性前提下，IVW提示每1 SD Omega-3升高对应跟腱炎OR＝1.18（P＝.020），多模型方向一致，敏感性分析稳健，但效应量小且FDR校正不显著，需更大样本及亚型分层验证。

讨论与结论总结

讨论指出本研究满足MR三基本假设：相关性（F＞10，弱工具偏倚低）、独立性（暴露与结局源自独立欧洲人群平台无重叠、多方法诊断无系统混杂模式）、排他性（Egger截距不显著、漏斗图对称、Cochran Q不显著、留一法排除单SNP驱动），内部效度较好。基因因果层面观察到基因升高Omega-3水平与跟腱炎风险呈微弱但方向一致、内部稳健的名义正相关，因未通过多重检验校正且未区分EPA/DHA/α-亚麻酸（ALA）亚型及通路分解，应视为需进一步验证的因果线索，其科学价值在于提示特定脂肪酸亚型、炎症—修复网络与肌腱退变间潜在联系，为后续机制研究与干预试验提供路线图。局限含仅适用于欧洲裔、未调整生活方式/运动/外伤史等，未来需更大样本、亚型分层、多变量MR（Multivariable MR, MVMR）纳入Omega-6/炎性标志物/血脂及FADS区共定位（colocalization）分析。

结论（翻译）

在欧洲裔人群中，研究人员用2SMR评估基因预测Omega-3水平（IEU: met-c-855）与跟腱炎（FinnGen R11: finngen_R11_M13_ACHILLESTEND）的因果关系。主要乘法随机效应IVW分析表明，基因决定的Omega-3每升高1 SD与跟腱炎风险名义上升高相关（OR＝1.18，P＝.020）。MR-Egger、加权中位数、加权模式及简单模式效应方向均一致为正；且Egger截距不显著、Cochran Q无显著异质性、漏斗图近似对称、留一法剔除任一变异不改变效应方向，支持结果内部稳健。BH FDR校正未达显著，提示若存在真实因果效应则其量小，需更多数据确认。综上，在遗传因果层面观察到基因高水平Omega-3与跟腱炎风险呈微弱、方向一致且内部稳健的正向关联；鉴于多重检验校正后不显著及尚未完成亚型与通路分解，该发现为可检验的因果假说与研究路线图，可通过更大样本、亚型分层、多变量MR及共定位分析进一步明确脂肪酸谱、炎症/修复网络与肌腱退变的关系。