自然同域河流中 Rhynchocypris lagowskii 与 R. oxycephala 的杂交格局：基于全基因组SNP的F1与回交鉴定及种界维持机制线索

《Ichthyological Research》：Hybridization patterns between Rhynchocypris lagowskii and R. oxycephala in a naturally sympatric river

　　

密切近缘物种间的杂交为理解生殖隔离与物种形成提供线索。研究人员调查了两种鲦类——Rhynchocypris lagowskii与Rhynchocypris oxycephala——在日本中部滋贺县Seri River自然同域范围内的杂交情况。线粒体DNA（mitochondrial DNA, mtDNA）系统发育证实，该河中两物种均为本地种群。利用MIG-seq进行的全基因组单核苷酸多态性（single nucleotide polymorphism, SNP）分析揭示两物种间存在清晰遗传分化，其中R. lagowskii在下游占优，R. oxycephala在上游占优。三角图分析鉴定出F₁杂交个体（1.1%）以及向R. lagowskii（4.8%）和R. oxycephala（12%）的回交个体，且均局限于中游。在本研究的采样设计与分析条件下，未检测到F₂或更晚世代个体及线粒体渐渗，提示潜在的合子后隔离（postzygotic isolation），而空间隔离可能促成合子前隔离（prezygotic isolation）。这些发现证实杂交在自然同域中发生但种界得以维持，使该系统成为研究物种形成的宝贵模型。

研究背景：密切近缘且分布重叠的物种是研究生殖隔离与物种形成的重要材料。Rhynchocypris lagowskii与R. oxycephala广泛分布于东亚，在日本本州中部自然分布重叠并可形成同域（sympatry）。两物种在分类上独立，但野生杂交长期被怀疑，依据是形态性状与mtDNA类型不一致，且多数涉及引入个体；自然同域河流中的同工酶（allozyme）研究结论冲突，一项报告频繁杂交，另一项未发现杂交证据。此前缺乏可靠基因组证据确认自然同域本地种群间是否发生杂交，以及强生殖隔离是否维持种界。因此，研究人员选择琵琶湖支流Seri River开展研究。该研究的意义在于确认自然同域中杂交发生但种界维持，为研究生殖隔离与物种形成提供模型系统。

研究人员开展的工作与结论：研究人员整合线粒体细胞色素b（Cyt b）系统发育与MIG-seq全基因组SNP，先确认Seri River种群本地性，再分析遗传结构、空间分布及杂交类别。结论是两物种在自然同域中发生杂交，但杂交频率低且种界未融合；杂交个体限于中游；未检测到F₂及以后世代和mtDNA渐渗；空间隔离可能促成合子前隔离，合子后隔离可能参与。其意义在于将该杂交系统确立为研究生殖隔离与物种形成的有效模型。

关键技术方法（≤250字）：样本队列来源为日本124个地点466个体；其中367个体来自94个地点用于Cyt b分析，293个体来自37个地点用于MIG-seq分析；Seri River详细分析192个体、6个站点（2022年10—11月，下游至上游），最终Seri River MIG-seq数据集为186样本。方法包括Cyt b的PCR与Sanger测序、723 bp比对及最大似然系统发育；MIG-seq建库测序，比对R. lagowskii参考基因组，Stacks调用SNP并行连锁不平衡（LD）剪枝；采用邻接网络（neighbor-net）、ADMIXTURE模型聚类（K = 1–6）和三角图分析（triangle plot analysis；主要阈值δ ≥0.9，并以δ ≥0.7、0.8、1.0进行敏感性分析）。

研究结果：
1. 遗传数据集概况（Overview of genetic datasets）：通过Cyt b获得364个目标样本序列，并整合GenBank中284条两物种序列及5条外群序列，其中日本149条序列来自92个地点；MIG-seq成功对37个地点287个样本进行基因分型。全数据集含360个SNP，经LD剪枝后216个SNP用于邻接网络；Seri River数据集含186个样本、6个站点、1,260个SNP，经LD剪枝后576个SNP用于模型聚类和杂交分析。
2. 线粒体DNA系统发育与Seri River种群本地性（Mitochondrial DNA phylogeny and native status of the Seri River population）：基于723 bp Cyt b最大似然系统发育，共识别257个单倍型；R. lagowskii与R. oxycephala单倍型总体形成两个高支持支系，但Rlag 3支系嵌套于R. oxycephala支系内。日本分布支系包括R. lagowskii的Rlag 2以及R. oxycephala的Roxy 3和Roxy 4。Seri River样本单倍型与琵琶湖周边单倍型匹配，提示该河Rhynchocypris鱼类发生人工引入的可能性低，支持其本地种群属性。
3. Seri River种群与其他地区种群的遗传关系（Genetic relationships between populations in the Seri River and other regions）：基于216个SNP的邻接网络显示，Seri River中存在R. lagowskii、R. oxycephala及遗传中间体；三大簇对应Rlag 2、Roxy 3和Roxy 4，异域样本形成清晰参考簇。类似中间体亦见于其他琵琶湖支流Kusano River和Tamura River。
4. Seri River沿河种群结构（Population structure along the Seri River）：ADMIXTURE在K = 2时交叉验证误差最小，两个祖先成分对应两物种。R. lagowskii个体倾向分布于下游（Sts. 1、2、4），R. oxycephala个体倾向分布于上游（Sts. 2–6）；遗传中间体仅确认于中游（Sts. 2–4）。该空间格局与既往关于两物种共存河流中R. lagowskii偏下游、R. oxycephala偏上游的发现一致。
5. Seri River种群内杂交（Hybridization within the Seri River populations）：三角图分析基于83个祖先信息SNP，强支持两物种在Seri River发生杂交。遗传中间体包括F₁杂交个体（杂交指数约0.50，类间杂合度接近1.00）及向任一亲本种的回交个体。F₁占1.1%（2/186），向R. lagowskii和R. oxycephala的回交分别占4.8%（9）和12%（23）；在本研究采样设计与分析条件下未观察到F₂及以后世代个体，且结果对δ阈值选择稳健。所有鉴定为纯R. lagowskii或向其回交的个体均携带R. lagowskii mtDNA单倍型；所有纯R. oxycephala或向其回交的个体均携带R. oxycephala mtDNA单倍型。

讨论与结论总结：研究人员认为，R. lagowskii与R. oxycephala即使在其自然同域范围内也会杂交，但种间杂交频率低，两物种未表现出融合。杂交个体仅见于河流中游，提示合子前隔离可能作为栖息地偏好和空间隔离的副产物建立。由于在本研究条件下未检测到F₂或更晚世代个体，也未观察到mtDNA渐渗，研究人员提出合子后隔离（如F₁雌性不育或线粒体—核不亲和）可能发挥作用；但这些机制仍属推测，需实验验证，尤其考虑到研究在时空采样、样本量及杂交类别判定方面的局限。其种界可能由这些生殖隔离机制组合维持。两物种在Seri River之外许多河流中共存，开展类似调查将有助于验证Seri River中观察到的杂交格局是否在河流系统间一致。R. lagowskii与R. oxycephala的杂交系统可作为研究生殖隔离与物种形成的有效模型。研究结论（译）：综合结果，R. lagowskii与R. oxycephala在自然同域中发生杂交，但杂交频率低且种界维持；杂交局限于中游，空间隔离可能促成合子前隔离，未检测到F₂及以后世代和mtDNA渐渗提示潜在合子后隔离；种界可能由多种生殖隔离机制组合维持，该系统可作为研究物种形成的宝贵模型。