IF31.7｜全基因组关联分析揭示扩张型心肌病分子病因

GWAS+单核转录组+罕见变异分析，锁定62个效应基因并阐明多基因背景对DCM风险的影响。

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

研究背景

扩张型心肌病（DCM）以左心室收缩功能障碍和扩张为特征，是心力衰竭和心脏移植的主要原因之一，患病率约1/250。传统观点认为DCM主要由单基因罕见变异引起，但近期证据表明常见多基因背景对DCM风险有直接和间接影响。然而，DCM的复杂遗传结构尚未被系统解析，尤其是常见变异如何与罕见变异协同作用、哪些细胞类型和通路驱动疾病发生，仍不清楚。

科学问题

卡在哪里： 既往DCM的GWAS样本量有限，发现的基因组风险位点较少；罕见变异研究多聚焦于已知孟德尔基因，对GWAS位点内效应基因的系统筛选不足；多基因背景对罕见变异携带者外显率的修饰作用缺乏量化证据。

本文要回答： 通过大规模GWAS和多性状分析鉴定DCM基因组风险位点，系统优先排序效应基因，并利用单核转录组和罕见变异分析阐明疾病分子机制。

技术路线

作者主要用了：GWAS meta分析、MTAG、精细定位、基因优先排序、单核RNA-seq、罕见变异负荷分析、PGS。

作者首先对16个队列的14,256例DCM和1,199,156例对照进行GWAS meta分析，并通过与左心室CMR性状的遗传相关性进行MTAG多性状分析以增加统计功效。随后，对80个风险位点进行功能精细定位和效应基因优先排序，整合八种预测方法。进一步利用单核RNA-seq数据解析细胞类型和通路，并通过罕见变异负荷分析验证候选基因。最后构建多基因评分，评估其在普通人群中的预测能力及对罕见变异携带者外显率的修饰作用。

图 1　Fig. 1。图1展示研究设计流程：从DCM表型定义（Narrow和Broad）出发，进行GWAS meta分析，评估与CMR性状的遗传相关性，纳入MTAG，下游包括基因优先排序、单核转录组、罕见变异分析和多基因评分。

核心亮点

亮点 1｜GWAS+MTAG鉴定80个DCM风险位点

GWAS meta分析发现26个位点达到全基因组显著（P<5×10-8），另有36个位点达到1% FDR。通过与左心室收缩末期容积、整体环向应变和射血分数三个CMR性状进行MTAG，额外发现18个位点，共80个位点用于下游分析。其中31个位点此前未报道与DCM或相关心脏性状关联。遗传力估计显示，严格定义DCM的h2SNP为20%，宽松定义为11%。

图 2　Fig. 2。图2曼哈顿图显示DCM GWAS和MTAG的关联信号。蓝色虚线为全基因组显著阈值（P<5×10-8），浅蓝色虚线为FDR 1%阈值。标注了达到阈值的位点及最近蛋白编码基因，共80个位点用于下游分析。

亮点 2｜系统优先排序锁定62个效应基因

对80个位点内的1,970个基因，结合最近基因、PoPS、V2G、编码变异、eQTL共定位、TWAS、ABC模型和孟德尔心肌病基因证据，为62个位点各确定一个高置信度效应基因。优先排序靠前的基因包括MYPN、HSPB8、ALPK3、ACTN2、SPATS2L和BAG3。通路富集显示这些基因富集于肌节和细胞骨架功能、细胞黏附、自噬以及Wnt和TGFβ信号通路。

图 3　Fig. 3。图3展示基因优先排序流程和通路富集结果。a：从1,970个基因中筛选380个候选基因，最终62个位点各确定一个优先基因；b：优先基因富集于肌节结构、细胞黏附、TGFβ和Wnt信号等通路。

亮点 3｜罕见变异分析发现新候选DCM基因

在UKB和GeL队列中对优先排序基因进行罕见蛋白截断变异（PTV）负荷分析，发现MAP3K7、NEDD4L和SSPN与DCM或左心室定量性状显著相关。MAP3K7 PTV在GeL中与DCM关联（OR=24.2，Padj=0.02），并在UKB中与左心室容积增加相关；NEDD4L PTV与DCM关联（OR=10.4，Padj=0.01）并在GeL中复制；SSPN PTV与左心室射血分数降低和收缩末期容积增加相关，且在局部DCM队列中富集。

图 4　Fig. 4。图4展示DCM遗传结构：常见变异（GWAS sentinel和PGS）与罕见PTV的效应量对比。罕见PTV在孟德尔DCM基因中的效应远大于常见变异，且MAP3K7、NEDD4L和SSPN的PTV与DCM或左心室性状关联。

亮点 4｜单核转录组揭示非心肌细胞关键作用

利用52例DCM和18例对照的左心室单核RNA-seq数据，发现心肌细胞、成纤维细胞和壁细胞（主要是周细胞）的疾病特异性差异表达基因富集遗传力。细胞间通讯分析显示DCM中整体信号增强，但心肌细胞-心肌细胞相互作用减弱；TGFβ通路相关配体-受体相互作用增加，COL4A1从成纤维细胞向心肌细胞的信号增强，BMP6-BMPR1A信号从心内膜细胞向心肌细胞和成纤维细胞增强。

图 5　Fig. 5。图5展示单核转录组分析结果：a-b：组织水平和细胞类型水平的遗传力富集，心肌细胞和壁细胞显著；c：优先基因在对照心肌细胞中高表达；d：SSPN、MAP3K7和NEDD4L在DCM心肌细胞中差异表达；e-f：COL4A1和BMP6信号增强。

亮点 5｜多基因评分预测风险并修饰外显率

基于541,841个SNP构建的多基因评分（PGSDCM）在UKB中与DCM显著关联（OR=1.76/标准差，AUROC=0.71）。最高百分位个体风险是最低百分位的7倍。在1,546名罕见致病变异携带者中，PGSDCM最高五分位DCM患病率为7.3%，最低五分位为1.7%（P=0.005），表明多基因背景修饰罕见变异的外显率。

生信可带走

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

• 方法栈： GWAS meta分析、MTAG、精细定位、基因优先排序、单核RNA-seq、罕见变异负荷分析、PGS
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事心血管遗传学、GWAS下游功能解析、单细胞转录组分析或罕见变异研究的生信和遗传学研究者。

带走一句

GWAS位点内的效应基因系统优先排序结合单核转录组和罕见变异分析，可高效发现新的疾病基因和机制；多基因背景对罕见变异外显率的修饰提示遗传咨询需整合常见变异信息。