🧬 biology

A highly contiguous reference genome for the Argentine hake (Merluccius hubbsi Marinni 1923, Gadiformes, Merluccidae)

本研究首次为具有重要经济与生态价值的阿根廷无须鳕（*Merluccius hubbsi*）提供了高质量的从头参考基因组，这不仅阐明了其进化历史与种群数量波动，还证实了尽管经历了历史性的气候变化和高强度的捕捞压力，该物种目前仍保持着强大的遗传多样性且不存在近期的近交现象。

原作者： Melisa E Magallanes Alba, Agustin Baricalla, Alejandro D’Anatro

发布于 2026-07-29

📖 1 分钟阅读☕ 轻松阅读

CC BY 4.0

原作者： Melisa E Magallanes Alba, Agustin Baricalla, Alejandro D’Anatro

原始论文采用 CC BY 4.0 许可（https://creativecommons.org/licenses/by/4.0/）。 ⚕️ 这是一篇未经同行评审的预印本的AI生成解释。这不是医疗建议。请勿根据此内容做出健康决定。阅读完整免责声明

想象一下，人体就像一座巨大的图书馆，包含了构建和运行一个人所需的所有指令。在生物学世界中，这座图书馆被称为基因组，其中的书是用由 A、C、G 和 T 四个字母组成的密码编写的。长期以来，科学家只能读取许多动物基因组中那些被撕碎的微小书页，这使得很难理解完整的故事。但最近，技术取得了进步，使我们能够将这些碎片重新缝合在一起，组成完整的、高分辨率的书籍。这就是我们所说的参考基因组。

为什么这很重要？把一个物种想象成一款复杂的电子游戏。如果你想修复一个漏洞（bug）、了解一个角色的移动方式，或者预测他们对新环境的反应，你需要完整的源代码。没有它，你只能靠猜测。对于像阿根廷鳕鱼（Argentine hake）这样人类食用的鱼类来说，拥有这份完整的“源代码”至关重要。它能帮助科学家研究鱼群是否健康、它们是如何在过去的冰河时代幸存下来的，以及如何保持捕捞的可持续性，从而避免“游戏崩溃”。这篇论文讲述了如何终于为这种非常重要的鱼类写下那份完整的源代码。

阿根廷鳕鱼的大书

认识一下阿根进制鳕鱼（Merluccius hubbsi）。它是西南大西洋海洋中的超级明星，这种鱼价值极高，以至于乌拉圭和阿根廷等国家都依赖它的渔业。但长期以来，科学家一直试图在黑暗的房间里用手电筒来阅读这条鱼的遗传故事。他们只有片段信息，却没有完整的图景。当一个研究小组决定为这个物种构建第一个高质量、“连续的”（意味着所有碎片都整齐地粘合在一起）参考基因组时，情况发生了改变。

把基因组组装想象成解决一个拥有数十亿块拼图的巨大拼图。研究人员并没有只使用一种类型的拼图块；他们使用了一种混合策略。他们将 PacBio 长读段（类似于可以跨越复杂、重复部分的细长丝带）与 NovaSeq 短读段（类似于修复任何拼写错误的超精确印章）结合起来。通过混合这两种技术，他们成功组装出了一个跨度约为 425.3 Mb（兆碱基对）的基因组。换句话说，如果这个基因组是一本书，它将是一本厚重的卷册，包含 28,998 个蛋白质编码基因——即构建鱼类身体部位和功能的实际指令。

遗传健康检查

一旦有了这本书，科学家就开始阅读它以检查鱼类的健康状况。他们担心多年的过度捕捞可能会伤害鱼群，导致近亲繁殖（就像小村庄里的表亲结婚一样），这通常会导致遗传健康受损。

他们寻找一种被称为**纯合性连续段（Runs of Homozygosity, ROH）**的现象。想象你的 DNA 是一串彩色的珠子。如果你在两条链上都有长段完全相同颜色的珠子，那就是近亲繁殖的迹象。研究人员发现了一些令人惊叹的事实：阿根廷鳕鱼几乎没有长段的相同颜色珠子。其近交系数小于 0.01%。这表明，尽管面临着巨大的捕捞压力，阿根廷鳕鱼的种群仍然非常庞大，就像一锅巨大的、健康的汤一样混合在一起，并没有遭受近期的遗传崩溃。这是种群具有强健性、基因自由混合的泛群现象（panmictic population）的迹象。

冰河时代的时光机

该团队还利用基因组作为时光机，观察鱼类种群在数千年间是如何变化的。通过使用一种称为 PSMC 的方法，他们追溯了“有效种群规模”（一种计算实际进行繁殖的鱼类数量的方法）的历史。

他们发现的故事非常戏剧化。大约在 10 万年前，种群达到了一个巨大的高峰，约有 12.5 万个有效个体。这很可能发生在气候温暖的时期，当时的海底栖息地非常完美。但随后，随着世界变冷并进入末次盛冰期（最后一次冰河时期的巅峰），种进行规模开始萎缩，到 2 万年前下降到了约 8 万个个体。这表明这些鱼在海洋冻结的过程中幸存了下来，但寒冷确实压缩了它们的数量。有趣的是，这种模式与它们的欧洲近亲不同，后者找到了安全的温暖避难所。阿根廷鳕鱼必须在巴塔哥尼亚陆架上熬过这场风暴。

家族树与大迁徙

研究人员还组装了鱼类的线粒体基因组（一种仅由母亲传递的微小、独立的环状 DNA）来构建整个鳕科（Merluccius）的家族树。他们发现，该属可能在大约 3400 万年前（晚始新世）起源于北大西洋。

家族树中最令人兴奋的部分是阿根廷鳕鱼所处的位置。研究证实，阿根廷鳕鱼是太平洋鳕鱼的“姐妹分类群”（即最亲近的亲戚）。这意味着，阿根廷鳕鱼的祖先是连接该家族在中新世时期（约 1500 万年前）从大西洋跨越到太平洋的桥梁。它们是其进化旅程中至关重要的踏脚石。

未来展望

虽然这个新基因组是一个巨大的飞跃，但作者谨慎地指出，这仅仅是个开始。他们使用的是单条鱼的 DNA，因此虽然它告诉了我们深层的历史和普遍的健康状况，但目前还无法绘制出今天不同特定区域鱼类之间的微小差异。此外，由于使用了混合测序方法，某些极具重复性的基因组部分可能与更新的超长读段技术相比被略微低估了。

然而，这份“源代码”现在已向所有人开放。它是一个基础工具，将帮助科学家更好地管理渔业，理解鳕鱼如何适应气候变化，并确保这一宝贵的资源能够世代留在我们的餐桌上。阿根廷鳕鱼终于在生命之书里发出了属于自己的声音。

技术摘要：阿根廷无刺尾鱼（Merluccius hubbsi）的高连续性参考基因组

问题陈述
阿根联无刺尾鱼（Merluccius hubbsi）是西南大西洋海洋生态系统中的核心底栖物种，也是乌拉圭和阿根廷的主要商业渔业资源。尽管其具有重要的生态和经济价值，但该物种一直承受着剧烈且持续的捕捞压力。由于缺乏高分辨率的基因组资源，有效的保护和种群评估一直受到阻碍。以往的研究依赖于低密度的分子标记或高度破碎的基因组组装，导致在理解该物种的基因组架构、进化历史及当前遗传健康状况方面存在空白。迫切需要从基础的遗传监测转向全面的基因组管理，以应对种群动态和适应能力问题。

方法论
本研究采用混合测序策略，生成了首个高质量的 M. hubbsi de novo 参考基因组。

测序： 使用 PacBio Sequel 长读段（目标插入片段为 10–20 kb）和 Illumina NovaSeq 短读段（151 bp 双端测序）对单个野生个体的超高分子量 DNA 进行测序，实现了约 112X 的覆盖度。
组装： 使用 MaSuRCA 构建了混合 de novo 草图组装。利用 Redundans 处理草图以减少单倍型冗余。使用 RagTag 进行支架搭建（scaffolding），利用亲缘关系较近的欧洲无刺尾鱼（M. merlucius）的染色体水平组装作为参考进行排序和定向。最后使用 Illumina 读段通过 POLCA 进行组装校正，以纠正碱基水平的错误。
注释： 结构注释采用了结合 BRAKER3（利用来自鳕鱼 Gadus morhua 和 OrthoDB 的蛋白质同源性）与 LiftOn（从 G. morhua 转移注释）的混合共识方法。功能注释则使用 eggNOG-mapper 完成。
群体与人口统计学分析：
- 杂合度与近交： 将 Illumina 读段回比对至参考基因组以进行 SNP 调用。通过识别纯合片段运行（Runs of Homozygosity, ROH）来计算基因组近交系数（ $F_{ROH}$ ）。
- 历史人口统计学： 应用成对序列马尔可夫合并（PSMC）模型，对二倍体共识序列进行分析，以重建历史有效种群大小（ $N_e$ ）的波动，并根据 3.5 年的世代时间和 $1.0 \times 10^{-8}$ 的突变率进行缩放。
系统发育学： 使用 GetOrganelle 组装了完整的线粒体基因组，并通过 MitoFish 进行注释。利用 BEEST v2.7.8，结合 13 个串联蛋白编码基因和三个鳕形目（Gadiformes）外类群，利用化石记录和生物地理事件进行时间定标的贝叶斯系统发育推断。通过离散系统地理学方法进行了祖先区域重建。

主要贡献与结果

基因组组装： 最终组装规模约为 425.3 Mb，包含 8,622 个支架（scaffolds）。其 N50 为 19.19 Mb，是迄今为止连续性第二高的 Merluccius 组装结果，仅次于 M. merlucius 的染色体水平组装。根据 BUSCUS 分析（actinopterygii_odb12），组装的完整度为 63%。
重复序列景观： 基因组中约 31.51% 由重复序列组成，其中以间插重复序列（24.53%）为主，特别是 LINEs 和 DNA 转座子。很大一部分（18.82%）的重复序列未被分类，这表明存在一个在公共数据库中未发现的谱系特异性转座元件库。
基因注释： 研究鉴定出 28,998 个蛋白编码基因。在过滤冗余后，保留了 27,276 个非冗余蛋白序列，其中 79.1% 成功分配了功能注释（GO term、KEGG 通路）。
基因组健康状况： 对纯合片段运行（ROH）的分析显示，几乎不存在近期近交现象（ $F_{ROH} < 0.01\%$ ），且总 ROH 长度仅为 3.7 kb。这表明尽管经历了数十年的高强度捕捞压力，该种群仍维持着大规模、泛性的基因库以及高水平的基因流。
历史人口统计学： PSMC 模型表明，在 1,000,000 至 300,000 年前，有效种群大小（ $N_e$ ）稳定在 50,000 至 60,000 之间。随后发生了种群扩张，在约 100,000 年前达到约 125,000 个体的峰值（可能对应于间冰期条件），随后在导致末次盛冰期（LGM）的降温期出现下降。
系统发育与生物地理学： 时间定标的系统发育分析将 Merluccius 属的起源定位于晚渐新世（约 34 MYA）。旧世界（欧非）与新世界（美洲）支系的分化发生在中新世中期（约 15 MYA）。分析确认 M. hubbsi（西南大西洋）是太平洋无刺尾鱼辐射演化的姐妹分类群，支持了“向西南大西洋的迁徙是入侵太平洋前的过渡步骤”这一假设。

意义
作者认为，这一高质量的参考基因组为未来关于种群评估、功能适应以及西南大西洋海洋渔业可持续管理的研究提供了宝贵的资源。研究结果表明，M. hubbsi 保持了强大的遗传多样性，并在经历剧烈的更新世气候波动时并未遭受近期的基因组侵蚀。然而，作者指出，虽然单体基因组分析提供了深刻的历史洞察，但无法捕捉当代精细的空间种群结构。因此，本文强调该资源是未来进行多个体、群体水平基因组筛选的基础框架，而这些筛选对于阐明面对气候变化和捕捞压力下的现代种群动态及局部适应至关重要。此外，该研究也强调了为非模式海洋硬骨鱼类生成 de novo 重复序列库的必要性，因为目前的公共数据库缺乏对其谱系特异性转座元件的充分代表性。