IF31.7｜罕见杂合致病变异对基因功能与人类性状的梯度效应

代谢组+外显子组测序，揭示罕见杂合变异如何模拟纯合敲除效应并影响身高与骨骼肌肉疾病。

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

研究背景

人类代谢由数千种酶和转运蛋白精密调控，其功能紊乱可导致先天性代谢缺陷（IEM）或常见代谢病。传统IEM研究受限于纯合突变患者极罕见，而全基因组关联研究（GWAS）发现的常见变异效应微弱且功能不明。基因水平的罕见变异聚合检验可在群体中评估罕见、有害变异的累积效应，但既往研究几乎仅关注循环代谢物，且未系统探索杂合变异与纯合敲除效应之间的关系。

科学问题

卡在哪里： 此前罕见变异代谢组研究主要局限于血浆代谢物，未充分利用尿液代谢组信息；且缺乏系统证据表明杂合有害变异能否模拟纯合敲除的代谢效应，以及代谢物水平能否作为人类复杂性状的中间表型。

本文要回答： 通过整合血浆和尿液代谢组与外显子组测序，发现新的基因-代谢物关联，并验证杂合变异能否提供与IEM一致的基因功能信息及剂量效应。

技术路线

作者主要用了：全外显子测序、非靶向代谢组、基因负担检验、全身体代谢模型、UK Biobank PheWAS。

作者利用德国慢性肾脏病（GCKD）研究中4737名参与者的全外显子组测序数据，对1294种血浆和1396种尿液代谢物进行基因水平的罕见变异聚合检验，使用两种互补的变异筛选掩模（LoF_mis和HI_mis）。随后通过UK Biobank独立复制、蛋白质组数据、计算建模（全身体代谢模型的虚拟IEM）以及体外功能实验进行多维度验证，并进一步在UK Biobank中探索代谢相关基因与人类性状和疾病的关联。

图 1　Fig. 1。图1展示研究整体设计：从GCKD队列获取WES和血浆/尿液代谢组数据，进行基因水平罕见变异聚合检验，随后通过UKB复制、蛋白质组、计算建模和实验验证，并关联人类性状和疾病。

核心亮点

亮点 1｜发现192个基因-代谢物关联

在血浆和尿液代谢物中分别以P<5.04×10⁻⁹和P<4.46×10⁻⁹为阈值，共发现192个显著基因-代谢物关联（涉及73个基因和179种代谢物），其中43个在血浆和尿液中均显著。氨基酸和脂质是主要代谢通路。与既往研究比较，69%的血浆关联为新发现。73个基因中38%已知为IEM致病基因，显著富集（OR=10.6，P=1.9×10⁻¹⁴）。

图 2　Fig. 2。图2显示全基因组范围内基因-代谢物关联的环形图。外圈为血浆代谢物关联（红色），内圈为尿液代谢物关联（蓝色），灰色标签为既往报道，黑色为新发现。氨基酸和脂质是主要代谢通路，血浆中脂质关联比例更高。

亮点 2｜多维度验证支持关联真实性

在UK Biobank中，对两个重叠代谢物（组氨酸和苯丙氨酸）使用相同QVs进行基因负担检验，效应量与GCKD高度一致。UKB蛋白质组数据显示，17个关联基因中15个的蛋白截断/罕见有害变异与循环蛋白水平降低相关（P<1×10⁻⁵），提示无义介导的mRNA降解。与既往血浆代谢组研究比较，变异水平和聚合水平的效应量相关系数>0.8。体外实验证实SLC6A19转运蛋白可转运新底物甲硫氨酸砜，且被抑制剂cinromide完全阻断。

图 3　Fig. 3。图3展示UKB中相同QVs的基因负担检验效应量与GCKD高度一致（a），UKB蛋白质组中15/17个关联基因的蛋白截断变异与循环蛋白水平降低相关（b），以及与既往血浆代谢组研究的效应量高度相关（c,d）。

亮点 3｜杂合变异反映剂量效应

X染色体基因TMLHE的男性半合子携带者尿液N6,N6,N6-三甲基赖氨酸水平显著高于女性杂合子，而产物羟基-N6,N6,N6-三甲基赖氨酸水平显著降低，提示半合子相当于纯合敲除。全身体代谢模型虚拟敲除24个基因，对60个基因-代谢物对的效应方向预测准确率达73.3%（男性）和76.7%（女性），显著优于随机。模型经人工修正后准确率提升至79.1%和83.58%。

图 4　Fig. 4。图4显示SLC6A19过表达CHO细胞对甲硫氨酸砜的转运活性呈浓度依赖性增加，与已知底物异亮氨酸相似（a），且被SLC6A19抑制剂cinromide完全阻断（b），证实甲硫氨酸砜是SLC6A19的新底物。

亮点 4｜个性化模型预测IEM生物标志物

以KYNU基因为例，利用569个微生物组个性化全身体模型进行in silico敲除，发现18种代谢物为潜在生物标志物，其中黄尿酸和3-羟基犬尿氨酸为已知标志物，而8-甲氧基犬尿氨酸为新发现。在1例纯合KYNU突变患者及其父母的尿液中，8-甲氧基犬尿氨酸水平显著升高，证实其为该IEM的新生物标志物。

图 5　Fig. 5。图5显示TMLHE基因QVs携带者尿液N6,N6,N6-三甲基赖氨酸水平：男性半合子携带者显著高于非携带者和女性杂合子，而产物羟基-N6,N6,N6-三甲基赖氨酸水平显著降低，支持剂量效应。

亮点 5｜硫酸盐转运变异与身高及骨骼肌肉疾病

SLC13A1和SLC26A1编码肾小管跨细胞硫酸盐重吸收的关键转运蛋白。6个功能性驱动变异在GCKD中与血浆硫酸盐水平显著相关，效应量在UKB中与坐高和站高呈正相关（r=0.57和0.70）。在UKB中，这些变异与多种骨折和肌肉骨骼疾病风险增加相关（如NaS1 p.Trp48*与颈部闭合性骨折OR=30.7，P=2.1×10⁻⁸）。携带多个变异或纯合子个体身高z-score呈剂量依赖性降低。

生信可带走

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

• 方法栈： 全外显子测序、非靶向代谢组、基因负担检验、全身体代谢模型、UK Biobank PheWAS
• 公开数据： 原文未给出公开组学登录号
• 代码： 公开仓库 https://github.com/DavisBrian/seqMeta、https://github.com/SysPsyHertel/CodeBase、https://github.com/VirtualMetabolicHuman/AGORA；本地已克隆：seqMeta、CodeBase
• 谁该点开原文： 从事代谢组学、罕见变异关联分析、复杂性状遗传学及药物靶点发现的研究者。

带走一句

罕见杂合致病变异可模拟纯合敲除的代谢效应，代谢物水平可作为连接基因型与复杂性状的中间表型，为药物靶点发现和疾病机制研究提供新思路。