IF31.7｜人类代谢遗传图谱：跨越等位基因频率谱

45万人、249种代谢表型，揭示代谢调控的罕见与常见变异全景。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

人类代谢调控的遗传基础研究长期受限于样本量和等位基因覆盖范围。近年来，高通量代谢组学与大规模基因组学的结合，使全基因组代谢组关联研究（mGWAS）成为可能，已发现众多与循环代谢物水平相关的遗传位点，并部分验证了已知生化通路。然而，既往研究多聚焦于常见变异，罕见变异对代谢表型的贡献尚缺乏系统评估。此外，代谢物高度多效性，同一遗传变异可影响多种代谢物，但多效性的模式及其对下游疾病风险的因果意义仍不清楚。

科学问题

卡在哪里： 现有mGWAS研究存在几个关键局限：一是每个代谢物关联的遗传位点数量有限，难以支撑可靠的因果推断；二是无法区分多效性变异是独立作用于不同代谢物（水平多效性）还是作为上游“根因”引发下游级联变化（垂直多效性）；三是在疾病风险位点共定位时，难以区分位点特异效应与代谢物丰度效应；四是新发现位点的效应基因注释置信度不足。

本文要回答： 本研究旨在整合常见与罕见遗传变异，系统绘制人类循环代谢物的遗传调控图谱，解析多效性模式，并评估代谢物与心血管疾病的因果关联。

技术路线

作者主要用了：NMR代谢组、GWAS、ExWAS、精细定位、机器学习效应基因注释、共定位、孟德尔随机化。

研究基于英国生物银行（UKB）约45万人的血浆样本，利用质子核磁共振（1H NMR）技术定量249种代谢表型，涵盖脂蛋白亚类、脂肪酸、氨基酸等。作者首先在三个主要祖先群体（欧洲、非洲、中南亚）中分别进行全基因组关联分析（GWAS），随后进行跨祖先荟萃分析以发现常见变异关联。同时，利用全外显子组测序数据进行罕见变异关联分析（ExWAS），包括基因负荷检验和单变异检验。为精细定位因果变异，采用SuSiE进行统计精细定位，并利用多祖先连锁不平衡差异进一步缩小可信集。效应基因注释通过机器学习模型整合功能基因组学资源与通路信息实现。最后，通过共定位分析和孟德尔随机化评估代谢物与心血管疾病的因果关联。

图 1　Fig. 1。图1展示了跨祖先荟萃分析的曼哈顿图，显示249个代谢表型的哨兵变异在染色体上的分布，不同颜色代表生化类别。b图比较了哨兵变异的等位基因频率与效应量，显示常见变异效应量较小，罕见变异效应量较大。c图显示效应量在三个祖先群体间高度一致，支持跨祖先遗传调控的保守性。

核心亮点

亮点 1｜跨祖先29,824个代谢关联

跨祖先荟萃分析在249个代谢表型中识别出29,824个区域哨兵变异-代谢物关联，映射到753个非重叠基因组区域。其中近半数区域（359个，47%）与超过10种代谢物相关，显示广泛多效性。大HDL颗粒特征关联区域最多（中位数166个），远超所有代谢物的中位数105个。效应量在三个祖先群体间高度一致，仅1.14%的位点显示显著异质性。非洲裔特异的CD36错义变异rs3211938与较低的血浆omega-3脂肪酸水平及15种其他代谢物显著相关，验证了CD36作为脂肪酸转运蛋白的功能。

图 2　Fig. 2。图2a比较了性别差异位点的效应量，显示两性间效应方向一致，但存在数量差异。b图展示了与临床脂质性状相比，NMR脂蛋白亚类测量在遗传关联中的排名分布，约25%的变异在对应NMR测量中排名前10%，表明精细脂蛋白谱能提供更特异的生物学信息。

亮点 2｜机器学习注释效应基因

通过训练随机森林分类器整合eQTL、代谢通路、OMIM和药物靶点等特征，成功为73.6%的欧洲裔精细定位mQTL（2,213个）分配了至少中等置信度的效应基因，其中28.2%为高置信度。例如，ELOVL6被优先为16种VLDL/HDL特征的效应基因，其编码的脂肪酸延长酶催化长链脂肪酸延长的限速步骤。在17q25.3位点，CYTH1被优先为5个独立变异、11种代谢物的效应基因，该基因通过激活ARF蛋白调节脂质囊泡运输，与脂蛋白代谢间接相关。

图 3　Fig. 3。图3通过散点图展示四种多效性模式：a通路多效性、b比例多效性、c不成比例多效性、d非特异性多效性。e图统计了各模式的数量，f图显示多效性模式与GWAS Catalog中表型多效性的关联，g图显示比例和非特异性多效性变异在代谢疾病、脂肪肝和动脉疾病等类别中显著富集。

亮点 3｜罕见变异揭示新调控基因

罕见变异分析（MAF≤0.05%）在209个基因中发现与代谢物显著关联，包括3,709个基因负荷关联和4,131个单变异关联。效应量显著大于常见变异，例如SLC13A5罕见功能缺失变异携带者血浆柠檬酸浓度升高1.41个标准差。SIDT2罕见变异与124种代谢物相关，主要影响HDL颗粒特征，支持其在肝脏脂质代谢和ApoA1分泌中的作用。JAK2 p.Val617Phe变异与73种代谢物关联，其与HDL参数的强逆相关提示JAK2可能通过ABCA1促进HDL介导的脂质清除。

图 4　Fig. 4。图4a比较了基因负荷和罕见外显子变异的效应量与MAF，显示罕见变异效应量更大。b图展示SIDT2罕见内含子变异与249种代谢物的关联，多数显著。c图显示罕见变异与常见可信集变异的基因组距离分布，d图展示等位基因系列证据，e和f图展示APOA1等位基因系列中不同变异对HDL颗粒特征的差异影响。

亮点 4｜APOA1等位基因系列异质性

在APOA1基因内发现由7个变异（5个罕见、2个常见）和罕见功能缺失负荷组成的等位基因系列，与87种代谢物差异关联。四个罕见错义变异（MAF≤0.03%）对HDL颗粒数量、大小和胆固醇含量的影响模式不同：p.Leu158Pro主要降低血清ApoA1浓度和HDL颗粒数，模拟功能缺失效应；p.Lys131del和p.Arg201Ser导致胆固醇从大HDL向小HDL转移，而p.Asp113Glu则相反。这种表型异质性提示同一基因内不同变异可能通过不同机制影响代谢。

图 5　Fig. 5。图5a总结了两样本孟德尔随机化分析结果，显示NMR代谢物与多种心血管疾病的因果效应估计。b图展示位点-疾病网络，突出显示与至少一种心血管疾病共定位的位点。c图展示HDL粒径与冠心病风险的关联，e图展示VEGFA位点的共定位证据。

亮点 5｜HDL粒径与冠心病风险

孟德尔随机化分析发现，遗传预测的较大HDL粒径与冠心病风险降低35%相关（OR=0.65，95% CI 0.50-0.83，Padj<0.007）。共定位分析显示VEGFA位点（rs4711750，PP=99%）与HDL粒径和冠心病风险共享遗传信号。VEGFA编码血管内皮生长因子A，已知参与血管生成，但本研究提示其可能通过激活SR-BI促进HDL胆固醇的跨内皮转运，从而影响冠心病风险。这一发现为超越LDL胆固醇降低的残余心血管风险管理提供了潜在靶点。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： NMR代谢组、GWAS、ExWAS、精细定位、机器学习效应基因注释、共定位、孟德尔随机化
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事代谢遗传学、心血管遗传学或药物靶点发现的研究者，尤其是关注多效性解析和罕见变异分析的生信人员。

带走一句

大规模代谢组遗传研究需整合常见与罕见变异，并系统解析多效性模式，才能准确识别因果基因和疾病相关通路。