南大西洋海龟种间基因流与潜在杂交衰退：基于种群基因组学的证据

《Molecular Ecology》：Interspecific Gene Flow and Potential Hybrid Breakdown in South Atlantic Sea Turtles

　　

种间杂交（interspecific hybridisation）可在进化中发挥重要作用，既可能产生新分类单元，也可能导致远交衰退（outbreeding depression），并可通过强化（reinforcement）引发行为生殖隔离（behavioural reproductive isolation）。然而，关于分歧时间如何与杂交速率相关及其对基因渗入（introgression）的影响，研究仍较少。海龟尽管具有较长的分歧时间以及行为、生态和形态差异，仍已知可发生杂交。若干遗传学研究试图确定七种现存海龟中的四种（Caretta caretta、Eretmochelys imbricata、Lepidochelys olivacea 和 Chelonia mydas）在巴西大西洋沿岸的杂交程度，该区域是这四种海龟的筑巢和觅食地。在本研究中，研究人员整合了这四种海龟及其杂交个体的新生成和已发表基因组数据，开展详细的种群基因组学（population genomics）研究。尽管杂交确实在巴西筑巢地发生于四种中的三种之间，但似乎仅限于第一代杂交个体和少量回交，提示种间存在某些合子后屏障（postzygotic barriers）和杂交衰退（hybrid breakdown）。此外，证据提示杂交个体间具有中等偏高的亲缘关系，表明独立杂交事件数量减少。尽管如此，研究人员发现涉及不同物种的回交证据，以及 Ca. caretta x E. imbricata F1 杂交个体之间的 F2 子代，提示杂交雄性可育。有限的种间基因流似乎是海龟长期分歧历史的一部分，因为推断出低水平的深时基因渗入。总之，本研究为大尺度分歧下的杂交速率和成功率提供了重要证据，并表明这似乎与低基因渗入率相对应，指向基因流的时间限制。

该研究聚焦南大西洋海龟的种间杂交、基因流及其长期进化后果。背景上，杂交在脊椎动物中并不罕见，估计可见于多达10%的物种；杂交后代可能携带有害等位基因组合而适合度降低或育性受损，成功回交才可能形成基因渗入（introgression，一个物种遗传物质进入另一物种基因组）。基因渗入取决于杂种育性、杂种与亲本物种的相容性及繁殖成功率；强生殖屏障下杂交也可能仅限第一代。深分歧谱系间杂交是否伴随显著基因渗入并改变物种进化轨迹尚不清楚。海龟现存7种，共同祖先约2.2亿年前出现，种间分歧约10—80 My；其生活史长、世代时间长、进化速率慢、基因组高度保守，且具有长距离迁徙和雌性高回归性。巴西东部海岸自20世纪80年代起实施保护管理，此前研究已记录 Ca. caretta、Ch. mydas、E. imbricata 和 L. olivacea 四种海龟间广泛杂交，但多为F1，少数回交和F2，且既往研究难以推断特定杂交组合和长期杂交历史。该研究旨在整合大规模基因组数据和高品质参考基因组，检验分歧时间与杂交/基因渗入的关系、杂交方向性及杂交是否等同于基因流。

关键技术方法：研究人员采用简化基因组测序（reduced-representation sequencing），新生成3RAD序列来自15尾 Ca. caretta、1尾 Ch. mydas、7尾 E. imbricata、23尾经分子鉴定杂交个体和46尾未鉴定个体，样本来自巴西六州、Eleva??o do Rio Grande（ERG）觅食区及Guadeloupe；结合已发表RADseq/ddRAD数据，总样本167尾，覆盖Martinique、French Guiana、Espírito Santo和Abrolhos等地点，并以 D. coriacea 为外群。序列比对 D. coriacea 参考基因组，经bcftools call和过滤构建All与LD_pruned数据集；使用PCA、ADMIXTURE和triangulaR分析混合与杂交世代，dadi拟合人口模型，Dsuite F_branch检测基因渗入。

研究结果：3.1 参考基因组比对与数据过滤：通过插入长度、覆盖位点重叠和映射深度评估数据整合可行性，发现RAD文库片段大小重叠，28,253,981个唯一位点（36%）重叠，167尾样本平均映射深度≥5，且不同参考基因组间深度无显著差异；过滤后All数据集含1,065,048个位点，193,233个双等位SNP，LD_pruned数据集含20,675个SNP。3.2 混合分析：PCA可分离四个物种，杂交个体位于亲本物种之间，PC1和PC2分别解释28.5%和19.9%方差；ADMIXTURE最佳K=5，K=4时四物种分离，K=5时 Ca. caretta 分为两个祖先群；共鉴定49尾杂交个体，未发现涉及 Ch. mydas 或ERG来源 Ca. caretta 的杂交。三角图显示 Ca. caretta 与 E. imbricata、Ca. caretta 与 L. olivacea 产生F1及第一代回交，Ca. caretta 与 E. imbricata 中还发现潜在晚期回交、F2及F2与 E. imbricata 的潜在回交；E. imbricata 与 L. olivacea 仅见F1；所有回交个体均为幼体。35尾杂交个体与至少另一尾样本亲缘系数≥0.25；纯样本中仅 L. olivacea 一个三样本家系相关。多样性以 Ch. mydas 最高，其次 E. imbricata、L. olivacea 和两个 Ca. caretta 群体；各物种Tajima's D均为负。F_ST在 Ca. caretta 与 E. imbricata 间低于姐妹种 Ca. caretta 与 L. olivacea，但D_XY更符合系统发育分歧。3.3 人口模型：通过dadi比较五种模型，所有物种对最高似然多为二次接触不对称基因流模型；连续对称基因流显著优于无基因流，连续不对称基因流除 Ch. mydas × L. olivacea 外优于简单模型。不对称模型中分歧时间与系统发育估计一致，迁移率最高为 Ca. caretta 与 L. olivacea，方向由 Ca. caretta 至 L. olivacea；所有有效迁移率N×m远小于1，提示基因流很受限；有效种群大小与多样性估计一致。3.4 基因渗入：Dsuite F_branch检测到 Ca. caretta 与 E. imbricata 间约1.6%基因渗入、E. imbricata 与 Ch. mydas 间约0.2%基因渗入，其他物种对不显著。

讨论部分总结：4.1 南大西洋海龟的杂交：研究支持南大西洋成年海龟中杂交频繁，涉及分歧约13.6—62.8 Mya的 Ca. caretta、E. imbricata 和 L. olivacea；Ca. caretta 分为两个遗传群，第二群见于ERG觅食区。杂交主要限于前几代，回交较少且均为幼体（部分死胎），提示杂交衰退；F2个体支持杂交雄性可育。Ca. caretta 与 E. imbricata 杂交最普遍，尤其在Bahia的Praia do Forte和Abrolhos，杂交雌性筑巢成功率较低；高杂交率若长期存在可能预期强化，但数据尚不能区分近期杂交或未发生强化。线粒体数据提示性别偏倚，但采样偏差不能排除；Ca. caretta 与 L. olivacea 的回交证据提示F1可育至少部分可育，缺乏成年F2或回交可能提示有害或采样偏差，需谨慎。4.2 古老基因流：人口模型显示四物种分化以来存在一定程度基因流，迁移率随分歧时间下降，可能反映基因流的时间限制；长期基因渗入仅见于 Ca. caretta 与 E. imbricata、E. imbricata 与 Ch. mydas，且量很低，提示屏障可能较快建立；姐妹种 Ca. caretta 与 L. olivacea 虽有F1和回交但无长期基因渗入证据，支持当前杂交衰退假说；样本量、简化基因组和低基因流模型选择需谨慎解释。4.3 遗传多样性的对比情景：保护努力下海龟种群恢复，反映于有效种群大小和Tajima's D；范围扩张可能增加杂交概率，若强化不发生则可能影响恢复。Ch. mydas 多样性最高，E. imbricata 次之；低多样性是长期特征还是近期基因组侵蚀尚需研究。4.4 局限性：简化基因组测序限制SNP数量并可能受等位基因丢失和SFS偏差影响；Ch. mydas 和 L. olivacea 样本量小，地理上Bahia和Sergipe代表性较强；海龟慢进化与基因组保守性限制结果外推。结论翻译：总体而言，这些结果指向不同海龟谱系间历史和当代杂交事件，并可能因杂交衰退降低已濒危种群的生存力；结果还显示杂交雄性与雌性间成功交配的潜在证据；这些发现对南大西洋海龟种群保护和管理具有区域意义，重叠筑巢区和种间杂交可能导致局部下降和遗传侵蚀。