南极棘鳍冰鱼的新分型组装为冰鱼辐射演化提供了新见解

A new phased assembly of the Antarctic spiny plunderfish provides novel insights into the evolution of the notothenioid radiation

作者信息Jacopo Martelossi, Ksenia Krasheninnikova, Amy Denton, Jonathan M D Wood, Thomas C Mathers, Richard Durbin, Nova Fong, David L Bentley, Melody S Clark, Iliana Bista
PMID42881419
发布时间2026-10-15
DOI10.1186/s13059-026-04293-7

实验完整度

高

包含新分型基因组组装、跨12种冰鱼比较基因组分析、afgp基因座单倍型解析、选择分析和重复元件分析等多个独立证据层级,并进行了机制层面的验证分析。

主要模型

Harpagifer antarcticus 分型基因组组装 12种冰鱼染色体水平基因组 H. antarcticus 鳃组织RNAseq样本

重点核对

H. antarcticus 基因组组装使用ONT长读长和Hi-C数据,两个单倍型总长度分别为1,069 Mb和1,222 Mb,BUSCO完整性99% afgp基因座位于16号染色体,Hap1和Hap2分别注释到17和16个afgp基因拷贝,以及各2个嵌合afgp/tlp基因 比较基因组分析使用12种冰鱼染色体水平基因组,包括新组装的H. antarcticus Hap1和11个公开组装 选择分析使用扩展数据集,额外加入7种温带或热带鱼类,测试冰鱼祖先分支上的多样化选择 LINE-L2分析基于Hap1和Hap2中长度大于1,000 bp的LINE-L2插入,提取逆转录酶结构域序列进行系统发育分析

摘要

背景:冰鱼是南极特有的、特征明确的物种群,也是研究基因组对极端寒冷适应的重要模式类群。一个新的参考组装和跨分支比较基因组分析,使我们能够研究冰鱼演化以及寒冷适应相关新功能的出现。结果:对模式冰鱼Harpagifer antarcticus的新分型组装显示,全基因组范围内单倍型变异水平较低。尽管如此,仍发现了来自多个LINE-L2分支的大量插入,提示持续转座并可能对物种形成有贡献。与预期相反,afgp基因座在单倍型之间高度相似,只是afgp基因存在大的长度等位变异。分析提示H. antarcticus中afgp基因座通过每次涉及两对afgp基因的片段串联重复进行扩张的模型。跨分支基因组的共线性重建展示了冰鱼中保守的大规模共线性关系以及类群特异性的染色体融合。对冰鱼多样化期间基因组获得和转座速率的量化显示,最初伴随历史温度下降出现了祖先的缓慢基因组扩张。随后出现了谱系特异性的基因组获得和转座活性的巨大峰值。最后,我们鉴定了一组在冰鱼出现期间经历了祖先多样化选择并获得新的保守非编码元件的基因。这些基因与抗氧化和蛋白质稳态相关,可能促进了冰鱼在南极的辐射。结论:多样化选择和与转座子活性相关的基因组获得是贯穿该分支的谱系特异性演化动态的主要贡献者,促进了对寒冷生活的适应。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
冰鱼辐射演化过程中,转座子活性、afgp基因座演化以及其他基因组特征如何促成对极端寒冷环境的适应?
核心机制
多样化选择与转座子活性相关的基因组获得是冰鱼谱系特异性演化动态的主要贡献者,祖先缓慢的转座子积累和近期谱系特异性转座爆发共同促进了寒冷适应;afgp基因座通过每次涉及两对基因的片段串联重复扩张。
主要证据
H. antarcticus分型基因组组装显示低单倍型变异但多个LINE-L2插入;afgp基因座在两个单倍型间结构几乎相同但存在长度等位变异;跨12种冰鱼比较基因组分析显示基因组获得远超丢失,且转座子活性与温度下降时间吻合;选择分析鉴定出334个在冰鱼祖先受选择且不在温带物种中受选择的基因,富集于氧化应激、脂质代谢和蛋白质稳态。
研究意义
该分型基因组组装为未来冷适应和气候变化适应性的功能研究提供了宝贵资源;研究表明基因组获得与转座子活性相关,并与南大洋温度下降相关,且冰鱼基因组在基因和控制区域均受到显著选择压力,使这些鱼能够在0°C以下生存。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

H. antarcticus 分型基因组组装与注释

获得首个冰鱼物种的单倍型解析染色体水平参考基因组,为后续比较基因组和基因座分析提供基础

使用ONT长读长和Hi-C数据组装H. antarcticus基因组,获得两个高质量单倍型;进行重复元件注释、基因注释、端粒和rDNA注释;评估组装完整性和杂合度

2

重复元件与LINE-L2家族分析

表征H. antarcticus基因组中转座子组成、活性历史以及LINE-L2家族的演化动态

使用RepeatModeler和RepeatMasker注释重复元件;提取LINE-L2逆转录酶结构域进行系统发育分析,鉴定高拷贝家族;重建家族共识序列并在12种冰鱼基因组中搜索;通过CpG校正的Kimura距离转换为绝对时间构建重复景观

3

单倍型变异与结构变异分析

评估H. antarcticus两个单倍型之间的序列和结构差异,识别近期活跃的转座子插入

使用Minimap2和Syri进行全基因组比对,鉴定SNP、结构变异和重排;使用GLMM分析SV与基因组特征的关联;通过ONT reads回贴和Sniffles2检测杂合SV;重叠RepeatMasker注释识别杂合TE插入

4

afgp基因座注释与演化分析

解析afgp基因座在H. antarcticus两个单倍型上的结构、基因拷贝数和演化历史

使用GeMoMa和BLASTN进行同源基因预测和手动注释;使用Klumpy验证注释完整性;基于exon-1和intron-1构建afgp和嵌合afgp/tlp基因的最大似然树;在多个冰鱼物种中自动注释afgp基因座

5

跨物种比较基因组与系统发育分析

重建冰鱼物种树和分化时间,量化基因组获得和丢失,分析共线性关系和染色体融合

使用12种冰鱼染色体水平基因组和7种温带/热带鱼类;提取单拷贝BUSCO基因进行系统发育分析;使用LSD2进行时间校准;使用GENESPACE和Progressive Cactus进行共线性分析和全基因组比对;使用HAL推断插入和缺失事件

6

选择分析和保守元件鉴定

鉴定冰鱼祖先分支上经历多样化选择的蛋白编码基因和新出现的保守非编码元件

使用OrthoFinder聚类直系同源组,DISCO分解多拷贝OG;使用PRANK进行密码子比对;使用aBSREL检测多样化选择;使用PhastCons预测保守元件;使用halAlignmentDepth确定CEs的起源;进行GO富集分析

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
Qubit 1X dsDNA高灵敏度检测试剂盒ThermoFisherQ33231
Qubit Flex荧光计ThermoFisherQ33327
Femto Pulse脉冲场电泳系统AgilentM5330AA
gDNA 165 kb分析试剂盒AgilentFP-1002-0275
连接法测序试剂盒V14--SQK-LSK114
Kapa Hyperprep UDI文库构建试剂盒KapaKK8544
无酚红L15培养基----
Nunclon组织培养板Nunclon--
TRI试剂----

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
基因组组装
ONT reads覆盖度约23.4X每单倍型;Hifiasm版本0.24.0-r702,使用--ont选项和Hi-C定相模式;Hi-C数据使用bwa-mem2比对,YaHS进行scaffolding
阅读提示:Methods中Genome assembly小节
RNAseq样本处理
三只H. antarcticus标本的鳃组织,在0°C和6°C分别孵育6天;每个温度一个样本;RNA提取使用TRI Reagent;1 µg rRNA depleted RNA用于文库构建;Novaseq PE150测序
阅读提示:Methods中Sequencing data generation小节
重复元件注释
RepeatModeler v.2.0.5使用-LTRStruct扩展;ProtExcluder v.1.2去除宿主基因;MCHelper v1.7.0默认参数;RepeatMasker sensitive模式
阅读提示:Methods中Repeat annotation小节
afgp基因座注释
使用D. mawsoni afgp基因座(accession HQ447059.1,haplotype 1)作为参考;GeMoMa v1.7.1默认参数;BLASTN v2.12.0手动注释;Klumpy验证;afgp基因座位于16号染色体约699 kb
阅读提示:Methods中Annotation of the afgp locus小节
比较基因组分析
使用12种冰鱼染色体水平基因组和7种温带/热带鱼类;单拷贝BUSCO基因使用MAFFT E-INS-i比对;TrimAl automated模式清理;1,000个比对随机选择;IQ-TREE v2.4.0约束树搜索;LSD2使用固定替换率进行时间校准,校准点:94 MYA、47 MYA、10.7 MYA
阅读提示:Methods中Species tree reconstruction and divergence time estimation小节
选择分析
OrthoFinder v2.5.4 ultra-sensitive模式聚类;DISCO分解多拷贝OG;PRANK v250331密码子模型比对;Gblocks和TrimAl过滤;OG比对长度>300密码子;aBSREL from HyPhy v2.5.8测试冰鱼祖先分支;排除在温带适应分支受选择的基因
阅读提示:Methods中Selection analyses小节
单倍型变异分析
Minimap2 v2.26全基因组比对;Syri v1.7.1检测SNP、SV和重排;GLMM使用lme4包,染色体为随机截距;ONT reads回贴使用Minimap2,Sniffles2检测杂合SV;SURVIVOR合并SV
阅读提示:Methods中Haplotype comparison小节