《Current Issues in Molecular Biology》：Integrative Multi-Omics Mendelian Randomization and Single-Cell Sequencing Prioritize Therapeutic Targets for Coronary Artery Disease
编辑推荐：
冠状动脉疾病（CAD）仍是全球重大健康负担，靶向治疗策略有限。研究人员将基于血液表达数量性状基因座（eQTLs）的两样本孟德尔随机化（MR）与贝叶斯共定位、外部GWAS验证、meta分析、单细胞RNA测序、细胞–细胞通讯分析和全表型组MR（Phe-MR）相整合
冠状动脉疾病（CAD）仍是全球重大健康负担，靶向治疗策略有限。研究人员将基于血液表达数量性状基因座（eQTLs）的两样本孟德尔随机化（MR）与贝叶斯共定位、外部GWAS验证、meta分析、单细胞RNA测序、细胞–细胞通讯分析和全表型组MR（Phe-MR）相整合。在2460个具有显著cis-eQTL的基因中，MR识别出26个靶点，其中6个（PTPN11、FES、VAMP8、LIPA、GGCX和ERP29）在两个队列及meta分析中均获重复。共定位支持4个（LIPA、GGCX、FES和VAMP8；两个队列中共享因果变异的后验概率（PP.H4）≥0.80；Tier A），而ERP29（两个队列均H3主导）和PTPN11（发现队列H3主导；外部队列不确定）未达到该阈值（Tier B，MR重复但无共定位支持）。LIPA、GGCX和ERP29与CAD风险正相关（比值比（ORs）1.05、1.12、1.28）；PTPN11、FES和VAMP8呈负相关（ORs 0.83、0.90、0.90）。这些基因在CAD组织中不同细胞簇间显示差异表达。覆盖4698个基因–表型组合的Phe-MR识别出24个Bonferroni显著和132个错误发现率（FDR）显著关联，作为潜在安全性信号报告。缺乏遗传关联并不能确立药物安全性，终生遗传预测表达也不能估计短暂药物暴露。本研究在两个明显不同的证据层级上为CAD的六个候选治疗靶点提供了遗传学证据。
研究背景：冠状动脉疾病（coronary artery disease, CAD；冠状动脉粥样硬化导致心肌缺血的疾病）仍是全球死亡和残疾的主要原因。现有药物治疗主要针对肾素–血管紧张素系统、交感神经系统、胆固醇代谢和抗血小板机制，但血运重建相关并发症及多因素发病机制使新靶点需求突出。全基因组关联研究（GWAS）、表达数量性状基因座（eQTL；与基因表达相关的遗传变异）和孟德尔随机化（Mendelian randomization, MR；以遗传变异为工具变量推断暴露与结局因果关联的方法）为靶点发现提供路径，单细胞RNA测序（single-cell RNA sequencing, scRNA-seq；解析细胞异质性和表达谱的技术）可补充细胞背景。研究人员据此整合遗传因果、共定位、细胞分布和安全性筛查，以优先识别CAD治疗靶点。
研究人员开展了两样本MR、外部GWAS重复、meta分析、贝叶斯共定位（Bayesian colocalization；评估暴露与结局是否共享因果变异的分析）、组织特异性MR、scRNA-seq和全表型组MR（phenome-wide MR, Phe-MR；系统筛查基因表达与多种表型的关联）。结论为六个候选靶点分属两个证据层级：LIPA、GGCX、FES和VAMP8为Tier A（两队列MR重复且共定位PP.H4≥0.80）；ERP29和PTPN11为Tier B（两队列MR重复但共定位未达阈值）。LIPA、GGCX和ERP29与CAD风险正相关，PTPN11、FES和VAMP8负相关。研究意义在于提供遗传学优先排序、细胞类型背景和潜在安全性信号，但原文强调其确立的是遗传关联和优先排序，而非药理调控的机制或因果性。
主要关键技术方法：基于eQTLGen全血cis-eQTL（31,684例欧洲血统样本）及DGIdb与Finan et al.候选基因集，采用两样本MR（Wald ratio/逆方差加权（inverse variance weighted, IVW））在FinnGen R11（51,098例/402,635对照）和CARDIoGRAMplusC4D（60,801例/123,504对照）中分析，并以meta分析合并；用贝叶斯共定位、Steiger过滤、MR-Egger、Cochran’s Q和MR-PRESSO评估敏感性；用GSE233870冠状动脉粥样硬化scRNA-seq（6例粥样硬化样本，4例对照被排除）进行分群、差异表达和CellPhoneDB通讯分析；用GTEx主动脉与皮下脂肪eQTL做组织特异MR；用UK Biobank 783个表型（n≤408,961，≥500例）做Phe-MR。
研究结果：
3.1 利用eQTL数据识别CAD潜在治疗基因：研究人员从DGIdb和Finan et al.获得4488个候选基因，与eQTLGen全血cis-eQTL交集并以p<5×10
?8、连锁不平衡（linkage disequilibrium, LD）r
2<0.03筛选后，2460个基因有显著cis-eQTL，2446个经LD clumping保留；2418个可对FinnGen R11 I9_CHD分析、2420个可对CARDIoGRAMplusC4D分析。MR识别26个Bonferroni显著基因，其中12个增加CAD风险、14个降低风险；风险增加的前三位为IL20RB、LRP8、PSMB10，风险降低的前三位为AGRP、PTPN11、GGT7。结论为这些基因表达与CAD风险存在遗传学关联，并作为后续验证候选。
3.2 潜在治疗基因的敏感性分析：对26个基因进行共定位和敏感性分析。以PP.H4≥0.80为标准，发现队列有8个基因获共享因果变异支持，其中LIPA、GGCX、VAMP8和FES在CARDIoGRAMplusC4D中亦达标；MAP3K8、HLA-DRB5、DSTYK和ADAM9仅在发现队列达标。ERP29在两队列均H3主导，PTPN11在发现队列H3主导、外部队列不确定。Steiger过滤未剔除变异；1995个多变异基因中369个（18.5%）Cochran’s Q显著，40个MR-Egger截距显著；LIPA在发现队列存在方向多效性证据，但关联方向与显著性未变。结论为四个靶点获较强共定位支持，ERP29和PTPN11的MR估计需保留工具变量可能不通过靶基因起作用的警示。
3.3 MR识别治疗靶点的外部验证：以CARDIoGRAMplusC4D为结局重复MR，10个基因Bonferroni显著，其中PTPN11、FES、VAMP8、LIPA、GGCX和ERP29在FinnGen中亦显著。meta分析合并后，LIPA（OR 1.05）、GGCX（1.13）和ERP29（1.28）与CAD风险正相关，PTPN11（0.83）、FES（0.90）和VAMP8（0.90）负相关。结论为六个靶点获两队列MR重复和meta分析支持，并按共定位结果分为Tier A与Tier B。
3.4 组织特异性MR分析：研究人员用GTEx主动脉和皮下脂肪eQTL重复MR，发现LIPA、GGCX和FES在血液、脂肪和动脉组织中的效应方向一致；VAMP8在血液中呈保护关联，但在主动脉和皮下脂肪中方向相反；PTPN11在血液中呈保护关联，在主动脉中方向相反。结论为部分靶点存在组织依赖的效应方向，组织特异结果不构成跨组织验证。
3.5 潜在治疗基因在不同细胞类型中的分布：基于GSE233870的6例冠状动脉粥样硬化样本，质控后保留57,196个细胞，鉴定50个Leiden簇并合并为13个标注细胞群。六基因呈细胞类型异质性表达；均在心肌细胞中相对低表达；VAMP8、ERP29和PTPN11在浆细胞样树突状细胞（plasmacytoid dendritic cells, pDC；仅65个细胞，0.11%，属探索性）中较高；LIPA在巨噬细胞中较高，FES在单核细胞和巨噬细胞中较高，GGCX在成纤维细胞中较高。自动标注的星形胶质细胞群217个细胞因表达完整脂肪细胞程序而被再标注为脂肪细胞。结论为六基因具有不同细胞分布，为后续机制研究提供细胞背景。
3.6 细胞–细胞通讯分析：CellPhoneDB分析显示心脏、血管、基质细胞与免疫细胞间存在密集相互作用，基质和内皮细胞信号活跃，巨噬细胞和pDC具较强发送与接收活性。HLA、淀粉样前体蛋白和Annexin通路作为示例展示；按显著配体–受体相互作用数排序，胶原/整合素黏附、WNT和趋化因子信号更突出。结论为冠状动脉粥样硬化中存在先天免疫–血管/基质细胞通讯网络，但结果为基于配体–受体数据库的模型推断，不证明体内相互作用。
3.7 CAD潜在治疗基因的全表型组MR：在4698个预设基因–表型组合中，4636个获得可用估计；24个关联达Bonferroni阈值（ERP29 17个、FES 5个、PTPN11 2个），132个达FDR<0.05。最强关联包括较高ERP29与甲状腺功能减退、较高FES与较低高胆固醇血症风险、较高PTPN11与较低原发性高血压风险；GGCX、LIPA和VAMP8无Bonferroni显著关联。结论为Phe-MR提供潜在表型后果和安全信号筛查结果，但不能确立药物安全性。
讨论部分总结：原文讨论围绕PTPN11（SHP2）、FES、LIPA、ERP29、GGCX和VAMP8展开，强调PTPN11和VAMP8的组织依赖方向、LIPA在发现队列的敏感性分析警示、ERP29和PTPN11的Tier B共定位局限，以及Phe-MR终生遗传预测与短暂药物暴露之间的差异。局限性包括仅分析cis-eQTL、组织eQTL可用性有限、scRNA-seq仅来自冠状动脉粥样硬化病变、共定位采用近似贝叶斯因子模型并假设单因果变异、部分基因仅单变异工具、潜在样本重叠及部分单细胞群较小等。
研究结论翻译：本研究在两种明显不同的证据层级上优先列出六个CAD候选治疗靶点。四个（LIPA、GGCX、FES和VAMP8；Tier A）在两个队列中均获共定位支持，并获MR和meta分析支持。其余两个（ERP29和PTPN11；Tier B）在两项MR分析中重复，但未达到PP.H4≥0.80：ERP29在两个队列中均H3主导，而PTPN11在发现队列H3主导、外部队列不确定。ERP29占24个Bonferroni显著全表型组关联中的17个，因此在其可被视为治疗靶点前需要显著更多证据。这些分析确立的是遗传关联和优先排序，而非机制或药理调控的因果性。组织特异性结果对两个靶点方向不一致，不构成跨组织验证；细胞–细胞通讯结果为数据库推断。所提示的实验方向包括Tier A靶点的组织相关扰动、Tier B结果所涉及的12q24位点的功能表征，以及利用疾病和对照组织进行独立单细胞验证。需要临床试验检验这些靶点是否具有治疗作用以及能否安全调控。
