Integrative Functional Genomics Identifies Candidate Brain and Pituitary Regulatory Mechanisms of Reproductive Aging

技术摘要：整合功能基因组学识别候选大脑与垂体调节机制，研究生殖衰老

问题陈述
自然绝经年龄（ANM）是一个复杂的生殖衰老性状，具有显著的遗传力（估计为 44–85%），并对女性健康（包括心血管代谢、神经系统和心理健康结果）具有重大影响。虽然以往的全基因组关联研究（GWAS）已经识别出许多与 ANM 相关的位点，但这些关联的功能解释仍然具有挑战性。调节变异通常以组织或细胞类型特异性的方式影响远端基因，而广泛的神经内分泌系统内的特定调节机制——特别是大脑和垂体中的机制——仍有待进一步表征。现有研究主要强调卵巢生物学，这导致在理解 ANM 相关变异如何在特定的神经和神经内分泌细胞群体中发挥作用方面存在空白。

方法论
作者开发了一个整合功能基因组学框架，利用了来自 173,4424 名欧洲血统女性的大规模 GWAS 汇总统计数据。分析工作流通过以下几个阶段进行：

1. 遗传架构与相关性： 使用连锁不平衡得分回归（LDSC）估算 SNP 遗传力和 ANM、重度抑郁障碍（MDD）以及焦虑症（ANX）之间的遗传相关性。
2. 位点定义与精细映射： 使用 FUMA 定义独立的基因组风险位点。在每个位点内进行精细映射，采用近似贝叶斯因子（ABF）方法（通过 coloc 实现），以优先排序具有后验包含概率（PIP）> 0.03 的代表性因果变异（位于 90% 置信集内）。
3. 基因优先排序与多层注释： 使用 FUMA 进行基因映射，并利用多层 xQTL 证据（eQTL、sQTL、pQTL 和单细胞 eQTL）进行注释。
4. 转录组范围关联分析（TWAS）： 使用 FUSION 在卵巢、垂体和下丘脑组织上进行联合/条件 TWAS，以优化基因表达关联，并考虑连锁不平衡及预测表达信号之间的相关性。
5. 孟德尔随机化（MR）： 进行基于顺式 eQTL 的 MR，以识别基因预测的基因表达与 ANM 之间的稳健关联。分析使用了来自 GTEx v8 的 15 种组织（包括 13 个脑区、卵巢和垂体）的条件独立顺式 eQTL。主要工具要求 $P_{eQTL} < 0.05$，而敏感性分析要求达到全基因组显著性水平（$P_{eQTL} < 5 \times 10^{-8}$）。
6. 共定位： 对 MR 优先排序的基因应用贝叶斯共定位分析（coloc.abf），以确定 ANM 和基因表达信号是否共享一个因果变异（定义为 $PP.H4 \geq 0.80$）。
7. 单细胞调节图谱： 将精细映射的变异与人类大脑和垂体的 SCENT（单细胞增强子靶基因映射）数据相结合，以识别细胞类型特异性的增强子-基因关系。
8. 跨性状与药物分析： 评估 ANM 优先排序基因与独立优先排序的 MDD 和 ANX 基因之间的重叠情况。使用药物-基因相互作用数据库（DGIdb）对药物-基因相互作用进行注释。

关键结果

• 遗传架构： 本研究识别出 230 个 ANM 基因组风险位点，其中 108 个为新发现。精细映射优先排序了 216 个代表性变异（其中 168 个为新发现），均位于 90% 置信集内。
• 遗传相关性： ANM 与 MDD ($r_g = 0.050$) 和 ANX ($r_g = 0.044$) 显示出轻微但具有统计学意义的正向遗传相关性。
• 基因优先排序： TWAS 在卵巢（如 GCKR, UIMC1）和垂体（如 TMEM106A）中发现了强关联。MR 识别出 352 个稳健的基因-组织关联，包括跨五个脑组织的 HSD17B6 以及小脑中的 NBR1。
• 共定位： 对于 43 个 MR 优先排序的基因，发现了存在共享因果变异的强有力证据（$PP.H4 \geq 0.80$）。其中七个基因（PRIM1, L3MBTL3, PNP, EBF3, CCDC88C, GBAS, TYMP）表现出特别强的区域共定位性（$PP.H4 > 0.99$）。
• 细胞类型特异性机制： 与 SCENT 的整合突显了四个具有特定调节背景的候选基因：
  • HSD17B6：与大脑中的抑制性神经元相关。
  • NBR1 和 HGS：与少突胶质前体细胞（OPCs）相关。
  • ERBB4：与垂体泌乳细胞相关。
• 跨性状重叠： 16 个 ANM 优先排序基因与 MDD 重叠，10 个与 ANX 重叠，且两个重叠集合之间没有共同基因。显著的重叠包括 ERBB4 (MDD) 和 GABBR1 (ANX)。药理学注释显示这些基因与包括抗精神病药物和抗肿瘤药物在内的已批准药物存在相互作用。

意义与主张
本文声称通过从仅识别位点转向优先排序候选因果变异、基因及特定调节机制，推进了对 ANM 遗传学的解释。研究强调，生殖衰老是在更广泛的神经内分泌系统中协调进行的，而非仅仅局限于卵巢。

作者强调，其整合方法为特定基因提供了互补证据：

• HSD17B6 受到多个脑组织 MR 结果以及细胞类型特异性（抑制性神经元）证据的支持。
• NBR1 和 HGS 受到 OPCs 证据的支持，暗示了少突胶质谱系细胞在生殖衰老生物学中的潜在作用。
• ERBB4 受到垂体 TWAS 和泌乳细胞特异性调节的支持，将其与神经内分泌功能及精神疾病的遗传架构联系起来。

研究谨慎地总结道，虽然这些发现为功能性研究提供了特定的位点和基因优先级，但并未建立治疗效果或因果机制。作者明确指出，需要实验研究来确定这些调节关联如何影响细胞功能，以及它们的效果是否随生殖或激素状态而变化。研究结果受限于欧洲血统女性，且作者承认，由于缺乏多样化的祖先数据以及使用基于汇总统计量的精细映射方法，精细映射的分辨率和普适性可能受到限制。