IF50.5｜减数分裂基因常见变异塑造人类重组与非整倍体风险

大规模PGT数据揭示SMC1B等减数分裂基因变异同时调控交叉重组与母源非整倍体。

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

研究背景

非整倍体是导致人类妊娠丢失的首要原因，常源于女性减数分裂中染色体分离错误。已知异常交叉重组会增加非整倍体风险，但缺乏大规模数据来解析两者的共同遗传基础。减数分裂交叉在胎儿期建立，与姐妹染色单体黏连共同维持同源染色体配对，直至排卵时分离。黏连蛋白复合物（如SMC1B）对染色体联会和交叉形成至关重要。此前研究多基于存活非整倍体个体或小样本配子，统计效力有限。

科学问题

卡在哪里： 此前最大规模的人类家系重组研究未纳入非整倍体参与者，无法直接检验交叉与非整倍体的遗传关联。单细胞测序虽可同时检测交叉和非整倍体，但样本量小，难以评估变异性和共享遗传结构。

本文要回答： 利用大规模PGT数据，系统绘制人类胚胎中的交叉重组和非整倍体图谱，定量检验两者关系，并发现调控这些表型的常见遗传变异。

技术路线

作者主要用了：SNP芯片PGT、karyoHMM、GWAS、TWAS、eQTL、EMSA、ARG推断。

作者收集139,416个囊胚期胚胎的SNP芯片PGT数据及父母DNA，开发karyoHMM隐马尔可夫模型追踪亲本单倍型传递，同时识别非整倍体和交叉事件。先验证非整倍体发生率与母源年龄关系，再比较整倍体与非整倍体胚胎的交叉计数，随后对母源减数分裂非整倍体进行GWAS，锁定SMC1B位点，并通过eQTL、EMSA等功能实验解析因果变异，最后结合TWAS和跨表型关联揭示减数分裂基因的广泛作用。

图 1　Fig. 1。图1展示数据来源和分析流程：a为PGT样本采集示意；b显示通过单倍型追踪识别交叉和非整倍体；c显示非整倍体主要涉及母源同源染色体，且在15、16、21、22号染色体富集；d显示母源年龄与非整倍体率的二次关系，父源年龄无显著效应；e显示交叉在基因组上的分布；f比较整倍体与非整倍体胚胎的二体染色体交叉数，显示非整倍体胚胎交叉数显著减少。

核心亮点

亮点 1｜非整倍体胚胎交叉数显著减少

在46,861个整倍体胚胎中鉴定出3,809,412个交叉，母源平均49.30个，父源31.99个，与既往家系研究一致。比较发现，非整倍体胚胎的二体染色体上交叉数显著低于整倍体胚胎（边际均值差异0.105，P=4.64×10⁻¹⁵⁰）。该结果支持交叉减少（尤其是缺失）会增加减数分裂非整倍体风险。

图 2　Fig. 2。图2展示母源减数分裂非整倍体GWAS结果：a曼哈顿图显示Chr.22上rs6006737达到全基因组显著；b显示不同基因型下母源年龄与非整倍体率的关系，风险等位基因携带者非整倍体率更高；c区域图显示关联信号位于包含SMC1B的约120kbp单倍型内。

亮点 2｜SMC1B常见单倍型关联非整倍体

GWAS发现Chr.22上rs6006737（P=2.21×10⁻⁸）与母源减数分裂非整倍体显著相关，并在独立测试集中复制（P=0.033）。风险等位基因在40岁患者中每拷贝增加约1.65%非整倍体风险，且效应随母源年龄增大。该单倍型跨越SMC1B等四个基因，SMC1B编码黏连蛋白复合物组分，对减数分裂染色体黏连和重组至关重要。

图 3　Fig. 3。图3展示SMC1B位点的功能注释：a为减数分裂黏连蛋白复合物示意图；b显示rs6006737风险等位基因与SMC1B表达降低相关（MAGE数据）；c显示rs6006737与SMC1B eQTL精细定位SNP的连锁不平衡；d显示rs2272804位于SMC1B启动子CpG岛内，并落在ATF1结合基序中。

亮点 3｜风险单倍型是SMC1B的eQTL

rs6006737风险等位基因与SMC1B表达降低显著相关（MAGE中β=-0.429，P=4.68×10⁻¹⁸）。精细定位锁定rs2272804（后验包含概率>0.99），位于SMC1B转录起始位点上游144bp。EMSA实验显示，含替代等位基因的DNA探针与ATF1转录因子的结合亲和力降低超过3倍（KD：参考型56.62nM vs 变异型173.39nM，P=2.60×10⁻⁴），支持非编码顺式调控机制。

图 4　Fig. 4。图4展示非整倍体和重组GWAS前导SNP与生殖衰老表型的多效性关联：颜色表示关联方向，深浅表示显著性。rs6006737风险等位基因与初潮年龄增大、绝经年龄提前相关；多个重组相关位点（如C14orf39、CCNB1IP1、RNF212）也显示与非整倍体的名义显著关联。

亮点 4｜TWAS揭示其他减数分裂基因

TWAS鉴定出C14orf39（P=1.65×10⁻⁷）与母源减数分裂非整倍体显著相关，该基因编码联会复合体组分，此前已发现罕见变异导致不孕，常见变异与重组表型关联。此外，NCAPD2（P=2.16×10⁻⁵）虽未达全转录组显著，但编码凝缩蛋白I复合物调控亚基，参与减数分裂染色体凝集。这些结果提示常见非编码调控变异通过影响减数分裂基因表达调节非整倍体风险。

图 5　Fig. 5。图5展示SMC1B风险等位基因的进化分析：a显示rs2272804衍生等位基因的频率轨迹，在非洲人群中较高，近期在非非洲人群中略有下降；b显示理论模型中适合度代理与实际适合度的缩放因子α与选择系数s的关系，表明只有极弱的适合度关联才能解释该等位基因的常见频率。

亮点 5｜多效性关联生殖衰老表型

非整倍体风险等位基因rs6006737同时与女性重组率降低（β=-0.033，P=0.002）、初潮年龄增大（β=0.021，P=3.82×10⁻¹²）和绝经年龄提前（β=-0.047，P=2.06×10⁻⁴）相关，提示缩短的女性生殖窗口。此外，三个重组GWAS位点（C14orf39、CCNB1IP1、RNF212）也显示与非整倍体的名义显著关联，方向一致，表明这些基因在重组和非整倍体风险中具有共享遗传基础。

亮点 6｜风险等位基因的进化悖论

尽管SMC1B风险等位基因降低适合度，但其在全球人群中频率较高（gnomAD全球AF=0.44），在非洲人群中更高（AF=0.71）。基于祖先重组图的推断显示该等位基因起源于约91万年前，频率轨迹基本符合中性进化。理论模型表明，只有当适合度代理与实际适合度的关联极弱（α≤0.01）时，该等位基因才能通过遗传漂变维持中间频率，解释了常见有害变异的维持机制。

生信可带走

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

• 方法栈： SNP芯片PGT、karyoHMM、GWAS、TWAS、eQTL、EMSA、ARG推断
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事生殖遗传学、GWAS/eQTL分析、减数分裂机制研究的生信与遗传学研究者。

带走一句

大规模PGT数据可同时解析重组与非整倍体，揭示减数分裂基因的共享遗传基础；分析时需注意非整倍体胚胎交叉检测的偏倚，并考虑多效性关联。