单细胞表达数量性状基因座孟德尔随机化揭示多囊卵巢综合征中免疫细胞特异性因果调控网络与可干预靶点

《Science Progress》：Single-cell expression quantitative trait locus Mendelian randomization reveals immune cell-specific causal regulatory networks and actionable targets in polycystic ovary syndrome

　　

目的：系统研究多囊卵巢综合征（polycystic ovary syndrome, PCOS）的发病是否与特定免疫细胞亚群中的基因表达失调存在因果关系，并评估这些因果基因作为可干预药物靶点的潜力。
方法：本研究采用两样本孟德尔随机化（Mendelian randomization, MR）框架，利用公开的全基因组关联研究（genome-wide association study, GWAS）汇总统计数据。参与者数据包括797例PCOS病例和140,558例对照（未直接招募患者）。工具变量来源于14种免疫细胞类型的高分辨率免疫细胞特异性单细胞表达数量性状基因座（single-cell expression quantitative trait locus, sc-eQTL）数据（OneK1K项目）。主要分析采用逆方差加权（inverse-variance weighted, IVW）法。共享因果变异通过贝叶斯共定位验证。进行了全表型关联分析（phenome-wide association analysis, PheWAS）、外部转录组数据集验证（GSE8157）和DrugBank数据库筛选，以评估多效性和药物重定位。
结果：MR分析显示，GLIPR1在非经典单核细胞（non-classical monocytes, Mono NC）和XBP1在CD4⁺效应记忆T细胞（CD4⁺ effector memory T cells, CD4 ET）中的表达与PCOS风险存在全基因组显著因果关联。GLIPR1表达升高与PCOS风险降低相关（OR = 0.669，P = 4.34×10^-6），而XBP1表达升高与风险增加相关（OR = 1.406，P = 9.53×10^-8）。共定位分析证实GLIPR1与PCOS共享一个因果变异（PP.H4 = 96.73%）。PheWAS和外部验证确认了GLIPR1的安全性特征和显著上调（P = 0.03）。药物重定位确定SOT-107（一种III期蛋白治疗药物）为GLIPR1的潜在相互作用药物。
结论：这项sc-eQTL MR研究揭示了PCOS中免疫细胞特异性因果调控网络。非经典单核细胞中的GLIPR1代表一个高置信度保护性靶点，而XBP1为免疫介导的发病机制提供了提示性证据。候选药物SOT-107突出了理论上的重定位机会，但仍需严格的临床前验证。

多囊卵巢综合征（polycystic ovary syndrome, PCOS）是一种异质性内分泌疾病，影响育龄女性，其管理依赖抗苗勒管激素（anti-Müllerian hormone, AMH）、雄激素及胰岛素抵抗稳态模型（homeostatic model assessment for insulin resistance, HOMA-IR）等多种标志物。雄激素过量破坏卵泡成熟，代谢紊乱与长期健康风险相关。近年研究发现免疫系统在PCOS发病中起核心作用，患者存在全身及卵巢局部免疫稳态失衡，T细胞、单核/巨噬细胞和树突状细胞等免疫细胞亚群的异常极化与浸润参与病理进展。然而，既往全基因组关联研究（genome-wide association study, GWAS）虽然鉴定出LHCGR、FSHR、INSR等风险位点，但难以确立关联背后的因果性。孟德尔随机化（Mendelian randomization, MR）利用遗传变异作为工具变量可推断基因表达与表型的因果效应，但基于组织bulk eQTL的研究无法捕捉细胞类型特异的调控机制。单细胞RNA测序（single-cell RNA sequencing, scRNA-seq）和单细胞表达数量性状基因座（single-cell expression quantitative trait locus, sc-eQTL）技术为解析免疫细胞亚群特异的调控网络提供了可能，但尚缺乏在PCOS中的系统性应用。为此，研究人员开展了本研究。

研究人员利用OneK1K项目14种免疫细胞类型的sc-eQTL数据作为暴露，以IEU Open GWAS数据库中PCOS GWAS（797例病例和140,558例对照）为结局，进行两样本MR、贝叶斯共定位、全表型关联分析、外部数据集验证及药物重定位分析，旨在揭示免疫细胞特异性因果调控网络并评估可干预靶点。该研究发表于《Science Progress》。结果发现非经典单核细胞中的GLIPR1（GLI pathogenesis related 1）为保护性靶点，CD4⁺效应记忆T细胞中的XBP1（X-box binding protein 1）为风险性靶点，并筛选出候选药物SOT-107。这些发现为理解PCOS的免疫-代谢失衡机制提供了新视角，并为精准靶向治疗奠定了基础。

主要技术方法包括：从OneK1K sc-eQTL数据（14种免疫细胞类型）筛选条件独立顺式eQTL（cis-eQTL）作为工具变量，阈值P < 5.0×10^-8、连锁不平衡（linkage disequilibrium, LD）聚簇r² < 0.001、F统计量>10；使用TwoSampleMR包（v0.6.6）进行两样本MR，单个SNP用Wald ratio，多个用逆方差加权法；贝叶斯共定位（coloc R包）评估共享因果变异；利用AstraZeneca PheWAS平台（约50万UK Biobank外显子组测序数据）评估多效性和脱靶效应；使用独立外周血全转录组数据集GSE8157（22例PCOS患者和22例对照）验证差异表达；通过DrugBank（v6.0）筛选药物-基因相互作用。所有数据均来自公开数据库，未直接招募患者。

基因表达对PCOS的免疫细胞类型特异性效应：通过两样本MR分析，在FDR校正阈值P < 6.78×10^-6下，鉴定出两个免疫细胞特异性eGenes（eQTL调控的基因）。在非经典单核细胞（Mono NC）中，GLIPR1表达每增加一个标准差，PCOS风险降低33.1%（OR = 0.669，P = 4.34×10^-6）；在CD4⁺效应记忆T细胞（CD4 ET）中，XBP1表达每增加一个标准差，PCOS风险增加40.6%（OR = 1.406，P = 9.53×10^-8）。其余12种免疫细胞类型中未达到全基因组显著性，所有工具变量F统计量均大于10（18.6–42.3），排除了弱工具偏倚。

共定位分析优先排序关键eGenes：贝叶斯共定位分析显示，GLIPR1在Mono NC中的表达与PCOS共享因果变异的后验概率PP.H4 = 96.73%，超过80%的强证据阈值，且先导变异rs7971504同时驱动eQTL和GWAS信号；XBP1在CD4 ET中的PP.H4 = 65.4%，超过50%的中度证据阈值，表明GLIPR1与PCOS共享因果变异的证据更强。

全表型关联分析评估潜在脱靶效应和多效性：利用AstraZeneca PheWAS平台分析，GLIPR1相关遗传变异在预设阈值（P < 1×10^-6）下与恶性肿瘤、自身免疫病及主要心血管不良事件无显著关联，但与多项代谢相关表型显著关联，最显著的是肌酐水平（P = 3.97×10^-7，效应值 = -4.34），并与极低密度脂蛋白（very-low-density lipoprotein, VLDL）相关脂质成分、支链氨基酸（branched-chain amino acid, BCAA）和糖酵解代谢物等NMR代谢组特征呈负相关，表明GLIPR1表达升高与更有利的代谢谱相关。

关键eGene表达差异的外部临床数据集验证：在独立外周血全转录组数据集GSE8157中，GLIPR1 mRNA表达在PCOS患者组显著高于正常对照组（log2表达值6.82 ± 0.45 vs 6.48 ± 0.52，Student's t-test P = 0.03），验证了GLIPR1在PCOS中上调，可能反映内源性代偿防御机制。

PCOS相关eGenes的药物重定位：通过DrugBank数据库（v6.0）筛选，鉴定出SOT-107（一种正处于脑癌III期临床试验的靶向蛋白治疗药物）为与GLIPR1相互作用的高置信候选化合物。SOT-107是转铁蛋白和白喉毒素的偶联物，目前无PCOS或其他生殖内分泌疾病的临床试验报告。数据库中未发现其他与GLIPR1或XBP1直接相互作用的已批准或临床阶段药物。

讨论部分指出，非经典单核细胞和CD4⁺记忆表型T细胞是PCOS相关免疫特征的主要贡献者。GLIPR1与肌酐及多种NMR代谢组特征关联，提示其可能作为整合全身代谢稳态、肾脏排泄机制和免疫调节的关键枢纽。在PCOS患者中观察到的GLIPR1上调可能是一种代偿性防御机制，旨在抑制异常炎症反应或维持局部代谢稳态。XBP1是未折叠蛋白反应（unfolded protein response, UPR）的主转录因子，其过表达可能通过促进Th17极化、抑制Treg分化以及驱动胆固醇合成等机制加速PCOS病理。先天免疫（非经典单核细胞）和适应性免疫（CD4⁺记忆T细胞）在PCOS发病中发挥独立而协同的作用。与近期多组学MR研究（如CBLN3、NRBP1、CDCP1等）相比，本研究的独特贡献在于单细胞分辨率、新机制维度以及通过PheWAS提供的安全性评估。药物重定位分析发现SOT-107可能通过增强GLIPR1生物学功能改善免疫代谢表型，但当前证据仅限于计算预测，需严格的临床前验证。

研究结论部分翻译如下：本研究利用sc-eQTL孟德尔随机化和共定位分析揭示了PCOS中免疫介导的遗传调控机制。非经典单核细胞中的GLIPR1和CD4⁺T细胞中的XBP1被确定为驱动免疫代谢失衡的关键因果基因。SOT-107等精准治疗候选药物为打破疾病循环提供了有前景的途径，为未来生殖内分泌领域的转化研究奠定了坚实基础。