小种群与大海种群刺鱼的重组率差异及关联选择效率：基于无缺口端粒到端粒参考基因组的多群体连锁图谱研究

《PLOS Genetics》：Recombination rate and efficiency of linked selection in small and large stickleback populations

　　

摘要专业翻译 种群遗传学理论预测，自然选择在大种群中比在小种群中更有效，因为遗传漂变降低了小种群中选择的效率。适应新环境的小种群也可预期进化出更高的重组率，以促进适应并解离和清除有害突变。研究人员通过（1）调查四个小淡水种群（平均Ne≈7 096）和四个大海种群（平均Ne=376 828）的九刺鱼（Pungitius pungitius）基因组中核苷酸多样性（π）与重组率之间关联强度的差异，以及（2）利用种群特异性连锁图谱比较小种群与大海种群之间的重组率，对上述假说进行了检验。研究人员发现，除最小的淡水种群外，所有种群均存在预测的π与重组率的正相关关系，表明即使在考虑了GC/CpG含量和基因密度变异之后，连锁选择仍然普遍存在。淡水种群与海洋种群之间的平均重组率无显著差异，唯独具有最小Ne的淡水种群表现出显著升高的重组率。探索性全基因组关联分析（GWAS）提示重组率可能具有多基因遗传基础。这些结果表明，连锁选择在降低低重组区域的π方面发挥重要作用，尤其是在大种群中。此外，正如理论所预测的，至少有一个小淡水种群似乎已进化出高于其海洋祖先的重组率。

论文解读文章

一、研究背景与科学问题

自然选择与遗传漂变是塑造基因组中性变异模式的两大核心进化力量。正向选择（选择性清除）和负向选择（背景选择）均可降低与选择靶标连锁区域的遗传多样性。由于选择效应在紧密连锁区域最为显著，理论上基因组中核苷酸多样性（π）应与重组率呈正相关。尽管这一预期在多种模式生物中已得到验证，但该相关性并非普遍存在，且针对非模式生物自然种群的研究仍然匮乏。尤其值得关注的是，小种群中遗传漂变会显著削弱选择效率，因此同一物种内小种群与大海种群之间的π-重组率关联强度差异，是检验选择效率理论预测的理想体系。

另一方面，重组率具有可遗传性且可进化。理论预测，适应新环境的小种群可能进化出更高的重组率，因为增加重组可产生更多新的等位基因或性状组合，并有利于解离和清除有害等位基因。然而，这种升高仅在特定条件下出现，包括小种群规模、强烈且持续的方向性选择、高有益突变率、顺式作用重组修饰因子以及选择位点间极低的重组等。这些条件在自然环境中是否满足尚不清楚，野外适应性重组率分化的经验证据极为稀缺。

此前同种不同种群间的重组率比较研究多为两两比较，且部分研究采用基于连锁不平衡的方法，易受种群历史和选择效应的混杂影响。因此，亟需开展包含多个种群、利用家系杂交连锁图谱（以避免种群历史和选择偏倚）的研究，系统比较不同有效种群大小（N_e）种群间的重组率差异。

本研究以九刺鱼（Pungitius pungitius）为研究对象，选取四个小N_e淡水种群和四个大N_e海洋种群，旨在检验两个核心预测：（1）π与重组率的正相关关系及其强度是否随N_e大小而变化；（2）小种群是否进化出高于大种群的重组率。该物种的小种群近期经历了强烈的方向性选择，增加了高重组率进化的可能性。研究利用高质量无缺口端粒到端粒（T2T）参考基因组以提高全基因组π和重组率估计的准确性。

二、主要关键技术方法

研究人员基于八个九刺鱼自然种群（四个淡水、四个海洋）的家系材料构建了十张高密度连锁图谱，涵盖158,980至3,838,162个SNP标记，基因组覆盖范围超过常染色体区域的99%。连锁图谱构建采用标准化Lep-MAP3流程，包括ParentCall2基因分型、系列过滤（dataTolerance=0.001、removeNonInformative=1、MAFLimit=0.05、missingLimit=0.1）、OrderMarkers2排序及遗传距离计算。对X染色体性别决定区（SDR）进行了母系信息标记的特殊处理。核苷酸多样性（π）基于F₀亲本SNP以1Mb非重叠窗口估计，使用VCFtools和Pixy两种方法交叉验证。基因密度和CpG含量基于T2T参考基因组注释和甲基化数据计算。重组率热点定义为重组率最高的前5%窗口，冷点定义为无重组检测的窗口。统计分析采用广义线性混合模型（GLMM）和广义线性模型（GLM），通过AIC准则筛选最优模型。全基因组关联分析（GWAS）采用rMVP软件中的GLM、MLM和FarmCPU三种方法，以亲本遗传变异与减数分裂交叉数为关联指标，前五个主成分作为协变量控制亲缘关系影响。为排除PYO种群低SNP密度对重组率估计的潜在影响，研究人员将其他七个种群的标记密度降至PYO水平并重复构建连锁图谱100次进行验证。研究样本涵盖四个淡水种群和四个海洋种群，其中TVA和POR种群的测序数据来自先前已发表研究。

三、研究结果

连锁图谱与雌性偏倚的异源交叉（heterochiasmy）

常染色体遗传图谱长度估计为1,556至1,815 cM。染色体水平上，性别平均遗传图谱长度与物理距离呈显著正相关（r_s=0.61, p=2.03e-16）。雌性遗传图谱长度在染色体水平上显著长于雄性（配对t检验，t₁₉₉=31.16, p=2.2e-16），淡水图谱平均雌雄比为1.73，海洋图谱为2.03。X染色体两性差异最为显著：雌性SDR遗传距离为106.5 cM（淡水）和119.7 cM（海洋），而雄性SDR无重组；PAR区域雌性图谱长度分别为雄性的2.00倍（淡水）和1.63倍（海洋）。

种群间全基因组重组率变异

单个后代检测到的交叉数从18到52不等，平均每次减数分裂母系21.1个、父系11.2个交叉。淡水种群平均交叉数显著高于海洋种群（t检验，t₂₅₂=-5.27, p=2.94e-07），其中PYO种群交叉数最高（均值40.30，95% CI: 38.90–41.70），显著高于其他所有种群（t检验，t_947.16=-8.97, p<2e-16），为交叉数最低的UME种群的1.36倍。种群水平交叉数与MSMC2估计的祖先N_e（10 kya）呈显著负相关（GLMM, β=-1.23, p=0.02），但与CurrentNe2估计的当代N_e无显著相关（GLMM, β=-0.92, p=0.098）。染色体重组率范围为1.22至7.18 cM/Mb，与染色体长度呈负相关（GLMM, β=-0.24, p=1.69e-11）。淡水与海洋种群的性别平均重组率无显著差异（MA=3.74, FW=3.91 cM/Mb; t检验, t_3.93=1.24, p=0.28）。淡水种群染色体水平重组率略高于海洋种群（MA=3.63, FW=4.07 cM/Mb; t检验, t_141.69=2.14, p=0.03），但去除PYO后差异不再显著（FW_non-PYO=3.72 cM/Mb, t检验, t_113.06=-0.61, p=0.55）。降采样验证表明，低标记密度仅使重组率估计膨胀约0.05 cM/Mb，远不足以解释PYO的高重组率。

局部重组率保守性

八个种群的遗传图谱在全基因组范围内变化模式高度相似，1Mb窗口重组率在两种生态型间紧密对齐（r₉₀₃=0.74, p<2.2e-16），两性间亦显著相关（r₉₀₄=0.54, p<2.2e-16）。重组热点倾向于分布在染色体端粒区域，这些区域CpG含量显著高于冷点（t_40.93=-14.39, p<2.2e-16）。重组率与CpG含量呈强正相关（负二项GLMM），且存在显著的性别特异性效应；基因密度高的区域重组事件较少，呈弱负相关。chr19上已知的倒位区域显著抑制了交叉的发生。

核苷酸多样性（π）与重组率的关系

基于种群和生态型特异性连锁图谱的最优GLM模型分别解释了π变异的88.77%和53.52%。由于海洋种群平均π高于淡水种群，生态型因素贡献了大部分方差（种群水平83.68%，生态型水平43.63%）。两种模型均表明局部π与性别平均重组率和CpG含量正相关，与基因密度弱负相关。π与重组率的正相关在淡水种群中显著更平缓（GLMM, β_{recom.rate:ecotype}=6.33e-05, p<2e-16），且本研究中最近交的PYO种群是唯一未表现出显著π-重组率关联的种群（β=2.91e-06, p=0.43）。按生态型分别建模时，两种生态型中重组率和CpG含量均为解释π变异的最重要固定效应，但淡水种群表现出远大于海洋种群的残差变异（FW: 34.6%, MA: 0.0%）。Pixy估计与VCFtools估计的π在各种群间高度相关（r=0.848至0.956）。

与交叉数相关的遗传位点

GWAS分析在不同方法中鉴定出约4至48个与全基因组交叉数相关的推定数量性状位点（QTL，p<1e-05）。FarmCPU方法鉴定出48个暗示性QTL，涉及26个功能未知的基因，以及多个已知候选基因如adprm、ncoa2、htra3a、mbtps1和slc25a26。六个QTL位于基因间区域。经Benjamini-Hochberg FDR校正（FDR<0.05）后，FarmCPU和GLM分析中各有三个QTL保持显著，MLM中无显著位点。

四、讨论与结论

本研究的核心发现包括：大海种群中核苷酸多样性与重组率之间存在一致且强烈的正相关，而小淡水种群中该正相关显著更平缓，PYO种群甚至完全缺失。这一模式符合遗传漂变削弱小种群选择效率的理论预期，得到多重佐证：淡水种群虽已清除大部分高有害突变负荷，但积累了相对较高的轻度有害突变负荷；淡水种群显著富集年轻活跃转座元件；淡水种群间分化程度和近交水平极高；其当代有效种群大小比海洋种群低千倍以上。然而，种群瓶颈也可能掩盖连锁选择信号，现有证据尚不能完全区分这两种解释。

PYO作为最小的淡水种群，其交叉数和重组率显著高于大海种群，提示该种群在约10,000年前定殖淡水环境后可能进化出了更高的重组率。这符合小种群适应新环境时通过增加重组产生新等位基因组合、促进有害等位基因解离和清除的理论预期。但其他淡水种群并未表现出升高的重组率，可能反映了种群遗传权衡：强遗传漂变降低选择效率，阻碍了重组率升高的适应性转变。值得注意的是，PYO基因组91%区域为纯合状态，在这种严重受损的遗传背景下，清除有害突变的迫切需求可能独特地克服了遗传漂变的限制。同时，长纯合片段可能降低重组事件检测概率，意味着淡水种群重组率可能被低估；若淡水种群基因组确实小于海洋种群，基于海洋参考基因组的物理距离换算还会进一步低估淡水种群的cM/Mb值。因此研究者在核心分析中采用不依赖参考基因组物理缩放的交叉数作为主要指标。

重组事件的分布受基因组结构制约。交叉干扰（CI）现象将交叉推向染色体末端，限制了每条染色体的最大交叉数，导致长染色体重组率较低。结构变异如chr19倒位也抑制交叉发生。重组率与CpG含量强正相关符合预期，因为九刺鱼的Prdm9基因已失去功能结构域KRAB和SSXRD，重组热点倾向于集中在富含CpG的启动子样区域，这与犬科和鸟类等prdm9非依赖物种的研究结果一致。GWAS结果进一步支持刺鱼重组率不受prdm9控制的观点：重组率呈多基因遗传特征，多数候选基因与哺乳动物中鉴定的基因不同。其中slc25a26直接参与甲基化过程和细胞增殖，klhdc3l在结构上类似于人类重组激活基因2（rag2）。

尽管研究采用了高质量T2T参考基因组和高深度测序数据，重复序列区域的短读段变异检测仍具挑战性，这些区域被排除在多数分析之外。未来整合长读段测序数据有望改进这些区域的变异检测，并提供转座元件和结构变异在影响重组模式和遗传多样性方面作用的新见解。

结论：利用高质量连锁图谱和改进的参考基因组，本研究在八个具有对比有效种群大小的刺鱼种群中揭示了重组率变异及其基因组基础。低重组区域核苷酸多样性显著降低，表明连锁选择在降低多样性方面发挥重要作用。小种群中π-重组率正相关弱于大种群，提示强漂变条件下选择效率降低。最小N_e种群表现出最高的重组率，可能为抵消遗传多样性降低和连锁不平衡升高的负面后果而进化出升高的重组率，但重组率差异与基因组大小和结构差异相关的可能性不能排除，需进一步研究。由于物理距离以较大的海洋参考基因组为基准，淡水种群的cM/Mb值如有偏差只会被低估，正确大小的淡水参考基因组将扩大而非缩小淡水-海洋差异。结果还表明刺鱼重组率变异具有多基因基础，但需进一步研究探索无prdm9物种中重组变异的详细分子机制和交叉分布的潜在分子驱动因素。