Analysis of morphological differences and discrimination between female and male Cololabis saira based on geometric morphometrics
-
摘要:
秋刀鱼 (Cololabis saira) 广泛分布于西北太平洋海域,是中国远洋渔业的主要捕捞对象之一。为探究其雌、雄个体间的形态差异并有效区分其性别,根据2022年5—11月西北太平洋公海采集的150尾秋刀鱼样本的鱼体图片和生物学数据,首次采用几何形态测量学的地标点法对雌、雄秋刀鱼进行形态学研究,并通过判别分析构建了性别鉴定模型。其中样本分为2批,第1批107尾用于形态分析和鉴定模型的构建,第2批43尾用于检验模型的实际应用情况。结果表明,在相对扭曲主成分分析中,第1和第2主成分分别解释了总变异的63.43%和11.79%,降维效果和散点图区分效果均较好。其中I型和II型地标点的累积贡献率分别为57.94%和41.83%,说明在秋刀鱼雌、雄个体的形态区分上,I型和II型地标点的作用较大,III型地标点的作用较小。薄板样条分析结果显示,秋刀鱼雌、雄个体的形态差异主要表现在眼部、躯干前部和尾部。第1批107尾秋刀鱼的判别分析和交叉验证结果显示,判别正确率分别为91.6%和88.8%。将该模型用于第2批43尾秋刀鱼的性别鉴定时,判别正确率为88.45%。综上所述,几何形态测量学分析可以作为秋刀鱼形态学研究和性别鉴定的有效方法。
Abstract:Pacific saury (Cololabis saira) is widely distributed in the Northwest Pacific Ocean and is one of the main targets of Chinese pelagic operations. To investigate the morphological differences between female and male individuals and differentiate the two sexes effectively, we applied geometric land mark method morphometrics for the first time on 150 Pacific saury samples which were collected from the Northwest Pacific high seas during May–November, 2022, so as to analyze their morphology and establish a sex differentiation model via discriminant analysis with fish images and biological data. The samples were divided into two batches: Batch 1 (107 individuals) was used for morphological analysis and model construction, while Batch 2 (43 individuals) was used to test the model's practical application. The results show that in the relative warp analysis, the first and second principal components explained 63.43% and 11.79% of the total variation, respectively, with good dimensionality reduction and scatter plot separation effects. Landmark I and II contributed cumulatively to 57.94% and 41.83% of the variance, respectively, indicating that these landmark types play a significant role in the discrimination between two sexes, whereas Landmark III had lesser influence. The results of the thin-plate spline analysis show that the morphological differences between them were mainly in the eyes, the front part of the trunk and the tail. The results of discriminant analysis and cross-validation of Batch 1 showed correct discrimination rates of 91.6% and 88.8%, respectively. When applying the model to Batch 2 for sex discrimination, the correct discrimination rate was 88.45%. Thus, geometric morphometrics analysis is an effective method for morphological studies and sex identification of C. saira.
-
带鱼 (Trichiurus japonicus) 隶属鲈形目、带鱼科、带鱼属,为广泛分布于温带、亚热带和热带水域的暖水性底层鱼类,是我国沿海的重要经济鱼种,也是南海北部底拖网和刺网的主要渔获对象[1]。朱江峰和邱永松[2]根据20世纪90年代调查数据对南海北部带鱼的生长、死亡等参数进行了估算,并利用动态综合模型模拟其资源开发状况,结果表明南海北部带鱼能承受较大的捕捞压力,但幼鱼遭到过度捕捞,资源利用不合理;王跃中等[3]对南海北部带鱼的渔获量、捕捞努力量及环境数据进行了综合分析,探讨了其渔获量对捕捞压力及气候变化的响应;颜云榕[4]通过对北部湾带鱼食物组成的定量分析,研究了带鱼的摄食习性和摄食生态。
鱼类生长、死亡和性成熟等生物学参数的评估是探索渔业资源变动规律和开展渔业资源评估的基础[5],尤其在数据缺乏条件下的渔业资源评估中具有重要意义[6]。随着捕捞压力增大和环境恶化加剧,鱼类的生长、死亡等生物学特征必将产生相应的变化,如渤、黄海的小黄鱼 (Larimichthys polyactis)[7]、东海鲐鱼 (Scomber japonius)[8]和北部湾的短尾大眼鲷 (Priacanthus macracanthus)[9]等均出现了个体规格小型化、性成熟提前和生长速度加快等现象。然而,近20年来南海北部带鱼的生物学特征是否发生了转变尚缺乏相应的研究。因此,本文根据2014—2015年南海北部底拖网调查的渔业生物学资料,对该时期带鱼的渔获率、群体结构进行分析,并估算其生长及死亡等参数,以期为南海北部带鱼资源的评估及合理利用提供基础数据。
1. 材料与方法
1.1 数据来源
研究样品数据来自2014—2015年南海北部底拖网调查资料,调查站点根据其资源按水深分布[10]的特点设计,站位图 (图1) 与蔡研聪等[11]的研究一致。调查按季节进行,分别是2014年夏季 (7—8月)、秋季 (10—11月) 与2015年冬季 (1—2月)、春季 (4—5月)。采集的带鱼样本冷冻保存后带回实验室测定肛长、体质量、性成熟度和胃饱满度等生物学数据。每个站位带鱼的尾数大于50尾时测定50尾,少于50尾时全部测定,4个航次调查共测定带鱼1 128尾 (表1)。
表 1 带鱼的群体结构Table 1. Population structure of T. japonicus采样时间
Sampling time尾数
Number性比 (♀/♂)
Sex ratio肛长 Anal length/mm 体质量 Body mass/g 范围
Range均值
Mean优势组
Dominant范围
Range均值
Mean优势组
Dominant春 Spring 243 2.96 70~451 185.0 130~170 11~1200 128.5 20~60; 80~90 夏 Summer 428 1.66 61~505 184.1 160~210 10~1462 101.7 30~50; 60~80 秋 Autumn 178 0.67 87~500 246.8 230~250 9~1400 232.4 190~200 冬 Winter 279 1.11 92~528 199.9 160~200 13~1650 136.8 50~90 1.2 群体结构分析
带鱼样本按肛长10 mm为组距分组,分别进行肛长频次分析,其中占样品总数8%及以上的肛长组定义为优势肛长组。采用非参数Kol-mogorov-Smirnov Z检验分析雌、雄个体肛长组成的显著性差异;采用单因素方差分析研究4个季度带鱼肛长均值的显著性差异。
1.3 肛长体质量关系
带鱼肛长、体质量的测定分别精确到1 mm和1 g。肛长与体质量关系由函数W=aLb进行拟合,式中W表示体质量,L为肛长,a为生长的条件因子,b为异速生长因子。运用t检验分析雌雄个体的肛长体质量关系是否存在显著性差异。将相对体质量rw (%) 的平均值表示为肥满度[12],计算公式为rw=(W/aLb)×100%。
1.4 性比与性腺发育
统计每个季度以及肛长组中带鱼的雌、雄个体数目及性腺发育状况,性腺按Ⅰ—Ⅵ期标准进行划分,规定Ⅲ期及以上的为性成熟个体,通过适合性卡方检验对不同季度和肛长组的雌雄比例是否为1∶1进行分析。根据带鱼的生殖习性,采用春、夏季雌性样品计算50%性成熟肛长 (L50),由最小二乘法 (Least square method, LSM) 拟合以下逻辑斯蒂模型得到[13]:
$${P_i} = \frac{1}{{1 + {{\rm e}^{[ - r({L_i} - {L_{50}})]}}}}$$ (1) 其中Pi为肛长组Li对应的成熟个体比例,r为模型参数。
1.5 生长、死亡参数估计
带鱼的生长关系可由von Bertalanffy生长方程[5]进行拟合:
$$\begin{split} &\;\qquad \qquad \qquad \\ &{L_t} = {L_\infty }\{ 1 - \exp [ - K(t - {t_0})]\} \end{split}$$ (2) 式中Lt为t龄的肛长,L∞表示渐近肛长,K为生长系数,t0表示肛长为0时的理论年龄。参数L∞和K通过FISAT II (FAO-ICLARM Stock Assessment Tool II) 软件中的ELEFAN I (Electronic Length Frequency Analysis I) 估算。t0及体质量生长拐点tp可用以下经验公式求得[5]:
$$ {\log _{10}}( - {t_0}) = - 0.392\;2 - 0.275{\log _{10}}{L_\infty } - 1.038{\log _{10}}K $$ (3) $${t_p} = \frac{{{\rm{ln}}b}}{K} + {t_0}$$ (4) 根据生长参数L∞和K的评估结果,总死亡系数 (Z) 由线性体长变换渔获曲线法求得[5],Z的95%置信区间利用线性回归计算。自然死亡系数(M)通过Pauly经验公式估算[14]:
$${\rm{ln}}(M)= - 0.152-0.279{\rm{ln}}({L_\infty })+0.654{\rm{ln}}K + 0.463{\rm{ln}}T$$ (5) 其中T为栖息海域的年平均水温 (℃)。本研究中平均水温为调查站点的实际测量数据,为23.4 ℃。
捕捞死亡系数(F)和开发率(E)分别通过以下公式计算[5]:
$$ F = Z - M $$ (6) $$E = F/Z$$ (7) 2. 结果
2.1 渔获率的季节变化
本研究中,南海北部带鱼主要分布在水深小于100 m的海域,且渔获率随季节变化明显 (图2)。夏季渔获率较高的区域集中在北部湾和海南岛东北部海域;春、秋季渔获率较高的区域分布广泛,从北部湾至广东省东南部沿海海域皆有分布;冬季带鱼群体渔获率较低。
2.2 群体结构特征
从2014—2015年南海北部带鱼肛长频率分布 (图3-a,表1) 分析可知,带鱼肛长频率分布为明显的单峰型,肛长介于61~528 mm,优势肛长组为160~190 mm,占29.2%,雌、雄个体肛长组成无显著差异 (P>0.05)。秋季的带鱼个体最大,肛长均值为246.8 mm,体质量均值为232.4 g,远超其他三季,其次是春季,夏季个体最小,单因素方差分析表明4个季节的肛长均值差异显著 (P<0.05)。
2.3 肛长体质量关系
用LSM方法评估2014—2015年带鱼肛长体质量关系为W=4.09×10−5L2.79,a的95%置信区间为 (3.30×10−5,4.88×10−5),b的95%置信区间为 (2.76, 2.82),经t检验,雌雄个体的肛长体质量关系无显著差异 (P>0.05)。由肥满度随肛长变化 (图3-c) 得出,肛长介于100~150 mm,肥满度波动较大,肛长大于150 mm时肥满度保持约100%,无较明显的变化;当肛长为380 mm时肥满度最低 (88%)。此外,由于未测定纯质量数据,仅用体质量进行分析的情况下,存在一定的偏差,肛长小于100 mm时肥满度较高 (最高超过200%)。
2.4 性比与性腺发育
除去151尾雌雄不分个体,2014—2015年南海北部带鱼雌、雄个体性比为1.54。从4个季节性比可看出,春季雌性个体占比较大,夏季雌性个体占比有所下降,秋季雌性个体数远低于雄性,冬季雌性个体占比回升 (表1)。适合性卡方检验分析表明,4个季节雌雄比与1∶1均差异显著 (P<0.05)。性比随肛长变化 (图4-a),肛长组60~70 mm、80~90 mm、380~390 mm、400~410 mm、430~440 mm、450~470 mm、490~500 mm、520~530 mm均为雄性,肛长组360~370 mm均为雌性。肛长组130~140 mm、260~270 mm、300~310 mm、330~340 mm、350~360 mm、420~430 mm、500~510 mm性比与1∶1无显著差异 (P>0.05) 外,其他肛长组性比与1∶1均有显著差异。
4个季节中秋季雄性的性成熟比例高于雌性,而春、夏、冬季雌性性成熟比例均高于雄性(图5)。雌性中夏季的性成熟比例高于其他三季,秋季性成熟比例最低;而雄性中秋季性成熟比例最高,春季性成熟比例最低。各季度中未达到性成熟的个体占比均超过50%,且雄性中无性腺成熟度达到Ⅵ期的个体。南海北部带鱼样品分析表明,2014—2015年南海北部带鱼的最小性成熟肛长为113 mm。逻辑斯蒂方程拟合结果 (图4-b) 表明,50%性成熟肛长为241.5 mm。
2.5 生长、死亡参数估计
由ELEFAN I技术估算得到带鱼von Bertalanffy生长参数L∞、K分别为585 mm和0.2 a−1,由经验公式 (3) 求得t0为0.37年。因此带鱼的生长方程为:
$${L_t} = 585 \times (1 - {{\rm{e}}^{ - 0.20 \times (t + 0.37)}})$$ (8) $${W_t} = 1\;763 \times {(1 - {{\rm{e}}^{ - 0.20 \times (t + 0.37)}})^{2.79}}$$ (9) Wt为t龄的体质量。根据体质量生长方程求得的体质量拐点年龄为4.76龄,对应拐点的体质量为722 g。线性体长变换渔获量曲线法估算的总死亡系数为1.172 (图6),95%置信区间为 (0.582, 1.762)。由经验公式 (5) 算得M为0.475,由公式 (6) 和 (7) 可算得F为0.697,E为0.59。
3. 讨论
3.1 与其他海域带鱼种群参数的比较
不同海区带鱼的生长、死亡和性成熟等生物学参数均存在较大差异 (表2)。由于不同海域的捕捞压力和生活环境相差较大,带鱼个体规格存在显著差异。福建海区带鱼个体的平均规格最小,印尼阿拉弗拉海的平均规格最大,而北部湾带鱼肛长的分布范围最广,最大肛长可达896 mm,最小仅20 mm。肛长体质量关系式中b与鱼类生活环境相关,本研究与其他海域相比,b差异较大,东海区带鱼b最大,印尼阿拉弗拉海带鱼的最小。这种不同海域的差异可能主要是由其营养条件所致[5]。L∞与K是描述鱼类个体生长规律的重要参数,分别代表了鱼类极限体长和生长曲线的平均曲率[5]。与本研究相比,伊朗阿曼湾海域带鱼极限肛长和生长参数均偏大,而东海海域带鱼极限肛长最小,生长参数最大。这种差异可能是不同海域带鱼的种群密度、栖息水温和饵料生物等不同所引起的。L50是带鱼渔业管理中重要的生物学参考点,由表2可知,各海域带鱼的L50差异显著,当前南海北部带鱼的最小,印尼阿拉弗拉海带鱼的最大。M与E是反映渔业资源开发现状的重要参数,以Gulland[23]提出的鱼类最适E为0.5作为判断依据,南海北部带鱼处于过度开发状态。
表 2 不同海区带鱼的种群参数Table 2. Population parameters of T. japonicus in different sea areas海区及采样时间
Sea area and sampling time肛长范围
Anal length range/mm平均
肛长
Mean anal
length/mm异速生
长因子
b渐近肛长
L∞/mm生长
系数
K/a50%性
成熟肛长
L50/mm自然死
亡系数
M总死亡
系数
E参数评估方法
Assessment method for parameters印尼阿拉弗拉海 (2008)[15]
Arafura Sea of Indonesia123~650 343.00 2.420 − − 368.0 − − 体长频率法
Length-frequency method缅甸外海 (2012—2013)[16]
Andaman Sea of Myanmar101~422 212.26 2.752 − − 242.7 − − 体长频率法
Length-frequency method伊朗阿曼湾 (2017—2018)[17]
Iran's Gulf of Oman119~580 − 2.914 610.4 0.220 242.3 − − 体长频率法
Length-frequency method福建闽南渔场 (2000—2002)[18]
Minnan Sea of Fujian52~332 160.00 2.529 − − − − − 体长频率法
Length-frequency method东海 (2002—2004)[19]
East China Sea− − 3.006 493.0 0.346 − − − 年龄鉴定法
Age-determination method东海近海 (2007—2008)[20]
Costal waters of East China Sea126~332 182.62 2.795 − − 164.7 − − 体长频率法
Length-frequency method北部湾 (2006—2007)[21]
Beibu Gulf20~896 − 2.970 − − − − − − 南海北部 (1981—1982)[22]
Northern South China Sea132~605 342.84 2.400 622.0 0.160 322.0 0.212 0.71 体长频率法
Length-frequency method南海北部 (1997—1999)[2]
Northern South China Sea20~670 187.60 3.028 700.0 0.270 276.0 0.390 0.87 体长频率法
Length-frequency method南海北部 (2014—2015) (本研究)
Northern South China Sea (This study)61~528 198.10 2.790 585.0 0.200 241.5 0.475 0.59 体长频率法
Length-frequency method注:−. 文献中未查找到相应的参数值或评估方法 Note: −. Values of parameter or assessment method are not found in the References. 3.2 南海北部带鱼生物学特征的年际变化
大多数鱼类的表型特征在经过较强的捕捞压力胁迫后会产生适应性响应,其中个体规格的变化较为明显[24]。3个时期南海北部带鱼肛长频率分布对比表明,在持续多年的捕捞胁迫下,南海北部带鱼的群体结构发生了较明显的转变,肛长频率分布从20世纪80年代初的多峰型变为90年代末的双峰型,至今已变为较明显的单峰型 (图7),渔获物肛长从342.84 mm下降至187.60 mm,虽然当前其平均肛长小幅度上升至198.1 mm (表2),但群体结构依然呈现简单化的趋势。需要注意的是,不同时期调查所使用的网具和调查船存在一定的差异,可能会对结果造成一定影响。相关渔业研究表明,当前南海带鱼群体组成以1龄鱼和当年鱼为主,渔获物趋于小型化,渔获质量显著下降[25-26]。从渐近肛长和50%性成熟肛长来看,L∞、L50分别从80年代的622和322.0 mm减小到现在的585和241.5 mm (表2),表明南海北部带鱼存在明显的小型化和性成熟提前现象,这可能是对捕捞压力的响应,并且与北部湾深水金线鱼 (Nemipterus bathybius)[27]、东海小黄鱼[28]等一致。K从20世纪80年代至今增加了0.04,表明南海北部带鱼生长加快,且b呈增大趋势,这可能是带鱼对资源结构性衰退的生物学适应性响应所致[29]。
在较高强度捕捞压力下,南海北部带鱼的生物学特征发生了明显变化,结构趋于简单,个体小型化,造成种群参数变化的原因可归结于2个:1) 高强度的捕捞压力使得种群密度降低,来自捕食及种间竞争的压力减小,剩余群体能够获得相对充足的饵料供给,从而使得小型个体迅速进入渔业[30],这种变化在降低捕捞压力的情况下是可逆的,如日本北海道的大麻哈鱼 (Oncorhynchus keta)[31];2) 捕捞驱动进化 (Fisheries-induced evolution, FIE),即长期的过度捕捞迫使种群的某种基因发生变化,该变化为不可逆。笔者分析认为南海北部带鱼表型特征变化可能是由第一种原因所引起。自1999年起,我国政府开始对南海区域的渔船数和渔船功率进行管控,调整了渔业生产结构,同时大力建设人工鱼礁并开展增殖放流,加上休渔期的实施和捕捞成本的不断攀升,使得南海渔业资源的开发强度有所降低[32]。有研究表明,南海近年来的渔业政策对鱼类生物学特征的变化存在积极影响,如北部湾二长棘犁齿鲷 (Evynnis cardinalis) 的平均体长、L50等有所恢复[33]。本研究中数据显示南海北部带鱼资源的E从20世纪80年代初的0.71上升至90年代末的0.87,至今降至0.59,捕捞强度明显减小,但仍处于过度捕捞状态。在捕捞强度下降期间,带鱼的平均肛长由187.6 mm恢复至198.1 mm,K则从0.27减小至0.20,其生长加快的现象有所缓解,这一变化趋势与北部湾蓝圆鲹 (Decapterus maruadsi) 类似[34]。因此,亟需控制带鱼的捕捞压力,探索开展总量管理和限额捕捞[35],并加大对网具及开捕规格的管理,帮助其资源的恢复。
由于本研究按季节采样,样本量有限,而基于肛长频率 (ELEFAN) 对生物学样本量和频率分布要求较高,因此在估算L∞和K时会产生一定误差;此外,解释捕捞压力对经济鱼类的影响及其种群产生的适应性变化还需要来自种群遗传学和室内生态受控实验的进一步证据。
-
表 1 秋刀鱼的采样信息
Table 1 Sampling information of C. saira
样本批次
Batch of samples性别
Sex样本数
Number of samples/尾体长 Body length/mm 体质量Body mass/g 范围
Range均值±标准差
Mean±SD范围
Range均值±标准差
Mean±SD第1批
Batch 1雌♀ 52 219~281 260.65±13.36 65.3~109.2 75.90±12.92 雄♂ 55 237~301 272.75±10.53 72.9~115.3 97.71±12.73 第2批
Batch 2雌♀ 25 235~297 258.40±10.61 79.2~111.6 92.56±10.14 雄♂ 18 244~286 266.86±11.37 76.3~117.9 94.12±13.62 表 2 地标点类型与定义
Table 2 Landmark types and definition
地标点类型
Landmark type定义
DefinitionI型地标点 Landmark I 4 眼中心 7 胸鳍起点 11 背鳍起点 12 尾鳍基部上端 14 尾鳍基部下端 15 臀鳍起点 16 消化道末端开口 17 腹鳍起点 II型地标点 Landmark II 2 口裂末端 3 眼前缘 5 眼后缘 6 鳃盖上缘 8 鳃盖后缘 9 以7点向上作垂线与鱼体的交点 13 尾柄末端 III型地标点 Landmark III 1 吻端 10 躯干前部最大弯曲处 表 3 每个地标点在相对扭曲主成分分析中的贡献率
Table 3 Contribution rate of each landmark on relative warps
地标点
Landmark贡献率
Contribution rate/%1 0.15 2 1.53 3 17.21 4 46.47 5 15.14 6 4.62 7 5.11 8 0.19 9 0.14 10 0.05 11 0.01 12 1.03 13 3.01 14 0.88 15 2.15 16 2.31 17 0.02 表 4 秋刀鱼鱼体的雌雄判别方程系数
Table 4 Coefficients of female and male discrimination functions for body of C. saira
性状
Characteristics方程1雌性
Function 1 for female方程2雄性
Function 2 for maleRW1 −30.086 26.344 RW4 95.254 −91.117 RW5 86.747 −54.486 RW6 122.796 −99.142 RW7 −59.675 54.684 RW8 −42.726 132.809 RW18 529.118 −580.991 RW25 385.023 −363.849 常数Constant −1.532 −1.466 注:RW表示相对扭曲得分。 Note: RW represents relative warps scores. 表 5 第1批秋刀鱼的雌雄判别结果
Table 5 Results of female and male discrimination in Batch 1 of C. saira
项目
Item实际
性别
Actual
sex预测性别
Predicted sex判别正确率
Identification
accuracy/%综合
正确率/%
Total
accuracy/%雌性
Female雄性
Male逐步判别分析
SDA雌♀ 46 6 88.5 91.6 雄♂ 3 52 94.5 交叉验证
CV雌♀ 44 8 86.4 88.8 雄♂ 4 51 92.7 表 6 判别方程对第2批秋刀鱼的雌雄判别结果
Table 6 Results of function for female and male discrimination in Batch 2 of C. saira
实际性别
Actual sex预测性别 Predicted sex 判别正确率
Identification
accuracy/%综合
正确率
Total
accuracy/%雌性
Female雄性
Male雌性 Female 22 3 88.0 88.45 雄性 Male 2 16 88.9 -
[1] 赵峰, 章龙珍, 庄平, 等. 鲟科鱼类性别鉴别技术的研究进展及其应用[J]. 海洋渔业, 2009, 31(2): 215-220. [2] 陈宇舒. 尖头塘鳢 (Eleotris oxycephala) 雌雄个体形态差异分析[J]. 渔业研究, 2023, 45(4): 346-355. [3] 韩伟国, 刘新富, 孟振, 等. 醋酸洋红染色技术在大菱鲆幼鱼生理性别鉴定中的应用[J]. 渔业科学进展, 2011, 32(6): 37-42. [4] BLYTHE B, HELFRICH L A, BEAL W E, et al. Determination of sex and maturational status of striped bass (Morone saxatilis) using ultrasonic imaging[J]. Aquaculture, 1994, 125(1/2): 175-184.
[5] FERGUSON A, TAGGART B J, PRODOHL A P, et al. The application of molecular markers to the study and conservation of fish populations, with special reference to Salmo[J]. J Fish Biol, 2006, 47(1): 103-126.
[6] 仇登高, 徐世宏, 刘鹰, 等. 超声成像技术在大西洋鲑早期性别及发育期鉴别的应用研究[J]. 海洋科学, 2016, 40(6): 23-29. [7] 毕宜慧, 张鑫宇, 李杰, 等. 暗纹东方鲀雌雄个体形态差异及其判别分析[J]. 海洋渔业, 2021, 43(5): 563-572. [8] 邵锋, 陈新军, 李纲, 等. 东黄海鲐鱼形态差异分析[J]. 上海水产大学学报, 2008, 17(2): 204-209. [9] 白明, 杨星科. 几何形态测量法在生物形态学研究中的应用[J]. 昆虫知识, 2007, 44(1): 143-147. [10] 韩霈武, 陈新军, 方舟, 等. 基于鱼体和耳石形态的东海两种鲐属鱼类判别分析[J]. 海洋渔业, 2020, 42(2): 161-169. [11] SLICE D E. Modern morphometrics in physical anthropology [M]. New York: Springer US, 2005: 325-330.
[12] BRENO M, LEIRS H, VAN DONGEN S. Traditional and geometric morphometrics for studying skull morphology during growth in Mastomys natalensis (Rodentia: Muridae)[J]. J Mammal, 2011, 92(6): 1395-1406. doi: 10.1644/10-MAMM-A-331.1
[13] FRED L, BOOKSETINFL. Size and shape spaces for landmark data in two dimensions[J]. Stat Sci, 1986, 1(2): 175-184.
[14] BOOKSTEINFL. Morphometric tools for landmark data: geometry and biology[M]. Cambridge: Cambridge University Press, 1991: 18, 435.
[15] ROHLF F J, SLICE D E. Extensions of the procrustes method for the optimal superimposition of landmarks[J]. Syst Zool, 1990, 39(1): 40-59. doi: 10.2307/2992207
[16] 朱清澄, 花传祥. 西北太平洋秋刀鱼渔业[M]. 北京: 海洋出版社, 2017: 1-2. [17] NAKAYA M, MORIOKA T, FUKUNAGA K, et al. Growth and maturation of Pacific saury Cololabis saira under laboratory conditions[J]. Fish Sci, 2010, 76(1): 45-53. doi: 10.1007/s12562-009-0179-9
[18] 张孝民, 石永闯, 李凡, 等. 基于MAXENT模型预测西北太平洋秋刀鱼潜在渔场[J]. 上海海洋大学学报, 2020, 29(2): 280-286. [19] KIMURA N, OKADA Y, MAHAPATRA K. Relationship between saury fishing ground and sea surface oceanographic features determined from satellite data along the northeastern coast of Japan[J]. J Mar Sci Technol, 2004, 2(2): 1-12.
[20] 朱文涛, 陈新军, 汪金涛, 等. 基于灰色系统的西北太平洋秋刀鱼资源丰度预测[J]. 广东海洋大学学报, 2018, 38(6): 13-17. [21] 石永闯, 朱清澄, 张衍栋, 等. 基于模型试验的秋刀鱼舷提网纲索张力性能研究[J]. 中国水产科学, 2016, 23(3): 704-712. [22] 花传祥, 高玉珍, 朱清澄, 等. 基于耳石微结构的西北太平洋秋刀鱼 (Cololabis saira) 年龄与生长研究[J]. 海洋学报, 2017, 39(10): 46-53. [23] 杨明树. 西北太平洋公海秋刀鱼繁殖生物学研究[D]. 上海: 上海海洋大学, 2017: 5-6. [24] 李非, 朱清澄, 花传祥, 等. 西北太平洋秋刀鱼主要体征及繁殖特性研究[J]. 水产科学, 2019, 38(2): 173-181. [25] 朱国平, 刘芳沁. 几何形态测量学及其在鱼类生态学研究中的应用进展[J]. 上海海洋大学学报, 2022, 31(5): 1180-1189. [26] ICHII T, NISHIKAWA H, MAHAPATRA K, et al. Oceanographic factors affecting interannual recruitment variability of Pacific saury (Cololabis saira) in the central and western North Pacific[J]. Fish Oceanogr, 2018, 27(5): 445-457. doi: 10.1111/fog.12265
[27] ROHLF F J, BOOKSTEIN F L. Proceedings of the Michigan morphometrics workshop[M]. New York: University of Michigan Museum of Zoology, 1990: 1-380.
[28] 张锋, 李萍. 几何形态测量学在古生物学研究中的应用概况[J]. 古生物学报, 2016, 55(4): 518-531. [29] VISCOSI V, CARDINI A. Leaf morphology, taxonomy and geometric morphometrics: a simplified protocol for beginners[J]. PLoS One, 2018, 6(10): e25630.
[30] 姜涛, 郑朝臣, 黄洪辉, 等. 基于地标点法的九龙江口和珠江口凤鲚和七丝鲚耳石形态学特征比较[J]. 南方水产科学, 2018, 14(6): 10-16. [31] SAUER F G, PFITZNER W P, JÖST H, et al. Using geometric wing morphometrics to distinguish Aedes japonicus japonicus and Aedes koreicus[J]. Parasit Vectors, 2023, 16(1): 417-418. doi: 10.1186/s13071-023-06032-4
[32] KEMP F. Modern applied statistics with S[J]. J R Stat Soc Ser D Statist, 2003, 52(4): 704-705.
[33] AKINOBU W. How many landmarks are enough to characterize shape and size variation?[J]. PlOS One, 2018, 13(6): e0198341. doi: 10.1371/journal.pone.0198341
[34] CARDINI A, SEETAH K, BARKER G. How many specimens do I need? Sampling error in geometric morphometrics: testing the sensitivity of means and variances in simple randomized selection experiments[J]. Zoomorphology, 2015, 134(2): 149-163. doi: 10.1007/s00435-015-0253-z
[35] 郁岳峰, 张勋, 黄洪亮, 等. 秋刀鱼舷提网集鱼方法的研究[J]. 浙江海洋学院学报(自然科学版), 2006(2): 154-156. [36] CRESPI-ABRIL A C, JOSEBARO P, MORSAN E M, et al. Analysis of the ontogenetic variation in body and beak shape of the Illex argentinus inner shelf spawning groups by geometric morphometrics[J]. J Mar Biol Assoc UK, 2010, 90(3): 547-553. doi: 10.1017/S0025315409990567
[37] ZISCHKE T M, GRIFFITHS P S, TIBBETTS R I, et al. Stock identification of wahoo (Acanthocybium solandri) in the Pacific and Indian Oceans using morphometrics and parasites[J]. ICES J Mar Sci, 2013, 70(1): 164-172. doi: 10.1093/icesjms/fss164
[38] 陈楠桦, 梁仁杰, 白义, 等. 基于几何形态测量学的四种滨螺形态差异与系统发生关系研究[J]. 海洋与湖沼, 2018, 49(6): 1365-1374. [39] 张秀霞, 朱巧莹, 赵俊. 利用几何形态测量学方法分析唐鱼群体的形态变异[J]. 水产学报, 2017, 41(9): 1365-1373. [40] 韩霈武, 王超, 方舟, 等. 基于几何形态测量的不同群体柔鱼个体形态生长变化[J]. 中国水产科学, 2022, 29(3): 483-493. [41] 侯刚, 王学锋, 朱立新, 等. 基于几何形态测量学的4种金线鱼矢耳石识别研究[J]. 海洋与湖沼, 2014, 45(3): 496-503. doi: 10.11693/hyhz20130300006 [42] 侯刚, 刘丹丹, 冯波, 等. 基于地标点几何形态测量法识别北部湾4种白姑鱼矢耳石形态[J]. 中国水产科学, 2013, 20(6): 1293-1302.