《Marine Genomics》：Genomic mechanisms of life history divergence shaped by demography and selection in the yellowtail kingfish species complex
编辑推荐：
摘要黄尾鰤（Seriola spp.）是上层硬骨鱼类（epipelagic teleosts），在全球水产养殖中的重要性日益增加。该物种复合体（species complex）包括东北太平洋的S. dorsalis、西北太平洋的S. aureovittata，
摘要黄尾鰤（Seriola spp.）是上层硬骨鱼类（epipelagic teleosts），在全球水产养殖中的重要性日益增加。该物种复合体（species complex）包括东北太平洋的S. dorsalis、西北太平洋的S. aureovittata，以及南太平洋和南部非洲的两个不同谱系的S. lalandi。虽然其遗传分化程度已有充分记录，但潜在的基因组机制及其与生活史变异的功能相关性仍知之甚少。因此，研究人员开展了一项比较基因组分析，其中包含一个组装的南非S. lalandi染色体水平参考基因组（678.3 Mb，21,796个蛋白编码基因，27.0%重复元件）。种群统计重建（demographic reconstructions）表明，南非在更新世（Pleistocene）期间充当了气候避难所（climate refugium），有助于保留支持活跃远洋生活方式（pelagic lifestyle）的祖先基因组架构。同样，来自中国的S. aureovittata基因组保留了若干与持续运动相关的基因，但其染色体重排（chromosomal rearrangements）数量增加，这可能归因于末次冰期（Last Glacial Period）期间的种群收缩。相比之下，澳大利亚的S. lalandi，尤其是加利福尼亚的S. dorsalis，表现出更高的基因家族周转（gene family turnover），表明生态分化更大。控制代谢、免疫和发育通路的基因中的选择信号（signatures of selection）表明，生活史的分化受对比鲜明的热和氧气条件以及病原体压力的影响。这些发现对资源管理（包括杂交育种潜力以及在持续人为变化下的适应能力）具有广泛意义。
**论文解读：黄尾鰤物种复合体生活史分化的基因组机制**
**研究背景、问题与目的**
黄尾鰤（Seriola spp.）是一类大型远洋性硬骨鱼类，在全球水产养殖中地位日益突出。该物种复合体包括东北太平洋的加州黄尾鰤（S. dorsalis）、西北太平洋的黄尾鰤（S. aureovittata）以及南太平洋和南部非洲的拉氏黄尾鰤（S. lalandi，含两个不同谱系）。尽管已有研究明确了这些谱系间的遗传分化程度，但对于背后的基因组机制及其与生活史变异的关联仍知之甚少。特别是，南非谱系缺乏参考基因组，限制了对整个复合体性状分化的比较。为此，该研究通过建立南非谱系染色体水平的参考基因组，并与太平洋谱系进行对比，旨在评估种群历史及基因组机制如何导致黄尾鰤间的表型分化。
**主要关键技术方法**
研究人员采用PacBio HiFi（17×）和Illumina双端测序（38×）对南非False Bay野生捕获的雌雄个体进行测序，利用hifiasm进行从头组装，并通过RagTag参考引导支架和Pilon碱基校正获得染色体级基因组。重复元件识别使用RepeatModeler和RepeatMasker；基因预测采用BRAKER v3.0.8，整合蛋白证据（OrthoDB及Seriola其他基因组）和RNA-seq（10个智利S. lalandi个体）。比较基因组分析包括：skani计算平均核苷酸同一性（ANI），MUM&Co鉴定结构变异，OrthoFinder进行直系同源聚类，CodeML分支位点模型检测正选择，PSMC进行种群历史重建。
**研究结果**
**3.1 基因组组装**：南非S. lalandi（SA）基因组最终大小为678.3 Mb，97.4%锚定至24条伪染色体，BUSCO完整度98.5%，MERQURY QV值47.0（HiFi）和41.3（Illumina）。与太平洋谱系相比，SA组装质量与S. aureovittata（CND、CNY）相当，优于短读组装的澳大利亚（AUS）和加州（USA）基因组。
**3.2 基因组注释**：SA基因组中重复元件占27%，以DNA转座子为主（14%）。预测21,796个蛋白编码基因，94.9%获得功能注释；AUS仅19,020个基因，蛋白完整度84.9%低于SA的96.2%。
**3.3 比较基因组分析**
**3.3.1 基因组核苷酸差异**：种间差异1.74%-3.05%，黄尾鰤复合体内，半球间差异（S. lalandi vs S. aureovittata）0.77%-0.85%，半球内差异（USA vs S. aureovittata）0.61%-0.65%，S. lalandi同种间0.38%，S. aureovittata间0.22%。
**3.3.2 密码子使用偏好**：所有基因组平均有效密码子数（ENC）约50，呈中度偏好。SA中非适应性组分压力（如GC偏向基因转换）对密码子使用的贡献最大（R2=0.484），而北太平洋谱系中翻译选择作用相对更强。
**3.3.3 直系同源聚类**：共鉴定20,344个蛋白簇，69.4%为所有基因组共享；149个谱系特异性簇（0.7%）。SA和AUS共享389个簇，北太平洋基因组共享771个簇。
**3.3.4 谱系特异性基因富集**：SA富集肌球蛋白重链基因（myh1、myh7、myo6）和运动活动；北太平洋基因组富集先天免疫相关NLR家族和TRIM基因，以及嗅觉受体；USA独特富集精子发生、心血管发育、低氧信号等基因；CNY富集肾脏发育基因。
**3.3.5 基因组结构变异**：以SA为参考，太平洋谱系结构变异数量为7,612-13,024，其中插入/缺失占95%，易位和倒位占0.3%。倒位长度在谱系间差异显著（p=0.038），北太平洋S. aureovittata的倒位和易位总数最多。
**3.3.6 染色体重排的功能注释**：倒位涉及免疫基因（nlrp12、nlrc3、gimap8）、热应激反应（hspa1）、低氧反应（hyou1）等；易位涉及细胞骨架动态和细胞粘附基因。
**3.4 系统基因组学、选择与种群历史**
**3.4.1 系统发育与分化时间**：基于11,090个单拷贝直系同源基因构建的超度量树显示，半球间分化最早发生在2.43百万年前（mya），S. lalandi谱系分离在1.21 mya，北太平洋谱系内CND与CNY在1.05 mya分离。
**3.4.2 基因家族动态与富集**：与S. quinqueradiata相比，黄尾鰤共同祖先经历显著收缩（+23/-343）。SA周转最小（+141/-270），主要集中在钙离子转运和血管舒张基因；USA周转最大（+411/-859），富集信号转导、细胞外基质和先天免疫基因。
**3.4.3 正选择基因**：6,195个单拷贝基因中12.9%显示正选择，92.1%为谱系特异性。各谱系富集类别包括：感官知觉、心血管功能、钙离子稳态、生殖过程、脂质代谢、免疫应答等。NLR、TRIM、GIMAP、MUC家族基因在多个谱系中受选择。
**3.4.4 种群历史**：PSMC重建显示，早更新世SA和CND有效种群大小较大，中更新世转折期所有谱系下降，末次冰期北太平洋谱系（CND、CNY）急剧下降，而USA持续增长，SA相对稳定。
**讨论与结论**
**讨论总结**：谱系分化主要受晚更新世种群事件驱动，南非作为气候避难所保留了祖先基因组架构，支持活跃的远洋生活方式；北太平洋谱系因末次冰期种群瓶颈积累了更多染色体重排，加州谱系通过快速扩张和局部适应导致基因家族扩张和功能创新。代谢、免疫、发育通路的选择信号表明生活史分化受热和氧气条件及病原体压力影响。结构变异（尤其是倒位）有助于保护适应性等位基因组合，但可能增加杂交衰退风险。南非谱系因稳定的种群历史和祖先性状保留，可能更适合作为杂交育种受体，但需标记辅助选择管理染色体不相容性。
**结论部分翻译**：比较基因组分析支持以下假说：黄尾鰤复合体间的分化主要是晚更新世种群事件的结果，这些事件差异性地调节了遗传变异，而持续的生态压力则影响生活史性状。可以推测，当代热和氧气条件的变化以及病原体暴露强烈影响了观察到的表型差异。南非谱系经历了更稳定的种群历史，保留了活跃的、泛化的远洋生活方式，预计能更稳健地应对更广泛的生态压力，并可能成为水产养殖杂交育种的理想受体。从分化谱系中渗入适应性性状可能提供利用杂种优势的途径，但需要标记辅助选择来管理染色体不相容性。此外，必须采取措施防止逃逸对当地适应性野生谱系遗传完整性的影响。然而，克服当前组装限制并纳入南太平洋其他同种个体将是充分阐明黄尾鰤复合体分化的必要条件。此外，已建立的南非参考基因组和分化基因目录为未来种群和数量遗传学研究提供了基础，有助于遗传改良和种群管理。
