Growth, maturity and mortality of Johnius distinctus and J. belangerii in offshore waters of southern Zhejiang Province
-
摘要: 鳞鳍叫姑鱼 (Johnius distinctus) 和皮氏叫姑鱼 (J. belangerii) 是我国东海的小型经济鱼类,且为生态系统中重要的饵料鱼类,研究分析其生物学特征对促进东海生态系统和渔业的科学管理具有重要意义。根据2015年11月至2018年11月浙江南部近海渔业资源季度调查数据,运用线性混合效应模型 (Linear mixed effects model, LMEM) 研究了这2种叫姑鱼的生长特征及其异质性;通过ELEFAN法估算了2种叫姑鱼的极限体长L∞、趋近极限体长的相对生长速率k、理论生长起始年龄t0、自然死亡系数M、总死亡系数Z、开发率E以及50%性成熟体长。结果表明,季节对2种叫姑鱼体长和体质量关系影响最为显著;鳞鳍叫姑鱼L∞=216.83 mm,k=0.84,t0=−0.18;皮氏叫姑鱼L∞=207.38 mm,k=0.63,t0=−0.25;鳞鳍叫姑鱼Z=3.80,M=1.53,E=0.60;皮氏叫姑鱼Z=1.88,M=1.24,E=0.34;鳞鳍叫姑鱼和皮氏叫姑鱼50%性成熟体长分别为133.1和128.9 mm。Abstract: Karut croaker (Johnius distinctus) and Belanger's croaker (J. belangerii) are common small economic fishes and important prey species in the ecosystem of the East China Sea. Studies on their biological characteristics are essential to understand the ecosystem and to improve the scientific management of these important fisheries in the East China Sea. Based on the seasonal surveys of fishery resources in the offshore waters of southern Zhejiang Province, we applied linear mixed effects model (LMEM) to study the growth characteristics and heterogeneity of the two species. The asymptotic body length (L∞), growth rate (k) and theoretical growth starting age (t0) were estimated by ELEFAN method. The natural mortality coefficient (M), total mortality coefficient (Z), development rate (E), and the 50% sexual mature body length of the two species were estimated. The LMEM results show that seasons have the most significant effect on their growth. For J. distinctus, L∞=216.83 mm, k=0.84, t0=−0.18; for J. belangerii, L∞=207.38 mm, k=0.63, t0=−0.25. For J. distinctus, Z=3.80, M=1.53, E=0.60; for J. belangerii, Z=1.88, M=1.24, E=0.34. The length of 50% maturity's estimates were 133.1 mm and 128.9 mm for J. distinctus and J. belangerii, respectively.
-
Keywords:
- Johnius distinctus /
- J. belangerii /
- Linear mixed effects model /
- Sex maturity /
- East China Sea
-
三疣梭子蟹 (Portunus trituberculatus) 隶属甲壳纲、十足目、梭子蟹科、梭子蟹属,广泛分布于西北太平洋及东南亚沿海,处于食物链较底层,生长速度快、适应能力强、分布范围广,是我国近海海洋生态系统中的优势物种,是海洋生态系统中物质传递和能量流动的重要载体,决定着海洋生态系统的许多重要生态过程。近年来,有关三疣梭子蟹的研究大多集中在生物学特性与养殖[1-5]、基因表达与免疫[6-9]等方面,对其年龄鉴定的相关研究甚少。但年龄鉴定一直是甲壳类生物研究的重要内容,有利于了解及掌握其种群结构、开展资源评估,为制订渔业管理措施提供依据,因此提供更加精确的直接年龄鉴定方法十分必要。
甲壳类生物因其特殊的生长方式 (蜕壳),其年龄鉴定尤其困难。长期以来大多数学者采用标志重捕法[10]和体长频率法[11]等间接方法鉴定其年龄,但因方法的局限性,鉴龄的精确度不高;放射性同位素分析法[12]虽可较为精准鉴龄,但成本较高。最新研究表明甲壳类硬组织结构 (眼柄和胃磨) 能够作为直接鉴定其年龄的方法[13]。Leland等[14-15]利用甲壳动物的胃磨研磨截面微结构进行年龄鉴定,为甲壳动物年龄鉴定研究提供了新方向。Kilada等[13,16]先后发现美国螯龙虾 (Homarus americanus) 与斑纹黄道蟹 (Cancer irroratus) 两者的胃磨生长纹与年龄有着一定关系。现有研究中对三疣梭子蟹进行年龄鉴定的较少,蒋瑞等[17]发现三疣梭子蟹眼柄生长纹与其生长有关,但未肯定宽带个数代表其年龄。本研究通过对三疣梭子蟹胃磨微结构进行观察,探讨利用胃磨微结构对三疣梭子蟹进行年龄鉴定的可行性,分析胃磨微结构与基础生物学数据及年龄间的关系,以期为今后深入开展三疣梭子蟹年龄鉴定和其资源评估提供参考依据。
1. 材料与方法
1.1 实验样本的采集
三疣梭子蟹样本共185只,分别采集于三疣梭子蟹南通养殖基地 (100只)、江苏南通近海 (45只) 和上海芦潮港附近 (40只)。
1.2 胃磨的制片
测量取样个体生物学数据 [甲长 (Carapace length,CL)、甲宽 (Carapace width,CW)、体质量 (Body mass,m) ] 后取出胃磨,放置于70%的乙醇溶液中运送回实验室,在实验室中将胃磨转移至乙醇甘油混合溶液中 [V (乙醇)∶V (甘油)∶V (水)=6∶3∶1];实验时使用解剖针剔除胃磨中附着的有机物,并使用蒸馏水冲洗干净。采用国内外胃磨结构划分的常规方法[18],可将胃磨结构分为5个部分,即前幽门骨、外幽门骨、尾贲门骨、轭贲门骨和翼贲门骨 (图1)。以往的甲壳类鉴龄研究表明,外幽门骨和尾贲门骨经过预处理并切片后能在显微镜下清晰地观察到类似鱼类耳石或甲壳类眼柄中生长纹结构[15],因此本研究采用这两种结构进行切片观察。
采用蒋瑞等[19]眼柄的制片方法,对于胃磨中外幽门骨和尾贲门骨在进行包埋前先使用无水乙醇脱水,再放入包埋盒中使用调配好的冷埋树脂溶液固定 (固化剂与亚克力粉按照体积比1.5∶1的比例调配,固化剂黏度为60 mPas)。包埋好的树脂样本放在阴凉通风处24 h待其硬化,后将包有胃磨的模具依次用120、240、600、1 200和2 500目的水磨砂纸打磨,在打磨过程中不断在显微镜下观察,直至能清晰观察到胃磨结构的平面,之后对此树脂薄片抛光,再清除其表面杂质,用吸水纸擦干表面水分。然后将抛光过的截面朝下,用热熔胶粘在载玻片上,使用同样的方法对另一面进行打磨,直至在显微镜下模具截面的生长纹清晰可见,再抛光得到切片。最后将做好的切片用密封袋密封保存,以备观察微结构时使用。
1.3 生长纹读取方法
生长纹的计数以胃磨中“暗纹”为计数标准。使用目视方法确定角质层边界,并且以该边界为计数起点。将胃磨切片照片依据照片的可读性分为3个等级。1) 清晰:沿着单个平面读数,可以连续且快速地读取所有生长纹的切片;2) 基本清晰:无法在单个平面内连续读取所有生长纹,但是通过平面的平移可以将生长纹全部读取的切片;3) 无法读取:无法辨别生长纹分界、明显脱钙、生长纹不完全及有明显缺陷的切片,这切片的读数不计入后续的数据分析。
每个样本的胃磨切片照片分别由3名不同实验室人员在相互独立的环境下判读其轮纹数量。采用变异系数 (CV) 以校对每个样本读数的准确性,公式为:
$$ {\rm{CV}} = \frac{{\sqrt {\frac{{{{\left( {{R_1} - R} \right)}^2} + {{\left( {{R_2} - R} \right)}^2} + {{\left( {{R_3} - R} \right)}^2}}}{2}} }}{R} \times 100{\text{%}} $$ (1) 其中R1、R2、R3分别表示3名不同人员对同一样本的读数,R表示3位人员读数的平均值,根据其他研究的经验[20-21],当每个样本读数的变异系数值不高于10%时认为该次读数可信。
1.4 蜕壳次数的估算
通过周期性的蜕壳实现生长是甲壳动物最明显的生长特点。掌握甲壳类动物的蜕壳规律是研究其年龄与生长等信息的重要依据。高天龙等[22]认为三疣梭子蟹蜕壳周期时间随着期数的增加而增大,但其个体生命史中只存在有限数量的蜕壳次数,同时引入了单位阶梯函数来表示蜕壳次数f(m) 与体质量m的关系,具体函数为:
当m <0,u(m) =0;当m >0,u(m) =1;
$$ \begin{array}{c} f\left( m \right) = u\left( m \right) + u\left( {m - 0.014} \right) + u\left( {m - 0.033} \right) + u( {m }-\\ 0.104 ) + u\left( {m - 0.411} \right) + u\left( {m - 1.4} \right) + u\left( {m - 3.89} \right) +\\ u\left( {m - 11.56} \right) +u\left( {m - 32.77} \right) + u\left( {m - 66.28} \right) + \\ u\left( {m - 114.02} \right) + u\left( {m - 177.84} \right) + u\left( {m - 256.61} \right) \end{array} $$ (2) 该阶梯函数对于蜕壳阶段在5—11期的三疣梭子蟹判对率高达91.38%,因此本文采用该阶梯函数来推算取样三疣梭子蟹的蜕壳次数,并以此来探索蜕壳次数与胃磨生长纹数量之间的关系。
1.5 生长方程模型
由于缺乏对其的精准鉴龄,有关其生长方程的研究也较少。故本研究采用目前推算水生生物常用的7种生长方程模型:线性模型、幂函数模型、指数模型、对数模型、Logistic模型、von Bertalanffy模型 (v-B模型) 和Gompertz模型来拟合三疣梭子蟹的生长方程,并开展最适模型的选择,各模型公式为:
线性方程,
$$ f\left( t \right) = {\rm{a}} + {\rm{b}}t $$ (3) 幂函数方程,
$$ f\left( t \right) = {\rm{a}}{t^{\rm{b}}} $$ (4) 指数方程,
$$ f\left( t \right) = {\rm{a}}{{\rm{e}}^{{\rm{b}}t}} $$ (5) 对数函数方程,
$$ f\left( t \right) = {\rm{a}}\ln t + {\rm{b}} $$ (6) Logistic生长方程,
$$ f\left( t \right) = \frac{{{{{L}}_\infty }}}{{1 + {{\rm{e}}^{ - {\rm{K}}\left( {t - {{\rm{t}}_0}} \right)}}}} $$ (7) von Bertalanffy生长方程,
$$ f\left( t \right) = {L_\infty } \times \left[ {1 - {{\rm{e}}^{ - {\rm{K}}\left( {t - {t_0}} \right)}}} \right] $$ (8) Gompertz生长方程,
$$ f\left( t \right) = {L_\infty } \times {{\rm{e}}^{\left[ {1 - {{\rm{e}}^{ - {\rm{K}}\left( {t - {t_0}} \right)}}} \right]}} $$ (9) 其中f (t) 表示三疣梭子蟹的甲长 (CL)、甲宽 (CW) 或体质量 (m),t表示胃磨生长纹个数,a、b、K、L∞和t0表示对应公式中常数项。对于Logistic方程、v-B方程和Gompertz方程,L∞表示无限年龄的三疣梭子蟹渐进的甲长、甲宽或体质量,t0表示理论上当甲长、甲宽或体质量为0时三疣梭子蟹的生长纹数量,在一般情况下,该值为负值。K值则表示常数。为了更加快速地拟合各公式,本研究中在迭代过程中起始值默认为所有样本中甲长、甲宽或体质量的最大值,而为了避免出现理论上的错误,t0和K的起始值则为0,并约束−t0和K值为正值。
在选取最适函数时,使用赤池信息量准则(Akaike information criterion, AIC),AIC为最小时,则认为该模型在常用7个生长模型当中为最适模型,选定最适函数后使用回归系数 (R2) 作为判断数据对该模型的拟合程度。
AIC计算方法为:
$$ {\rm{AIC}} = 2k + n{\rm{ln}}\left( {\frac{{{\rm{RSS}}}}{n}} \right) $$ (10) 式中k为模型的参数数量,n为样本数量,RSS为所有样本量的残差平方和。
2. 结果
2.1 胃磨的微结构
胃磨位于虾蟹类的消化系统中,大部分虾蟹类食道通入一个能够收缩的三角形囊状胃,胃内部角质层发达,与胃壁共同形成的嵴构成胃磨。由于胃磨形状的不规则性及胃磨由多块相对独立的部分组合而成,因此将胃磨在显微镜下所观察的微结构分为两部分:1) 在肉眼下能观察到类似与动物牙齿的齿板部分;2) 不包含齿板的剩余白色部分 (柄部及翼部)。生长纹主要分布于胃磨各部分的柄部和翼端,本研究中胃磨微结构特指能够观察到生长纹的柄部切面微结构。
通过显微镜下观察三疣梭子蟹胃磨微结构后发现 (图2),胃磨微结构由4个部分组成:蜡质层、外角质层、钙化层和膜层。蜡质层紧邻外角质层,外角质层位于胃磨靠向外一侧,边缘规则且光滑,但因内部含有少量色素透光度低,但随着切片厚度减小,透光度有所提升,主要由几丁质构成,在显微镜下观察不到任何的周期性生长纹结构。膜层在胃磨内侧与肌肉或其他组织接触部分,在显微镜下可以观察到黑色部分,透光度较高,钙化层是胃磨微结构中占比最大的一部分,钙化程度较高,在显微镜下能观察明暗相间的生长纹结构,与膜层为平移关系,同时钙化层中还存在垂直生长纹。
与同样作为甲壳类直接年龄鉴定材料的眼柄微结构 (图2-c) 相比,胃磨微结构与其略有不同。眼柄切片中色素层充满了红色色素,而胃磨中表层较薄,色素层呈现黑色,且色素层占比比眼柄切片小。胃磨钙化层所占比例要大于眼柄中钙化层。在胃磨和眼柄切片中均可观察到密集且间距均匀的细纹,还可以观察到颜色较深、排列稀疏的宽带,一个宽带内可以看到包含数条或十数条细纹。
2.2 生长纹读数的精确性
生长纹读数的准确度会直接影响到最终的实验结果,因此在利用生长纹个数拟合生长方程之前,要判断对样本生长纹读数的精确度。Kilada等[23]利用眼柄和胃磨直接鉴定南极磷虾 (Euphausis superba) 年龄发现,当变异系数CV小于10%,该组数据可被使用。在本研究中,由3名实验人员单独对同一样本的胃磨生长纹读数,经计算得出平均变异系数为8.6% (n=185),因此该生长纹数据准确可靠,可用于拟合有关生长方程的研究。
2.3 生长模型建立
2.3.1 甲长、甲宽与体质量生长模型
通过基本生物学测量分析后,185只样本中包含雌蟹120只,雄蟹65只。总体样本中甲长为20~99 mm [平均值 (57.27±1.30) mm],甲宽介于43~190 mm [平均值 (116.25±2.22) mm],体质量介于3.70~365.89 g [平均值 (105.86±5.92) g]。雌雄三疣梭子蟹基础生物学数据见表1,在平均值上雌蟹数据要略高于雄蟹。分别对所有的三疣梭子蟹甲长、甲宽和体质量数据进行分析,得出以下三疣梭子蟹甲长CL与甲宽CW的关系 (图3、图4、图5)。
表 1 三疣梭子蟹样本个体基础生物学数据Table 1. Basic biological data of gastric abrasion growth samples of P. trituberculatus项目
Item雄蟹 Male 雌蟹 Female 甲长 CL/mm 甲宽 CW/mm 体质量 m/g 甲长 CL/mm 甲宽 CW/mm 体质量 m/g 样本量 Sample number 65 65 65 120 120 120 最大值 Maximum 99 190 365.89 97 189 363.17 最小值 Minimum 20 43 3.7 24 60 7.29 平均值 Mean 50.05 105.47 85.33 61.18 122.08 166.97 标准误 Standard error 2.27 4.17 10.16 1.47 2.67 7.10 雄性三疣梭子蟹:
$$ {\rm{CW}} = 1.77{\rm{CL}} + 16.94\;(n = 65,{R^2} = 0.93,P < 0.01) $$ (11) 雌性三疣梭子蟹:
$$ {\rm{CW}} = 1.72{\rm{CL}} + 16.93\;(n = 120,{R^2} = 0.89,P < 0.01) $$ (12) 185只三疣梭子蟹:
$$ {\rm{CW}} = 1.72{\rm{CL}} + 17.93\;(n = 185,{R^2} = 0.89,P < 0.01) $$ (13) 三疣梭子蟹甲宽CW与体质量m的关系:
雄性三疣梭子蟹:
$$ m = 8.28 \times {10^{ - 5}}{\rm{C}}{{\rm{W}}^{2.91}}(n = 65,{R^2} = 0.87,P < 0.01) $$ (14) 雌性三疣梭子蟹:
$$ m = 4.92 \times {10^{ - 5}}{\rm{C}}{{\rm{W}}^{3.02}}(n = 120,{R^2} = 0.92,P < 0.01) $$ (15) 185只三疣梭子蟹:
$$ m = 6.31 \times {10^{ - 5}}{\rm{C}}{{\rm{W}}^{2.97}}(n = 185,{R^2} = 0.92,P < 0.01) $$ (16) 2.3.2 甲长、甲宽、体质量与生长纹的关系
根据最小AIC法原则 (表2),185个样本中,甲长、甲宽与胃磨生长纹个数的关系均适合v-B生长模型来表示 (图6),关系式分别为:
表 2 7个生长模型函数拟合结果Table 2. Seven growth model function fitting results模型
Model甲长-生长纹个数
CL-t甲宽-生长纹个数
CW-t体质量-生长纹个数
m-tR2 AIC R2 AIC R2 AIC 线性 Linear 0.884 473.77 0.922 539.78 0.868 710.89 幂函数 Power 0.834 504.01 0.886 571.06 0.910 678.42 指数 Exponent 0.865 486.36 0.906 555.32 0.913 675.90 对数 Logarithm 0.841 500.07 0.897 562.54 0.698 780.75 逻辑斯蒂 Logistic 0.897 465.69 0.953 499.09 0.907 683.98 v-B 0.910 449.19 0.995 4 498.26 0.914 677.15 冈珀茨 Gompertz 0.870 484.63 0.949 504.56 0.922 668.93 注:加粗表示拟合结果最好 Note: Bold indicates the best fit result. $$ {\rm{CL}}=103.42 \times ( {1-{{\rm{e}}^{ - 0.027\;7t-1.149}}} )( {n=185,{R^2}=0.91,P < 0.01} ) $$ (17) $$ {\rm{CW}} = 202.96 \times ( {1 - {{\rm{e}}^{ - 0.032\;2t}}} )(n = 185,{R^2} = 0.90,P < 0.01) $$ (18) 体质量与胃磨生长纹个数的关系适合使用Gompertz模型来表示 (图6),且关系式为:
$$ m = 356.27 \times {{\rm{e}}^{\left[ {1 - {{\rm{e}}^{ - 0.027\;1\left( {t - 59.6} \right)}}} \right]}}\left( {n = 185,{R^2} = 0.92,P < 0.01} \right) $$ (19) 同样,雄性三疣梭子蟹甲长、甲宽、体质量与其胃磨生长纹数目的关系均适于用Gompertz模型描述,表达式分别为:
$$ {\rm{CL}} = 51.296 \times {{\rm{e}}^{\left[ {1 - {{\rm{e}}^{0.028\left( {t - 26.014} \right)}}} \right]}}( {n = 65,{R^2} = 0.98,{{P}} < 0.01} ) $$ (20) $$ {\rm{CW}} = 87.853 \times {{\rm{e}}^{\left[ {1 - {{\rm{e}}^{ - 0.033\left( {t - 18.471} \right)}}} \right]}}( {n = 65,{R^2} = 0.99,{{P }}< 0.01} ) $$ (21) $$ m = 2{\rm{}}512.410 \times {{\rm{e}}^{\left[ {1 - {{\rm{e}}^{0.015\left( {t - 129.841} \right)}}} \right]}}( {n = 65,{R^2} = 0.96,{{P}}< 0.01} ) $$ (22) 对于雌性三疣梭子蟹,其甲宽以及体质量与胃磨生长纹数目的关系均适于用幂函数生长模型来描述,表达式分别为:
$$ {\rm{CW}} = 17.704{t^{0.573}}(n = 120,{R^2} = 0.89,P < 0.01) $$ (23) $$ m = 0.357{t^{1.684}}(n = 120,{R^2} = 0.81,P < 0.01) $$ (24) 雌性甲长与胃磨生长纹数目的关系则适于用v-B生长模型描述,表达式为:
$$ {\rm{CL}}=98.007\times\left[ {1-{{\rm{e}}^{ - 0.037{\rm{}}\left( {t - 1.201} \right){\rm{}}}}} \right](n=120,{R^2}=0.84,P < 0.01) $$ (25) 2.4 生长纹与预估蜕壳次数关系
在本研究中选取70个胃磨切片清晰的样本,利用单位阶梯函数的判定,得到其样本的预估蜕壳次数。通过对胃磨生长纹和预估蜕壳次数的拟合,两者存在幂函数关系 (图7),关系式为u=0.207 1t2.105 2(n=70,R2=0.77)。
3. 讨论
Leland等[14]在研究红鳌鳌虾 (Cherax quadricarinatus) 胃磨时发现胃磨研磨截面展示出与鱼类耳石相近的同心生长纹,这开启了虾蟹类年龄鉴定研究的新方向,而后又使用钙黄绿素染色追踪红鳌鳌虾胃磨生长后,在胃磨切片中能明显看到生长痕迹,这也是目前证明胃磨能够作为甲壳类直接年龄鉴定最有力的证据[15]。Sheridan等[24]和Kilada等[13]在实验中采用外幽门骨或尾贲门骨作为甲壳类年龄鉴定,结果显示两者组织切片在显微镜下能够观察到清晰的生长纹,并且其生长纹与年龄信息显著相关。Sheridan等[24]利用硬组织生长纹计数法对普通黄道蟹 (Cancer pagurus) 鉴龄时,其轭贲门骨中生长纹清晰可见,但未在眼柄中观察到清晰的生长纹,这不排除是样本制备所导致。本研究结果显示,胃磨微结构中存在与眼柄相似的生长纹,且胃磨微结构中钙化层比例较大,生长纹清晰可见。对于甲壳类动物来说,目前为止没有哪一种方法可以独立对其进行精确鉴龄,因此可以同时对眼柄及胃磨生长纹分析,两者互为辅助,并与多种鉴龄方法相结合,交叉验证,以提高鉴龄的精确性。对甲壳类年龄鉴定最大的问题就是所使用的材料是否会在甲壳类独特的生长过程 (周期性蜕壳) 中被蜕去。如果在蜕壳过程中所使用的年龄鉴定材料随着老化的外壳丢失,那么该鉴龄材料则只能记录两次蜕壳间的生长信息。但是,Kilada等[13]在对雪蟹 (Chionoecetes opilio) 蜕壳的观察中发现,其眼柄不会随着蜕壳过程丢失,这使得眼柄可以成为甲壳类鉴龄材料。同样地,Leland等[15]也证明了在水晶蟹 (Chaceon albus) 蜕壳过程中,其胃磨也不会随着蜕壳过程而丢失。Krafft等[25]研究发现,胃磨当中的生长纹与已知能够作为甲壳类年龄鉴定材料的眼柄生长纹数量呈函数关系。因此胃磨中生长纹可以作为年龄鉴定材料。三疣梭子蟹通过蜕壳实现生长,在整个生长周期内,蜕壳次数有限。王小刚[26]认为三疣梭子蟹的龄期 (每蜕壳1次,称作1个龄期) 与其蜕壳次数呈正相关关系。故本研究引入高天龙等[22]建立的通过体质量判断期龄的判定函数,以此来探索蜕壳次数与胃磨生长纹数量之间的关系。结果发现,蜕壳次数与生长纹数目之间呈幂函数关系。即在幼蟹时期,三疣梭子蟹每次蜕壳间胃磨所累积的生长纹数目要远少于成年后三疣梭子蟹每次蜕壳间胃磨所累积的生长纹数目。
国内外渔业生物学研究者在对甲壳类生物研究中,各甲壳类甲宽 (头胸甲长) 与体质量关系的研究均采用幂函数,其中幂函数系数b均介于2.5~3.5。通过对所采集的样本分析,本研究中总样本的幂指数系数为3.05,雄蟹为2.91,雌蟹为3.02。即可认为本研究所使用样本符合大部分三疣梭子蟹的生长规律,即等速生长,说明三疣梭子蟹个体在整体生活史中体型随着质量增加等比增大。同时本研究中雌雄三疣梭子蟹幂指数系数相差不大,说明三疣梭子蟹甲宽与体质量关系函数在性别上没有显著差异。与王小刚等[27]使用v-B生长模型所计算得出的渐近甲宽 (L∞=210 mm) 相比,本研究中通过胃磨生长纹代替实际年龄所拟合出的函数其渐近甲宽 (L∞=202.96 mm) 要小。同时从以往的研究来看,在20世纪后期林景祺[28]在中国海洋渔业资源中指出中国近海三疣梭子蟹甲宽范围在30~250 mm,优势甲宽为120~190 mm;张健等[29]研究表明,在中国东海近海的三疣梭子蟹甲宽范围在63.3~201.5 mm,优势甲宽仅为110~160 mm。说明中国近海部分三疣梭子蟹偏向于小型化,且相对于早期个体有明显差距,这也是本研究中渐近甲宽、甲长和体质量要略小于早几年研究的原因。三疣梭子蟹甲宽与甲长之间为线性关系,这也和大多数学者研究结果一致。此外,逻辑斯蒂生长方程、v-B生长模型和Gompertz方程均可用于描述水生生物生长,而在本研究中对于三疣梭子蟹胃磨生长纹个数与其基础生物学数据的关系中,甲长与甲宽和胃磨生长纹数目关系适于用v-B生长模型描述,而体质量和胃磨生长纹数目则更适于用Gompertz模型进行描述。Kilada等[30]在雌性欧洲龙虾 (Homarus gammarus) 和雌性斑纹黄道蟹 (Cancer irroratus) 的轭贲门骨切片中分别观察到了4个和9个宽生长纹,且证实宽带个数与两者实际年龄相近。Kilada和Acuna[16]也在黄蹲龙虾 (Curvimunida johni) 的中位骨切片中观察到了宽带,且认为此宽带是以年为单位形成的。Ismail等[31]研究得出榄绿青蟹 (Scylla oliviacea) 轭贲门骨切片中的宽带数目可以作为估计其年龄的依据。本研究在三疣梭子蟹外贲门骨切片中也观察到了宽带,依据前人研究结果,推测可以用宽带数目估计其年龄。
本研究表明,采用胃磨对三疣梭子蟹鉴龄是可行的,但本实验采用的均为未知年龄的三疣梭子蟹样本,因此只能依据前人研究成果推测胃磨切片中宽带个数与实际年龄相近,后续应使用已知年龄的样本对该推测进行验证。在探讨三疣梭子蟹胃磨生长纹与其年龄的相关性时,笔者尚无法对胃磨生长纹的形成机制作出合理解释。未来的研究重点应落在以下两个方向:1) 延长实验样本 (已知年龄) 采集周期并缩短采集间隔,以验证更加精确的周期性实验,提高三疣梭子蟹胃磨鉴龄的精确性;2) 采用材料学方法 (电子显微镜、荧光染色和同位素分析等) 探究三疣梭子蟹胃磨的形成机制。
-
表 1 2015—2018年浙江南部近海鳞鳍叫姑鱼和皮氏叫姑鱼的样本数量
Table 1 Sample size of J. distinctus and J. belangerii in offshore waters of southern Zhejiang Province from 2015 to 2018
种类
Species年份
Year春季
Spring夏季
Summer秋季
Autumn冬季
Winter总计
Total鳞鳍叫姑鱼
J. distinctus2015 — — 125 — 125 2016 54 123 160 1 338 2017 123 21 162 1 307 2018 120 28 99 3 250 总计 Total 297 172 546 5 1 020 皮氏叫姑鱼
J. belangerii2015 — — 69 — 69 2016 124 54 124 17 319 2017 37 33 150 35 255 2018 21 13 101 96 231 总计 Total 182 100 444 148 874 表 2 浙江南部近海鳞鳍叫姑鱼和皮氏叫姑鱼体长体质量关系的拟合模型及其评价指标
Table 2 Candidate models to describe length-mass relationship for J. distinctus and J. belangerii in offshore waters of southern Zhejiang Province and their AIC and RMSE values
模型编号
Model ID截距的随机效应
Random effects for loga斜率的随机效应
Random effects for b鳞鳍叫姑鱼 J. distinctus 皮氏叫姑鱼 J. belangerii 赤池信息量准则
AIC均方根误差
RMSE/%赤池信息量准则
AIC均方根误差
RMSE/%1 — — −759 16.60 −298 20.33 2 性别 — −757 16.53 −302 20.34 3 — 性别 −757 16.52 −301 20.34 4 性别 性别 −795 16.85 −298 20.33 5 年份 — −777 16.33 −298 20.26 6 — 年份 −775 16.36 −298 20.26 7 年份 年份 −812 15.96 −294 20.23 8 季节 — −848 15.75 −320 19.91 9 — 季节 −848 15.75 −320 19.91 10 季节 季节 −844 15.75 −317 19.86 11 性别、年份 — −775 16.34 −301 20.34 12 — 性别、年份 −773 16.32 −301 20.33 13 性别、年份 性别、年份 −813 16.24 −295 20.27 14 季节、性别 — −846 16.20 −322 19.54 15 — 季节、性别 −846 16.20 −322 19.54 16 季节、性别 季节、性别 −841 16.19 −316 19.63 17 季节、年份 — −860 15.58 −318 19.91 18 — 季节、年份 −860 15.58 −318 19.91 19 季节、年份 季节、年份 −868 15.39 −312 19.80 20 季节、性别、年份 — −858 15.42 −320 19.61 21 — 季节、性别、年份 −858 15.43 −320 19.61 22 季节、性别、年份 季节、性别、年份 −862 15.48 −311 19.59 表 3 不同海域皮氏叫姑鱼生长、死亡参数
Table 3 Growth and mortality parameters of J. belangerii in different China's seas
生长条件
因子
a异速生长
因子
b极限
体长
L∞相对
生长速率
k生长起始
年龄
t0总死亡
系数
Z自然死亡
系数
M捕捞死亡
系数
F开发率
E马鞍列岛[17] Ma'an Archipelago 1.35×10–5 3.09 福建沿岸[38] Fujian Coastal 8.47×10–6 3.2 200 0.56 –0.29# 2.01 0.63 1.38 0.69 台湾西南部*[20] Southwestern Coast of Taiwan 4.44×10–5 2.83 闽江口*[21] Min River Estuary 1.02×10–5 3.16 厦门海域[19] Xiamen Sea Area 3.44×10–5 3.08 220.5 0.56 –0.61 2.25 0.63 1.62 0.72 厦门湾[28] Xiamen Bay 1.14×10–5 3.11 210 0.81 –0.19 1.09 0.61 0.48 0.44 流沙湾*[18] Liusha Gulf 9.02×10–6 3.17 珠江口[23] Pearl River Estuary 3.44×10–5 2.88 北部湾*[22] Beibu Gulf 1.67×10–5 3.05 注:*. 原文计算单位为cm,本研究进行对比时a已按mm换算;#. 原文中无t0,本研究根据公式 (4) 计算得出t0。 Note: *. The unit in the references is cm, and the value of a has been converted in this study; #. There is no t0 estimate in the references, and t0 was estimated according to Formula (4) in this study. -
[1] FISHBASE. Johnius distinctus [DB/OL]. [2021-01-03]. www.fishbase.se/summary/Johnius-distinctus.html.
[2] FISHBASE. Johnius belangerii [DB/OL]. [2021-01-03]. www.fishbase.se/summary/Johnius-belangerii.html.
[3] 台湾鱼类资料库. Johnius distinctus [DB/OL]. [2021-01-03]. http://fishdb.sinica.edu.tw/mobi/species.php?id=382470.html. [4] 台湾鱼类资料库. Johnius belangerii [DB/OL]. [2021-01-03]. http://fishdb.sinica.edu.tw/mobi/species.php?id=382469.html. [5] 赵盛龙, 徐汉祥, 钟俊生, 等. 浙江沿海鱼类志[M]. 杭州: 浙江科学技术出版社, 2017: 693-696. [6] 刘鸿, 叶振江, 李增光, 等. 黄海中部近岸春夏季鱼卵、仔稚鱼群落结构特征[J]. 生态学报, 2016, 36(12): 3775-3784. [7] 李凡, 徐炳庆, 吕振波, 等. 莱州湾鱼类群落优势种生态位[J]. 生态学报, 2018, 38(14): 5195-5205. [8] 丁丽丽, 杨红, 戴小杰, 等. 浙江瓯飞围填海区渔业资源优势种生态位分析[J]. 上海海洋大学学报, 2016, 25(1): 116-123. doi: 10.12024/jsou.20150401399 [9] 袁华荣, 陈丕茂, 秦传新, 等. 南海柘林湾鱼类群落结构季节变动的研究[J]. 南方水产科学, 2017, 13(2): 26-35. doi: 10.3969/j.issn.2095-0780.2017.02.004 [10] 谢仰杰, 李军, 黄良敏, 等. 2006~2007年福建沿岸海域石首鱼类资源量时空变化[J]. 台湾海峡, 2012, 31(3): 403-411. [11] 王东旭, 陈国宝, 汤勇, 等. 大亚湾南部海域渔业资源水声学评估[J]. 安徽农业科学, 2017, 45(6): 95-98, 159. doi: 10.3969/j.issn.0517-6611.2017.06.033 [12] 徐兆礼. 春夏季闽江口和兴化湾鱼类数量特征的研究[J]. 水产学报, 2010, 34(9): 1395-1403. [13] 张云雷, 徐宾铎, 张崇良, 等. 海州湾春季皮氏叫姑鱼摄食习性的空间异质性及其索饵场分布[J]. 中国水产科学, 2020, 27(3): 315-326. [14] YE Z J, ZHANG C, PANHWAR S K, et al. Ageing Belanger's croaker, Johnius belangerii (Cuvier, 1830), based on otolith shape analysis[J]. J Appl Ichthyol, 2015, 31(1): 27-31. doi: 10.1111/jai.12451
[15] ALIREZA S, MOHAMMAD T R, FAZEL A M, et al. Heavy metal concentration in Belanger's croaker fish, Johnius belangerii from petrochemical waste receiving estuary in the Persian Gulf. [C]//Proceedings of International Conference on Environmental Engineering and Applications. Asia-Pacific Chemical, Biological & Environmental Engineering Society, 2011: 5.
[16] AZMAT R, TALAT R, AHMED K. The length-weight relationship, condition factor and impact of fluorideconcentration in Johnius belangerii of Arabian Sea[J]. J Environ Toxicol, 2007, 1(3): 138-143. doi: 10.3923/rjet.2007.138.143
[17] 王凯, 章守宇, 汪振华, 等. 马鞍列岛海域皮氏叫姑鱼渔业生物学初步研究[J]. 水产学报, 2012, 36(2): 228-237. [18] 郭慧, 徐明, 申玉春, 等. 流沙湾海域皮氏叫姑鱼生物学特性的初步研究[J]. 海南热带海洋学院学报, 2017, 24(2): 7-12. [19] 申世常, 黄良敏, 王家樵, 等. 厦门海域皮氏叫姑鱼渔业生物学特性的初步研究[J]. 海洋湖沼通报, 2020(1): 129-135. [20] CHU W S, WANG J P, HOU Y Y, et al. Length-weight relationships for fishes off the southwestern coast of Taiwan[J]. Afr J Biotechnol, 2011, 10(19): 3945-3950.
[21] WANG J Q, HUANG L M, LI J, et al. Length-weight relationships of 45 fish species in the Min River Estuary, East China Sea[J]. J Appl Ichthyol, 2016, 32(1): 131-133. doi: 10.1111/jai.12910
[22] WANG X H, QIU Y S, ZHU G P, et al. Length-weight relationships of 69 fish species in the Beibu Gulf, northern South China Sea[J]. J Appl Ichthyol, 2011, 27(3): 959-961. doi: 10.1111/j.1439-0426.2010.01624.x
[23] 李永振, 陈国宝, 孙典荣. 珠江口鱼类组成分析[J]. 水产学报, 2000, 24(4): 312-317. [24] 孙典荣, 李渊, 林昭进, 等. 海南岛近岸海域鱼类群落结构研究[J]. 中国海洋大学学报(自然科学版), 2011, 41(4): 33-38. [25] 乔延龙, 林昭进, 邱永松. 北部湾秋、冬季渔业生物群落结构特征的变化[J]. 广西师范大学学报(自然科学版), 2008, 26 (1): 100-104. [26] 徐兆礼. 兴化湾和闽江口海域春夏季鱼类区系特征[J]. 生物多样性, 2011, 19(1): 79-84. [27] 于海成. 长江口及邻近海域鱼类群落结构分析[D]. 北京: 中国科学院研究生院, 2008: 26-45. [28] 黄良敏. 厦门湾渔业资源与生态环境[M]. 北京: 中国农业出版社, 2018: 67. [29] 黄良敏, 谢仰杰, 张雅芝, 等. 厦门海域渔业资源现存量评析[J]. 集美大学学报(自然科学版), 2010, 15(2): 81-87. [30] 戴小杰, 杨志金, 田思泉, 等. 浙江南部近海鱼类分类多样性研究[J]. 海洋学报, 2019, 41(8): 43-51. [31] 陈新军. 渔业资源生物学[M]. 北京: 科学出版社, 2017: 95-96. [32] 詹秉义. 渔业资源评估[M]. 北京: 中国农业出版社, 1995: 27. [33] PAULY D. ELEFAN: users instruction and program listings[M]. Manila: ICLARM, 1980: 12-71.
[34] PAULY D. On the interrelationships between natural mortality, growth parameters and mean environmental temperature in 175 fish stocks[J]. ICES J Mar Sci, 1980, 39(2): 175-192. doi: 10.1093/icesjms/39.2.175
[35] VELDE T D, GRIFFITHS S P, FRY G C. Reproductive biology of the commercially and recreationally important cobia Rachycentron canadum in northeastern Australia[J]. Fish Sci, 2010, 76: 33-43.
[36] 张云雷, 薛莹, 于华明, 等. 海州湾春季皮氏叫姑鱼栖息地适宜性研究[J]. 海洋学报, 2018, 40(6): 83-91. [37] RICKER W E. Computation and interpretation of biological statistics of fish population[J]. Bull Fish Res Bd Can, 1975, 191: 1-382.
[38] 王家樵, 张雅芝, 黄良敏, 等. 福建沿岸海域主要经济鱼类生物学研究[J]. 集美大学学报(自然科学版), 2011, 16(3): 161-166. [39] GULLAND J A. Fish stock assessment: a manual of basic methods[M]. New York: Wiley, 1983: 87-125.
[40] 任一平. 胶州湾渔业资源与栖息环境[M]. 北京: 中国农业出版社, 2018: 68. [41] 严利平, 刘尊雷, 张辉, 等. 小黄鱼生物学特征与资源数量的演变[J]. 海洋渔业, 2014, 36(6): 481-488. doi: 10.3969/j.issn.1004-2490.2014.06.001 [42] 张学健, 程家骅, 沈伟, 等. 黄鮟繁殖生物学研究[J]. 中国水产科学, 2011, 18(2): 290-298. [43] 余文娟, 沈建忠, 龚江, 等. 长江中游贝氏䱗繁殖生物学研究[J]. 淡水渔业, 2018, 48(3): 53-60. doi: 10.3969/j.issn.1000-6907.2018.03.009 -
期刊类型引用(3)
1. 廖东荣,黎傲雪,陈宁,吕少梁,王学锋. 湛江港皮氏叫姑鱼种群生物学特征和资源开发状态. 南方水产科学. 2024(03): 27-35 . 本站查看
2. Wen Yang,Wenjia Hu,Bin Chen,Hongjian Tan,Shangke Su,Like Ding,Peng Dong,Weiwei Yu,Jianguo Du. Impact of climate change on potential habitat distribution of Sciaenidae in the coastal waters of China. Acta Oceanologica Sinica. 2023(04): 59-71 . 必应学术
3. 杨雯,胡文佳,陈彬,谭红建,董鹏,俞炜炜,杜建国. 基于MaxEnt模型的中国近海主要石首鱼科鱼类潜在适生区. 生态学杂志. 2022(09): 1825-1834 . 百度学术
其他类型引用(5)