IF50.5｜卵巢衰老罕见变异与癌症风险及新生突变率遗传关联

罕见编码变异揭示卵巢衰老新基因，并提示早绝经遗传风险可能增加子代新生突变。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

女性生殖寿命受卵巢储备初始大小和卵泡丢失速率影响，绝经年龄个体差异大，且与生育力及多种疾病风险相关。既往大规模全基因组关联研究（GWAS）已鉴定约300个常见变异位点，解释10-12%的绝经年龄变异，其中三分之二涉及DNA损伤应答（DDR）通路，提示卵母细胞对DNA损伤的敏感性是卵巢衰老的核心机制。然而，常见变异效应较小，且主要位于非编码区，对具体致病基因和生物学机制的解析有限。罕见蛋白编码变异往往具有更大效应，能直接指向基因功能，但此前基于外显子组测序的研究仅发现少数基因，且样本量和分析方法存在局限。

科学问题

卡在哪里： 此前对绝经年龄的遗传研究主要关注常见变异，罕见编码变异分析不足，发现的基因数量有限，且效应大小和生物学机制尚不清楚。此外，卵巢衰老与癌症易感性及生殖细胞突变率之间的遗传联系缺乏系统研究。

本文要回答： 利用大规模外显子组测序数据系统鉴定影响绝经年龄的罕见蛋白编码变异，并探索这些变异与癌症风险及子代新生突变率的关联。

技术路线

作者主要用了：全外显子组测序、基因负担检验、BOLT-LMM、PheWAS、多基因评分、新生突变分析。

作者利用英国生物银行106,973名欧洲裔绝经后女性的全外显子组测序数据，通过基因负担检验（BOLT-LMM）系统评估罕见蛋白截短变异（PTV）和有害错义变异对绝经年龄的影响。随后在deCODE和BRIDGES队列中验证，并结合GWAS数据、癌症表型组关联分析以及100kGP和deCODE的三人家系全基因组测序数据，探讨这些基因的疾病关联及对子代新生突变率的影响。

图 1　Fig. 1。图1展示了基因负担检验的Manhattan图和QQ图。Manhattan图显示9个基因达到外显子组显著性（P<1.08×10-6），其中ZNF518A、ETAA1等为新发现。QQ图显示λ值接近1（1.047-1.050），表明统计模型校准良好，无系统性偏差。

核心亮点

亮点 1｜发现五个新的卵巢衰老基因

外显子组负担分析鉴定出9个与绝经年龄显著相关的基因，其中5个为首次报道：ETAA1、ZNF518A、PNPLA8、PALB2和SAMHD1。效应量远超常见变异，如ZNF518A蛋白截短变异携带者绝经年龄提前5.61年，SAMHD1有害变异携带者绝经年龄推迟1.35年。这些基因中7个参与DNA损伤修复，支持DDR在卵巢衰老中的核心作用。

图 2　Fig. 2。图2展示了9个显著基因的效应量及置信区间。ZNF518A蛋白截短变异效应最大（-5.61年），SAMHD1有害变异效应为+1.35年。不同变异类别（HC-PTV、CADD≥25错义、DMG）的效应方向一致，支持这些基因的真实关联。

亮点 2｜ZNF518A缩短生殖寿命

ZNF518A罕见蛋白截短变异与绝经年龄提前5.61年、初潮年龄推迟0.56年相关，表明其功能缺失缩短女性生殖寿命。ZNF518A编码一个C2H2锌指转录因子，与PRC2和G9A-GLP抑制复合物相互作用。ChIP-seq数据富集分析显示，绝经年龄相关常见变异显著富集于ZNF518A结合位点，提示其通过转录抑制调控生殖衰老相关基因。

图 3　Fig. 3。图3展示了SAMHD1等基因与癌症风险的性别分层关联。SAMHD1有害变异在男性和女性中均显著增加全癌风险，且在激素敏感性癌症中效应更明显。图中显示Bonferroni校正后的显著性阈值，SAMHD1在多个癌症类型中达到显著。

亮点 3｜SAMHD1变异增加癌症风险

SAMHD1有害变异与男性和女性全癌风险增加相关（男性OR=2.12，女性OR=1.61），并在deCODE队列中成功复制。癌症类型分析显示其增加前列腺癌、间皮瘤等多种癌症风险。SAMHD1变异同时与绝经年龄推迟相关，类似于CHEK2，提示其可能通过干扰DNA损伤感知和凋亡，减缓卵巢储备耗竭，同时增加癌症易感性。

图 4　Fig. 4。图4展示了9个绝经年龄相关基因与90种癌症的PheWAS结果。SAMHD1与多种癌症风险增加相关，包括前列腺癌、间皮瘤等。图中标注了达到外显子组显著性的关联，y轴截断在-log10(P)=30，完整统计量在补充表中。

亮点 4｜早绝经遗传风险与子代新生突变率

在100kGP的8,089个三人家系中，母亲绝经年龄多基因评分（PGS）与子代母源新生突变率显著负相关（β=-0.082 DNMs/s.d., P=0.00033），即遗传上倾向于早绝经的母亲其子代携带更多母源新生突变。但在deCODE队列中未能复制（P=0.53），meta分析仍显著（P=0.016）。孟德尔随机化分析支持该关联。提示卵巢衰老相关DDR基因变异可能影响卵母细胞基因组稳定性。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： 全外显子组测序、基因负担检验、BOLT-LMM、PheWAS、多基因评分、新生突变分析
• 公开数据： GSE86146、GSE107746
• 代码： 公开仓库 https://github.com/hadley/reshape、https://github.com/rgcgithub/regenie；本地已克隆：reshape、regenie
• 谁该点开原文： 从事生殖遗传学、癌症遗传学或新生突变研究的生信与遗传学研究者。

带走一句

罕见变异负担分析能发现常见GWAS遗漏的大效应基因，且与癌症表型组和新生突变数据整合可揭示卵巢衰老的跨表型遗传效应。