立陶宛淡水生态系统中拟鲤（Rutilus rutilus）和河鲈（Perca fluviatilis）的比较种群遗传学

《Animals》：Comparative Population Genetics of Rutilus rutilus and Perca fluviatilis in Lithuanian Freshwater Ecosystems

　　

理解共存淡水物种对同一水文景观的响应方式，对于解读遗传多样性和种群连通性具有重要意义。研究人员比较了两种广泛分布的淡水鱼类——拟鲤（Rutilus rutilus (Linnaeus, 1758)）和河鲈（Perca fluviatilis Linnaeus, 1758）——在立陶宛主要河流流域中的种群遗传结构，所用标记为线粒体ATP₆（ATP₆）和D-loop序列。两个物种均表现出较高的单倍型多样性（H_d），但遗传结构差异显著。R. rutilus显示出弱分化、流域间广泛的单倍型共享，以及符合种群扩张的人口统计学信号。相比之下，P. fluviatilis表现出更强、空间上更异质的结构化特征，人口统计学信号一致性较低，分化更具局部性。与预期相反，两个物种的遗传变异均未与主要流域边界持续关联。相反，观察到的模式表明，物种特异性生态性状（尤其是扩散能力和栖息地利用的差异）可能比当代水文边界单独发挥更大的作用，从而塑造种群结构。这些发现表明，即使共享同一水文系统，共存的淡水鱼类也可能表现出对比鲜明的遗传多样性模式。尽管核基因或全基因组标记将增强这些线粒体推断的可靠性，但本研究凸显了比较方法在理解种群连通性以及为淡水生物多样性保护和管理提供信息方面的价值。

解读文章：

**研究背景与科学问题**

淡水生态系统是全球破碎化程度最高的生态系统之一，其物种遗传多样性受水文连通性与历史景观格局的共同塑造。北欧地区在更新世经历多次冰川作用，水系网络和水生生境被反复重塑，这为许多淡水类群的遗传结构留下了长期印记。因此，当前遗传多样性格局既反映历史过程，也反映当代环境条件。然而，如何区分这两类因素的相对贡献，一直是淡水种群遗传学研究的核心难题。一种有效策略是比较同一水文景观中共存但生态性状不同的物种，因为共存物种经历了相似的环境和水文历史，若它们表现出遗传结构差异，则可归因于物种特异性生态而非共享历史。

立陶宛位于波罗的海流域，境内包含涅曼河（Nemunas）、利耶卢佩河（Lielupe）、文塔河（Venta）和道加瓦河（Daugava）四个主要流域，这些流域在较小地理范围内构成相对独立的单元，为同时评估水文边界的效应提供了理想设计。然而，以往对拟鲤（Rutilus rutilus）和河鲈（Perca fluviatilis）的研究多为单物种研究，在同一景观下用相同取样设计和分子标记直接比较二者的情况十分罕见，尤其缺乏波罗的海地区的线粒体数据。这一空白使得难以区分物种特异性生态与共享水文的作用。为此，研究人员利用线粒体ATP₆基因和控制区（D-loop），比较了两种鱼类在立陶宛主要流域的种群遗传结构，旨在检验两个互补假说：其一，若结构由共享水文驱动，则两个物种的遗传分化应与流域边界显著关联；其二，若物种特异性生态占主导，则即使共享同一水文景观，两个物种的遗传分化程度和空间模式也应出现差异。该研究发表于《Animals》。

**关键技术与方法**

研究人员从立陶宛各地淡水系统（包括河流、湖泊、水库及库尔斯泻湖）采集了209尾R. rutilus和191尾P. fluviatilis，样本涵盖涅曼河、道加瓦河、利耶卢佩河等流域，采集时间为2009至2025年。对肌肉或鳍组织提取DNA后，扩增并测序线粒体ATP₆和D-loop片段。分析手段包括：基于DnaSP计算遗传多样性指数，利用NETWORK构建中位连接网络，MEGA推断最大似然系统发育树，Arlequin计算成对F_ST和层级分子方差分析（AMOVA），并开展中性检验、错配分布分析和贝叶斯天际线图（BSP）重建种群历史动态。

**研究结果**

**3.1 序列特征与多态性**

通过对拼接的ATP₆–D-loop序列分析，R. rutilus在209条序列中鉴定出37个多态位点（含2个插入缺失），P. fluviatilis在191条序列中鉴定出33个多态位点（含4个插入缺失）。所有插入缺失均位于D-loop区域，P. fluviatilis中发现的异质长度变异（10–20 bp插入）因无法可靠比对而被排除。结果表明两个物种均以不变位点为主，R. rutilus的多态位点略多于P. fluviatilis。

**3.2 种群遗传学与结构**

遗传多样性指数显示，R. rutilus整体单倍型多样性为0.760 ± 0.031，核苷酸多样性为0.0017 ± 0.0001；P. fluviatilis整体单倍型多样性为0.892 ± 0.012，核苷酸多样性为0.0032 ± 0.0001。P. fluviatilis的多样性和平均核苷酸差异数均高于R. rutilus。成对F_ST分析中，R. rutilus的F_ST范围为0–0.269，中位数为0.072，多数比较为低至中度分化；P. fluviatilis的F_ST范围为0–0.816，约40%的比较超过0.15，分化范围更广。层级AMOVA显示，两个物种的流域间变异组分均不显著（R. rutilus：F_CT = 0.026，p = 0.188；P. fluviatilis：F_CT = 0.099，p = 0.086），但种群内变异显著，且P. fluviatilis的总体分化（F_ST = 0.285）远强于R. rutilus（F_ST = 0.110）。这说明水文边界并不能解释两个物种的遗传分化，且P. fluviatilis的结构化程度更高。

**3.3 单倍型网络与系统发育关系**

中位连接网络显示，R. rutilus呈明显的星状结构，以单一优势单倍型H_6为中心，该单倍型包含98条序列并见于所有采样点及GenBank记录，其他单倍型多通过一步突变与之相连，缺乏地理结构。P. fluviatilis的网络则缺乏单一优势中心，存在多个中等频率单倍型，且单倍型在不同采样点间的共享程度较低，网络结构更异质。最大似然系统发育树在两个物种中均为浅层拓扑，枝长较短，未支持任何深支系或明确谱系，网络划分的簇（A–D）仅为描述性分组，不代表系统发育支持的单系群。

**3.4 种群历史动态**

中性检验中，R. rutilus的Tajima's D为负（?0.835，p = 0.273）但不显著，Fu's Fs显著为负（?3.291，p = 0.002）；P. fluviatilis的Tajima's D（?0.039）和Fu's Fs（?0.798）均为负但不显著。错配分布分析中，R. rutilus呈单峰分布，但最小二乘法拟合未收敛；P. fluviatilis呈双峰分布，对扩张模型的拟合未遭拒绝（SSD和粗糙指数均不显著）。贝叶斯天际线图显示，两个物种均在晚更新世（约5万至10万年前）出现有效种群大小增长，其中P. fluviatilis的增长幅度更明显，R. rutilus增长更平缓，但95% HPD区间都很宽。综合来看，R. rutilus的种群扩张信号更强，而P. fluviatilis的信号较弱且更不具一致性。

**讨论与结论**

讨论部分指出，尽管两个物种共享同一水文景观，但遗传结构差异显著：R. rutilus因活动能力强、栖息地利用广泛且存在季节性迁移，可能促进基因流，因此分化弱、单倍型共享广泛；P. fluviatilis位点忠诚度高、迁移距离短，更易形成局部分化。研究还排除了流域边界的决定性作用，并强调人为放流、历史冰期避难所、有效种群大小差异等可能影响遗传格局，但现有单亲遗传线粒体标记不能完全区分这些过程的相对贡献。建议未来结合核基因组标记、更大范围采样和模型化分析，以更准确界定管理单元。

研究结论可概括为：本研究以相同标记和分析框架比较了同一水文景观中两种生态差异明显的共存鱼类，成功分离了共享水系配置与物种特异性生态的影响。在R. rutilus和P. fluviatilis中，遗传分化均未与主要流域边界一致对应，说明当代水文划分本身不决定遗传变异。相反，两个物种在分化的程度和空间模式上均存在差异：P. fluviatilis结构更强、空间更异质，而R. rutilus结构弱、单倍型广泛共享，这与基因流或近期扩张相符。这种种间对比表明，扩散能力和生态差异而非水系配置是这些系统中种群结构的主要驱动因素。该种间差异应视为一种总体趋势而非绝对界线，并且推断受限于单一母系遗传线粒体位点和有限的地理采样。未来需整合核基因组数据、扩大地理覆盖并采用基于模型的种群历史分析，以确证这些模式并解析塑造波罗的海淡水鱼类遗传多样性的过程。