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.
-
神经肽是从神经系统和内分泌系统产生和释放，参与调节生物个体的生理活动，维持生物体内环境稳态的一类重要信号因子，是NEI系统中最重要的调节因子之一，在脊椎和无脊椎动物中均有免疫调节功能。其中FMRFamide (Phe-Met-Arg-Pheamide) 是一种酰胺化四肽，是类 FMRF酰胺肽 (FMRFarmide-related peptides, Farps或FMRFamide-like peptide, FLPs) 家族的成员。其在神经系统中合成并释放，能通过神经元直接传递信号，并作为神经递质、神经调质和神经激素发挥作用[1]，是无脊椎动物中含量最丰富的一类神经肽[2]。此外，该神经肽也能通过内分泌系统调节摄食和运动行为[3]。近年来，FMRFamide的生理功能受到普遍关注，如Marciniak等[4]研究发现其具有特异性的心脏活性，能有效调节甲虫内脏肌肉的内源性收缩活性。Kim等[5]研究皱纹盘鲍 (Haliotis discus hannai) 发现，FMRFamide神经肽对海洋腹足类动物生殖调节有重要促进作用。FMRFamide也被发现能够通过抑制有丝分裂来保护细胞免受凋亡[6]。而作为无脊椎动物中特有的神经肽，一些研究也表明其可能参与了海洋无脊椎动物的免疫应答。如Li等[7]通过对太平洋牡蛎 (Crassostrea gigas) 炎症刺激后的FMRFanide多肽治疗，发现FMRFamide多肽不仅由神经内分泌系统合成，而且还在血细胞中合成加工，首次证实了该多肽存在于血细胞中。Guan等[8]对加利福尼亚海兔 (Aplysia california) 的研究发现FMRFamide可以通过激活P38丝裂原激活蛋白 (MAP) 激酶来调节突触的可塑性和调节免疫效应物。
瘤背石磺 (Onchidium reevesii) 隶属于软体动物门、石磺科，又名土鸡、海赖子、涂龟、土海参等[9-10]，常栖息在滩涂上的石缝中、植被的根部和阴暗潮湿的角落，主要分布于我国东部沿海地区[11-12]。瘤背石磺的中枢神经系统 (Central nervous system, CNS) 是一个简单的环状结构，体积大且易获得[13]。瘤背石磺不仅是研究无脊椎动物神经系统功能的模式种，而且由于其生活在潮间带，有良好的应激耐受机制，能适应波动的生存环境、各种重金属离子和病原体[14]，对于研究无脊椎动物应对外界环境变化和病原入侵的免疫反应机制也有显著优势。为此，本研究在实验室瘤背石磺高通量转录组测序的基础上克隆得到了OrFMRFamide基因的全长，分析该基因的分子特征及其mRNA和多肽在组织中的分布，研究其在炎症刺激后的反应，从而为探究FMRFamide基因在瘤背石磺免疫机制中的作用提供了理论依据，也为进一步理解生活在潮间带的两栖类动物对恶劣环境的适应机理奠定基础。
1. 材料与方法
1.1 实验材料
2020年6月初在江苏省盐城市沿海滩涂采集成体鲜活瘤背石磺，带回实验室后暂养于70 cm×120 cm×50 cm的塑料养殖箱内，养殖箱底部覆盖5~10 cm厚的海泥，种植滩涂采集的植被，放上瓦片，放养滩涂捕捉的螃蟹在泥中打洞，为石磺提供休息躲避的巢穴，尽可能地模拟自然生态环境，消除实验外干扰。每天早晚定时投喂玉米粉和喷洒海水保持土壤湿润，定期检查瘤背石磺生长状况，以确保其较高的成活率。暂养1周后，取5只体态均匀、健康有活力的瘤背石磺 [ 体长 (43.55±1.32) mm，体质量 (15.73±0.78) g]，用超纯水清洗表面的泥沙，取血细胞、神经节、肝胰腺、肌肉、性腺、皮肤和腹足等组织用液氮速冻后放入–80 ℃冰箱备用，用于组织cDNA模板的制备。
1.2 实验方法
1.2.1 总RNA的提取和cDNA的合成
实验所用器皿均经180 ℃高温杀菌处理3 h或使用DEPC处理过后进行高压灭菌。实验开始前将所有器材预冷备用，全程在冰上操作，取–80 ℃保存的瘤背石磺各组织，按照Trizol (Invitrogen, 美国) 说明书提取各组织的总RNA。配制质量浓度为10 mg·mL−1的琼脂糖凝胶，135 V、电泳15 min检测RNA提取质量，在核酸检测仪上标定提取RNA的浓度和纯度。将神经节RNA按照PrimeScript II 1st Strand cDNA Synthesis Kit (TaKaRa) 说明书合成cDNA第一链，作为后续实验的RACE (Rapid-amplification of cDNA ends) 模板。将提取所得的各组织RNA按照PrimeScript™ RT reagent Kit with gDNA Eraser (TaKaRa) 试剂盒说明书进行反转录，将获得的cDNA于–20 ℃保存备用。
1.2.2 OrFMRFamide基因的克隆
根据笔者实验室测得的转录组数据库 (未发表) 筛选得到目的基因的部分基因片段，使用Primer Premier 5.0软件设计2对引物，进行PCR扩增反应，反应产物用10 mg·mL−1的琼脂糖凝胶电泳检测，使用D2000 Ladder作为Maker，产物条带单一且清晰，送至苏州金唯智生物科技有限公司测序，测序结果在NCBI上BLAST比对，确认为FMRFamide基因片段，然后按照SMARTER® RACE5'/3'Kit及3'Full RACE Core Set With PrimerScriptTM RTase试剂盒的说明书在Primer Premier 5.0软件上设计5'和3' RACE特异性引物，通过RACE技术得到产物，产物经过1.0%琼脂糖凝胶电泳检测合格后送往苏州金唯智生物科技有限公司测序，将所得到的全部测序结果通过Sequebcher 5.0软件拼接得到OrFMRFamide基因的序列全长。
1.2.3 OrFMRFamide基因的序列分析
把拼接得到的基因序列全长在NCBI上用ORF Finder (https://www.ncbi.nlm.nih.gov/orffinder/) 在线工具对OrFMRFamide基因开放阅读框 (Open reading frame, ORF) 及编码氨基酸进行预测，得到预测的编码氨基酸后在NCBI上用在线工具BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi) 进行同源性比较，使用NetNGlyc 1.0 Server (http://www.cbs.dtu.dk/services/NetNGlyc/) 在线工具确定糖基化位点，理化性质用ProtParam (http://web.expasy.org/protparam/) 在线工具分析，使用NetPhos 2.0 Server (http://www.cbs.dtu.dk/services/NetPhos–2.0/) 在线工具查找磷酸化位点，使用Signal 5.0 (http://www.cbs.dtu.dk/services/SignalP/) 在线工具查找信号肽。用MEGA-X软件基于比邻法，以1 000次Bootstraps构建瘤背石磺与其他物种FMRFamide基因的系统进化树。
1.2.4 OrFMRFamide基因在不同组织中的表达
通过在线工具Genscrip (https://www.genscript.com/tools/real–time–pcr–taqman–primer–design–tool) 根据OrFMRFamide基因全长序列设计荧光定量引物 (表1)，以18S作为内参基因，以上述验证过的瘤背石磺各组织cDNA为模板，使用 ChamQ Universal SYBR qPCR Master Mix (Vazyme，美国) 试剂盒在CFX-96 (Bio-Rad，美国) 进行荧光定量反应，每个样品和内参都设置3个平行对照，以Rnase-free双蒸水 (ddH2O) 为阴性对照，实时荧光定量PCR数据用2–ΔΔCT法进行分析。
表 1 实验中所用引物序列Table 1. Primers used in this experiment引物
Primer引物序列 (5'–3')
Primer sequence用途
FunctionOrFMRFamide-F1 CTTAGGAGTGGGAACAGC 验证目的片段 OrFMRFamide-R1 CGGCTGGAGATTTGATTG OrFMRFamide-F2 CGGACCAGTACGACCAAC OrFMRFamide-R2 GTTCAGTCCGCCCTAATG OrFMRFamide-3'RACE-inter ACTGGTTTGGGTAGCA 3'RACE 特异性引物 OrFMRFamide-3'RACE-outer ATGGCAACAATGTCTTTCG OrFMRFamide-5'RACE-inter GAGCAGAAGATGGCGT 5'RACE 特异性引物 OrFMRFamide-5'RACE-outer TCCTCGCTTTGCCTCA OrFMRFamide-RT-F AGCTGGAGGACACACTGAGGCA Real-time RT-PCR OrFMRFamide-RT-R TGCCACATCGCCCTCATCGG 18S-F TCCGCAGGAGTTGCTTCGAT 18S-R ATTAAGCCGCAGGCTCCACT 1.2.5 免疫组织化学分析
免疫组织化学取养殖箱暂养未经处理的5只瘤背石磺进行实验。采集未经处理的新鲜的神经节、肝胰腺、性腺、腹足、肌肉和皮肤，用4%的多聚甲醛固定6 h，转移至20%的蔗糖溶液固定过夜，然后转移至40%的蔗糖溶液固定到组织沉底，使用冰冻切片机 (莱卡CM1950) 切片后贴在黏附载玻片上，按照一步法免疫组化试剂盒 (凯基) 说明书完成免疫组化实验，用显微镜观察拍照。
1.2.6 LPS刺激
在养殖箱中暂养7 d后随机选取100只瘤背石磺，平均分为2组，暂养在2个养殖箱中，分别给每只瘤背石磺注射100 μL磷酸盐缓冲溶液 (Phosphate buffered saline, PBS) 和100 μL脂多糖 (LPS, 0.5 mg·mL–1，溶于PBS)，经处理的瘤背石磺分别放回养殖箱暂养，分别在注射后第0、第6、第12、第24、第48小时取样，每次随机抽取9只石磺，将3只的血细胞、神经节、肝胰腺、肌肉、性腺、皮肤和腹足合并为1个样本，每个时间点采集3个平行样本。血细胞的采集是使用1 mL注射器从瘤背石磺腹部抽取血液，在800 ×g、4 ℃条件下离心10 min，获得血细胞后将其溶解在trizol中，用液氮速冻，于–80 ℃保存。其他组织取样后放入液氮速冻，于–80 ℃保存。用trizol法提取各组织RNA，用HiScript Q RT SuperMix for qPCR (+gDNA wiper, Vazyme，美国) 试剂盒进行反转录，得到的cDNA于–20 ℃备用。
1.2.7 OrFMRFamide基因在炎症刺激后的时空表达
按照上述方法，以各组织各个时期刺激前后的cDNA为模板，18S为内参，进行实时荧光定量RT-PCR扩增实验，数据使用2–ΔΔCT法进行分析。
1.2.8 统计学分析
使用SPSS 26.0软件对荧光定量的数据进行处理，单因素方差分析结果以“平均值±标准差 (
$ \overline X \pm {\rm{SD}} $ )”表示。P<0.05为差异显著，P<0.01为差异极显著。2. 结果
2.1 OrFMRFamide的cDNA序列特征
根据瘤背石磺转录组测序结果获得的OrFMRFamide基因片段，通过3'和5'-RACE扩增得到该基因全长序列为2 618 bp，开放阅读框 (Open reading frame, ORF)为882 bp，基因序列末端加尾信号AATAA与Poly (A) 之间有17个碱基 (图1)。ProtParam推测编码1个含有293个氨基酸的蛋白，分子量为33847.88 kD，理论等电点 (pI) 为9.25，不稳定指数为57.00，预测为不稳定蛋白，溶脂指数为55.67，平均疏水指数为–0.890。Server预测该蛋白不存在糖基位点。NetPhos 2.0 Server预测该蛋白含有12个磷酸化位点 [ 丝氨酸 (Ser)：4；苏氨酸 (Thr)：2；酪氨酸 (Tyr)：6]。信号肽软件预测该基因有1个N端信号肽 (1~24)。
![]()
图 1 FMRFamide基因的氨基酸序列及cDNA全长浅灰色阴影部分是该基因编码信号肽部分，起始密码子和终止密码子用方框标出，深灰色阴影是磷酸化位点，末端加尾信号AATAAA用下划线标出，poly (A) 用斜体表示Figure 1. Amino acid sequence and full length cDNA of FMRFamide geneThe light grey shaded portion is the signal peptide portion of the gene; the start and stop codons are boxed; the dark grey shaded portion is the phosphorylation site; the end plus tail signal AATAAA is underlined, and poly (A) is in italics.2.2 系统进化树分析
将预测得到的OrFMRFamide基因的氨基酸序列放到NCBI上进行BLAST比对，发现其与其他物种的FMRFamide基因有较高的同源性，下载这些物种的氨基酸序列，使用MEGA-X软件构建进化树 (图2)。结果显示，瘤背石磺和静水椎实螺 (Lymnaea stagnalis) 聚为一支，再与加利福尼亚海兔、棕蜗牛 (Cornu aspersum) 等腹足纲软体动物聚为一大支，与双壳纲海洋贝类聚为一大支，与环节动物和节肢动物各聚为一支。瘤背石磺与静水椎实螺进化关系最为接近。
![]()
图 2 利用MEGA-X软件基于NJ法构建的FMRFamide系统进化树Figure 2. NJ phylogenetic tree of FMRFamide by MEGA-X2.3 OrFMRFamide在不同组织中的表达分析
荧光定量RT-PCR结果显示，OrFMRFamide基因在神经节中的表达量极显著高于在其他组织的(P<0.01)，此外，该基因也在血细胞、肝胰腺、皮肤、性腺组织中表达，几乎不在肌肉和腹足中表达 (图3)。以性腺的表达量为参照，血细胞、神经节、肝胰腺和皮肤的表达量分别约为性腺的2.33、1 025.47、36.96和6.86倍。
![]()
图 3 OrFMRFamide在各组织的相对表达量**. 极显著性差异 (P<0.01)Figure 3. mRNA relative expression of OrFMRFamide in different tissues**. Very significant at 0.01 level (P<0.01)2.4 FMRFamide多肽免疫组织化学分析
利用免疫组织化学的方法分析瘤背石磺FMRFamide神经肽在神经节、肝胰腺、肌肉、性腺、皮肤和腹足中的分布情况 (图4)，结果显示OrFMRFamide基因在不同组织中的表达模式和FMRFamide神经肽的分布水平相一致，在神经节中大量表达，在性腺、皮肤和肝胰腺中微量表达，在腹足和肌肉中几乎不表达。
![]()
图 4 免疫组化结果a. 神经节；b. 腹足；c. 肝胰腺；d. 肌肉；e. 皮肤；f. 性腺；箭头表示FMRFamide多肽Figure 4. Immunohistochemistry resultsa. Ganglia; b. Pleopod; c. Hepatopancreas; d. Muscle; e. Skin; f. Gonad; arrows represent FMRFamide polypeptides.2.5 OrFMRFamide在LPS刺激后的时空表达分析
基于瘤背石磺OrFMRFamide基因在不同组织中的相对表达量，研究该基因在LPS刺激后48 h内的表达变化 (图5)。LPS组肝胰腺中OrFMRFamide基因的表达量在注射后第12小时达最高值，与PBS组差异显著 (P<0.05)，皮肤中OrFMRFamide基因的表达量随时间变化差异不大，但LPS组整体的表达量显著高于PBS对照组 (P<0.05)，神经节中该基因的表达量在第12 小时达最高值，整体表达水平显著高于PBS对照组 (P<0.05)，血细胞中该基因的相对表达量在第12小时达最高值，性腺中的表达量随时间变化的差异不显著，且LPS组和PBS组差异也不显著 (P>0.05)。
![]()
图 5 OrFMRFamide基因在炎症刺激后在不同组织中的相对表达量图中不同字母表示差异显著 (P<0.05)Figure 5. Relative expression of OrFMRFamide gene in different tissues after inflammatory stimulationDifferent letters in the figure indicate significant difference (P<0.05).3. 讨论
本研究克隆了瘤背石磺FMRFamide基因，系统进化树显示其与静水椎实螺、棕蜗牛聚为一支，说明其与软体动物门腹足纲贝类的亲缘关系非常近，与软体动物双壳纲贝类明显区分开来。
神经节、肝胰腺和血细胞属于NEI (神经节、肝胰腺、血细胞) 系统，有报道指出NEI系统在应激和感染过程中能维持机体内环境稳定，有效清除病原体，调节机体平衡，减少对宿主的伤害[15]。FMRFamide作为软体动物神经系统中的重要组成部分[16]，实时荧光定量结果显示OrFMRFamide在神经节、血细胞和肝胰腺组织中均有较高的表达，在神经节中表达量最高，可能因为其参与神经系统中神经递质的合成和释放及神经内分泌系统的可塑性调节。杨金龙等[17]在研究厚壳贻贝 (Mytilus coruscus) 早期幼虫发育时也发现，FMRFamide大量存在于神经系统中。在其他软体动物如双壳纲[7,18]、腹足纲[19]、头足纲[20]中，FMRFamide多肽也被报道主要存在于神经系统、消化系统和血细胞中。FMRFamide在瘤背石磺NEI系统中的高表达说明该基因可能对瘤背石磺的稳态起重要作用，推测FMRFamide可能在神经节、肝胰腺和血细胞中合成和加工，然后通过自分泌/旁分泌信号通路释放到血细胞中来调节机体。荧光定量结果显示瘤背石磺FMRFamide还存在于皮肤中，高凤娟等[21]和Allan等[22]对花背蟾蜍 (Bufo raddei) 和北美木蛙 (Rana sylvatica) 的研究表明，两栖类动物皮肤具有先天免疫功能，说明瘤背石磺作为无贝壳的两栖型贝类，为适应潮间带恶劣的生存环境，FMRFamide多肽可能参与了其先天免疫。FMRFamide应用免疫组化学进一步证实了OrFMRFamide在不同组织中的分布和FMRFamide多肽的分布一致，这与Terrnina等[23]对猫肝吸虫 (Opisthorchis felineus) 和Allan等[22]对星虫 (Themiste lageniformis) 中FMRFamide的定位结果相似。
炎症刺激实验结果显示，瘤背石磺神经节、肝胰腺和血细胞中FMRFamide基因相对表达量变化趋势为先上升后降低，在注射脂多糖12 h后的表达量达最大值，说明FMRFamide可能参与了免疫防御反应。在无脊椎动物中，血细胞被认为是宿主在免疫防御反应中起重要作用的免疫相关细胞[24]，更有研究表明免疫细胞可以通过自分泌信号途径重新合成神经肽来调节肌体免疫反应[25-27]。Wang等[28]用寄生虫感染光滑双脐螺 (Biomphalaria glabrata)，发现FMRFamide相对表达量高于对照组，FMRFamide表达量上升以抵抗寄生虫感染，这和本研究结果一致，可以推测FMRFamide多肽能保护机体免受侵害，响应免疫刺激，在维持机体稳态中发挥重要作用。炎症刺激后瘤背石磺性腺中FMRFamide相对表达量变化不显著，Hada等[29]研究发现FMRFamide还能影响根结线虫 (Meloidogyne graminicola) 的繁殖因子，Kerbl等[30]的研究也证实FMRFamide能作用于生殖腺，影响动物的交配行为，推测瘤背石磺性腺中的FMRFamide可能对瘤背石磺的繁殖有影响，但不参与维持自身稳态。
-
![]()
图 3 秋刀鱼鱼体所有地标点的平均形和叠加示意图
Figure 3. Diagrams of mean shape and superimposition of all landmarks on body of C. saira
![]()
图 4 不同性别秋刀鱼形态特征的第1和第2主成分散点示意图
注：实线和虚线分别代表雌性和雄性的置信区间。
Figure 4. Scatter plots of PC1 and PC2 of different sexes of C. saira
Note: Solid and dashed lines represent confidence intervals for females and males, respectively.
![]()
图 5 不同性别秋刀鱼的网格变形图和变异可视化(变异扩大3倍)
Figure 5. Grid deformation and variation visualization of different sexes of C. saira (Variation are enlarged 3 times)
表 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 
下载: 导出CSV
表 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 躯干前部最大弯曲处
下载: 导出CSV
表 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
下载: 导出CSV
表 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.
下载: 导出CSV
表 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
下载: 导出CSV
表 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
下载: 导出CSV
-
[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. -
期刊类型引用(1)
1. QIU Jiayin，SUN Lianlian，LI Shuang，ZHOU Xu，CHI Changfeng，ZHENG Libing. Identification and Expression Profile of a Neuropeptide LFRFamide-Like Gene During Different Stages of Gonadal Development in the Cephalopod Sepia pharaonis. Journal of Ocean University of China. 2024(02): 499-508 . 
必应学术
其他类型引用(0)
计量
- 文章访问数: 161
- HTML全文浏览量: 54
- PDF下载量: 49
- 被引次数: 1
粤公网安备 44010502001741号
