Effect of dietary protein level on hepatic function and antioxidant capacity of juvenile common carp (Cyprinus carpio)
-
摘要:
研究了不同蛋白水平饲料对鲤(Cyprinus carpio)幼鱼肝功能和抗氧化能力的影响。选用体质量为(10±1.2)g的健康鲤幼鱼, 随机分为5组, 每组设3个重复, 每个重复35尾鱼。放养在15个规格相同的白色水箱中, 分别投喂蛋白水平为30%、30.5%、31%、31.5%和32%的饲料(分别记作Diet 1、Diet 2、Diet 3、Diet 4和Diet 5), 饲养实验时间为8周, 每隔4周对相关生理生化指标进行测定。结果表明, 养殖第4和第8周时, 31%蛋白水平组血清和肝胰脏的谷丙转氨酶(ALT)和谷草转氨酶(AST)活性均显著降低(P < 0.05);血清、肝胰脏和肾脏中超氧化物歧化酶(SOD)、过氧化氢酶(CAT)和溶菌酶(LZM)活性升高(P < 0.05), 丙二醛(MDA)浓度显著降低(P < 0.05)。根据实验结果, 从饲料配制的成本考虑, 31%蛋白水平饲料能够达到维持鲤幼鱼机体正常的抗氧化能力以及保护肝功能的营养需求。
Abstract:Common carps (Cyprinus carpio) with average body weight of (10±1.2) g were randomly divided into five treatments with three replicates (35 fish in each treatment), to investigate the effect of dietary protein level on hepatic function and antioxidant capacity by feeding diets at five protein levels (30%, 30.5%, 31%, 31.5% and 32%).Every four weeks, the related physiological and biochemical indices were determined.At the 4th and 8th week, the results show:1) Activities of alanine transaminase (ALT) and aspartate transaminase (AST) in serum and hepatopancreas in 31% protein group were significantly higher (P < 0.05).2) Compared with the fish fed with the other treatment groups, those fed with 31% protein group had significantly higher activities of superoxide dismutase (SOD), catalase (CAT) and lysozyme (LZM) (P < 0.05).3) Contents of malonaldehyde (MDA) of fish which were fed with 31% protein group decreased significantly (P < 0.05).Considering the cost of feed preparation for common carps, dietary protein of 31% meets the requirement for hepatic function and antioxidant capacity for the juveniles.
-
Keywords:
- common carp /
- dietary protein level /
- hepatic function /
- antioxidant capacity
-
全球气候变化对生态系统产生深远影响，制约着物种的生态习性和地理分布，是影响海洋生物分布格局的重要因素之一[1]。随着气候变化引起的海水酸化和温度上升等进程的持续，大多数海洋生物的新陈代谢、个体发育和生活史特征等在一定程度上受到了影响，导致群落结构发生变化，同时分布格局也随气候变化产生生境迁移[2-4]。这些改变可能导致生物面临种群数量显著下降或加速灭绝的风险，进而引起生态系统结构和功能紊乱、生物多样性丧失[5]。作为海洋生物的典型代表，长棘海星 (Acanthaster planci) 和褐拟鳞鲀 (Balistoides viridescens) 对气候变化的响应受到了广泛关注。了解它们对全球气候变化的响应并确定影响其生境分布的关键因素，对维持生物多样性的稳定以及珊瑚礁生态系统保护与恢复管理措施的制定具有至关重要的作用。
长棘海星隶属于棘皮动物门、海星纲、有棘目、长棘海星科、长棘海星属，分布于印度-太平洋热带珊瑚礁海域[6]，是著名的“珊瑚杀手”。据统计，1963年以来澳大利亚大堡礁共经历了3次长棘海星大规模暴发，导致大量珊瑚白化，珊瑚礁生态系统严重退化、生物多样性显著降低[7]。1985—2012年间对珊瑚礁的广泛调查发现，全球珊瑚礁的活珊瑚覆盖率从28.0%下降至13.8%，其中由长棘海星造成的损害高达42%，对珊瑚礁生态系统造成了严重危害[8]。由于人类对其天敌 [如法螺 (Charonia tritonis) 等] 的过度捕捞、营养盐大量流入等影响，导致长棘海星易出现周期性暴发[9]。了解长棘海星在未来不同气候情景下的潜在分布变化可为保护和修复珊瑚礁生态系统提供理论依据和科学指导。遗憾的是我国对于长棘海星的研究仍处于起步阶段，相关研究鲜有报道。褐拟鳞鲀是长棘海星的重要捕食者，隶属于辐鳍鱼纲、鲀形目、鳞鲀科，分布在热带和温带的淡水和海洋环境中[10]，具有非常高的经济价值。在印度洋-太平洋的海域中，褐拟鳞鲀会捕食成年长棘海星，与长棘海星的种群动态存在一定关联[11]。此外，气候变化对两个物种潜在适宜生境分布的影响也尚未得到评估。
物种分布模型 (Species distribution models, SDMs) 是研究气候变化对物种潜在生境分布格局影响及生物多样性响应的有效工具[12]。随着统计学和生态学理论的不断完善及计算机应用技术的深入发展，SDMs技术也得到了改进，并广泛应用于陆地、淡水和海洋物种的地理分布预测、生境气候适宜性评价及野生资源保护与管理等方面[13]。Feitosa等[14]利用SDMs并结合公开数据对小尾真鲨 (Carcharhinus porosus) 当前和未来的潜在分布做出预测，确定南美北部海岸为保护该物种最有潜力的区域；Zhang等[15]利用SDMs预测了未来不同情景下中国、韩国和日本沿海的日本少鳞鱚 (Sillago japonica) 的潜在生境分布情况，指出未来管理战略需考虑其适宜生境范围的缩小及北移，为制定管理策略提供了重要指导。可见，SDMs在模拟未来气候变化情景下物种的生境变化具有广泛应用。
目前关于长棘海星的研究主要集中在其生物累积与分布特征、暴发原因探讨、危害与防治等方面[16-17]，也有少数利用单一物种分布模型对长棘海星的研究[18]。而褐拟鳞鲀的相关报道则集中在其生物学特征与遗传方面[10,19]。本研究结合长棘海星、褐拟鳞鲀的现有分布数据，利用组合物种分布模型对两者的适宜生境分布进行模拟，揭示影响分布的关键环境因子，并绘制两个物种当前和未来不同气候情景下的适宜生境，为维持珊瑚礁生态可持续发展及制定未来气候适应管理决策提供科学参考。
1. 材料与方法
1.1 物种分布数据
本研究主要针对长棘海星及其捕食者褐拟鳞鲀在全球范围内的分布情况。通过查询全球生物多样性信息网络 (Global biodiversity information facility, GBIF, http://www.Gbif.org) 获取2000年以来长棘海星及褐拟鳞鲀的发生记录。为减少采样带来的偏差并与环境数据分辨率相匹配，筛选掉重复及坐标不精确的样本，保证每个栅格数据内随机选取1个样本，最终获得1 094个长棘海星和1 793个褐拟鳞鲀的数据用于模型构建 (图1)。建模过程中部分算法有时需要环境条件来与发现物种的环境进行对比，需要物种的实际和缺席分布数据[20]。由于缺乏两物种不存在点数据，因此运行模型时在全球地理范围内随机生成10 000个坐标点作为伪不存在点，以提高模型性能。
1.2 环境变量
综合考虑数据可用性和环境相关性，本研究以全球海洋环境数据为分析依据。当前和未来 (2040—2050年) 的海洋遥感数据来源于BIO-Oracle v2.1 (https://bio-oracle.org/) 数据库及Global Marine Environment Datasets数据库 (http://gmed.auckland.ac.nz/index.html)。共下载57个环境变量，经过处理后空间分辨率为0.083 333 333。考虑到环境因子间的自相关性会影响模型的准确度，因此计算每对变量间的Pearson相关系数，若 >|0.7| 则认为相关性较强，保留其中1个变量[21]，最终得到表1的12个环境变量。未来海洋环境数据的科学假设来自于4种排放情景，即“典型浓度目标 (Representative concentration pathways, RCP)”，分别代表不同CO2排放标准下的未来情景，包括RCP26 (最低温室气体排放)、RCP45 (中等温室气体排放)、RCP60 (较高温室气体排放) 和RCP85 (最高温室气体排放)[22]。本研究选择RCP26、RCP60和RCP85情景来模拟未来环境，用对应变量建模并对未来气候进行投影。
表 1 环境变量说明Table 1. Description of environmental variables变量 Variable 单位 Unit 最大叶绿素质量浓度 Maximum chlorophyll mg·m−3 最大流速 Maximum current velocity m·s−1 最小流速 Minimum current velocity m·s−1 平均溶解氧浓度Mean dissolved oxygen mol·m−3 最小铁离子浓度 Minimum iron μmol·m−3 最小硝酸盐浓度 Minimum nitrate mol·m−3 最小浮游动物浓度 Minimum phytoplankton μmol·m−3 最小初级生产力浓度
Minimum primary productivityg·m−3·d−1 平均盐度 Mean salinity ‰ 最小温度 Minimum temperature ℃ 深度 Depth m 离岸距离 Land distance km 1.3 物种分布模型构建
本研究基于R语言biomod2程序包对两个物种当前和未来的潜在分布区进行预测。选取广义线性模型 (Generalized Linear Model, GLM)、推进式回归树模型 (Generalized Boosting Model, GBM)、分类树分析 (Classification Tree Analysis, CTA)、人工神经网络 (Artificial Neural Network, ANN)、表面分布区分室 (Surface Range Envelope, SRE) 模型、柔性判别分析 (Flexible Discriminant Analysis, FDA) 以及随机森林 (Random Forest, RF) 7种单一模型。为构建更多单一模型结果，减小单次模型运行带来的误差，每种模型算法运行10次，模型参数设置为默认。而后对单一模型结果进行集成产生组合模型用于后续的结果分析，以提高模型精度。
1.4 模型评估方法
模型运行时，将数据随机拆分，80%作为训练集样本，20%作为测试集样本，对模型精度进行交叉验证。模型评价使用真实技巧统计 (True skill statistic, TSS) 和受试者操作特征曲线 (Receiver operating curve, ROC) 2个统计学参数[23]。ROC曲线以假正率为横坐标，真正率为纵坐标。TSS值综合敏感度和特异性指标，考虑了疏忽性和替代性的错误，其取值范围为 [0,1]，>0.8时拟合效果极好，0.6~0.8拟合效果好，0.4~0.6拟合效果一般[24]。
1.5 模型集成方法
综合考虑所有模型的运算结果，保留TSS和ROC值均 >0.9的单模型构成最终组合模型的结果。采用投票平均 (Committee Averaging, CA) 和概率加权平均 (Weighted Mean of Probabilities, WM) 对所有单一模型进行集成产生组合模型 (Ensemble Model)，并对模型的不确定性进行评估[25]。
2. 结果
2.1 模型性能评价
模型交叉验证结果显示，两个物种的模型性能结果存在一定的相似性 (图2)，表现最好的模型均为随机森林模型 (RF)，长棘海星RF的TSS和ROC值分别为0.96和0.99 (图2-a)，褐拟鳞鲀RF的TSS和ROC值分别为0.97和0.99 (图2-b)，表明模型的预测精度高，整体性能较优。然而SRE模型表现出了最差的性能 (图2)，两物种SRE模型的TSS值均为0.62、ROC值均为0.81，结果显示了不同单一模型间的不确定性。长棘海星与褐拟鳞鲀投票平均和概率加权平均组合模型的TSS、ROC值分别为0.96、0.99与0.97、0.99，评估结果表明，组合模型能够一定程度上减少不同单一模型的不确定性，提高预测的准确性。
![]()
图 2 基于TSS和ROC的7种物种分布模型统计结果注：a. 长棘海星；b. 褐拟鳞鲀；GLM. 广义线性模型；GBM. 推进式回归树模型；CTA. 分类树分析模型；ANN. 人工神经网络模型；SRE. 表面分布区分室模型；FDA. 柔性判别分析模型；RF. 随机森林模型。Figure 2. Results of statistical analyses of seven species distribution models based on analyses of true skill statistics and receiver operating curveNote: a. A. planci; b. B. viridescens; GLM. Generalized Linear Model; GBM. Generalized Boosting Model; CTA. Classification Tree Analysis; ANN. Artificial Neural Network; SRE. Surface Range Envelope; FDA. Flexible Discriminant Analysis; RF. Random Forest.2.2 环境变量重要性评估
筛选后参与建模的环境变量共12个 (表1)，其Pearson相关系数均 <|0.7|，表明变量间的相关性较弱，可提高预测准确度。对于长棘海星，最小温度和离岸距离的重要性最高，两者重要性之和超过70% (表2)，表明这两个变量是长棘海星空间分布的主要影响因子。对于褐拟鳞鲀，最小温度、平均溶解氧浓度和离岸距离的重要性最高，三者重要性之和占70%以上，因此，这3个变量可能是决定褐拟鳞鲀空间分布的主要影响因子 (表3)。
表 2 长棘海星影响因子重要性排序结果Table 2. Importance ranking of factors influencing distribution of A. planci影响因子
Environmental factor重要性
Importance/%排序
Sorting最小温度 Minimum temperature 43.34 1 离岸距离 Land distance 28.54 2 最小铁离子浓度 Minimum iron 8.03 3 平均溶解氧浓度 Mean dissolved oxygen 7.61 4 平均盐度 Mean salinity 2.75 5 最小初级生产力浓度
Minimum primary productivity1.90 6 深度 Depth 1.90 6 最小硝酸盐浓度 Minimum nitrate 1.48 8 最大流速 Maximum current velocity 1.27 9 最小浮游动物浓度
Minimum phytoplankton1.27 9 最大叶绿素浓度 Maximum chlorophyll 1.06 11 最小流速 Minimum current velocity 0.85 12 表 3 褐拟鳞鲀影响因子重要性排序结果Table 3. Importance ranking of factors influencing distribution of B. viridescens影响因子
Environmental factor重要性
Importance/%排序
Sorting最小温度 Minimum temperature 48.82 1 平均溶解氧浓度 Mean dissolved oxygen 15.94 2 离岸距离 Land distance 11.81 3 最小铁离子浓度 Minimum iron 8.86 4 平均盐度 Mean salinity 4.13 5 最小初级生产力
Minimum primary productivity2.36 6 最大流速 Maximum current velocity 2.17 7 最小硝酸盐浓度 Minimum nitrate 1.57 8 深度 Depth 1.57 8 最小浮游动物浓度
Minimum phytoplankton1.18 10 最大叶绿素浓度 Maximum chlorophyll 0.98 11 最小流速 Minimum current velocity 0.59 12 环境变量的响应曲线表示的是物种存在概率与对应环境变量之间的关系，可从中解析长棘海星和褐拟鳞鲀在不同环境下的生态位，从而确定其适宜生境。长棘海星的存在概率与温度整体上呈非线性关系，多数模型预测其存在概率随温度增加而逐渐增大并趋于稳定 (图3)。离岸距离响应曲线的结果显示出了模型间的差异，整体而言，离岸距离越小，长棘海星的生境适宜性越高。不同模型的铁离子浓度响应曲线也表现出一定的差异，长棘海星生境适宜性随铁离子浓度的增加而降低 (图3)。
![]()
图 3 预测长棘海星与褐拟鳞鲀发生概率随环境变量的响应曲线注：GLM. 广义线性模型；GBM. 推进式回归树模型；CTA. 分类树分析模型；ANN. 人工神经网络模型；SRE. 表面分布区分室模型；FDA. 柔性判别分析模型；RF. 随机森林模型。Figure 3. Response curves of predicted occurrence probability of A. planci and B. viridescens against environment variablesNote: GLM. Generalized Linear Model; GBM. Generalized Boosting Model; CTA. Classification Tree Analysis; ANN. Artificial Neural Network; SRE. Surface Range Envelope; FDA. Flexible Discriminant Analysis; RF. Random Forest.由图3可知，当温度低于6 ℃时，褐拟鳞鲀存在的概率很低，随着温度的增加，其存在的概率增大，当温度高于16 ℃时，其生境适宜性陡然增加，达0.8以上。而多数模型溶解氧浓度响应曲线的结果相似，其溶解氧浓度的适宜范围约为100~200 mol·m−3。褐拟鳞鲀的离岸距离响应曲线结果与长棘海星相似，不同模型间的差异较大。总体而言，离岸距离越小，生境适宜性越高。
2.3 长棘海星、褐拟鳞鲀适宜生境分布预测
2.3.1 长棘海星、褐拟鳞鲀当前在世界范围内的适宜生境分布预测
当前长棘海星潜在分布区主要集中在北印度洋和西太平洋海域，其极适宜生境广泛分布于澳大利亚北部海岸、印度尼西亚海岸以及中国南海一直延伸至东海海域，美洲中部海域及红海海域也有少量分布 (图4-a—4-b)。从褐拟鳞鲀的预测结果来看，其分布的主要区域与长棘海星大致重叠。其极适宜生境同样分布在澳大利亚北部海岸、印度尼西亚海岸以及中国南海海域 (图4-c—4-d)。与长棘海星不同的是，褐拟鳞鲀在非洲南部海岸、美洲中部海岸的潜在分布明显低于长棘海星，而在红海海域的分布明显高于长棘海星。两物种的预测结果总体表现为极高适宜生境与高适宜生境分布在热带、亚热带等低纬度地区，而低、中适宜生境向高纬度地区延伸扩散。
![]()
图 4 当前环境下长棘海星与褐拟鳞鲀的潜在分布适生区注：a、b. 长棘海星；c、d. 褐拟鳞鲀；WM. 概率加权平均组合模型；CA. 投票平均组合模型。Figure 4. Potential suitable habitat of A. planci and B. viridescens in current environmentNote: a, b. A. planci; c, d. B. viridescens; WM. Weighted Mean of Probabilities Ensemble Model; CA. Committee Averaging Ensemble Model.长棘海星潜在适宜生境面积达3 027.80万km2，其中极适宜生境面积为460.90万km2 (约占全球海洋面积的1.2%)，占适宜生境面积的15.22% (图4-a，表4)，非适宜生境面积占绝大部分，达37 039.68万km2 (92.4%)。褐拟鳞鲀的潜在适宜生境比长棘海星更广，其面积达3 813.84万km2，其中极适宜生境面积为590.00万km2 (约占全球海洋面积的1.5%)，占适宜生境面积的15.47% (图4-c，表4)，非适宜生境面积达36 253.64万km2 (90.5%)。
表 4 长棘海星与褐拟鳞鲀的适宜生境占比Table 4. Proportions of suitable habitat for A. planci and B. viridescens物种
Species气候情景
Climate scenario低适宜生境占比
Proportion of
low suitability
habitat/%中适宜生境占比
Proportion of
moderately suitable
habitat/%高适宜生境占比
Proportion of
highly suitable
habitat/%极适宜生境占比
Proportion of
extremely suitable
habitat/%长棘海星
A. planci现在Current 40.23 24.33 20.22 15.22 RCP26 40.35 24.11 20.12 15.42 RCP60 40.18 24.06 20.09 15.67 RCP85 40.13 24.05 20.14 15.68 褐拟鳞鲀
B. viridescens现在Current 42.22 24.07 18.24 15.47 RCP26 41.82 24.10 18.31 15.77 RCP60 41.86 23.99 18.18 15.97 RCP85 41.86 23.90 18.24 16.00 注：非适宜生境不在计算范畴内。 Note: Unsuitable habitat is not included in the calculation. 2.3.2 未来不同气候情景下长棘海星和褐拟鳞鲀的适宜生境预测
两种不同组合模型 (WM、CA) 预测的长棘海星的适宜生境面积存在差异，但预测的总体趋势均为增大，不同气候情景下增幅不同。WM模型预测未来长棘海星在RCP26、RCP60和RCP85情景下总适宜生境面积分别增加78.87万、74.70万和102.41万km2，增幅分别为2.60%、2.47%和3.38%，其中极适宜生境增幅分别为3.97%、5.49%和6.49% (图5-a、5-c、5-e)。而CA模型预测RCP26、RCP60和RCP85情景下总适宜生境面积分别增加84.14万、75.83万和107.71万km2，增幅为3.50%、3.15%和4.48%，其中极适宜生境分别增加3.41%、3.82%和4.60% (图5-b、5-d、5-f)。不同气候情景下长棘海星未来的潜在分布面积均会扩大，且CA模型下的增幅略高于WM模型，但WM模型预测的极适宜生境面积增幅明显高于CA模型。
![]()
图 5 长棘海星未来不同气候情景下的潜在分布注：WM. 概率加权平均组合模型；CA. 投票平均组合模型；RCP26. 最低温室气体排放情景；RCP60. 较高温室气体排放情景；RCP85. 最高温室气体排放情景。Figure 5. Potential future distribution of A. planci under different climate scenariosNote: WM. Weighted Mean of Probabilities Ensemble Model; CA. Committee Averaging Ensemble Model; RCP26. Minimum greenhouse gas emission scenarios; RCP60. Higher greenhouse gas emission scenarios; RCP85. Maximum greenhouse gas emission scenarios.褐拟鳞鲀全球潜在适宜生境在未来不同气候情景下也均表现为增加，WM模型预测RCP26、RCP60和RCP85情景下其总适宜生境面积分别增加90.75万、90.64万和112.09万km2，增幅分别为2.38%、2.37%和2.94%，其中极适宜生境分别增加4.38%、5.69%和6.49% (图6-a、6-c、6-e)。CA模型预测RCP26、RCP60和RCP85情景下总适宜生境面积分别增加111.13万、106.98万和145.99万km2，增幅为3.71%、3.58%和4.88%，其中极适宜生境分别增加5.02%、5.77%和6.78% (图6-b、6-d、6-f)。与长棘海星结果不同的是，CA模型预测褐拟鳞鲀的未来潜在分布总面积与极适宜生境面积增幅均高于WM模型，表现出模型对于不同物种的差异。除此之外，两物种的潜在分布区域均有向高纬度地区移动的趋势。
![]()
图 6 褐拟鳞鲀未来不同气候情景下的潜在分布注：WM. 概率加权平均组合模型；CA. 投票平均组合模型；RCP26. 最低温室气体排放情景；RCP60. 较高温室气体排放情景；RCP85. 最高温室气体排放情景。Figure 6. Potential future distribution of B. viridescens under different climate scenariosNote: WM. Weighted Mean of Probabilities Ensemble Model; CA. Committee Averaging Ensemble Model; RCP26. Minimum greenhouse gas emission scenarios; RCP60. Higher greenhouse gas emission scenarios; RCP85. Maximum greenhouse gas emission scenarios.3. 讨论
3.1 物种分布组合模型在适宜生境研究中的应用
多数研究表明，使用物种分布模型对海洋生物的潜在分布进行模拟时，各类单一模型的运行过程及算法存在不确定性和差异性，会降低预测的精度[26]。采用本研究构建的组合模型评估当前及未来3种不同气候情景下长棘海星与褐拟鳞鲀的潜在适生区动态变化，可有效降低单一模型造成的偏差。真实技巧统计值与受试者工作特征曲线是评价模型常用的指标[27]。研究结果显示，长棘海星与褐拟鳞鲀组合模型的TSS和ROC值高于大多数单一模型的结果，表明组合模型对两个物种潜在适生区的模拟结果准确。Weinert等[28]使用Biomod2程序包的6种单模型构建组合模型研究北海底栖生物的分布变化情况，结果也表明组合模型可降低模型的不确定性。本研究的模型预测结果与已知的长棘海星和褐拟鳞鲀的分布范围大致符合 (主要分布在北印度洋、西太平洋海域)。但本研究使用的建模方式仍存在一定的局限性。首先，物种扩散能力与种群动态息息相关，长棘海星与褐拟鳞鲀的扩散受到洋流及自身游泳能力的限制，但本研究并未考虑扩散的影响[29]。其次，人为因素对海洋生物的影响加剧[30]，本研究只考虑了部分重要环境因素的影响。此外，受限于未来环境变量的获取，以及在实际情况下存在各种限制因素，本研究可能会高估长棘海星和褐拟鳞鲀的分布范围。生物间的相互作用也会对物种产生重要影响，褐拟鳞鲀作为长棘海星的重要捕食者，对长棘海星的种群动态也会产生干扰，并未纳入此次研究的范围。
3.2 影响长棘海星与褐拟鳞鲀适宜生境分布的主要因素
长棘海星与褐拟鳞鲀均生活在珊瑚礁生境，其种群动态与珊瑚礁生态系统密切相关。模型预测影响长棘海星空间分布的主要因素是温度和离岸距离，其中温度是最重要的因子。对褐拟鳞鲀分布有重要影响的环境变量为温度、溶解氧浓度与离岸距离，同样温度为最重要的因子。温度与海洋生物的新陈代谢有关，从而对生物的种群动态产生影响，模型结果显示，当温度高于15 ℃，长棘海星的生境适宜性增加。温度对长棘海星的产卵及幼虫成活率有重要影响，海水温度达27 ℃阈值时，长棘海星才开始产卵。其幼虫最佳发育温度在不同发育阶段有所差异，胚囊幼虫期为13~34 ℃，羽腕幼虫期为14.5~32 ℃，中晚期幼虫为26~30 ℃，当高于32 ℃时幼虫会死亡，低于25 ℃时则停止发育并退回早期发育阶段[31-32]。海水温度对中上层鱼类的水平和垂直移动有重要影响[33]。海温、叶绿素与上升流是影响鱼类空间分布的主要驱动因素，如鳀鱼 (Engraulis ringens) 等[34]。由此可见，温度可能是最主要和最普遍的鱼类分布驱动因子。离岸距离影响长棘海星分布，其原因可能与珊瑚礁形成有关，但未见相关的研究报道。长棘海星以珊瑚为食，在很大程度上珊瑚决定长棘海星的空间分布，而珊瑚对于生长环境的要求比较严格，主要分布在热带浅水海域，离海岸的距离较近。陆架边缘鱼类分布广泛、种类多，许多鱼类会在繁殖季节回到近岸海域产卵，因此离岸距离是影响鱼类栖息地的重要环境因子。一些鱼类的分布可能与特定生境的距离高度相关[35]。此外，溶解氧是海洋鱼类获得氧气的主要来源[36]，对褐拟鳞鲀的栖息地分布起着重要作用。由于温室气体大量排放，未来全球气候变暖极有可能使当前适宜生存的环境因素超出物种的耐受阈值，诱发长棘海星和褐拟鳞鲀迁移到符合其生境偏好的新水域或在当前适宜生境做出适应性改变或急剧减少，从而引起两个物种的生境范围发生变化[37]。
3.3 不同气候情景下长棘海星与褐拟鳞鲀潜在适生区的变化
气候变暖导致全球海洋生物的丰度呈现向极地方向增加且向赤道方向减少的趋势[38]。本研究结果表明，当前环境下长棘海星的适宜栖息地主要分布于印度-太平洋海域、中国南海到澳大利亚海岸线、红海和阿拉伯海，与Sill和Dawson[18]的研究结果相似。本研究对长棘海星和褐拟鳞鲀当前和未来的栖息地分布做出了较为可靠的预测，可以为制定长棘海星暴发管理的策略提供一些参考。例如可针对长棘海星未来生境适宜性较高的海域，定期开展长棘海星的跟踪监测，对可能出现长棘海星暴发的海域进行预警[39]。当前环境下褐拟鳞鲀极适宜生境分布在40°N—40°S，其分布区域大部分覆盖了长棘海星的分布范围，可能是由于褐拟鳞鲀会捕食长棘海星，长棘海星存在的地方可以为它们提供部分食物来源，适合其生存[40]。未来不同气候情景下，长棘海星和褐拟鳞鲀的共性是两者的总体适生区范围均呈现出一定程度向高纬度地区扩大的趋势，模型预测温度是影响两个物种空间分布的最重要因素，未来全球气候变暖致使水温升高，促进了长棘海星与褐拟鳞鲀的生长发育。长棘海星具有逃离天敌存在区域的生物特性[41]，因此可在其可能暴发的海域大量投放成本较低的长棘海星幼体和成体捕食者，从而减少长棘海星暴发的风险。在升温2 ℃的情况下，长棘海星幼虫的成活率可能会提高240%，在食物供应不变的情况下，发育速度会提高30%[42]，从而导致其分布范围发生变化。本研究结果也符合未来条件下物种生存范围的变化趋势。研究显示气候变化改变了127种海洋生物的分布范围，其中大部分物种的潜在分布向高纬度海域移动[43]。在非洲南部海岸、美洲中部海岸，长棘海星的适宜生境明显比褐拟鳞鲀更广泛，这可能是因为在不适合褐拟鳞鲀栖息的地区，由于捕食者压力的潜在释放，会出现更为丰富的长棘海星种群[44]。到2050s，2种生物的潜在适宜生境在热带及亚热带海域的部分缩减，推测可能因气候变化使得海洋温度上升，海洋生物可能会向更适宜生存的极地迁徙[45]。在本研究记录数据以外的部分非洲和美洲沿岸海域也发现了褐拟鳞鲀的适宜栖息地，可能是由于气候变暖导致这些海域的理化性质改变，从而变得适合褐拟鳞鲀的生长。此外，控制海洋环境改变以及减少人为活动的干扰对减缓长棘海星暴发周期也具有重要作用[46]，通过控制陆源污染以改善珊瑚礁区水质、设置禁渔或保护区以保护鱼类种群丰度和多样性，同时减少对长棘海星天敌的过度捕捞等都是维护珊瑚礁生态系统稳定和可持续发展的必要措施。
4. 结论
本研究基于组合物种分布模型，预测了不同气候情景下长棘海星与褐拟鳞鲀全球范围内的潜在分布格局，并分别探讨了2个物种分布的主要影响因子。结果表明，组合模型的结果优于单一模型，一定程度上能够提高预测的准确性。影响长棘海星空间分布的主要因素是温度和离岸距离，而影响褐拟鳞鲀空间分布的主要因素是温度、溶解氧浓度和离岸距离。长棘海星与褐拟鳞鲀在未来气候情景下的潜在分布均有扩大且向高纬度地区移动的趋势。本研究结果将为长棘海星生态管理的应用、维护生态系统多样性、制定未来气候适应管理决策提供科学依据。
-
表 1 实验饲料组成及营养水平(风干基础)
Table 1 Composition and nutrient levels of experimental diets (air-dry basis)
% 项目 item 饲料编号 diet No. Diet 1 Diet 2 Diet 3 Diet 4 Diet 5 原料 ingredient 鱼粉 fish meal 1 2 3 4 5 河虫 clam worm 4 4 4 4 4 花生粕 peanut meal 8 8 8 8 8 豆粕 soybean meal 32 32 32 32 32 棉粕 cotton meal 8 8 8 8 8 菜籽粕 rapeseed 6 6 6 6 6 麸皮 bran 18 18 18 18 18 次粉 wheat middling 22 21 20 19 18 预混料 premix 1 1 1 1 1 合计 total 100 100 100 100 100 营养水平 nutrient level 粗蛋白质 crude protein 30.0 30.5 31.0 31.5 32.0 粗脂肪 ether extract 2.5 2.6 2.7 2.8 2.9 总能/kJ·g-1 gross energy 11.90 11.88 11.86 11.84 11.82 
下载: 导出CSV
表 2 饲料蛋白水平对鲤幼鱼肝胰脏转氨酶活性的影响
Table 2 Effect of dietary protein level on activitiy of aminotransferase in hepatopancreas of juvenile common carps
U·g-1 饲料编号
diet No.谷丙转氨酶 ALT 谷草转氨酶 AST 第4周 第8周 第4周 第8周 Diet 1 3.16±1.30a 1.34±0.04a 20.34±5.20a 4.98±0.54 Diet 2 1.96±0.04a 0.77±0.08bc 12.38±1.32b 4.41±0.46 Diet 3 0.73±0.22b 0.62±0.04c 10.14±0.36b 4.03±0.30 Diet 4 2.78±0.78a 1.16±0.07a 26.30±0.34a 4.41±0.22 Diet 5 3.78±1.30a 0.94±0.08b 26.87±0.34a 4.64±0.36 注:表中数据为平均数及3个重复的标准差; 同列数据肩标不同小写字母表示差异显著(P < 0.05), 后表同此。
Note: The data are the average values and standard deviation of three replicates; the data in the same column with different superscript lowercase letters indicate significant difference (P < 0.05).The same case in the following tables.
下载: 导出CSV
表 3 饲料蛋白水平对鲤幼鱼血清转氨酶活性的影响
Table 3 Effect of dietary protein level on activity of aminotransferase in serum of juvenile common carps
U·g-1 饲料编号
diet No.谷丙转氨酶 ALT 谷草转氨酶 AST 第4周 第8周 第4周 第8周 Diet 1 0.44±0.12a 6.26±0.58a 1.63±0.06a 10.58±2.50a Diet 2 0.12±0.01b 2.63±0.56b 0.74±0.01b 7.22±1.08ab Diet 3 0.11±0.03b 1.08±0.13c 0.22±0.05c 4.92±0.72b Diet 4 0.12±0.01b 2.48±0.16bc 0.70±0.12b 5.34±1.24b Diet 5 0.29±0.08ab 6.53±0.70a 0.70±0.35b 8.78±0.80ab
下载: 导出CSV
表 4 饲料蛋白水平对鲤幼鱼超氧化物歧化酶活性的影响
Table 4 Effect of dietary protein level on activity of superoxide dismutase of juvenile common carps
U·g-1 饲料编号
diet No.血清 serum 肝胰脏 hepatopancreas 肾脏 kidney 第4周 第8周 第4周 第8周 第4周 第8周 Diet 1 18.54±0.32c 277.99±58.88b 9.70±0.33b 12.75±1.42c 8.54±0.31c 163.45±17.18c Diet 2 20.08±0.49bc 437.56±57.64ab 14.28±2.50a 17.44±0.34bc 10.08±0.48bc 238.40±13.94b Diet 3 24.34±0.65a 509.12±40.84a 14.80±0.26a 30.34±2.20a 14.34±0.65a 264.23±0.40a Diet 4 21.59±0.94b 424.90±48.86ab 11.74±0.69ab 20.58±2.59b 11.58±0.94b 162.94±8.10c Diet 5 15.16±0.57d 469.86±55.43a 13.39±0.74ab 22.78±3.88b 5.16±0.58d 159.88±6.84c
下载: 导出CSV
表 5 饲料蛋白水平对鲤幼鱼过氧化氢酶活性的影响
Table 5 Effect of dietary protein level on activity of catalase of juvenile common carps
U·g-1 饲料编号
diet No.血清 serum 肝胰脏 hepatopancreas 肾脏 kidney 第4周 第8周 第4周 第8周 第4周 第8周 Diet 1 9.94±0.88a 5.32±0.27ab 0.37±0.06b 9.10±1.84 1.27±0.26b 10.08±0.64b Diet 2 11.65±0.70a 6.26±0.42a 0.46±0.02b 9.68±0.84 1.74±0.24b 11.30±1.05b Diet 3 12.20±0.48a 6.56±0.33a 3.98±0.76a 15.66±5.56 4.52±0.45a 14.34±1.50a Diet 4 5.96±0.44b 4.00±1.15b 0.62±0.22b 12.36±1.24 4.61±0.30a 11.20±0.16b Diet 5 11.78±1.02a 1.07±0.04c 0.47±0.18b 9.02±2.24 1.20±0.40b 10.17±0.62b
下载: 导出CSV
表 6 饲料蛋白水平对鲤幼鱼丙二醛质量摩尔浓度的影响
Table 6 Effect of dietary protein level on concentration of malonaldehyde of juvenile common carps
nmol · mg-1 饲料编号
diet No.血清 serum 肝胰脏 hepatopancreas 肾脏 kidney 第4周 第8周 第4周 第8周 第4周 第8周 Diet 1 15.81±0.46a 19.48±0.83bc 1.74±0.02a 0.60±0.06a 19.84±2.18b 5.11±0.25ab Diet 2 11.76±1.62b 19.70±0.66b 1.14±0.18b 0.48±0.06ab 25.92±1.78a 3.75±0.18bc Diet 3 7.84±0.94c 15.62±0.30c 0.82±0.01b 0.43±0.02b 16.08±1.06b 3.00±0.09c Diet 4 10.94±1.30bc 17.70±0.52bc 1.10±0.08b 0.49±0.04ab 28.56±1.16a 4.53±0.45b Diet 5 9.59±0.75bc 24.44±2.53a 1.52±0.19a 0.46±0.02ab 18.38±1.18b 5.61±0.50a
下载: 导出CSV
表 7 饲料蛋白水平对鲤幼鱼溶菌酶活性的影响
Table 7 Effect of dietary protein level on activity of lysozyme of juvenile common carps
U·mg-1 饲料编号
diet No.血清 serum 肝胰脏 hepatopancreas 肾脏 kidney 第4周 第8周 第4周 第8周 第4周 第8周 Diet 1 8.16±0.44a 0.76±0.24ab 11.02±0.07a 4.72±0.28ab 28.20±3.62b 5.38±0.36b Diet 2 8.18±0.27a 0.97±0.12a 11.34±0.51a 5.58±0.27ab 33.89±0.54a 6.98±0.24ab Diet 3 8.46±0.20a 1.07±0.04a 11.76±0.99a 5.94±1.23a 36.08±0.30a 8.07±0.13a Diet 4 7.80±0.26a 0.72±0.18ab 10.07±0.14b 5.47±0.36ab 35.16±2.65a 7.56±0.13a Diet 5 7.48±0.28b 0.34±0.14b 9.02±0.22b 3.67±0.34b 27.00±2.69b 5.74±0.85b
下载: 导出CSV
-
[1] 李爱杰.水产动物营养与饲料学[M].北京:中国农业出版社, 1996:26-36. [2] 张加润, 林黑着, 黄忠, 等.饲料中用混合植物蛋白并添加氨基酸替代鱼粉对斑节对虾生长及免疫力的影响[J].南方水产科学, 2013, 9(5):44-50. http://www.schinafish.cn/CN/abstract/abstract9068.shtml [3] 艾春香, 陶青燕.鱼粉替代—鱼粉高价运行下水产配合饲料研发的技术对策[J].饲料工业, 2013, 34(10):1-7. http://www.wenkuxiazai.com/doc/6bb4ab7331b765ce050814aa-2.html [4] 叶元土, 蔡春芳, 张宝彤, 等.鱼类营养与饲料配制[M].北京:化学工业出版社, 2013:116. [5] ZEITLER M H, KIRCHGEBNER M, SCHWARZ F J.Effects of different protein and energy supplies on carcass composition of carp (Cyprinus carpio L.)[J].Aquaculture, 1983, 36:37-48. http://www.sciencedirect.com/science/article/pii/0044848684900528
[6] 李贵锋, 蒋广震, 刘文斌, 等.不同蛋白质和能量水平对建鲤幼鱼生长性能、体组成和消化酶活性的影响[J].上海海洋大学学报, 2012, 21(2):225-232. http://www.cnki.com.cn/Article/CJFDTOTAL-SSDB201202011.htm [7] 李爱杰, 徐玮, 孙鹤田, 等.鲤鱼营养需要的研究[J].水利渔业, 1999, 19(5):18-20. http://www.cnki.com.cn/Article/CJFDTOTAL-SCAN199905010.htm [8] FAN Z, LI J H, CHENG Z Y, et al.Protein sparing effect of lipid diets for common carp (Cyprinus carpio)[J].Adv Engin Res, 2015, 45:357-368. http://www.academia.edu/13056658/Differential_gene_expression_after_total...
[9] 郭立, 井润珍, 程镇燕, 等.降低鲤鱼饲料蛋白质的初步研究[J].饲料工业, 2013, 34(8):41-45. http://www.cnki.com.cn/Article/CJFDTOTAL-FEED201308010.htm [10] 范泽, 李静辉, 王安琪, 等.饲料中木薯淀粉添加水平对鲤生长性能、消化能力及糖代谢的影响[J].动物营养学报, 2016, 28(4):4404-4053. http://www.cnki.com.cn/Article/CJFDTOTAL-DWYX201612040.htm [11] 喻礼怀, 张军, 魏宗友, 等.饲粮色氨酸水平对扬州鹅蛋白质代谢的影响[J].动物营养学报, 2014, 26(1):81-89. http://www.chinajan.com/CN/abstract/abstract11217.shtml [12] 王桂芹, 李子平, 牛小天, 等.饲料能量和维生素B6对乌鳢生长和蛋白质代谢酶活性的影响[J].中山大学学报(自然科学版), 2011, 50(4):96-99, 105. http://lprapp14.fao.org/XML_Output/2010/CN/CN0912.xml [13] 范泽, 孙金辉, 程镇燕, 等.降低饲料中植物性蛋白水平对鲤幼鱼生长和肌肉营养成分及转氨酶活力的影响[J].中国畜牧杂志, 2016, 52(21):59-64. doi: 10.3969/j.issn.0258-7033.2016.21.014 [14] 孙瑞健, 张文兵, 徐玮, 等.饲料蛋白质水平与投喂频率对大黄鱼生长、体组成及蛋白质代谢的影响[J].水生生物学报, 2013, 37(2):281-289. doi: 10.7541/2013.15 [15] 林淑琴, 成琳, 张文兵, 等.饲料蛋白质水平对大黄鱼体内蛋白质的沉积和代谢反应的影响[J].中国海洋大学学报(自然科学版), 2014, 44(5):32-39. http://www.cnki.com.cn/Article/CJFDTOTAL-QDHY201405005.htm [16] BARCELLOS L J G, KREUTZ L C, SOUZA C, et al.Hematological changes in jundia (Rhamlda qnelen Quoy and Gaimard pimelodadae) after acute and chronic stress caused by usual aquacultural management, with emphasis on immunosuppressive effects [J].Aquaculture, 2004, 237(1/2/3/4):220-236. http://europepmc.org/abstract/AGR/IND43637873
[17] 强俊, 杨弘, 王辉, 等.饲料蛋白水平对低温应激下吉富罗非鱼血清生化指标和HSP70 mRNA表达的影响[J].水生生物学报, 2013, 37(3):434-443. doi: 10.7541/2013.40 [18] 李彬. 饲料蛋白质水平对三种规格草鱼生长、饲料利用以及氮代谢的影响[D]. 武汉: 华中农业大学, 2014: 15, 19. http://cdmd.cnki.com.cn/Article/CDMD-10504-1014231425.htm [19] KUCUKBAY F Z, YAZLAK H, KARACA L, et al.The effects of dietary organic or inorganic selenium in rainbow trout (Oncorhynchus mykiss) under crowding conditions [J].Aquacult Nutr, 2009, 15:569-576. doi: 10.1111/anu.2009.15.issue-6
[20] 杨弘, 徐起群, 乐贻荣, 等.饲料蛋白质水平对尼罗罗非鱼幼鱼生长性能、体组成、血液学指标和肝脏非特异性免疫指标的影响[J].动物营养学报, 2012, 24(12):2384-2392. doi: 10.3969/j.issn.1006-267x.2012.12.013 [21] 唐媛媛. 饲料蛋白质水平对卵形鲳鲹消化酶、非特异性免疫和肠道内容物细菌的影响[D]. 厦门: 厦门大学, 2014: 44-57. http://cdmd.cnki.com.cn/Article/CDMD-10384-1014227609.htm [22] 吴庆元, 蒋玫, 李磊, 等.低盐度胁迫对鲻鱼(Mugil cephalus)幼鱼鳃丝、肌肉、肠Na +-K +-ATP酶活性和MDA浓度的影响[J].生态与农村环境学报, 2014, 30(4):481-487. http://www.cnki.com.cn/Article/CJFDTOTAL-NCST201404012.htm [23] 周立斌, 张伟, 王安利, 等.饲料中锌对花鲈(Lateolabrax japonicus)幼鱼生长、免疫和组织积累的影响[J].海洋与湖沼, 2009, 40(1):42-47. doi: 10.11693/hyhz200901007007 [24] SAHIN K, YAZLAK H, ORHAN C, et al.The effect of lycopene on antioxidant status in rainbow trout (Oncorhynchus mykiss)reared under high stockingdensity [J].Aquaculture, 2014, 418/419:132-138. doi: 10.1016/j.aquaculture.2013.10.009
[25] HUANG J F, XU Q Y, CHANG Y M.Effects of temperature and dietary protein on gene expression of Hsp70 and Wap65 and of juvenile mirror carp (Cyprinus carpio)[J].Aquacult Res, 2015, 46(11):2776-2788. doi: 10.1111/are.2015.46.issue-11
[26] 肖温温. 饲料中脂肪与蛋白质水平对大口黑鲈生长、体组成、非特异性免疫和血液学的影响[D]. 上海: 上海海洋大学, 2012: 45-46. http://cdmd.cnki.com.cn/Article/CDMD-10264-1013000335.htm [27] 乐贻荣, 杨弘, 徐起群, 等.饲料蛋白水平对奥尼罗非鱼(Oreochromis niloticus×O.aureus)生长、免疫功能以及抗病力的影响[J].海洋与湖沼, 2013, 44(2):493-498. doi: 10.11693/hyhz201302036036
计量
- 文章访问数: 3116
- HTML全文浏览量: 1026
- PDF下载量: 610
粤公网安备 44010502001741号
