《Diversity and Distributions》：Genomics Reveals Anthropogenic Fragmentation and Introgression in a Highly Endangered New Zealand Freshwater Fish
编辑推荐：
研究目标：新西兰特有的淡水鱼类区系受到排水系统改造和鲑科鱼类（salmonid）捕食者引入的严重影响。本研究旨在检验人为种群破碎化对一种极度濒危特有鱼类的遗传多样性的影响，为该物种未来的保护管理提供信息。
研究地点：新西兰南岛 Teviot 河及周边水系。
研
研究目标：新西兰特有的淡水鱼类区系受到排水系统改造和鲑科鱼类（salmonid）捕食者引入的严重影响。本研究旨在检验人为种群破碎化对一种极度濒危特有鱼类的遗传多样性的影响，为该物种未来的保护管理提供信息。
研究地点：新西兰南岛 Teviot 河及周边水系。
研究方法：研究人员生成并分析了包含275个 Galaxias 标本、共47,754个 SNP（单核苷酸多态性）的全基因组数据集，以评估极度濒危淡水鱼类 G. 'Teviot' 全分布范围内的遗传多样性。研究人员进行了系统发育、STRUCTURE（群体遗传结构分析）和遗传多样性分析，以评估 Teviot 河种群在该区域重大人为河流排水改造背景下的进化关系。
研究结果：分析揭示了与人为种群破碎化相关的遗传分化，其中 G. 'Teviot' 内部最深的谱系地理分裂与第一座 Lake Onslow 大坝（建于1890年）相关。研究还表明，随后的排水改造（水渠连接）与一个新识别候选分类单元 G. 'Manor Burn' 的等位基因明显渐渗进入 Teviot 河上游种群有关。与直觉相反的是，Lake Onslow 大坝似乎对 Teviot 河下游种群免受这种人为基因渐渗起到了缓冲作用。
主要结论：研究人员确定了 G. 'Teviot' 的关键遗传多样性丰富且未发生渐渗的种群，可作为潜在易地保护（translocation）的来源，以从遗传层面保障这一极度濒危鱼类物种。
论文解读文章
一、研究背景与科学问题
淡水鱼类辐射演化常常在极小的空间尺度上表现出高度的生物多样性，这一现象通常归因于淡水生态系统的地质与生态隔离。然而，随着人类对淡水系统的压力不断增加，监测和缓解人为变化对生物多样性的影响已成为迫切需求。全球范围内，淡水生物多样性尤其受到种群灭绝和破碎化的威胁，其驱动因素包括栖息地退化、人为扩散屏障以及外来类群的引入。
新西兰特有的淡水鱼类区系以南半球广布属 Galaxias（南乳鱼属，胡瓜鱼目）为主。其中，淡水局限型的 Galaxias vulgaris 复合群包含12个目前公认的淡水局限型分类单元，特有种分布于新西兰南岛各河流。该复合群中多个窄域特有分类单元面临重大保护关注。特别是濒危的 Teviot 平头南乳鱼（G. 'Teviot'）主要局限于 Clutha 河水系内 Teviot 河的少数小型支流中。该分类单元的狭窄地理分布范围因引入鲑科捕食者、淡水蓄水工程及其他排水改造而进一步破碎化，被列为"国家极危"等级。
基因组方法在淡水生物多样性的保护与管理中日益重要，对于量化淡水种群遗传多样性、连通性以及应对人为基因渐渗和生物多样性下降至关重要。本研究利用基因分型测序（GBS，Genotyping-by-Sequencing）技术，对 G. 'Teviot' 全分布范围进行全基因组分析，旨在表征该物种的种群遗传多样性，以识别具有潜在恢复力和多样性的种群，为保护优先级排序和资源分配提供依据。
二、研究内容与主要结论
新西兰环保部在 Teviot 河水系及邻近的 Manor Burn 排水区共12个采样点采集了95个 Galaxias 标本，每个位点通常采集10个个体，少数种群较小的位点样本量较少。研究人员将新样本与 Campbell 等人（2022）的多物种 Galaxias 数据集整合，以 Galaxias brevipinnis（短鳍南乳鱼）参考基因组为比对参考，最终获得包含275个标本、47,754个 SNP 的全基因组数据集。
系统发育分析表明，Teviot 河和 Manor Burn 样本分为三个不同的系统发育簇：Teviot 河下游 Old Hut Creek 样本与 G. pullus 紧密聚类；其余大部分 Teviot 样本与先前测序的 G. 'Teviot' 样本聚类；而 Manor Burn 系统样本及 Teviot 位点 T4 的一个样本与先前测序的该 Clutha 亚流域及邻近水域样本聚类，构成"Manor Burn"支（G. 'Manor Burn'），该支是 G. 'Teviot' 的姐妹群。
STRUCTURE 分析在 K=2 和 K=3 时均识别出强烈的遗传分组，三个主要基因型簇与系统发育分析检测到的三个支系相对应。分析揭示了 G. 'Manor Burn' 等位基因向 G. 'Teviot' 的核基因组渐渗信号，但该渐渗仅限于 Teviot 上游区域（North Branch 样本 T5、T6 及 Lake Onslow 的一个小支流 T4），且总体上程度有限。其中一个 Lake Onslow 支流样本（T4）表现出强烈的杂交特征（约75% G. 'Manor Burn' 祖先成分），可能代表"纯" G. 'Manor Burn' 与 F1 杂交个体的回交后代。相比之下，G. pullus 和 Teviot 系统下游的 G. 'Teviot' 种群几乎没有渐渗证据。
系统发育分析揭示了 G. 'Teviot' 内部浅层但显著的谱系地理结构，样本通常按采样点聚类。单个支流内相邻样本（相距约1-2公里）存在精细尺度结构。更重要的是，被人为屏障隔离的 G. 'Teviot' 样本间检测到显著更深的结构，最深的分裂对应于1890年 Lake Onslow 大坝的隔离效应。该人为屏障也界定了 G. 'Manor Burn' 渐渗进入 G. 'Teviot' 的推断地理界限。在上游区域，一个相对较浅的谱系地理分裂可能与1982年 Lake Onslow 的扩建有关，该扩建终止了 Lake Onslow 支流 T4 与 Teviot North Branch 之间的河流连通。然而，整个 Teviot 系统的 Mantel 检验未检测到显著的隔离距离效应（r2=0.01；p=0.54）。
全基因组杂合度估计在不同 G. 'Teviot' 种群间存在差异，范围约为0.10至0.15。部分情况下，同一溪流内"上游"（源头）与"下游"样本间检测到局部多样性差异，T2U 和 T3U 的杂合度分别低于 T2L 和 T3L。小型源头种群 T6 携带的多样性显著低于 T5。值得注意的是，发生渐渗的种群 T4 和 T5 属于遗传多样性最高的种群之列，未发生渐渗的 T2L 同样如此。
三、主要技术方法
研究人员采用以下关键技术方法开展研究：1）基因分型测序（GBS）技术，按照 Elshire 等人（2011）的 PstI 方案制备文库，在 Illumina NextSeq 2000 平台上进行双端150 bp测序；2）变异 calling 流程，使用 BCFtools 进行 SNP 检测，并以 Galaxias brevipinnis 参考基因组（GCA_057930015.1）进行比对，最终保留至少在80%个体中存在的双等位 SNP，最小等位基因频率（MAF）>0.05；3）系统发育分析采用 IQ-TREE v2.1.1 的最大似然框架，以 ModelFinder 选择最佳模型，5000次超快自举评估节点支持率；4）STRUCTURE v2.3.4 进行群体遗传结构分析，采用混合模型和相关性等位基因频率，K=1至K=5各运行5次独立重复，以 ΔK 方法确定最佳聚类数；5）利用 R 包 adegenet 计算各采样点平均观测杂合度，ARLEQUIN 3.5 计算种群间两两 FST，GENEPOP 4.7.5 进行 Mantel 检验评估地理距离与遗传距离的关联。样本来源为新西兰环保部采集的 Teviot 河水系及邻近 Manor Burn 排水区共12个采样点的95个标本。
四、研究结果
4.1 系统发育关系
基于47,754个全基因组 SNP 的系统发育分析获得了分辨率良好的最大似然树，大多数主要节点获得100%自举支持。Teviot 河和 Manor Burn 样本分为三个不同的系统发育簇：Teviot 河下游 Old Hut Creek 样本与 G. pullus 聚类；其余大部分 Teviot 样本与 G. 'Teviot' 聚类；Manor Burn 系统样本及 Teviot 位点 T4 的一个样本构成"Manor Burn"支（G. 'Manor Burn'），该支是 G. 'Teviot' 的姐妹群。
4.2 结构、杂交与渐渗
STRUCTURE 重复运行在各 K 值内结果一致，K=2 和 K=3 时均识别出强烈遗传分组，三个主要基因型簇与系统发育分析检测到的三个支系对应。分析揭示了 G. 'Manor Burn' 等位基因向 G. 'Teviot' 的核基因组渐渗，但仅限于 Teviot 上游区域（T4-T6），总体上程度有限。T4 位点一个样本表现出约75% G. 'Manor Burn' 祖先成分的强烈杂交特征，可能代表回交个体。Teviot 下游种群和 G. pullus 几乎没有渐渗证据。
4.3 种群结构与人为屏障
G. 'Teviot' 种群系统发育分析揭示了浅层但显著的谱系地理结构，样本通常按采样点聚类。被人为屏障隔离的样本间检测到显著更深的结构，最深的分裂对应于1890年 Lake Onslow 大坝的隔离效应，该屏障也界定了 G. 'Manor Burn' 渐渗的地理界限。1982年 Lake Onslow 扩建可能导致了 T4 与 T5-T6 之间的浅层分裂。整个 Teviot 系统未检测到显著的隔离距离效应。
4.4 遗传多样性
全基因组杂合度估计在不同 G. 'Teviot' 种群间存在差异，范围约0.10至0.15。部分溪流内"上游"样本（T2U、T3U）杂合度低于"下游"样本（T2L、T3L），小型源头种群 T6 多样性显著低于 T5。发生渐渗的种群 T4 和 T5 属于遗传多样性最高的种群之列，未发生渐渗的 T2L 同样具有较高多样性。
五、讨论与结论
人为种群破碎化：研究表明 G. 'Teviot' 现存少数种群间的遗传结构与 Teviot 河水系内的人为屏障密切相关。1890年 Lake Onslow 大坝对应于 G. 'Teviot' 内部显著的谱系地理分裂，是该系统遗传破碎化的关键驱动因素。系统中缺乏遗传隔离距离效应可能反映了这些小型且近期破碎化种群中遗传漂变的主导作用。此外，引入的鲑科鱼类通过竞争和捕食进一步加剧了种群破碎化。研究识别出未发生渐渗且遗传多样性相对丰富的种群（T2L）可作为潜在易地保护来源，该种群也是目前数量最丰富的剩余种群。随着破碎化加剧，源头种群间的人工易地保护可能成为最大化遗传多样性的重要工具。
人为易位与遗传侵蚀：人为排水改造可将物种跨自然流域转移，使先前隔离的谱系发生二次接触，导致杂交和特有分类群的遗传淹没。19世纪金矿开采常涉及建造连接相邻水系的人工水渠，导致无意中的本地鱼类易位。遗传分析揭示了 Manor Burn-Teviot 分水岭历史水渠连接的明显影响，G. 'Manor Burn' 基因渐渗进入 G. 'Teviot' 的三个种群（T4-T6）。与直觉相反的是，1890年 Lake Onslow 大坝可能在限制 G. 'Manor Burn' 等位基因向上游 Teviot 区域传播中发挥了意外作用——下游种群（T1-T3）未受渐渗影响，表明该渐渗被大坝有效"阻断"。据此推断，G. 'Manor Burn' 基因向 Teviot 系统的传播发生在1890年大坝建成之后、1982年湖泊进一步抬升之前。
隐秘多样性：G. 'Manor Burn' 谱系具有与已正式描述 Galaxias 分类单元相当的遗传分化，但缺乏独特的线粒体谱系，可能表明完全的线粒体"捕获"。这种核质不一致在该 Galaxias 辐射演化中广泛存在，分类多样性常被线粒体渐渗所掩盖。T4 位点发现一个近乎纯 G. 'Manor Burn' 基因型样本，暗示这两个谱系在二次接触中可能存在一定程度的生殖隔离，但需进一步采样验证。
结论：研究揭示人为屏障与极度濒危 G. 'Teviot' 现存少数种群的遗传破碎化密切相关。同时，这些物理屏障反而限制了人为排水改造后渐渗等位基因向该物种的传播。研究进一步识别出未发生渐渗且遗传多样性丰富的种群，可作为潜在易地保护来源，以协助这一高度濒危物种的保护管理。
