Study on relationship between albacore catch rate and water depth and temperature in South Indian Ocean
-
摘要: 印度洋金枪鱼延绳钓渔业是我国远洋渔业的重要组成部分,海洋不同深度的水温影响到长鳍金枪鱼 (Thunnus alalunga) 延绳钓渔获率。文章利用2008—2017年延绳钓生产作业数据,并结合Argo浮标水温数据,采用广义加性模型 (Generalized additive model, GAM) 分析长鳍金枪鱼空间分布与不同深度水温之间的关系。结果表明,海表面 (0 m)、200和400 m 3个水层的温度显著影响长鳍金枪鱼的空间分布,最优的GAM模型对渔获率 (单位捕捞努力量渔获量,Catch per unit effort, CPUE) 的方差解释率为53.3%,模型拟合的决定系数为0.527。长鳍金枪鱼渔获率与所选取的3个水层温度均呈非线性关系,高渔获区集中分布于17~30 ℃的表层海域,17~20 ℃的200 m层海域,9~15 ℃的400 m层海域,以及他们的交集海线。文章初步得出了南印度洋长鳍金枪鱼空间分布与水深断面温度的关系,可为指导长鳍金枪鱼的合理生产提供技术支撑。Abstract: Longline tuna fisheries in the Indian Ocean are commerically important for Chinese offshore fisheries. Water temperature at different depths can affect longline albacore catch rates. In this paper, the generalized additive model (GAM) was used to analyze the relationship between the the catch rate of albacore (Thunnus alalunga) and water temperature at different depths based on the fishery data and corresponding Argo buoy data during 2008−2017. The results show that the catch rate of albacore tuna catch rate was significantly affected by the sea surface temperature (0 m), water temperature at depth of 200 and 400 m. The optimal GAM model explained the variance of catch rate (Catch per unit effort) by 53.3%, and the determining coefficient of model was 0.527. The catch rates of albacore tuna had a nonlinear relationship with the temperature of the three selected water layers. High catch areas were concentrated at the surface layer of 17−30 ℃, 200 m deep sea area of 17−20 ℃, and 400 m deep sea area of 9−15 ℃, and their intersection areas. We have derived the relationship between the spatial distribution of albacore tuna catch rate in the South Indian Ocean and the temperature of three water depth sections for the first time, and the results provide technical support for guiding the rational production of albacore tuna in the Indian Ocean.
-
Keywords:
- Thunnus alalunga /
- Generalized additive model /
- Catch rate /
- South Indian Ocean
-
弧菌广泛存在于养殖水体,是海水养殖动物的主要病原菌,弧菌病是目前为止对养殖鱼虾危害最大,造成损失最严重的细菌性疾病之一。弧菌病可以发生在养殖的各个时期,因此防治对虾弧菌病是养殖成败的关键因素。防治细菌性疾病最常用的方法是使用抗菌素,但频繁无节制的使用药物,导致耐药性菌株的产生并且造成环境的污染,因此,抗菌素的有效替代品的研究显得尤为迫切,应用生物及生态法防治疾病日益受到重视。大量的研究表明[1-7],乳酸菌能够调节机体肠道正常菌群,保持微生态平衡,提高食物消化率和生物价,增强机体的免疫功能,提高机体的抗病能力,产生抑菌活性代谢产物,如乳酸菌肽、细菌素、乳酸、过氧化氢、乙酸等,对许多革兰氏阳性菌李斯特氏菌、芽孢菌、梭菌等及革兰氏阴性菌大肠杆菌等有强烈的抑制作用,可抑制肠道内腐败菌的生长繁殖和腐败产物的产生,乳酸菌被作为饲料添加剂而广泛应用于禽畜养殖中,防治腹泻、下痢、肠炎等肠道功能紊乱的许多疾病。乳酸菌益生素作为鱼、虾饲料添加剂也受到广泛的研究,有提高养殖动物的免疫力,抵御病原菌的侵袭,提高养殖成活率的效果。体外拮抗实验是筛选益生菌的重要步骤,本实验通过乳酸菌体外对致病弧菌的拮抗作用研究,旨在筛选乳酸菌有益菌株,为进一步在养殖生产中应用益生素产品防治病害提供理论基础及水产养殖动物的病害防治提供一种生物防治方法。
1. 材料与方法
1.1 试验菌株
测试菌:乳酸杆菌L1是经点接种法初步筛选出的对弧菌有拮抗作用的菌株。
指示菌:溶藻弧菌T1是由本所鱼病室提供的从患病军曹鱼中分离,并经回归感染确认有致病性的菌株。鲨鱼弧菌T2也是从患病军曹鱼中分离。
1.2 培养基
乳酸菌培养基为改良的MRS培养基:蛋白胨10 g,酵母膏5 g,牛肉膏10 g,葡萄糖20 g,无水乙酸钠3 g,柠檬酸三铵2 g,K2HPO4 2 g,MgSO4 · 7H2O 0.2 g,MnSO4 · H2O 0.05 g,水1 000 mL,pH 6.8。
弧菌培养基为2216E培养基。
检测弧菌培养基为TCBS培养基。
1.3 测定方法
抑菌活性检测采用平板打孔抑菌圈测定法。
平板打孔抑菌圈测定法:制作2216E培养基,无菌操作倒平板,每个平板的培养基厚度为6 mm,取指示菌弧菌菌悬液0.1 mL,涂布于2216E培养基表面,选择一定的位置,在无菌条件下,用无菌打孔器打孔,孔径为10 mm,将乳酸杆菌发酵液注入孔内(不能溢出),于30℃培养24 h,检测抑菌圈直径的大小。以相同pH值的HCl水作对照。
2. 结果
2.1 乳酸杆菌及其代谢产物对弧菌的抑制作用
表 1为L1接种MRS培养液,经24 h培养后的发酵液及发酵液经5 000 rpm,15 min离心后取得的上清液,采用平板打孔法测定抑菌圈的大小。
表 1 L1菌发酵液及离心后上清液对指示菌T1、T2的抑菌圈大小Table 1 Size of zone plate inhibitory of fermentation liquid and supernatant liquid of L1 strain to indicator strains T1 and T2(mm) 指示菌
indicator strains发酵液(pH 3.5)
fermentation liquid上清液(pH 3.5)
supernatant liquid对照HCl水组(pH 3.5)
HCl diluted liquidT1 23 20 0 T2 18 16 0 从表 1试验可看出发酵液与上清液对指示菌T1、T2都有抑菌圈,而相同pH值的HCl水却没有,可见发酵液与上清液的抑菌效果并不是pH值低造成的。试验结果表明发酵液比上清液对弧菌的抑菌圈大,说明发酵液的抑菌活性强于上清液,乳酸杆菌及其代谢产物对弧菌有协同抑制作用。发酵液对弧菌T1、T2的抑菌圈大小不同,上清液也是如此,表明L1发酵液的代谢产物对不同种类弧菌的抑菌活性不同。
2.2 不同稀释倍数发酵液及离心后上清液的抑菌活性
L1接种MRS,30℃培养24 h,发酵液用无菌生理盐水释稀1、2、3倍,其稀释液对弧菌的抑菌能力见表 2。其发酵液经5 000 rpm,15 min离心,取得上清液,上清液用无菌生理盐水释稀1、2、3倍,其稀释液对弧菌的抑菌能力见表 3。
表 2 不同稀释倍数L1菌发酵液对指示菌T1、T2的抑菌圈大小Table 2 Size of zone plate inhibitory of diluted fermentation liquid of L1 strain to indictor strains T1 and T2(mm) 指示菌
indicator strains发酵液稀释倍数 times of dilution for fermentation liquid 0 1 2 3 T1 25 23 21.5 13.5 T2 23 21 19 无 表 3 不同稀释倍数L1菌发酵液对指示菌T1、T2的抑菌圈大小Table 3 Size of zone plate inhibitory of diluted supernatant liquid of L1 strain to indictor strains T1 and T2(mm) 指示菌
indicator strains上清液稀释倍数 times of dilution for supernatant liquid 0 1 2 3 T1 18 14.5 无 无 T2 15 无 无 无 从表 2、3可看出发酵液经3倍稀释对弧菌T1仍有抑制作用,上清液经1倍稀释后对弧菌T1的抑菌圈大小与发酵液3倍稀释液抑菌圈的大小相近。
2.3 不同培养时间发酵液的抑菌活性
L1接种MRS培养基,30℃培养,培养18 h、24 h、36 h、48 h、72 h、5 d的发酵液的抑菌圈大小见表 4。
表 4 不同培养时间的L1菌发酵液对指示菌T1、T2的抑菌圈大小Table 4 Size of zone plate inhibitory to indictor strains T1 and T2 of fermentation liquid of L1 strain under different time(mm) 指示菌
indicator strains不同培养时间的发酵液 fermentation liquid of different culture time 18 h 24 h 36 h 48 h 72 h 5 d T1 无 22.5 24 24 25 28 T2 无 20 22 23 24 25 从表 4可看出,在菌生长的对数期18 h时的发酵液对两种指示菌都无抑菌圈,随培养时间的延长,抑菌圈越来越大,生长衰退期(36 h)比稳定期(24 h)的抑菌圈大,5 d的陈培养液的抑菌圈最大。原因可能是随着培养时间的延长,乳酸杆菌的生长及代谢产物抑菌物质的积累提高了其杀菌作用;陈培养物和衰退期的发酵液的抑菌圈大于生长期及稳定期的抑菌圈,可能是一些抗菌活性物质在乳酸菌死亡后释放出来。
2.4 发酵液抑菌活性物质的耐热性试验
将MRS 24 h的发酵液,于60、80℃恒温水浴中保温15 min,及沸水浴中保温5、10、15 min,与原发酵液一起作抑菌活性试验,试验结果见表 5。
表 5 不同温度及时间处理的L菌发酵液对指示菌T1、T2的抑菌圈大小Table 5 Size of zone plate inhibitory to indictor strains T1 and T2 of fermentation liquid of L1 strain treated under different temperature and time(mm) 指示菌
indictor strains60℃ 80℃ 100℃ 未处理组
control group15 min 15 min 5 min 10 min 15 min T1 23 23 23 23 23 23 T2 22 22 22 22 22 22 从表 5的试验结果可见,发酵液经60、80℃恒温水浴处理15 min,及沸水浴中处理5、10、15 min,其抑菌圈大小与未作温度处理的发酵液的抑菌圈大小无差异,可见代谢抑菌物质具有很好的耐热性。对虾饲料制作过程有高温制粒这一步骤,瞬间温度可达100℃,益生菌代谢抑菌物质具有良好的耐热性,可避免对虾饲料制作过程瞬间高温的破坏,而保持其活性。
3. 讨论
体外抑菌试验是筛选益生菌的手段,通过L1菌株体外对弧菌抑制试验,结果表明L1菌株对致病弧菌T1、T2有较强的抑制作用,而且生长速度快,有望成为微生态制剂的候选菌种,但能否成为益生菌,还需做一系列的的工作如进一步确定其对宿主有否致病作用,对被选菌株进行致病性评估;对宿主潜在效果的评价等。
黄沧海等报导乳酸杆菌的代谢产物对不同血清型的大肠杆菌的抑制作用存在一定的差别,本试验的结果也表明,乳酸杆菌L1对不同的弧菌有不同的抑菌活性。水产养殖动物不同的疾病因病原不同,因此,实际在应用益生素制剂时,如果没有了解清楚益生菌的适用对象和范围就可能就会有不同的效果。Gatesoupe(1999)[8]通过每天添加乳杆菌在大菱鲆幼体活体食物的轮虫的培养基中,大大提高了鱼的成活率,当病原性弧菌侵袭幼体时,添加的乳杆菌能够大大降低幼体在9 d前的死亡率。因此认为乳杆菌可以防御病原菌-弧菌入侵大比目鱼的幼体。而Gildberg等(1997)[9-10]用添加了产乳酸细菌的饲料喂养大西洋鳕,将它们与经腹膜内感染了气单胞菌(Aeromonas salmonicida)的鱼苗一起养殖,在以后4周时间内记录鱼的死亡率,结果表明产乳酸细菌作为鱼苗饲料的添加成分可以促进肠微生物的定植,但未出现防止气单胞菌感染的现象,与预想相反,在饲料中添加了产乳酸细菌的鱼苗的死亡率最高。
本试验的结果表明乳酸杆菌与代谢产物有协同抑菌作用。因此在评价乳酸菌的抑菌效果时,以乳酸杆菌与代谢产物的协同抑菌效果作为衡量指标应当是更科学合理。在使用乳酸菌益生素产品时,不应当只是利用其菌体,菌体及其代谢产物能一起使用效果可能会更好。
-
表 1 广义加性模型统计参数表
Table 1 Statistical parameters of GAM model
公式
Formula赤池准则
信息
AIC方差
解释率
Deviance
explained决定
系数
R2广义交叉
验证
GCVlog(CPUE)=NULL 15023.60 0 0 2.4984 log(CPUE)=s(Year) 15005.51 0.8 0.00642 2.4871 log(CPUE)=s(Year)+s(Month) 14690.30 8.6 0.0831 2.2988 log(CPUE)=s(Year)+s(Month)+s(Lat, k=5) 13173.10 37.6 0.373 1.5735 log(CPUE)=s(Year)+s(Month)+s(Lat, k=5)+s(Lon) 12602.84 46.1 0.458 1.3645 log(CPUE)=s(Year)+s(Month)+s(Lat, k=5)+s(Lon)+s(Temp_0) 12156.86 52.0 0.516 1.2207 log(CPUE)=s(Year)+s(Month)+s(Lat, k=5)+s(Lon)+s(Temp_0)+s(Temp_50) 12153.27 52.1 0.516 1.2196 log(CPUE)=s(Year)+s(Month)+s(Lat, k=5)+s(Lon)+s(Temp_0)+s(Temp_50)+s(Temp_100) 12063.03 53.4 0.528 1.1925 log(CPUE)=s(Year)+s(Month)+s(Lat, k=5)+s(Lon)+s(Temp_0)+s(Temp_50)+s(Temp_100)+s(Temp_200) 12050.30 53.4 0.528 1.1925 log(CPUE)=s(Year)+s(Month)+s(Lat, k=5)+s(Lon)+s(Temp_0)+
s(Temp_50)+s(Temp_100)+s(Temp_200)+s(Temp_300)12007.99 54.3 0.536 1.1763 log(CPUE)=s(Year)+s(Month)+s(Lat, k=5)+s(Lon)+s(Temp_0)+s(Temp_50)+s(Temp_100)+
s(Temp_200)+s(Temp_300)+s(Temp_400)11994.66 54.5 0.538 1.1724 表 2 广义加性模型模型分析结果
Table 2 Analysis results of GAM model
变量
Variable自由度
dfF P Year 8.588 6.594 0.0003 Month 7.3419 7.332 <0.0001 Lat 4.5165 56.543 <0.0001 Lon 9.1298 37.893 <0.0001 Temp_0 9.2925 14.781 <0.0001 Temp_50 0.9521 2.840 0.0210 Temp_100 11.516 0 9.446 <0.0001 Temp_200 4.8355 3.394 0.0051 Temp_300 7.8455 1.301 0.1981 Temp_400 4.8479 3.261 0.0019 -
[1] 马璐璐. 印度洋长鳍金枪鱼资源的时空分布和资源评估[D]. 上海: 上海海洋大学, 2017: 2. [2] 马璐璐, 朱江峰, 耿喆, 等. 运用生物量动态模型评估印度洋长鳍金枪鱼资源[J]. 上海海洋大学学报, 2018, 27(2): 259-264. [3] 张亚男, 官文江, 李阳东. 印度洋长鳍金枪鱼栖息地指数模型的构建与验证[J]. 上海海洋大学学报, 2020, 29(2): 268-279. [4] DHURMEEA Z, ZUDAIRE I, CHASSOT E, et al. Reproductive biology of albacore tuna (Thunnus alalunga) in the Western Indian Ocean[J]. PLoS One, 2016, 11(12): e168605.
[5] 张艳波. 东南太平洋长鳍金枪鱼垂直分布与栖息地指数研究[D]. 上海: 上海海洋大学, 2015: 26. [6] 张嘉容, 杨晓明, 田思泉. 基于最大熵模型的南太平洋长鳍金枪鱼栖息地预测[J]. 中国水产科学, 2020, 27(10): 1222-1233. [7] 闫敏, 张衡, 樊伟, 等. 南太平洋长鳍金枪鱼渔场CPUE时空分布及其与关键海洋环境因子的关系[J]. 生态学杂志, 2015, 34(11): 3191-3197. [8] 唐衍力, 徐湛. 中西太平洋长鳍金枪鱼 (Thunnus alalunga) 延绳钓渔场分布的时空变化[J]. 渔业科学进展, 2015, 36(3): 22-29. doi: 10.11758/yykxjz.20150304 [9] DOMOKOS R, SEKI M P, POLOVINA J J, et al. Oceanographic investigation of the American Samoa albacore (Thunnus alalunga) habitat and longline fishing grounds[J]. Fish Oceanogr, 2007, 16(6): 555-572. doi: 10.1111/j.1365-2419.2007.00451.x
[10] 郭刚刚, 张胜茂, 樊伟, 等. 南太平洋长鳍金枪鱼垂直活动水层空间分析[J]. 南方水产科学, 2016, 12(5): 123-130. doi: 10.3969/j.issn.2095-0780.2016.05.016 [11] 杨胜龙, 张忭忭, 唐宝军, 等. 基于GAM模型分析水温垂直结构对热带大西洋大眼金枪鱼渔获率的影响[J]. 中国水产科学, 2017, 24(4): 875-883. [12] 张嘉容, 杨晓明, 戴小杰, 等. 南太平洋长鳍金枪鱼延绳钓渔获率与环境因子的关系研究[J]. 南方水产科学, 2020, 16(1): 69-77. [13] WANG J, CHEN X, STAPLES K W, et al. The skipjack tuna fishery in the west-central Pacific Ocean: applying neural networks to detect habitat preferences[J]. Fish Sci, 2018, 84(2): 309-321. doi: 10.1007/s12562-017-1161-6
[14] GUISAN A, EDWARDS T C, HASTIE T. Generalized linear and generalized additive models in studies of species distributions: setting the scene[J]. Ecol Model, 2002, 157(2): 89-100.
[15] 武胜男, 陈新军, 刘祝楠. 基于GAM的西北太平洋日本鲭资源丰度预测模型建立[J]. 海洋学报 (中文版), 2019, 41(8): 36-42. [16] 张云雷, 徐宾铎, 张崇良, 等. 基于Tweedie-GAM模型研究海州湾小黄鱼资源丰度与栖息环境的关系[J]. 海洋学报, 2019, 41(12): 78-89. [17] 王晓晴, 林宇, 吴锦仁. 东南印度洋长鳍金枪鱼延绳钓渔场探捕浅析[J]. 渔业信息与战略, 2020, 35(3): 198-207. [18] 陈雪忠, 樊伟, 崔雪森, 等. 基于随机森林的印度洋长鳍金枪鱼渔场预报[J]. 海洋学报 (中文版), 2013, 35(1): 158-164. [19] 郭爱, 陈新军. ENSO与中西太平洋金枪鱼围网资源丰度及其渔场变动的关系[J]. 海洋渔业, 2005(4): 338-342. doi: 10.3969/j.issn.1004-2490.2005.04.015 [20] 陈新军, 郑波. 中西太平洋金枪鱼围网渔业鲣鱼资源的时空分布[J]. 海洋学研究, 2007(2): 13-22. doi: 10.3969/j.issn.1001-909X.2007.02.002 [21] 史珩瑜, 张祖强, 任宏利. 近百年来ENSO强度的变化特征[J]. 气候变化研究进展, 2017, 13(1): 1-10. doi: 10.12006/j.issn.1673-1719.2016.117 [22] 孙诗. 南印度洋长鳍金枪鱼黄鳍金枪鱼生物学特性及金枪鱼延绳钓捕捞技术初步分析[D]. 舟山: 浙江海洋大学, 2019: 54. [23] 范永超, 陈新军, 汪金涛. 基于多因子栖息地指数模型的南太平洋长鳍金枪鱼渔场预报[J]. 海洋湖沼通报, 2015(2): 36-44. [24] 刘洪生, 蒋汉凌, 戴小杰. 中西太平洋长鳍金枪鱼渔场与海温的关系[J]. 上海海洋大学学报, 2014, 23(4): 602-607. [25] 范江涛, 陈新军, 钱卫国, 等. 瓦努阿图周边海域长鳍金枪鱼渔场分布及其与表温关系[J]. 海洋湖沼通报, 2011(1): 71-78. doi: 10.3969/j.issn.1003-6482.2011.01.010 [26] MACDONALD J I, FARLEY J H, CLEAR N P, et al. Insights into mixing and movement of South Pacific albacore Thunnus alalunga derived from trace elements in otoliths[J]. Fish Res, 2013, 148: 56-63. doi: 10.1016/j.fishres.2013.08.004
[27] SAITO S. Studies on fishing of albacore, Thunnus alalunga (Bonnaterre) by experimental deep-sea tuna long-line[J]. Mem Fac Fish Hokkaido Univ, 1973, 21(2): 107-184.
[28] 陈雪忠, 杨胜龙, 张禹, 等. 热带印度洋大眼金枪鱼垂直分布空间分析[J]. 中国水产科学, 2013, 20(3): 660-671. [29] 储宇航, 戴小杰, 田思泉, 等. 南太平洋长鳍金枪鱼栖息水层深度的时间变化与空间分布[J]. 海洋通报, 2016, 35(2): 216-224. doi: 10.11840/j.issn.1001-6392.2016.02.013 [30] 曹晓怡, 周为峰, 樊伟, 等. 大眼金枪鱼渔场与环境关系的研究进展[J]. 海洋渔业, 2008(2): 176-182. doi: 10.3969/j.issn.1004-2490.2008.02.012 [31] 翟天晨, 戴小杰, 朱江峰. 东太平洋长鳍金枪鱼个体大小与钓获深度的关系[J]. 海洋渔业, 2015, 37(1): 10-16. doi: 10.3969/j.issn.1004-2490.2015.01.002 [32] 宋婷婷, 樊伟, 伍玉梅. 卫星遥感海面高度数据在渔场分析中的应用综述[J]. 海洋通报, 2013, 32(4): 474-480. doi: 10.11840/j.issn.1001-6392.2013.04.017 -
期刊类型引用(0)
其他类型引用(1)