《International Journal of Women's Health》：Single-Cell Immune Cis-eQTL Mendelian Randomization Prioritizes Cell-Type-Resolved Candidate Causal Genes for Endometriosis
编辑推荐：
背景：子宫内膜异位症具有显著的免疫成分，但遗传调控疾病易感性的细胞环境仍未完全阐明。
方法：研究人员利用来自20种OneK1K外周血细胞类型的单细胞顺式表达数量性状位点（cis-eQTL）工具变量，针对芬兰基因数据库第12版（FinnGen R12）中的子宫内
背景：子宫内膜异位症具有显著的免疫成分，但遗传调控疾病易感性的细胞环境仍未完全阐明。
方法：研究人员利用来自20种OneK1K外周血细胞类型的单细胞顺式表达数量性状位点（cis-eQTL）工具变量，针对芬兰基因数据库第12版（FinnGen R12）中的子宫内膜异位症数据（20,190例病例；130,160例对照），采用单变异孟德尔随机化（MR）、共定位分析和敏感性分析。本杰明尼-耶库蒂埃利（BY）校正将保守性发现与探索性的本杰明尼-霍赫贝格（BH）筛选区分开来。方向一致性、全血孟德尔随机化和病灶表达提供了背景比较。
结果：在9,157项独特的基因-细胞类型检验中，61个基因的135个关联达到了BH错误发现率（BH-FDR）< 0.05。默认先验下的共定位分析优先识别了跨越28对组合的10个基因（其中6个编码蛋白质）。LOC102723878和RUVBL1同样通过了BY稳健性检验；自然杀伤（NK）细胞中的胰岛素样生长因子1型受体（IGF1R）得到BH支持但未通过BY稳健性检验（后验概率PP4 = 0.89；每标准差表达增加的比值比OR为1.11）。其NK细胞归因反映了工具变量的可用性，并不确立排他性。共定位结果对先验设定敏感。在部分重叠样本中的方向一致性（8/8）仅为描述性。全血孟德尔随机化在6个可检验的蛋白编码基因中提名了3个。IGF1R在病灶中的表达呈现非显著上升趋势（P = 0.101）。
结论：该分析优先识别了用于功能研究的免疫环境候选基因。IGF1R仍是一个探索性的、机制尚未阐明的提名。遗传支持和蛋白质可药性并不能确定治疗获益或生殖安全性；需要实验验证。
研究背景与现存问题
子宫内膜异位症是一种影响约10%育龄女性的常见妇科疾病，以慢性盆腔痛和不孕为主要特征，但其诊断常被延误多年。目前主流治疗方案依赖激素抑制疗法，该方法具有避孕效果、疗效不完全且长期耐受性差，因此亟需基于非激素、作用机制的创新治疗假说。虽然大多数女性均存在经血逆流现象，但仅有少数发展为子宫内膜异位症，这提示免疫清除功能的缺陷在疾病建立中起关键作用。大量证据表明患者外周和腹腔自然杀伤细胞的细胞毒性降低，腹腔巨噬细胞极化为促病变表型，T细胞亚群和B细胞相关反应也发生改变，这些现象将子宫内膜异位症部分定性为一种免疫监视紊乱。然而，现有观察未能识别究竟是哪些基因、在何种免疫细胞类型中因果性地驱动了疾病风险。
人类全基因组关联研究虽已鉴定数十个风险位点并提示激素和炎症通路参与，但绝大多数关联变异位于非编码区，推测它们通过改变基因表达发挥作用。如何将一个非编码关联转化为特定的效应基因及其作用的细胞类型，仍是核心瓶颈。孟德尔随机化结合共定位分析为解决这一问题提供了框架，但应用于整体组织时会平均不同细胞类型的信息，从而掩盖调节关联的细胞环境。大规模单细胞分辨率表达数量性状位点资源的出现，使得在特定免疫细胞环境中进行因果推断成为可能。
研究设计与主要结论
研究人员整合了来自OneK1K队列（982名供体）经eQTL Catalogue统一处理的20种外周血免疫细胞亚群的单细胞顺式表达数量性状位点数据，作为暴露工具变量，并以FinnGen R12数据库中20,190例子宫内膜异位症病例和130,160例女性对照的全基因组关联汇总统计作为结局，开展了双样本、细胞类型分辨的顺式表达数量性状位点孟德尔随机化研究，辅以共定位分析和多项敏感性验证。主要结论如下：在9,157项基因-细胞类型检验中，61个基因的135个关联达到BH-FDR < 0.05的探索性阈值；进一步共定位分析在默认先验下优先识别了10个基因（28对基因-细胞类型组合），其中6个为蛋白编码基因（KLF12、ATG7、RUVBL1、IGF1R、USP53、SZT2），4个为非编码转录本。LOC102723878和RUVBL1同时通过了最严格的BY校正，属于第一层级；其余8个基因为第二层级，包括自然杀伤细胞中的IGF1R（BH-FDR = 0.0085，BY-FDR = 0.082）。IGF1R的NK细胞归因受限于工具变量在该细胞类型中的可获得性，不代表生物学排他性。方向一致性检验在部分重叠样本中显示8/8个可用基因方向一致。整体血液孟德尔随机化基准比较表明，对于这6个蛋白编码基因，细胞类型分辨率改变了用于工具化的变异选择，进而影响了哪些基因被提名。在子宫肌瘤交叉疾病比较中，7/10个基因未达到显著性。在卵巢病灶与配对在位内膜的比较中，4/6个蛋白编码基因差异表达显著，而IGF1R呈非显著上升趋势（P = 0.101）。本研究的重要意义在于，首次在单细胞免疫细胞分辨率下系统性地优先识别了子宫内膜异位症的候选因果基因，并将遗传证据、细胞类型背景和蛋白质可药性注释相结合，为后续功能研究和靶向治疗假说的提出提供了聚焦资源。
关键技术方法
研究人员采用了以下主要技术方法：①从eQTL Catalogue获取OneK1K队列（982名欧洲裔供体）的20种外周血免疫细胞亚型的顺式表达数量性状位点汇总统计，并经SuSiE精细定位获得可信集；②对每个基因在每个细胞类型中选取最大后验包含概率变异作为单一工具变量，限定常染色体、暴露P值 < 5×10-8且F统计量 > 10，并通过GRCh38染色体、位置和等位基因匹配进行工具变量与结局数据的协调；③采用Wald比率法进行单变异孟德尔随机化估计，以一阶delta法计算标准误；④对所有BH显著基因-细胞类型对及其它可工具化的细胞类型，使用coloc.abf在完整±500 kb窗口上进行共定位分析，定义PP4 > 0.8为共定位阳性；⑤采用BH校正作为探索性多重检验控制，BY校正作为任意依赖性下的保守性检验，并进行了变异去重BH和基因内Bonferroni后接BH两种补充分析；⑥通过近似Steiger方向性检验、GWAS Catalog表型组扫描和Open Targets平台可药性注释进行辅助评估；⑦在部分重叠的Rahmioglu等荟萃分析中进行方向一致性检验，在eQTLGen全血顺式表达数量性状位点中进行工具交叉验证，并在GSE7305数据集（10名患者的卵巢病灶与配对在位内膜）中检测差异表达。
研究结果
仪器强度与显着性阈值下的统计功效：经协调后共纳入9,157项独特基因-细胞类型检验，涉及3,821个基因。F统计量中位数为64.9，近似表达方差解释率中位数约为6.2%。在名义α = 0.05水平下，80%功效的最小可检测比值比中位数为1.09；在严格阈值α = 0.05/9,157下，中位数和最大可检测比值比分别为1.18和1.27，35.3%的工具变量在此保守阈值下对OR = 1.15达到80%功效。IGF1R在NK细胞中的F值为137.8，近似R2为0.12。
跨免疫区室的细胞类型分辨孟德尔随机化筛查：总计135个关联（61个基因）达到BH-FDR < 0.05。信号分布于多种免疫谱系而非局限于某一亚群，数量最多的为初始CD4+ T细胞（31个）、中央记忆CD4+ T细胞（27个）、NK细胞（18个）和效应记忆CD8+ T细胞（16个）。
检验间依赖性与依赖性感知显着性层级：9,157次检验使用了7,110个不同的先导变异身份，另有2,047次检验在不同基因和/或细胞类型间重复使用变异。基因内成对检验的符号Z相关性为0.80，随机染色体间比较为-0.002。变异去重BH支持136个基因-细胞关联（61个基因），基因内Bonferroni后接BH支持46个基因。BY校正保留了7个基因的15个关联。IGF1R/NK属于第二层级（BH-FDR = 0.0085；基因水平FDR = 0.012；BY-FDR = 0.082）。
完整窗口共定位优先识别10个基因：对123个BH显著基因-细胞类型对及可工具化细胞类型进行共定位测试（排除MHC区域）。10个基因在PP4 > 0.8水平上共定位，覆盖28对基因-细胞类型组合，所有共定位对均通过近似方向性检验。6个为蛋白编码基因（KLF12、ATG7、RUVBL1、IGF1R、USP53、SZT2），4个为非编码转录本。LOC102723878和RUVBL1为第一层级，其余为第二层级。条件后验概率PP4/(PP3+PP4)在所有共定位对中均≥ 0.83，且在±250 kb和±100 kb窗口上稳定性良好。ESR1产生最强边际MR关联（FDR = 6.6×10-17）但未共定位（PP4 = 0.00）。
NK细胞IGF1R作为探索性候选基因：在工具选择标准下，IGF1R仅在NK细胞中保留了合格且协调的工具变量。其NK细胞顺式表达数量性状位点与子宫内膜异位症在完整窗口中共定位（PP4 = 0.89；PP4/(PP3+PP4) = 0.94）。MR工具rs2311766为SuSiE可信集中最大后验包含概率变异，位于IGF1R转录起始位点下游3 kb处。较高遗传预测的NK细胞IGF1R表达与较高的子宫内膜异位症风险相关（每SD增加OR为1.11）。该结果受四项限定：属第二层级；NK细胞归因描述的是信号可检测的位置而非必须作用的位置；蛋白质可药性不代表治疗效果；局部GWAS信号未达全基因组显著性（最小P≈3.8×10-6）。
非编码共定位转录本的注释：四个非编码转录本包括LOC102723878（长链非编码RNA，距CCDC88B 5'端6 bp）、LINC-PINT（p53诱导的基因间长链非编码RNA）、ENSG00000278022（长链非编码RNA，IGF1R内含子内）和MRPL42P1（转录加工假基因，位于USP53顺式区域内）。两个与共定位蛋白编码候选基因共享先导变异或位于250 kb范围内。
蛋白编码基因的可药性与多效性：IGF1R蛋白携带获批药物类别，ATG7和RUVBL1具有小分子先例，USP53具有基于定位的抗体验证注释，KLF12和SZT2无可药性类别。GWAS Catalog查询发现少量报道关联，如IGF1R与冒险行为，RUVBL1与盆腔器官脱垂等。LOC102723878的先导变异与多种免疫介导疾病相关。
方向一致性与工具交叉验证：在Rahmioglu等荟萃分析中，8/8个可用基因方向与发现一致（3个名义P < 0.05），包括NK细胞IGF1R（OR 1.05，P = 0.087）。在eQTLGen全血中，5/6个蛋白编码先导变异存在顺式表达数量性状位点且方向一致。
全血基准、跨疾病比较与病灶表达：以最强eQTLGen全血顺式表达数量性状位点工具化时，6个基因中3个名义显著，RUVBL1、IGF1R和SZT2不显著。对于IGF1R，与子宫内膜异位症共定位的变异在全血中仅为弱表达数量性状位点，而最强的全血IGF1R表达数量性状位点位于另一无子宫内膜异位症关联的变异。在子宫肌瘤比较中，7/10个基因未达到BH-FDR < 0.05，IGF1R的子宫肌瘤FDR = 0.96。在卵巢病灶vs配对在位内膜中，RUVBL1、USP53和ATG7在病灶中表达降低，KLF12升高；IGF1R呈非显著上升趋势（log2倍变化+0.30，P = 0.101）。
讨论总结与结论翻译
讨论部分强调，该筛查在免疫细胞环境中优先识别了候选基因，并明确区分了BY稳健性共定位基因（LOC102723878、RUVBL1）与探索性BH支持的候选基因（包括IGF1R）。细胞类型归因仅反映工具变量和支持证据的可获得性，不能自行确立排他性或细胞间差异。可药性是蛋白质的属性，不证明其在子宫内膜异位症中的有效性和安全性。
对于IGF1R信号的解读需特别谨慎。MR估计表明较高的NK细胞IGF1R表达与较高的子宫内膜异位症风险相关，按药物靶标孟德尔随机化的逻辑指向抑制IGF1R为降低风险的方向。但这与已知的NK细胞生物学不符——IGF-1/IGF1R信号经典地促进而非抑制NK细胞发育和细胞毒性，因此“高IGF1R损害NK清除异位细胞”的直观叙述缺乏文献支持。该正向关联最好被视为遗传支持但机制未解的观察：NK细胞IGF1R表达数量性状位点可能通过非细胞毒性的NK程序起作用，可能标记而非驱动某种NK状态，或其主要疾病相关效应位于血液中无法测量的效应功能之外。此外，小鼠颗粒细胞特异性Igf1r敲除导致不育，替妥木单抗临床试验报告的高血糖和听力损伤也是额外安全考量。因此，NK细胞IGF1R被呈现为一个优先的候选基因，其治疗方向、细胞机制和安全性均有待解决，明确不作为重定位靶点。
其余蛋白编码基因提出了可供追踪的假说：ATG7（自噬必需酶）暗示细胞存活程序；KLF12（Krüppel样转录抑制因子）在记忆B细胞中保护性共定位且在病灶组织中更高，已有实验联系但本研究的B细胞遗传关联增加了不同细胞环境；RUVBL1（染色质相关AAA+ ATP酶）关联增殖程序；SZT2（KICSTOR复合物组分，营养应激下约束mTORC1信号）在CD4+ T细胞中保护性共定位，提示T细胞代谢检查点假说；USP53（去泛素化酶样蛋白）在B和T细胞亚群中共定位，尚无明确免疫功能，是机制研究的候选。
方法学上，该研究展示了将孟德尔随机化与共定位在细胞类型分辨率下配对的优势和局限。ESR1产生最强边际MR信号却未通过共定位，若仅凭MR显著性报告则不会满足联合证据标准。修订过程中进行的修正对解读至关重要：共定位窗口曾因数据访问接口分页而被截断，所有结果均在完整±500 kb窗口上重新计算，改变了共定位集合。此前版本将一致性样本描述为欧洲裔，现更正为99.3%欧洲裔含少量日本样本。全血基准比较更准确地表述为：细胞类型分辨率改变了哪些变异被用于工具化基因，从而改变了哪些基因被提名。
若干局限性限制了解读：每个暴露仅有一个工具变量，无法应用多工具稳健估计量和MR-Egger/MR-PRESSO分析，无法排除水平多效性；coloc.abf假设每个性状在每个窗口内只有一个因果变异；默认先验支持对较小的p12敏感；方向一致性样本与发现样本存在重叠；外周血欧洲裔表达数量性状位点可能不代表病灶驻留细胞或其他祖先群体；名义功效估计可能高估严格筛选阈值下的检测能力；去标识的整体病灶表达受组织和细胞组成混杂影响，无法验证细胞特异性遗传效应的方向。所有提名均为计算假说，需功能验证和生殖安全性评估。
结论：单细胞免疫顺式表达数量性状位点孟德尔随机化结合完整窗口共定位分析，在一个层级化框架中优先识别了10个子宫内膜异位症候选因果基因，其中6个为蛋白编码基因。该框架将BY稳健性提名（LOC102723878、RUVBL1）与BH和基因水平标准支持的提名区分开来，并将工具变量可用性与筛查证据的变异性加以区别，两者均不能单独确立生物学排他性。NK细胞IGF1R位居第二层级之首，是一个优先的、产生假说的候选基因，其治疗方向、机制和生殖安全性均未解决。这些基因均非治疗就绪靶点；它们的价值在于将实验工作聚焦于特定免疫环境中的特定基因，服务于这个至今仍主要依靠激素抑制治疗的疾病领域。
