摘要

目的 通过双向两样本孟德尔随机化(MR)分析探索血浆代谢物与冠心病(CHD)的因果关系。方法 采用双向两样本MR分析方法，正向MR分析时将血浆代谢物作为暴露因素，CHD作为研究结局；反向MR分析时则将CHD作为暴露因素，血浆代谢物作为研究结局。从GWAS Catalog数据库获取1 400种血浆代谢物，其均为欧洲人群样本，样本量为8 299例，单核苷酸多态性(SNP)数量约1 540万个。从IEU OpenGWAS project网站获取CHD数据集，其为欧洲人群样本，样本量为86 995例，SNP数量为2 420 361个。使用R 4.3.3软件、采用Two Sample MR包(版本号：0.4.22)在Windows 11系统上进行MR分析，以逆方差加权法(IVW)分析结果为主，以MR-Egger回归、加权中位数法(WME)、简单模型、加权模型分析结果为补充。采用Cochran's Q检验判断SNP间是否存在统计学异质性，采用MR-Egger回归的截距项、MR-PRESSO Outlier检验、MR-PRESSO Global检验分析SNP的水平多效性，采用留一法分析单个SNP对IVW分析结果的影响，绘制漏斗图以评估SNP的潜在偏倚。结果 共筛选出62个与血浆代谢物强相关的SNP，15个与CHD强相关的SNP。最终共纳入6种血浆代谢物。正向IVW分析结果显示：1-棕榈酰-GPG(16∶0)升高是CHD的保护因素〔OR=0.842，95%CI(0.766～0.926)〕，N-乙酰腐胺〔OR=1.125，95%CI(1.067～1.186)〕、N-乙酰天冬氨酸〔OR=1.068，95%CI(1.029～1.109)〕、5-乙酰氨基-6-甲酰氨基-3-甲基尿嘧啶〔OR=1.095，95%CI(1.048～1.144)〕、N-乙酰瓜氨酸〔OR=1.079，95%CI(1.034～1.126)〕、花生四烯酸(20∶4n6)至油酸至反式十八碳一烯酸的比率〔OR=1.316，95%CI(1.136～1.524)〕升高是CHD的危险因素(P<0.05)；反向IVW分析结果显示：CHD是1-棕榈酰-GPG(16∶0)〔OR=1.089，95%CI(1.014～1.169)〕、花生四烯酸(20∶4n6)至油酸至反式十八碳一烯酸的比率〔OR=1.094，95%CI(1.013～1.182)〕升高的危险因素(P<0.05)，CHD与N-乙酰腐胺、N-乙酰天冬氨酸、5-乙酰氨基-6-甲酰氨基-3-甲基尿嘧啶、N-乙酰瓜氨酸无因果关系(P>0.05)。MR-Egger回归、WME、简单模型、加权模型分析的β值均与正向和反向IVW的β值方向一致。Cochran's Q检验分析结果显示，与上述6种血浆代谢物、CHD强相关的SNP间不存在统计学异质性(P>0.05)。MR-Egger回归的截距项、MR-PRESSO Outlier检验、MR-PRESSO Global检验分析结果显示，与上述6种血浆代谢物、CHD强相关的SNP不存在水平多效性(P值均>0.05)。留一法分析结果显示，剔除单个SNP后，MR分析结果无明显改变。漏斗图分析结果显示，与上述6种血浆代谢物、CHD强相关的SNP基本对称分布，无潜在偏倚。结论 1-棕榈酰-GPG(16∶0)、花生四烯酸(20∶4n6)至油酸至反式十八碳一烯酸的比率与CHD互为因果关系，N-乙酰腐胺、N-乙酰天冬氨酸、5-乙酰氨基-6-甲酰氨基-3-甲基尿嘧啶、N-乙酰瓜氨酸升高是CHD的危险因素。

关键词