扁吻鱼胚胎发育及仔稚鱼形态与生长观察研究

韩军军, 贺疆滔, 陈朋, 胡江伟, 祁峰, 封永辉, 蔡林钢, 时春明, 张人铭

韩军军, 贺疆滔, 陈朋, 胡江伟, 祁峰, 封永辉, 蔡林钢, 时春明, 张人铭. 扁吻鱼胚胎发育及仔稚鱼形态与生长观察研究[J]. 南方水产科学, 2021, 17(1): 59-66. DOI: 10.12131/20200177
引用本文: 韩军军, 贺疆滔, 陈朋, 胡江伟, 祁峰, 封永辉, 蔡林钢, 时春明, 张人铭. 扁吻鱼胚胎发育及仔稚鱼形态与生长观察研究[J]. 南方水产科学, 2021, 17(1): 59-66. DOI: 10.12131/20200177
HAN Junjun, HE Jiangtao, CHEN Peng, HU Jiangwei, QI Feng, FENG Yonghui, CAI Lingang, SHI Chunming, ZHANG Renming. Observation on embryonic development, morphology and growth of larvae and juveniles of Aspiorhynchus laticeps[J]. South China Fisheries Science, 2021, 17(1): 59-66. DOI: 10.12131/20200177
Citation: HAN Junjun, HE Jiangtao, CHEN Peng, HU Jiangwei, QI Feng, FENG Yonghui, CAI Lingang, SHI Chunming, ZHANG Renming. Observation on embryonic development, morphology and growth of larvae and juveniles of Aspiorhynchus laticeps[J]. South China Fisheries Science, 2021, 17(1): 59-66. DOI: 10.12131/20200177

扁吻鱼胚胎发育及仔稚鱼形态与生长观察研究

基金项目: 农业农村部财政专项“西北地区重点水域渔业资源与环境调查”
详细信息
    作者简介:

    韩军军(1990—),男,硕士,工程师,从事水生生物繁育和生态修复研究。E-mail: hanjunjun1990@163.com

    通讯作者:

    张人铭 (1965—),男,研究员,从事渔业资源及鱼类人工繁殖技术研究。E-mail: xj01zrm@163.com

  • 中图分类号: S 961.2+3

Observation on embryonic development, morphology and growth of larvae and juveniles of Aspiorhynchus laticeps

  • 摘要:

    文章采用干法受精获取受精卵的方法,对扁吻鱼(Aspiorhynchus laticeps)胚胎发育和仔稚鱼的形态变化进行观察。研究结果表明,扁吻鱼卵为沉性卵,卵细胞呈青灰色、微黏性。该次孵化共获取成熟卵21.45×104粒,受精率88.95%,孵化率93.27%。扁吻鱼成熟卵卵径 (2.46±0.13) mm,受精卵吸水后膨胀,卵径 (3.42±0.15) mm。根据胚胎发育的外部形态和典型特征将其分为受精卵、胚盘形成、卵裂、囊胚、原肠、神经胚、器官分化和孵化共8个阶段,整个胚胎发育耗时115.72 h,所需积温3 765.42 ℃·h。初孵仔鱼全长 (7.58±0.31) mm,从初孵到卵黄吸收完全历时9 d,全长特定生长率为 (3.63±0.68)%。仔鱼后期历时36 d,全长特定生长率为 (3.51±2.10)%。稚鱼期历时63 d,全长特定生长率为 (1.32±0.17)%,与仔鱼前期和仔鱼后期差异显著 (P<0.05)。1年后幼鱼全长 (15.71±0.89) cm,体质量 (28.29±6.15) g,体色和体型接近成鱼。

    Abstract:

    In this study, the embryonic development and morphological changes of larvae and juveniles of Aspiorhynchus laticeps were observed by dry fertilization. The results show that the eggs were sinking eggs, and the color of egg cells was bluish gray and slightly sticky. A total of 21.45×104 mature eggs were obtained and the average of fertilization rate and hatching rate were 88.95% and 93.27%, respectively. The diameter of mature eggs was (2.46±0.13) mm, and that of fertilized eggs was (3.42±0.15) mm. According to the external morphology and typical characteristics of embryo development, it can be divided into eight stages: fertilization, cleavage, blastocyst, gastrula, neuroembryo, organ differentiation and hatching. It took 115.5 h for the whole embryo to develop, and the total accumulative temperature was 3 765.42 ℃·h. The total length of newly hatched larvae was (7.58±0.31) mm. The development of prelarva from hatching to yolk absorption took approximately 9 d, and the specific growth rate (SGR) of total length was (3.63±0.68)% and that of total length was (3.51±2.10)%. The juvenile period lasted 63 d and the SGR of total length was (1.32±0.17)%, significantly different from those of early and late larval stages (P<0.05). The total length and body mass of A. laticeps young fish was (15.71±0.89) cm and (28.29±6.15) g, respectively, and the body color and body shape were similar to adult fish.

  • 南海又称南中国海或中国南海,是世界第三大边缘海[1],拥有丰富的海洋生物资源。同时,南海也是海上“丝绸之路”的起点和国家海洋强国战略的“主战场”,具有重要的生态、经济和政治意义。但20世纪80年代以来,南海渔业资源受海洋污染和过度捕捞的双重影响而不断衰退[2]。因此自1999年起,我国政府开始对南海区域渔船数量等进行管控,调整渔业生产结构,以降低南海渔业资源的开发强度[3],保护海洋渔业资源。

    鱼类是海洋中最主要的动物类群,现存种类主要包括硬骨鱼纲、软骨鱼纲、盲鳗纲和头甲纲4个纲[4]。张俊等[5]发现南海中南部中层鱼中硬骨鱼纲的种类最多,是南海海洋生态系统的主要物种。鳃是鱼类等海洋生物的主要呼吸器官,除了具有呼吸作用,还具有滤食、渗透压调节、pH平衡调节和氨排泄等功能[6]。同时,鳃也是海洋环境微生物进入宿主体内的第一道关卡,对宿主内生的微生物群落具有重要调控作用。海洋鱼类鳃组织中的微生物类群是一个重要的免疫屏障,也是鱼类抵御海洋环境中潜在病原体的第一道防线。研究发现,鱼类鳃组织中的细菌与宿主的疾病具有重要的相关性,Toenshoff等[7]通过对大西洋鲑 (Salmon salar) 鳃组织菌群研究发现了一种新型的β-proteobacterial会引起宿主鱼体上皮囊炎等疾病。因此,对鳃组织上附着的微生物群落进行研究,对更好地理解海洋动物在海洋环境中的适应性作用、以及病原入侵与鱼类病害的发生机制具有重要意义。近年来,海洋鱼类微生物研究大多关注在各种鱼类的肠道组织,如日本花鲈 (Lateolabrax japonicus)[8]、大西洋鳕 (Gadus morhua)[9]、大西洋鲑[10]、大菱鲆 (Scophthalmus maximus)[11]和大黄鱼 (Larimichthys crocea)[12]等经济性物种,而围绕鳃组织微生物多样性及其功能的研究匮乏。王小彤和许强华[13]报道了4种南极鱼类鳃组织微生物群落多样性,Pratte等[14]对珊瑚礁不同鱼类鳃组织微生物群落多样性进行了报道。但相较于种类繁多的海洋鱼类资源,鳃组织微生物的研究数据还十分匮乏。

    16S rDNA (或16S rRNA) 基因作为细菌的保守特征序列是进行细菌系统发育和分类鉴定的标签基因[15]。随着高通量测序成本的降低和扩增子 (Amplicon) 技术的日臻完善,16S rDNA扩增子已逐渐成为微生物多样性研究的主流方法[16]。本研究基于16S扩增子技术,对从南海近岸海域7个站点采集的45种海洋硬骨鱼纲生物的鳃组织样品进行了细菌群落分析,以期为了解影响海洋动物鳃组织微生物群落的主要因素、探索海洋动物与海洋环境微生物的相互依存关系提供参考。

    2019年3月,在7个站点共采集了52份鱼类鳃组织样品 (分类见表1),分属于11目33科40属的45种硬骨鱼纲鱼类。7个站点的位置见图1,其中D3、D8、F3、F8、H3、H8和I9站点距离海岸的距离分别为60、215、10、175、20、130和165 km。采集到的鲜活样品迅速装于无菌密封袋中,于−80 ℃暂存。上岸后用无菌解剖工具完整取下鳃组织并于液氮中速冻后,置于−80 ℃长期保存。本实验将同一站点采集的同种鱼类的鳃组织,根据鱼体大小,进行1~3尾的混合而作为一个样品,以保证文库构建和测序数据的质量。

    表  1  样品分类信息
    Table  1  Taxonomy information of samples
    站点
    Sampling station
    样品编号
    Sample ID

    Order
    物种
    Species
    D3D3.1S 鲈形目 条纹䱨 Terapon theraps
    D3D3.2S 鲈形目 日本金线鱼 Nemipterus japonicus
    D3D3C.2S 鲈形目 列牙䱨 Pelates quadrilineatus
    F3F3.2S 鲈形目 鹿斑鲾 Leiognathus ruconius
    F3F3.4S 鲽形目 少牙斑鲆 Pseudorhombus oligodon
    F3F3.5S 鲈形目 带鱼 Trichiurus lepturus
    D8D8.1S 鲈形目 蓝圆鲹 Decapterus maruadsi
    D8D8.3S 鳗鲡目 网纹裸胸鳝 Gymnothorax reticularis
    D8D8.5S 鮟目 黑鮟 Lophiomus setigerus
    D8D8.6S 鲉形目 环纹蓑鲥 Pterois lunulata
    D8D8.8S 鲈形目 深水金线鱼 Nemipterus bathybius
    D8D8.12S 鲉形目 深海红娘鱼 Lepidotrigla abyssalis
    F8F8.1S 鲉形目 环纹蓑鲥 Pterois lunulata
    F8F8.4S 鲉形目 单棘豹鲂鮄 Daicocus peterseni
    F8F8.8S 鲈形目 二长棘鲷 Parargyrops edita
    F8F8.11S 鮟目 棘茄鱼 Halieutaea stellata
    F8F8.16S 鲉形目 日本红娘鱼 Lepidotrigla japonice
    F8F8.17S 鲉形目 拟蓑鲉 Parapterois heterurus
    F8F8.18S 海龙鱼目 鳞烟管鱼 Fistularia petimba
    H3H3.3S 鲈形目 翼红娘鱼 Lepidotrigla alata
    H3H3.4S 鲈形目 中华䲢 Uranoscopus chinensis
    H3H3.5S 鳗鲡目 大鳞鳞头鲉 Sebastapistes megalepis
    H3H3.7S 鮟目 无备虎鲉 Minous inermis
    H8H8.1S 鲉形目 大鳞短额鲆 Engyprosopon grandisquama
    H8H8.2S 鲈形目 鳞烟管鱼 Fistularia petimba
    H8H8.3S 鲉形目 繁星鲆 Bothus myriaster
    H8H8.7S 鲉形目 孔鰕虎鱼 Trypauchen vagina
    H8H8.10S 鲉形目 斑鰶 Konosirus punctatus
    H8H8.12S 鲈形目 皮氏叫姑鱼 Johnius belangerii
    H8H8.15S 鮟目 短吻鲾 Leiognathus brevirostris
    H8H8.24S 鲉形目 大头狗母鱼 Trachinocephalus myops
    H8H8.25S 鲉形目 大头狗母鱼 Trachinocephalus myops
    I9I9.1S 海龙鱼目 大鳞鳞头鲉 Sebastapistes megalepis
    I9I9.5S 鲉形目 黄纹拟鲈 Parapercis xanthozona
    I9I9.7S 鲈形目 寿鱼 Banjos banjos
    I9I9.8S 鲉形目 棕斑宽吻鲀 Amblyrhynchotes rufopunctatus
    I9I9.9S 鲉形目 单棘豹鲂鮄 Daicocus peterseni
    I9I9.10S 鲽形目 松球鱼 Monocentris japonica
    I9I9.11S 海龙鱼目 新平鲉 Neosebastes entaxis
    I9I9.12S 鲽形目 二长棘鲷 Parargyrops edita
    I9I9.13S 鲈形目 横带眶棘鲈 Scolopsis inermis
    I9I9.14S 鲱形目 盔蓑鲉 Ebosia bleekeri
    I9I9.15S 鲈形目 冠鲽 Samaris cristatus
    I9I9.16S 鲈形目 叉斑狗母鱼 Synodus macrops
    I9I9.17S 海龙鱼目 美尾䲗 Calliurichthys japonicus
    I9I9.18S 鲉形目 瑞氏红鲂鮄 Satyrichthys rieffeli
    I9I9.20S 鲉形目 虻鲉 Erisphex pottii
    I9I9.21S 鲉形目 红鲬 Bembras japonicus
    I9I9.22S 鲈形目 六带拟鲈 Parapercis sexfasciata
    I9I9.23S 鲈形目 黄纹拟鲈 Parapercis xanthozona
    I9I9.24S 鲈形目 尖牙鲈 Synagrops japonicus
    I9I9.25S 鲈形目 水珍鱼 Argentina kagoshimae
    下载: 导出CSV 
    | 显示表格
    图  1  广东近岸南海海域采样站点示意图
    Fig. 1  Diagram of sampling stations in Guangdong coastal waters in South China Sea

    分别按照采样站位、宿主来源可将52个样品进行两种分组编排。依站点不同可分为7组,依宿主目水平可分为3组,其中鲈形目样品最多 (19个),鲉形目16个,其他小目统一编为一组 (17个)。

    使用HiPure Bacterial DNA kit (Magen,广州) 试剂盒,参照说明书进行总DNA提取,利用Thermo NanoDrop One (Thermo Fisher Scientific Inc., 上海) 检测DNA的纯度和浓度。针对16S rDNA基因的V4区进行巢套式扩增,使用带barcode的特异引物及TaKaRa Premix Taq® Version 2.0 (TaKaRa Biotechnology Co., 大连) 进行PCR扩增。所用一轮引物为27F (5'-AGAGTTTGATCCTGGCTCAG-3') 和1492R (5'-TACGGYTACCTTGTTAYGACTT-3'),二轮引物为515F (5'-GTGCCAGCMGCCGCGGTAA-3') 和806R (5'-GGACTACHVGGGTWTCTAAT-3'),扩增区域长度为310 bp。一轮PCR反应体系为 (20 μL):Premix TaqTM (Ex TaqTM Version 2.0 plus dye)/Mix 10 μL,上下游引物 (10 μmol·L−1) 各0.4 μL,模板DNA (40~60 ng) 各1 μL,Nuclease-free water 8 μL。二轮PCR反应体系 (50 μL) 为:Premix TaqTM (Ex TaqTM Version 2.0 plus dye)/Mix 25 μL,上下游引物 (10 μmol·L−1) 各1 μL,模板DNA (40~60 ng) 1 μL,Nuclease-free water 22 μL。PCR反应条件为94 ℃ 5 min;94 ℃ 30 s,53 ℃ 30 s,72 ℃ 30 s,30个循环;72 ℃ 8 min。DNA提取、文库构建及测序工作委托广东美格基因科技有限公司进行。

    采用Illumina Nova 6000平台对构建的扩增子文库进行PE250测序 (广东美格基因科技有限公司, 广州),测序完成后得到原始的下机数据:双端原始序列 (Paired-end Raw Reads),利用fastp[17] (An ultra-fast all-in-one FASTQ preprocessor, V 0.14.1,https://github.com/OpenGene/fastp) 分别对两端的原始序列进行滑窗 (sliding window) 质量剪裁,同时根据序列首尾两端的引物信息利用cutadapt软件 (https://github.com/marcelm/cutadapt/) 去除引物,得到质控和过滤后高质量的双端优化序列 (Clean Reads)。对于双端优化序列,根据reads之间的overlap关系,利用usearch-fastq_mergepairs (V10, http://www.drive5.com/usearch/预设参数包含最小overlap长度设置为16 bp,拼接序列的overlap区允许的最大错配5 bp等) 过滤不符合的Tags,获得原始的拼接序列 (Raw Tags),利用fastp对拼接序列数据进行滑窗质量剪裁,得到优化的拼接序列 (Clean Tags)。利用UPARSE[18]方法进行OTU聚类。利用usearch-sintax将每个OTU的代表序列与SILVA (16S)、RDP (16S)、Greengenes (16S) 数据库进行比对获得物种注释信息,以达到了解所有序列物种来源的目的,最终用SILVA (16S) 数据库注释得到的OTU结果进行后续分析。

    使用过滤后获得的最终OTU丰度表进一步进行物种群落组成分析、Alpha多样性分析、Beta多样性分析。Alpha多样性使用物种丰富度或均匀度来表征,其中Chao1[19]、 ACE和Observed species指数用于估算群落物种丰富度[20]。而Shannon[21]和Simpson指数用于反映群落物种多样性。Beta多样性通过非度量多维尺度分析 (Nonmetric Multidimensional Scaling, NMDS) 来表征,NMDS 可以基于进化关系或数量距离矩阵对比样本之间的差异。

    使用PICRUSt软件[22]预测OTU的功能,首先消除16S marker gene在物种基因组中的拷贝数的影响,然后通过每个OTU所对应的greengene id,分别比对到京都基因和基因组百科全书 (Kyoto Encyclopedia of Genes and Genomes, KEGG) 数据库和COG (Clusters of Orthologous Groups) 数据库,获得功能预测信息,根据OTU丰度计算各功能类别的丰度,分别得到各样本在不同分类水平的功能丰度表。

    使用在线微生物组分析平台MicrobiomeAnalyst (https://www.microbiomeanalyst.ca) 进行微生物数据的多样性和功能分析[23]

    对52份16S rDNA文库进行测序,共获得有效拼接序列2 952 366条,平均每个文库56 776条;共产生10 617个OTU。图2表明各样本稀释曲线随着测序深度的增加逐渐趋于平缓,所有样品均达到平台期,表明本研究的测序深度足够。由稀释曲线的位置可知,I9站点样品的OTU数量相对较高,D8站点样品较低。

    图  2  各个样本的稀释曲线
    Fig. 2  Rarefaction curve of each sample

    对不同站点样品的细菌门水平分类组成和相对丰度情况进行了统计 (图3)。在总体水平上,平均相对丰度≥1%的优势菌门依次为变形菌门 (Proteobacteria, 71.3%)、厚壁菌门 (Firmicutes, 10.3%)、浮霉菌门 (Planctomycetes, 8.3%)、放线菌门 (Actinobacteria, 3.4%)、拟杆菌门 (Bacteroidetes, 2.1%) 和蓝细菌 (Cyanobacteria, 1%)。不同站点的优势菌门丰度不同,但变形菌门均为所有站点中丰度最高的优势类群,其中以F3和F8站点最高,D3站点最低。厚壁菌门主要分布于D3、D8、I9站点,而浮霉菌门主要分布于D3、H8、I9站点,其他门在不同站点间丰度比例略有差异。另外,F3和F8站点的细菌种类较少,而H8和I9站点的细菌种类相对较多。

    图  3  鱼鳃细菌门水平物种组成相对丰度
    Fig. 3  Relative abundance of taxonomy composition of bacteria in fish gill at phylum level

    进一步在属水平上对各站点样品的分类组成和相对丰度进行统计分析。在所有站点中,相对丰度≥2%的优势菌属依次为不动杆菌属 (Acinetobacter, 17.2%)、青枯菌属 (Ralstonia, 7.5%)、嗜冷杆菌属 (Psychrobacter, 7.5%)、鞘氨醇单胞菌 (Sphingomonas, 3.5%)、BD1_7_clade (3.3%) 和Clostridium_sensu_stricto_1 (3%,图4)。不同菌属在不同站点的丰度占比不同,但不动杆菌属在所有站点中均为优势菌属。青枯菌属 (26.5%) 和鞘氨醇单胞菌属 (20.5%) 主要分布于F3站点,BD1_7_clade在F8站点占比最高 (20.3%),而嗜冷杆菌属在H3 (20.7%) 和I9站点 (17.9%) 丰度较高。总体来看,D8、H8、I9站点的细菌在属水平上的种类相对其他几个站点较多,各个站点优势菌属的相对丰度相差较大。

    图  4  鱼鳃细菌属水平物种组成相对丰度
    Fig. 4  Relative abundance of taxonomy composition of bacteria in fish gill at genus level

    Alpha多样性指数可用于反映样品内部微生物群落的相对丰度及多样性。本研究计算了各站点分组样品的Alpha多样性指数 (图5)。结果表明,H3、H8和I9站点的物种丰富度,包括Chao1 (图5-a)、ACE (图5-b) 和Observed species指数 (图5-c) 均极显著高于D3、D8、F3和F8站点 (P<0.01),尤以H8和I9最高。各站点的群落多样性,包括Simpson (图5-d) 和Shannon指数 (图5-e),也表现出了显著差异 (P<0.05),其中D3站点最高,F3站点最低。然而,基于宿主目水平的分组结果显示,不同物种分组间的Alpha指数均无显著差异 (P>0.05,图5-f,以Chao1指数为例)。上述结果表明,不同站点细菌群落多样性和丰富度指数均存在显著差异,而基于宿主目水平的分组之间未表现出显著差异。

    图  5  Alpha多样性指数箱线图
    注:a—e. 基于站点分组的Chao1、ACE、Observed、Simpson、Shannon 多样性指数箱线图,不同颜色代表不同站点的样本;f. 基于样本宿主(目水平) 分类分组的Chao1 指数箱线图。
    Fig. 5  Boxplot of Alpha diversity index
    Note: a–e. Boxplot of diversity indexes of Chao1, ACE, Observed, Simpson and Shannon based on site grouping, and different colors represent samples of different sites; f. Boxplot of Chao1 index based on sample host (Order level) classification.

    Beta多样性可用于描述不同样本微生物群落之间物种组成的差异。我们基于不同分组对各个样本进行NMDS分析 (图6)。按照站点分组结果显示,各个站点样本的Beta多样性具有极显著差异 (P<0.001),且不同站点样本分别聚集在一起 (图6-a)。总体来看,与α多样性结果一致,即H3、H8、I9站点与D3、D8、F3、F8站点明显分为两个主要类群。其中D3和F3、D8 和F8 站点的样本距离更近,H3、H8和I9站点的10个样本聚集在一起,I9站点一部分样本单独分群(图6-a:I9-2)。然而,基于样本目水平的分组未发现明显分群,不同的目分组间不存在显著差异 (P>0.05,图6-b)。由此可见,地理站位对硬骨鱼鳃组织内微生物可能具有重要影响,而样本宿主遗传因素可能不是影响鱼鳃细菌群落分布的主要因素。

    图  6  基于不同分组方式的NMDS分析 [非参数多元方差分析 (PerMANOVA)]
    注:a. 基于站点分组的 NMDS 分析; b. 基于样本的目水平分组。图中两个样本点越接近,表示两个样本物种组成越相似。横纵坐标表示样本间的相对距离,无实际意义。图形中的每一个点代表一个样本,不同颜色代表不同的样本分组信息。NMDS 结果的优劣用胁迫系数 (stress) 来衡量,此值越小越好,当小于 0.2 时表示可以用 NMDS 的二维点图表示组间或组内差异。
    Fig. 6  NMDS analysis based on different grouping modes [Nonparametric multivariate analysis of variance (PerMANOVA)]
    Note: a. NMDS analysis based on site grouping; b. Grouping based on the order level of the sample. The closer the two sample points are,the more similar the composition of the two sample species are. The abscissa and ordinate indicate the relative distance between samples, which is of no practical significance. Each point in the graph represents a sample, and different colors represent different sample grouping information. The advantages and disadvantages of testing NMDS results are measured by stress coefficients, the smaller, the better, and the value less than 0.2 indicates that the two-dimensional dot plot of NMDS can be used to represent differences between or within groups.

    基于KEGG比对共获得8个一级层级和41个二级层级功能分类。各样本在一级功能层级组成上的差异不明显 (数据未展示),占比最高的均为基因与新陈代谢相关功能 (Metabolism, 50%),其次是遗传信息处理 (Genetic information processing, 16%)、环境信息处理 (Environmental information processing, 13%)。图7为KEGG数据库比对到的22个相对丰度较高 (>1%) 的二级功能层级热图。总体来看,二级功能层级中膜转运 (Membrane transport)功能丰度最高 (11.02%),其次为氨基酸代谢 (Amino acid metabolism, 10.5%)。但不同功能在不同样本中丰度差异明显,膜转运 (14%)、碳水化合物代谢 (Carbohydrate metabolism, 10%)、氨基酸代谢 (9%)、复制和修复 (Replication and repair, 9%) 在I9.11S中丰度最高,在D8.6S中次之。样本聚类显示鱼鳃中细菌功能基因分布及丰度与站点没有明显的关联性。另外,基于COG数据库预测结果与KEGG的一级功能结果类似 (数据未展示),各样本中的功能基因比例总体差别不大。

    图  7  22个KEGG-L2功能层级相对丰度热图
    注:横坐标代表不同的样本,纵坐标代表不同的功能层级,图片上方色条表示不同分组,并按照样本聚类。
    Fig. 7  Heatmap of relative abundance of 22 KEGG-L2 functional levels
    Note: The abscissa represents different samples; the ordinate represents different functional levels; the color bars indicate different groups, which are clustered according to samples.

    两种功能数据库预测结果均表明,鱼鳃细菌群落可能主要参与了宿主营养物质转运及代谢等相关进程。

    本研究发现,硬骨鱼鳃组织细菌群落中的优势菌门包括变形菌门、厚壁菌门、浮霉菌门、放线菌门、拟杆菌门和蓝细菌,这与其他学者对海洋环境菌群统计分析所得结果一致,如穆大帅等[24]对我国分离鉴定的577种海洋细菌进行多样性分析,发现在门分类水平上分布最多的是变形菌门,其次是拟杆菌门、放线菌门和厚壁菌门。白洁等[25]采用高通量测序方法分析了南海南部海域浮游细菌丰度,发现南海南部海域浮游细菌的优势类群为变形菌门、蓝藻门和拟杆菌门,优势亚群为γ-变形菌纲、α-变形菌纲、蓝藻菌纲和黄杆菌纲。徐新亚等[26]报道广西北部湾分离鉴定到1 843种海洋细菌,分别隶属于放线菌门、拟杆菌门、厚壁菌门、浮霉菌门和变形菌门。本研究结果表明海洋环境与鱼类鳃中的微生物群落至少在门水平上具有较高的一致性。变形菌门作为海洋中主要的微生物类群之一,在海洋生态系统的碳循环、氮循环等多种物质循环过程中均发挥重要作用[27],有研究证实海洋表面的变形菌门含有丰富的固氮种群[28]。厚壁菌门最主要的功能是降解纤维类物质,拟杆菌门主要降解非纤维物质[29-30],而上述3类细菌在鱼鳃组织上的富集可能与宿主在海洋环境中的代谢需求有关。此外,Legrand等[31]认为,黄尾鰤 (Seriola lalandi) 鳃中丰度较高的浮霉菌门可能在海洋鱼类体内固氮作用中发挥功能,有助于促进氮元素在海洋环境与海洋生物之间的垂直传播。由于本研究是通过16S的物种组成反推预测群落功能,并未直接对宏基因组和功能基因进行测序。因此,功能分析结果会与真实情况有一定差距。未来研究还需进一步结合宏基因组测序及海洋环境有机物质的测定,才能更好地解读鳃组织微生物在鱼体环境适应性方面发挥的功能。

    近年来,国内外学者在海洋细菌多样性研究领域贡献了大量成果,但多数只关注水体中游离的细菌,对于宿主,尤其是鳃组织上,共附生细菌群落潜在功能的研究相对较少。本研究预测到鱼类鳃中细菌多参与宿主的新陈代谢过程,积极参与氨基酸代谢反应和碳水化合物代谢反应,表明鱼鳃中细菌可能在碳循环中发挥作用。Kessel等[32]研究了鲤 (Cyprinus carpio) 和斑马鱼 (Danio rerio) 鳃中微生物在氮循环中的作用,结果发现鳃中含有一个活跃的N-循环细菌群,属于β-变形杆菌纲成员,可能对许多鱼类降低氨浓度、防止高浓度氨中毒有益。本研究发现所有样本中均包含丰度不一的β-变形菌纲成员,提示海水鱼类鳃组织中可能也具有类似功能的菌群。Pratte等[14]对于法国Moorea珊瑚礁采集的鱼鳃和肠道等样本的微生物群落进行研究,发现成鱼和幼鱼鳃中的微生物群落组成与肠道有显著差异。鳃中最丰富的细菌是希瓦氏菌属和Endozoicimonaceae成员,且微生物群落与环境样品共享的OTU百分比高达75%,表明鳃微生物组可能受到周围浮游微生物群落变化的影响。本研究中D8站点的网纹裸胸鳝 (Gymnothorax reticularis)、黑鮟 (Lophiomus setigerus)、环纹蓑鲥 (Pterois lunulata) 和深水金线鱼 (Nemipterus bathybius) 中的希瓦氏菌属细菌占比同样较高,与该结果吻合。然而,其他多数样品中则是变形菌门的不动杆菌属细菌比例较高,是所有站点的优势菌属。该结果说明不同站点或不同宿主物种的微生物组成可能差异很大。有研究报道,不动杆菌属细菌作为非发酵革兰阴性菌,具有代谢多样性,可有效降解海水中己内酰胺、石油油烃等海洋污染成分[33]。不动杆菌属在海水鱼鳃部的富集或可有效促进环境污染物等的代谢,从而对宿主产生潜在的保护作用。除了D3站点,其他6个站点的优势菌属还包括嗜冷杆菌属。该属细菌多分布于潮湿、含盐或寒冷的生境中,且能分泌适应低温的酶,可在海洋低温下发挥其生物活性[34],且具有较高的催化效率。然而,本研究发现的广泛存在于鱼鳃组织中的嗜冷杆菌是否也具有类似功能值得深入研究。

    本研究发现站点对鱼鳃细菌群落的影响远高于宿主物种遗传因素的影响。近年来,其他研究者也发现地理位置对海洋环境中的微生物群落产生影响。翟万营等[35]发现地理位置对南极中山站附近海水微生物群落组成有重要影响;Zhang等[36]认为除理化性质外,地理位置对南极中山站、长城站周边土壤样品和近海沉积物中的细菌群落组成产生一定影响。与土壤等环境样本不同,生物体内的微生物群落除受到环境影响外,还会受到饮食、遗传、发育阶段、健康等因素的影响。本研究发现海洋鱼类鳃中的菌群也受到地理位置的影响,可能由不同站点的水质环境、浮游微生物和营养物质与鳃的直接接触造成。海洋动物鳃组织与肠道等内生组织不同,鳃暴露于水体环境中,没有额外的免疫屏障,可能导致鱼鳃细菌群落受到环境的影响尤为显著。鉴于海洋鱼类具有自主、灵活的游泳能力,由此推测,其在海洋微生物传播和海洋微生物群落调节中可能发挥着重要作用。

    此外,本研究还发现I9站点的20个样本明显分为两个群,这可能有两个原因:1) 两个群落在宿主物种水平上存在差异,第一个群落包括大鳞鳞头鲉 (Sebastapistes megalepis)、黄纹拟鲈 (Parapercis xanthozona)、寿鱼 (Banjos banjos) 和棕斑宽吻鲀 (Amblyrhynchotes rufopunctatus) 等物种;第二个群落包括冠鲽 (Samaris cristatus)、叉斑狗母鱼 (Synodus macrops)、美尾䲗 (Calliurichthys japonicus)、虻鲉 (Erisphex pottii) 和红鲬 (Bembras japonicus) 等物种。通过调查采样记录发现,第一个类群部分物种鱼类营养级高于第二个类群,而第二个类群的样本鱼体较小,由此推测,营养级和个体大小 (发育阶段) 或许是造成该批次群落差异的一个原因。2) 扩增、建库、测序等试验环节批次的不同也有可能会对上述群落差异造成一定影响,具体原因还需要更多试验证实。

    综上所述,南海具有丰富的海洋微生物资源,本研究发现南海近岸硬骨鱼鳃组织中细菌群落的组成非常丰富。采样站点相比宿主来源对鳃组织上细菌群落结构及多样性具有更主要的影响,它们可能在辅助宿主营养物质转运及代谢方面发挥积极作用。然而本研究也存在很多不足如:样本和站点数量较少,未对海水样品进行对照研究等。未来还需进一步扩大样品的种类、数量,并结合水质环境、理化因子和宿主的生理生长指标及生活方式等多维度指标,同时探讨海洋动物鳃、肠道等多个器官中海洋微生物的分布、富集和传播模式,为研究海洋动物宿主微生物与环境浮游微生物群落之间的互作关系及生态影响提供更深入、全面的数据支持。

  • 图  1   扁吻鱼胚胎发育

    a. 成熟卵;b. 受精卵;c. 1细胞期;d. 2细胞期;e. 4细胞期;f. 8细胞期;g. 16细胞期;h. 桑葚期;i. 囊胚早期;j. 囊胚中期;k. 囊胚晚期;l. 原肠早期;m. 原肠中期;n. 原肠晚期;o. 神经胚期;p. 胚孔闭合期;q. 视囊出现期;r. 肌肉效应期;s. 耳石形成期;t. 孵化前期

    Figure  1.   Embryonic development of A. laticeps

    a. Mature egg; b. Fertilized egg; c. 1-cell stage; d. 2-cell stage; e. 4-cell stage; f. 8-cell stage; g. 16-cell stage; h. Morula stage; i. Early blastula stage; j. Middle blastula stage; k. Late blastula stage; l. Early gastrula stage; m. Mid-gastrula stage; n. Late gastrula stage; o. Neurula stage; p. Blastopore closed stage; q. Optic capsule stage; r. Muscular effect stage; s. Otolith formation stage; t. Prehatching stage

    图  2   扁吻鱼仔鱼卵黄囊体积与出膜天数的关系

    Figure  2.   Relationship between yolk sac volume and days of hatching of A. laticeps

    图  3   扁吻鱼仔稚鱼发育

    a. 初孵仔鱼;b. 出膜2 d;c. 出膜9 d;d. 出膜26 d;e. 出膜45 d;f. 出膜86 d;g. 1龄幼鱼;h. 2龄幼鱼

    Figure  3.   Delopment of A. laticeps larva and juvenile

    a. Newly hatched larva; b. 2 d post hatching; c. 9 d post hatching; d. 26 d post hatching; e. 45 d post hatching; f. 86 d post hatching; g. yearling; h. Two-year-old of young fish

    图  4   扁吻鱼前期仔鱼全长

    Figure  4.   Total length of A. laticeps prelarva

    图  5   后期仔鱼、稚鱼全长

    Figure  5.   Total length of A. laticeps postlarva and juvenile

    表  1   扁吻鱼胚胎发育时序

    Table  1   Process of embryonic development of A. laticeps

    发育时期
    Developmental stage
    耗时
    Time duration/h
    积温
    Accumulative temperature/℃·h
    图序
    Plate
    受精卵 Fertilized egg stage 0 0 图1-a~b
    胚盘形成期 Blastodem stage 1.25 20.34 图1-c
    卵裂期 Cleavage stage 5.12 83.25 图1-d~h
    囊胚期 Blastula stage 9.55 155.38 图1-i~k
    原肠期 Gastrula stage 7.68 125.01 图1-l~n
    神经胚期 Neural stage 17.18 279.57 图1-o~p
    器官形成期 Organogenesis stage 47.67 775.54 图1-q~s
    出膜前期 Prehatching stage 27.27 443.63 图1-t
    下载: 导出CSV

    表  2   其他裂腹鱼类发育特点比较

    Table  2   Comparison of developmental characteristics between other Schizothoracinaes

    种类
    Species
    卵径
    Egg diameter/mm
    卵色
    Egg color
    初孵仔鱼长度
    Length at
    hatching/mm
    孵化温度
    Rearing
    temperature/℃
    文献
    References
    吸水前
    Before hydration
    吸水后
    After hydration
    塔里木裂腹鱼 S. biddulphi 1.70±0.06 2.90±0.07 7.29±0.18 16~19 [7]
    金黄色、青灰色 7.0~8.0 14~19 [9]
    齐口裂腹鱼 S. prenanti 2.9~3.0 4.2 金黄色 11 16.8 [10]
    斑重唇鱼 D. maculates 金黄色 1.0~1.2 4.5 [11]
    光唇裂腹鱼 S. lissolabiatus 2.2±0.13 3.2±0.21 金黄色 8.5 15~17 [12]
    尖裸鲤 O. stewartii 2.57±0.07 3.22±0.07 黄色 10.27±0.15 9.5~11.8 [13]
    拉萨裸裂尻鱼 S. younghusbandi 2.50±0.08 3.54±0.07 金黄色 10.86±0.33 9.5~11.1
    拉萨裂腹鱼 S. waltoni 2.95±0.09 4.03±0.06 棕色 10.67±0.17 10.0~12.0
    短须裂腹鱼 S. wangchiachii 2.36 3.68 淡黄色 8.7 14±1 [14]
    长丝裂腹鱼 S. dolichonema 2.89±0.30 4.17±0.09 淡黄色 10.91 15.9±0.6 [15]
    小裂腹鱼 S. parvus 1.9 3 淡黄色 7.5 14.5 [16]
    双须叶须鱼 P. dipogon 3.7~3.9 5.1~5.3 淡黄色 12.4 10.0 [17]
    异齿裂腹鱼 S. o'connori 2.92±0.16 3.73±0.31 淡黄色 8.97±0.31 13.0~14.5 [18]
    重口裂腹鱼 S. davidi 2.4±0.1 橘黄色 11.67±0.38 15 [22]
    下载: 导出CSV
  • [1] 郭焱, 张人铭, 蔡林刚, 等. 新疆鱼类志[M]. 乌鲁木齐: 新疆科学技术出版社, 2012: 122.
    [2] 马燕武, 张人铭, 吐尔逊, 等. 温度和水流对扁吻鱼胚胎发育的影响[J]. 淡水渔业, 2008, 38(1): 57-61. doi: 10.3969/j.issn.1000-6907.2008.01.014
    [3] 马燕武, 张人铭, 李红, 等. 扁吻鱼的栖息地及其群落保护生物学研究[J]. 水生态学杂志, 2010, 31(1): 38-42.
    [4]

    EINBINDER S, PEREIBERG A, BEN-SHAPRU O, et al. Effects of artificial reefs on fish grazing in their vicinity: evidence from algae presentation experiments[J]. Mar Environ Res, 2006, 61(1): 110-119. doi: 10.1016/j.marenvres.2005.07.001

    [5] 殷名称. 鱼类早期生活史研究与其进展[J]. 水产学报, 1995, 19(4): 335-342.
    [6] 邵建春, 刘春雷, 秦芳, 等. 汉江地区翘嘴鲌胚胎及仔鱼发育观察[J]. 华中农业大学学报, 2016, 35(6): 111-116.
    [7] 龚小玲, 崔忠凯, 吴敏芝, 等. 塔里木裂腹鱼胚胎和仔鱼的发育与生长[J]. 上海海洋大学学报, 2013, 22(6): 827-834.
    [8] 王亚利, 刘玥, 田佳佳, 等. 四川华鳊仔稚鱼生长与形态发育研究[J]. 渔业科学进展, 2020, 41(4): 58-67.
    [9] 谢春刚, 张人铭, 马燕武, 等. 塔里木裂腹鱼人工繁殖技术初步研究[J]. 干旱区研究, 2010, 27(5): 734-737.
    [10] 吴青, 王强, 蔡礼明, 等. 齐口裂腹鱼的胚胎发育和仔鱼的早期发育[J]. 大连水产学院学报, 2004(3): 218-221.
    [11] 牛建功, 张涛, 刘鸿, 等. 斑重唇鱼人工驯养繁育技术研究[J]. 科学养鱼, 2016(12): 10-11.
    [12] 申安华, 李光华, 赵树海, 等. 光唇裂腹鱼胚胎发育与仔鱼早期发育的研究[J]. 水生态学杂志, 2013, 34(6): 76-80. doi: 10.3969/j.issn.1674-3075.2013.06.013
    [13] 许静. 雅鲁藏布江四种特有裂腹鱼类早期发育的研究[D]. 武汉: 华中农业大学, 2011: 39-44.
    [14] 左鹏翔, 李光华, 冷云, 等. 短须裂腹鱼胚胎与仔鱼早期发育特性研究[J]. 水生态学杂志, 2015, 36(3): 77-82.
    [15] 刘小帅, 王红梅, 甘维熊, 等. 雅砻江长丝裂腹鱼胚胎形态发育及仔鱼生长研究[J]. 安徽农业科学, 2017, 45(31): 118-121. doi: 10.3969/j.issn.0517-6611.2017.31.037
    [16] 冷云, 徐伟毅, 刘跃天, 等. 小裂腹鱼胚胎发育的观察[J]. 水利渔业, 2006, 26(1): 32-33. doi: 10.3969/j.issn.1003-1278.2006.01.015
    [17] 刘海平, 刘孟君, 牟振波, 等. 西藏双须叶须鱼早期发育特征[J]. 水生生物学报, 2019, 43(5): 1041-1055. doi: 10.7541/2019.123
    [18] 张良松. 异齿裂腹鱼胚胎发育与仔鱼早期发育的研究[J]. 大连海洋大学学报, 2011, 26(3): 238-242. doi: 10.3969/j.issn.1000-9957.2011.03.009
    [19] 任波, 任慕莲, 郭焱, 等. 扁吻鱼的生物学研究[J]. 水产学杂志, 2006, 19(2): 9-22. doi: 10.3969/j.issn.1005-3832.2006.02.002
    [20] 任波, 任慕莲, 郭焱, 等. 扁吻鱼胚胎及仔鱼发育的形态学观察[J]. 大连水产学院学报, 2007, 22(6): 397-402.
    [21] 张人铭, 马燕武, 吐尔逊, 等. 新疆扁吻鱼的胚胎发育和仔鱼发育的初步观察[J]. 干旱区研究, 2008, 25(2): 190-195.
    [22] 彭淇, 吴彬, 陈斌, 等. 野生重口裂腹鱼[Schizothorax (Racoma) davidi (Sauvage)]的性腺发育观察与人工繁殖研究[J]. 海洋与湖沼, 2013, 44(3): 651-655. doi: 10.11693/hyhz201303017017
    [23] 金丹璐, 张清科, 王友发, 等. 鲤科经济鱼类马口鱼 (Opsariichthys bidens) 胚胎发育及仔稚鱼形态与生长观察研究[J]. 海洋与湖沼, 2017, 48(4): 838-847.
    [24]

    BOBE J, LABBÉ C. Egg and sperm quality in fish[J]. Gen Comp Endocrinol, 2010, 165(3): 535-548. doi: 10.1016/j.ygcen.2009.02.011

    [25] 甘维熊, 王红梅, 邓龙君, 等. 雅砻江短须裂腹鱼胚胎和卵黄囊仔鱼的形态发育[J]. 动物学杂志, 2016, 51(2): 253-260.
    [26] 蔡瑞钰, 赵健蓉, 黄静, 等. 云南盘仔稚鱼发育的初步观察[J]. 南方水产科学, 2018, 14(3): 120-125. doi: 10.3969/j.issn.2095-0780.2018.03.015
    [27] 徐滨, 朱祥云, 魏开金, 等. 拉萨裂腹鱼的胚胎及仔稚鱼发育特征[J]. 大连海洋大学学报, 2020, 35(5): 663-670.
    [28] 严太明, 周翠萍, 李忠利, 等. 实验室饲养的重口裂腹鱼仔鱼的形态发育与生长[J]. 四川农业大学学报, 2007, 25(4): 493-497, 525. doi: 10.3969/j.issn.1000-2650.2007.04.025
    [29] 骆豫江, 朱新平, 潘德博, 等. 高体革䱨仔稚鱼的生长和发育[J]. 水产学报, 2008, 32(5): 697-702.
    [30] 黄俊, 朱挺兵, 杨德国, 等. 短须裂腹鱼仔稚鱼发育及生长特性的初步研究[J]. 水生态学杂志, 2019, 40(6): 99-105.
推荐阅读
Application of fishery acoustic frequency difference technology in fishery resource assessment of marine ranching in southern sea area of yintan, guangxi
CHEN Guobao et al., SOUTH CHINA FISHERIES SCIENCE, 2025
Trophic niches ofceratoscopelus warmingiiandbolinichthys longipesin adjacent waters of cold seep area in south china sea
SHI Juan et al., SOUTH CHINA FISHERIES SCIENCE, 2024
Estimation on catches by guangxi fishing vessels in beibu gulf
GONG Binhao et al., SOUTH CHINA FISHERIES SCIENCE, 2024
Assessment of fishery resources in southern sea area of yintan marine ranching, guangxi province
NIU Lulian et al., SOUTH CHINA FISHERIES SCIENCE, 2024
Study on community structure and seasonal dynamics of zooplankton and macrobenthos in rice-fish symbiotic system in qingtian, zhejiang province
ZHENG Sidi et al., JOURNAL OF FISHERIES OF CHINA, 2025
Difference in the fishery resource density between the bottom and surface layers and an analysis of multiple types of related factor importance in the northern south china sea
SUN Mingshuai et al., JOURNAL OF FISHERY SCIENCES OF CHINA, 2024
The diversity and community composition of three plants' rhizosphere fungi in kaolin mining areas
Xiao, Wenqi et al., JOURNAL OF FUNGI, 2024
Spatio-temporal distribution of harmful algal blooms and their correlations with marine hydrological elements in offshore areas, china
Chen, Chao et al., OCEAN & COASTAL MANAGEMENT, 2023
Morphological and phylogenetical analyses of pathogenic hypomyces perniciosus isolates from agaricus bisporus causing wet bubble disease in china
PHYTOTAXA, 2021
Exploring water quality as a determinant of the existence of soft shell crab (scylla olivacea) in different habitats of the large marine ecosystem in the bay of bengal.
WATER CONSERVATION AND MANAGEMENT, 2024
Powered by
图(5)  /  表(2)
计量
  • 文章访问数:  1551
  • HTML全文浏览量:  611
  • PDF下载量:  123
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-08-22
  • 修回日期:  2020-11-01
  • 录用日期:  2020-11-18
  • 网络出版日期:  2020-12-04
  • 刊出日期:  2021-02-04

目录

/

返回文章
返回