IF31.7｜大规模GWAS揭示扩张型心肌病遗传机制与心肌韧性

70个新位点、63个优先基因，揭示心肌细胞收缩装置在DCM中的核心作用。

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

研究背景

扩张型心肌病（DCM）以左心室扩大和收缩功能下降为特征，是心衰和猝死的重要原因。约25%的先证者可检出罕见致病变异，但多数病例缺乏已知的单基因病因。近年来，全基因组关联研究（GWAS）已证实DCM存在多基因遗传基础，但既往研究样本量有限，仅发现少数显著位点，对疾病机制的理解仍不充分。心肌细胞功能障碍被认为是DCM的核心，但具体致病基因和通路尚未完全阐明。此外，DCM与收缩性心衰的关系、常见变异与罕见变异的相互作用、以及可干预的危险因素等问题均有待深入探索。

科学问题

卡在哪里： 既往DCM的GWAS样本量小，发现的显著位点少，统计功效不足，难以全面揭示遗传结构。细胞类型特异性富集分析结果存在争议，不同研究对心肌细胞和其他细胞类型的贡献结论不一。此外，DCM的因果危险因素、多基因风险评分在不同人群和临床场景中的预测价值，以及常见变异与罕见致病变异的叠加效应均缺乏系统评估。

本文要回答： 通过大规模GWAS和MTAG分析，系统鉴定DCM相关遗传位点和基因，阐明关键细胞类型和通路，评估多基因风险评分的预测性能，并利用孟德尔随机化识别可干预的因果危险因素。

技术路线

作者主要用了：GWAS、MTAG、LDSC、eQTL/pQTL共定位、PoPS、单细胞RNA-seq、孟德尔随机化、PRS。

作者首先整合六个队列（三个临床DCM病例对照和三个生物样本库）进行GWAS荟萃分析，得到GWAS-DCM。随后利用与DCM强遗传相关的两个心脏MRI性状（整体周向应变Ecc和左室收缩末期容积指数LVESVi）进行MTAG分析，以提升发现功效。接着通过连锁不平衡评分回归（LDSC）评估遗传力和细胞类型富集，利用多种方法（包括PoPS、eQTL/pQTL共定位、蛋白质改变变异等）进行变异到基因的映射和优先排序。进一步整合三个已发表的单细胞RNA测序数据集，分析优先基因在左室各细胞类型中的表达模式。最后，通过双向孟德尔随机化筛选DCM的因果危险因素和后果，并在多个独立队列中构建和验证多基因风险评分（PRS）。

图 1　Fig. 1。图1展示了研究设计：a部分显示GWAS在生物样本库队列（NICM和NI-DCM）和临床DCM队列中进行，随后将NI-DCM和临床DCM结果进行荟萃分析得到GWAS-DCM，再与心脏MRI性状（Ecc和LVESVi）结合进行MTAG。b部分列出下游分析，包括组织/细胞类型富集、基因优先排序、单细胞表达、孟德尔随机化和PRS分析。该图概括了从遗传发现到功能解读的完整流程。

核心亮点

亮点 1｜GWAS与MTAG发现70个显著位点

GWAS-DCM纳入9,365例病例和946,368例对照，发现38个全基因组显著位点，其中27个为新位点。MTAG-DCM进一步发现65个显著位点，其中50个为新位点。综合两者共鉴定出70个全基因组显著位点。在独立样本（HERMES、MVP、All of Us）中，GWAS-DCM位点92%复制，MTAG-DCM位点88%复制，证实了结果的稳健性。这些位点广泛涉及心肌细胞功能和收缩装置相关基因，为DCM的遗传结构提供了全面视图。

图 2　Fig. 2。图2为迈阿密图，上方为GWAS-DCM结果（9,365例/946,368例对照），下方为MTAG-DCM结果。横轴为基因组位置，纵轴为-log10(P值)。GWAS-DCM显示38个显著位点，MTAG-DCM显示65个显著位点。图中可见BAG3、HSPB7等已知位点以及多个新位点。MTAG显著提升了发现功效，新增大量位点。

亮点 2｜心肌细胞是DCM遗传力的主要贡献细胞类型

利用两个已发表的左室单核RNA测序数据集进行分层LDSC分析，发现仅心肌细胞特异性基因显著富集DCM遗传力（P < 3×10^-7），而其他细胞类型未显示稳健富集。这一结果强调了心肌细胞功能障碍在DCM发病中的核心地位，与既往部分研究结论不同，可能源于统计方法差异。该发现为后续基因优先排序和机制研究提供了明确的细胞类型先验。

图 3　Fig. 3。图3为气泡热图，展示63个优先基因在三个已发表的DCM和对照左室单细胞RNA测序数据集中的表达。横轴为协调后的左室细胞类型，纵轴为优先基因。点的大小代表表达该基因的细胞核/细胞百分比，颜色代表相对表达水平。许多基因在心肌细胞中高表达，如TTN、OBSCN、ACTN2等，支持收缩装置在DCM中的核心作用。

亮点 3｜63个优先基因揭示收缩装置与热休克通路

通过整合PoPS、eQTL/pQTL共定位、蛋白质改变变异、最近基因和V2F五种方法，对GWAS-DCM和MTAG-DCM位点进行基因优先排序，最终得到63个独特优先基因。这些基因包括已知的孟德尔心肌病基因（如BAG3、TTN、FHOD3）以及心脏中功能未知或较少描述的基因（如CRIM1、MLF1、HSPA4、ERBB4等）。基因集富集分析显示‘细胞热应激反应’通路显著富集，此外还涉及心肌发育、ERBB信号和细胞骨架组织等通路。单细胞表达分析表明许多优先基因在心肌细胞中高表达或优先表达，包括结构基因SVIL和PDLIM5，以及转录因子MITF和MLIP。

图 4　Fig. 4。图4展示孟德尔随机化筛选结果。a部分为气泡图，左面板为DCM作为结局，右面板为DCM作为暴露。纵轴为带符号的-log10(P值)，气泡代表不同疾病/性状。红色和蓝色线表示Bonferroni校正阈值。b部分为MR散点图，显示DCM易感性对心衰风险的影响，37个遗传工具变量，黑色线为加权中位数法，虚线为CAUSE方法。

亮点 4｜孟德尔随机化支持体重和收缩压为因果危险因素

双向孟德尔随机化筛选73种疾病和数量性状，发现体重、BMI、房颤、收缩压和身高为DCM的潜在因果危险因素，而心衰和平均血小板体积为DCM的潜在后果。多变量MR显示体重、收缩压和房颤为独立危险因素。CAUSE分析进一步支持BMI、体重和收缩压增加DCM风险。这些结果提示早期血压控制和体重管理可能是DCM预防和管理的合理策略。

图 5　Fig. 5。图5展示PRS分析结果。a部分为MR散点图，x轴为GWAS-DCM中37个遗传工具变量的beta系数，y轴为对一般心衰的beta系数，黑色线为WM方法，虚线为CAUSE。b部分为森林图，显示DCM PRS在All of Us数据集中与NI-DCM、收缩性心衰及去除DCM/NICM后的收缩性心衰的关联，均显著。

亮点 5｜多基因风险评分预测DCM及收缩性心衰

基于GWAS-DCM和MTAG-DCM构建的多基因风险评分（PRS）在All of Us、Amsterdam UMC和UKB数据集中均显著关联DCM。MTAG-DCM PRS表现最佳，在欧洲裔（OR=1.73）和非洲裔（OR=1.61）中均显示强关联。在Amsterdam队列中，携带罕见致病变异的DCM患者PRS显著低于非携带者，但两组PRS均高于对照，表明常见变异在有无罕见变异的情况下均贡献疾病风险。此外，DCM PRS与收缩性心衰显著相关（OR=1.30），且在去除DCM/NICM病例后仍显著，提示其捕获了内在心肌功能障碍的易感性。

生信可带走

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

• 方法栈： GWAS、MTAG、LDSC、eQTL/pQTL共定位、PoPS、单细胞RNA-seq、孟德尔随机化、PRS
• 公开数据： GSE183852
• 代码： 公开仓库 https://github.com/FinucaneLab/pops、https://github.com/JonJala/mtag、https://github.com/OliverVoogd/edgeR；本地已克隆：DESeq2、cause、seurat、METAL、ldsc、regenie、coloc、edgeR、mtag、PRScs、MR_pipeline_sjj
• 谁该点开原文： 从事心血管遗传学、GWAS后功能解读、多基因风险评分或单细胞组学整合分析的研究者。

带走一句

大规模GWAS结合MTAG和单细胞数据，能有效揭示复杂心脏病的遗传机制；PRS在跨种族和临床场景中具有预测价值，提示多基因风险评分在心肌病精准管理中的潜力。