IF50.5｜精子测序揭示男性生殖系广泛正选择

精子中大量驱动突变在老年男性中富集，显著增加后代致病突变风险。

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

研究背景

人类细胞在生命周期中不断积累突变。在增殖组织中，获得性驱动突变可赋予克隆竞争优势，促进特定干细胞克隆扩增。癌症中此类选择已被广泛研究，但正常组织衰老过程中的克隆选择近年才受到关注。精原干细胞是唯一能将突变传递给后代的成体增殖细胞，其突变率极低，但靶向测序已发现精原干细胞可发生驱动突变并沿生精小管扩增，导致突变克隆在精子中比例升高。已知的男性生殖系驱动基因均为激活型错义热点突变，且多与严重发育障碍相关，可导致散发疾病出生率显著升高并与父亲年龄相关。

科学问题

卡在哪里： 由于睾丸和精子的多克隆性及极低突变率，此前研究仅限于少数已知基因，无法全面刻画男性生殖系正选择的范围和机制。

本文要回答： 利用高精度双链测序技术NanoSeq，在全外显子组水平系统鉴定男性生殖系正选择基因，并量化其对精子致病突变负荷的影响。

技术路线

作者主要用了：全基因组NanoSeq、全外显子组NanoSeq、靶向NanoSeq、dNdScv选择分析。

作者对81份精子样本进行全基因组NanoSeq测序，确定突变负担和特征；随后对38份样本进行全外显子组NanoSeq测序，对81份样本进行靶向NanoSeq测序，深度检测编码区突变；利用改进的dNdScv模型计算dN/dS比值，在全基因组范围筛选正选择基因和热点；最后通过整合变异等位基因频率估计精子中致病突变的比例。

图 1　Fig. 1。图1展示了精子与血液中突变负担的对比。a和b显示精子中替换和插入缺失随年龄线性增加，斜率与父源新生突变和睾丸数据一致。c和d显示血液中突变负担远高于精子。e显示血液与精子突变率的比值，替换约7.6倍，插入缺失约6.3倍。f和g显示突变特征，精子主要为SBS1和SBS5，血液还包含SBS19。

核心亮点

亮点 1｜精子突变率远低于血液

全基因组NanoSeq测序显示，精子中每个单倍体基因组每年积累1.67个替换突变（95% CI 1.41-1.92），与家系父源新生突变和睾丸精原细胞估计一致。血液中每个二倍体基因组每年积累19.9个替换突变，约为精子的7.6倍。精子突变特征主要为SBS1和SBS5，与年龄相关，而血液还包含SBS19。这些结果表明生殖系受到保护，突变率显著低于体细胞。

图 2　Fig. 2。图2展示了正选择基因的鉴定结果。a和b显示基因和热点水平的显著性，c显示突变类型富集和基因功能注释，d和e展示CUL3和SMAD4的突变分布。f显示dN/dS比值随年龄增加的趋势。g显示高表达于精原干细胞的基因dN/dS比值更高。h显示精子突变与癌症和发育障碍数据库的重叠。

亮点 2｜发现31个新正选择基因

通过dN/dS分析，共鉴定出40个显著正选择基因，其中31个为新发现。这些基因不仅包括激活型错义突变，还包括大量功能丧失（LOF）突变，涉及RAS-MAPK、WNT、TGFβ-BMP信号通路以及表观遗传修饰和RNA代谢等过程。新基因中16个为癌症基因，27个与单基因发育障碍相关。

图 3　Fig. 3。图3展示了致病突变和驱动突变在精子中的比例。a显示致病突变中由中性模型、驱动突变和未知因素解释的比例。b显示驱动突变比例随年龄增加。c显示致病突变比例随年龄增加且高于预期。d显示少数基因贡献了大部分致病或驱动突变。

亮点 3｜正选择增加致病突变负荷

精子中致病突变比例随年龄增加，从30岁的约2%上升到70岁的约4.5%，是中性模型预期的2-3倍。其中约三分之一由已知驱动突变解释，三分之一由中性突变解释，剩余三分之一尚未找到驱动基因。驱动突变比例也从30岁的0.5%上升到70岁的2.6%。

图 4　Fig. 4。图4展示了不同来源生殖系变异的dN/dS比值比较。a显示精子、健康家系、发育障碍家系和gnomAD不同等位基因频率的dN/dS比值，精子中正选择信号强，群体中负选择占主导。b显示正选择基因在gnomAD中的LOF z分数，大多数基因LOF缺失，但四个基因LOF富集。

亮点 4｜生殖系选择与体细胞选择共享机制

精子中富集的突变类型与癌症和发育障碍中的热点高度重叠，COSMIC高频突变富集11倍，发育障碍高频突变富集66倍。但SMAD4在精子中的热点是Myhre综合征相关，在癌症中少见。

亮点 5｜群体水平负选择占主导

比较不同来源的dN/dS比值显示，精子中正选择信号强，而gnomAD群体中高频变异dN/dS比值下降，表明负选择在世代间占主导。四个基因（SMAD6、MIB1、LZTR1、SSX1）在群体中LOF富集，可能反映生殖系正选择输入增加。

生信可带走

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

• 方法栈： 全基因组NanoSeq、全外显子组NanoSeq、靶向NanoSeq、dNdScv选择分析
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事生殖遗传学、体细胞进化或罕见病研究的科研人员。

带走一句

生殖系突变分析需考虑正选择影响，否则会扭曲群体遗传学推断和疾病基因发现。