参考引导的染色体规模基因组组装及大西洋巨型石斑鱼（Epinephelus itajara）种群基因组学研究——哥伦比亚罗萨里奥群岛

《Ecology and Evolution》：Reference-Guided Chromosome-Scale Genome Assembly With Insights on Population Genomics of the Atlantic Goliath Grouper (Epinephelus itajara), Islas del Rosario, Colombia

　　

大西洋巨型石斑鱼（Epinephelus itajara）是西北大西洋最大的石斑鱼物种，处于极度濒危状态。该物种具有重要的生态、文化和经济作用，并已成为哥伦比亚罗萨里奥群岛海洋馆（Oceanario of the Rosario Islands）圈养繁殖工作的焦点。然而，尽管其具有生态和保护重要性，E. itajara的基因组资源和种群基因组数据仍然稀缺，尤其是在哥伦比亚加勒比地区。本研究利用PacBio HiFi测序和Illumina技术，对该物种进行了参考引导的染色体规模基因组组装，并分析了其种群基因组结构。组装的基因组总大小为1.12 Gb，contig N50为42.69 Mb，scaffold N50为46.30 Mb。在重复元件约占46%的基因组进行掩蔽后，共鉴定出22,692个蛋白质编码基因。比较基因组分析揭示了其与近缘石斑鱼物种的高度共线性，并确定E. lanceolatus为其最近亲缘种，支持该属内近期分化及保守的基因组结构。此外，研究人员利用7,706个高质量SNP进行了种群基因组分析，以评估圈养种群的基因组结构。结果显示存在四个不同的基因组谱系，样本个体间具有中度遗传分化。在哥伦比亚加勒比地区，研究人员识别出两个独有谱系，分别与西斯帕塔湾（Bahía Cispatá）和巴尔沃科阿斯湾（Bahía Barbacoas）相关联，提示可能存在地理隔离。这些基因组资源为了解E. itajara的基因组多样性、进化历史和繁殖机制提供了宝贵工具和新机遇。此外，它们还为保护策略奠定了基础，包括旨在提高圈养种群基因组多样性的选择性育种计划，以及指导其自然栖息地的恢复工作。

**论文解读：大西洋巨型石斑鱼参考基因组及种群基因组学研究**

大西洋巨型石斑鱼（Epinephelus itajara）是西北大西洋最大的石斑鱼类，最大体长可达2米，体重约400公斤，具有生长缓慢、性成熟晚（雄性4–6年，雌性6–7年）、对出生地高度忠诚等特点。其幼鱼依赖红树林作为育幼场，成鱼则迁移至珊瑚礁或沉船等结构。该物种在生态链中处于顶级捕食者地位，并支撑着游钓、 subsistence 渔业及小型商业渔业，但过度捕捞导致其种群急剧下降。尽管国际自然保护联盟（IUCN）红色名录目前将其列为“易危”（vulnerable），其保护状况仍受到关注。此前，基于线粒体DNA和微卫星的研究显示，该物种在巴西、美国等分布区存在显著遗传分化和较低的遗传多样性，但这些传统标记仅能反映母系遗传或受限于突变模式和无效等位基因，无法在全基因组水平提供高分辨率信息。此外，哥伦比亚加勒比地区此前没有任何E. itajara基因组组装，当地种群的遗传组成几乎未知。因此，构建高质量参考基因组并开展种群基因组分析，对保护和管理该物种至关重要。

该研究由Jácome Castro等人完成，发表于《Ecology and Evolution》。研究人员利用PacBio HiFi测序和Illumina测序技术，对来自哥伦比亚罗萨里奥群岛海洋馆的18个个体及1个公共佛罗里达样本进行了全基因组分析，构建了参考引导的染色体规模基因组组装，并基于SNP解析了圈养种群的基因组结构。研究获得的关键结果为：基因组总长度1.12 Gb，contig N50为42.69 Mb，scaffold N50为46.30 Mb；共注释22,692个蛋白质编码基因；重复序列占基因组的46.4%。比较基因组分析表明，E. itajara与E. lanceolatus亲缘关系最近，两者间同义替代率（Ks）约0.01，推测分化时间约4百万年前。种群结构分析识别出四个遗传谱系，其中哥伦比亚加勒比地区存在两个特有谱系，且样本间呈现中度遗传分化（F_ST = 0.1823）。这项研究填补了该物种在哥伦比亚地区的基因组资源空白，为理解其进化历史、遗传多样性及繁殖机制提供了重要工具，并可直接服务于圈养繁殖的遗传管理。

**主要技术方法**
研究样本包括哥伦比亚罗萨里奥群岛海洋馆的18个个体，其中13个来自西斯帕塔湾（Bahía Cispatá），3个来自巴尔沃科阿斯湾（Bahía Barbacoas），1个来自阿拉尔卡（Ararca），以及1个圈养F2后代；另纳入1个来自美国佛罗里达州卡纳维拉尔港的公共测序数据。核心方法包括：PacBio HiFi与Illumina测序；使用Hifiasm进行初步组装，NGSEP进行错误校正，NextPolish2进行抛光；利用RagTag基于E. lanceolatus基因组进行参考引导的染色体支架装配；EDTA注释转座元件，MAKER注释基因；BUSCO、Merqury、QUAST等评估组装质量；NGSEP GenomesAligner和SynVisio进行共线性及比较基因组分析；STRUCTURE、主成分分析（PCA）和层次聚类解析种群结构；自定义Python脚本计算观测杂合度（H_o）、期望杂合度（H_e）、核苷酸多样性（π）、固定指数（F_ST）及亲缘系数（φ）。

**研究结果**
**3.1 大西洋石斑鱼的新基因组组装（A New Genome Assembly for the Goliath Grouper）**
通过单个体（ID 61517）的PacBio HiFi（约81×深度）和Illumina（约94×深度）测序，经Hifiasm组装获得220个contig，参考引导支架后获得168个scaffold，其中76个contig被合并为24条假染色体。BUSCO评估显示99.1%的完整放线鳍鱼单拷贝基因，Merqury一致性质量值（QV）为53.98，k-mer完整性98.72%，PacBio回比率为100%，表明组装具有高准确性和完整性。重复序列注释显示转座子占基组组的46.4%，其中DNA转座子最丰富（28.44%）。线粒体基因组全长16,752 bp，包含13个蛋白质编码基因、2个核糖体RNA基因和22个转运RNA。

**3.2 比较基因组学（Comparative Genomics）**
将E. itajara的基因组与8种石斑鱼基因组进行比较，发现染色体数目保守且多数物种呈现高度共线性。E. lanceolatus表现出最高的共线性，其Ks值约0.01，支持两物种最近亲缘关系。在染色体4、5和19上检测到重复出现的推定倒位信号，其中染色体4的倒位在全部配对比较中一致存在。基于Ks估算的物种分化时间显示，E. itajara与四个物种（E. fuscoguttatus、E. moara、E. polyphekadion、E. tukula）的分化约发生在10百万年前，而与另三个物种（E. akaara、E. awoara、E. cyanopodus）的分化约在20百万年前。对20,353个直系同源群的分析显示，其中8,685个为单拷贝基因家族，且存在多拷贝扩张事件。

**3.3 圈养大西洋石斑鱼的基因组多样性（Genomic Diversity Within Captive Goliath Groupers）**
基于19个个体共7,706个高质量SNP，STRUCTURE分析（Puechmaille估计量）确定最优K=4，PCA进一步支持该群体结构。结果显示四个遗传谱系，其中哥伦比亚样本中橙色和绿色谱系分别与西斯帕塔湾和巴尔沃科阿斯湾相关，但受限于各采样点样本量较小（如巴尔沃科阿斯湾仅3个样本），这些谱系的地理独特性仅被视为初步发现。遗传多样性指标显示，蓝色谱系（n=11）的期望杂合度（H_e=0.2361）和核苷酸多样性（π=0.2474）高于橙色谱系（n=5，H_e=0.1593，π=0.1770），表明蓝色谱系保留了更多遗传变异。橙色谱系内亲缘系数较高（φ=0.2048），接近一级亲缘关系预期值，提示该谱系内个体可能密切相关。蓝色谱系内个体间亲缘系数接近零，符合自然随机交配种群特征。两个主要谱系间F_ST为0.1823，呈中至高度分化，且各谱系拥有大量私有等位基因（蓝色谱系2,096个，橙色谱系371个），进一步支持限制性基因流的存在。

**讨论与结论**
讨论部分指出，E. itajara基因组大小（1.12 Gb）与E. malabaricus和E. lanceolatus相近，但蛋白编码基因数略少。共线性分析显示与E. lanceolatus高度保守，尽管参考引导支架可能在一定程度上增加保守性，但仅76个contig构建24条染色体表明结构信息可靠。种群结构方面，哥伦比亚样本内存在多个遗传谱系，暗示可能存在地理隔离或有限的连通性，但因采样不平衡，尚不能确定谱系是否地理专属。研究强调传统标记的局限性，而全基因组SNP提供了更精细的种群差异信息。保护策略上，研究人员建议通过被动声学监测和遥测技术识别繁殖聚群，建立受保护的捕鱼禁区和迁徙廊道；同时建议设定最小捕捞尺寸（如至少140厘米），并恢复该物种在IUCN红色名录中的“极危”（critically endangered）等级。该基因组框架为圈养繁殖提供了遗传指导，例如通过亲缘关系和多样性指标筛选配对方可最大化保留遗传多样性，并为未来野外恢复和遗传拯救项目奠定基础。总之，该研究首次提供了E. itajara的染色体规模参考基因组及哥伦比亚加勒比地区的种群基因组数据，揭示了圈养种群的遗传结构和多样性差异，为保护生物学和养殖遗传管理提供了关键资源。