非洲鳢Parachanna obscura沿物种形成连续体的基因组证据：种群分化与隐性多样性

《Global Ecology and Conservation》：Genomic evidence along the speciation continuum of Parachanna obscura (Channidae): population differentiation and cryptic diversity

　　

物种形成是一个进化过程，祖先物种通过外在（环境）和内在（如遗传）因素的逐渐积累而分裂成两个或多个独立谱系，这些因素可能最终阻止亚群间的基因流并导致生殖隔离。传统上根据地理隔离程度分为异域、邻域和同域模式，现越来越多地被视为由空间、生态和生殖过程相互作用塑造的物种形成连续体。淡水鱼类由于河流环境的碎片化和流域间潜在的隔离，为研究这些过程提供了特别合适的系统。非洲蛇头鱼Parachanna obscura（鳢科）广泛分布于西非和中非多个水文孤立流域，是探索这些过程的理想模型。本研究利用全基因组重测序技术鉴定的基因组范围变异，调查了P. obscura在西非12个河流流域的种群结构、进化历史和多样化格局。研究人员发现了五个遗传分化的种群，其遗传分化受限于异构水文网络间的有限连通性、地理距离和环境异质性。人口统计建模支持中更新世时期河流网络逐步破碎化驱动的分歧。此外，伊莫河（Imo River）的一个深度分歧谱系表明P. obscura复合体内可能存在一个推定的隐性物种。这些发现强调了历史水文气候过程在塑造西非河流系统生物多样性格局中的作用，并为该地区淡水鱼类的未来保护和进化研究提供了基因组框架。

研究背景

物种形成是进化生物学中的核心议题，传统上依据地理隔离程度被划分为异域、邻域和同域三种模式。然而，现代观点倾向于将其视为一个由空间、生态和生殖过程共同作用的“物种形成连续体”。在这一过程中，地理隔离长期被认为是驱动物种形成的首要因素。当种群被物理屏障分隔后，基因流中断，遗传和表型差异逐渐累积，直至达到生殖隔离。在空间结构化的种群中，“距离隔离”（IBD）可使连续种群最终碎裂为生殖隔离的群体；而在异质性景观中，局部适应会进一步加剧这种分化并促进物种形成。值得注意的是，物种形成并不总是需要完全隔离，在邻域物种形成中，相邻种群在不同选择压力下发生分歧的同时仍可维持低水平基因交流。

淡水鱼类因其栖息地——河流环境的碎片化和流域间的潜在隔离特性，成为研究上述过程的理想系统。流域常被集水区边界、水流变化（包括高速急流段）、瀑布和季节性干涸等物理屏障分隔，从而限制扩散和基因流。这些条件促进了异域和邻域分歧，而人为改变（如水坝建设、土地利用、水资源开采和污染）则可能进一步加剧碎片化并改变选择压力，加速分化进程。在进化时间尺度上，这些相互作用的力量塑造了精细尺度的种群结构，并可导致隐性物种形成。

非洲淡水生态系统拥有约3300种鱼类，是世界上物种最丰富的鱼类区系之一。这一显著生物多样性分布在一个动态的水文景观之上，其特点是内陆三角洲、古湖泊、急流和河流袭夺事件。上新世和更新世期间的气候波动，以湿润和干旱期的交替为特征，深刻重塑了非洲的地貌和水生生境。这些环境波动导致了反复的种群隔离、局部灭绝和重新定殖，被认为加速了遗传分化并推动了物种形成。然而，尽管有这一总体框架，很少有研究明确探讨历史气候振荡如何塑造西非河流系统中种群的碎片化模式和随后的异域分歧。该区域拥有超过550种淡水鱼类，其特征是众多半隔离的“生物地理岛屿”状流域，受流域面积、水文连通性和历史障碍的共同制约。

非洲蛇头鱼Parachanna obscura（鳢科）是探索这些过程的合适模型。该物种广泛分布于西非和中非，但不分布在北非。在西非，它出现在多个水文孤立的流域中。其广泛的分布跨越众多水文孤立的河流盆地，自然限制了基因流，而流域间的环境变异则创造了多样的选择景观。这些条件可能促进了距离隔离和环境隔离，并结合了上新世和更新世非洲历史气候事件中可能的避难所隔离。因此，可以预期强烈的种群结构甚至出现反映隐性物种的遗传隔离种群。一项近期研究揭示了主要来自科特迪瓦采样点的P. obscura种群间存在显著的遗传分歧，可能受到地理隔离、空间距离和当地环境条件的影响。然而，该研究使用的微卫星标记数量有限以及样本来源的空间范围受限，限制了其在更广空间尺度上解析种群结构和进化历史的分辨率。

研究内容与结论

为克服上述局限性，研究人员利用全基因组重测序发现的基因组范围变异，调查了P. obscura在西非12个河流流域的种群结构、进化历史、种群统计学历史和多样化格局。研究人员首先通过估算遗传分化、基因流和种群混合，并评估地理距离和环境异质性在塑造遗传分歧中的作用，研究了种群分化和结构的模式。研究人员假设西非复杂的地质历史和水文网络促进了流域间的强烈分歧，并得到当地环境条件的进一步加强。随后，研究人员利用人口统计建模重建了种群统计学历史，假设西非所有P. obscura种群起源于一个共同的祖先种群，该祖先通过历史上的水文连接殖民了不同的流域，最终形成了现今孤立的种群。最后，研究人员通过估算谱系分歧指数（GDI）来评估遗传簇间的物种界定，探索了潜在物种形成的信号。研究人员假设长期隔离极大地限制了西非河流复杂水文网络中种群间的基因流，可能推动新物种的出现。

研究结果表明，种群内基因组多样性普遍较低，而流域间遗传分化通常较高，但在水文相连或地理邻近的系统间分化较低。基因流估计值证实了这一模式，在地理接近或相连的地点间观察到更高水平的基因流。距离隔离分析显示地理距离对遗传分化有显著影响，而环境隔离测试表明环境异质性也通过空间梯度对种群分化有所贡献。人口统计分析揭示了一个层级式的分歧历史，涉及祖先种群的收缩，随后是西非河流系统的渐进式碎片化，所有这些都发生在距今约16.1万至26.8万年的中更新世。这一模式支持了阶梯式碎片化而非定向定殖的情景，即一个最初广泛分布的种群逐渐被细分到孤立的河流谱系中。最后，一个种群与其他种群高度分化，被识别为一个推定的隐性物种，表明P. obscura实际上可能代表一个物种复合体。

主要关键技术方法

研究人员采集了来自西非五国12个河流流域30个采样点的127份P. obscura标本。采用MGI DNBseq平台进行全基因组鸟枪法测序。将测序读段比对至P. obscura染色体级别参考基因组（GCA_060808335.1，来源于Bia河的KRIN种群）。使用Platypus进行SNP鉴定，经滤获得3,439,341个高质量变异（含2,873,386个SNP和565,955个Indel）。同时提取完整的线粒体基因组序列（16,114 bp）用于独立分析。使用Popcluster进行种群结构分析，PLINK进行主成分分析（PCA），dartR计算遗传分化指数（F_ST）和基因流（Nm）。运用Mantel检验评估距离隔离（IBD）和环境隔离（IBE）。利用fastsimcoal2基于多维位点频谱（SFS）推断种群统计学历史。使用IQ-TREE2构建最大似然系统发育树，并通过贝叶斯系统发育与生物地理学（BPP）程序计算谱系分歧指数（GDI）以评估物种界定。

研究结果

3.1 测序、变异鉴定与过滤

鸟枪法重测序结果显示高质量数据，目标测序深度为10-20X。每个样本的干净读段数约为3600万至9300万以上，采用双端150 bp（PE150）测序。Q20大多高于99%，Q30介于96.2%至97.4%。GC含量相对一致，范围为41.2%至43.0%。比对效率高，大多数样本比对率约96%至99%。初始变异数据集包含10,636,855个变异（8,038,526个SNP和2,598,329个Indel），最终数据集包含3,439,341个高质量变异（2,873,386个SNP和565,955个Indel）。

3.2 种群结构与分化

30个采样点的观测杂合度（Ho）和期望杂合度（He）显示变异性，但遗传多样性水平较低。Ho从0.05（MEDI）到0.14（DANGBO），He介于0.04（MEDI）和0.14（GPOPO）之间。近交系数（F_IS）从-0.66到0.20不等。多数种群F_IS为负值。大多数地点偏离哈代-温伯格平衡（HWE）不显著，但DOUN、KOUBI、KRIN、MEDI和IMO等地表现出大量偏离。基于Popcluster分析，所有个体被划分为五个遗传簇（K=5）。簇1包括来自科特迪瓦北部和布基纳法索西南部的个体，涵盖尼日尔河下游和沃尔特河上游流域。簇2包含科特迪瓦东南部、中部和东部阿涅比河、班达马河和科莫埃河流域的样本。簇3包括贝宁和多哥的韦梅河、莫诺河和沃尔特河下游流域的样本。簇4由科特迪瓦西部和西南部的个体组成，大部分采集于萨桑德拉河流域及其支流或布约湖。簇5几乎完全由尼日利亚伊莫河的个体组成。PCA分析与Popcluster结果一致，其中簇5（尼日利亚）在遗传上与其他所有群体最为不同。采样点间的遗传分化相对较高（>0.15），但最高分化（F_ST>0.20）出现在包含高度分化的IMO位点的配对中。同一流域内的种群间遗传分化低（F_ST<0.05）。分子方差分析（AMOVA）显示，37.0%的遗传变异源于流域间差异，22.4%源于流域内位点间差异，40.5%源于位点内差异。基因流估计显示，同一流域内或地理上接近但无水文连接的位点间基因流较高。

3.3 地理距离和局部环境条件对遗传结构的影响

距离隔离（IBD）得到确认，地理距离与遗传距离之间存在强正相关关系（Mantel r = 0.49, p = 0.001）。环境隔离（IBE）检测显示局部环境条件对遗传分化有显著贡献（Mantel r = 0.21, p = 0.008）。然而，控制地理距离后的偏Mantel检验显示无显著相关性（r = 0.006, p = 0.46），表明观察到的IBE信号主要由IBD驱动。

3.4 P. obscura的进化历史与轨迹

在30个采样点和河流流域中鉴定出120个mtDNA单倍型。大多数地点单倍型多样性高（24个地点Hd=1.0），但核苷酸多样性普遍较低。单倍型网络显示一个高度互联的核心结构，中心单倍型源自萨桑德拉流域，可能是祖先节点。伊莫河流域的单倍型形成最独特和分歧最大的谱系。在fastsimcoal2测试的人口统计场景中，最佳支持模型显示一个层级式的分歧历史，涉及祖先种群收缩，随后是西非河流系统的渐进式碎片化。推断的种群多样化发生在约16.1万至26.8万年前。最早的分歧对应伊莫河种群与祖先种群的分裂，约发生在268,481年前。后续分歧事件发生在16.1万至16.5万年前的狭窄时间区间内。

3.5 物种形成信号的调查

身份矩阵显示，先前定义为簇5的个体间相似性最低（约98%），表明这些种群可能代表一个不同的物种。基于mtDNA和核基因组数据的最大似然系统发育树均支持簇5形成最长分支，遗传距离最大。谱系分歧指数（GDI）密度图显示，连接（簇1+簇2+簇4+簇3）与簇5的节点GDI值最高，均值达0.98，峰值接近1.0，强烈表明簇5（伊莫河）代表一个不同的物种。相比之下，连接簇1至4的节点GDI值大多低于0.15，表明它们属于同一物种。与GenBank中其他Parachanna物种的COI序列比较分析显示，伊莫河样本形成一个独立分支，在系统发育上最接近P. obscura，但遗传上与之不同。

讨论与结论

讨论部分总结了研究的主要发现：P. obscura种群内基因组多样性低，流域间遗传分化高，地理距离和环境异质性共同塑造了遗传结构。人口统计分析支持中更新世时期河流系统逐步碎片化驱动的层级式分歧历史。最重要的是，伊莫河种群被确认为一个推定的隐性物种，GDI、身份矩阵和系统发育树均支持其为独立的进化谱系。这一发现与此前基于COI基因的研究一致，并得到了核基因组数据的进一步证实。伊莫河流域与邻近水系间的地质屏障（如Udi–Okigwe–Awgu陡崖）以及其独特的生态和生物地理特征，共同促进了该谱系的长期隔离和独立进化。研究还讨论了种群结构模式与西非古地理历史（如苏鲁河的古连通性）以及沿海泻湖环境可能提供的季节性基因流之间的联系。

研究结论指出：这项研究强调了历史水文气候过程在塑造西非河流系统生物多样性格局中的作用，并为该地区淡水鱼类的未来保护和进化研究提供了基因组框架。大多数地点观察到的低遗传多样性表明，迫切需要适当的管理策略和保护行动以确保该物种的长期可持续性。检测到的高水平遗传分化，特别是来自不同或地理孤立流域的种群间分化，表明了持续的异域进化，凸显了将这些遗传上独特的种群视为进化重要单元（ESUs）并加以保护的重要性。重要的是，结果还揭示了P. obscura复合体内存在一个未表征的隐性谱系（伊莫河），这可能需要提升到物种级别，强调需要进行全面的分类学修订，结合基因组、形态学和生物地理学数据来阐明Parachanna属内的物种界限。