水文连通性对比人为破碎化：拉脱维亚与立陶宛河流域中欧亚河鲈（Perca fluviatilis）的种群遗传结构

《Diversity》：Hydrological Connectivity Versus Anthropogenic Fragmentation: Population Genetic Structure of Eurasian Perch Across Latvian and Lithuanian River Basins

　　

淡水生态系统日益受到生境改变与连通性下降的影响，这会影响水生物种的遗传结构。研究人员利用9个微卫星位点（microsatellite loci）调查了拉脱维亚道加瓦河（Daugava River）以及立陶宛尼曼河（Nemunas River）与内里斯河（Neris River）等主要水系中 Eurasian perch（Perca fluviatilis L.）的种群遗传结构。研究人员通过标准种群遗传学分析评估遗传多样性数量参数，并基于贝叶斯聚类（Bayesian clustering）与系统发育重建评价种群遗传结构。所有受检种群均表现出较高遗传多样性，每位点等位基因数（NA）为4至24，观测杂合度（Ho）为0.536至0.716。种群间遗传分化低至中等但显著，道加瓦河流域内 FST 为0.021至0.079，尼曼河与内里斯河间为0.034。贝叶斯聚类与邻接树（Neighbor-Joining）分析揭示区域和局部尺度的遗传结构，包括道加瓦河因三座水电站而被分割后出现的明显亚结构。上述模式表明种群连通性在研究区并不均一，可能反映历史与当代过程的共同影响。这些发现增进了对波罗的海河流系统中河鲈种群结构的理解，并为淡水鱼类遗传资源保护与管理提供相关信息。

研究背景方面，淡水生态系统对维持生物多样性与生态稳定具有基础作用，但受水电调节、城市化和水坝建设等人为压力影响突出。河流系统尤其容易发生生境转化与纵向连通性中断，进而限制扩散和基因流，导致遗传分化增强、遗传多样性下降和适应潜力减弱。波罗的海地区拥有破碎化程度不同的河流系统，其中拉脱维亚道加瓦河被普利亚维尼亚斯、凯古姆斯和里加水电站梯级改造，而立陶宛尼曼河与内里斯河水系的大型水电屏障较少。已有研究多关注沿海、过渡带或湖泊种群，利用线粒体DNA（mtDNA，matochondrial DNA）或基因组手段开展广泛系统地理分析，但针对拉脱维亚与立陶宛主要河流、使用核标记开展跨境比较的研究仍然有限。因此，研究人员以生态可塑性强、兼具定居与迁移行为、可作为连通性与破碎化指示种的 Eurasian perch（河鲈）为对象，比较不同破碎化水平下河流种群的遗传结构。研究人员假设不同流域种群存在显著遗传分化，道加瓦河内部分化强于立陶宛水系，且水电大坝影响河段分化更明显，同时历史流域结构与冰后期定殖历史也会塑造遗传变异。论文成果发表于《Diversity》。

关键技术方法方面，研究人员共采集立陶宛尼曼河23尾、内里斯河10尾，以及拉脱维亚道加瓦河斯克里韦里24尾、列普卡尔尼22尾、普利亚维尼亚斯20尾样本；其中道加瓦河样本分别位于普利亚维尼亚斯大坝上、下游不同河段。研究人员采用盐提法提取肌肉组织核DNA（nuclear DNA），使用9个微卫星位点（其中6个二核苷酸、2个三核苷酸、1个四核苷酸位点）进行基因分型。数据分析依靠 MICRO-CHECKER 检查零等位基因（null alleles）与分型错误，POPGENE 与 GenAlEx 计算等位基因数、观测/期望杂合度（Ho/He）、私有等位基因丰富度（NRPA）与基于 AMOVA 的 F_ST；HP-RARE 做等量样本稀化下的等位基因丰富度（N_RA）；BOTTLENECK 检测近期瓶颈；STRUCTURE 结合 LOCPRIOR 模型进行贝叶斯聚类，并用 Evanno 的 ΔK 判断最适簇数 K；Populations 软件基于 Nei 遗传距离构建邻接树；GenAlEx 实施 Mantel 检验与主坐标分析（PCoA，Principal Coordinates Analysis）。研究人员还将本研究5个种群与此前研究的里加河 Voleri、利耶卢佩河 Jelgava 与 Griva 样本合并为8个种群再做区域尺度分析。

研究结果方面，质量控制在基因型质控部分，MICRO-CHECKER 提示 YP60、PflaL6、PflaL10、PflaL2 与 YP111 等位点可能存在零等位基因并出现杂合子缺失；剔除 PflaL6 与 PflaL10 后，簇数与种群结构结论不变，说明主要遗传模式不是分型假象。遗传变异跨位点部分，9个微卫星位点多态性高，N_A 从4到24，He 为0.449至0.901，Ho 为0.386至0.902；总体每位点 F_ST 为0.028至0.143，多数位点间接估计的迁移个体数（Nm）较高，但 YP111 显示较低基因流。种群内遗传多样性部分，N_A 为60至77，N_RA 为5.55至5.87，NRPA 以普利亚维尼亚斯0.92与内里斯河0.85最高、尼曼河0.40最低；仅内里斯河等位频率分布出现偏离L形的模式偏移，提示可能近期瓶颈。Hardy–Weinberg 平衡（HWE）检验中，道加瓦河部分种群在 PflaL6、PflaL2、YP78、PflaL10 等位点显著偏离，其他种群基本符合平衡。种群间遗传分化部分，成对 F_ST 为0.021至0.079，全部显著；普利亚维尼亚斯与斯克里韦里间达0.079，尼曼河与内里斯河间为0.034。STRUCTURE 在 K=2 时将尼曼河、内里斯河与斯克里韦里归为一簇，列普卡尔尼与普利亚维尼亚斯归为另一簇；更高 K 值进一步拆分道加瓦河内部与立陶宛种群。邻接树支持尼曼河—内里斯河簇及列普卡尔尼—普利亚维尼亚斯高支持簇，斯克里韦里成独立分支。Mantel 检验显示遗传距离与地理距离几乎无相关，Rxy 约0.02，不存在隔离随距离现象。PCoA 前两轴解释67.38%变异，尼曼河、内里斯河与斯克里韦里接近，普利亚维尼亚斯与列普卡尔尼位置更孤立。整体遗传结构部分，研究表明流域间分化与道加瓦河内部亚结构共存，连通性在空间上不均一。更广区域尺度部分，8种群 STRUCTURE 在 K=2 将 Voleri、Jelgava、Griva 与其余种群分开；邻接树中列普卡尔尼—普利亚维尼亚斯亚簇支持率75%，Voleri—Griva 亚簇支持率87%；PCoA 前两轴解释57.07%，此前加入样本占据较孤立位置。

讨论部分总结，研究人员指出高遗传多样性与低至中等 F_ST 并存符合淡水鱼类中受限扩散与持续基因流共同作用的特征。道加瓦河内部差异化与长期无过鱼设施的水电梯级一致，但不能仅凭本数据量化大坝贡献，也可能来自冰后期定殖、河网拓扑和河鲈产卵地忠实性等自然过程。尼曼河与内里斯河水系微卫星分化较弱，提示当代连通性较高，但已有 mtDNA 研究显示该区存在历史谱系分组，说明母系遗传、雌体扩散偏向或历史连通性可能与核标记模式不同。内里斯河的模式偏移提示有效种群大小下降，但瓶颈信号短暂且方法敏感，需谨慎解释。管理含义上，不同流域甚至同流域不同河段可考虑作为独立管理单元；放流与移植若忽视结构可能造成遗传均质化与本地适应丧失；道加瓦河应优先监测河段连通性并评估屏障削减措施对河鲈基因流的效果。

结论部分翻译为：总之，本研究揭示了拉脱维亚与立陶宛河流系统中 Eurasian perch（Perca fluviatilis L.）的显著遗传分化与种群结构。虽然高水平遗传多样性表明种群内仍维持大量遗传变异，但所观察到的流域水平分化与河流内部亚结构说明研究区种群连通性并不均一。所检测到的遗传模式与历史及当代过程的影响一致，包括流域配置与水文连通性，而人为破碎化的潜在贡献仍需进一步研究。结果应结合若干局限来理解，包括不同河流系统采样年份不同、样本量不等、部分流域空间覆盖有限、部分位点可能存在零等位基因，以及无法评估群体历史、未记录移植或放流活动的影响。尽管存在这些局限，本研究仍为波罗的海河流系统的种群遗传结构提供了新认识，并凸显在淡水种群遗传研究中考虑多空间尺度的重要性。上述模式可为由其发展的未来保护与渔业管理策略提供参考，但需纳入更广空间采样、基因组标记与环境数据的后续研究，才能提出确定性管理建议。