Spatio-temporal variation of Ommastrephes bartramii resources (winter & spring groups) in Northwest Pacific under different climate modes
-
摘要:
柔鱼 (Ommastrephes bartramii) 是大洋性短生命周期物种，其生活史和生物量受环境和气候因子影响明显。基于2004—2015年西北太平洋柔鱼渔捞日志、海表面温度 (Sea surface temperature，SST)、叶绿素 (Chlorophyll a，Chl a) 浓度及太平洋年代际震荡 (Pacific decadal oscillation，PDO) 数据，利用空间自相关统计方法、热点分析和小波分析法研究PDO冷期与暖期两种气候模态下西北太平洋柔鱼渔场的时空变化。结果表明，PDO指数 (Pacific Decadal Oscillation Index，PDOI) 与柔鱼单位捕捞努力量渔获量 (Catch per unit effort，CPUE) 呈正相关，且CPUE滞后PDO 8个月 (R2=0.548，P<0.05)，CPUE与PDO指数的共轭周期为2~4个月。不同气候模态下的渔场热点分布特征为，暖期时渔场重心向高纬方向分布较为明显，空间上的集聚性强；冷期时向低纬分布较为明显，空间上的集聚性相对较差。PDO冷暖时期对西北柔鱼时空分布影响显著，该研究对柔鱼资源的可持续开发具有一定的科学意义。
Abstract:Ommastrephes bartramii is an oceanic species with short life cycle, whose life history and biomass are affected by environmental and climatic factors significantly. Based on the fishery data, sea surface temperature (SST), chlorophyll a (Chl a) and Pacific decadal oscillation (PDO) from 2004 to 2015 in the Northwest Pacific, we applied patial autocorrelation statistics, hot spot analysis and wavelet analysis to study the spatio-temporal variation of O. bartramii in that area under two PDO models with cold and warm periods. The results show that the PDO index was positively correlated with the catch per unit effort (CPUE) and the CPUE lagged PDO for eight months. However, the conjugate period was 2−4 months under the influence of CPUE and PDO index. The distribution of fishery hotspots under two PDO modes is characterized by the fact that the gravity centers of fishing ground is obviously distributed in high and low latitude direction with strong and weak spatial agglomeration in warm and cold periods, respectively. PDO index affects the spatio-temporal distribution of O. bartramii in different periods significantly, which has positive implications for the sustainable development of O. bartramii in the Northwest Pacific Ocean.
-
中华绒螯蟹 (Eriocheir sinensis) 简称河蟹，属于甲壳纲、十足目、方蟹科、绒螯蟹属，主要分布于东亚地区[1]。由于其味道鲜美、营养均衡，广受消费者喜爱[2-3]。近几年，中国年养殖总产量近80×104 t，产值超过500×108元[4]。然而，产业高速发展也导致一系列问题出现，例如水系之间的盲目引种和随意杂交，亲本规格的降低和近亲交配，养殖密度的不断增加和病害的频发等，造成长江水系部分中华绒螯蟹养殖群体种质退化[5]。遗传选育是一种重要的种质改良方案。现有的5个中华绒螯蟹新品种[6] (长江1号、长江2号、江海21号、光合1号和诺亚1号) 和长荡湖1号等新品系[7]均主要是以生长速度和规格为指标进行群体选育而产生的，在一定程度上推动了中华绒螯蟹养殖产业的可持续发展。然而，越来越多的消费者开始追求差异化的特色产品，现有的种质资源尚不能满足中华绒螯蟹养殖产业对良种经济性状多样化的需求。
壳色是甲壳动物重要的经济性状，直接影响着消费者的喜好和市场价格[8]。本课题组于2017年以壳色为指标，在国际上首次选育出一种白壳色中华绒螯蟹新品系，暂定为“海大玉蟹”，迄今已经繁育至G3代。现有研究表明，白壳新品系营养均衡且具有生长速度快的优势 (另文发表)，对其进行规模化人工繁育和种质推广可能具有重要意义。生殖性能和胚胎质量评估是开展虾蟹类遗传选育工作的首要步骤，生殖性能评估主要包括亲本抱卵速度、抱卵率、亲本成活率等，而胚胎质量主要为卵径、卵质量、色泽、类胡萝卜素含量以及生化组成等指标[9-12]。以往对中华绒螯蟹生殖性能和胚胎质量的研究主要集中在不同亲本水系来源[10-11]、不同亲本饲料添加[13-14]、亲本营养强化[15-16]、亲本成熟时间[17]和亲本规格[18]，有关壳色的研究仅见于紫壳中华绒螯蟹[19] (色泽不可遗传)。研究表明，紫壳雌蟹亲本的成活率明显高于绿壳雌蟹，虽然紫壳亲本抱卵率略低，但其胚胎中花生四烯酸 (C20:4 n-6, ARA) 和二十二碳五烯酸(C22:5 n-3, EPA) 含量显著高于绿壳[19]。中华绒螯蟹白色壳品系的生殖性能与正常绿壳亲本间是否存在一定差异，至今尚不清楚，这在一定程度上，阻碍了白壳新品系中华绒螯蟹的种质资源开发和推广。鉴于此，本研究测定了中华绒螯蟹白壳新品系和绿壳品系亲本生殖性能与卵质量、常规生化、脂肪酸和类胡萝卜素组成，以期为河蟹基于壳色的多样化育种提供理论基础和实践参考。
1. 材料与方法
1.1 亲本蟹的来源
中华绒螯蟹白壳品系为本课题组2017年初通过家系选育获得。2020年11月初，从上海海洋大学崇明蟹类育种基地室外同一土池中挑选完成生殖蜕壳、体质量相近、肢体完整和体质健壮的白壳奇数年G2代品系和正常绿壳品系 (雌蟹各50只和雄蟹各20只，白壳和绿壳品系雌体亲本的外观照片见图1)。将两种壳色蟹在崇明基地室内循环水槽中按照雌雄分开原则进行暂养，每日17:00投喂适量的冰冻杂鱼和鲜活缢蛏 (Sinonovacula constricta)。
![]()
图 1 中华绒螯蟹白壳 (a) 和绿壳 (b) 品系亲本雌体外观照片Figure 1. Photos of white (a) and green (b) carapace stains of female E. sinensis broodstock1.2 亲本交配
2020年1月14日，在上海海洋大学崇明蟹类育种基地的室内循环水系统中实施亲本交配。白壳和绿壳组各用3个规格和尺寸一致的黑色PVC桶 (Φ108×H120 cm) 养殖，每个桶底部铺设一层约10 cm的细沙供蟹抱卵，同时投放4片PVC塑料瓦片作为遮蔽物。每个桶投放雌蟹10只和雄蟹3只 (生物学指数见表1，同一性别两种壳色各指标间无显著差异)。海水为人工配置的半咸水，盐度17，水深80 cm。每日17:00对室内所有亲本蟹投喂适量冰冻杂鱼和鲜活缢蛏，每日9:00和21:00记录室内交配桶内水温并检查抱卵蟹抱卵和死亡情况。交配7~10 d后剔除雄蟹以防影响雌蟹抱卵。雌体亲本的成活率 (Survival rate, RS, %) 和抱卵率 (Ovigerous rate, RO, %) 计算公式为：
表 1 白壳和绿壳雌体和雄体亲本的生物学性状Table 1. Biological characteristics of white and green carapace strains of E. sinensis broodstock类目 Item 雌体 Female (N=30) 雄体 Male (N=9) 白壳 White carapace 绿壳 Green carapace 白壳 White carapace 绿壳 Green carapace 体质量
Body mass/g161.85±12.39 152.45±10.13 222.57±3.31 215.65±12.09 甲壳长
Carapace length/cm62.12±1.44 60.52±0.70 67.34±0.56 66.63±0.96 甲壳宽
Carapace width/cm68.51±1.45 67.52±1.14 74.79±1.79 75.07±1.16 $$ R_{\mathrm{S}}=\left(S_{\mathrm{i}} / S_{\mathrm{t}}\right) \times 100 {\text{%}} $$ (1) $$ R_{\mathrm{O}}=\left(O_{\mathrm{i}} / O_{\mathrm{t}}\right) \times 100 {\text{%}}$$ (2) 式中：Si和St分别为存活雌体亲本个体数和投放雌体亲本总数；Oi和Ot分别为雌体抱卵蟹数量和交配后雌体存活总亲本数量。
1.3 亲本生殖性能测定
当胚胎发育至原肠胚期，从白壳和绿壳组各取6只抱卵蟹 (每个桶各取2只)，参照李清清等[10]方法测定亲本生殖性能，具体操作为用干毛巾将体表擦干并用镊子将腹部卵块全部取出，剔除携卵的刚毛和绒毛。用海水清洗卵块并用纸巾吸干卵块水分，对吸干后的卵块和去卵后的抱卵蟹精确称质量。每个样品取10 mg胚胎装入称量杯中，采用百万分之一电子天平 (赛多利斯CPA225D，0.001 mg，德国) 精确称质量并用计数器计数，据此计算单卵湿质量。每个蟹卵块随机测定3次，取平均值为该蟹单卵湿质量。随后将卵块装袋密封后置于−80 ℃冰箱保存。计算抱卵量 (Egg production, EP, 104 粒·只−1)、生殖力 (Fecundity, F, 粒·g −1) 和生殖指数 (Reproductive effort, RE, %)，计算公式为：
$$ \mathrm{EP}=W_{\mathrm{e}} / W_{\mathrm{i}} $$ (3) $$ F=\mathrm{EP} / W_{0} $$ (4) $$ \mathrm{RE}=W_{\mathrm{e}} / W_{\mathrm{t}} $$ (5) 式中：We和Wi分别表示卵块质量 (g) 和单卵湿质量 (g)；W0为抱卵蟹净体质量 (去卵后质量，g)；Wt为抱卵蟹体质量 (g)。
1.4 卵径、单卵湿质量和干质量测定
采用ToupView显微图像测量分析软件在显微镜下测定单个卵的卵径，每个抱卵蟹随机测定60粒卵，取平均值为该蟹的卵径。参考李清清等[10]方法进行单卵湿质量和干质量测定，其中，单卵湿质量与生殖性能中测定方法一致；单卵干质量具体操作为取称质量后的湿卵，用计数器精确计数后放入锡箔称量杯中，用百万分之一天平准确称量总质量，然后放置于65 ℃烘箱中烘至恒质量，再对烘干后的卵块和锡箔纸总质量准确称量以计算单卵干质量。每只抱卵蟹随机取样3次进行测定，最终取平均值为该抱卵蟹的单卵干质量。
1.5 胚胎的色泽和类胡萝卜素测定
每个抱卵蟹取少量鲜卵装入透明自封袋，用Lovibond-RT200色差仪测定鲜卵色泽。干卵色泽测定之前，将冻干后的蟹卵粉碎过60目筛，随机取样倒入小培养皿 (直径50 mm) 中，用色差仪测量表面层色度。每个抱卵蟹胚胎鲜样和冻干样均随机取6个点测量，平均值为该蟹卵的最终色泽参数，具体参数为亮度 (L*)、红度 (a*) 和黄度 (b*)，具体参考李清清等[10]的方法。
每个抱卵蟹取0.3~0.4 g冻干后的卵进行总类胡萝卜素的提取和测定，主要参考孙秋凤等[11]方法：采用80%丙酮对样品进行首次提取，用超声波在遮光环境下浸提15~20 min，离心取出上清液，剩余样品继续加入100%丙酮 (C3H6O)，反复提取3~4次，直至提取的上清液基本无色，将提取液用100%丙酮定容至15 mL后，使用分光光度计在474 nm波长下测定上清液的光密度 (OD)。采用100%丙酮稀释虾青素标准品至不同浓度梯度并测定对应梯度吸光值，制备标准曲线，计算各样品中总类胡萝卜素含量。样品经 0.22 μm过滤器过滤后进行高效液相色谱分析。采用安捷伦1260高效液相色谱系统 (配备YMC™类胡萝卜素C30柱，4.6×150 mm，包装直径=3 μm，YMC) 对提取液进行主要类胡萝卜素含量测定和分析。参数设置为流速1.0 mL·min−1，柱温30 ℃，进样量20 μL。样品中类胡萝卜素种类的鉴定依据虾青素、叶黄素、玉米黄素、角黄素和β-胡萝卜素5种商业标准品出峰时间，样品中类胡萝卜素含量根据5种标准品制备的标准曲线和样品出峰面积计算。
1.6 胚胎的常规生化和脂肪酸测定
称量0.2 g左右湿胚胎于105 ℃烘干至恒质量测定水分，按照AOAC[20]标准方法即凯氏定氮法测定冻干后胚胎中粗蛋白，总脂测定采用经典的氯仿/甲醇 (V/V=2/1) 抽提方法[21]。脂肪酸的测定参考Wu等[22]的方法，采用三氟化硼-甲醇法对提取的总脂进行甲酯化，用正己烷溶解装入2 mL的棕色进样瓶，使用安捷伦 (GC-7890B, MS-5977A) 气相色谱仪 (色谱柱为Omegawax-320毛细管柱，30 m×0.25 mm，美国Supelco公司) 进行脂肪酸测定和分析。脂肪酸质量分数的计算采用面积百分比法。
1.7 数据分析
数据以“平均值±标准误 (
$\overline { X}\pm { \rm {SE}} $ )”表示，统计分析采用SPSS 19.0软件，作图采用Excel 2019软件，通过Levene法对实验数据进行方差齐性检验，两组数据之间采用独立样本t检验进行显著性分析，P<0.05为差异显著。2. 结果
2.1 生殖性能
白壳和绿壳品系雌体亲本交配后第4天均出现抱卵蟹，交配后4~16 d，白壳雌体的抱卵率均显著高于绿壳雌体 (P<0.05，图2)，这说明白壳品系交配后抱卵速度较快，两品系雌体的最终抱卵率无显著性差异 (P>0.05)。抱卵蟹养殖期间，白壳品系雌体的成活率为93.33%，绿壳为86.67%，但差异不显著 (P>0.05，表2)。白壳品系雌体的抱卵量、生殖力和生殖指数与绿壳品系较为接近，两者无显著性差异 (P>0.05)。白壳和绿壳胚胎卵径几乎一致，分别为369.99和370.04 μm，无显著性差异 (P>0.05，表3)。白壳组单卵湿质量为45.41 μg，绿壳组为42.92 μg，白壳组略高于绿壳组，但差异不显著 (P>0.05)。白壳组单卵干质量为12.06 μg，绿壳组为11.42 μg，白壳组也略高于绿壳组，但差异不显著 (P>0.05)。
![]()
图 2 白壳和绿壳中华绒螯蟹雌蟹抱卵率及水温变化Figure 2. Ovigerous rate and water temperature of white and green carapace strains of E. sinensis broodstock表 2 中华绒螯蟹白壳和绿壳品系雌蟹亲本成活率、抱卵率和生殖性能比较Table 2. Survival rate, ovigerous and reproductive performance of white and green carapace strains of E. sinensis broodstock类目
Item样本数
Sample number白壳
White carapace绿壳
Green carapace成活率 survival rate/% 30 93.33±6.67 86.67±8.82 抱卵率 Ovigerous rate/% 30 92.50±3.82 83.74±6.18 抱卵量 Egg production/(104 粒·只−1) 6 43.33±8.38 44.88±3.73 生殖力 Fecundity/(粒·g−1) 6 3 083.68±429.67 3 361.96±150.18 生殖指数 Reproductive effort/% 6 11.71±1.12 12.57±0.40 表 3 中华绒螯蟹白壳和绿壳品系卵径和单卵质量比较 (N=6)Table 3. Comparison of egg diameter and egg mass of white and green carapace strains of E. sinensis broodstock类目
Item白壳
White carapace绿壳
Green carapace卵径 Egg diameter/μm 369.99±3.82 370.04±2.76 单卵湿质量 Egg wet mass/μg 45.41±1.65 42.92±1.10 单卵干质量 Egg dry mass/μg 12.06±0.67 11.42±0.35 2.2 胚胎色泽和类胡萝卜素质量分数
白壳组胚胎鲜样a*显著高于绿壳组 (P<0.05)，L*和b*也高于绿壳组，但差异不显著(P>0.05，表4)；而干样色泽参数中，白壳组L*略高于绿壳组，a*和b*略低于绿壳组，但差异均不显著 (P>0.05)，与外观肉眼观察结果类似 (图3)。两组壳色胚胎胡萝卜素组成中，白壳组虾青素和叶黄素质量分数显著低于绿壳组 (P<0.05)，但两组间总类胡萝卜素、玉米黄素、角黄素以及β-胡萝卜素质量分数无显著性差异 (P>0.05)。
表 4 中华绒螯蟹白壳和绿壳品系雌体胚胎色泽参数和总类胡萝卜素质量分数比较 (湿质量，N=6)Table 4. Color parameters and carotenoids content of eggs produced by white and green carapace strains of E. sinensis broodstock (wet mass, N=6) mg·g−1类目
Item白壳
White carapace绿壳
Green carapace湿样 Fresh 亮度 L* 33.17±0.91 30.69±0.83 红度 a* 5.14±0.36a 3.67±0.42b 黄度 b* 5.72±0.77 4.07±0.47 干样 Dried 亮度 L* 53.79±1.19 52.19±0.88 红度 a* 14.96±0.84 16.80±1.20 黄度 b* 18.63±1.04 19.81±1.61 总类胡萝卜素 Total carotenoids 291.41±40.45 400.47±60.36 虾青素 Astaxanthin 11.78±1.45a 19.29±2.19b 叶黄素 Lutein 16.53±2.21a 30.05±4.77b 玉米黄素 Zeaxanthin 3.43±0.60 3.61±0.98 角黄素 Canthaxanthin 3.02±0.15 2.71±0.22 β-胡萝卜素 β-carotene 22.56±4.10 44.55±9.77 注：同行标有不同小写字母表示组间差异显著 (P<0.05)，表5同此。 Note: Values with different small letters within the same row have significant difference (P<0.05). The same case in the Table 5. ![]()
图 3 白壳和绿壳中华绒螯蟹原肠胚期蟹卵冻干后比较Figure 3. Comparison of dried egg in gastrula stage produced by white and green carapace strains of E. sinensis broodstock2.3 胚胎中常规生化和脂肪酸组成
就常规生化而言，白壳和绿壳组胚胎中水分质量分数无显著性差异 (P>0.05，表5)；白壳组胚胎的粗脂肪质量分数 (5.58%) 和绿壳组 (5.66%) 接近，无显著性差异 (P>0.05)；白壳组粗蛋白质量分数 (6.24%) 略高于绿壳组 (5.96%)，但无显著性差异 (P>0.05)。就脂肪酸而言，白壳和绿壳组胚胎中的总饱和脂肪酸 (ΣSFA) 和总单不饱和脂肪酸 (ΣMUFA) 以及总多不饱和脂肪酸 (ΣPUFA) 质量分数无显著差异 (P>0.05，表5)，但是白壳组顺-9-十八碳一烯酸 (C18:1 n-9) 质量分数显著高于绿壳组 (P<0.05)。虽然白壳组中ARA略高于绿壳组，EPA 和二十二碳六烯酸 (C22:6 n-3, DHA) 略低于绿壳组，但差异均不显著 (P>0.05)。白壳组和绿壳组n-3长链多不饱和脂肪酸/n-6长链多不饱和脂肪酸 (n-3/n-6) 比例分别为1.18和1.49，无显著性差异 (P>0.05)。
表 5 中华绒螯蟹白壳和绿壳品系雌体卵常规生化和脂肪酸组成比较 (湿质量)Table 5. Proximate composition and fatty acid composition of eggs produced by white and green E. sinensis broodstock (wet mass) %类目
Item白壳
White carapace绿壳
Green carapace类目
Item白壳
White carapace绿壳
Green carapace水分 Moisture 74.88±0.46 76.45±0.52 总单饱和脂肪酸 ΣMUFA 30.33±1.40 29.67±0.79 脂肪 Lipid 5.58±0.15 5.66±0.18 亚油酸 C18:2 n-6 14.21±1.56 12.18±2.00 蛋白 Protein 6.24±0.08 5.96±0.15 α-亚麻酸 C18:3 n-3 4.21±0.54 4.38±0.30 肉豆蔻酸 C14:0 0.69±0.05 0.96±0.16 二十碳二烯酸 C20:2 n-6 2.26±0.20 1.81±0.18 十五碳酸 C15:0 0.43±0.04 0.5±0.02 花生四烯酸 C20:4 n-6 3.85±0.46 3.45±0.29 棕榈酸 C16:0 10.12±0.22 10.01±0.11 二十碳三烯酸 C20:3 n-3 0.28±0.02 0.41±0.05 珠光脂酸 C17:0 0.73±0.05 0.91±0.07 二十碳五烯酸 C20:5 n-3 (EPA) 7.05±0.44 7.74±0.64 硬脂酸 C18:0 8.03±0.72 8.24±0.43 二十二碳五烯酸 C22:5 n-3 0.82±0.04 0.75±0.04 花生酸 C20:0 0.49±0.02 0.44±0.12 二十二碳六烯酸 C22:6 n-3 (DHA) 8.50±1.51 9.33±1.52 总饱和脂肪酸 ΣSFA 20.48±0.69 21.05±0.23 总多不饱和脂肪酸 ΣPUFA 41.18±0.90 40.05±0.67 肉豆蔻烯酸 C14:1 n-5 0.08±0.01 0.09±0.01 总长链多不饱和脂肪酸 ΣLC-PUFA 22.76±1.39 23.50±1.84 棕榈油酸 C16:1 n-7 6.24±0.84 6.39±0.59 总n-3长链多不饱和脂肪酸 Σn-3 PUFA 21.14±1.34 23.02±1.88 十七碳一烯酸 C17:1 n-7 0.93±0.07 1.09±0.02 总n-6长链多不饱和脂肪酸 Σn-6 PUFA 20.32±2.09 17.44±2.37 顺-9-十八碳一烯酸 C18:1 n-9 18.57±0.51a 16.28±0.55b n-3长链多不饱和脂肪酸/n-6长链
多不饱和脂肪酸 n-3/n-61.18±0.28 1.49±0.31 顺-7-十八碳一烯酸 C18:1 n-7 2.26±0.50 2.07±0.54 n-3长链多不饱和脂肪酸/n-6长链
多不饱和脂肪酸 n-3/n-61.18±0.28 1.49±0.31 二十碳一烯酸 C20:1 n-9 1.70±0.13 2.12±0.26 二十二碳六烯酸/二十碳五烯酸
DHA/EPA1.20±0.17 1.18±0.11 芥酸 C22:1 n-9 0.54±0.13 1.05±0.20 3. 讨论
3.1 两种壳色品系河蟹的生殖性能和胚胎质量
甲壳动物亲本的生殖性能主要受到亲本营养状况、交配环境、遗传选育和抱卵频次等影响[23-26]。本研究中，亲本均来自于室外同一养殖池塘，亲本营养状况接近，暂养和交配环境也相似，且均为第一次抱卵。所以排除了亲本营养状况、配对环境和抱卵次数因素，一定程度可以反映不同壳色与生殖性能之间的关系。本研究中白壳组亲本成活率和抱卵率明显高于绿壳组，但差异不显著，在其他条件基本一致的前提下，这种差异很可能是遗传因素导致的[27-28]。Li等[19]研究发现，紫壳与绿壳抱卵蟹的成活率和抱卵率也存在一定差异，但不显著，结合本研究结果可说明壳色对抱卵蟹的成活率和抱卵率存在一定影响。抱卵量是评估甲壳动物生殖性能的重要指标之一，中华绒螯蟹的抱卵量随着规格的增大而增加[18]。本研究中，白壳组亲本平均体质量略高于绿壳组，白壳品系的平均抱卵量反而略低于绿壳组，这可能是由于白壳品系雌体的抱卵量差异较大造成的。本研究中白壳品系亲本的生殖力和生殖指数略低于绿壳组，可能是因为白壳品系雌体平均体质量略大于绿壳品系。有研究表明，中华绒螯蟹的生殖力和生殖指数均随着雌体体质量的增加而下降[18]。本研究中不同壳色品系的胚胎卵径、单卵湿质量和干质量均无显著差异，这说明壳色品系对其胚胎参数影响较小；也有研究表明，中华绒螯蟹胚胎的卵径和单卵质量相对稳定，受亲本营养状况、色泽和亲本规格的影响较小[15,18-19]。值得一提的是，本研究发现白壳蟹抱卵速度明显优于绿壳蟹，笔者推测可能是绿壳蟹在黑色养殖桶中更容易伪装，而白壳外观呈鲜亮的白色与黑色养殖桶存在强烈的颜色反差，使得白壳蟹在水中更容易被异性发现且相遇，从而交配机会更大且频次更多，所以抱卵速度会更快。甲壳动物的体表颜色对交配机会的影响在椰子蟹 (Birgus latro)[29]、招潮蟹 (Austruca mjoebergi)[30]、相手蟹 (Parasesarma peninsulare)[31]和蓝蟹 (Callinectes sapidus)[32]中均有相关报道，不同物种影响程度也存在一定差异。
3.2 两种壳色品系河蟹的胚胎色泽和类胡萝卜素含量
类胡萝卜素有利于甲壳动物胚胎发育、幼虫早期色素沉积和眼柄发育[33-34]。有研究表明，中华绒螯蟹卵巢、肝胰腺和蟹壳中，a*与虾青素含量高度相关[35]。本研究中，白壳组胚胎干样中a*和b*明显低于绿壳组，这主要是因为其胚胎中总类胡萝卜素、虾青素、叶黄素以及玉米黄素含量也低于绿壳组，尤其虾青素和叶黄素含量显著低于绿壳组。推测白壳和绿壳组胚胎中虾青素和叶黄素含量存在显著差异可能有2个原因。1) 甲壳动物由于缺乏从头合成类胡萝卜素的能力，所以其体内的类胡萝卜素主要从食物中获取[34,36]。本研究中白壳和绿壳亲本来自同一室外养殖池塘，白壳和绿壳亲本可能对饵料中类胡萝卜素的吸收和保留效率不同，从而导致胚胎中类胡萝卜素含量存在差异。2) 尽管甲壳动物自身缺乏从头合成类胡萝卜素的能力，但具有一定的将体内β-胡萝卜素转化为虾青素的能力，不同物种和品系的这种转化能力可能存在一定差异[37-38]。本课题组的研究表明，白壳品系的体色具有可遗传性，且其蟹壳中的虾青素含量显著低于绿壳品系 (另文发表)；本研究结果表明白壳品系胚胎中的虾青素含量也低于绿壳品系，这暗示白壳品系将β-胡萝卜素转化为虾青素的能力可能弱于绿壳品系。
3.3 两种壳色品系河蟹的胚胎常规生化和脂肪酸组成
脂类和蛋白质是甲壳动物胚胎发育过程中重要的能源物质[39-40]。因此，对胚胎生化组成的研究有利于了解个体或组织对能源物质的储存和利用情况[40-41]。研究表明，甲壳动物胚胎中的营养物质主要来自亲本成熟卵母细胞，因此亲本的营养状况直接影响着胚胎中的生化组成[41]。本研究中，白壳和绿壳品系胚胎中水分、脂肪和蛋白含量无显著性差异，说明他们胚胎中积累的常规营养成分含量接近。Li等[19]研究表明，中华绒螯蟹紫壳和绿壳品系胚胎中水分质量分数约为61.11%~63.33%，脂肪质量分数约为10.95%~11.08%，蛋白质量分数约为18.20%~18.56%；本研究中白壳和绿壳品系胚胎中的水分质量分数高于前者，但是脂肪和蛋白质量分数低于前者，这主要因为本研究中的胚胎样品处于原肠胚期，而Li等[19]研究中的胚胎样品处于囊胚期。甲壳动物胚胎发育过程中，胚胎中水分含量一般呈逐渐上升趋势，而脂肪和蛋白呈下降趋势[40,42-44]。
脂肪酸在甲壳动物胚胎发育中扮演着重要角色，饱和脂肪酸和单不饱和脂肪酸是重要的能量型脂肪酸，高度不饱和脂肪酸 (HUFA) 主要影响甲壳动物繁殖性能、胚胎发育和幼体质量[14,45-47]。本研究中，白壳组中C18:1 n-9百分含量显著高于绿壳组 (P<0.05)，说明白壳组亲本中可能积累更多的C18:1 n-9，造成这一差异的具体原因有待进一步研究，可能与其胚胎中的类胡萝卜素含量不同有关。有研究发现亲本饵料和卵中的HUFA质量分数与卵的孵化率、胚胎发育速度和初孵幼体质量相关，在适宜质量分数条件下，胚胎中的HUFA含量越高，胚胎质量越好[45-46]。也有研究发现，亲本饲料中HUFA质量分数在一定范围内 (0.14%、0.65%、1.15%、1.64%和2.07%) 对亲本繁殖性能影响不显著[48]。本研究中，白壳和绿壳胚胎中HUFA组成和质量分数无显著差异，这提示白壳和绿壳组胚胎孵化率和初孵幼体质量可能接近。
综上所述，中华绒螯蟹白壳品系和常规绿壳品系具有比较相似的生殖性能和胚胎质量，但白壳品系交配后的抱卵速度显著高于常规绿壳品系。白壳品系胚胎中的虾青素和叶黄素含量显著低于绿壳品系，但其C18:1 n-9含量显著高于绿壳品系。这些研究结果可为白壳品系中华绒螯蟹的遗传选育、生殖性能评价和推广应用提供基础资料。
-
![]()
图 1 2004—2015年7—11月的月平均单位捕捞努力量渔获量和太平洋年代际震荡指数时间序列
Figure 1. Time series of monthly average CPUE and PDO index in July to November during 2004−2015
![]()
图 2 太平洋年代际震荡指数与单位捕捞努力量渔获量交相关分析
当滞后性为负值时，太平洋年代际震荡指数滞后单位捕捞努力量渔获量；当滞后性为正值时，单位捕捞努力量渔获量滞后于太平洋年代际震荡指数；蓝色虚线表示95%置信区间
Figure 2. Cross correlation function study on PDOI and CPUE
When time lag is negative, PDO index lags CPUE; when the time lag is positive, PDO index leads CPUE. ACF is automatic correlationfunction; blue dotted line indicates 95% confidence interval.
![]()
图 3 柔鱼单位捕捞努力量渔获量与渔场环境因子同太平洋年代际震荡指数的交叉小波功率谱
白色实线圈闭合区域为功率谱值通过置信水平为95%的标准背景谱检验，锥形范围表示有效值；箭头表示位关系，→表示单位捕捞努力量渔获量/海洋表面温度/叶绿素与太平洋年代际震荡同位相，说明两者为正相关关系；↗、↙表示单位捕捞努力量渔获量/海表温度/叶绿素滞后太平洋年代际震荡且分别为正相关关系、负相关关系，↘、↖表示单位捕捞努力量渔获量/海洋表面温度/叶绿素超前太平洋年代际震荡且分别为正相关关系、负相关关系
Figure 3. Cross-wavelet power spectrum of O. bartramii CPUE and fishery environmental factors and PDO index
The white solid coil closed area is the standard background spectrum test with a power spectrum value of 95% confidence level, and the cone range represents the effective value; arrow indicates bit relationship: ‘→’ indicates that the CPUE/SST/Chl a is in phase with PDO, indicating that the both are positively correlated; ‘↗, ↙’ indicate that CPUE/SST/Chl a lags PDO and is positively correlated and negatively correlated; ‘↘, ↖’indicate that CPUE/SST/Chl a leads PDO and is positively correlated and negatively correlated.
![]()
图 4 太平洋年代际震荡冷期与暖期柔鱼渔场热点重心位置分布图
a, b 分别为柔鱼资源2004—2015年的热点和冷点的空间分布；c. 太平洋年代际震荡暖期与冷期柔鱼空间热点分布重心在空间上的变化方向和距离；d. 太平洋年代际震荡冷期与暖期柔鱼资源每个月的热点重心位置在空间上的变化方向与距离
Figure 4. Distribution of center of gravity of hot fishing ground in cold and warm periods of PDO
a–b. The spatial distribution of hot and cold spots of squid resources from 2004 to 2015, respectively; c. The direction and distance of spatial hot spot distribution center of PDO in warm and cold periods; d. The spatial variation direction and distance of the hot spots of soft fish resources in cold and warm periods of PDO
![]()
图 5 暖期和冷期模态下柔鱼中心渔场重心与海洋表面温度的关系
Figure 5. Relationship between central gravity point of fishing ground and SST in warm and cold periods
表 1 西北太平洋柔鱼冬春生群资源样本统计参数与全局空间自相关分析
Table 1 Statistical parameters and global spatial autocorrelation analysis of O. bartramii resources (winter & spring groups) in Northwest Pacific
年份
Year最大值
Max.均值
Mean标准差S
Std. deviation偏态
Skewness峰态
Kurtosis变异值
Cv=S/m集聚系数
S2/mZ得分
Z-scoreP 全局相关系数*
Gearyc*2004 10.67 2.64 1.81 0.78 0.76 0.69 1.24 10.07 0.00 0.13 2005 8.00 2.60 1.68 0.78 0.69 0.53 1.09 7.43 0.00 0.31 2006 5.67 1.88 1.43 0.33 0.33 0.66 1.09 8.47 0.00 0.36 2007 15.00 4.28 2.71 0.85 1.31 0.63 1.71 4.96 0.00 0.22 2008 9.23 2.39 1.85 0.78 0.29 0.77 1.43 9.12 0.00 0.10 2009 3.33 1.02 1.81 0.59 0.19 0.79 3.19 1.94 0.00 0.34 2010 5.00 1.60 1.33 0.53 0.18 0.65 1.11 11.83 0.00 0.34 2011 16.50 1.90 1.56 3.58 22.96 0.82 1.28 1.94 0.00 0.12 2012 8.50 2.24 1.90 1.44 1.56 0.85 1.61 3.60 0.00 0.53 2013 13.74 1.80 1.41 3.29 18.98 0.73 1.11 1.89 0.00 0.13 2014 5.47 1.31 1.87 1.20 2.68 0.66 2.66 6.03 0.00 0.22 2015 7.00 2.60 1.82 0.78 0.69 0.53 1.27 3.87 0.00 0.28 平均值 Average 9.01 2.19 1.76 1.24 4.13 0.69 1.57 5.93 0.00 0.26 注：*. 当S2/m>1 表示集聚分布，S2/m=1表示随机分布，S2/m<1表示均匀分布；*. 大于1表示负相关，等于1表示不相关，而小于1表示正相关 Note: *. S2/m>1 indicates the agglomeration distribution; S2/m=1 indicates a random distribution; S2/m<1 indicates a uniform distribution; *. Greater than 1 indicates a negative correlation, and that equal to 1 indicates no correlation, but that less than 1 indicates a positive correlation. 
下载: 导出CSV
-
[1] 余为, 陈新军. 西北太平洋柔鱼冬春生群体栖息地的变化研究[J]. 海洋学报, 2018, 40(3): 86-94. [2] 魏广恩. 北太平洋柔鱼渔场的时空分析与资源丰度的预测[D]. 上海: 上海海洋大学, 2018: 6-8. [3] YU W, CHEN X J, YI Q, et al. A review of interaction between neon flying squid (Ommastrephes bartramii) and oceanographic variability in the north Pacific Ocean[J]. J Ocean Univ Chin, 2015, 14(4): 739-748. doi: 10.1007/s11802-015-2562-8
[4] BOWER J R, ICHII T. The red flying squid (Ommastrephes bartramii): a review of recent research and the fishery in Japan[J]. Fish Res, 2005, 76(1): 39-55. doi: 10.1016/j.fishres.2005.05.009
[5] CHEN X J, CHEN Y, TIAN S Q. An assessment of the west winter–spring cohort of neon flying squid (Ommastrephes bartramii) in the northwest Pacific Ocean[J]. Fish Res, 2008, 92(2/3): 221-230. doi: 10.1016/j.fishres.2008.01.011
[6] 唐峰华, 樊伟, 伍玉梅, 等. 北太平洋柔鱼渔场资源与海洋环境关系的季节性变化[J]. 农业资源与环境学报, 2015, 32(3): 242-249. [7] 崔雪森, 唐峰华, 张衡, 等. 基于朴素贝叶斯的西北太平洋柔鱼渔场预报模型的建立[J]. 中国海洋大学学报(自然科学版), 2015, 45(2): 37-43. [8] 余为, 陈新军, 易倩. 不同气候模态下西北太平洋柔鱼渔场环境特征分析[J]. 水产学报, 2017, 41(4): 525-534. [9] 余为. 西北太平洋柔鱼冬春生群对气候与环境变化的响应机制研究[D]. 上海: 上海海洋大学, 2016: 70-88. [10] 刘冀, 董晓华, 李帅. 基于小波消噪的径流趋势变化对比分析[J]. 水电能源科学, 2009, 27(5): 4-7. doi: 10.3969/j.issn.1000-7709.2009.05.002 [11] 杨勇, 普宗朝, 张山清, 等. 近48年新疆博州地区气候变化趋势的分析[J]. 石河子大学学报(自然科学版), 2010, 28(4): 446-452. [12] WANG X L, YAN F J, WANG T R, et al. Changing features of sunshine duration and effects on greenhouse production in Shouguang City [J]. Meteoral Environl Ves, 2016, 7(6): 76-79.
[13] 邱海军, 曹明明, 曾彬. 基于小波分析的西安降水时间序列的变化特征[J]. 中国农业气象, 2011, 32(1): 23-27. doi: 10.3969/j.issn.1000-6362.2011.01.005 [14] 杨晓明, 周应祺, 陈新军, 等. 西北印度洋延绳钓大眼金枪鱼渔获率的时间序列小波特征[J]. 中国水产科学, 2009, 16(6): 941-948. doi: 10.3321/j.issn:1005-8737.2009.06.015 [15] 蔡建堤, 陈方平, 吴建绍. 渔场重心信度测算及渔场重心修正理论构建——以闽南-台湾浅滩渔场二长棘鲷为例[J]. 生态学报, 2015, 35(6): 1929-1937. [16] 魏广恩, 陈新军, 李纲. 西北太平洋柔鱼洄游重心年际变化及预测[J]. 上海海洋大学学报, 2018, 27(4): 573-583. doi: 10.12024/jsou.20171102171 [17] 冯永玖, 陈新军, 杨铭霞, 等. 基于ESDA的西北太平洋柔鱼资源空间热点区域及其变动研究[J]. 生态学报, 2014, 34(7): 1841-1850. [18] 崔丽. 空间多尺度下西北太平洋柔鱼空间模式及其与海洋环境的关系[D]. 上海: 上海海洋大学, 2018: 17-20. [19] CHEN X J, ZHAO X H, CHEN Y. Influence of El Niño/La Niña on the western winter-spring cohort of neon flying squid (Ommastrephes bartramii) in the northwestern Pacific Ocean[J]. ICES J Mar Sci, 2007, 64(6): 1152-1160. doi: 10.1093/icesjms/fsm103
[20] 杨铭霞, 陈新军, 冯永玖, 等. 中小尺度下西北太平洋柔鱼资源丰度的空间变异[J]. 生态学报, 2013, 33(20): 6427-6435. [21] GENCAY R, SELCUK F, WITCHER B. An introduction to wavelets and other filtering methods in finance and economics[J]. Wave Random Media, 2002, 12(3): 399-399. doi: 10.1088/0959-7174/12/3/701
[22] HAN M, NAICHANG Y, CHI H. Signal analysis by time-frequency distributions: applications to impulse radar scattering echoes of plate and dihedral corner[J]. Int J Infrared Millim Wave, 1998, 19(7): 1017-1025. doi: 10.1023/A:1022635710654
[23] 牟萍, 杨胜发, 童思陈. 基于Morlet小波分析的嘉陵江流域年径流周期研究[J]. 重庆交通大学学报 (自然科学版), 2016, 35(1): 76-79. [24] 田俊, 黄淑娥, 祝必琴, 等. 江西双季早稻气候适宜度小波分析[J]. 江西农业大学学报, 2012, 34(4): 646-651, 670. doi: 10.3969/j.issn.1000-2286.2012.04.005 [25] 杨晓明, 戴小杰, 田思泉, 等. 中西太平洋鲣鱼围网渔业资源的热点分析和空间异质性[J]. 生态学报, 2014, 34(13): 3771-3778. [26] DORMANN C. Effects of incorporating spatial autocorrelation into the analysis of species distribution data[J]. Glob Ecol Biogeogr, 2010, 16(2): 129-138.
[27] MARTINO F D, SESSA S. A fuzzy particle swarm optimization algorithm and its application to hotspot events in spatial analysis[J]. J Ambient Intell Hum Comput, 2013, 4(1): 85-97. doi: 10.1007/s12652-011-0096-5
[28] TSENG C T, SUN C L, YEH S Z, et al. Spatio-temporal distributions of tuna species and potential habitats in the Western and Central Pacific Ocean derived from multi-satellite data[J]. Int J Remote Sens, 2010, 31(17/18): 4543-4558.
[29] 田思泉, 陈新军, 冯波, 等. 西北太平洋柔鱼资源丰度与栖息[J]. 上海海洋大学学报, 2009, 18(5): 586-592. [30] 崔雪森, 唐峰华, 周为峰, 等. 基于支持向量机的西北太平洋柔鱼渔场预报模型构建[J]. 南方水产科学, 2016, 12(5): 1-7. doi: 10.3969/j.issn.2095-0780.2016.05.001 [31] 马金, 陈新军, 刘必林, 等. 北太平洋柔鱼渔业生物学研究进展[J]. 上海海洋大学学报, 2011, 20(4): 563-570. [32] 武胜男, 余为, 陈新军, 等. 太平洋褶柔鱼秋生群产卵场环境变化及对资源丰度的影响[J]. 海洋渔业, 2018, 40(2): 4-13. [33] YU W, CHEN X J, YI Q, et al. Variability of suitable habitat of western winter-spring cohort for neon flying squid in the Northwest Pacific under anomalous environments[J]. PLoS One, 2015, 10(4): e0122997. doi: 10.1371/journal.pone.0122997
[34] 陈峰, 陈新军, 钱卫国, 等. 水温变动对2009年西北太平洋柔鱼产量下降的影响[J]. 广东海洋大学学报, 2010, 30(1): 65-71. doi: 10.3969/j.issn.1673-9159.2010.01.013 [35] 陈新军, 田思泉. 西北太平洋海域柔鱼的产量分布及作业渔场与表温的关系研究[J]. 中国海洋大学学报(自然科学版), 2005, 35(1): 101-107. -
期刊类型引用(3)
1. 马晨夕，赵金良，宋银都，赵亮亮. 不同营养水平对鳜幼鱼糖异生和糖原分解的影响. 水产科学. 2024(05): 707-716 . 
百度学术
2. 王崇懿，王伦，刘建勇，傅学丽. 凡纳滨对虾(Litopenaeus vannamei)不同品系生长与耐高盐性状配合力及杂种优势分析. 海洋与湖沼. 2022(01): 161-167 . 百度学术
3. 张瑞祺，赵金良，郝月月，宋银都. 鳜颅部侧线系统的胚后发育. 南方水产科学. 2020(06): 57-66 . 本站查看
其他类型引用(3)
计量
- 文章访问数: 5250
- HTML全文浏览量: 2827
- PDF下载量: 50
- 被引次数: 6
粤公网安备 44010502001741号
