Chromosome-level genome and characteristic analysis of Platax teira
-
摘要:
尖翅燕鱼 (Platax teira) 具有生长速度快、肉质鲜美、营养价值高等特点,其外形奇特,尤其幼鱼更为特殊,可作为观赏鱼,是南海发展网箱养殖的重要潜力鱼类之一。由于目前缺乏尖翅燕鱼的基因组信息,其多数功能基因未被挖掘,已成为制约其遗传育种的重要因素。运用三代测序技术和组装获得染色体水平的尖翅燕鱼高质量基因组图谱,通过基因组注释获得尖翅燕鱼基因组序列的基本生物学信息。结果显示,尖翅燕鱼基因组大小为697.98 Mb,Hi-C挂载至24条染色体,挂载率为99.26%,含重复序列177.79 Mb,占整个基因组的25.47%,注释到22 851个蛋白编码基因。系统进化分析显示尖翅燕鱼与波纹唇鱼 (Cheilinus undulatus) 亲缘关系最近,分化时间距今约82.89 Ma。正选择基因主要富集在与离子通道和心脏功能有关的通路中,扩张的基因家族主要富集在嗅觉传导和氮代谢通路上,揭示了其在特定环境中的生存、适应依据和生态适应策略。
Abstract:Platax teira has the characteristics of fast growth rate, delicious meat and high nutritional content, and its strange appearance, especially for young fish, makes it an ornamental fish, being one of the important potential fishes for the cage culture development in the South China Sea. Due to the lack of genomic information, most of the functional genes of P. teira have not been explored, which has become an important factor of restricting its genetic breeding. We utilized triple sequencing technology and assembly to obtain a high-quality genome map of P. teira on chromosome level, and obtained the basic biological information of P. teira genome sequence through genome annotation. The results show that the genome size of P. teira was 697.98 Mb, assembling into 24 chromosomes with an assembly rate of 99.26%. P. teira genome contained 177.79 Mb of repetitive sequences, accounting for 25.47% of the total genome and encoding 22 851 genes. Comparative genomic analysis with 11 other fish species reveals that P. teira shared the closest relationship with Cheilinus undulatus, and the differentiation time was about 82.89 Ma. Genes under positive selection in P. teira were enriched in pathways related to ion channels and cardiac function, while the expanded gene families were enriched in pathways related to olfactory transmission and nitrogen metabolism, which reveals its survival, adaptation basis and ecological adaptation strategies in specific environment.
-
Keywords:
- Platax teira /
- Chromosome /
- Genome /
- Comparative genomics /
- Ephippidae
-
河口盐沼湿地初级生产力高、生境异质性强的生态环境为众多水生生物提供了丰富的饵料资源和栖息生境[1],因此被认为是鱼类重要的育幼场所[2-4]。长江口湿地属于河口型冲击沙洲,在陆海物质交汇、咸淡水混合、径流和潮汐的相互作用下,形成了特殊而又多样的生态环境,主要有潮沟、盐沼植被和光滩3种亚生境类型[5],它们为水生生物提供了产卵场、育幼场、索饵场和洄游通道[6],且对维持流域和近海鱼类生物多样性、保障渔业资源有效补充方面发挥着重要作用。近年来长江口盐沼受到了众多人为因素的干扰并已造成严重影响,如水域污染引起鱼类死亡、渔获量下降以及饵料生物减少;滩涂围垦等各种水利工程建设使栖息生境和育幼场所破碎化;过度捕捞和毁灭性渔业生产方式对亲鱼和幼鱼资源造成了严重破坏[7]。中国花鲈 (Lateolabrax maculatus) 也是深受其害的鱼类之一。
中国花鲈广泛分布于我国沿海及河口,是长江口的重要经济鱼类[7],对河口进行选择性利用。夏、秋季活跃于河口近岸水域,11月至翌年1月到外海产卵[8]。国内外学者已对中国花鲈的分类学、人工繁育、养殖、营养学、食性等进行了大量相关研究[9-14],但有关其如何选择、利用育幼生境的研究尚不多见。
生态位宽度是指一个种群在一个群落中所利用的各种资源的总合[15]。时间生态位宽度反映物种在时间尺度上对资源利用的情况;空间生态位宽度反映物种在空间维度上的分布范围,其大小反映物种对不同空间资源利用的能力大小[16]。本研究根据在长江口九段沙湿地3种亚生境的调查资料,分析中国花鲈的生态位、栖息分布以及种群特征的变化情况,旨为丰富中国花鲈对长江口盐沼湿地的利用提供基础资料,也为中国花鲈资源保护和可持续利用提供科学依据,并为长江口盐沼湿地生态修复提供新视角。
1. 材料与方法
1.1 样品采集与测量
研究地点位于上海九段沙湿地国家级自然保护区 (121°50'E—122°05'E、31°03'N—31°38'N,图1),受亚热带季风气候的影响,四季分明,冬季寒冷干燥,夏季暖热多雨,年均降水量为1 145 mm,年平均气温为15.7 ℃,1月最低 (平均4.2 ℃),7月最高 (平均27.3 ℃)。九段沙附近潮汐属于非正规半日浅海潮,平均潮差为2.67 m,最大潮差为4.62 m。冬季为枯水季节 (11月—翌年1月),被高盐度咸水控制,夏季为丰水季节 (7—9月),被低盐度江水控制[17]。
本研究于2019年5月 (春季)、8月 (夏季)、11月 (秋季) 和2020年1月 (冬季) 利用插网在光滩、海三棱藨草 (Scirpus mariqueter) 植被区和潮沟中各设一站点进行调查。插网由网目为15 mm的网翼和网囊组成,其中网翼长5 m、高1 m,网口1 m×1 m,网囊长6 m。每次调查插网连续收集3 d渔获物,当场冰上冷藏并带回实验室,分类鉴定到种。对渔获物进行生物学测量,包括全长、体长 (精确至1 mm)、体质量 (精确至0.01 g) 等。根据孙帼英等[8]研究并结合解剖结果辨别中国花鲈的生长发育情况。
1.2 数据分析
采用Pinkas相对重要性指数 (Index of relative importance, IRI) 来分析优势种[18],公式为:
$$ {\rm{IRI}}=(N{\text{%}}+W{\text{%}})\times F{\text{%}} $$ (1) 式中
$ N{\text{%}} $ 为某种类的尾数在渔获物总尾数所占百分比;$ W{\text{%}} $ 为某种类的质量在渔获总质量所占百分比;$ F{\text{%}}$ 为某种类出现的站位数占调查站位总数百分比,将IRI>500定义为优势种[19]。应用Shannon-Winner公式计算生态位宽度[20]:
$$\begin{split} &\qquad \qquad \qquad \qquad \qquad \;\\ &{B}_{{i}}=-{\sum }_{j=1}^{r}\left({P}_{i_{\dot{J}}}\mathrm{ln}{P}_{i{j}}\right) \end{split}$$ (2) 式中
${P}_{i_{\dot{J}}}$ =$ {n}_{ij}/{N}_{ij} $ 和$ {P}_{kj}={n}_{kj}/{N}_{kj} $ 分别代表种$ i $ 和$ k $ 在第$ j $ 个资源状态下的个体数占该种所有个体数的比例,$ j $ 表示某个调查航次。生态位宽度反映物种对资源的利用能力,空间生态位宽度可反映某物种分布区域广泛与否,时间生态位宽度可反映物种时间上出现的长短[21]。采用摄食等级来描述摄食强度[22-23]。2. 结果
2.1 研究区域的环境特征
调查期间,九段沙3种亚生境的水温随季节发生明显变化,水温夏季最高、冬季最低,3种亚生境间水温无显著性差异 (P>0.05);盐度则相反,夏季最低,冬季最高,3种亚生境间盐度无显著性差异 (P>0.05);3种亚生境最大淹没水深在夏季均最高,其他月份相差不大,潮沟比海三棱藨草、光滩的更高且差异显著 (P<0.05,表1)。春季至夏季海三棱藨草生长旺盛,秋季至冬季海三棱藨草逐渐枯萎、消亡。
表 1 九段沙湿地不同亚生境水温、盐度和水深的季节变化Table 1. Seasonal variation of water temperature, salinity and water depth in different subhabitats of Jiuduansha Wetland季节
Season水温
Water temperature/℃盐度
Salinity最大水深
Maximum water depth/m潮沟
Tidal creek海三棱藨草
S. mariqueter光滩
Mudflat潮沟
Tidal creek海三棱藨草
S. mariqueter光滩
Mudflat潮沟
Tidal creek海三棱藨草
S. mariqueter光滩
Mudflat春季 Spring 13.5 16.1 13.6 1.83 1.46 1.49 3.5 1.6 1.7 夏季 Summer 28.2 28.3 28.5 1.04 0.43 0.81 4.3 2.1 2.2 秋季 Autumn 14.7 14.8 14.5 7.24 7.12 7.77 3.4 1.5 1.6 冬季 Winter 7.6 8.1 7.9 14.88 12.34 14.51 3.1 1.4 1.6 2.2 长江口九段沙盐沼湿地主要种及其生态位
调查期间共采集鱼类样本8目12科,29种,1 346尾。3种亚生境数量百分比N>1%的鱼类有5种:中国花鲈、斑尾刺虾虎鱼 (Acanthogobius ommaturus)、鮻 (Liza haematocheila)、小黄鱼 (Larimichthys polyactis)、黄姑鱼 (Nibea albiflora)。不同物种在3种亚生境的优势度有差异,中国花鲈、斑尾刺虾虎鱼在3种亚生境中均为优势种 (IRI>500),鮻、小黄鱼、黄姑鱼在部分亚生境中为优势种 (表2)。
表 2 九段沙3种亚生境主要种的相对重要性指数Table 2. Relative importance index of main species in three subhabitats of Jiuduansha Wetland种类
Species潮沟
Tidal creek海三棱藨草
S. mariqueter光滩
MudflatN/% W/% F/% IRI N/% W/% F/% IRI N/% W/% F/% IRI 中国花鲈 L. maculatus 26.5 24.8 100 5 130 35.9 40.5 100 7 645 17.7 43.5 100 6 119 斑尾刺虾虎鱼 A. ommaturus 23.8 63.2 100 8 698 8.3 7.9 75 1 216 2.8 9 100 1 180 鰉 L. haematocheila 1.5 4.9 75 480 20.5 40.7 75 6 117 6.1 28.3 100 3 444 小黄鱼 L. polyactis 40 1.6 25 1 040 12.8 0.2 25 326.6 17.1 0.8 50 897 黄姑鱼 N. jordan 1.6 0.1 25 43 1.9 0.02 25 48.6 22.3 5.7 25 700 中国花鲈的空间生态位宽度和时间生态位宽度变化较大,空间生态位宽度介于0.16~0.88,在春、夏、秋季依次递减,且春季和夏季时空间生态位宽度比其他优势种的大,春、夏季中国花鲈在3种亚生境的分布较秋、冬季更均匀;中国花鲈的时间生态位宽度介于0.39~0.78,在时间尺度上中国花鲈对潮沟的利用程度大于海三棱藨草植被区和光滩亚生境 (表3)。
表 3 九段沙3种亚生境主要鱼种生态位宽度Table 3. Niche breadth of main fish species in three subhabitats of Jiuduansha Wetland种类
Species空间生态位宽度
Spatial niche breadth时间生态位宽度
Temporal niche breadth春季
Spring夏季
Summer秋季
Autumn冬季
Winter潮沟
Tidal creek海三棱藨草
S. mariqueter光滩
Mudflat中国花鲈 L. maculatus 0.88 0.52 0.16 0.26 0.78 0.66 0.39 斑尾刺虾虎鱼 A. ommaturus 0.17 0.41 0.45 0.51 0.97 0.22 0.18 鰉 L. haematocheila 0.57 0.21 0.45 0.89 0.27 0.90 0.72 小黄鱼 L. polyactis 0.57 0 0.01 0 0.16 0.14 0.28 黄姑鱼 N. albiflora 0 0.57 0 0 0.29 0.11 0.17 2.3 长江口中国花鲈分布特征
本研究共采集到中国花鲈449尾,3种亚生境数量差异明显,其中潮沟266尾,海三棱藨草植被区121尾,光滩62尾。各亚生境中国花鲈的数量在不同季节间变化显著,春、夏、秋和冬季各亚生境中国花鲈的数量先逐渐减少、后增加,春季最多,秋季最少 (图2)。
中国花鲈在3种亚生境的数量与总数量的比例随季节发生明显变化,潮沟中的数量占比在春、夏、秋、冬季逐渐减小;海三棱藨草和光滩的数量占比在春、夏、秋、冬季逐渐增加 (图3)。
2.4 中国花鲈的体长特征
中国花鲈性成熟时体长约350 mm,本研究捕获的中国花鲈体长最大为339 mm,最小为20 mm,经对比和解剖发现所捕获的中国花鲈均为幼鱼。体长小于70 mm的个体最多,占比94.24%,且该体长规格的个体数在不同亚生境中的比例存在差异,占比大小依次为潮沟>海三棱藨草>光滩。
中国花鲈的平均体长在3种亚生境均表现出春、夏、秋季逐渐增大的趋势,海三棱藨草亚生境中的平均体长在春季与夏季之间的变化幅度明显大于光滩和潮沟 (图4)。
各季节各亚生境中国花鲈平均体长的标准偏差较大,个体间的体长差异大。为进一步分析各亚生境间的体长差异,本研究将中国花鲈的体长数据进行合理分组。结果显示,春季中国花鲈体长主要集中在10~50 mm体长组,对该体长区间的3种亚生境体长数据进行单因子方差分析发现,潮沟、海三棱藨草和光滩亚生境中国花鲈的平均体长分别为34.13、29.27和32.15 mm,海三棱藨草亚生境中国花鲈的体长与其他亚生境有显著性差异 (P<0.05),潮沟亚生境与光滩亚生境的中国花鲈体长无显著性差异 (P>0.05);夏季各亚生境体长为91~130 mm的个体占比均最大,此体长组潮沟、海三棱藨草和光滩亚生境的平均体长分别为101.85、113.17和112.4 mm,潮沟亚生境中的中国花鲈体长与海三棱藨草亚生境和光滩亚生境差异显著 (P<0.05);秋季各亚生境体长分布较为分散,未发现体长小于90 mm的个体,各亚生境中的中国花鲈体长无显著性差异 (P>0.05);冬季中国花鲈体长集中在131~210 mm,未发现体长大于251 mm的个体 (图5)。
春、夏、秋季除了占比较大的小规格个体外,还存在大于251 mm的大规格个体,且此规格的个体在海三棱藨草和光滩亚生境出现的比例更高,但冬季3种亚生境中均无此规格的个体出现。
2.5 各亚生境中国花鲈平均摄食强度变化
中国花鲈的平均摄食等级在潮沟亚生境呈现从春季到冬季逐渐增大的趋势;在光滩亚生境呈现从春季到秋季逐渐增大而冬季骤减的趋势;在海三棱藨草亚生境呈现夏季最大、秋冬季逐渐减小的趋势。
中国花鲈平均摄食等级在春、夏、秋季表现为海三棱藨草>潮沟>光滩;冬季较其他季节特殊,平均摄食等级表现为潮沟>海三棱藨草>光滩 (图6)。
3. 讨论
3.1 优势种与生态位
本研究发现长江口盐沼湿地优势种为中国花鲈、斑尾刺虾虎鱼、鮻和小黄鱼,且中国花鲈优势度比其他优势种大,说明中国花鲈在九段沙湿地的群落中起主导和控制作用。本研究与杨刚等[24]在长江口的调查结果——优势种为中国花鲈、鮻和窄体舌鳎 (Cynoglossus gracilis)及庄平等[25]的调查结果——优势种为长蛇 (Saurogobio dumerili)、中国花鲈、窄体舌鳎、刀鲚 (Coilia nasus)、鮻和鲻 (Mugil cephalus) 略有差异。导致这种差异主要有2个原因:1) 时间。本研究进行4次采样,对应春、夏、秋、冬4个季节,较全面地反映了中国花鲈在九段沙湿地的动态变化,杨刚等[24]采样时间是2010年5月、8月、11月,庄平等[25]采样时间是2004—2006年每年的5—8月,与本研究采样时间略有差异。2) 生境。本研究采样地点包括3种不同的亚生境,能更客观地反映九段沙盐沼湿地鱼类群落特征,而杨刚等[24]、庄平等[25]对采样地点未作详细划分,可能因某些物种的选择偏好性而低估了盐沼湿地鱼类群落的物种丰度[26]。
本研究中中国花鲈的空间生态位宽度在春、夏、秋季逐渐减小,秋季的较小是因为光滩亚生境中国花鲈的数量偏少,为潮沟亚生境和海三棱藨草亚生境中数量的1/2,说明中国花鲈随季节变化在3种亚生境中的分布差异增大,这可能是因为季节变化引起水温、水位、盐度等环境因素的变化,进而影响各亚生境中的环境和饵料生物。春、夏季中国花鲈的空间生态位较秋、冬季大,说明春、夏季中国花鲈在3种亚生境的分布较秋、冬季更均匀,此时中国花鲈对生境的选择利用差异较小,侧面反映了春、夏季更适合中国花鲈在河口盐沼湿地育幼和生长。
时间生态位宽度反映物种在时间尺度上对资源利用的情况[16,27],中国花鲈时间生态位宽度在潮沟、海三棱藨草、光滩亚生境间逐渐减小。光滩亚生境中国花鲈的时间生态位最小,说明光滩亚生境中的中国花鲈在时间上较其他亚生境的分布差异大。潮沟亚生境中国花鲈的时间生态位宽度较其他亚生境大,说明潮沟亚生境中的中国花鲈在时间尺度上较其他亚生境的分布差异小,这可能是因为潮沟亚生境环境复杂,相较其他亚生境,季节变化对潮沟的干扰较小。
3.2 九段沙湿地中国花鲈的时空分布差异
河口洄游性鱼类一般受到自身生活史及生物学特性的影响,其分布具有明显的时空交替变化特征。本研究结果表明中国花鲈作为近海咸淡水鱼类,在不同亚生境和不同季节均存在数量差异 (潮沟的数量>海三棱藨草>光滩),春季数量最多。中国花鲈为肉食性鱼类,具有主动捕食的特征,其季节分布差异受自身生物学特性和饵料生物的影响[28-29]。作为洄游性鱼类,中国花鲈在秋季至冬季性成熟,洄游到水位深、盐度高的水域产卵,随后仔、稚鱼在早春游到近岸低盐度水域索饵[7],因此,本研究调查到的数量表现出春季最多、随后逐渐减小、冬季有所回升的趋势。春季,大量中国花鲈幼体集中游向河口盐沼湿地索饵、躲避敌害,此区域的低盐度和丰富的适口饵料生物适合其在此栖息[30-31];夏季,中国花鲈逐渐长大,游泳能力增强,食性也发生转变,为满足生长需求会游向其他水域以获取更多的食物[14,32-33];秋季,成体和亚成体需游向深水位、高盐度水域为繁育做准备,因此秋季在盐沼湿地捕到的中国花鲈个体很少;冬季,产过卵的中国花鲈会游到河口盐沼湿地索饵,这可能是冬季捕获个体数量多于秋季的原因。有研究表明长江口潮间带是中国花鲈的育肥场所[8],中国花鲈在长江口3种亚生境的空间分布差异,可能与不同亚生境的生境特质有关。潮沟和海三棱藨草植被区比光滩亚生境有更丰富的底栖生物,3种亚生境浮游生物优势种和丰度也存在差异[34],且本研究发现3种亚生境中中国花鲈的摄食等级有明显差异,因此,中国花鲈数量在长江口盐沼湿地的时空变化可能是基于其自身生物学特性以及生境间的差异、季节性利用河口育幼功能的结果。
3.3 中国花鲈的体长特征
国内外的一些研究表明,盐沼湿地是一些近岸渔业物种的重要栖息生境[4,35]。张衡和朱国平[36]通过在长江口的调查推测长江口盐沼湿地是重要的育幼场,金斌松[37]在崇明东滩和九段沙湿地调查发现潮沟中存在大量的中国花鲈幼鱼,平均体长为24 mm。本研究春季中国花鲈的平均体长为43.6 mm,经解剖并与孙帼英等[8]研究对比发现,本研究所捕中国花鲈均为幼鱼;3种亚生境中国花鲈体长小于50 mm数量占比均达86%以上,幼体比例较高,说明长江口盐沼湿地是中国花鲈的育幼场。中国花鲈在3种亚生境的体长存在差异,尤其在春季和夏季。相较于潮沟,大规格个体在海三棱藨草、光滩出现的比例更高,这可能因为海三棱藨草、光滩亚生境较为开阔,更利于中国花鲈追捕游泳动物。
-
表 1 尖翅燕鱼基因组组装和校正后组装结果统计
Table 1 Statistical analysis of genome assembly and corrected assembly results of P. teira
类型
TypeNextDenovo软件 NextPolish软件 长度
Length/bp数量
Count长度
Length/bp数量
CountN50 26 172 673 12 26 178 423 12 N90 8 653 897 29 8 656 320 29 最小长度
Min. length11 814 11 779 最大长度
Max. length33 665 020 33 669 695 平均长度
Ave. length5 207 983 5 208 839 总长度
Total length697 869 769 134 697 984 514 134 表 2 尖翅燕鱼BWA比对结果
Table 2 Results of BWA comparison of P. teira
总读长
Total reads匹配读长
Map reads匹配率
Map rate/%双端读长
Paired reads双端比对读长
Paired map reads真正的读长
Proper paired reads正确匹配率
Properly map rate/%758 790 194 755 705 560 99.59 701 824 260 698 739 626 671 793 148 95.72 表 3 尖翅燕鱼基因组BUSCO评估
Table 3 Genome BUSCO evaluation of P. teira
类型
Type数量
Number百分比
Percentage/%完整的BUSCOs
Complete BUSCOs3 595 98.8 完整的单拷贝BUSCOs
Complete and single-copy BUSCOs3 571 98.1 完整的重复序列BUSCOs
Complete and duplicated BUSCOs24 0.7 碎片BUSCOs
Fragmented BUSCOs9 0.2 未比对上的BUSCOs
Missing BUSCOs36 1.0 总的BUSCO
Total BUSCO3 640 100.0 表 4 尖翅燕鱼各染色体长度
Table 4 Length of each chromosome of P. teira
染色体
Chromosome长度
Length/bp染色体
Chromosome长度
Length/bpHiC_scaffold_1 32 115 465 HiC_scaffold_13 22 937 044 HiC_scaffold_2 28 296 500 HiC_scaffold_14 29 051 425 HiC_scaffold_3 30 192 628 HiC_scaffold_15 28 201 383 HiC_scaffold_4 24 662 710 HiC_scaffold_16 23 075 928 HiC_scaffold_5 16 755 790 HiC_scaffold_17 25 711 309 HiC_scaffold_6 29 747 999 HiC_scaffold_18 34 711 268 HiC_scaffold_7 31 770 170 HiC_scaffold_19 31 217 840 HiC_scaffold_8 30 909 847 HiC_scaffold_20 34 223 832 HiC_scaffold_9 27 971 542 HiC_scaffold_21 32 267 482 HiC_scaffold_10 28 320 680 HiC_scaffold_22 30 609 981 HiC_scaffold_11 34 468 000 HiC_scaffold_23 27 290 755 HiC_scaffold_12 26 316 500 HiC_scaffold_24 32 023 099 合计 Total 692 849 177 表 5 尖翅燕鱼重复序列注释结果
Table 5 Results of repeated sequence annotation of P. teira
元件类型
Elements type元件数量
Number of
elements长度
Length/bp百分比
Percentage/%逆转录因子
Retroelements137 738 39 456 150 5.65 DNA转座子
DNA transposons266 011 55 441 011 7.94 环状
DNA Rolling-circles2 389 682 132 0.10 无分类
Unclassified366 703 66 929 139 9.59 总夹杂的重复序列
Total interspersed
repeats161 826 300 23.18 小RNA
Small RNA700 80 729 0.01 卫星DNA
Satellites DNA660 189 494 0.03 简单的重复序列
Simple repeats317 455 12 696 083 1.82 低复杂性
Low complexity42 640 2 317 723 0.33 总的重复
Total repeats177 792 461 25.47 表 6 尖翅燕鱼结构注释结果BUSCO评估
Table 6 BUSCO evaluation of structural annotation results of P. teira
类型
Type数量
Number百分比
Percentage/%完整的BUSCOs
Complete BUSCOs3 213 88.3 完整的单拷贝BUSCOs
Complete and single-copy BUSCOs3 176 87.3 完整的重复序列BUSCOs
Complete and duplicated BUSCOs37 1.0 碎片BUSCOs
Fragmented BUSCOs155 4.3 未比对上的BUSCOs
Missing BUSCOs272 7.4 总的BUSCO
Total BUSCO3 640 100.0 表 7 尖翅燕鱼功能注释结果
Table 7 Functional annotation results of P. teira
数据库
Database数量
Number百分比
Percentage/%NR 21 016 91.97 SwissProt 18 945 82.91 KEGG 17 511 76.63 COG 5 787 25.32 Trembl 20 996 91.88 功能注释基因
Functional annotation21 104 92.35 总计Total 22 851 100.00 表 8 12个物种基因家族聚类结果
Table 8 Results of gene family clustering in 12 species
物种
Species总的
Total单拷贝
Single特异性
Specific未聚类
Unclustered眼斑双锯鱼A. ocellaris 23 035 4 787 226 185 泰国斗鱼B. splendens 22 791 4 787 245 227 大白鲨C. carcharias 19 440 4 787 1 252 1 273 波纹唇鱼C .undulatus 23 316 4 787 381 184 斑马鱼D. rerio 32 717 4 787 3 503 1 659 云纹石斑鱼E .moara 23 735 4 787 112 336 食蚊鱼G. affinis 23 135 4 787 98 166 青鳉O. latipes 22 071 4 787 460 241 尖翅燕鱼P .teria 22 851 4 787 689 2 727 深裂眶锯雀鲷S. partitus 22 589 4 787 70 363 红箭鱼X. helleri 23 921 4 787 98 138 花斑剑尾鱼X. maculatus 23 238 4 787 34 126 -
[1] BRAY R A, CRIBB T H. Lepocreadiidae (Digenea) from the batfish of the genus Platax Cuvier (Teleostei: Ephippidae) from the southern Great Barrier Reef, Queensland, Australia[J]. Syst Parasitol, 2003, 55(1): 1-9. doi: 10.1023/A:1023974022432
[2] MARIMUTHU N, WILSON J J, KUMARAGURU A K. Teira batrish, Platax teira (Forsskal, 1775) in Pudhumadam coastal waters, drifted due to the tsunami of 26 December 2004[J]. Current Sci, 2005, 89(8): 1310-1312.
[3] GOLANI D, SONIN O, EDELIST D. Second records of the Lessepsian fish migrants Priacanthus sagittarius and Platax teira and distribution extension of Tylerius spinosissimus in the Mediterranean[J]. Aquat Invasions, 2011, 6(S1): S7-S11.
[4] 刘明鉴, 郭华阳, 高杰, 等. 尖翅燕鱼早期胚胎发育及仔稚鱼形态观察[J]. 南方水产科学, 2022, 18(4): 103-111. doi: 10.12131/20210251 [5] BILECENOGLU M, KAYA M. A new alien fish in the Mediterranean Sea-Platax teira (Forsskål, 1775) (Osteichthyes: Ephippidae)[J]. Aquat Invasions, 2006, 1(2): 80-83. doi: 10.3391/ai.2006.1.2.5
[6] LEIS J M, HAY A C, HOWARTH G J. Ontogeny of in situ behaviours relevant to dispersal and population connectivity in larvae of coral-reef fishes[J]. Mar Ecol Prog Ser, 2009, 379: 163-179. doi: 10.3354/meps07904
[7] LIU M J, GAO J, GUO H Y, et al. Transcriptomics reveal the effects of breeding temperature on growth and metabolism in the early developmental stage of Platax teira[J]. Biology, 2023, 12(9): 1161. doi: 10.3390/biology12091161
[8] 陈松林, 徐文腾, 刘洋. 鱼类基因组研究十年回顾与展望[J]. 水产学报, 2019, 43(1): 1-14. [9] APARICIO S, CHAPMAN J, STUPKA E, et al. Whole-genomeshotgun assembly and analysis of the genome of Fugu rubripes[J]. Science, 2002, 297(5585): 1301-1310.
[10] LIU D, WANG X Y, GUO H Y, et al. Chromosome-level genome assembly of the endangered humphead wrasse Cheilinus undulatus: insight into the expansion of opsin genes in fishes[J]. Mol Ecol Resour, 2021, 21(7): 2388-2406.
[11] ZHENG S Q, SHAO F, TAO W J, et al. Chromosome-level assembly of southern catfish (Silurus meridionalis) provides insights into visual adaptation to nocturnal and benthic lifestyles[J]. Mol Ecol Resour, 2021, 21(5): 1575-1592. doi: 10.1111/1755-0998.13338
[12] 廖静. 人工养殖尖翅燕鱼性价比高[J]. 海洋与渔业, 2018(11): 62-63. [13] LIU B, GUO H Y, ZHU K C, et al. Nutritional compositions in different parts of muscle in the longfin batfish, Platax teira (Forsskål, 1775)[J]. J Appl Anim Res, 2019, 47(1): 403-407. doi: 10.1080/09712119.2019.1649680
[14] CHEN S F, ZHOU Y Q, CHEN Y R, et al. fastp: an ultra-fast all-in-one FASTQ preprocessor[J]. Bioinformatics, 2018, 34(17): i884-i890. doi: 10.1093/bioinformatics/bty560
[15] ROBERTS R J, CARNEIRO M O, SCHATZ M C. The advantages of SMRT sequencing[J]. Genome Biol, 2013, 14(7): 405. doi: 10.1186/gb-2013-14-6-405
[16] HU J, WANG Z, SUN Z Y, et al. NextDenovo: an efficient error correction and accurate assembly tool for noisy long reads[J]. Genome Biol, 2024, 25(1): 107. doi: 10.1186/s13059-024-03252-4
[17] HU J, FAN J P, SUN Z Y, et al. NextPolish: a fast and efficient genome polishing tool for long-read assembly[J]. Bioinformatics, 2020, 36(7): 2253-2255. doi: 10.1093/bioinformatics/btz891
[18] SEPPEY M, MANNI M, ZDOBNOV E M. BUSCO: assessing genome assembly and annotation completeness[J]. Methods Mol Biol, 2019, 1962: 227-245.
[19] FLYNN J M, HUBLEY R, GOUBERT C, et al. RepeatModeler2 for automated genomic discovery of transposable element families[J]. PNAS, 2020, 117(17): 9451-9457. doi: 10.1073/pnas.1921046117
[20] HAAS B J, SALZBERG S L, ZHU W, et al. Automated eukaryotic gene structure annotation using EVidenceModeler and the Program to Assemble Spliced Alignments[J]. Genome Biol, 2008, 9(1): R7. doi: 10.1186/gb-2008-9-1-r7
[21] BIRNEY E, CLAMP M, DURBIN R. GeneWise and Genomewise[J]. Genome Res, 2004, 14(5): 988-995. doi: 10.1101/gr.1865504
[22] KIM D, LANGMEAD B, SALZBERG S L. HISAT: a fast spliced aligner with low memory requirements[J]. Nat Methods, 2015, 12(4): 357-360. doi: 10.1038/nmeth.3317
[23] PERTEA M, PERTEA G M, ANTONESCU C M, et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads[J]. Nat Biotechnol, 2015, 33(3): 290-295. doi: 10.1038/nbt.3122
[24] STANKE M, DIEKHANS M, BAERTSCH R, et al. Using native and syntenically mapped cDNA alignments to improve de novo gene finding[J]. Bioinformatics, 2008, 24(5): 637-644. doi: 10.1093/bioinformatics/btn013
[25] EMMS D M, KELLY S. OrthoFinder: phylogenetic orthology inference for comparative genomics[J]. Genome Biol, 2019, 20(1): 238. doi: 10.1186/s13059-019-1832-y
[26] KATOH K, STANDLEY D M. MAFFT Multiple Sequence Alignment Software Version 7: improvements in performance and usability[J]. Mol Biol Evol, 2013, 30(4): 772-780. doi: 10.1093/molbev/mst010
[27] KUMAR S, SULESKI M, CRAIG J M, et al. TimeTree 5: an expanded resource for species divergence times[J]. Mol Biol Evol, 2022, 39(8): msac174. doi: 10.1093/molbev/msac174
[28] YANG Z H. PAML 4: phylogenetic analysis by maximum likelihood[J]. Mol Biol Evol, 2007, 24(8): 1586-1591. doi: 10.1093/molbev/msm088
[29] MENDES F K, VANDERPOOL D, FULTON B, et al. CAFE 5 models variation in evolutionary rates among gene families[J]. Bioinformatics, 2021, 36(22/23): 5516-5518.
[30] EDDY S R. Accelerated profile HMM searches[J]. PLoS Comput Biol, 2011, 7(10): e1002195. doi: 10.1371/journal.pcbi.1002195
[31] YU G C, WANG L G, HAN Y Y, et al. clusterProfiler: an R package for comparing biological themes among gene clusters[J]. OMICS, 2012, 16(5): 284-287. doi: 10.1089/omi.2011.0118
[32] 高杰, 郭华阳, 刘明鉴, 等. 尖翅燕鱼染色体核型分析[J]. 海洋渔业, 2022, 44(5): 535-542. [33] OSHIUMI H , TSUJITA T , SHIDA K, et al. Prediction of the prototype of the human Toll-like receptor gene family from the pufferfish, Fugu rubripes, genome[J]. Immunogenetics, 2003, 54: 791-800.
[34] HU Y C, TAN R H, ZHU X, et al. Genome-wide identification, phylogeny and expressional profile of the Dmrt gene family in Chinese sturgeon (Acipenser sinensis)[J]. Sci Rep, 2024, 14(1): 4231. doi: 10.1038/s41598-024-54899-9
[35] ZHU K C, ZHANG N, LIU B S, et al. A chromosome-level genome assembly of the yellowfin seabream (Acanthopagrus latus; Hottuyn, 1782) provides insights into its osmoregulation and sex reversal[J]. Genomics, 2021, 113(4): 1617-1627. doi: 10.1016/j.ygeno.2021.04.017
[36] ZHANG D C, GUO L, GUO H Y, et al. Chromosome-level genome assembly of golden pompano (Trachinotus ovatus) in the family Carangidae[J]. Sci Data, 2019, 6(1): 216. doi: 10.1038/s41597-019-0238-8
[37] LIANG Y, XIAN L, PAN J M, et al. De Novo genome assembly of the whitespot parrotfish (Scarus forsteni): a valuable scaridae genomic resource[J]. Genes (Basel), 2024, 15(2): 249. doi: 10.3390/genes15020249
[38] ZHOU Q, GUO X Y, HUANG Y, et al. De novo sequencing and chromosomal-scale genome assembly of leopard coral grouper, Plectropomus leopardus[J]. Mol Ecol Resour, 2020, 20(5): 1403-1413. doi: 10.1111/1755-0998.13207
[39] CHEN X H, ZHONG L Q, BIAN C, et al. High-quality genome assembly of channel catfish, Ictalurus punctatus[J]. GigaScience, 2016, 5(1): 39. doi: 10.1186/s13742-016-0142-5
[40] LI J, BIAN C, HU Y C, et al. A chromosome-level genome assembly of the Asian arowana, Scleropages formosus[J]. Sci Data, 2016, 3: 160105.
[41] LI S S, XIE Z Z, CHEN P, et al. The complete mitochondrial genome of the Platax teira (Osteichthyes: Ephippidae)[J]. Mitochondrial DNA A DNA Mapp Seq Anal, 2016, 27(2): 796-797.
[42] HUGHES L C, ORTÍ G, HUANG Y, et al. Comprehensive phylogeny of ray-finned fishes (Actinopterygii) based on transcriptomic and genomic data[J]. PNAS, 2018, 115(24): 6249-6254. doi: 10.1073/pnas.1719358115
[43] HE S, LI L, LYU L Y, et al. Mandarin fish (Sinipercidae) genomes provide insights into innate predatory feeding[J]. Commun Biol, 2020, 3(1): 361. doi: 10.1038/s42003-020-1094-y
[44] MARSHALL H D, COULSON M W, CARR S M. Near neutrality, rate heterogeneity, and linkage govern mitochondrial genome evolution in Atlantic cod (Gadus morhua) and other gadine fish[J]. Mol Biol Evol, 2009, 26(3): 579-589.
[45] TEREKHANOVA N V, LOGACHEVA M D, PENIN A A, et al. Fast evolution from precast bricks: genomics of young freshwater populations of threespine stickleback Gasterosteus aculeatus[J]. PLoS Genet, 2014, 10(10): e1004696. doi: 10.1371/journal.pgen.1004696
[46] AO J Q, MU Y N, XIANG L X, et al. Genome sequencing of the perciform fish Larimichthys crocea provides insights into molecular and genetic mechanisms of stress adaptation[J]. PLoS Genet, 2015, 11(4): e1005118. doi: 10.1371/journal.pgen.1005118
[47] NIEDERRITER A R, DAVIS E E, GOLZIO C, et al. In vivo modeling of the morbid human genome using Danio rerio[J]. J Vis Exp, 2013(78): e50338.
[48] DAVIDSON W S, KOOP B F, JONES S J M, et al. Sequencing the genome of the Atlantic salmon (Salmo salar)[J]. Genome Biol, 2010, 11: 403. doi: 10.1186/gb-2010-11-9-403
[49] LEVANTI M, RANDAZZO B, VIÑA E, et al. Acid-sensing ion channels and transient-receptor potential ion channels in zebrafish taste buds[J]. Ann Anat, 2016, 207: 32-37. doi: 10.1016/j.aanat.2016.06.006
[50] MOHAMED N A, SAAD M F, SHUKRY M, et al. Physiological and ion changes of Nile tilapia (Oreochromis niloticus) under the effect of salinity stress[J]. Aquac Rep, 2021, 19: 100567. doi: 10.1016/j.aqrep.2020.100567
[51] IP Y K, CHEW S F. Ammonia production, excretion, toxicity, and defense in fish: a review[J]. Front Physiol, 2010, 1: 134.
[52] RANDALL D J, TSUI T K N. Ammonia toxicity in fish[J]. Mar Pollut Bull, 2002, 45(1/2/3/4/5/6/7/8/9/10/11/12): 17-23.
[53] ARILLO A, MARGIOCCO C, MELODIA F, et al. Ammonia toxicity mechanism in fish: studies on rainbow trout (Salmo gairdneri Rich.)[J]. Ecotoxicol Environ Saf, 1981, 5(3): 316-328. doi: 10.1016/0147-6513(81)90006-3