基于基因组学解析濒危塔斯马尼亚淡水鱼天鹅南乳鱼（Galaxias fontanus）自然与迁地保护群体的遗传多样性与保护管理策略

《Aquatic Conservation: Marine and Freshwater Ecosystems》：Conservation Genomics of Natural and Translocated Populations of the Swan galaxias (Galaxias fontanus), an Endangered Tasmanian Freshwater Fish

　　

摘要：对濒危物种残存于破碎化种群中的基因组分析，可揭示遗传多样性的分布格局，厘清历史连通性，并为各群体的适应潜力提供见解。此类分析还可直接为保护性迁地等管理措施提供科学依据。天鹅南乳鱼（Galaxias fontanus）是一种溪流栖息鱼类，局限于澳大利亚塔斯马尼亚州两个河流流域的少数源头支流中，是包括迁地保护在内的保护管理重点对象。研究人员采用双酶切限制性位点关联DNA测序（ddRAD sequencing）技术，基于利用PacBio HiFi长读长测序生成的612-Mb重叠群水平基因组组装（contig-level genome assembly），对7个核心天鹅南乳鱼群体（3个自然群体和4个迁地群体）的遗传多样性进行了表征，并测定了单核苷酸多态性（SNP）变异。所有群体均表现出显著的遗传分化，表明历史连通性有限，其中流域间的遗传分化最为强烈。麦夸里河（Macquarie River）流域两个自然群体的核苷酸多样性（nucleotide diversity, π）是斯旺河（Swan River）流域仅存自然群体的两倍。迁地群体保持了与其源群体相当的遗传多样性，但等位基因频率发生了明显偏移，这与迁地后30年间两个群体均经历遗传漂变（genetic drift）的预期一致。杜克斯河（Dukes River）迁地群体保留了上斯旺河主河道中已灭绝鱼群的遗传变异。在赛格尼特河（Cygnet River）这一两物种发生接触的迁地地点，未发现与短鳍南乳鱼（G. brevipinnis）杂交的证据。这些基因组数据直接指导了两次新的多来源混合迁地保护的设计与实施，并为天鹅南乳鱼的持续保护管理奠定了坚实基础。

论文解读：濒危鱼类天鹅南乳鱼（Galaxias fontanus）的保护基因组学研究

一、研究背景与科学问题

全球生物多样性危机日益严峻，淡水鱼类及其栖息地正承受着不成比例的灭绝压力。多数淡水鱼类被不可栖息环境分隔为离散种群，类似于岛屿陆生物种，因而表现出高度的种内遗传结构。维持遗传多样性长期被视为保护生物学的核心目标，基因组技术的进步使得生成可直接指导管理决策的精细数据集成为可能。然而，遗传多样性指标在预测保护成效方面的价值近期受到质疑，这促使研究者重新审视遗传变异在种群存续中的作用。小种群规模导致遗传变异丧失，进而通过近交衰退、遗传漂变负荷和适应潜力降低而削弱种群存续能力，这一基本原理已得到再次确认。

天鹅南乳鱼是一种小型（体长小于135毫米）濒危淡水限制性鱼类，1978年发现于塔斯马尼亚东部上斯旺河，是塔斯马尼亚特有种。其自然分布跨越两个独立排水系统：东海岸的斯旺河流域和经南埃斯克河汇入北海岸的麦夸里河流域。该物种被列入澳大利亚政府2022–2032年优先受威胁物种行动计划，凸显其危急状态。几乎所有已知种群均位于某种自然屏障之上，使其栖息地不受入侵性褐鳟（Salmo trutta）影响。褐鳟是该物种的关键威胁，其幼鱼在开阔水域集群的习性使其极易成为捕食目标。1990年，斯旺河模式产地瀑布上方发现褐鳟入侵，数年内该地天鹅南乳鱼种群几乎灭绝，仅存个体局限于不适合褐鳟的小支流。目前，近10年内仅观察到3个自然种群存续：蓝阶溪（高数量）、乳品溪（中数量）和斯旺河支流（低数量）。为应对保护挑战，1989至1995年间实施了10次保护性迁地，其中4个迁地群体保持相对稳定：杜克斯河（高数量）、赛格尼特河（中数量）、罗卡河（低数量）和圣保罗河（低数量）。赛格尼特河迁地群体建立后，短鳍南乳鱼于2011年自然定殖该地，两物种形态相似导致监测中难以可靠区分，引发杂交可能性的推测。该物种恢复计划提出了两项保护遗传学优先事项：表征遗传结构和多样性以保护遗传上独特的管理单元，以及评估迁地群体是否保留自然群体的遗传多样性。本研究正是针对这些优先事项展开。

二、研究内容与主要结论

研究人员对7个天鹅南乳鱼高丰度群体（3个自然群体和4个迁地群体）进行了基因组分析。首先利用PacBio HiFi长读长测序生成了该物种的参考基因组，随后采用ddRAD测序量化群体水平遗传变异。对自然群体，研究评估了遗传独特性和多样性，并比较了流域内和流域间的差异以推断历史连通性；对迁地群体，比较了其与源群体的遗传多样性；最后检验了赛格尼特河迁地地点天鹅南乳鱼与短鳍南乳鱼间是否存在杂交。研究结果显示：所有群体均遗传独特，流域间分化最强（平均F_ST = 0.338），麦夸里流域两个自然群体间也存在强分化（F_ST = 0.169）；麦夸里流域自然群体核苷酸多样性（蓝阶溪π = 0.001089，乳品溪π = 0.000813）约为斯旺流域（斯旺河支流π = 0.000364）的两倍；迁地群体大体保留了源群体的遗传多样性水平，但等位基因频率发生偏移，符合30年间遗传漂变预期；杜克斯河群体遗传多样性显著高于斯旺河支流，推测源自现已灭绝的上斯旺河主河道种群；下圣保罗河鱼群与杜克斯河群体遗传上高度相似（F_ST = 0.019），含杜克斯河特征标记；赛格尼特河群体未检测到短鳍南乳鱼等位基因渗入。这些基因组数据直接指导了2023年实施的两次新迁地保护：德林溪迁地汇集了蓝阶溪和杜克斯河来源的鱼类，涵盖两个流域最具多样性的群体；埃弗克里奇河迁地则整合了杜克斯河/下圣保罗河与赛格尼特河来源的鱼类。

三、关键技术与方法

研究采用的主要技术方法包括：PacBio HiFi长读长测序生成612-Mb重叠群水平基因组组装，BUSCO评估显示97.47%完整基因；整合ab initio预测、蛋白同源映射和转录组组装的多策略基因注释流程；ddRAD测序对45个样本（5个2013年采集冷冻样本，40个2020–2021年乙醇保存样本）进行基因分型，比对至参考基因组后经Stacks流程鉴定SNP；经严格过滤获得57,429个SNP的全数据集和8,517个SNP的均匀间隔数据集；采用邻接法构建个体树、主成分分析（PCA）评估群体结构；利用JeDi流程估计核苷酸多样性和Hudson's F_ST；通过物种诊断性SNP检验杂交。

四、研究结果

基因组组装与注释：利用98.10-Gb PacBio HiFi读长生成611.56-Mb基因组组装，含2,270个重叠群，N50为875.5 kb，BUSCO评估显示97.47%基因完整。注释获得157.87 Mb重复序列（占25.81%）和20,531个蛋白编码基因。

遗传多样性与群体结构：初始数据集含176,136个SNP，过滤后全数据集含57,429个SNP，均匀间隔数据集含8,517个SNP。自然群体中，麦夸里流域蓝阶溪和乳品溪个体杂合度显著高于斯旺流域斯旺河支流。群体核苷酸多样性在麦夸里流域约为斯旺流域两倍。流域间遗传分化极强，麦夸里流域内两自然群体分化也显著。罗卡河迁地群体与源群体蓝阶溪遗传多样性相似（平均uHe分别为0.000892和0.000990），分化较小（F_ST = 0.072）。赛格尼特河群体遗传多样性最低（平均uHe = 0.000329），与其假定源群体斯旺河支流（平均uHe = 0.000342）相似，但群体间分化中等（F_ST = 0.117）。杜克斯河群体遗传多样性较高（平均uHe = 0.000613），下圣保罗河群体与之匹配（平均uHe = 0.000582）并含杜克斯河特征标记。树分析和PCA均先按流域、再按自然群体、最后按迁地群体与源群体关系分离遗传变异。均匀间隔数据集中1.17%的SNP在流域间固定。

杂交评估：结合34个塔斯马尼亚短鳍南乳鱼样本和34个天鹅南乳鱼样本，鉴定出84,469个物种诊断性SNP（占40.2%）。赛格尼特河天鹅南乳鱼群体中99.92%的等位基因与其它地点天鹅南乳鱼固定等位基因匹配，未发现短鳍南乳鱼等位基因渗入。

五、讨论与结论

自然群体遗传多样性与结构方面，麦夸里流域群体遗传多样性较高，表明该水系历史上曾存在更大的种群；斯旺流域仅存一个隔离区域，历史上可能曾有更高多样性的种群。流域间最大遗传距离最可能源于历史跨分水岭事件，而非独立陆封事件。麦夸里流域内两自然群体间约150公里低梯度河道无明显自然屏障，但显著遗传分化表明历史基因流有限，过去200年间褐鳟入侵和河流改造进一步限制了近期移动。斯旺流域内杜克斯河群体与现斯旺河支流群体的遗传差异表明历史上存在隔离，可能源于陡峭梯度或建群瓶颈。

迁地群体遗传学方面，两对源-迁地群体（蓝阶溪至罗卡河，斯旺河支流至赛格尼特河）杂合度相似，表明迁地过程中未发生极端瓶颈。但群体间存在中等到强分化（F_ST = 0.072和0.117），符合30年间遗传漂变驱动的等位基因频率分化。杜克斯河与上圣保罗河两次迁地事件中，杜克斯河群体遗传多样性较高且丰度高，而上圣保罗河群体丰度低，可能反映栖息地质量差异或遗传多样性差异导致的适合度差异。下圣保罗河鱼群与杜克斯河群体遗传相似，未显示来自上圣保罗群体的下游扩散遗传信号。

杂交方面，赛格尼特河样本中未检测到杂交个体，但样本量小（n = 11），未采样杂交个体可能存在。两物种在该地共存提供了理解其相互作用和共存潜力的重要机会。

保护管理方面，研究识别出四个遗传集群。从遗传角度，乳品溪和杜克斯河群体是最高优先级的新迁地来源，因其遗传谱系各仅存在于单一地点。蓝阶溪（或罗卡河）来源的新迁地也有益于维持这一目前仅限两个地点的多样性谱系。斯旺河支流和赛格尼特河群体因共享且有限的遗传多样性而优先级较低。混合来源迁地虽可能破坏共适应基因复合体，但在迁入新栖息地时能最大化适应潜力。2023年实施的两次新迁地已纳入遗传考量：德林溪迁地汇集蓝阶溪和杜克斯河来源鱼类，埃弗克里奇河迁地整合杜克斯河/下圣保罗河与赛格尼特河来源鱼类。未来迁地应优先纳入乳品溪及两个未纳入当前遗传分析的自然群体（格林斯溪和塔特花园溪）的鱼类。通过揭示遗传多样性在自然和迁地群体间的分布格局，本研究为天鹅南乳鱼的持续恢复提供了坚实的基因组基础，持续将基因组见解融入保护管理将有助于保护该物种的进化潜力。