Effects on growth and hepatic glucose metabolism of grass carp fed with high dietary carbohydrates
-
摘要:
为了从生长性状、器官健康和营养代谢水平方面探究草食性鱼类草鱼 (Ctenopharyngodon idella) 对高糖饲料的适应和利用能力,以初始体质量 (132.01±16.43) g的草鱼为研究对象,分别配制正常糖水平 (15%,对照组) 和高糖水平 (40%,高糖组) 的等氮等能饲料分组投喂草鱼140 d,每组设3个重复,每个重复300 尾鱼。结果显示,相较于对照组,第60、第80 和第140天时高糖组草鱼的体质量显著增加 (P<0.05),第140天时,高糖组草鱼的脏体比、肠系膜脂肪系数、中肠绒毛高度和肝糖原含量均显著提高且肝脏和中肠肠道组织出现损伤。随实验时间的延长,高糖组葡萄糖激酶 (gk)、丙酮酸激酶 (pk)、糖原合酶 (gys2)、脂肪酸合成酶 (fas) 基因表达量呈下降趋势,磷酸烯醇式丙酮酸激酶 (pepck)、葡萄糖-6-磷酸酶 (g6pase) 基因表达量呈上升趋势,对照组gys2、fas、pepck和g6pase基因表达量变化趋势与高糖组一致,但gk、pk基因表达量在各时间点均无显著性变化。第140天时,与对照组相比,高糖组pk基因表达量显著下调,g6pase、pepck和fas基因表达量显著上调。综上,采用高糖饲料投喂草鱼可促进其体质量增加,但长期饲喂会对其生理指标造成一定负面影响,肝脏糖代谢相关基因表达变化显著。
Abstract:In order to explore the growth performance, organic health and nutritional metabolism level of Ctenopharyngodon idella, a fed herbivorous fish-grass carp, with high dietary carbohydrates, we divided 1 800 individuals of (132.01±16.43) g into normal carbohydrate diet group (Control group) and high carbohydrate diet group (Group H), feeding them with 15% and 45% level carbohydrate diets for 140 d. Each group included three replications and each replication included 300 individuals. The results show that compared with the control group, the body mass of Group H was higher significantly on 60th, 80th and 140th day (P<0.05). The visceral body ratio, mesenteric fat coefficient, villus height and hepatic glycogen content of Group H increased significantly on 140th day, and damage was observed in liver and intestinal tract. The expression of glucose kinase gene (gk), pyruvate kinase gene (pk), glycogen synthase 2 gene (gys2) and fatty acid synthetase gene (fas) showed a decreasing trend in Group H along with time, while phosphoenolpyruvate carboxykinase gene (pepck) and glucose-6-phosphatase gene (g6pase) showed an increasing trend. The expressions of gys2, fas, pepck and g6pase had the same trend in Group L as that in Group H, while no significant variation trend was observed in gk and pk gene at each time. Compared with the control group, the expression of pk was significantly lower, while those of g6pase, pepck and fas were significantly higher in Group H. In conclusion, 40% level dietary carbohydrate can improve the weight gain rate of grass carp, but high dietary carbohydrates have a negative impact on physiological indicators and significantly affect glucose metabolism at molecular level.
-
Keywords:
- Ctenopharyngodon idella /
- Glucose metabolism /
- Liver /
- Intestinal structure
-
增殖放流是目前有效恢复渔业资源的方法之一[1],掌握放流种类的生态学特征是实施有效增殖放流的基本要求,选择适宜放流海区是确保放流物种成活率的关键。其中,物种的生境适宜性分析是保障物种增殖放流获得最大效益的重要环节。该过程需分析目标物种与所处环境之间的关系,分析目标物种生境中关键因子的需求,构建目标物种栖息地选择模型,最后将其潜在适生地理分布区直观地表达在数字化地形图或遥感图上[2]。头足类具有市场价值高、生长周期短和生长速度快等优点[3],对头足类增殖放流时进行物种的生境适宜性分析有利于进一步实现渔业资源可持续发展和蓝色经济增长的目标。
最大熵 (Maximum Entropy, MaxEnt) 模型是基于物种已知的分布数据和相关环境因子对物种分布及其生态需求进行模拟,从而分析得出物种的潜在适生区的物种分布模型 (Species Distribution Models, SDMs)[4]。在头足类研究领域,目前物种分布模型的应用以广义加性模型(Generalized Additive Model, GAM) 为主[5-7]。但与GAM模型相比,MaxEnt模型的预测结果较为准确[8-9]。MaxEnt模型的预测精度在一定程度上受所选生态因子的影响,但通过筛选生态因子和参数优化后,模型可在研究样本量较小的情况下获得较高的预测精度[10-11]。该模型不仅能有效处理环境因子间的相互作用,还可定量揭示环境因子对物种分布的影响[12]。适用于不同的生境类型、空间尺度和分类类群[13],在海洋生物保护方面极具发展前景[14],常被用于分析物种栖息地需求、气候变化对物种分布的影响、物种入侵监测和自然保护区的选择和保护[15]。例如,MaxEnt模型被应用于预测南太平洋长鳍金枪鱼 (Thunnus alalunga) 的栖息地[16],模拟水温变化对太平洋柔鱼 (Ommastrephes bartramii) 潜在栖息地的影响[17],分析底层鱼类在阿拉斯加湾的生境适宜性等[18]。
曼氏无针乌贼 (Sepiella maindroni) 是我国重要的头足类生物之一,隶属于软体动物门、头足纲、二鳃亚纲、十腕目、乌贼科、无针乌贼属。其可食部分比重大 (约占总体质量的92%)、蛋白质含量高且肉质富含钙 (Ca)、磷 (P)、铁 (Fe) 等元素[19],具有较高的营养和经济价值。曼氏无针乌贼曾是东海的优势物种,在浙江省的年产量曾高达6万吨[20],但受过度捕捞、亲体过度开发以及产卵场破坏等影响,其产量自1981年后持续下降[21],直至20世纪90年代初基本绝迹[22]。开展大规模增殖放流后,其产量在浙江海域呈上升趋势,至2008年增殖效果初显[21,23],但目前其产量仍不能满足人们的需求[24]。舟山海域岛屿众多,适宜曼氏无针乌贼卵附着的藻类丰富[25],是其重要的栖息地和产卵地。该海域内曼氏无针乌贼种群胴长范围为29~141 mm,体质量范围为5.1~341.8 g[26],目前对该海域内曼氏无针乌贼生境适宜性的研究较少。本文以舟山海域曼氏无针乌贼为研究对象,收集其分布数据和舟山海域的海洋环境数据,结合ArcGIS 10.8软件和MaxEnt模型对曼氏无针乌贼进行生境适宜性评价,以期为曼氏无针乌贼的增殖放流选址提供参考,并为其资源养护提供科学支撑。
1. 材料与方法
1.1 物种分布数据
构建MaxEnt模型需要物种分布点的经纬度数据。本研究的物种分布数据来源于2017—2021年浙江近海的底拖网调查数据、2021—2022年自主航次的底拖网调查数据和全球生物多样性信息机构GBIF (Global Biodiversity Information Facility, https://www.gbif.org/)。2017—2021年调查范围为120°E—131°E、23°N—33°N,站位采用均匀网格设置,每隔经纬度15'设置1个采样站点。2021—2022年自主航次调查范围为122°30'E—123°20'E、29°36'N—30°30'N。调查船为“浙普渔68712”单船底拖网渔船,使用单拖网,网目尺寸为80 mm,网囊最小网目尺寸为20 mm,每网拖曳1 h,平均拖速3 kn。调查时记录捕获曼氏无针乌贼时的经纬度信息。GBIF是一个由世界各国政府资助的国际开放基础数据库,数据来源包括物种名录、书籍、相关论文、动植物考察和观察资料以及数字化的在线数据集,目前已有5 000多万条分布点的数据记录。为扩大样本量,利用该数据库及相关文献资料,在底拖网调查数据的基础上补充曼氏无针乌贼的分布点数据。为避免因群集效应造成的取样偏差,利用ENMTools (http://purl.oclc.org/enmtools) 剔除冗余数据[27]。最终共得到曼氏无针乌贼有效分布样点167个 (图1),通过Excel 2020软件转换为csv格式。采用ArcGIS 10.8 软件绘制地图,数据来源于全国地理信息资源目录服务系统 [https://www.webmap.cn/,1∶100万公众版基础地理信息数据 (2021)]。
1.2 海洋环境数据
海洋环境数据集来源于Bio-ORACLE (https://bio-oracle.org) 和全球海洋环境数据库 (Global Marine Environment Datasets, GMED)。研究表明,海表面温度 (Sea surface temperature, SST)、海表面盐度 (Sea surface salinity, SSS)、pH、光照均可影响曼氏无针乌贼的生长和活动[28-30]。在此基础上补充选取影响海洋生物分布的相关生态因子 (海洋流速、初级生产力、硝酸盐含量等) 共同构成模型预实验的生态因子集合。
相关性过高的生态因子会影响模型的运行结果[31],为提高预测结果的准确性,对预实验中的生态因子进行筛选。利用SPSS 20.0软件中的相关性分析法和MaxEnt模型中的刀切法 (Jackknife),以贡献率小且相关性绝对值 |r|>0.8为标准,剔除易导致模型过拟合的生态因子[32],最终获得与曼氏无针乌贼相关的关键生态因子 (表1)。
表 1 筛选出的关键生态因子Table 1. Selected key ecological factors生态因子
Ecological factor单位
Unit海流流速最大值 Maximum currents velocity m·s−1 海流流速最小值 Minimum currents velocity m·s−1 溶解氧质量浓度最小值 Minimum dissolved oxygen mg·L−1 底层光照强度最大值 Maximum light at bottom μmol·(m2·s)−1 硝酸盐质量浓度最大值 Maximum nitrate mg·L−1 磷酸盐质量浓度最大值 Maximum phytoplankton mg·L−1 初级生产力最小值 Minimum primary productivity g·(m3·day)−1 海表盐度最小值 Minimum sea surface salinity ‰ 海表温度最大值 Maximum sea surface temperature ℃ 海表温度最小值 Minimum sea surface temperature ℃ 酸碱度 pH 1.3 模型参数优化及运行
导入数据至MaxEnt模型,运用由Cobos等[33]发明的kuenm软件包进行MaxEnt模型的参数优化。其利用了R语言和MaxEnt模型的功能,帮助实现校准、创建和评估最终的模型。为最大化获得参数优化效果,设置调控倍频 (Regularization multiplier, RM) 区间为0.5~4.0,间距为0.5 (0.5、1.0、1.5、2.0、2.5、3.0、3.5、4.0),其他设置为系统默认值。优先选择AIC (Akaike information criterion) 值最小和遗漏率最低的调控系数。AIC值是衡量统计模型拟合优良性的一种标准,能权衡所估计模型拟合数据的优良性,遗漏率是用于所用的数据在特定阈值上的评估,可分析建模时的数据是否存在空间自相关性[34]。根据运行结果,设置调控系数为3.0。
导入处理后的数据至MaxEnt 3.4.4 (https://biodiversityinformatics.amnh.org/open_source/maxent/)运行,勾选MaxEnt模型中Jackknife工具,获得主要生态因子对曼氏无针乌贼在舟山海域分布影响的权重值。随机选取75%的分布点用于建立模型,其余25%的分布点用于模型验证,且为了保证模型稳定性,采用10次重复交叉验证[35]。模型其他参数设置为默认值,结果以Logistic格式的栅格数据输出。
1.4 划分物种适生区等级
参考IPCC报告有关评估可能性的划分方法[36],根据曼氏无针乌贼分布的概率值,将其适宜区分为5个等级,分别为P<0.2、0.2≤P<0.4、0.4≤P<0.6、0.6≤P<0.8和0.8≤P<1,依次为非适生区 (Non-adaptive areas)、低适生区 (Low-adaptive areas)、中适生区 (Middle-adaptive areas)、次高适生区 (Sub-superior suitable area) 和高适生区 (High-adaptive areas)。
2. 结果
2.1 模型结果精度及检验
本文使用受试者工作特征曲线 (Receiver operating characteristic curve, ROC) 验证最大熵模型的预测精度。ROC曲线下面积 (Area under curve, AUC) 是评估预测物种分布的模型准确性的主要工具[37]。AUC值越高、越接近1则表示模型的性能越好,当AUC值大于0.9时,表示结果准确度较高[38]。应用表1中的生态因子进行建模,AUC值为0.995 (图2)。证明选择表1中的生态因子建立MaxEnt模型进行预测的结果准确度较高,可用于曼氏无针乌贼在舟山海域的生境适宜性分析。
2.2 物种对生态因子的响应
表2反映了生态因子对MaxEnt模型预测的相对贡献率,其中SST最大值、SSS最小值和pH是影响曼氏无针乌贼在舟山海域分布的主要生态因子。三者对模型的相对贡献率分别为27.8%、25.5%和20.3%,均超过20%,并且其相对贡献率之和占所有生态因子相对贡献率之和的一半以上。初级生产力最小值、磷酸盐质量浓度最大值、SST最小值和硝酸盐质量浓度最大值对模型的相对贡献率介于1%~15%,对曼氏无针乌贼分布的影响小于SST最大值、SSS最小值和pH。溶解氧质量浓度最小值、底层光照强度最大值、最大海流流速和最小海流流速的相对贡献率较小,均小于1%。
表 2 生态因子对模型的相对贡献率Table 2. Percent contributions of ecological factors to MaxEnt model生态因子
Ecological factor相对贡献率
Relative
contribution
rate/%海表温度最大值 Maximum sea surface temperature 27.8 海表盐度最小值 Minimum sea surface salinity 25.5 酸碱度 pH 20.3 初级生产力最小值 Minimum primary productivity 12.2 磷酸盐质量浓度最大值 Maximum phytoplankton 8.7 海表温度最小值 Minimum sea surface temperature 2.5 硝酸盐质量浓度最大值 Maximum nitrate 2.2 溶解氧质量浓度最小值 Minimum dissolved oxygen 0.4 底层光照强度最大值 Maximum light at bottom 0.2 海流流速最大值 Maximum currents velocity 0.2 海流流速最小值 Minimum currents velocity 0.1 生态因子重要性结果显示 (图3),仅包含SST最大值、SSS最小值和pH等生态因子建模时,模型获得的增益值较大。排除相对贡献率较大的生态因子 (如:SST最大值、SSS最小值和pH) 后,分别利用除这类生态因子外的其他生态因子建模,此时模型所获得的增益值较小。表明模型预测过程中SST最大值、SSS最小值和pH包含有价值的信息,无法被其他生态因子替代。
生态因子响应曲线可反映生态因子与曼氏无针乌贼存在概率之间的关系 (图4)。因相对贡献率较高的生态因子对曼氏无针乌贼分布的影响较大,故重点分析曼氏无针乌贼对SST最大值、SSS最小值、pH等生态因子的响应。以0.5概率值作为分界点,结果显示,舟山海域内SST最大值介于25.5~28.5 ℃,SSS最小值介于20‰~30‰,pH介于7.8~8.3时,曼氏无针乌贼存在的概率较大。
在SST最大值介于25.5~26.6 ℃时,曼氏无针乌贼的存在概率随温度升高而增加,温度达到26.6 ℃后其存在概率随温度升高而减少。SSS最小值响应曲线呈先增加后减少的变化,SSS最小值达25‰时曼氏无针乌贼物种存在的概率最大。在pH的响应曲线中,pH为8.2时曼氏无针乌贼的存在概率达到最大值。pH介于7.1~8.2时,其存在概率随pH升高呈上升趋势;pH介于8.2~8.4时,存在概率随pH升高呈下降趋势。
2.3 物种的适生区分布及检验
利用ArcGIS 10.8软件对模型结果进行重分类,获得曼氏无针乌贼在舟山海域的适生区分布图 (图5)。经计算,曼氏无针乌贼在舟山海域的适生区面积为30 213.40 km²,占研究区域总面积的75.11%。适生区分布范围较广,呈南北方向分布,北至马鞍列岛海域,南至舟山本岛东侧海域。高适生区面积为4 600.78 km²,主要分布在嵊泗列岛、马鞍列岛和中街山列岛海域,少量分散在马鞍列岛东侧。该海域内曼氏无针乌贼次高适生区面积为8 205.52 km²,中适生区面积为9 723.31 km²,低适生区面积为7 683.78 km²,三者面积相近 (表3)。高适生区、次高适生区分布范围内包含了舟山市内的3个海洋特别保护区 (图6)。
表 3 适生区面积及其占海域面积百分比Table 3. Area of suitable habitat and its percentage of sea area区域
Area概率
Probability面积
Square/km2占比
Percentage/%高适生区 High-adaptive area 0.8≤P<1 4 600.78 11.44 次高适生区 Sub-superior suitable area 0.6≤P<0.8 8 205.52 20.40 中适生区 Middle-adaptive area 0.4≤P<0.6 9 723.31 24.17 低适生区 Low-adaptive area 0.2≤P<0.4 7 683.78 19.10 非适生区 Non-adaptive area P<0.2 10 007.89 24.88 研究区域 Research area 40 221.29 100.00 为检验模型结果,采取访问调查的方式走访舟山市各区县典型渔村,开展曼氏无针乌贼捕捞出现海域调研。范围包括定海区 (金塘镇、长白乡)、普陀区 (沈家门、朱家尖、六横岛、桃花镇、东极岛)、岱山县 (衢山镇、长涂镇) 和嵊泗县 (嵊山镇、枸杞乡、洋山镇)。据2021年舟山市统计年鉴,2020年定海区、普陀区、岱山县和嵊泗县渔业人口分别为12 575、73 596、56 961和44 007人。综合考虑人力物力等客观因素及调研结果准确度,2022年共访问调查了1 285人,其中定海区87人、普陀区505人、岱山县391人和嵊泗县302人。所调查渔民中70.3%为专业从业人员,19.2%为兼业从业人员,10.5%为临时从业人员。调查统计结果表明,较多渔民在中街山列岛海域、嵊泗列岛海域和浪岗山列岛海域捕获曼氏无针乌贼,渔民数量分别占调查总体的32.84%、24.51%和15.02% (图7),该结果与模型预测结果总体相符。
3. 讨论
3.1 影响曼氏无针乌贼分布的主要生态因子
海表温度和海表盐度与头足类分布联系紧密[39],因此在头足类分布与生态因子的关系研究中常被作为重要生态因子[40-41]。同时,溶解氧、叶绿素、深度、悬浮物和混合层深度也常被作为分析头足类分布的重要环境变量。目前针对曼氏无针乌贼分布与环境变量的关系,主要探讨了水温、盐度、叶绿素、水深和溶解氧与其栖息地选择之间的偏好关系。pH、初级生产力、营养盐含量、海流流速以及光照强度也与海洋生物分布关系密切。但针对这部分生态因子对曼氏无针乌贼分布影响的研究较少,因此本文在其他学者研究[42-44]的基础上进行补充研究。该过程涉及的生态因子较多,为避免主观选取生态因子导致模型预测结果产生误差,文章采取相关性分析法和刀切法筛选生态因子。模型运行结果表明,SST最大值、SSS最小值和pH对模型的贡献率较大,是影响曼氏无针乌贼的主要生态因子,该结果与目前已有研究的结果总体一致[5,42]。初级生产力、营养盐含量对曼氏无针乌贼分布具有一定影响,而溶解氧浓度、底层光照强度和海流流速对其分布的影响较小。
SST最大值、SSS最小值和pH生态因子的响应曲线表明,曼氏无针乌贼偏向在海表温度较高、盐度较大的环境活动,且能接受弱碱性的环境,舟山海域内SST最大值介于25.5~28.5 ℃、SSS最小值介于20‰~30‰、pH介于7.8~8.3时,曼氏无针乌贼存在的概率较大。该结论与实验室中观察不同温度、盐度、pH下曼氏无针乌贼活力的实验结果相符[45],同时进一步验证了MaxEnt模型结果的可靠性。在舟山海域中,适宜曼氏无针乌贼生长的SST最大值介于25.5~28.5 ℃,且温度达到 (26±0.5) ℃时其存在概率可达0.75以上。符合在最佳温度范围内曼氏无针乌贼具有高生长速度和高成活率的特点[46]。Peng等[47]发现水温在25~27 ℃时,曼氏无针乌贼的瞬时生长速率、摄食速率、食物转化率和体长增加长度最佳。排除Peng等[47]在实验过程中所选种群 (浙南—闽东种群) 与舟山海域内曼氏无针乌贼 (浙北种群) 之间的差异,模型对曼氏无针乌贼在海表温度阈值的分析较为准确。SST最小值对模型的贡献率小于SST最大值 (表2),表明曼氏无针乌贼分布受高温变化的影响较大,且在低温条件下曼氏无针乌贼的活动和摄食频率降低,即曼氏无针乌贼更适宜在温度较高的环境下生存。盐度有利于海洋生物维持渗透压,对曼氏无针乌贼幼体而言,稳定的盐度环境有利于其发育和成熟[48]。该特性解释了SSS最小值对曼氏无针乌贼分布的影响较大,且在海水盐度较低的条件下曼氏无针乌贼存在概率较小的原因。此外,SSS最小值达到25‰ 时,曼氏无针乌贼存在概率最大,表明曼氏无针乌贼对盐度的耐受性较强。曼氏无针乌贼存在概率与pH变化关系表明,曼氏无针乌贼对pH的耐受范围较窄,仅在pH为7.8~8.3时的存在概率较大。面对敌害或外界环境条件突变时,曼氏无针乌贼会出现喷墨现象,该现象在pH过高时也会发生[28]。结合该现象和物种在pH达8.2后其存在概率降幅较大的现象可知,曼氏无针乌贼对pH的变化敏感。因此后续对曼氏无针乌贼潜在适生区进行保护时需注意海洋水质及其酸碱度。
3.2 曼氏无针乌贼适生区分布及其修复策略
MaxEnt模型预测的物种适生区分布是模型在一定的研究区域和时间尺度下根据已有的约束条件 (目标物种的存在数据及其所处环境信息) 而运算出的物种时空分布拟合[49],具有科学性和直观性,有利于为物种保护区规划研究和物种增殖放流选址提供参考。根据曼氏无针乌贼的适生区分布图及面积大小,可获得其在研究区域的分布特征。首先,曼氏无针乌贼在研究区域内的高适生区总体分布以嵊泗列岛、马鞍列岛和中街山列岛为核心。该特征与曼氏无针乌贼的繁殖特点有关,浙江北部的曼氏无针乌贼种群在需要进行交配和产卵时会游向浅海区和岛屿周围。舟山海域内岛屿众多,诸如嵊泗列岛、马鞍列岛和中街山列岛等岛屿为曼氏无针乌贼繁殖创造了优良条件。其次,曼氏无针乌贼适生区分布面积广,占研究区总面积的75.11%,大部分高适生区集中在舟山本岛的东北方海域,从马鞍列岛到舟山本岛东侧海域呈南北方向延伸,具有一定的连通性。且高适生区与次高适生区包含该海域的3个海洋特别保护区,其中在高适生区内包含2个国家级海洋特别保护区,分别为嵊泗马鞍列岛海洋保护区和普陀中街山列岛海洋生态特别保护区。次高适生区则包含了舟山市东部省级海洋保护区。利用该特征,加强对舟山海域内海洋特别保护区的保护,维持其生态系统平衡,有利于保护曼氏无针乌贼潜在适生区,从而间接促进其资源恢复。除自然条件外,资金和政策的支持也同样给予了曼氏无针乌贼资源恢复优势。2017年舟山市为海洋保护区提供法律支撑,出台了第一部国家级海洋特别保护区地方性法规《舟山国家级海洋特别保护区管理条例》。2018—2022年,为生态修复累计投入资金超过1亿元,陆续投放了近50万立方米的人工鱼礁[50]。因此,鉴于曼氏无针乌贼在舟山海域内适生区分布面积广、适生区保护程度高、资金投入充足及相关保护政策完善,可得出结论:舟山海域曼氏无针乌贼资源恢复的潜力较大。
根据适生区分布情况可采取相应策略进行曼氏无针乌贼增殖放流,从而促进其资源恢复。首先,可在曼氏无针乌贼产卵期选择高适生区投放曼氏无针乌贼成体,该措施有利于较大程度提高曼氏无针乌贼产量。原因在于,与其他海洋生物相比,头足类动物的生命周期相对较短,物种的数量变化受产卵数量和孵化数量的影响较大。并且曼氏无针乌贼产卵期较短,一般为4月中旬至6月下旬,亲体繁殖后不久便死亡[51]。其次,次高适生区、中适生区和低适生区3个适生等级的面积大小相近。从生境保护的角度出发,提高曼氏无针乌贼资源有两种策略,一是扩大曼氏无针乌贼高适生区面积,二是将次高适生区、中适生区和低适生区向更高一级的适生区转化。在面对次高适生区、中适生区和低适生区面积大小相近的情况时,选择次高适生区进行生境修复获得的效益较大。
3.3 MaxEnt模型在头足类中的应用
MaxEnt模型在头足类中的应用逐渐增加,本文以头足类曼氏无针乌贼为研究对象展现了MaxEnt模型在多元生态因子条件下模拟预测物种适生区分布及物种生境需求分析上的适用性。同时,与ArcGIS软件相结合获得了较为直观的结果,可为曼氏无针乌贼增殖放流选址提供参考。但该模型仍具有一定的局限性,在适生区等级划分上仍然具有一定的主观性。此外,头足类生物活动范围广,其分布不具有明确的分界线,而MaxEnt模型结果的准确性受样本量和研究范围的影响[52],因此需考虑研究区的边界范围。综上,应用MaxEnt模型分析头足类时,应注重模型参数优化和模型预测结果的精度检验。大部分头足类具有洄游分布的特点,根据其洄游规律划分时间尺度,再利用MaxEnt模型进行物种生境适宜性分析可获得更为精准的物种生境适宜区域分布情况。由于调查过程中受疫情影响,部分月份收集统计到的曼氏无针乌贼分布点较少。为了保证模型的预测精度,本文未将时间尺度划分至月份后再进行曼氏无针乌贼的生境适宜性分析。在今后的研究中,可进一步补充相关数据,探索不同时空尺度下曼氏无针乌贼的生境适宜性,或选择不同物种分布模型与MaxEnt模型的预测结果进行对比,多方位、多角度地为曼氏无针乌贼的资源养护提供科学参考和理论支持。
-
图 1 草鱼肝脏组织结构
注:a. 对照组草鱼肝脏 (400×);b. 高糖组草鱼肝脏 (400×);c. 对照组草鱼肝脏 (1 000×);d. 高糖组草鱼肝脏 (1 000×)。
Figure 1. Tissue slices of liver of grass carp
Note: a. Liver of grass carp (400×); b. Liver in Group H of grass carp (400×); c. Liver in control group of grass carp (1 000×); d. Liver in Group H of grass carp (1 000×).
图 2 草鱼中肠组织结构
注:a. 对照组草鱼中肠 (100×);b. 高糖组草鱼中肠 (100×);c. 对照组草鱼中肠 (400×);d. 高糖组草鱼中肠 (400×)。
Figure 2. Tissue slices of midgut of grass carp
Note: a. Midgut in Control group grass carp (100 ×); b. Midgut in Group H of grass carp (100×); c. Midgut in control group of grass carp (400×); d. Midgut in Group H of grass carp (400×).
图 4 草鱼肝脏糖代谢相关基因的表达
注:方柱上的不同大写字母表示同一时间点高糖组、对照组的相对表达量差异显著,不同小写字母表示同一组在不同时间点的相对表达量差异显著 (P<0.05)。
Figure 4. Relative expression of genes related to liver glucose metabolism in grass carp
Note: Different uppercase letters on the square columns indicate significant differences in relative expression levels between the high glucose group and the control group at the same time point, while different lowercase letters on the square columns indicate significant differences in relative expression levels of the same group at different time points (P<0.05).
表 1 实验饲料配方及常规成分分析
Table 1 Formulation and proximate composition of diets
主要成分
Essential component普通饲料
Normal
diet/%高糖饲料
High-carbohydrate
diet/%鸡肉粉 Poultry meal 5 5 普通蒸汽级鱼粉 Ordinary fish meal 15 15 豆粕 Bean meal 33.57 33.57 木薯生淀粉 Raw cassava starch 15 40 豆油 Soybean oil 3 3 细统糠 Fine bran 25 0 多维多矿预混料① Premix (1%) 1 1 磷酸二氢钙 Ca(H2PO4)2 2 2 氯化胆碱 Choline chloride (50%) 0.2 0.2 食盐 Common salt 0.2 0.2 抗氧化剂 Antioxidant 0.03 0.03 合计 Total 100 100 营养成分 Nutritional ingredient 粗蛋白 Primary protein 28.55 28.55 粗脂肪 Crude fat 5.22 5.32 钙 Ca 1.35 1.35 磷 P 1.41 1.41 赖氨酸 Lys 1.83 1.83 蛋氨酸 Met 0.53 0.53 注:① 多维多矿预混料为商用预混料,由淮安市禾丰饲料有限公司提供。 Note: ① Multi-mineral premix is a commercial premix provided by Huai'an Hefeng Feed Co. Ltd.. 表 2 草鱼糖代谢相关基因 qRT-PCR 引物
Table 2 qRT-PCR primers of grass carp glycometabolism genes
基因
Gene引物序列 (5'—3')
Primer sequence (5'–3')GenBank 登录号
GenBank ID参考文献
Referencegk F: GGACAATGGTGAACTGGAGGACTTC
R: TAGCAGCACAAGACGCACAAGCGU065314 [19] pk F: GCCGAGAAAGTCTTCATCGCACAG
R: CGTCCAGAACCGCATTAGCCACJQ951928 [20] pepck F: ATCGTCACGGAGAACCAA
R: CCTGAACACCAAACTTAGCAJQ898294.1 [20] g6pase F: AAGGACAGCAGGTAGAAGAGG
R: ACGGAAAACAAGAAGAGCAGKY742725.1 gys2 F: CCTCCAGTAACAACTCACAACA
R: CAGATAGATTGGTGGTTACGCJQ792167 [19] fas F: GATGGGTCTAAGCCTGATGG
R: GACACCCTGTGGACATTGAGCMG241310 [19] β-actin F: GCTATGTGGCTCTTGACTTCG
R: GGGCACCTGAACCTCTCATTM25013.1 [21] 表 3 高糖饲料对草鱼生长性能的影响
Table 3 Effects of dietary carbohydrate on growth performance of grass carp
项目
Item样本量/组
N/Group养殖时间
Culture time/d对照组
Control group高糖组
Group H初始体质量 Initial body mass/g 60 0 130.61±14.79 132.83±13.36 终末体质量 Final body mass/g 60 40 146.08±27.97 143.88±20.43 60 153.18±29.40 165.54±36.64* 80 168.71±36.73 195.53±48.70* 140 361.18±67.61 388.57±59.02** 体质量增长率 RF/% 60 40 109.98±15.64 110.18±9.33 60 121.24±4.21 133.28±4.09* 80 129.82±0.23 147.20±9.68* 140 271.78±13.07 297.51±15.85 体长 Body length/cm 60 40 20.74±1.29 20.13±1.05 60 20.34±1.21 20.75±1.41 80 21.21±1.71 21.58±2.16 140 29.27±0.82 29.22±1.09 体高 Body height/cm 60 40 4.90±0.37 4.96±0.37 60 4.88±0.44 5.20±0.51** 80 5.34±0.59 5.49±0.70 140 6.54±0.15 6.87±0.32** 肥满度 RF 60 40 1.61±0.13 1.76±0.32** 60 1.78±0.28 1.85±0.23 80 1.77±0.19 1.92±0.28** 140 1.48±0.13 1.72±0.17** 肝体比 RHI/% 30 40 1.94±0.47 2.19±0.62 60 1.68±0.61 2.14±0.60 80 2.17±0.27 2.29±0.49 140 1.23±0.24 1.39±0.14 脏体比 RVI/% 30 140 6.64±1.29 8.12±0.71** 肠系膜脂肪系数 RMFI/% 30 140 1.16±0.54 2.20±0.64** 空壳质量 WE/g 30 140 345.71±20.67 383.23±27.46** 注:*. 同一时期下高糖组与对照组相比差异显著 (P<0.05);**. 差异极显著 (P<0.01)。后表同此。 Note: *. There are significant differences between the high glucose group and the control group during the same period (P<0.05); **. Extremely significant differences (P<0.01). The same case in the following table. 表 4 高糖饲料对草鱼血糖、胰岛素、糖原的影响
Table 4 Effects of dietary carbohydrate on plasma glucose, insulin and glycogen of grass carp
项目
Item养殖时间
Culture time/d对照组
Control group高糖组
Group H血糖浓度 Blood glucose/(mmol·L−1) 40 10.95±0.99 10.99±1.17 60 11.46±2.15 13.65±1.63 80 9.77±1.77 10.79±1.78 140 13.39±0.63 12.78±1.00 胰岛素水平 Insulin level/(mIU·L−1) 140 9.23±1.18 10.31±1.17 肝糖原质量分数 Liver glycogen/(mg·g−1) 140 24.23±13.03 36.29±10.41* 肌糖原质量分数 Muscle glycogen/(mg·g−1) 140 0.15±0.06 0.24±0.08 -
[1] 草鱼产业发展报告[J]. 中国水产, 2021(2): 27-37. [2] WILSON R P. Utilization of dietary carbohydrate by fish[J]. Aquaculture, 1994, 124(1): 67-80.
[3] 许霄霄, 刘伟, 文华, 等. 高糖饲料对吉富罗非鱼生长性能、饲料利用和糖脂代谢的影响[J]. 南方水产科学, 2017, 13(5): 94-102. doi: 10.3969/j.issn.2095-0780.2017.05.013 [4] 齐野, 孙向军, 于刚, 等. 饲料可消化糖水平对宝石鲈餐后代谢的影响[J]. 南方水产, 2010, 6(2): 59-65. [5] PAN M Z, LIU D N, LIU J M, et al. Biotin alleviates hepatic and intestinal inflammation and apoptosis induced by high dietary carbohydrate in juvenile turbot (Scophthalmus maximus L.)[J]. Fish Shellfish Immunol, 2022, 130(8): 560-571.
[6] 陈少莲, 刘肖芳. 我国淡水优质草食性鱼类的营养和能量学研究: 草鱼、团头鲂对七种水生高等植物的最大摄食量和消化率的测定[J]. 水生生物学报, 1993, 17(1): 1-12. [7] 张蕾, 章文琪, 吴仁福, 等. 中华绒螯蟹成蟹养殖池塘常用水草的营养成分比较[J]. 浙江海洋学院学报(自然科学版), 2016, 35(2): 113-121. [8] TIAN L X, LIU Y J, YANG H J, et al. Effects of different dietary wheat starch levels on growth, feed efficiency and digestibility in grass carp (Ctenopharyngodon idella)[J]. Aquac Int, 2012, 20(2): 283-293. doi: 10.1007/s10499-011-9456-6
[9] SU J Z, MEI L Y, XI L W, et al. Responses of glycolysis, glycogen accumulation and glucose-induced lipogenesis in grass carp and Chinese longsnout catfish fed high-carbohydrate diet[J]. Aquaculture, 2021, 533: 736146. doi: 10.1016/j.aquaculture.2020.736146
[10] CAI W J, LIANG X F, YUAN X C, et al. Different strategies of grass carp (Ctenopharyngodon idella) responding to insufficient or excessive dietary carbohydrate[J]. Aquaculture, 2018, 497: 292-298. doi: 10.1016/j.aquaculture.2018.07.042
[11] LI X S, ZHU X M, HAN D, et al. Carbohydrate utilization by herbivorous and omnivorous freshwater fish species: a comparative study on gibel carp (Carassius auratus gibelio. var CAS III) and grass carp (Ctenopharyngodon idellus)[J]. Aquac Res, 2016, 47(1): 128-139. doi: 10.1111/are.12476
[12] 曹俊明, 关国强, 刘永坚, 等. 饲料蛋白质、脂肪、碳水化合物水平对草鱼生长和组织营养成分组成的影响[J]. 水产科技情报, 1997(2): 8-12. doi: 10.16446/j.cnki.1001-1994.1997.02.002 [13] YUAN X C, LIANG X F, LI A X, et al. The feedback regulation of carbohydrates intake on food intake and appetite in grass carp (Ctenopharyngodon idella)[J]. Fish Physiol Biochem, 2021, 47(5): 1395-1403. doi: 10.1007/s10695-020-00914-5
[14] 樊佳佳, 唐小红, 白俊杰, 等. 草鱼PKMa基因SNPs筛选及与耐糖性状的关联分析[J]. 农业生物技术学报, 2019, 27(6): 1072-1080. [15] 孙雪. 草鱼生长相关SNPs标记的筛选及优势基因型的聚合效果分析[D]. 上海: 上海海洋大学, 2020: 15. [16] MARANDEL L, SEILIEZ I, VERON V, et al. New insights into the nutritional regulation of gluconeogenesis in carnivorous rainbow trout (Oncorhynchus mykiss): a gene duplication trail[J]. Physiol Genomics, 2015, 47(7): 253-263. doi: 10.1152/physiolgenomics.00026.2015
[17] BOONANUNTANASARN S, KUMKHONG S, YOOHAT K, et al. Molecular responses of Nile tilapia (Oreochromis niloticus) to different levels of dietary carbohydrates[J]. Aquaculture, 2018, 482: 117-123. doi: 10.1016/j.aquaculture.2017.09.032
[18] 陈团, 胡毅, 张德洪, 等. 不同糖源膨化饲料对大规格草鱼生长、越冬及血清部分生化指标的影响[J]. 水产学报, 2019, 43(4): 1069-1079. [19] FANG L, GUO X, LIANG X F. First feeding of grass carp (Ctenopharyngodon idellus) with a high-carbohydrate diet: the effect on glucose metabolism in juveniles[J]. Aquac Rep, 2021, 21: 100830. doi: 10.1016/j.aqrep.2021.100830
[20] LI A X, YUAN X C, LIANG X F, et al. Adaptations of lipid metabolism and food intake in response to low and high fat diets in juvenile grass carp (Ctenopharyngodon idellus)[J]. Aquaculture, 2016, 457: 43-49. doi: 10.1016/j.aquaculture.2016.01.014
[21] GONG L, HE H C, LI D J, et al. A new isolate of Pediococcus pentosaceus (SL001) with antibacterial activity against fish pathogens and potency in facilitating the immunity and growth performance of grass carps[J]. Front Microbiol, 2019, 10: 1384. doi: 10.3389/fmicb.2019.01384
[22] GAO W, LIU Y J, TIAN L X, et al. Effect of dietary carbohydrate-to-lipid ratios on growth performance, body composition, nutrient utilization and hepatic enzymes activities of herbivorous grass carp (Ctenopharyngodon idella)[J]. Aquac Nutr, 2010, 16(3): 327-333.
[23] 胡毅, 陈云飞, 张德洪, 等. 不同碳水化合物和蛋白质水平膨化饲料对大规格草鱼生长、肠道消化酶及血清指标的影响[J]. 水产学报, 2018, 42(5): 777-786. [24] SONG X R, LIU H K, HAN D, et al. Two strains of gibel carp (Carassius gibelio) exhibit diverse responses to carbohydrates in a low-lipid diet[J]. Aquac Nutr, 2022, 2022: 9417331. DOI: 10.1155/2022/9417331.
[25] SONG X R, MARANDEL L, SKIBA-CASSY S, et al. Regulation by dietary carbohydrates of intermediary metabolism in liver and muscle of two isogenic lines of rainbow trout[J]. Front Physiol, 2018, 9: 1579. doi: 10.3389/fphys.2018.01579
[26] YUAN X C, ZHOU Y, LIANG X F, et al. Molecular cloning, expression and activity of pyruvate kinase in grass carp Ctenopharyngodon idella: Effects of dietary carbohydrate level[J]. Aquaculture, 2013, 410/411: 32-40. doi: 10.1016/j.aquaculture.2013.06.009
[27] WANG J T, LI X Y, HAN T, et al. Effects of different dietary carbohydrate levels on growth, feed utilization and body composition of juvenile grouper Epinephelus akaara[J]. Aquaculture, 2016, 459: 143-147. doi: 10.1016/j.aquaculture.2016.03.034
[28] LI X F, LIU W B, LU K L, et al. Dietary carbohydrate/lipid ratios affect stress, oxidative status and non-specific immune responses of fingerling blunt snout bream, Megalobrama amblycephala[J]. Fish Shellfish Immunol, 2012, 33(2): 316-323. doi: 10.1016/j.fsi.2012.05.007
[29] ABIMORAD E G, CARNEIRO D J, URBINATI E C. Growth and metabolism of pacu (Piaractus mesopotamicus Holmberg 1887) juveniles fed diets containing different protein, lipid and carbohydrate levels[J]. Aquac Res, 2007, 38(1): 36-44. doi: 10.1111/j.1365-2109.2006.01621.x
[30] VASQUEZ-TORRES W, ARIAS-CASTELLANOS J A. Effect of dietary carbohydrates and lipids on growth in cachama (Piaractus brachypomus)[J]. Aquac Res, 2013, 44(11): 1768-1776.
[31] YE W J, TAN X Y, CHEN Y D, et al. Effects of dietary protein to carbohydrate ratios on growth and body composition of juvenile yellow catfish, Pelteobagrus fulvidraco (Siluriformes, Bagridae, Pelteobagrus)[J]. Aquac Res, 2009, 40(12): 1410-1418. doi: 10.1111/j.1365-2109.2009.02239.x
[32] TIAN J, WU F, YU L J, et al. The effects of high-macronutrient (protein, fat and carbohydrate) diets on growth performance and muscular metabolic responses in grass carp[J]. Aquac Nutr, 2020, 26(6): 2135-2146. doi: 10.1111/anu.13152
[33] WANG J L, LU R H, SUN J J, et al. Differential expression of lipid metabolism-related genes and miRNAs in Ctenopharyngodon idella liver in relation to fatty liver induced by high non-protein energy diets[J]. Aquac Nutr, 2017, 48(8): 4070-4085.
[34] 付兵, 郁二蒙, 王广军, 等. 投喂蚕豆对草鱼肝脏脂肪蓄积及脂肪代谢的影响[J]. 上海海洋大学学报, 2020, 29(1): 45-54. doi: 10.12024/jsou.20190302546 [35] CAHU C L, GISBERT E, VILLENEUVE L, et al. Influence of dietary phospholipids on early ontogenesis of fish[J]. Aquac Res, 2009, 40(9): 989-999. doi: 10.1111/j.1365-2109.2009.02190.x
[36] LI X F, XU C, ZHANG D D, et al. Molecular characterization and expression analysis of glucokinase from herbivorous fish Megalobrama amblycephala subjected to a glucose load after the adaption to dietary carbohydrate levels[J]. Aquaculture, 2016, 459: 89-98. doi: 10.1016/j.aquaculture.2016.03.035
[37] SONG X R, MARANDEL L, DUPONT-NIVET M, et al. Hepatic glucose metabolic responses to digestible dietary carbohydrates in two isogenic lines of rainbow trout[J]. Biol Open, 2018, 7(6): bio032896.
[38] 徐晶, 梁旭方, 蔡文静, 等. 草鱼3个6-磷酸葡萄糖酶催化亚基的基因表达分析及高糖饲料对其表达的影响[J]. 中国水产科学, 2020, 27(1): 24-34. [39] QIANG J, YANG H, HE J, et al. Comparative study of the effects of two high-carbohydrate diets on growth and hepatic carbohydrate metabolic enzyme responses in juvenile GIFT tilapia (Oreochromis niloticus)[J]. Turkish J Fish Aquat, 2014, 14(2): 515-525.
[40] SHI H J, LIU W B, XU C, et al. Transcriptional regulation of the AMP-activated protein kinase and glycolipid metabolism-related genes by insulin and glucagon in blunt snout bream (Megalobrama amblycephala): a comparative study[J]. Aquaculture, 2020, 515: 734553. doi: 10.1016/j.aquaculture.2019.734553
-
期刊类型引用(3)
1. 吴春梅,杨柳明,王维奇,罗欣婕. 闽江河口湿地围垦后土壤微生物残体碳特征. 环境科学学报. 2024(12): 360-368 . 百度学术
2. 洪妍,杨平,仝川,赵光辉,李玲,唐晨,张彦,谭莹莹,阮蔓菁. 亚热带河口区对虾养殖池塘浮游植物初级生产力变化. 湖泊科学. 2022(03): 881-893 . 百度学术
3. 林于蓝,陈钰,尹晓雷,曾庆淞,王维奇. 围垦养殖与退塘还湿对闽江河口湿地土壤铁碳结合特征的影响. 环境科学学报. 2022(07): 466-477 . 百度学术
其他类型引用(4)