富山武装乌贼 (Enoploteuthis chunii) 最早发现于日本海富山湾，是一种小型头足类，其活动范围很广，昼深夜浅^[1]，在中国、韩国海域内也有分布^[2]。富山武装乌贼作为西太平洋海洋食物网的重要一环，是海洋哺乳动物和鱼类的重要食饵，常被大量发现于条纹原海豚 (Stenella caeru-leoalba) ^[3]、太平洋蓝鳍金枪鱼 (Thunnus orientalis)^[4]的胃中。目前国内外学者对其相关研究较少，Hong^[2]研究了富山武装乌贼的形态信息，但针对其生长规律、摄食习惯等其他信息的研究仍是空白。头足类生长发育过程中角质颚的形态大小、色素沉积等表型特征的变化与其摄食习性的转变息息相关^[5-7]。因此，了解和掌握富山武装乌贼角质颚的生长发育规律，不仅有利于掌握其个体发育期摄食生态的变化，也有助于进一步研究其在海洋生态系统中的地位和作用。

角质颚是研究头足类的重要硬组织，作为头足类的摄食器官，有着结构稳定、不易腐蚀等优点^[8-11]。除此之外，角质颚还储藏了头足类丰富的个体生长信息 ^[12-14]。刘必林等^[12]研究了智利外海茎柔鱼 (Dosidicus gigas) 角质颚微结构，探讨了其生长规律；林静远等^[13]分析了剑尖乌贼 (Uro-teuthis edulis) 的日龄及生长；方舟等^[14]研究了北太平洋柔鱼 (Ommastrephes bartramii) 的角质颚形态，认为环境因素也会对角质颚生长产生影响。角质颚的色素沉积情况也是研究生长的重要依据，陆化杰等^[15]分析了中国南海西沙群岛海域鸢乌贼 (Sthenoeuthis oualaniensis) 色素沉积等级与胴长、体质量、性腺成熟度和角质颚形态参数的生长关系，陈炫妤等^[16]也通过此方法对西北印度洋鸢乌贼的色素沉积进行了研究。且由于富山武装乌贼昼夜垂直活动的生活习性为西太平洋众多生物提供了饵料，因此不同水层生活的鱼类胃中均能发现其角质颚，研究其角质颚的表型特征与生长的关系有助于研究西北太平洋海洋生物间的食物网关系。本研究将富山武装乌贼角质颚形态参数及色素沉积面积变化与胴长结合，通过研究其表型变化分析生长发育规律。

本研究根据2019年上海海洋大学“淞航号”在西太平洋海域进行渔业资源调查时采集的样本，通过研究富山武装乌贼角质颚的各部分形态指标及色素沉积面积的变化特点，结合胴长分析其不同个体发育时期的生长特点及规律，可为后续研究其捕食者的摄食规律提供相关资料，也可根据其捕食者胃含物中富山武装乌贼的数量来估算其在西太平洋的资源量^[17-18]。另外，本研究还结合富山武装乌贼摄食习性的变化分析了不同胴长组角质颚表型变化差异的原因。

1.   材料与方法

1.1   材料

富山武装乌贼样本由上海海洋大学“淞航号”远洋渔业资源调查船于2019年3月和9月采自西北太平洋海域 (148°00'E—148°10'E, 31°99'N—38°76'N)，共91尾 (表1)。样品在实验室进行解冻及清洗处理后，使用游标卡尺测量其胴长 (Mantel length, ML)，精确到1 mm。从头部提取角质颚，清洗后存放于含75%乙醇的离心管中^[18]。

表  1  富山武装乌贼采样信息

Table  1.  Sampling information of E. chunii

采样海域 Sampling sea area	采样日期 Sampling date	样品数量 Sample size	胴长范围 Mantle length range/mm
148°00'E, 33°00'N	2019年3月21日	25	32~62
148°00'E, 31°99'N	2019年3月23日	4	43~49
148°10'E, 38°76'N	2019年9月5日	62	18~45

下载: 导出CSV 

| 显示表格

1.2   数据处理方法

1.2.1   拍照及图片处理

使用体视显微镜OLYMPUSSZ61对角质颚的侧视图进行拍照，拍取的照片用 Photoshop CS 6.0 (Adobe Systems Inc.，2019) 对图像进行亮度、对比度和去噪声等处理。

1.2.2   角质颚长度测定

使用生物图像处理软件Digimizer测量角质颚形态学指标，包括上头盖长 (Upper hood length, UHL)、上脊突长 (Upper crest length, UCL)、上喙长 (Upper rostrum length, URL)、上侧壁长 (Upper lateral wall length, ULWL)、上翼长 (Upper wing length, UWL)、下头盖长 (Lower hood length, LHL)、下脊突长 (Lower crest length, LCL)、下喙长 (Lower rostrum length, LRL)、下侧壁长 (Lower lateral wall length, LLWL)、下翼长 (Lower wing length, LWL) (图1)，测量结果精确到0.01 mm^[19]。

图  1  角质颚形态参数示意图

A. 头盖长 (HL)；B. 脊突长 (CL)；C. 喙长 (RL)；D. 侧壁长 (LWL)；E. 翼长 (WL)。

Figure  1.  Scheme of morphometric measurements for beak

A. Hood length (HL); B. Crest length (CL); C. Rostrum length (RL); D. Lateral wall length (LWL); E. Wing length (WL).

下载: 全尺寸图片 幻灯片

1.2.3   角质颚色素面积计算

根据角质颚色素沉积情况，将其分为2个区域，色素沉积较深的部分为Ⅰ区，色素沉积较浅的部分为Ⅱ区 (图2)。由于通过拍照将角质颚转换为图片形式进行测量，因此其色素沉积面积均为实际的投影面积，使用Photoshop测量角质颚色素沉积面积，计算其占角质颚总面积的比例。

图  2  角质颚色素沉积分区示意图

Figure  2.  Partition of pigmentation in beak

下载: 全尺寸图片 幻灯片

1.2.4   统计分析

分别采用线性、指数、幂函数、对数关系对富山武装乌贼角质颚长度与胴长进行拟合，根据最小AIC原则选择最适方程^[20]。采用单因素方差分析法 (One-way ANOVA) 对不同胴长组富山武装乌贼角质颚长度和色素沉积进行差异性分析，显著水平P=0.05 ^[21]。

2.   结果

2.1   角质颚长度与胴长总体关系

富山武装乌贼角质颚各项形态参数随着胴长增加逐渐增大 (表2)。根据AIC原则分析显示，角质颚长度与胴长呈显著的线性关系 (图3)，其方程式为：

表  2  富山武装乌贼角质颚形态参数值

Table  2.  Morphometric parameters of beak of E. chunii mm

形态指标 Morphologic index	最大值 Maximum	最小值 Minimum	均值 Mean	标准差 Standard deviation
上头盖长 UHL	6.16	3.34	4.48	0.86
上脊突长 UCL	8.68	4.57	6.14	1.18
上喙长 URL	2.18	0.87	1.45	0.31
上侧壁长 ULWL	7.09	3.56	4.95	0.98
上翼长 UWL	3.36	1.55	2.38	0.50
下头盖长 LHL	3.25	1.87	1.47	0.32
下脊突长 LCL	3.43	1.69	2.32	0.44
下喙长 LRL	2.60	1.01	1.60	0.38
下侧壁长 LLWL	6.45	3.33	4.56	0.92
下翼长 LWL	4.12	1.90	2.60	0.52

下载: 导出CSV 

| 显示表格

图  3  角质颚各项形态参数与胴长关系

Figure  3.  Relationship between each beak morphometric parameter and mantle length

下载: 全尺寸图片 幻灯片

Y_ML=0.072 3X_UHL−2.008 (R²=0.730 3, N=91)

Y_ML=0.099 8X_UCL−2.721 1 (R²=0.735 8, N=91)

Y_ML=0.023 4X_URL−0.649 (R²=0.602 5, N=91)

Y_ML=0.083 5X_ULWL−2.094 (R²=0.743 3, N=91)

Y_ML=0.042 4X_UWL−0.928 2 (R²=0.746, N=91)

Y_ML=0.026 1X_LHL−0.570 7 (R²=0.703 8, N=91)

Y_ML=0.038 2X_LCL−1.013 3 (R²=0.772 7, N=91)

Y_ML=0.031 2X_LRL−0.529 7 (R²=0.697 2, N=91)

Y_ML=0.079X_LLWL−1.855 2 (R²=0.752 8, N=91)

Y_ML=0.042 7X_LWL−1.142 1 (R²=0.692 2, N=91)

2.2   角质颚色素沉积面积与胴长的关系

根据最小AIC原则分析显示，角质颚色素沉积面积随胴长总体呈线性增长趋势 (图4)。

图  4  角质颚色素沉积面积占比与胴长的关系

Figure  4.  Relationship between proportion of beak pigmentation area and mantle length

下载: 全尺寸图片 幻灯片

上颚色素沉积面积与胴长关系式为：

Y_上1=0.002 4X_ML−0.112 2 (R²=0.395 4, N=91)

Y_上2=0.006 1X_ML+0.112 9 (R²=0.594 9, N=91)

下颚色素沉积面积与胴长关系式为：

Y_下1=0.005 2X_ML+0.121 1 (R²=0.422 4, N=91)

Y_下2=0.002 9X_ML+0.106 2 (R²=0.396 6, N=91)

2.3   角质颚表型变化分析

2.3.1   角质颚各项形态参数变化

富山武装乌贼角质颚各项形态参数随着角质颚生长总体呈上升趋势，各项形态指标在31、41、62 mm处发生明显的跃变 (图5)。ANOVA进一步分析发现，角质颚各项形态指标在18~31 mm、32~41 mm、>42 mm均存在显著的组间差异 (P<0.05)。

图  5  角质颚各项形态参数随胴长变化图

Figure  5.  Variation of beak morphometrics variation with mantle length

下载: 全尺寸图片 幻灯片

2.3.2   角质颚色素沉积面积占比变化

富山武装乌贼角质颚色素沉积由喙部开始向侧壁延伸，呈依次减少的趋势，由深变浅，且随着胴长增加其色素覆盖面积也逐步增加；上颚头盖、下颚头盖和下颚侧壁在胴长32 mm处色素沉积发生跃变，上颚侧壁在胴长43 mm处色素沉积突变明显 (图6)。ANOVA进一步分析发现，角质颚上颚头盖Ⅰ区、上颚Ⅰ、Ⅱ区色素沉积在18~31 mm、>42 mm两个胴长组间存在显著性差异 (P<0.05)，上颚侧壁Ⅰ、Ⅱ区色素沉积在18~41 mm、>42 mm两个胴长组间存在显著性差异 (P<0.05)。上颚头盖Ⅱ区色素沉积随胴长增加不断变化，胴长组间无明显差异，随着角质颚生长，在胴长介于31~40 mm，上、下颚色素沉积面积增加明显。

图  6  角质颚各部分色素面积占比随胴长变化图

Figure  6.  Variation of proportion of beak pigmentation area with mantle length

下载: 全尺寸图片 幻灯片

2.3.3   个体发育期角质颚表型变化

根据角质颚长度和色素沉积度可按胴长将富山武装乌贼分为A组(18~31 mm)、B组(32~41 mm)、C组(>42 mm)共3组 (表3)。胴长18~31 mm时，上颚喙部相对较短，上、下颚侧壁色素沉积面积较小，色素沉积程度较浅，下颚翼部无色素沉积，下颚侧壁脊不明显；胴长31~41 mm时，上颚喙部长度增加且变尖，上下颚侧壁及头盖部分色素沉积明显增加，下颚翼部有轻微色素沉积，下颚侧壁脊略明显；胴长>42 mm时，上颚喙部尖且长，上、下颚色素沉积面积变大，程度加深，下颚翼部有明显色素沉积，下颚侧壁脊明显 (图7)。

表  3  角质颚表型变化信息

Table  3.  Information of beak phenotype change

胴长组 Mantle length group/mm	上颚 Upper beak			下颚 Lower beak
	喙部 Beak	侧壁色素沉积面积 Pigmentation area of lateral wall		翼部色素沉积面积 Pigmentation area of wing	侧壁脊 Ridge of lateral wall	侧壁色素沉积面积 Pigmentation area of lateral wall
18~31	较长	24%		无沉积	不明显	42%
32~41	尖且长	36%		较小	略明显	64%
42~62	尖且长	51%		较大	明显	65%

下载: 导出CSV 

| 显示表格

图  7  角质颚生长变化图

Figure  7.  Growth and change of beak

下载: 全尺寸图片 幻灯片

3.   讨论

3.1   角质颚大小、色素沉积与胴长关系

最小AIC分析显示，角质颚各项形态参数和胴长呈显著的线性相关，这与南海北部秋季杜氏枪乌贼 (U. duvaucelii) ^[22]、印度洋西北海域鸢乌贼^[23]结果相似。分析关系式斜率发现，角质颚各部分生长也存在一定差异，斜率越小就代表其生长越快。富山武装乌贼UCL、ULWL、LLWL斜率分别为0.099 8、0.083 5、0.079，因此上颚脊突部、侧壁部及下颚侧壁部生长较快。角质颚不同部位的生长速度与其摄食习惯相关，上、下颚的脊突部与侧壁部在富山武装乌贼摄食活动中主要起支撑作用，快速生长可方便其咬合食物^[24]。

角质颚色素沉积面积占比随胴长总体也呈增长趋势，这与阿根廷滑柔鱼 (Illex argentinus) 和秘鲁外海茎柔鱼的研究结果相似。方舟等^[25]分析阿根廷滑柔鱼角质颚色素变化及其与个体生长等因素的关系，发现色素沉积的平均等级总体上随着其生长而增加；胡贯宇等^[26]发现上喙长和下翼长是影响角质颚色素沉积等级的主要因素，且角质颚各形态参数与角质颚色素沉积等级呈正相关关系。

3.2   个体发育期角质表型变化

在个体发育的不同阶段，角质颚的增长速度也有不同。研究结果显示胴长18~31 mm角质颚生长较缓慢，胴长32~41 mm角质颚生长较快，胴长>42 mm角质颚生长速度又逐渐减缓 (图5)，这与柔鱼从仔鱼到稚鱼再到成鱼的生长规律相符^[27]。根据已有的研究资料显示，19 mm富山武装乌贼为仔鱼幼体 ^[4]；胴长38 mm的希氏武装乌贼 (E. higginsi) 为雌性未成熟个体，胴长62 mm的琼氏武装乌贼 (E. jonesi) 为成熟个体^[28]。希氏武装乌贼、琼氏武装乌贼与富山武装乌贼同为武装乌贼属，因此本研究胴长18~31 mm富山武装乌贼处于仔稚鱼期，胴长32~41 mm处于未成熟的幼鱼期，胴长>42 mm慢慢处于成鱼期。仔鱼期时，头足类摄食小型浮游动物、营养需求有限，因而角质颚发育较慢；稚鱼期时，个体快速生长需提高摄食量，以获取足够丰富的营养；成鱼期时，个体发育逐渐变大，发育逐渐趋于成熟，生长速度也逐渐趋于稳定，摄食习惯相对固定，因此其角质颚变化速度逐渐下降^[29-31]。

在个体发育的不同阶段，角质颚的色素沉积也不同。研究结果显示，富山武装乌贼角质颚上颚头盖、下颚头盖和下颚侧壁在胴长32 mm处色素沉积发生跃变，上颚侧壁在胴长43 mm处色素沉积跃变明显 (图6)。结合个体发育期分析认为，第一次色素沉积的跃变发生在仔稚鱼向未成熟幼体过渡期，这个阶段往往食性会发生改变，食物由小型浮游动物转变成较大的甲壳类；第二次色素沉积的跃变发生在幼体向成体过渡期，此时进入性腺成熟阶段，食物转向具有较高能量的鱼类^[32-34]。在富山武装乌贼生长过程中，由于其摄食习惯的改变，导致其色素沉积逐渐加深，角质颚的硬度也逐渐增加 ^[34-37]，进而有助于其更好地撕碎食物^[32]，提高进食效率。

综合角质颚大小及色素沉积变化，本研究将富山武装乌贼角质颚表型变化分为仔稚鱼期 (胴长18~31 mm)、幼鱼期 (胴长32~41 mm) 和成鱼期 (胴长>42 mm) 3个阶段 (图7)。仔稚鱼期，较短的上颚喙部和较浅的色素沉积适合捕食游泳能力弱的浮游动物类；幼鱼期，尖长的上颚喙部和较深的色素沉积适合捕食有壳的甲壳类；成鱼期，更加尖长的上颚喙部和更深的色素沉积适合捕食游泳迅速的鱼类^[10]。

3.3   展望与不足

富山武装乌贼具有昼夜垂直移动的习性，一般是昼深夜浅，这为西太平洋众多生物提供了丰富的饵料。角质颚因具有耐腐蚀特性长存于捕食者的胃中^[3-4]。因此，在掌握个体发育角质颚的生长变化基础上，有助于通过分析捕食者胃含物中残留的角质颚来评估富山武装乌贼的资源量，揭示其在海洋生态系统中的作用^[18]。目前，有关富山武装乌贼的研究较少，仅针对特定海域内以及一些海洋生物胃含物中残留个体有零星报道^[38]。本研究分析了富山武装乌贼不同生长发育时期角质颚长度及色素沉积面积的变化，并结合其摄食、生活习性等分析其表型变化的原因和特点。本研究认为，富山武装乌贼角质颚表型的变化主要与其个体发育期的食性转变有关。然而，由于样本量偏少，并没有开展其食性方面的研究。因此，今后希望采用胃含物、脂肪酸及碳、氮稳定同位素等分析技术，来准确解析富山武装乌贼的摄食习性，以此来佐证本文的相关研究结论，同时为研究其他海洋生物的摄食习性提供参考。