Sensitivity analysis of CMSY based on varying parameter priorassumptions
-
摘要:
受数据不足的限制，大部分渔业难以采用传统的渔业资源评估方法为渔业管理提供建议。因此，基于数据缺乏方法的评估模型逐渐受到关注，其中渔获量-最大可持续产量 (Catch-maximum sustainable yield, CMSY) 是目前国际上应用广泛的数据缺乏评估方法。然而，现有研究表明，CMSY模型的评估结果对先验假设具有高度依赖性，且不同参数先验对模型输出结果的具体影响尚不明确。为探究CMSY参数先验假设对评估结果的影响，从RAM传统资源评估数据库 (RAM Legacy Stock Assessment Database) 中随机选取200个来自不同海域的鱼类及无脊椎动物种群对CMSY敏感性进行分析。系统分析了内禀增长率 (r) 以及开始、中间、最后年份的资源量水平 (Bstart/K、Bint/K、Bend/K，B为生物量，K为环境容纳量) 的上下限先验假设变化对模型估算结果生物量轨迹、最大可持续产量 (MSY)、最大可持续产量对应的生物 (BMSY)、最大可持续产量对应的捕捞死亡率 (FMSY)、相对生物量 (B/BMSY)、相对捕捞死亡率 (F/FMSY)]的影响。结果表明：1) 生物量轨迹、BMSY、FMSY受r下限的影响大，B/BMSY、F/FMSY受Bend/K上限的影响大，而MSY受各参数先验的影响较小。2) 目级分类差异对CMSY方法的敏感性无显著性影响。3) 参数先验的改变易导致Kobe图中的种群资源状况发生象限迁移，其中Bend/K上限对其影响最大。研究表明，CMSY模型对参数先验设置具有较高的敏感性，建议使用该模型时谨慎设定参数和进行结果分析。
Abstract:Due to data limitation, most fisheries, most fisheries find it difficult to use traditional fisheries resource assessment methods to provide recommendations for fisheries management. As a result, data-limited assessment models have gainedincreasing attention, among which catch-maximum sustainable yield (CMSY) has become a widely used method internationally. However, studies indicate that CMSY is highly dependent on prior assumptions, and the specific impacts of different parameter priors on model outputs remain uncertain. To address this problem, we conducted a sensitivity analysis of CMSY by randomly selecting 200 fish and invertebrate populations from different sea areas from the RAM Legacy Stock Assessment Database. The model estimation results, including biomass trajectory, MSY, BMSY, FMSY, B/BMSY and F/FMSY were systematically analyzed to examine their variations by incrementally adjusting the upper and lower bounds of the intrinsic growth rate (r), and the biomass levels at the start, middle and end years (Bstart/K, Bint/K, Bend/K). The results indicate that: 1) Biomass trajectory, BMSY, and FMSY were significantly influenced by the lower bound of r, while B/BMSY and F/FMSY were more sensitive to the upper bound of Bend/K. MSY was less affected by the prior settings of the parameters. 2) Taxonomic classification at the order level exhibited no significant influence on the sensitivity of the CMSY method. 3) Variations in parameter priors were observed to induce quadrant transitions in the population stock status depicted in the Kobe plot, among which the upper bound of Bend/K showed the most significant effect. This study reveals that the CMSY model exhibited high sensitivity to parameter prior settings. It is suggested that careful parameter configuration and result interpretation are essential when utilizing this model.
-
Keywords:
- CMSY /
- Sensitivity analysis /
- Stock assessment /
- Data-limited
-
无乳链球菌 (Streptococcus agalactiae) 也称B族链球菌 (Group B Streptococcus, GBS)，是一种人、鱼共患的革兰氏阳性菌[1]。与许多其他致病菌一样，无乳链球菌具有许多毒力因子可提高病原体感染或破坏宿主的能力[2]，其毒力因子可通过影响宿主细胞的黏附和侵袭以逃避宿主免疫[3]，从而有效削弱机体免疫系统和新陈代谢过程中诱发的生理变化[4]，严重时会降低机体的先天防御力，进而诱导机体死亡。作为一种由病原菌引起的常见疾病，链球菌病在罗非鱼养殖中具有破坏性影响[5]。2019年，无乳链球菌和海豚链球菌 (S. iniae) 给罗非鱼养殖业造成的经济损失约达25亿元[6]。此外，无乳链球菌感染还可导致斑马鱼 (Danio rerio)[7]在内的多种鱼体患病并大量死亡。硬骨鱼感染无乳链球菌的典型症状包括鱼体表变黑、眼球突出或混浊、发白、出血、腹部斑点和鳃盖内侧出血等[8]。目前已有学者从菌株的分离鉴定和分型、防治链球菌感染的药物筛选、疫苗研制等方面进行了无乳链球菌感染鱼体的研究报道[9]。关于无乳链球菌感染卵形鲳鲹 (Trachinotus ovatus) 的发病机理和调节机制方面的研究仍有不足。
卵形鲳鲹作为我国南方深远海养殖的主要品种，年产量近20万吨[10]。近年来，随着养殖环境变化以及养殖密度增加，导致卵形鲳鲹病害频发，给养殖产业造成严重的经济损失。Cai等[11]2016年首次报道了无乳链球菌感染卵形鲳鲹，感染后可导致其大量死亡，死亡率最高可达每天2%[12]。有研究报道，无乳链球菌侵袭生物体主要通过血流和主要器官在全身传播。作为鱼类的重要免疫调节器官，肝脏、脾脏和肾脏是病原菌入侵鱼体时的广泛研究对象。例如，罗非鱼感染无乳链球菌后可导致其肝脏、肠、鳃、脾、头肾和脑发生不同程度的病变，其中在肾脏和脾脏中检测到较高的细菌密度[5]。杂交蛇头鱼 (Channa maculata♀ × C. argus♂) 脾脏中存在的免疫和凋亡相关通路在舒伯特气单胞菌 (Aeromonas schubert) 感染中发挥重要作用[13]。爱德华氏菌 (Edwardsiella) 感染可诱导许氏平鲉 (Sebastes schlegeli) 脾、肝和头肾中钙网蛋白的显著上调，从而减少病原体在其体内的传播和复制[14]。然而，对卵形鲳鲹脾脏在无乳链球菌入侵后的调控机制研究仍有不足。
因此，本研究在已有研究基础上，通过酶活性测定、组织病理学观察和实时定量PCR探究了卵形鲳鲹脾脏在无乳链球菌感染后的生理生化反应和凋亡诱导机制，以期为卵形鲳鲹抗病育种及长期健康养殖提供参考依据。
1. 材料与方法
1.1 实验鱼和细菌
实验所用卵形鲳鲹均来自中国水产科学研究院南海水产研究所深圳试验基地，平均体质量为31.15 g。实验开始前，选取无病原感染、健壮的500尾卵形鲳鲹，转入与实验条件一致的实验桶进行暂养，养殖水温为 (27±0.5) ℃，盐度25‰，溶解氧质量浓度保持在5.5 mg·L−1以上。每天投喂2次，饲喂量约占体质量的4%。选用的无乳链球菌菌株2021年分离自中国水产科学研究院南海水产研究所深圳试验基地的患病卵形鲳鲹。
1.2 感染预实验
实验开始前，将细菌接种于BHI液体培养基中，并在180 r·min−1、27 ℃的摇床上孵育24 h。将BHI液体培养物在8 000 r·min−1[15]的条件下离心8 min，收集沉淀物。然后用无菌磷酸盐缓冲溶液 (PBS) 洗涤沉淀4次，用不同浓度 (1.0×1010、1.0×109、1.0×108、1.0×107和1.0×106 CFU·mL−1) 的无乳链球菌注射感染卵形鲳鲹。感染后，观察鱼的行为变化，记录不同感染浓度下死亡率与时间的关系，根据死亡率最终得到半致死浓度为1.0×108 CFU·mL−1。
1.3 人工感染实验
感染实验在150 L的水族箱中进行，水体140 L。选取300尾健康的卵形鲳鲹，随机平均分为6组，其中3个处理组，3个对照组。按照预实验结果，将鱼麻醉后处理组每尾鱼腹腔注射200 μL无乳链球菌 (1.0×108 CFU·mL−1)，对照组注射相同剂量的无菌PBS溶液。
注射后的第0、第6、第12、第24、第48、第72、第96和第120小时取样，每个时间点取3尾鱼，用40 mg·L−1丁香酚麻醉。鱼体经体积分数为75%乙醇消毒后，用1.5 mL无菌注射器静脉取血，将3尾鱼的血液混合置于2.0 mL离心管中静置5 h，离心 (1 000 r·min−1、20 min) 分离获得血清，置于1.5 mL冻存管中，于−80 ℃冰箱中保存用于分析酶活性等相关指标。另外用肝素钠 (100 IU·mL−1) 润洗后的离心管分别收集第0和第120小时血液用于血常规检测。取新鲜的脾脏组织液氮速冻后置于−80 ℃冰箱中备用；另取新鲜脾脏组织通过生理盐水清洗表面血液后用4%多聚甲醛固定用于组织学检查。
1.4 病理观察
对鱼样品进行宏观和组织学检查，记录外部和内部形态变化。将固定在多聚甲醛中的脾脏组织用体积分数为70%的乙醇清洗，经过脱水后使用常规技术进行石蜡包埋、切片，切片厚度5 μm，并用苏木精-伊红染色，使用正置荧光显微镜 (NIKON ECLIPSE C1) 获得切片图像[16]。
1.5 血常规指标检测
通过全自动血液细胞分析 (深圳迈瑞生物医疗电子股份有限公司，BC-5000Vet) 检测血液中白细胞 (WBC)、嗜中性粒细胞 (NE)、淋巴细胞 (LYM)、红细胞 (RBC)、血小板 (PLT)、血红蛋白浓度 (HGB)、红细胞比容 (HCT)、平均红细胞体积 (MCV)、平均红细胞血红蛋白浓度 (MCHC)、平均血小板体积 (MPV) 和血小板压积 (PCT) 的变化。
1.6 生化指标检测
对照组和实验组卵形鲳鲹血清中乳酸脱氢酶 (LDH)，脾脏中酸性磷酸酶 (ACP) 和碱性磷酸酶 (ALP) 的活性采用北京华英生物技术研究所试剂盒和酶标仪 (华卫德朗DR-200BS) 进行测定。
1.7 Caspase基因表达
按照RNA提取试剂盒 (广州美基生物科技有限公司) 说明书提取卵形鲳鲹脾脏总RNA，质量分数为1%琼脂糖凝胶电泳和NanoDrop 2000 (Thermo Fisher，美国) 检测其质量和浓度，使用 PrimeScript™ RT试剂盒和gDNA Eraser合成cDNA，于−20 ℃冰箱中储存备用。
使用Primer Premier 5设计Caspase-3、Caspase-8和Caspase-9的引物序列 (表1)，由于EF-1α不受无乳链球菌感染的影响，因此选择EF-1α作为内参基因。通过实时荧光定量PCR仪 (Roche Light Cycler® 480 II，罗氏诊断产品有限公司，上海) 进行实时定量PCR (qRT-PCR)，反应体系为12.5 μL。反应条件为95 ℃预变性30 s；95 ℃变性5 s，60 ℃退火30 s，72 ℃延伸 30 s，共40个循环。每个样品重复3次实验，在获得每个样品的阈值之后，使用2−ΔΔt方法计算Caspase-3、Caspase-8和Caspase-9 mRNA的相对表达水平，使用SPSS 20.0软件进行统计分析。数据以“平均值±标准差 (
$\overline { X}\pm { \rm {SD}} $ )”表示。所有统计分析均使用单因素方差分析 (One-way ANOVA) 进行比较，P<0.05表示差异显著。表 1 引物信息及序列Table 1. Primers and sequences information applied in this study引物名称
Primer name引物序列 (5'—3')
Primer sequence (5'–3')应用
ApplicationCaspase-3-F GCTGCTCTACTGCTTCTGCCTGATG qRT-PCR Caspase-3-R TGGCTGAGGATTGTGATGTTGCTG Caspase-8-F GCAACAAAACAGCCATCCA qRT-PCR Caspase-8-R GCAGGGGTAAAGGGTCATT Caspase-9-F GAATGGCGTCCGTCTGGTCATC qRT-PCR Caspase-9-R GGCAGCACGTCTCAGTTCAGC EF-1α-F AAGCCAGGTATGGTTGTCAACTTT qRT-PCR EF-1α-R CGTGGTGCATCTCCACAGACT 2. 结果
2.1 病理学观察
人工感染无乳链球菌24 h后，卵形鲳鲹游动不定、食欲不振、嗜睡、运动不协调，且明显观察到角膜混浊、突眼、眼出血、内脏充血和内鳃盖充血等典型病变 (图1)，并在感染后48 h出现大量死亡。组织病理学分析显示，脾脏发生较为严重的病变损伤，典型特征为白髓中B淋巴细胞和T淋巴细胞的坏死和丢失；轻度至中度病变表现为脾炎，白髓区域显著减少，淋巴细胞数目降低。在严重病变中，脾脏白髓弥漫性坏死至白髓区域完全消失，伴有网状内皮增生，可见大量被HE染成棕黄色的“小结”(图2)。未感染的卵形鲳鲹脾脏具有正常的红髓和白髓。
![]()
图 1 卵形鲳鲹感染无乳链球菌后内部和外部病变注：a. 角膜混浊、突眼、眼出血；b. 内脏充血；c. 内鳃盖充血 (箭头所指)。Figure 1. Internal and external lesions recorded on T. ovatus specimens affected by S. agalactiaNote: a. Corneal opacity, exophthalmia and eye hemorrhage; b. Gut congestion; c. Congestion in the inner gill cap (Arrow).![]()
图 2 卵形鲳鲹脾脏组织病理学特征注：a. 健康脾脏组织切片；b—c. 感染无乳链球菌120 h 后脾脏组织切片。Figure 2. Histopathological characteristics of spleen of T. ovatusNote: a. Healthy spleen tissue section; b−c. Spleen tissue section 120 h after infection with S. agalactiae.2.2 血液指标参数
卵形鲳鲹感染无乳链球菌后血液参数见表2。与对照组相比，感染个体的红细胞数量、血红蛋白浓度和红细胞比容明显降低 (P<0.05)，白细胞、嗜中性粒细胞、淋巴细胞数量则出现相反的趋势 (P<0.05)，与对照组相比嗜中性粒细胞和淋巴细胞数量增加近3倍。同时平均红细胞体积明显增大，其他指标未见明显变化。
表 2 卵形鲳鲹感染无乳链球菌对血液学指标的影响Table 2. Effect of S. agalactiae infection on hematological parameters of T. ovatus项目
Item对照组
Control感染组
Infection白细胞数量 WBC/(109·L−1) 8.32±0.94 13.67±0.58* 嗜中性粒细胞数量 NE/(109·L−1) 0.18±0.09 0.54±0.12* 淋巴细胞数量 LYM/(109·L−1) 1.01±0.35 3.04±0.75* 红细胞数量 RBC/(109·L−1) 5.25±1.07 3.08±0.25* 血红蛋白浓度 HGB/(g·L−1) 179.67±10.25 129.58±9.58* 红细胞比容 HCT/% 35.69±0.55 23.37±0.17* 平均红细胞体积 MCV/fL 160.63±1.27 190.31±0.79* 平均红细胞血红蛋白量
MCH/pg39.75±2.47 41.11±3.73 平均红细胞血红蛋白质量浓度
MCHC/(g·L−1)289.21±2.99 292.36±4.13 血小板数量 PLT/(109·L−1) 19.68±0.25 13.43±0.86 平均血小板体积 MPV/fL 7.91±1.39 7.41±1.05 血小板压积 PCT/% 0.70±0.17 0.75±0.11 注：*. 与对照组相比存在显著性差异 (P<0.05)。 Note: *. Significant difference compared with the control group (P<0.05). 2.3 血清、脾脏生化指标结果
卵形鲳鲹感染无乳链球菌后的血清和脾脏生化指标变化见图3，相较于对照组，感染组个体的血清LDH和脾脏中ALP、ACP活性显著增加 (P<0.05)，并于第6小时达到峰值 (P<0.01)；同时随着感染时间延长，其活性均逐渐降低但ACP与对照组相比仍有极显著差异 (P<0.01)。其中LDH活性在感染6 h后升高最为显著 (P<0.01)，活性约为对照组的10倍。
![]()
图 3 卵形鲳鲹感染无乳链球菌后生化指标变化注：a. 不同时间血浆乳酸脱氢酶变化水平；b. 不同时间脾脏碱性磷酸酶变化水平；c. 不同时间脾脏酸性磷酸酶变化水平；**. 差异极显著 (P<0.01)；*. 差异显著 (P<0.05)。Figure 3. Changes of biochemical indexes after infection of T. ovatus at different timeNote: a. The levels of plasma lactate dehydrogenase at different time before and after the challenge; b. The levels of spleen alkaline phosphatase at different time before and after the challenge; c. The levels of spleen acid phosphatase at different time before and after the challenge; **. Very significant difference (P<0.01); *. Significant difference (P<0.05).2.4 脾脏中凋亡基因表达
卵形鲳鲹感染无乳链球菌后凋亡基因表达谱见图4。在第0—第120小时的实验期内，Caspase-3、Caspase-8和Caspase-9的mRNA表达水平随时间推移均不断升高。其中Caspase-8基因表达水平在感染第6小时显著高于对照组 (P<0.05)，并在第120小时达到峰值 (P<0.01)；Caspase-3在感染后第120小时达到峰值；Caspase-9在感染后第12小时表达量显著升高 (P<0.05)，并在第72小时达到峰值后趋于稳定。
![]()
图 4 卵形鲳鲹感染无乳链球菌后脾脏凋亡基因相对表达量注：**. 差异极显著(P<0.01)；*. 差异显著 (P<0.05)。Figure 4. Relative expression of apoptosis genes in spleen after T. ovatus infection with S. agalactiaeNote: **. Very significant difference (P<0.01); *. Significant difference (P<0.05).3. 讨论
本研究通过无乳链球菌体外注射实验感染鱼体后，自患病鱼中分离得到的菌株经形态分析和16S RNA测序比对后均与无乳链球菌具有高度一致性，因此确定其为无乳链球菌感染。经腹腔注射200 μL浓度为1.0×108 CFU·mL−1的病原菌后，在48 h后实验组死亡率开始增加，这与罗非鱼人工注射相似剂量无乳链球菌后的死亡高峰期接近[17]。观察发现实验组卵形鲳鲹出现与自然感染无乳链球菌后相似的临床症状[18]，组织病理学分析显示，卵形鲳鲹在无乳链球菌感染后表现出脾脏炎症病变并伴有含铁血黄素沉积。有研究表明，机体对铁元素利用受阻导致脾脏中含铁血黄素的积累或红细胞被大量破坏是鱼类对环境应激的重要信号[19]，因此，脾脏中含铁血黄素的检测对卵形鲳鲹的健康监测具有重要意义。本研究中，脾脏作为卵形鲳鲹的主要淋巴器官在无乳链球菌感染后损伤严重，造成无法区分红髓和白髓以及淋巴细胞减少，同时可见产生大量炎症细胞，这与罗非鱼在自然状态下感染无乳链球菌后的组织切片结果相似[3]。
血清生化参数是组织或器官在临床病理学研究中的重要指标[20]。研究表明，高于或低于基线水平的酶活性通常反映组织或器官代谢紊乱，严重时会引起细胞凋亡[21]。细胞凋亡是感染的最终结果，主要包括细胞坏死、凋亡和自噬[22]。LDH作为糖酵解途径中的末端酶，在辅酶NADH协助下可将乳酸转化为丙酮酸。当机体受到外界压力时，LDH可以进行逆反应并将丙酮酸 (糖酵解的最终产物) 转化为乳酸[23]。而细胞内LDH的释放是细胞坏死的特征之一。本研究中感染无乳链球菌后卵形鲳鲹血清中LDH显著升高，表明机体可能存在细胞凋亡情况并伴有病理性坏死反应，机体通过凋亡过程清除体内坏死的细胞并参与病原体入侵后的免疫调控过程。此外，ALP和ACP作为溶酶体中的主要水解酶参与生物体的先天免疫过程，参与包括生长和细胞分化等一系列生理代谢活动，常被作为病原菌感染的敏感指标，在评估鱼类病原菌引发的疾病中起重要作用[23-25]。Sharkoori等[26]研究指出ALP和ACP活性可作为细胞坏死的指标。本研究发现，无乳链球菌感染可使卵形鲳鲹脾脏中ALP和ACP活性在短时间内显著升高，这可能是由于脾脏组织功能活性增加促使其在短时间内大量合成，以降低病原菌感染对机体的损害[27]。此外，本研究中ACP活性变化与南亚野鲮 (Labeo rohita)感染嗜水气单胞菌 (A. hydrophila) 后的变化相似[28]，推测可能是机体内溶酶体在抵抗无乳链球菌感染过程中大量增殖，促使溶酶体膜破裂，致使膜内保持潜伏状态的ACP被激活。
血常规检验对于病原体感染的监测具有重要意义[29]。病原菌感染鱼体后可通过细胞膜扩散到血液系统，并对机体血液氧气输送能力和血液电解质平衡产生负面影响[30]。血液中红细胞比容、红细胞和血红蛋白浓度是评估鱼类在暴露于各种环境压力、化学毒性和细菌感染后健康状况的重要指标[25]。本研究中，实验组红细胞比容、红细胞和血红蛋白浓度均显著低于对照组，与草鱼 (Ctenopharyngodon idella)[31]人工感染草鱼呼肠孤病毒 (Grass carp reovirus) 后的变化相似，因此认为血液生化特性可用作卵形鲳鲹对无乳链球菌感染的敏感指标。此外，有研究发现，红细胞比容、红细胞和血红蛋白浓度降低会造成贫血，这可能是细菌感染动物后抑制红细胞生成或是细菌感染导致红细胞生命周期缩短而引起贫血状态[32]。同时本研究中实验组平均红细胞体积较大，推测其可能伴有红细胞增多症。以往研究发现，由于鱼类与哺乳动物不同，幼体红细胞比成体红细胞小，细胞质较少，当外源微生物入侵后可使鱼平均红细胞体积增大[33]。此外，有研究报道，血液中白细胞、嗜中性粒细胞和淋巴细胞数量常因病原菌感染而增多[34]。本研究发现，受感染鱼白细胞、嗜中性粒细胞和淋巴细胞的含量显著升高，可能是病原菌感染诱导机体产生可逆性造血干细胞损伤后出现的侵袭性炎症病变或血液早期再增殖，机体在试图通过发展获得性免疫而增强对病原体的抗性[35]。
当先天免疫被破坏导致应激反应受阻时，细胞会遭受严重损伤并启动凋亡或坏死。细胞凋亡是由基因编码的自杀程序的作用引起的，该程序触发了一系列特征性的形态和生化变化，Caspase是鱼类检测细胞凋亡的关键指标[36]。细胞凋亡通常通过外在途径 (受体-凋亡途径) 和内在途径 (线粒体-凋亡途径) 两种主要途径发生[37]，其中外在途径通过细胞外信号调节死亡受体，募集相关死亡结构域和Caspase-8相关的蛋白质，随后Caspase-8的激活直接或间接地激活Caspase-3。内在途径是由细胞色素C与凋亡蛋白酶激活因子 (Apaf-1) 结合，激活下游效应子Caspase-9后进一步激活Caspase-3[38]。河豚 (Takifugu obscurus)[39]感染嗜水气单胞菌后头肾中Caspase-3和Caspase-8的mRNA表达水平显著上调，促使细胞凋亡。此外，已有研究报道无乳链球菌诱导的细胞凋亡可以改变水生动物凋亡相关基因的表达水平[40]。本研究发现，无乳链球菌可诱导卵形鲳鲹脾脏中Caspase-3、Caspase-8和Caspase-9基因的转录水平显著增加，表明机体的内外凋亡途径可同时参与对病原体入侵的免疫调控过程。
-
![]()
图 1 改变参数上限导致Kobe图中种群资源状况发生象限移动的示例
注：PCEE、SWHI和ATKA三个种群代码来源于RAM数据库，分别代表小孔羽鼬鳚 (Genypterus blacodes)、弗氏鱚 (Sillago flindersi) 和单鳍多线鱼 (Pleurogrammus monopterygius)。黑色圆点的-0编号代表其初始参数设置情景，黑色三角形的-8编号代表增加Bend/K参数上限的情景。箭头代表了这些种群先验参数改变后资源状况移动的方向。
Figure 1. Example of quadrant shift in stock status in Kobe plot caused by changing parameter upper limit
Note: The stock codes PCEE, SWHI, and ATKA originated from the RAM database and represented Genypterus blacodes, Sillago flindersi and Pleurogrammus monopterygius, respectively. The black circles with a -0 designation represent the initial parameter setting scenarios, while the black triangles with a -8 designation represent scenarios with an increased upper limit for Bend/K parameters. The arrows represent the directional changes in the stock status from adjustments to the prior parameter settings.
![]()
图 2 每个先验参数增加10%时生物量轨迹(a)、MSY (b)、BMSY (c)、FMSY (d)、B/BMSY (e)、F/FMSY (f)与参考值的差异
注：箱型图内部灰色点为均值；图中未显示离群点。横轴坐标：1. r下限；2. r上限；3. 开始年份资源量水平下限；4. 开始年份资源量水平上限；5. 中间年份资源量水平下限；6. 中间年份资源量水平上限；7. 最后年份资源量水平下限；8. 最后年份资源量水平上限。“*”代表显著性差异，p<0.05；“ns”代表非显著性差异，p≥0.05。黑色虚线为0值参考线。
Figure 2. Difference in biomass trajectory (a), MSY (b), BMSY (c), FMSY (d), B/BMSY (e), F/FMSY (f) relative to reference values when each prior parameter increased by 10%
Note: Grey points inside the box plot are means; outliers are not shown in the plot. The horizontal coordinates are: 1. r-lower; 2. r-upper; 3. Bstart/K-lower; 4. Bstart/K-upper; 5. Bint/K-lower; 6. Bint/K-upper; 7. Bend/K-lower; 8. Bend/K-upper. "*" represents a significant difference with p< 0.05; "ns" represents a non-significant difference with p≥ 0.05. The black dashed line is the zero value reference line.
![]()
图 3 不同目的各个种群参数先验增加10%时结果 (生物量轨迹、MSY、BMSY、FMSY、B/BMSY、F/FMSY) 与参考值的差异
注：横坐标中：1. 生物量轨迹；2. 最大可持续产量(MSY)；3. 最大可持续产量对应的生物量(BMSY)；4. 最大可持续产量对应的捕捞死亡系数(FMSY)；5. 相对生物量(B/BMSY)；6. 相对捕捞死亡系数(F/FMSY)。蓝色虚线为0值参考线。
Figure 3. Difference in results (Biomass trajectory, MSY, BMSY, FMSY, B/BMSY, F/FMSY) relative to reference values when prior parameters for stocks from different orders increased by 10%
Noet: The horizontal coordinates are: 1. Biomass trajectory; 2. MSY; 3. BMSY; 4. FMSY; 5. B/BMSY; 6. F/FMSY. The blue dashed line is the zero value.
表 1 基于恢复力分类的r参数先验范围
Table 1 Prior range of r parameter based on resilience classification
恢复力
Resiliencer先验范围
Prior range of r高 High 0.600~1.500 中 Medium 0.200~0.800 低 Low 0.050~0.500 非常低 Very low 0.015~0.100 
下载: 导出CSV
表 2 58个种群中不同参数先验改变对种群资源状况的影响
Table 2 Impact of prior parameter changes on population stock status
先验参数
Prior parameterB/BMSY F/FMSY 向右越过基准线
(个数)
Rightward shift
across baseline
(Number)向左越过基准线
(个数)
Leftward shift
across baseline
(Number)合计
Sum向上越过基准线
(个数)
Upward shift
across baseline
(Number)向下越过基准线
(个数)
Downward shift
across baseline
(Number)合计
Sumr-下限 r-lower 6 1 7 3 10 13 r-上限 r-upper 2 1 3 6 9 15 Bstart/K-下限 Bstart/K-lower 4 2 6 2 8 10 Bstart/K-上限 Bstart/K-upper 4 2 6 3 7 10 Bint/K-下限 Bint/K-lower 3 2 5 4 10 14 Bint/K-上限 Bint/K-upper 5 2 7 2 11 13 Bend/K-下限 Bend/K- lower 4 2 6 4 12 16 Bend/K-上限 Bend/K- upper 30 1 31 1 19 20 合计 Sum 58 13 71 25 86 111
下载: 导出CSV
-
[1] 石永闯, 樊伟, 张衡, 等. 适用于数据缺乏渔业的资源评估方法研究进展[J]. 中国水产科学, 2021, 28(5): 673-691. doi: 10.12264/JFSC2020-0309 [2] COSTELLO C, OVANDO D. Status, institutions, and prospects for global capture fisheries[J]. Annu Rev Environ Resour, 2019, 44(1): 177-200. doi: 10.1146/annurev-environ-101718-033310
[3] 石永闯. 西北太平洋秋刀鱼 (Cololabis saira) 资源评估研究[D]. 上海: 上海海洋大学, 2020: 48. [4] COSTELLO C, OVANDO D, HILBORN R, et al. Status and solutions for the world’s unassessed fisheries[J]. Science, 2012, 338(6106): 517-520. doi: 10.1126/science.1223389
[5] ROSENBERG A A, FOGARTY M J, COOPER A B, et al. Developing new approaches to global stock status assessment and fishery production potential of the seas[R]. Rome: FAO, 2014: 0_1, I, IV, X, XI, 1, 3, 5-25, 27-45, 47-49, 51-81, 83-89, 91-103, 105-175.
[6] FROESE R, DEMIREL N, CORO G, et al. Estimating fisheriesreference points from catch and resilience[J]. Fish Fish, 2017, 18(3): 506-526. doi: 10.1111/faf.12190
[7] VASCONCELLOS M, COCHRANE K. Overview of world status of data-limited fisheries: inferences from landings statistics[C]//Anchorage: University of Alaska Fairbanks. Fisheries assessment and management in data-limited situations, 2005: 1-20.
[8] THORSON J T, MINTO C, MINTE-VERA C V, et al. A new role for effort dynamics in the theory of harvested populations and data-poor stock assessment[J]. Can J Fish Aquat Sci, 2013, 70(12): 1829-1844. doi: 10.1139/cjfas-2013-0280
[9] MACCALL A D. Depletion-corrected average catch: a simple formula for estimating sustainable yields in data-poor situations[J]. ICES J Mar Sci, 2009, 66(10): 2267-2271. doi: 10.1093/icesjms/fsp209
[10] DICK E, MACCALL A D. Depletion-based stock reduction analysis: a catch-based method for determining sustainable yields for data-poor fish stocks[J]. Fish Res, 2011, 110(2): 331-341. doi: 10.1016/j.fishres.2011.05.007
[11] COPE J M. Implementing a statistical catch-at-age model (Stock Synthesis) as a tool for deriving overfishing limits in data-limited situations[J]. Fish Res, 2013, 142: 3-14. doi: 10.1016/j.fishres.2012.03.006
[12] PONS M, COPE J M, KELL L T. Comparing performance of catch-based and length-based stock assessment methods in data-limited fisheries[J]. Can J Fish Aquat Sci, 2020, 77(6): 1026-1037. doi: 10.1139/cjfas-2019-0276
[13] CHEUNG W W, PALACIOS-ABRANTES J, FRÖLICHER T L, et al. Rebuilding fish biomass for the world's marine ecoregionsunder climate change[J]. Glob Change Biol, 2022, 28(21): 6254-6267. doi: 10.1111/gcb.16368
[14] FROESE R, WINKER H, CORO G, et al. Status and rebuilding of European fisheries[J]. Mar Pol, 2018, 93: 159-170. doi: 10.1016/j.marpol.2018.04.018
[15] ZHANG K, ZHANG J, XU Y W, et al. Application of a catch-based method for stock assessment of three important fisheries in the East China Sea[J]. Acta Oceanol Sin, 2018, 37: 102-109.
[16] FROESE R, WINKER H, CORO G, et al. New developments in the analysis of catch time series as the basis for fish stock assessments: the CMSY plus plus method[J]. Acta Ichthyol Piscat, 2023, 53: 173-189. doi: 10.3897/aiep.53.e105910
[17] OVANDO D, FREE C M, JENSEN O P, et al. A history and evaluation of catch-only stock assessment models[J]. Fish Fish, 2022, 23(3): 616-630. doi: 10.1111/faf.12637
[18] LIANG C, XIAN W W, PAULY D. Assessments of 15 exploited fish stocks in Chinese, South Korean and Japanese waters using the CMSY and BSM methods[J]. Front Mar Sci, 2020, 7: 623. doi: 10.3389/fmars.2020.00623
[19] FALSONE F, SCANNELLA D, GERACI M L, et al. How fishery collapses: the case of Lepidopus caudatus (Pisces: Trichiuridae) in the Strait of Sicily (Central Mediterranean)[J]. Front Mar Sci, 2021, 7: 584601. doi: 10.3389/fmars.2020.584601
[20] BOUCH P, MINTO C, REID D G, et al. Comparative performance of data-poor CMSY and data-moderate SPiCT stock assessment methods when applied to data-rich, real-world stocks[J]. ICES J Mar Sci, 2021, 78(1): 264-276. doi: 10.1093/icesjms/fsaa220
[21] 陈镝安, 朱江峰, 耿喆. 应用3种数据有限方法对印度洋近海金枪鱼类和马鲛类资源评估[J]. 上海海洋大学学报, 2024, 33(3): 728-740. [22] FROESE R, DEMIREL N, CORO G, et al. User guide for CMSY++[R]. Kiel: GEOMAR, 2021: 1-17.
[23] RICARD D, MINTO C, JENSEN O P, et al. Examining the knowledge base and status of commercially exploited marine species with the RAM Legacy Stock Assessment Database[J]. Fish Fish, 2012, 13(4): 380-398. doi: 10.1111/j.1467-2979.2011.00435.x
[24] 莫奇龙, 胡世森, 谢君君, 等. 南海北部海鳗资源利用状况评估[J]. 渔业信息与战略, 2022, 37(1): 27-33. [25] MAUNDER M N, AIRE-DA-SILVA A. Evaluation of the kobe plot and strategy matrix and their application to tuna in the EPO[R]. La Jolla, USA: IATTC, 2012.
[26] DE BRUYN P. Report of the 2017 ICCAT shortfin mako assessment meeting[R]. Madrid: ICCAT, 2017.
[27] IOTC SECRETARIAT. Report of the 26th session of the IOTC working party on tropical tunas[R]. Seychelles: IOTC, 2024.
[28] FROESE R, PALOMARES M L D, PAULY D. Estimation of life history key facts[M]//FROESE R, PAULY D. FishBase 2000: Concepts, Design and Data Sources. Philippines: ICLARM, 2000.
[29] 石永闯, 陈云云, 樊伟, 等. 基于Catch-MSY模型的西北太平洋日本鲭资源评估及敏感性分析[J]. 海洋渔业, 2023, 45(1): 26-38. doi: 10.3969/j.issn.1004-2490.2023.01.003 [30] BOX G E P. Robustness in the strategy of scientific model building[M]//Robustness in statistics. Newyork: Academic Press, 1979: 201-236.
[31] WINKER H, CARVALHO F, SHARMA R, et al. Initial results for North and South Atlantic shortfin mako (Isurus oxyrinchus) stock assessments using the Bayesian surplus production model JABBA and the catch-resilience method CMSY[R]. Madrid: ICCAT, 2017: 1836-1866.
[32] FITZGERALD C J, DELANTY K, SHEPHARD S. Inland fish stock assessment: Applying data-poor methods from marine systems[J]. Fish Manag Ecol, 2018, 25(4): 240-252. doi: 10.1111/fme.12284
[33] FREE C M, JENSEN O P, ANDERSON S C, et al. Blood from a stone: performance of catch-only methods in estimating stock biomass status[J]. Fish Res, 2020, 223: 105452. doi: 10.1016/j.fishres.2019.105452
-
期刊类型引用(12)
1. 徐志华，任娣，刘崇万，刘熠，葛筱琴，胡月，朱晓华. 江苏某地区稻虾蟹综合种养模式下养殖水和产品中农药残留与风险评价. 农业环境科学学报. 2024(01): 37-47 . 
百度学术
2. 杨秋红，皇甫鑫，郑瑜，董靖，艾晓辉. 湖北省稻虾养殖模式下小龙虾农药残留与风险评估. 华中农业大学学报. 2024(06): 297-306 . 百度学术
3. 吴松，黄晓丽，陶月，覃东立，陈中祥，郝其睿，高磊，王海涛，孙言春，王鹏，张颖. 氯虫苯甲酰胺对中华绒螯蟹的急性毒性研究. 水产学杂志. 2023(03): 46-51 . 百度学术
4. 曾斯丽，高磊，王鹏，陈中祥，黄晓丽，覃东立. 淡水环境及水产品中酰胺类除草剂残留现状研究进展. 中国渔业质量与标准. 2023(03): 43-52 . 百度学术
5. 黄晓丽，高磊，汤施展，黄丽，霍堂斌，王鹏，陈中祥，覃东立. 东北地区养殖淡水鱼有机氯和拟除虫菊酯农药残留及健康风险评价. 中国渔业质量与标准. 2022(02): 1-9 . 百度学术
6. 卜媛媛，何小燕，李翱，金华，金诗语，罗思. 稻虾综合种养模式下小龙虾的膳食风险讨论. 食品安全导刊. 2022(09): 130-133 . 百度学术
7. 黄丽，高磊，郝其睿，吴松，王鹏，陈中祥，黄晓丽，覃东立. 哈尔滨郊区池塘养殖鱼类六六六和滴滴涕残留特征及食用健康风险. 食品安全质量检测学报. 2022(15): 4867-4875 . 百度学术
8. 姚晶晶，吴旋，李葳，王明锐，崔文文，张惠贤，易甜. 湖北省某地区不同种植模式下土壤与稻米中重金属含量和污染风险评价. 湖北农业科学. 2022(22): 146-151 . 百度学术
9. 黄晓丽，吴松，黄丽，王鹏，陈中祥，覃东立. 苯噻酰草胺胁迫对中华绒螯蟹幼蟹的毒性效应. 中国渔业质量与标准. 2021(04): 1-8 . 百度学术
10. 彭中校，张华威，黄会，宫向红，张秀珍. 异丙甲草胺对水生生物的毒性效应研究进展. 中国渔业质量与标准. 2021(04): 62-68 . 百度学术
11. 周军，周刚，李旭光，邓燕飞，许郑超，陆全平. 中华绒螯蟹质量安全风险研究. 中国渔业质量与标准. 2020(01): 13-28 . 百度学术
12. 唐建军，李巍，吕修涛，王岳钧，丁雪燕，蒋军，汤亚斌，李坚明，张金保，杜军，游宇，李晓东，李斌，成永旭，窦志，高辉，陈欣. 中国稻渔综合种养产业的发展现状与若干思考. 中国稻米. 2020(05): 1-10 . 百度学术
其他类型引用(5)
计量
- 文章访问数: 18
- HTML全文浏览量: 4
- PDF下载量: 1
- 被引次数: 17
粤公网安备 44010502001741号
