多组织eQTL整合用于子宫平滑肌瘤药物靶点发现：12个潜在候选基因的遗传学证据与药物重定位评估

《Intelligent Pharmacy》：Multi-tissue eQTL Integration for Druggable Target Discovery in Uterine Leiomyoma: Genetic Evidence and Drug Repurposing Assessment of 12 Potential Candidates

　　

摘要 背景： 子宫平滑肌瘤（uterine leiomyoma, UL）是女性生殖系统最常见的良性肿瘤，现有药物治疗方案因副作用和生育力保护不足而受限。本研究旨在通过整合多组织表达数量性状位点（expression quantitative trait locus, eQTL）与全基因组关联研究（genome-wide association study, GWAS）数据，识别并优先排序UL的潜在药物靶点。 方法： 研究人员首先从Drug-Gene Interaction数据库（DGIdb v4.0）和Finan等人发表的数据集中整合了可成药基因。随后，利用来自GTEx V.8的全血、卵巢和子宫组织的顺式eQTL数据，结合FinnGen R10的UL GWAS汇总统计数据（34,422例病例，195,888例对照），进行了基于汇总数据的孟德尔随机化（summary-data-based Mendelian randomization, SMR）分析和转录组全关联研究（transcriptome-wide association study, TWAS）。通过贝叶斯共定位分析评估基因表达信号与UL风险是否共享同一因果变异。此外，还进行了全表型组孟德尔随机化（phenome-wide Mendelian randomization, Phe-MR）分析以评估潜在不良效应，并通过检索多个药物数据库评估药物的可开发性。 结果： 经SMR、TWAS和共定位分析三重验证后，研究人员确定了12个高置信度的候选可成药基因：ESR2、BACE1、PIM3、CCHCR1、THBS3、SLC5A6、REV1、SEMA4D、MICU1、ATRAID、NT5C3A和CHRNB1。其中，CCHCR1和THBS3在全血和卵巢组织中均被识别。Phe-MR分析显示，ESR2和CCHCR1表现出显著的多样性关联，而其余10个基因未发现与其他疾病的显著关联。药物评估结果显示，针对ESR2的选择性激动剂（如Erteberel和DPN）、针对PIM3的泛PIM抑制剂（如SGI-1776和TP-3654）以及针对MICU1的临床前小分子抑制剂（如MCU-i4和MCU-i11）具有潜在治疗价值，但需进一步验证其疗效和生殖安全性。CCHCR1、THBS3、REV1、ATRAID、NT5C3A、SLC5A6和CHRNB1则代表了全新的治疗靶点。 结论： 本研究通过多组织eQTL整合分析，成功识别了12个与UL风险相关的可成药基因，并评估了其现有药物开发状态。这些发现为UL的药物开发和药物重定位提供了重要的遗传学依据和优先顺序。

论文解读文章

1. 研究背景与问题提出

子宫平滑肌瘤（uterine leiomyoma, UL）是女性生殖系统最常见的良性肿瘤，在育龄女性中患病率高达20%至40%。UL常导致不孕、盆腔疼痛、月经过多及不良妊娠结局，严重降低患者生活质量。目前，手术干预仍是UL的主要治疗手段，但可能伴随术中出血、血栓栓塞等并发症，且常无法充分保留生育功能。药物治疗方面，促性腺激素释放激素激动剂（Gonadotropin-Releasing Hormone agonist, GnRHa）和米非司酮虽能缓解症状并抑制肌瘤生长，但多数药物会干扰受孕，长期使用可引起不良反应，停药后症状和肌瘤生长易复发，因此不适用于有积极备孕需求的UL患者。近年来，虽然天然化合物如表没食子儿茶素没食子酸酯（epigallocatechin gallate, EGCG）、合成药物如他汀类以及辅助补充剂如维生素D等新兴治疗策略已被探索，但药物干预很少能完全消除肌瘤，且可能存在显著副作用或长期使用限制。因此，迫切需要为育龄期UL女性开发更有效、更安全、更精准的药物治疗方案，这对于生育力保护和改善生殖结局至关重要。

将表达数量性状位点（expression quantitative trait locus, eQTL）分析与全基因组关联研究（genome-wide association study, GWAS）相结合的最新进展，使研究人员能够在特定组织或细胞类型中识别新的、与疾病相关的药物靶点基因。例如，整合eQTL和GWAS数据已成功识别出调控乳腺癌、子宫内膜癌和高级别浆液性上皮性卵巢癌发病机制中多个靶点基因的遗传变异。这些研究为应用药物靶点孟德尔随机化进行复杂疾病的精准医学提供了重要见解，为探索子宫平滑肌瘤的治疗靶点和药物重定位奠定了基础。

2. 研究方法概述

本研究采用了多阶段、多层次的因果推断框架。首先，研究人员整合了来自Drug-Gene Interaction数据库（DGIdb v4.0）和Finan等人数据集的可成药基因。其次，利用来自GTEx V.8项目（欧洲裔人群）的全血（n=670）、卵巢（n=167）和子宫（n=129）三个组织的顺式eQTL数据，以及FinnGen R10的UL GWAS汇总统计数据（34,422例病例，195,888例对照，欧洲裔人群），进行了基于汇总数据的孟德尔随机化（summary-data-based Mendelian randomization, SMR）分析（辅以HEIDI检验排除多效性或连锁不平衡）和转录组全关联研究（transcriptome-wide association study, TWAS）（辅以条件与联合（conditional and joint, COJO）分析识别独立遗传信号）。接着，采用贝叶斯共定位分析（COLOC包）评估基因表达与UL风险是否共享同一因果变异（以PPH4>80%为强共定位证据）。此外，利用来自英国生物银行（UK Biobank）408,961名欧洲裔参与者的SAIGE GWAS数据库中的783种表型，进行了全表型组孟德尔随机化（phenome-wide Mendelian randomization, Phe-MR）分析以评估潜在不良效应。最后，通过检索DrugBank、ChEMBL、ClinicalTrials和UniProt等多个药物数据库，评估了候选基因的药物开发状态和药物重定位机会。

3. 研究结果

3.1 可成药基因组

研究人员从DGIdb v4.0数据库获得3,952个潜在可成药基因，从Finan等人的综述中获得4,463个基因，整合去重后最终得到5,883个独特的可成药基因用于后续分析。

3.2 UL的候选可成药基因

在SMR分析中，研究人员利用全血eQTL数据集检测到17个可成药基因，利用卵巢eQTL数据集检测到10个可成药基因（FDR校正P<0.05；HEIDI P>0.01），但基于子宫eQTL数据集未发现任何可成药基因。在TWAS分析中，研究人员在全血中检测到44个显著基因，卵巢中24个，子宫中18个（FDR校正P<0.05；P-joint<0.05）。通过交集SMR和TWAS结果，研究人员在全血中识别出11个可成药基因，卵巢中9个，共计18个独特基因，其中CCHCR1和THBS3在两个组织中均被识别。随后的共定位分析进一步支持了这些发现，在全血中识别出9个候选基因，卵巢中5个，其中CCHCR1和THBS3在两个组织中共享，显示出基因表达与UL风险之间的强共定位证据（PPH4>0.8）。具体而言，CCHCR1和THBS3在全血和卵巢中均观察到强共定位；ESR2、BACE1、PIM3、SLC5A6、REV1、SEMA4D和MICU1在全血中；ATRAID、NT5C3A和CHRNB1在卵巢中。基于SMR估计，这些基因与UL风险的关联方向各异：全血中ESR2、BACE1、CCHCR1、THBS3、REV1和SEMA4D的高表达，以及卵巢中CCHCR1、THBS3、ATRAID和NT5C3A的高表达与UL风险降低相关；相反，全血中PIM3、SLC5A6和MICU1的低表达，以及卵巢中CHRNB1的低表达也与较低的UL风险相关。

3.3 潜在多效性效应和安全性信号的评估

通过整合SMR和Phe-MR分析，曼哈顿图显示，在12个候选基因中，仅ESR2（全血）和CCHCR1（全血和卵巢）表现出显著的多效性关联。具体而言，ESR2高表达与心律失常风险降低相关。CCHCR1表现出显著的组织特异性效应：其在卵巢中的高表达降低了2型糖尿病和烟草使用障碍的风险，但显著增加了甲状腺功能减退症的风险；而在全血中的高表达则与骨关节炎和其他眼部疾病风险降低相关。其余10个候选基因（BACE1、SLC5A6、THBS3、MICU1、ATRAID、NT5C3A、CHRNB1、REV1、SEMA4D和PIM3）未显示出与其他疾病的显著关联。

3.4 候选可成药基因的现有药物和潜在治疗策略评估

研究人员对12个候选基因的药物开发进展进行了系统评估。对于需要上调以降低UL风险的基因（ESR2、BACE1、CCHCR1、THBS3、REV1、SEMA4D、ATRAID、NT5C3A）：ESR2的现有激活剂选择性低，但选择性ESR2激动剂（如Erteberel和DPN）和天然化合物（如liquiritigenin和MF101）具有潜力；BACE1和SEMA4D的现有药物均为抑制剂，与所需的上调方向相反；CCHCR1、THBS3、REV1、ATRAID和NT5C3A尚无上市或在研药物。对于需要下调以降低UL风险的基因（PIM3、MICU1、SLC5A6、CHRNB1）：PIM3的获批药物fostamatinib存在生殖毒性，但泛PIM抑制剂和天然化合物（如coptisine和timosaponin A3）值得关注；MICU1的临床前抑制剂（MCU-i4和MCU-i11）需进一步验证体内疗效和生殖安全性；CHRNB1的现有非选择性神经肌肉阻断剂不适合精准治疗；SLC5A6尚无直接抑制剂。

4. 讨论与结论总结

讨论部分总结： 本研究通过全面的全基因组跨性状分析，为识别UL的潜在可成药基因靶点提供了宝贵见解。研究发现强调了组织特异性基因表达在UL治疗中的关键作用，例如CCHCR1和THBS3在两种组织中均发挥保护作用，而其他10个基因仅在单一组织中表现出保护性关联。Phe-MR结果进一步凸显了组织特异性的重要性：CCHCR1在卵巢中的高表达既带来治疗益处也伴随潜在副作用，而在全血中的高表达则未观察到显著不良结局，这为设计组织选择性药物提供了关键框架。研究还为ESR2与UL之间的既往争议提供了新的遗传学证据，表明更高ESR2表达具有潜在保护作用，这可能解释了以往观察性研究中发现的局部肿瘤ERβ水平升高是一种代偿性生物学反应而非病因。通过对现有药物数据库的系统检索，研究发现大多数现有药物因选择性低、效应方向不一致或生殖安全性问题而不适用于UL治疗，但若干候选药物和天然化合物值得优先验证。研究的优势包括聚焦于预验证的可成药基因以增强转化潜力、整合三种组织的eQTL数据以观察组织特异性效应、采用多阶段因果推断框架减少假阳性，以及通过Phe-MR系统评估脱靶效应。局限性包括eQTL和全表型组数据集包含男性参与者可能引入性别偏倚、GWAS数据主要来自欧洲裔人群限制了普适性、全基因组跨性状分析不能替代实验验证或随机对照试验（randomized controlled trial, RCT）、子宫组织SMR分析未产生显著靶点可能由于样本量较小，以及候选药物必须经过严格的功能测试才能进入临床转化。

结论部分翻译： 总之，通过整合多组织eQTL和UL GWAS数据，并采用包含SMR、TWAS和贝叶斯共定位的三重验证策略，本研究识别了12个具有遗传学支持的与UL风险相关的可成药基因。药物评估结果表明，尽管ESR2、PIM3和MICU1已有候选药物，但这些药物受到低选择性、生殖毒性或体内验证不足等限制，因此开发高选择性的新型药物应成为未来研究的优先事项。对于BACE1和SEMA4D，当前的药物研发工作仅集中于抑制剂，而本研究的发现提示靶点激活在UL中的潜在治疗作用值得进一步探究。此外，CCHCR1、THBS3、REV1、ATRAID、NT5C3A、SLC5A6和CHRNB1目前缺乏潜在的候选药物，作为新型治疗靶点具有巨大潜力。未来的药物开发应整合实验研究和随机对照试验来验证这些遗传学发现，阐明潜在的生物学机制，并严格评估靶向干预的安全性和临床疗效。