Effect of mosaic disease on skin mucosa microbial community structure of Chinese sturgeon (Acipenser sinensis)
-
摘要:
花斑病是近年来中华鲟 (Acipenser sinensis) 养殖过程中出现的一种新型疾病。为探索用于监测中华鲟健康状况的菌群标志物，通过高通量测序，分析比较了健康和患花斑病中华鲟幼鱼皮肤黏膜层微生物菌群的结构特征。结果显示，与健康中华鲟相比，患病中华鲟皮肤黏膜层微生物多样性和丰富度均显著下降 (p<0.001)。在属分类水平上，健康中华鲟皮肤黏膜层的优势菌属为乳杆菌属 (Lactobacillus)、不动杆菌属 (Acinetobacter)、norank_f_Muribaculaceae、明串珠菌属 (Leuconostoc) 和假单胞菌属 (Pseudomonas)；患病中华鲟皮肤黏膜层的优势菌属为黄杆菌属 (Flavobacterium) 和不动杆菌属，其中黄杆菌属占比78.38%。通过差异显著性分析，获得了5种反映中华鲟健康状态的敏感菌群。研究表明，与健康中华鲟相比，花斑病破坏了中华鲟皮肤黏膜层正常的微生态稳态结构，优势菌群由乳杆菌属、不动杆菌属等转变为黄杆菌属；因此黄杆菌属可作为监测中华鲟花斑病病原的菌群标志物，其相对丰度可用于评估中华鲟的患病风险。
Abstract:Mosaic disease is a novel disease that has emerged during the culture of Chinese sturgeon (Acipenser sinensis) in recent years. In order to study the flora markers for monitoring the health status of A. sinensis, we analyzed the microbial community structure of the skin mucosa of healthy and diseased juveniles by using high-throughput sequencing technology. The results show that compared with healthy fish, the diversity and richness of diseased fish reduced significantly (p<0.001). On genus level, the skin mucosa of healthy fish was composed of Lactobacillus, Acinetobacter, norank_f-Muribaculaceae, Leuconostoc and Pseudomonas, while the skin mucosa of diseased fish was composed of Flavobacterium and Acinetobacter, and the proportion of Flavobacterium was 78.38%. The significant difference analysis shows that five sensntive taxa were screened, reflecting the health status of Chinese sturgeon. This study indicates that compared with healthy fish, mosaic disease destroyed the normal homeostasis of microbiota in skin mucosa. The dominant bacterial community changed from Lactobacillus and Acinetobacter to Flavobacterium. Therefore, Flavobacterium have been identified as a marker to monitor the pathogen of mosaic disease in A. sinensis. Its relative abundance can be utilized to evaluate the risk of disease in A. sinensis.
-
粒径谱 (Size spectrum) 是以能量生态学为基础，把各类复杂的水域生态系统、食物循环、种群变动和影响机制等简化为谱图的一种方法[1]。它不依赖于传统的生物学分类知识[2]，具有数据需求少、易操作和结果直观等特点，避免了生物特性鉴定不准确或调查资源有限引起的误差，在总体上、宏观地研究不同生态系统和群落动态及其变化机制和控制因素[1]，可以用来比较不同历史时期和环境条件下生态系统的特点和差异[3]，从某一粒级估计其他粒级的生物量或产量[1-2]，评估捕捞[4]或环境对生态系统的压力[5-6]，为研究水生生态学提供了新思路。1988年王荣等[7]首次把粒径谱的概念引入国内，如今我国粒径谱研究的发展及应用越来越广泛，涉及海湾[8]、湿地[9]和河口[10]等生境，在浮游生物[8,10]、底栖生物[9]中已有一定的应用和研究成果，而粒径谱应用于鱼类群落研究[11]近年才开始起步。
浙江省南部近海位于暖温带海区的温州和台州近海，有多条江河入海，受到浙江沿岸水、台湾暖流和黑潮次表层水等水团影响，悬浮物和营养物充足，饵料生物丰富[12]，为多种鱼类提供了必要的栖息、索饵和繁衍场所。近年来，由于过度捕捞、环境污染等问题，近海水域的渔业资源严重衰退，渔业群落结构发生较大变动[13]。目前，国内已有学者对浙江省南部近海鱼类群落进行了物种生态优势度[14]、多样性[15]、生态位[16]等方面的研究，但受到气候变化和捕捞压力的持续影响，该海域鱼类群落结构发生着动态变化，需要进一步掌握其群落特征。本文根据2019年浙江南部近海4个季节的渔业资源调查数据，从时间和空间上进行鱼类群落粒径谱分析，探究该海域鱼类群落粒径谱特征，为鱼类资源保护与生态可持续发展提供参考。
1. 材料与方法
1.1 采样与测定
研究区域为浙江南部近岸海域 (120.93°E—122.95°E, 27.21°N—28.97°N)，采样时间为2019年2月26日—3月2日 (冬季)、5月23—27日 (春季)、8月1—7日 (夏季) 和11月1—4日 (秋季)，共4个航次的渔业资源和环境调查 (图1)。
![]()
图 1 浙江南部近海采样站位Figure 1. Sampling stations in offshore waters of southern Zhejiang Province调查船“浙洞渔10109”，主机马力540 HP，总吨位305 t，取样网具为网口宽40 m、网高7.5 m、网囊网目2 cm的底拖网。其中冬、夏季各27个站位，春季26个站位，秋季24个站位 (图1)。调查时每站位拖曳1次，平均拖速3 kn，每次作业时间1 h。同时，使用温盐深测量仪 (Conductivity temperature depth, CTD) 于每个站位测定盐度、温度、深度和pH等环境数据 (表1)。现场采样中样品的采集、保存和测定等均按照《海洋监测规范 (GB 17378)》和《海洋调查规范 (GB/T 12763)》执行。
表 1 浙江南部近海环境指标统计值Table 1. Statistical values of environmental indicators in offshore waters of southern Zhejiang Province季节
Season环境因子
Environmental factor最小值
Minimum value最大值
Maximum value平均值±标准差
$\overline { X}\pm { \rm {SD}} $春季 Spring 水深 Water depth/m 20.00 66.50 47.66±14.32 水温 Temperature/℃ 21.20 24.50 22.46±0.79 盐度 Salinity 26.30 33.10 28.77±1.74 pH 8.23 8.71 8.44±0.14 夏季 Summer 水深 Water depth/m 19.90 65.20 48.50±14.18 水温 Temperature/℃ 28.00 30.20 28.94±0.45 盐度 Salinity 31.70 34.00 33.36±0.71 pH 8.09 8.32 8.17±0.06 秋季 Autumn 水深 Water depth/m 16.70 63.70 51.21±12.25 水温 Temperature/℃ 21.60 23.90 22.96±0.58 盐度 Salinity 28.40 33.20 31.64±1.40 pH 8.15 8.27 8.22±0.04 冬季 Winter 水深 Water depth/m 21.50 66.00 48.26±13.63 水温 Temperature/℃ 10.40 15.80 13.21±1.80 盐度 Salinity 28.80 34.30 32.24±1.60 pH 8.18 8.34 8.26±0.04 1.2 研究方法
1.2.1 Sheldon鱼类粒径谱
鱼类个体生物量 (g) 经过log2转换后划分粒级，以最小生物量为起点，每一个粒级的上限值是下限值的2倍[17]，即粒径间距以公比2成等比形式增长，作为横轴；单位面积 (km2)上对应粒级的生物量经过log2转换后作为纵轴，以此构建鱼类Sheldon粒径谱。谱上的散点用曲线连接后形成的峰型可以反映鱼类群落的结构特征[18-19]。
1.2.2 标准化鱼类粒径谱
鱼类标准化粒径谱的横轴是粒径级上限值经过log2转换后划分粒级作为横轴；每个单位面积上的鱼类总生物量与该粒级上鱼类个体生物量的变化幅度的比值，经log2对数化后作为纵轴，并对谱线进行相关性和回归分析，实质上代表了鱼类群落的粒级丰度[18-20]。
$$ \mathrm{\beta }\left({w}_{i}\right)={\rm{log}}_{2}\left[\frac{{B}\left({{w}}_{i}\right)}{\Delta {{w}}_{i}}\right] $$ (1) 式中：
$ \beta \left({w}_{i}\right) $ 为第$ i $ 粒级所对应的纵坐标；$ B\left({w}_{i}\right) $ 为单位面积上第$ i $ 粒级的总生物量；$ \Delta {w}_{i} $ 为第$ i $ 粒级上个体生物量的变化幅度。在理想的稳定状态下，标准化粒径谱为线性，且斜率为−1[2]。当鱼类群落受外界干扰打破平衡时，标准化粒径谱为抛物线型[11,21]，曲线随捕捞强度增大而更加陡峭[18]，且曲率或斜率大小受到当年生幼鱼和鱼类粒径大小[22]、人类捕捞活动强度[4]、栖息地环境[18]等因素影响。
1.2.3 多维尺度法
为了揭示不同站点鱼类群落的相似度结构，筛选出同时在3个或以上站点出现且丰度占1个站位1%以上的物种用于分析[23]，以平衡稀有种、偶见种对整个鱼类群落的影响，再用该数据计算Bray-Curtis相似矩阵，在该矩阵的基础上，用多维尺度法 (Multidimensional Scaling, MDS) 对浙江省南部近岸海域的鱼类群落空间特征进行分析。使用胁强系数 (Stress) 检验MDS结果，一般情况下Stress<0.2时，可用MDS的二维散点图来表示，其图形有一定意义[24]。为判断分组的统计学意义，用相似性分析 (Analysis of similarities, ANOSIM) 检验组间的差异显著性，Global R>0为有效分组，表示组间差异大于组内差异。
1.3 数据处理与分析
鱼类生物学原始数据进行拖网时间 (1 h) 和拖网速度 (3 kn) 标准化，在R 4.0.2中完成鱼类粒径谱建立、显著性检验以及环境因子与标准化粒径谱曲率的相关性分析，鱼类群落分组用Primer 6.15进行MDS分析，采用Primer 6.15的ANOSIM和SPSS 19.0的单因素方差分析 (One-way ANOVA) 检验粒径谱特征差异和各分组间差异的显著性。
2. 结果
2.1 全年鱼类粒径谱
全年Sheldon鱼类粒径谱形状呈不规则的锯齿状 (图2-a)，大体为单峰模式，粒级介于−2~11 (表2)。在−2~6粒级单位面积的生物量呈上升趋势，6~7粒级生物量增加趋势放缓，在7~9粒级呈下降趋势，9~11粒级缓慢增加。生物量最高值出现在6~7粒级，主要由带鱼 (Trichiurus lepturus)、刺鲳 (Psenopsis anomala)、镰鲳 (Pampus echinogaster) 等组成 (表2)。其中带鱼的生物量和数量占整个粒级上生物量和数量百分比的36.04%和36.57%，平均个体生物量为69 g；刺鲳的生物量和数量占比分别为12.02%和11.33%，平均个体生物量为75 g；镰鲳的生物量和数量占比分别为10.49%和9.70%，平均个体生物量为76 g。
![]()
图 2 浙江南部近海全年鱼类粒径谱a. Sheldon粒径谱；b. 标准化粒径谱。Figure 2. Annual fish size spectrum in offshore waters of southern Zhejiang Provincea. Sheldon size spectrum; b. Normalized size spectrum.表 2 浙江南部近海鱼类群落组成及营养级状况Table 2. Composition of fish community and trophic level in offshore waters of southern Zhejiang Province粒径
Size range粒径级上限值
Upper limit of
size class/g总生物量
Total biomass/
(g·km−2)主要鱼类物种
Main fish species营养级
Trophic level
$\overline X \pm {\rm{SD}}$−3~−2 −2 1 麦氏犀鳕 Bregmaceros mcclellandi 3.30±0.42* −2~−1 −1 2 麦氏犀鳕 B. mcclellandi 3.30±0.42* 粗吻海龙 Trachyrhamphus serratus — −1~0 0 9 麦氏犀鳕 B. mcclellandi 3.30±0.42* 六丝钝尾虾虎鱼 Amblychaeturichthys hexanema 3.65±0.05[25] 舒氏冠海龙 Corythoichthys schultzi 3.80±0.50* 0~1 1 40 麦氏犀鳕 B. mcclellandi 3.30±0.42* 六丝钝尾虾虎鱼 A. hexanema 3.65±0.05[25] 细条天竺鲷 Jaydia lineata 3.70±0.50* 赤鼻棱鳀 Thryssa kammalensis 3.91±0.31[25] 鳄齿鱼 Champsodon capensis 4.20±0.73* 1~2 2 163 六丝钝尾虾虎鱼 A. hexanema 3.65±0.05[25] 细条天竺鲷 J. lineata 3.70±0.50* 麦氏犀鳕 B. mcclellandi 3.30±0.42* 拉氏狼牙虾虎鱼 Odontamblyopus lacepedii — 2~3 3 613 拉氏狼牙虾虎鱼 O. lacepedii — 龙头鱼 Harpadon nehereus 3.62±0.20[25] 发光鲷 Acropoma japonicum − 细条天竺鲷 J. lineata 3.70±0.50* 大头白姑鱼 Pennahia macrocephalus 4.10 ±0.64* 蓝圆鲹 Decapterus maruadsi 3.63±0.06[25] 带鱼 T. lepturus 3.76±0.13[25] 赤鼻棱鳀 Thryssa kammalensis 3.91±0.31[25] 3~4 4 2 778 赤鼻棱鳀 T. kammalensis 3.91±0.31[25] 大头白姑鱼 P. macrocephalus 4.10 ±0.64* 白姑鱼 Pennahia argentata 4.10 ±0.70* 发光鲷 A. japonicum — 龙头鱼 H. nehereus 3.62±0.20[25] 4~5 5 8 507 黄鲫 Setipinna tenuifilis 3.74±0.15[25] 黑姑鱼 Atrobucca nibe 4.00±0.69* 大头白姑鱼 P. macrocephalus 4.10 ±0.64* 龙头鱼 H. nehereus 3.62±0.20[25] 5~6 6 17 971 小黄鱼 Larimichthys polyactis 3.70±0.40* 龙头鱼 H. nehereus 3.62±0.20[25] 带鱼 T. lepturus 3.76±0.13[25] 刺鲳 P. anomala 4.00±0.28* 6~7 7 19 326 带鱼 T. lepturus 3.76±0.13[25] 刺鲳 P. anomala 4.00±0.28* 镰鲳 P. echinogaster 3.50±0.20* 7~8 8 10 651 带鱼 T. lepturus 3.76±0.13[25] 镰鲳 P. echinogaster 3.50±0.20* 刺鲳 P. anomala 4.00±0.28* 8~9 9 4 018 镰鲳 P. echinogaster 3.50±0.20* 蓝圆鲹 D. maruadsi 3.63±0.06[25] 绿鳍鱼 Chelidonichthys kumu 3.70±0.57* 黄鳍东方鲀 Takifugu xanthopterus — 9~10 10 5 234 黄鮟 Lophius litulon 4.17±0.09[25] 中国花鲈 Lateolabrax maculatus — 10~11 11 5 867 黄鮟 L. litulon 4.17±0.09[25] 注：*. 数据来源于FishBase (www.fishbase.org)；—. 无数据。
Note: *. Data from FishBase (www.fishbase.org); —. No data.全年标准化鱼类粒径谱 (图2-b) 的拟合曲线呈“穹顶”抛物线状，回归系数R2为0.88 (P<0.01)，曲线穹顶位于3~7粒级，曲率为−0.15，表明浙江南部近海鱼类群落总体处于受干扰状态。
2.2 鱼类粒径谱的季节差异
2.2.1 Sheldon鱼类粒径谱
浙江南部近海4个季节Sheldon型鱼类生物量粒径谱形状基本呈单峰模式，季节差异不显著 (P=0.93，图3)。4个季节鱼类粒级介于−2~11，其中冬季粒级介于−1~11，春季为0~10，夏季为−2~10，秋季为0~9；鱼类总生物量以夏季 (8月) 最大，春季 (5月) 最小。冬季，Sheldon峰值位于10~11粒级，主要由黄鮟 (Lophius litulon)组成，其渔获重量百分比和数量百分比分别占该粒级的68.95%和66.67%，生物量最高值位于5~6粒级，主要由小黄鱼 (Larimichthys polyactis)、黄鲫 (Setipinna tenuifilis)构成，小黄鱼的渔获重量百分比和数量百分比分别为34.51%和28.34%，黄鲫的相应占比分别为18.93%和22.81%。春季，鱼类生物量最高值位于5~6粒级，主要由带鱼、龙头鱼 (Harpadon nehereus)构成，其中带鱼的渔获质量百分比和数量百分比分别为38.18%和35.55%，龙头鱼的相应占比分别为18.57%和20.09%；夏季，鱼类生物量最高值位于6~7粒级，主要由带鱼和日本鲭 (Scomber japonicus) 组成，其中带鱼的渔获质量百分比和数量百分比分别为50.60%和49.73%，日本鲭的相应占比分别为15.55%和14.32%；秋季，鱼类生物量最高值位于6~7粒级，主要由刺鲳、镰鲳等组成，其中刺鲳的渔获质量百分比和数量百分比分别为28.27%和24.23%，镰鲳的相应占比分别为18.97%和18.37%。
![]()
图 3 浙江南部近海的各季节Sheldon鱼类粒径谱Figure 3. Sheldon fish size spectra by seasons in offshore waters of southern Zhejiang Province2.2.2 标准化鱼类粒径谱
浙江省南部近岸海域4个季节标准化鱼类粒径谱的拟合曲线均呈“穹顶”抛物线型 (图4)，所有回归结果均极显著 (P<0.01)。回归系数R2以冬季最大 (R2=0.90)，夏季最小 (R2=0.75，表3)；曲率介于−0.18~−0.09，以冬季最大 (−0.09)，夏秋次之，春季最小 (−0.18，表3)，且冬季曲线相对平缓，春季曲线较为陡峭，这表明春季浙江南部近海鱼类生物量受外界干扰程度最大，而冬季相对稳定。
![]()
图 4 浙江南部近海的各季节标准化鱼类粒径谱Figure 4. Normalized fish size spectra in offshore waters of southern Zhejiang Province in different seasons表 3 标准化鱼类粒径谱各季拟合曲线的参数Table 3. Parameters of seasonal fitted curve of normalized fish size spectra季节
Season拟合方程
Fitted
equation曲率
Curvature回归系数R2
Regression
coefficient冬季 Winter y=−0.09x2+0.83x+4.54 −0.09 0.90 春季 Spring y=−0.18x2+1.70x+2.63 −0.18 0.87 夏季 Summer y=−0.13x2+1.07x+5.08 −0.13 0.75 秋季 Autumn y=−0.16x2+1.51x+3.59 −0.16 0.88 2.3 鱼类群落空间差异
对浙江省南部近海四季各站位鱼类丰度进行多维尺度分析，得出四季MDS图 (图5)。Stress值分别为0.13、0.17、0.18和0.13，均在0.2以下，表明MDS结果可以反映各站位间的关系。冬季27个站位以69%的群落相似性将鱼类群落分为组Ⅰ、组Ⅱ、组Ⅲ、组Ⅳ，春季以57%的相似性分为4组，夏季以48%的相似性分为3组，秋季以64%的相似性分为3组。其中冬季组Ⅰ(站位21)、春季组Ⅰ(站位42) 和夏季组Ⅰ(站位04) 在各自季节自成一组，聚类相似性接近0。ANOSIM结果显示，冬季Global R=0.83 (P<0.01)，春季Global R=0.74 (P<0.01)，夏季Global R=0.69 (P<0.01)，秋季Global R=0.85 (P<0.01)，说明各分组间差异性极显著。
2.4 鱼类粒径谱的空间差异
根据MDS结果，构建四季14个群落组的Sheldon粒径谱 (图6)，每个群落组均呈不规则的锯齿状。除夏季外，远岸站位比近岸站位的粒级区间大、生物量高。
![]()
图 6 浙江南部近海不同鱼类群落Sheldon粒径谱Figure 6. Sheldon fish size spectra of different fish communities in offshore waters of southern Zhejiang Province冬季，组Ⅰ (站位21) 粒级介于4~10，主要由带鱼等鱼种组成；组Ⅱ(站位33、34) 粒级介于1~7，主要以六丝钝尾虾虎鱼 (Amblychaeturichthys hexanema)等为主；组Ⅲ (01、06、11等近岸站位) 粒级介于−1~10，以丰度高的拉氏狼牙虾虎鱼 (Odontamblyopus lacepedii)、黄鲫等小粒级鱼类为主；组Ⅳ (03、04、13等远岸站位) 粒级介于0~11，出现黄鮟等生物量较高鱼类。
春季，组Ⅰ(站位42) 粒级介于4~9，以龙头鱼等种类为主；组Ⅱ (32、33、23、24等远岸站位) 粒级介于0~10，总生物量最大，主要由带鱼等组成；组Ⅲ (39、35等近岸站位) 粒级介于2~9，峰值出现在5~7粒级上，以小黄鱼等为主；组Ⅳ (01、21等近岸站位) 粒级介于2~9，以黄鲫等鱼类为主。
夏季，组Ⅰ(站位04) 粒级介于4~8，在5~6粒级上捕获大量白姑鱼 (Pennahia argentata)，生物量较高形成峰值，而在6~7粒级上捕获1尾短尾大眼鲷 (Priacanthus macracanthus)，生物量低于5~6和7~8粒级上的生物量，致使波谷出现；组Ⅱ(33、14、42等远岸站位) 粒级介于−2~9，最小粒级主要由麦氏犀鳕 (Bregmaceros mcclellandi)、六丝钝尾虾虎鱼等低生物量鱼类组成；组Ⅲ (01、31、35等近岸站位) 粒级介于−1~10，捕获1尾海鳗 (Muraenesox cinereus) 使9~10粒级上生物量较高。
秋季，组Ⅰ(站位31、35) 粒级介于3~8，以龙头鱼居多；组Ⅲ (24、32、42等远岸站位) 粒级介于0~9，较群落组Ⅱ (01、06、16等近岸站位) 粒级 (1~9) 更完整，峰值出现在粒级7~8上，主要由刺鲳、带鱼、镰鲳等组成。
14个鱼类群落组的标准化粒径谱 (图7) 均存在显著差异 (P<0.01)，除冬季群落组Ⅰ近似线性外，其他各组拟合曲线均呈穹顶型。冬季，拟合曲线的曲率为−0.17~−0.01，R2介于0.68~0.89 (表4)，其中组Ⅰ拟合曲线的曲率最大且近似直线，其次为组Ⅳ和组Ⅲ，组Ⅱ最小，生物量由组Ⅳ向组Ⅰ依次减少。春季，拟合曲线的曲率为−0.29~−0.19，R2介于0.31~0.80 (表4)，其中组Ⅱ拟合曲线的曲率最大，其次为组Ⅲ和组Ⅰ，组Ⅳ最小；组Ⅱ生物量最大，其次是组Ⅲ、Ⅳ、Ⅰ。夏季，拟合曲线的曲率为−0.38~−0.10，R2介于0.57~0.82 (表4)，其中组Ⅲ拟合曲线的曲率最大，其次为组Ⅱ，组Ⅰ最小；组Ⅲ生物量最大，组Ⅰ生物量最小。秋季，拟合曲线的曲率为−0.24~−0.15，R2介于0.26~0.88 (表4)，其中组Ⅰ拟合曲线的曲率最大，其次为组Ⅲ，组Ⅱ最小；组Ⅲ生物量最大，组Ⅰ生物量最小。四季各群落组标准化粒径谱的曲率大小表明，除夏季外，近岸鱼类群落整体上比远岸鱼类群落受到更大的干扰。
![]()
图 7 浙江南部近海不同鱼类群落标准化粒径谱Figure 7. Normalized fish size spectra of different fish communities in offshore waters of southern Zhejiang Province表 4 各季节分组标准化鱼类粒径谱拟合曲线参数Table 4. Parameters of seasonal fitted curve of normalized fish size spectra in different groups季节
Season分组
Group拟合方程
Fitted equation曲率
Curvature回归系数R2
Regression coefficient生物量
Biomass/(g·km−2)冬季 Winter 组Ⅰ y=−0.01x2+0.64x+0.92 −0.01 0.89 8 004 组Ⅱ y=−0.17x2+1.33x+3.44 −0.17 0.68 4 291 组Ⅲ y=−0.10x2+0.93x+4.27 −0.10 0.85 17 109 组Ⅳ y=−0.08x2+0.70x+4.99 −0.08 0.75 28 858 春季 Spring 组Ⅰ y=−0.28x2+3.83x−8.03 −0.28 0.31 5 232 组Ⅱ y=−0.19x2+1.80x+3.10 −0.19 0.74 24 650 组Ⅲ y=−0.28x2+2.78x+0.19 −0.28 0.80 19 059 组Ⅳ y=−0.29x2+2.89x−0.57 −0.29 0.78 9 667 夏季 Summer 组Ⅰ y=−0.38x2+3.21x+0.49 −0.38 0.82 8 076 组Ⅱ y=−0.13x2+1.12x+4.72 −0.13 0.65 29 358 组Ⅲ y=−0.10x2+0.92x+5.27 −0.10 0.57 34 790 秋季 Autumn 组Ⅰ y=−0.15x2+1.60x+1.05 −0.15 0.26 4 069 组Ⅱ y=−0.24x2+2.32x+1.78 −0.24 0.88 16 406 组Ⅲ y=−0.16x2+1.31x+4.03 −0.16 0.59 20 425 3. 讨论
3.1 鱼类群落特征
Sheldon鱼类粒径谱的峰型主要反映海域生态系统对外界干扰的适应性响应[26]。Warwick[27]与Jung和Houde[19]的研究表明，粒径谱图应为“双峰模式”，小粒级物种和大粒级物种分别形成波峰，并于交汇处形成波谷。本研究显示，浙江南部近海2019年整体粒级介于−2~11，但仅在6~7粒级出现最高生物量 (19 326 g·km−2)，由带鱼、刺鲳、镰鲳等数量多的低生物量鱼类组成，在8~9粒级处捕获的镰鲳、蓝圆鲹等生物量 (4 018 g·km−2) 最低，低于7~8 (10 651 g·km−2) 和9~10 (5 234 g·km−2) 粒级，致使于8~9粒级处形成波谷，而不存在第二个峰值，表明该海域鱼类粒径处于不平衡状态，鱼类组成主要以小型鱼类为主，大型经济鱼类较少，这可能与过度捕捞有关。
标准化粒径谱拟合曲线的曲率、回归系数可指示群落的健康和稳定状况[28-29]。本研究中，浙江南部近海鱼类标准化粒径谱拟合曲线呈“穹顶”抛物线型，曲率为−0.15，高于稳态理论值−1，曲线较陡，说明小型鱼类居于有利地位，不易被大型鱼类捕食，向大型鱼类能量传递效率偏低，表明鱼类群落中的营养循环效率低。标准化粒径谱拟合曲线R2 (0.88) 接近1，说明各个粒级之间的能量传递效率近似。
本研究发现，东海传统“四大海产”中，鱼类除小黄鱼、带鱼能保持一定的采样率外，大黄鱼 (L. crocea) 等采样率极低。浙江省南部近岸海域鱼类群落主要以六丝钝尾虾虎鱼、麦氏犀鳕、拉氏狼牙虾虎鱼、鳄齿鱼 (Champsodon capensis) 等饵料性鱼类为主，而刀鲚 (Coilia nasus)、大黄鱼、中国花鲈等高营养级的经济鱼类少，几乎无法形成利于捕捞的渔汛期，这可能是因为不合理的捕捞活动造成亲体资源量剧减，加之近年海岸带开发、环境污染等原因，也可能与调查采样的海域范围有关。
3.2 鱼类粒径谱季节变动及其影响因素
鱼类粒径谱反映鱼类群落的粒径结构[18,30]及鱼类大小与丰度的关系[31]。本研究中，浙江南部近海冬季鱼类的粒级介于−1~11，春季为0~10，夏季为−2~10，秋季为0~9，各季节粒级大小存在季节变动。本次调查中冬季平均水温为13.21 ℃，最大粒级 (10~11) 捕获冷温性鱼类黄鮟、远东海鲂 (Zeus faber) 和中国花鲈，最小粒级 (−2~−1) 捕获虻鲉 (Erisphex pottii)、舒氏冠海龙 (Corythoichthys schultzi)、粗吻海龙 (Trachyrhamphus serratus)、六丝钝尾虾虎鱼等生物量低的鱼类和幼鱼；春季水温上升，平均值22.46 ℃，黄鮟等大粒级的冷温性鱼种游向水温较低处使粒级区间小于冬季，最小粒级 (−1~0) 捕获眼镜鱼 (Mene maculata)、翼红娘鱼 (Lepidotrigla alata)、蓝圆鲹的幼鱼，而高粒级只捕获1尾乌鲹 (Parastromateus niger)；夏季正值禁渔期，平均水温为28.94 ℃，营养盐丰富，带鱼等暖温性鱼类洄游至近海进行育幼，以幼鱼数量居多，最小粒级 (−3~−2) 以带鱼为主；秋季平均水温为22.96 ℃，在伏季休渔期结束后，人类捕捞强度逐渐增大导致高粒级经济鱼类剧减，鱼类粒级区间缩小，最小粒级 (−1~0) 捕获发光鲷、六丝钝尾虾虎鱼、丝棘虎鲉 (Minous pusillus)、麦氏犀鳕、棕腹刺鲀 (Gastrophysus spadiceus)，最大粒级有带鱼、鮸 (Miichthys miiuy)、鲶 (Silurus asotus)、棕腹刺鲀。
同时，粒径谱图形的曲率受到各季节最小粒级上鱼类生物量[9]、当年生幼鱼数量、洄游性与定栖性鱼类生物量[26]、捕捞强度[32]等影响。本研究发现，冬季浙江南部近海鱼类群落相对稳定，最小粒级上捕获数量较多的小个体鱼类，使得该粒级总生物量较高。冷温性鱼种黄鮟等个体较大的成鱼在冷水系沿岸流影响下出现，这或许是冬季存在高粒级 (10~11) 的原因，从而使得冬季鱼类群落的粒径结构较其余季节好，标准化粒径谱的曲率相对平缓。春、夏季，最小粒级上鱼类的数量和生物量均较少，秋季虽在最小粒级上捕获较多麦氏犀鳕，但随着捕捞强度的增大，高粒级鱼类缺乏，使标准化鱼类粒径谱的曲率相对陡峭。此外，小粒径的丝棘虎鲉、粗吻海龙、六丝钝尾虾虎鱼等定栖性鱼类的存在，发光鲷、麦氏犀鳕、蓝圆鲹、眼镜鱼等洄游性鱼类的季节性出现，以及带鱼等当年生鱼类的出现，都会使得谱图曲率产生季节变化。冬季到春季曲率降低与位于春季最小粒级 (0~1) 上的1尾眼镜鱼、蓝圆鲹、翼红娘鱼有关，春季生物量为3 g·km−2，低于冬季生物量 (7 g·km−2)，且粒级区间小于冬季；春季到夏季曲率升高与夏季最小粒级 (−2~0) 捕获到较多麦氏犀鳕、带鱼幼鱼、六丝钝尾虾虎鱼有关，且0~1粒级上夏季生物量 (9 g·km−2) 高于春季 (4 g·km−2)；夏季到秋季曲率下降是由于秋季最小粒级 (0~1) 上的麦氏犀鳕、发光鲷等造成，生物量 (6 g·km−2) 低于夏季生物量 (9 g·km−2)，且粒级区间较秋季小；秋季到冬季曲率上升与冬季最小粒级 (−1~0) 捕获粗吻海龙、六丝钝尾虾虎鱼、丝棘虎鲉等，少数大粒级 (10~11) 黄鮟有关，且0~1粒级生物量 (7 g·km−2) 高于秋季生物量 (6 g·km−2)。
因此，浙江南部近海鱼类粒径谱季节变化明显，主要受环境变化、最小粒径级的鱼类生物量变化、洄游性与定栖性鱼类、当年生幼鱼影响。禁渔期的季节性、人类捕捞活动的季节性强度不同也会对鱼类群落结构和生物量产生影响，进而使得鱼类粒径谱结构发生季节变化。徐姗楠等[11,26]对大亚湾和胶州湾的研究也得到类似结论。
3.3 鱼类粒径谱空间特征
鱼类群落的空间分布特征受鱼类生活习性、丰度分布、栖息地环境等影响[33-34]。对比浙江南部近海四季14个群落组的相似性发现，冬、春季在69%和57%的相似性水平上分别聚为4组，夏、秋季在48%和64%的相似性水平上分别聚为3组，即春、夏季群落相似性均低于冬、秋季，这可能与鱼类生活习性、栖息地环境等有关。秋、冬季鱼类分散索饵或越冬洄游，分布较均匀，各组结构较相似；春、夏季部分鱼类在沿岸产卵、索饵，因环境需求不同，呈分群、斑块状分布，导致各组差异较明显[33]，带鱼、龙头鱼、黄鲫等物种并非只出现在某一个站位或群聚中，这可能由于该部分鱼类适温、适盐性强，不存在空间上的相互隔离，即由鱼类丰度分布的空间差异导致[34]，本研究进一步证实了这点。
同时，鱼类群落的空间分布也与鱼类生态特征、过度捕捞、环境因素等有关[35]。综合分析14个组发现，各组间均有极显著性差异 (P<0.01)，除单个站位组的冬季组Ⅰ (站位21)、春季组Ⅰ (站位42)、夏季组Ⅰ (站位04) 外，近岸站位 (如站位06、16、26等) 和远岸站位 (如站位23、32、13等) 基本不会划分成同一组 (图5)，即近岸站位和远岸站位群落结构具有较大差异，这可能与鱼类生态特性有关，与梁海[34]研究结果一致。各组Sheldon粒径谱中，冬季最高值在组Ⅳ，春季为组Ⅱ，夏季为组Ⅲ，秋季为组Ⅲ，这4组粒级区间大、生物量高且均包含大部分远岸站位。从种类组成上看，远岸站点包含了带鱼、中国花鲈、海鳗等生物量稍高的鱼类种类，近岸站点主要由黄鲫、镰鲳、龙头鱼、黑姑鱼等丰度大的小型鱼类组成。标准化粒径谱中，冬季群落组Ⅳ、春季群落组Ⅱ和秋季群落组Ⅲ的曲率分别大于冬季组Ⅲ、春季组Ⅳ和秋季组Ⅱ，整体上，近岸比远岸受到更大干扰，鱼类群落结构发生变化。而夏季正值禁渔期，致使组Ⅱ曲率小于组Ⅲ，证实禁捕政策取得一定成效，近岸站位渔业资源得到一定恢复。研究表明，人类捕捞影响鱼类粒级结构[18]，过度捕捞导致大粒级鱼类缺乏[36]，标准化粒径谱的曲率随捕捞强度的增加而减小[18,37]。同时，本研究发现标准化粒径谱曲率与环境因子有显著相关性 (图8)，其中曲率大小与海表温度 (R=−0.27) 呈极显著负相关，与盐度 (R=0.59) 呈非常显著正相关，与pH (R=−0.46) 呈非常显著负相关，与水深深度 (R=0.11) 无明显相关性，即曲率与海域盐度、pH和温度的相关性显著。有研究认为，蓝圆鲹等鱼类受温度影响大，带鱼、龙头鱼等鱼类的越冬洄游与外侧较深海域的盐度有关[38]，环境因子通过影响鱼类生长全过程导致各粒径级对应物种组成变化，从而影响整个鱼类粒径谱形状[18]。此外，宋伦等[8,39]研究表明，海域富营养化对粒径谱的粒径结构和曲率也具有影响。
![]()
图 8 标准化粒径谱曲率与环境因子的相关性分析*. P<0.05，显著相关；**. P<0.01，极显著相关；***. P<0.001，非常显著相关。Figure 8. Correlation analysis between curvature of normalized size spectra and environmental indicators*. P<0.05, significant correlation; **. P<0.01, extremely significant correlation; ***. P<0.001, very significant correlation.3.4 展望
根据鱼类粒径谱的结构形态，对偏离理论稳定状态的系统制定修复计划，如生物操纵策略中的增殖放流，根据实际情况在群落/生态系统中补充低于拟合值的生物。小范围海域可通过定向捕捞、投放幼鱼或成鱼等措施；在大范围海域，可先将幼鱼投放在邻近较为稳定的环境中，以确保鱼类资源的亲体补充量，增强群落稳定性和抗干扰能力。本研究从粒径谱角度出发，探索了鱼类群落的粒径谱结构和组成，并从季节和空间两方面分析了鱼类群落的空间差异和变动，为浙江南部近岸海域渔业资源管理和生态可持续发展提供了重要的参考依据。
致谢：对上海海洋大学渔业资源和生态系统量化评估与管理研究室、浙江省海洋水产养殖研究所的老师和同学在样品采集和生物学分析中的付出，戴黎斌、王士聪、薛铭华等提供的帮助，谨致谢忱！
-
![]()
图 1 健康和患病中华鲟的特征和病症
Figure 1. Features and symptoms of healthy and diseased A. sinensis
![]()
图 2 健康和患病中华鲟皮肤黏膜层OTUs数量分布韦恩图
Figure 2. Venn diagram of quantity distribution of OTUs among skin mucosa from healthy and diseased A. sinensis
![]()
图 3 健康和患病中华鲟皮肤黏膜层微生物结构组成分析
Figure 3. Dominant species of bacteria in skin mucosa from healthy and diseased A. sinensis
![]()
图 4 健康和患病中华鲟皮肤黏膜层微生物α多样性分析
注：*. p<0.05；**. p<0.01；***. p<0.001。
Figure 4. Alpha diversity analysis of skin mucosa from healthy and diseased A. sinensis
Note: *. p<0.05; **. p<0.01; ***. p<0.001.
![]()
图 5 健康和患病中华鲟皮肤黏膜层微生物PCoA图
Figure 5. PCoA of skin mucosa from healthy and diseased A. sinensis
-
[1] 张艳敏, 杨国坤, 李克克, 等. 鱼类黏膜层微生物研究进展[J]. 水产学报, 2022, 46(6): 1117-1127. [2] 翟万营. 南极鱼亚目鱼类肠道、皮肤黏液及其生存水环境微生物结构和功能研究[D]. 上海: 上海海洋大学, 2020: 7-12. [3] BLETA M C, KELLY M, SABINO-PINTO J, et al. Disruption of skin microbiota contributes to salamander disease[J]. Proc Biol Sci, 2018, 285(1885): 20180758.
[4] MENG K F, DING L G, WU S, et al. Interactions between commensal microbiota and mucosal immunity in teleost fish during viral infection with SVCV[J]. Front Immunol, 2021, 12: 654758. doi: 10.3389/fimmu.2021.654758
[5] 詹梦婷. 虹鳟感染IHNV后皮肤黏膜免疫及其表面微生物变化的研究[D]. 武汉: 华中农业大学, 2023: 47-48. [6] 许峻旗, 谢建平, 王志坚. 大鲵皮肤溃疡死亡相关的微生物区系[J]. 水产学报, 2022, 46(11): 2186-2195. [7] 许兰馨, 周亮, 危起伟. 长江鲟类资源现状及保护[J]. 水产学报, 2023, 47(2): 62-72. [8] 田甜, 张德志, 杜合军. 中华鲟主要细菌性疾病及其临床诊断与防控[J]. 水产科学, 2018, 37(6): 800-805. [9] 张书环, 聂品, 舒少武, 等. 子二代中华鲟分枝杆菌感染及血液生理生化指标的变化[J]. 中国水产科学, 2017, 24(1): 136-145. [10] DI J, ZHANG S H, HUANG J, et al. Isolation and identification of pathogens causing haemorrhagic septicaemia in cultured Chinese sturgeon (Acipenser sinensis)[J]. Aquac Res, 2018, 49(11): 3624-3633. doi: 10.1111/are.13830
[11] 邸军, 张书环, 黄君, 等. 中华鲟脑膜败血伊丽莎白菌的分离鉴定及药敏特性[J]. 水产学报, 2018, 42(1): 120-130. [12] 田甜, 张建明, 杜合军. 中华鲟源维氏气单胞菌的分离鉴定及其药敏特性[J]. 安徽农业大学学报, 2017, 44(6): 1016-1017. [13] BAY L, RING H C. Human skin microbiota in health and disease[J]. APMIS, 2021, 130(12): 706-718.
[14] 文海燕, 胡晓, 周鑫, 等. 高通量测序技术在呼吸道病原检测中的应用[J]. 口岸卫生控制, 2024, 29(2): 29-32. doi: 10.3969/j.issn.1008-5777.2024.02.008 [15] CAMARA-RUIZ M, CEREZO I M, GUARDIOLA F A, et al. Alteration of the immune response and the microbiota of the skin during a natural infection by Vibrio harveyi in European seabass (Dicentrarchus labrax)[J]. Microorganisms, 2021, 9: 964. doi: 10.3390/microorganisms9050964
[16] 张雪萍, 徐莉, 符宏高, 等. 基于16sRNA扩增子测序分析患烂皮病棘胸蛙皮肤的细菌多样性[J]. 福建农业科技, 2022, 53(11): 1-9. [17] HUANG H Y, ZHOU P J, CHEN P, et al. Alteration of the gut microbiome and immune factors of grass carp infected with Aeromonas veronii and screening of an antagonistic bacterial strain (Streptomyces flavotricini)[J]. Microb Pathogenesis, 2020, 143: 104092. doi: 10.1016/j.micpath.2020.104092
[18] LI W F, YANG H Y, ZHAO Q, et al. Polyphenol-rich loquat fruit extract prevents fructose-induced nonalcoholic fatty liver disease by modulating glycometabolism, lipometabolism, oxidative stress, inflammation, the intestinal barrier, and the gut microbiota in mice[J]. J Agr Food Chem, 2019, 67(27): 7726-7737. doi: 10.1021/acs.jafc.9b02523
[19] PARK J, ROH H, LEE Y, et al. Characterization and pathogenicity of Flavobacterium psychrophilum isolated from rainbow trout (Oncorhynchus mykiss) in Korea[J]. Microorganisms, 2023, 11(10): 2546. doi: 10.3390/microorganisms11102546
[20] BRUCE T J, MA J, JONES E M, et al. Assessment of Flavobac terium psychrophilum-associated mortality in Atlantic salmon (Salmo salar) and brook trout (Salvelinus fontinalis)[J]. J Fish Dis, 2021, 44(5): 645-653. doi: 10.1111/jfd.13349
[21] CHINCHILLA B, VAZQUEZ-FERNANDEZ E, REBOLLADA-MERINO A, et al. First detection of Flavobacterium psychrophilum in juvenile Siberian sturgeons (Acipenser baerii) and description of the pathological findings[J]. J Fish Dis, 2023, 46(8): 887-894. doi: 10.1111/jfd.13801
[22] 许莹萱. 罗非鱼柱状黄杆菌感染的组织载菌量和鳃转录组分析[D]. 海口: 海南大学, 2023: 26-27. [23] 柴静茹, 王荻, 卢彤岩, 等. 虹鳟源嗜冷黄杆菌的分离鉴定及致病性[J]. 水产学报, 2021, 45(6): 929-939. [24] 刘丽娟, 王红莉, 王亚军, 等. 笋壳鱼源柱状黄杆菌的分离鉴定及致病性和药物敏感性[J]. 南方农业学报, 2023, 54(2): 638-646. doi: 10.3969/j.issn.2095-1191.2023.02.032 [25] 田甜, 张建明, 张德志, 等. 健康和患病中华鲟皮肤黏膜层微生物菌群结构特征[J]. 中国水产科学, 2023, 30(10): 1259-1270. [26] 徐杨, 江昀, 侯金亮, 等. 发育不良和歪头病对变态末期棘胸蛙皮肤菌群的影响[J]. 水生生物学报, 2024, 48(5): 839-848. doi: 10.7541/2024.2023.0303 [27] 张雪晨. 健康和患病对虾肠道细菌组成比较及不同生长阶段溯源分析[D]. 宁波: 宁波大学, 2021: 18-20. [28] WANG C, SUN G X, LI S S, et al. Intestinal microbiota of healthy and unhealthy Atlantic salmon Salmo salar L. in a recirculating aquaculture system[J]. J Oceanol Limnol, 2018, 36(2): 414-426. doi: 10.1007/s00343-017-6203-5
[29] 柴静茹, 王荻, 卢彤岩, 等. 嗜冷黄杆菌及细菌性冷水病的研究进展[J]. 大连海洋大学学报, 2020, 35(5): 755-761. [30] 张玉蕾. 柱状黄杆菌病原特性与耐药机制研究[D]. 武汉: 华中农业大学, 2017: 1-8. -
期刊类型引用(2)
1. 朱海晨，汤建华，吴磊，施金金，王储庆，王燕平，闫欣，葛慧. 黄海南部海域鱼类群落结构及生物量粒径谱特征. 中国水产科学. 2023(09): 1102-1111 . 
百度学术
2. 王寅，赵静，高春霞，王士聪，叶深. 浙江南部近海春季鱼类群落结构特征及其年际变化. 海洋与湖沼. 2022(06): 1455-1466 . 百度学术
其他类型引用(2)
下载:
计量
- 文章访问数: 111
- HTML全文浏览量: 28
- PDF下载量: 56
- 被引次数: 4
粤公网安备 44010502001741号
