IF31.7｜主动脉瓣狭窄多组学分析揭示新治疗靶点与预测模型

大规模多祖先GWAS与瓣膜eQTL整合，发现261个风险位点并验证CMKLR1和LTBP4功能。

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

研究背景

主动脉瓣狭窄（AS）是一种常见心脏瓣膜病，全球患者超900万，但缺乏有效药物，中重度未治疗患者四年死亡率高达34%和45%。多数AS由钙化性主动脉瓣疾病（CAVD）引起，遗传度估计高达49%。既往GWAS已发现49个风险位点，主要涉及脂代谢、炎症、钙化等通路，其中LPA基因编码的脂蛋白(a)是最强遗传风险因子，相关降Lp(a)临床试验正在进行。然而，AS的遗传结构仍大部分未解释，且既往研究主要基于欧洲血统人群，缺乏瓣膜组织特异性表达调控信息。

科学问题

卡在哪里： 此前AS的GWAS样本量有限，发现的位点较少，且缺乏X染色体、性别分层和非洲血统等非欧洲人群的分析。更重要的是，风险位点如何通过主动脉瓣组织特异性基因表达影响AS尚不清楚，缺少整合瓣膜eQTL的TWAS和功能验证。

本文要回答： 通过大规模多祖先GWAS和瓣膜eQTL整合，全面解析AS遗传结构，发现新风险位点、性别和血统特异性位点，并利用TWAS和基因沉默实验鉴定潜在因果基因和通路。

技术路线

作者主要用了：多祖先GWAS、X染色体分析、性别分层、TWAS、eQTL共定位、PRS、siRNA基因沉默、免疫荧光。

作者首先汇总30个队列共285万人的基因型和AS表型，进行多祖先GWAS及X染色体、性别和血统分层分析，鉴定独立风险位点。随后利用484例人主动脉瓣RNA-seq数据构建表达预测模型，进行TWAS和eQTL共定位，筛选瓣膜表达相关基因。接着整合多种注释方法对每个位点进行因果基因优先级排序，并开发新的AS多基因风险评分（PRS）在多个独立队列验证。最后对5个优先基因在瓣膜间质细胞中进行siRNA敲低实验，检测钙化表型，并用免疫荧光验证关键基因在钙化瓣膜中的表达定位。

图 1　Fig. 1。图1展示了多祖先AS GWAS的环形曼哈顿图，内圈红线为全基因组显著性阈值（P=5×10^-8）。外圈标注了优先基因，红色为新位点，黑色为已复制位点。该图直观呈现了241个常染色体独立风险位点的染色体分布，新位点占绝大多数，显示本研究大幅扩展了AS遗传图谱。

核心亮点

亮点 1｜发现261个独立风险位点

多祖先GWAS共发现241个常染色体独立风险位点，其中203个为新位点；X染色体分析新增3个位点；欧洲血统分层新增11个，非洲血统新增1个，男性分层新增5个，总计261个独立位点，其中223个为新发现。5个位点效应较大（OR>2），多为罕见变异，如APOB p.R3527Q。遗传力估计为0.087，lead变异解释5.5%方差。这些结果大幅扩展了AS的遗传图谱。

图 2　Fig. 2。图2为喇叭图，展示lead变异的效应量和效应等位基因频率。绿色为多祖先发现，灰色为既往位点，红色为非洲血统特异，蓝色为欧洲血统特异，黄色为男性特异。高效应变异被标注最近基因。图中可见少数罕见变异具有较大OR（>2），如APOB p.R3527Q，提示罕见变异在AS中的贡献。

亮点 2｜TWAS揭示54个新基因

利用欧洲AS GWAS和人主动脉瓣eQTL进行TWAS，发现192个基因达到Bonferroni校正显著性，其中66个基因与AS风险共定位。54个基因为新发现，48个位于全基因组显著位点。8个基因具有瓣膜表达特异性（ESS>0.1），包括COMP、LTBP2、ACAN等。这些基因富集于肌动蛋白细胞骨架、TGFβ结合等通路，提示瓣膜局部生物学在AS中的关键作用。

图 3　Fig. 3。图3为TWAS的曼哈顿图，使用欧洲AS GWAS汇总统计和人主动脉瓣eQTL数据。Z分数表示基因表达与AS风险的关联方向和强度，红色为新TWAS关联。图中可见多个基因达到Bonferroni显著性，其中许多为新发现，且多数位于全基因组显著位点，表明瓣膜局部表达调控在AS中的重要性。

亮点 3｜CMKLR1和LTBP4功能验证

对5个优先基因进行siRNA敲低实验，发现沉默CMKLR1和LTBP4显著降低瓣膜间质细胞钙化（CMKLR1平均染色强度差-0.19，P=0.028；LTBP4差-0.31，P=0.049）。免疫荧光显示LTBP4定位于钙化结节和肌样细胞，CMKLR1与钙化细胞共定位并在微血管表达。这些结果支持CMKLR1（多不饱和脂肪酸信号）和LTBP4（TGFβ信号）在AS钙化中的功能作用。

图 4　Fig. 4。图4为示意图，展示至少三种方法支持的基因优先排序结果。绿色代表多祖先GWAS，蓝色代表欧洲血统。共154个多祖先lead变异、8个欧洲特异、5个男性特异、1个非洲特异和3个X染色体lead变异因支持方法少于三种而未显示。该图强调多证据整合提高了因果基因的可信度。

亮点 4｜新PRS预测性能翻倍

基于LDpred2构建的AS PRS在UKB、TIMI试验和UCLA ATLAS中验证，每标准差HR分别为1.92、1.80和1.23-1.38，较既往PRS关联强度提高约两倍。在UKB中，高遗传风险（前20%）的HR高于除年龄外的所有临床风险因素。加入PRS后C-index从0.85提升至0.87，连续NRI为29%。在非欧洲人群中也有显著关联，显示跨血统预测潜力。

图 5　Fig. 5。图5为森林图，a部分展示AS PRS在UKB、TIMI试验和UCLA ATLAS中的验证结果，每标准差HR/OR及95%置信区间，并与既往PRS比较。b部分比较UKB中遗传风险类别与个体临床风险因素的HR。结果显示新PRS关联强度约为既往的两倍，且高遗传风险HR高于除年龄外的所有临床因素。

亮点 5｜性别与血统特异性信号

性别分层分析发现IL6R位点rs7512646在女性中效应更强（女性OR=1.09，男性OR=1.03，异质性P=2.9×10^-6），提示炎症通路存在性别差异。非洲血统GWAS发现一个新位点rs148649124（LOC102723966），为长非编码RNA。LPA区域在非欧洲人群中未达全基因组显著，但rs10455872方向一致，提示跨血统证据仍有限。

生信可带走

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

• 方法栈： 多祖先GWAS、X染色体分析、性别分层、TWAS、eQTL共定位、PRS、siRNA基因沉默、免疫荧光
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事心血管遗传学、GWAS后功能验证或多基因风险评分研究的科研人员。

带走一句

大规模多祖先GWAS结合瓣膜eQTL和功能验证，可有效发现新风险基因并提升PRS预测力，为复杂疾病遗传研究提供范式。