《Frontiers in Genetics》：Multi-omics mendelian randomization integrating GWAS and eQTL data revealed potential drug target for irritable bowel syndrome
编辑推荐：
摘要
引言：肠易激综合征（Irritable Bowel Syndrome, IBS）是一种主要影响年轻人群及女性的常见胃肠道疾病，现有治疗手段有限，因此亟需识别新的药物靶点。本研究旨在通过利用大规模人类遗传数据，鉴定并优先排序IBS的新型、经遗传学验证的药物
摘要
引言：肠易激综合征（Irritable Bowel Syndrome, IBS）是一种主要影响年轻人群及女性的常见胃肠道疾病，现有治疗手段有限，因此亟需识别新的药物靶点。本研究旨在通过利用大规模人类遗传数据，鉴定并优先排序IBS的新型、经遗传学验证的药物靶点。
方法：研究人员开展了一项系统的可成药基因组范围孟德尔随机化（Mendelian Randomization, MR）分析，以评估5,642个潜在可成药蛋白编码基因对IBS风险的因果效应。该分析整合了当前最大规模可用的IBS全基因组关联研究（Genome-Wide Association Study, GWAS）汇总统计数据（包括53,400例病例与433,201例对照）以及全面的血液表达数量性状位点（expression Quantitative Trait Loci, eQTL）数据。显著发现进一步采用共定位分析（colocalization analysis）进行验证。通过表型组范围关联研究（Phenome-Wide Association Study, PheWAS）评估潜在的靶向不良反应风险。最后，利用药物特征数据库（Drug Signatures Database, DSigDB）和分子对接预测潜在的治疗化合物。
结果：MR分析确定了八个与IBS存在潜在因果关联的基因。经过严格的共定位分析验证，EP300和P2RY14成为最有希望的候选靶点。遗传学预测的较高EP300表达（比值比[Odds Ratio, OR]: 1.128，95%置信区间[Confidence Interval, CI]: 1.079-1.180）和P2RY14表达（OR: 1.118，95% CI: 1.067-1.172）均与IBS风险增加存在提示性因果关联。PheWAS分析表明，EP300和P2RY14未与其他任何表型呈现全基因组显著关联。此外，分子对接预测卡托普利（captopril）和甲萘醌（menadione）等现有化合物可有效结合EP300蛋白。
结论：本研究提供了确立EP300和P2RY14作为IBS治疗潜力药物靶点的遗传学证据，为未来的药物开发与老药新用工作奠定了基础。
研究背景
肠易激综合征（IBS）是一种常见的慢性功能性胃肠病，主要表现为无器质性病变证据的反复发作性腹痛或不适，并伴有排便习惯改变。该病全球患病率为1.1%至35.5%，女性发病率约为男性的1.5至2倍，且在儿童中亦属高发，给医疗系统带来了沉重负担。尽管临床常采用饮食调整、生活方式干预及解痉药、止泻药等对症治疗，但疗效有限且患者复发率高，加之IBS亚型异质性大，导致“一刀切”疗法难以奏效。因此，寻找新型治疗靶点迫在眉睫。近年来，将人类遗传学证据融入药物研发被证明是提高临床试验成功率的有效策略，特别是利用孟德尔随机化（MR）整合GWAS与eQTL数据，能够推断基因表达终生变化对疾病风险的因果效应。本研究正是基于此背景，旨在通过系统的遗传学分析挖掘IBS的潜在药物靶点。
主要技术方法
研究人员首先整合了Drug-Gene Interaction Database（DGIdb）及Finan等人的综述构建了包含5,642个可成药基因的列表。暴露变量选用eQTLGen联盟的血液cis-eQTL数据，结局变量选用GWAS Catalog中最大的IBS GWAS汇总数据（含53,400例病例与433,201例对照，均为欧洲血统）。采用两样本MR分析，利用逆方差加权法（IVW）及Wald比率法估计因果效应，并通过Cochran’s Q检验、MR-Egger回归及Steiger方向性检验进行敏感性分析。随后进行贝叶斯共定位分析（使用coloc包）验证共享因果变异。安全性评估通过UK Biobank的PheWAS分析完成。靶点表达验证结合了GTEx数据库、临床样本RT-qPCR及公共数据集GSE166869和GSE168759的批量RNA测序。下游机制探索采用GRNBoost2算法构建基因调控网络，并通过单细胞转录组数据（Gut Cell Atlas）解析细胞特异性。药物预测利用DSigDB及AutoDock Vina进行分子对接。
研究结果
3.1 可成药基因的工具变量
通过整合数据库并排除假基因，共获得5,642个可成药基因，匹配cis-eQTL后最终筛选出26,101个合格的单核苷酸多态性（SNP）作为工具变量。
3.2 候选可成药基因
MR分析初步筛选出8个与IBS存在潜在因果关联的cis-eQTL基因（FDR<0.05）。其中EP300（OR: 1.128）和P2RY14（OR: 1.118）的表达增加与IBS风险升高相关，而HLA-C、PAM、PRMT5则呈保护性关联。敏感性分析未发现显著异质性与水平多效性，且支持从基因表达到IBS风险的因果方向。
3.3 共定位分析
鉴于工具变量可能存在连锁不平衡，研究人员进行了共定位分析。结果显示，在8个候选基因中，EP300和P2RY14的eQTL信号与IBS GWAS信号存在强共定位证据（PP.H4 > 0.8），表明两者共享相同的因果SNP，排除了连锁不平衡导致的假阳性。
3.4 表型组范围关联分析
PheWAS分析显示，EP300和P2RY14在严格的全基因组显著性阈值（P < 10-8）下未与其他任何表型产生关联。这表明靶向这两个基因可能引发的脱靶效应较少，具备较好的安全性前景。
3.5 基因调控与组织表达谱分析
利用GTEx数据分析发现，调控EP300和P2RY14表达的SNP在不同组织中均表现出显著的基因型-表达量关联。此外，这两个基因在人类小肠组织中均呈中度表达，提示其可能在肠道局部发挥生物学功能。
3.6 靶向药物预测与分子对接
DSigDB预测有四个小分子化合物可能下调EP300表达。分子对接结果显示，柔红霉素（daunorubicin）与EP300结合能最强（-12.47 kcal/mol），卡托普利、甲萘醌和伏立诺他（vorinostat）也显示出良好的结合活性。
3.7 EP300及其潜在调控靶点在IBS中上调
在两个独立RNA-seq数据集（GSE166869和GSE168759）中，IBS患者肠道组织的EP300表达均显著升高。临床样本RT-qPCR验证也证实，IBS患者血液中EP300表达较健康对照上调约1.8倍。通过整合数据集并构建EP300调控网络（Regulon），研究人员鉴定出382个直接靶基因，其中有92个在IBS中显著上调。功能富集分析显示这些靶基因参与了皮质微管组织、细胞命运决定以及Notch、Rap1和mTOR信号通路。
3.8 IBS相关EP300调控靶点在肠上皮富集
利用人肠道细胞图谱（Gut Cell Atlas）的单细胞数据，研究发现IBS相关的EP300调控网络活性高度富集于肠上皮谱系（包括吸收性肠上皮细胞、分泌细胞和祖细胞），而在免疫细胞（B/T细胞）中活性较低。EWCE分析进一步确认该富集主要位于BEST4+上皮细胞、肠上皮细胞及簇状细胞中，提示EP300可能通过破坏肠上皮稳态促进IBS发病。
讨论与结论
讨论部分指出，P2RY14作为一种嘌呤能受体，在细胞应激时被激活，可能通过放大受损黏膜的信号传导参与IBS的低度炎症。EP300作为关键的转录共激活因子，其过表达可能通过改变组蛋白修饰（如H3K27ac）触发炎症相关基因表达，破坏肠上皮屏障功能。研究还探讨了老药新用的可能性，卡托普利和甲萘醌因具备抗炎特性且安全性较好，成为潜在的候选药物。尽管本研究样本主要基于欧洲血统且eQTL来自血液，但仍为IBS的药物研发提供了坚实的遗传学基础。
结论翻译
本研究通过整合多组学孟德尔随机化分析，提供了确立EP300和P2RY14作为肠易激综合征治疗潜力药物靶点的遗传学证据。这一发现不仅揭示了IBS潜在的新生物学通路，也为未来的药物开发与老药新用工作奠定了数据驱动的基础，有望加速针对IBS更有效疗法的研发进程。该研究发表于《Frontiers in Genetics》。
