Modeling habitat of skipjack tuna of free swimming school in Western and Central Pacific Ocean based on MaxEnt model
-
摘要:
由于漂流人工集鱼装置 (Fish aggregating device, FAD) 的大量使用对金枪鱼种群带来的负面效应，金枪鱼围网渔业转向捕捞自由群成为发展趋势，因此开展鲣 (Katsuwonus pelamis) 自由群的栖息地利用研究非常必要。使用2016—2020年中西太平洋渔业委员会 (Western and Central Pacific Fisheries Commission, WCPFC) 统计的月度金枪鱼渔业数据和不同层的水温 (SST、Temp200)、海表盐度 (SSS)、溶解氧浓度 (DO0、DO50、DO200)、东西向海流速度 (EV)、南北向海流速度 (NV)、混合层深度 (MLD)、叶绿素 a 浓度 (CHL0、CHL50、CHL100、CHL200) 共13个环境变量，通过最大熵模型 (Maximum Entropy Model, MaxEnt) 模拟鲣自由群的栖息地分布及其月变化规律。结果表明：模型测试和训练集的AUC值及灵敏度值均大于0.90，真实技巧统计值大于0.80，模型具有很强的预测能力，可用于鲣的栖息地适宜性模拟；SST和DO200是影响鲣自由群栖息地偏好的关键因子，最适范围分别为30~31 ℃、114~153 mmol·m−3。研究期内，鲣自由群高适宜栖息地主要靠近巴布亚新几内亚和所罗门群岛海域，不同时期向东延伸的范围有较大变化，经度差达到6°。研究结果可为中国金枪鱼围网船队进行鲣自由群中心渔场的预报提供参考。
Abstract:Due to the negative effects of extensive use of drifting artificial fish aggregating devices (FADs) on tuna stocks, tuna purse seine fishing has become a development trend towards catching free swimming school, so it is necessary to specifically study the habitat use of free swimming school of skipjack tuna (Katsuwonus pelamis). In this study, we used monthly tuna fishery data from the Western and Central Pacific Fisheries Commission (WCPFC) from 2016 to 2020, different layers of water temperature (SST, Temp200), sea surface salinity (SSS), dissolved oxygen concentration (DO0, DO50, DO200), east-west current velocity (EV), north-south current velocity (NV), mixed layer depth (MLD), chlorophyll a concentration (CHL0, CHL50, CHL100, CHL200), and a total of 13 environmental variables by Maximum Entropy (MaxEnt) model to simulate the habitat distribution of the free swimming school and their monthly variation patterns. The results show that the AUC and sensitivity values of both the test and training set of the model were greater than 0.90, and the true skill statistics values were greater than 0.80, indicating that the model has strong predictive ability and can be used for the habitat suitability modeling of skipjack tuna. SST and DO200 were the key factors affecting the habitat preference of free swimming school, with the optimal ranges of 30−31 ℃ and 114−153 mmol·m−3, respectively. During the survey period, the highly suitable habitat for free swimming school was mainly near the waters of Papua New Guinea and Solomon Islands, with a large variation in the range extending eastward in different periods, and the difference in longitude reached 6°. The results provide references for the prediction of the central fishing ground of free swimming school of skipjack tuna by Chinese tuna purse seine fleet.
-
几丁质由N-乙酰-D-氨基葡萄糖 (GlcNAc) 通过β-1,4-糖苷键连接而成[1]，在自然界中广泛分布于虾和蟹的外壳[2]、真菌细胞壁、昆虫外骨骼以及一些绿藻中。仅水生生态系统中，几丁质每年的产量就超过1×1011 t[3]。天然几丁质由于分子内和分子间强大的氢键作用而性质稳定，难溶于水、稀酸、碱，仅溶于浓盐酸、硝酸等强酸以及一些离子液体 (如咪唑基离子液体等)[4]，所以很难得到有效利用，造成了严重的资源浪费。然而，几丁质的水解产物N-乙酰壳寡糖却具有水溶性、细胞膜渗透性和易吸收性等理化特点；还具有抗菌、抗病毒、抗肿瘤、抗氧化、免疫调节、降血压和降低胆固醇等多种生物活性[5]，广泛应用于生物化工、生物医药、农产品、化妆品、保健品等领域[6]。其中，N, N-二乙酰基壳二糖 [(GlcNAc)2] 具有抗肿瘤、抗菌和免疫刺激作用等生物活性[7]。(GlcNAc)2在医疗方面，能够防止术后组织黏连[8]；在农业方面，能够促进植物生长和诱导植物先天免疫[9]。此外，(GlcNAc)2还能诱导几丁质酶在链霉菌、蘑菇等生物中的表达[10-11]。
降解几丁质的方法包括物理法、化学法和生物酶法。迄今为止，采用物理法制备N-乙酰壳寡糖的研究较少，在工业上尚未形成规模化生产。传统的化学法采用强酸处理几丁质，反应非常复杂，难以控制，并会产生各种二级化合物，难以去除。生物酶法作为绿色生产技术，可克服物理和化学法的缺点[12]，其反应条件温和、产率高、对环境无污染，因而受到越来越多的关注[13]。
几丁质酶 (Chitinase, EC 3.2.1.14) 属于糖苷水解酶，通过水解几丁质的β-(1→4) 糖苷键将其降解为N-乙酰壳寡糖或GlcNAc[14]。在碳水化合物-活性酶数据库 (CAZy) 中，根据催化结构域氨基酸序列的保守性，被分为糖苷水解酶18家族 (GH18)、19家族 (GH19) 和20家族 (GH20)。GH18家族几丁质酶主要由细菌、真菌、昆虫、植物和动物产生，GH19家族几丁质酶主要来源于高等植物和某些细菌[15]，而GH20家族几丁质酶来源于细菌和人类。不同来源的几丁质酶已在大肠杆菌 (Escherichia coli) 和枯草芽孢杆菌 (Bacillus subtilis) 等异源宿主系统中以重组蛋白的形式表达[16-17]。几丁质酶有许多潜在的应用，除了用于制备N-乙酰壳寡糖[18-19]外，还可以用作致病菌的抗真菌剂[20]，利用和清理几丁质废弃物并将其转化为生物乙醇[21]。
海洋环境中含有丰富的几丁质，是筛选产几丁质酶微生物的重要来源地之一。在前期研究中，本课题组从红树林虾场土壤中筛选到一株能够利用几丁质的放线菌橙黄小单孢菌(Micromonospora aurantiaca)。迄今为止，已有诸多放线菌来源的几丁质酶被分离和鉴定[22-23]，但是橙黄小单孢菌来源几丁酶的相关研究尚未见报道。本研究以橙黄小单孢菌的基因组DNA为模板，克隆到了一条几丁质酶基因Machi3，并成功在大肠杆菌中表达。之后，对重组酶进行纯化及酶学特性的研究，采用扫描电镜 (SEM) 进一步验证了该重组酶对几丁质底物的水解活性。最后，利用该重组酶水解胶体几丁质制备 (GlcNAc)2，为几丁质资源的开发和利用奠定了基础。
1. 材料与方法
1.1 菌株与质粒
橙黄小单孢菌：于海南省海口市东寨港红树林虾场土壤中筛选得到，保存于本实验室；pMD-19T载体购于宝生物工程 (大连) 有限公司；E. coli DH5α、E. coli BL21 (DE3) 和表达载体pET-28a购于Novagen公司。
1.2 主要材料与仪器
细菌基因组DNA提取试剂盒、质粒提取试剂盒 (上海捷瑞生物工程有限公司)；PCR产物纯化试剂盒、DNA胶回收试剂盒、蛋白Marker [生工生物工程 (上海) 股份有限公司]；T4 DNA连接酶、DNA Marker [宝生物工程 (大连) 有限公司]；限制性内切酶 (赛默飞世尔科技公司)；GlcNAc (上海源叶生物科技有限公司)；(GlcNAc)2-6标准品 (山东青岛博智汇力有限公司)；牛血清白蛋白 (武汉赛维尔生物科技有限公司)；其他试剂均为国产分析纯。
DYCZ-24DN型双板垂直电泳槽、DYPC-31DN型水平电泳槽 (北京六一生物科技有限公司)；A300型PCR扩增仪、LG2020D型凝胶成像系统 (杭州朗基科学仪器有限公司)；高效液相色谱柱：Shodex Asahipak NH2P-50 4E (4.6 mm×250 mm, 5 μm) (北京迪科马科技有限公司)；K6600B型全波长酶标仪 (北京凯奥科技发展有限公司)；NRY-200型全温培养摇床 (上海南荣实验室设备有限公司)；Phenom ProX型SEM (飞纳科学仪器有限公司)。
1.3 实验方法
1.3.1 目的基因的克隆和生物信息学分析
本研究根据美国国家生物信息中心 (NCBI) 已公布的橙黄小单孢菌ATCC 27029的全基因组DNA序列信息 (NCBI登陆号CP002162)，从中挖掘出一条可能编码几丁质酶的基因序列，根据该基因序列采用Primer Premier 5.0软件设计特异性引物Primer-F (CATATGGTGGTCCTCCCGAACAACTTC，酶切位点为Nde I) 和Primer-R (AAGCTTTCAGGCCGGCAGCCGC，酶切位点为Hind III)，委托生工生物工程 (上海) 股份有限公司合成。PCR反应体系 (50 µL)：TaKaRa LA Taq 0.5 µL，Primer-F (20 µmol·L−1) 1 µL，Primer-R (20 µmol·L−1) 1 µL，基因组DNA 1 µL，dNTP Mixture 8 μL，2 × GC Buffer I 25 μL，ddH2O 13.5 µL。PCR扩增条件：94 ℃预变性 1 min；94 ℃变性 30 s，58 ℃退火 30 s，72 ℃ 延伸100 s，30个循环；72 ℃ 延伸 7 min。PCR产物经1% (质量分数) 琼脂糖凝胶电泳验证后，根据DNA胶回收试剂盒说明书回收目的产物片段。
生物信息学分析软件及网站：氨基酸序列分析采用NCBI的BLAST功能，信号肽序列分析用SignalP 5.0在线服务器 (http://www.cbs.dtu.dk/services/SignalP/)，蛋白质同源序列比对采用DNAMAN软件，保守域的分析采用NCBI的CD-search工具 (https://www.ncbi.nlm.nih.gov/cdd/)，理论分子量和等电点的预测采用ExPASy的ProtParam工具 (http://web.expasy.org/protparam/)。
1.3.2 重组表达载体的构建
将目的基因连接至pMD19-T载体，之后转化至E. coli DH5α感受态细胞中。根据菌落PCR鉴定结果，挑选阳性克隆子送至北京擎科生物科技有限公司测定目的基因序列。将测序正确的基因序列和pET-28a表达载体分别用Nde I和Hind III双酶切，回收目的基因片段和酶切质粒，用T4 DNA连接酶连接，构建重组表达质粒，并转化至E. coli BL21 (DE3) 感受态细胞中。
1.3.3 重组几丁质酶的诱导表达和分离纯化
将含有目的基因的重组菌接种到含50 μg·mL−1卡那霉素的Luria-Bertani (LB) 液体培基中，在37 ℃、180 r·min−1下振荡培养至菌液OD600 nm为0.8时，加入终浓度为0.2 mmol·L−1的异丙基-β-D-硫代半乳糖苷 (IPTG)，在16 ℃、180 r·min−1下诱导培养16 h。诱导培养结束后，离心 (4 ℃, 6 000 r·min−1, 5 min) 收集菌体，用0.85% (质量分数) 的生理盐水清洗后用pH 7.0的磷酸盐缓冲液重悬，之后置于冰水浴中超声破碎 (功率350 W，工作 2 s，间隔5 s)。超声破碎结束后在4 ℃、6000 r·min−1下离心15 min，上清液即为粗酶液。将粗酶液上样于镍亲和层析色谱柱进行分离纯化，并利用SDS-PAGE分析蛋白的纯度和分子量。
1.3.4 胶体几丁质的制备
参照Pan等[24]的方法制备胶体几丁质，略有改动。称取10 g α-几丁质粉末，缓慢加入到200 mL预冷的浓盐酸中 (边加边搅拌)，在4 ℃冰箱中静置至溶液呈黄色。之后，向上述混合物中加入2 L预冷的50% (体积分数) 乙醇，快速搅拌后于4 ℃下放置过夜。然后在4 ℃、4 000 r·min−1下离心15 min，并收集沉淀部分。用大量去离子水反复冲洗沉淀部分至中性，离心 (10 000 r·min−1, 10 min) 后获得的沉淀物即为胶体几丁质。
1.3.5 重组几丁质酶的酶活力测定
以胶体几丁质作为底物，采用3,5-二硝基水杨酸 (DNS) 法测定还原糖产量[25]。用0~5 μmol·L−1的GlcNAc制作还原糖标准曲线，得到标准曲线方程y=0.122 6x−0.030 9，R2=0.990 8。酶活测定方法如下：将胶体几丁质加入至pH 7.0的Tris-HCl缓冲液中，混合均匀后制成终质量浓度为50 mg·mL−1的胶体几丁质混悬液。取0.5 mL上述胶体几丁质混悬液与0.5 mL酶液混合，于50 ℃水浴中磁力搅拌反应30 min。反应结束后离心 (10 000 r·min−1, 2 min)，获得上清液。取0.5 mL上清液，加入0.5 mL DNS溶液，摇匀后煮沸5 min进行显色反应，用酶标仪于540 nm处测定吸光值，根据还原糖标准曲线计算还原糖产量。同时，以加热灭活的酶液作为空白对照组。酶活力单位定义为：在上述反应条件下，每分钟催化胶体几丁质产生1 μmol GlcNAc所需要的酶量，用U表示。
1.3.6 重组几丁质酶的蛋白浓度测定
采用Bradford法[26]测定重组酶的蛋白浓度。以0~0.1 mg·mL−1的牛血清白蛋白制作蛋白浓度标准曲线，得到标准曲线方程y=4.579 1x−0.002，R2=0.990 7。根据吸光值和标准曲线方程计算待测样的蛋白浓度。
1.3.7 重组几丁质酶的酶学性质研究
1) 最适反应温度及热稳定性研究。取等量的酶液，在不同反应温度 (25~65 ℃) 下按照标准方法测定重组几丁质酶的酶活。将最适反应温度下的相对酶活定义为100%。为了测定重组几丁质酶的热稳定性，将酶液置于不同温度 (40~60 ℃) 下孵育1 h后取样，按标准方法，在最适反应温度下测定残余酶活。将孵育前几丁质酶的酶活定义为100%。
2) 最适反应pH及pH稳定性研究。取等量的酶液，在不同pH的缓冲液 (50 mmol·L−1，pH 3.0~5.0为柠檬酸-柠檬酸钠缓冲液；pH 5.0~8.0为磷酸盐缓冲液；pH 7.0~9.0为Tris-HCl缓冲液；pH 9.0~10.0为甘氨酸-NaOH缓冲液) 中按照标准方法测定重组几丁质酶的酶活。将最适反应pH下的相对酶活力定义为100%。为了测定重组几丁质酶的pH稳定性，将酶液加入到不同pH的缓冲液 (50 mmol·L−1，pH 3.0~4.0为柠檬酸-柠檬酸钠缓冲液；pH 5.0~6.0为磷酸盐缓冲液；pH 7.0~8.0为Tris-HCl缓冲液；pH 9.0~10.0为甘氨酸-NaOH缓冲液) 中，在4 ℃下孵育1 h。之后，在最适反应温度和pH条件下，按照标准方法测定残余酶活。孵育前重组几丁质酶的酶活被定义为100%。
3) 化学试剂和金属离子对酶活力的影响。为了研究化学试剂对重组几丁质酶酶活力的影响，向反应体系中分别加入不同的化学试剂，包括尿素 (终浓度200 mmol·L−1)，十二烷基硫酸钠 (SDS，终浓度10 mmol·L−1)，乙二胺四乙酸 (EDTA，终浓度10 mmol·L−1)，吐温 (Tween 20、40、60、80，终体积分数1%) 和曲拉通X-100 (终体积分数1%)，在最适反应温度和pH条件下，按照标准方法测定重组酶MaChi3的酶活力。为了研究金属离子对重组几丁质酶酶活力的影响，向反应体系中分别加入终浓度1 mmol·L−1的钾离子 (K+)、镁离子 (Mg2+)、亚铁离子 (Fe2+)、钙离子 (Ca2+)、钴离子 (Co2+)、铁离子 (Fe3+)、钡离子 (Ba2+)、锌离子 (Zn2+)、铜离子 (Cu2+) 和银离子 (Ag+)，在最适反应温度和pH条件下，按照标准方法测定重组酶MaChi3的酶活力。将未添加任何化学试剂和金属离子时的相对酶活定义为100%。
4) 底物特异性及动力学参数测定。将酶液分别与终质量浓度为2 mg·mL−1的α-几丁质、胶体几丁质、虾壳粉、淀粉、纤维素、不同脱乙酰度的壳聚糖 (脱乙酰度50%、70%、80%、85%、90%、95%) 进行反应，在最适反应温度和pH条件下，测定重组酶MaChi3的酶活力，研究该重组酶对不同底物的特异性。
在含有不同质量浓度 (1~8 mg·mL−1) 胶体几丁质的反应体系中，在最适反应温度和pH条件下测定重组酶MaChi3的酶活力。采用双倒数作图法，利用Origin 8.5软件进行线性拟合，根据拟合公式计算MaChi3的动力学参数米氏常数Km和最大反应速率Vmax值。
1.3.8 几丁质的SEM表征
将10 mg α-几丁质和胶体几丁质分别加入到50 mmol·L−1 pH 7.0的Tris-HCl缓冲液中，均匀分散后加入2 mL酶液，于55 ℃下振荡孵育24 h。孵育结束后，在10 000 r·min−1下离心10 min，收集底物样品，在60 ℃下干燥24 h。采用SEM对孵育前后的α-几丁质和胶体几丁质的表面微观形态进行观察。观察前用导电双面胶带将样品固定在样品台上，在真空环境下进行镀金处理，物料表面被一层薄薄的金覆盖形成导电膜以获得导电性。使用SEM的颗粒表面分析系统观察样品的表面形态，加速电压为10 kV，放大倍数为2500倍[27]。
1.3.9 N, N-二乙酰基壳二糖的制备
将10 mg胶体几丁质加入到50 mmol·L−1 pH 7.0的Tris-HCl缓冲液中，均匀分散后加入5 mL酶液，于55 ℃下振荡反应。定时从反应瓶中取样 (0、2、4、6、8、10、12、24 h)，在沸水浴中孵育10 min后离心 (10 000 r·min−1, 5 min)，取上清液用0.22 μm滤头过滤，滤液用于HPLC检测分析。检测条件：安捷伦1260高效液相色谱仪，配备Asahipak NH2P-504E (250 mm×4.6 mm，Shodex，日本) 色谱柱。流动相为75% (体积分数) 乙腈，流速为0.5 mL·min−1，柱温箱温度设定为25 ℃，检测器温度为25 ℃，进样量为10 μL，检测波长为210 nm。(GlcNAc)1-6的保留时间依次为13.2、18.0、25.3、34.1、49.7和77.3 min。
1.3.10 数据处理
采用Origin Pro 2018软件对数据进行统计分析，每组实验做3次平行。
2. 结果
2.1 几丁质酶基因的克隆及生物信息学分析
以本实验室保存的橙黄小单孢菌的基因组DNA为模板，利用特异性引物进行PCR扩增，获得一条大约1 400 bp的基因条带 (图1-a)。将测序结果提交至NCBI数据库进行BLAST分析，发现其与查找到的几丁质酶基因序列一致性达到100% (NCBI登录号WP_013285221)，表明已成功从M. aurantiaca的基因组DNA中克隆到几丁质酶基因，并将其记为Machi3。Machi3基因的全长为1 371 bp，编码456个氨基酸，预测的蛋白理论相对分子质量和pI分别为47.44 kD和7.55。信号肽预测结果表明，在MaChi3的N端含有信号肽，信号肽的切割位点在第33—第34个氨基酸之间。MaChi3含有一个Glyco_18结构域 (位于40—305位氨基酸残基) 和一个具有碳水化合物结合作用的RICIN结构域 (位于331—455位氨基酸残基，图1-b)。将MaChi3的氨基酸序列提交至NCBI数据库进行BLAST分析，发现其与来源于那拉提瓦小单孢菌 (M. narathiwatensis) (NCBI登录号WP_091195038)、邱氏疣孢菌 (Verrucosispora sioxanthis) (NCBI登录号WP_164447165) 和石渣球孢菌 (Allocatelliglobosispora scoriae) (NCBI登录号WP_184833644) 的几丁质酶序列一致性分别为87.3%、83.5%和78.3%。但是，上述几丁质酶均为未经验证的假定几丁质酶。因此，将MaChi3的氨基酸序列提交至PDB数据库，结果表明该酶与来源于约氏黄杆菌 (Flavobacterium johnsoniae) UW101的GH18家族几丁质酶具有32%的一致性。GH18家族几丁质酶具有保守序列DxDxE和SxGG，DxDxE基序为几丁质酶的催化中心，而SxGG基序负责酶与底物的结合。在催化过程中，SxGG基序中的丝氨酸对DxDxE基序中第一个天冬氨酸上产生的暂时过剩的负电荷起稳定作用[28-29]。如图1-c所示，在MaChi3的催化域中含有符合上述保守基序的145DMDWE149序列，但是SxGG基序中的丝氨酸残基被丙氨酸残基 (103AVGG106) 所代替，这一发现与Vaaje-Kolstad等[30]报道的来源于乳酸乳球菌乳酸亚种 (Lactococcus lactis ssp. Lactis)的LlChi18A相一致。事实上，除了SxGG基序中的丝氨酸外，酪氨酸也具有稳定电荷的作用。尽管MaChi3缺乏丝氨酸残基，但却含有酪氨酸残基 (Y45，相当于LlChi18A的Y48)。另外，Vaaje-Kolstad等[30]研究发现，在与LlChi18A催化域相似性较高的50个不同来源的GH18家族几丁质酶的催化域中，约有一半蛋白质的丝氨酸或酪氨酸残基被替代，但是这两个保守氨基酸残基不会被同时替代。总之，以上序列分析结果表明，MaChi3是一种新型的GH18家族几丁质酶。
![]()
图 1 MaChi3 基因 PCR 扩增产物 (a)、保守域分析 (b) 和部分催化域氨基酸多序列比对 (c)注：(a) M为DNA分子量标准；(c)图中列出的几丁质酶氨基酸序列来源于B. circulans WL-12 (PDB ID : 1ITX)，C. shinanonensis (PDB ID : 6KXL)，B. thuringiensis (PDB ID : 6BT9)，Paenibacillus sp. FPU-7 (PDB ID : 5GZU) 和L. lactis ssp. Lactis (NCBI登录号AAK06048)。保守基序SxGG和DxDxE在底部用带下划线的红色字母标识，保守的酪氨酸残基在上部用※标识。Figure 1. PCR products (a), conserved domain analysis (b) and sequence alignment of partial predicted catalytic domain (c) of MaChi3Note: (a) M is DNA molecular mass standard. (c) The listed sequences include the chitinases from B. circulans WL-12 (PDB ID: 1ITX), C. shinanonensis (PDB ID: 6KXL), B. thuringiensis (PDB ID: 6BT9), Paenibacillus sp. FPU-7 (PDB ID: 5GZU) and L. lactis ssp. Lactis (NCBI accession NO.: AAK06048). The conserved sequence motifs SxGG and DxDxE are showed at the bottom with red color underlined, and the conserved tyrosine residue is indicated by ※.2.2 MaChi3的诱导表达
重组菌经诱导表达后，利用镍亲和层析色谱柱对粗酶液进行纯化，在大约47 kD处出现一条蛋白条带 (图2)，与推测的蛋白相对分子质量 (47.44 kD) 基本吻合，表明重组酶MaChi3已成功在大肠杆菌中表达。
![]()
图 2 重组酶 MaChi3 的 SDS-PAGE 分析Figure 2. SDS-PAGE analysis of purified recombinant chitinase MaChi32.3 MaChi3的酶学性质研究
2.3.1 MaChi3的最适反应温度及热稳定性
为了测定重组几丁质酶MaChi3的最适反应温度，以胶体几丁质作为底物，在25~65 ℃下进行反应。结果显示，MaChi3的最适反应温度为55 ℃，并且在50~60 ℃范围内具有较高的相对酶活力 (>90.4%)。当反应温度提高至65 ℃时，相对酶活力急剧下降，仅为23.5% (图3-a)。为了测定MaChi3的温度稳定性，将其置于40~60 ℃水浴中孵育1 h，然后在最适反应温度下测定其残余酶活力。由图3-b可知，MaChi3的温度稳定性良好，在50 ℃下孵育1 h后残余酶活力超过70%，在其最适反应温度55 ℃下孵育1 h后残余酶活力保持在67.3%。
![]()
图 3 重组几丁质酶MaChi3的最适反应温度 (a) 及热稳定性 (b)Figure 3. Optimal temperature (a) and thermal stability (b) of recombinant chitinase MaChi32.3.2 MaChi3的最适反应pH及pH稳定性
重组几丁质酶MaChi3的最适反应pH为7.0，并且在pH 6.0~9.0范围内表现出较好的酶活力 ，相对酶活力高于64% (图4-a)。为了测定MaChi3的pH稳定性，将该酶置于pH 3.0~10.0的缓冲液中孵育，之后在最适反应温度和pH条件下测定其残余酶活力。MaChi3在pH 6.0~9.0范围内孵育1 h后残余酶活力保持在80%以上，展现出了较好的稳定性 (图4-b)。此外，在其最适反应pH (pH 7.0) 下孵育1 h后对其酶活力基本无影响，残余酶活力为98%。
![]()
图 4 重组几丁质酶 MaChi3 的最适反应pH (a) 及pH稳定性 (b)Figure 4. Optimal pH (a) and pH stability (b) of recombinant chitinase MaChi32.3.3 化学试剂和金属离子对MaChi3的影响
如表1所示，Triton X-100和EDTA对MaChi3的酶活力无显著影响。尿素、SDS、吐温 20和吐温60对MaChi3具有不同程度的抑制作用，相对酶活力为58.1%~90.4%。吐温40和吐温80对MaChi3的酶活力有轻微的促进作用，相对酶活力分别为104.1%和108.4%。不同的金属离子对MaChi3表现为不同程度的激活或抑制作用。Mg2+和Ca2+对其酶活力具有激活作用，相对酶活力分别为117.7%和107.3%。Zn2+和Ba2+对MaChi3的酶活力几乎无影响。而K+、Ag+、Co2+、Fe2+、Cu2+和Fe3+对MaChi3的酶活力呈现出不同程度的抑制作用，相对酶活力降低至43.3%~87.2% (表1)。
表 1 化学试剂和金属离子对重组酶 MaChi3 酶活力的影响Table 1. Effects of metal ions and chemical reagents on activity of recombinant chitinase MaChi3化学试剂
Chemical reagent终浓度/体积分数
Final concentration/
Volume fraction相对酶活力
Relative activity/
%对照 Control 100 尿素 Urea 200 mmol·L−1 70.9±0.8 乙二胺四乙酸 EDTA 10 mmol·L−1 94.2±1.6 十二烷基硫酸钠 SDS 10 mmol·L−1 58.1±0. 5 吐温20 Tween 20 1% 90.4±0.8 吐温40 Tween 40 1% 104.1±1.6 吐温60 Tween 60 1% 79.6±1.0 吐温80 Tween 80 1% 108.4±0.8 曲拉通 X-100 Triton X-100 1% 101.0±0.9 钾离子 K+ 1 mmol·L−1 86.4±0.7 银离子 Ag+ 1 mmol·L−1 43.3±0.4 镁离子 Mg2+ 1 mmol·L−1 117.7±0.7 亚铁离子 Fe2+ 1 mmol·L−1 65.4±0.7 钙离子 Ca2+ 1 mmol·L−1 107.3±0.4 钴离子 Co2+ 1 mmol·L−1 87.2±0.7 钡离子 Ba2+ 1 mmol·L−1 99.5±1.3 锌离子 Zn2+ 1 mmol·L−1 97.6±1.1 铜离子 Cu2+ 1 mmol·L−1 45.0±0.7 铁离子 Fe3+ 1 mmol·L−1 59.3±1.1 2.3.4 MaChi3的底物特异性和动力学参数
分别以α-几丁质、胶体几丁质、淀粉、纤维素、虾壳粉和不同脱乙酰度的壳聚糖作为底物，研究重组几丁质酶MaChi3的底物特异性。如图5所示，MaChi3在催化胶体几丁质水解时表现出了最高的酶活力 (2.24 U·mg−1)。但是，以α-几丁质粉末作为底物时，酶活力仅为0.59 U·mg−1，前者约是后者的3.8倍。当MaChi3以虾壳粉为底物时，酶活力为0.84 U·mg−1，与其催化α-几丁质水解的活性相当。以50%脱乙酰度的壳聚糖为底物时，MaChi3展现了较高的水解能力 (1.66 U·mg−1)。但是，随着脱乙酰度的增加水解效率逐渐降低。以淀粉和纤维素为底物时，酶活力均较低，分别为0.40和0.39 U·mg−1。此外，以胶体几丁质作为底物，MaChi3的Km和Vmax值分别为5.93 mg·mL−1和8.58 μmol·(min·mg)−1。
2.3.5 SEM分析
为了进一步验证重组几丁质酶MaChi3对几丁质底物的水解能力，对酶解前及酶解后的α-几丁质和胶体几丁质进行SEM观察。α-几丁质的纤维结构紧密、表面略光滑 (图6-a)。经MaChi3酶解后，表面微观结构被破坏，变得粗糙、不平整，并出现大小不均匀的碎片附着在表面 (图6-b)。α-几丁质被制成胶体几丁质后，纤维结构变得松散 (图6-c)，暴露出更多的酶切位点，因而经MaChi3酶解后，表面变得更加粗糙和褶皱，有更多的小碎片附着在表面 (图6-d)。以上结果表明，重组几丁质酶MaChi3对几丁质具有水解活力，并对胶体几丁质的水解活力更高，该结论与底物特异性结果 (图5) 一致。
![]()
图 6 α-几丁质 (a) 和胶体几丁质 (b) 以及MaChi3酶解后的α-几丁质 (c) 和胶体几丁质 (d) SEM图Figure 6. SEM images of α-chitin (a), colloidal chitin (b), MaChi3-hydrolyzed α-chitin (c) and MaChi3-hydrolyzed colloidal chitin (d)2.3.6 N, N-二乙酰基壳二糖的制备
如图7-a所示，在反应初始阶段，重组几丁质酶MaChi3快速地将胶体几丁质水解为 (GlcNAc)2，同时生成了少量的GlcNAc。反应24 h后，(GlcNAc)2的产率基本保持稳定，为285.5 mg·g−1几丁质，GlcNAc的产率为39.8 mg·g−1几丁质，考虑到MaChi3的热稳定性 (图3-b)，此时几丁质酶可能已经失去了活性。
![]()
图 7 胶体几丁质水解产物分析 (a) 和 (GlcNAc)1-6 标准品及水解产物的HPLC图 (b)Figure 7. Hydrolysis products of colloidal chitin (a) and HPLC spectra of (GlcNAc)1-6 stands and hydrolysis products (b)3. 讨论
虾、蟹等甲壳类动物在生产加工过程中会生成大量的废弃物，其中，外壳占废物干质量的20%~58%。这些外壳中富含蛋白质 (20%~40%)、矿物质 (30%~60%) 和几丁质 (20%~40%) 等成分。在工业上，虾、蟹的外壳是生产几丁质的主要原料。几丁质经水解后得到的N-乙酰壳寡糖和单糖具有多种生物活性，已被应用于多个领域。因此，从甲壳类动物的壳中提取几丁质并将其转化为更具有应用价值的水解产物 [ 如(GlcNAc)2] 是缓解废弃物污染问题的有效途径。
在原核生物中，放线菌是最好的几丁质分解者之一，其体内含有丰富的几丁质降解酶系[31]，可以利用几丁质作为唯一的碳源和氮源。本研究从放线菌橙黄小单孢菌的基因组DNA中克隆到了一条几丁质酶基因Machi3，经过生物信息学分析得知，MaChi3属于GH18家族几丁质酶。为了简化酶的分离纯化过程，并提高酶的产量[32]，本研究将几丁质酶基因Machi3在E. coli BL21 (DE3) 中重组表达，经过纯化获得了电泳纯的重组几丁质酶。
不同微生物来源的几丁质酶对温度表现出不同的敏感性，如来源于橙色假交替单胞菌(Pseudoalteromonas aurantia) 的几丁质酶最适反应温度为60 ℃[33]，而来源于冷喜几丁质节杆菌(Chitiniphilus shinanonensis) 的几丁质酶最适反应温度为50 ℃[34]。Huo等[35]克隆表达了来源于维氏气单胞菌 (Aeromonas veronii) B565的几丁质酶，该重组酶的最适反应温度为40 ℃，在50 ℃下孵育1 h后残余酶活力仅为20%。相比之下，本研究克隆表达的重组几丁质酶MaChi3具有较高的最适反应温度和良好的热稳定性，为拓宽几丁质酶的应用奠定了良好的基础。同时，不同来源的几丁质酶对酸碱适应性差异较大，已报道的几丁质酶最佳pH值范围为3.0~10.0[36-37]。由于几丁质酶在不同的氢离子(H+)浓度下解离程度不同，其高级结构中S-S键、疏水键及氢键受到影响，从而影响几丁质酶构象的变化而表现出不同的活性。Wang等[38]克隆、表达了假交替单胞菌 (Pseudoalteromonas sp.) DL-6来源的几丁质酶ChiA，其最适反应pH为9.0，在碱性条件下 (pH 8.0~12.0) 较稳定。Chen等[39]研究发现铜绿假单胞菌 (Pseudomonas aeruginosa) 几丁质酶的最适pH为4.5，在酸性条件下 (pH 4.5~6.0) 较稳定。而本研究的重组几丁质酶MaChi3是一种中性酶，并且在中性pH附近 (pH 6.0~9.0) 具有较好的稳定性。由表1可知，金属离子螯合剂EDTA对MaChi3的酶活性无显著影响，说明MaChi3为非金属依赖性蛋白酶。除了SDS外，MaChi3对所测试的化学试剂表现出了较强的抵抗力。据报道，SDS能够破坏蛋白质分子间和分子内的氢键，从而破坏蛋白质的二级和三级结构[40]。MaChi3与大多数几丁质酶相似，表现出严格的底物特异性[41-42]，对胶体几丁质的酶活力高于其他底物，这可能与底物的物理形态和脱乙酰度有关。此外，本研究对MaChi3降解胶体几丁质的产物进行了分析。在反应过程中 (GlcNAc)2为主要产物，GlcNAc为唯一副产物，产物类型简单，有利于后续的分离纯化。(GlcNAc)2已被广泛用作药用生物活性物质以及抗菌剂和螯合剂。因此，MaChi3在制备 (GlcNAc)2 方面具有潜在应用价值。
4. 结论
本研究从海洋微生物橙黄小单孢菌的基因组DNA中克隆到了一条几丁质酶基因Machi3，并成功在大肠杆菌中表达。MaChi3的相对分子质量为47 kD，为新型的GH18家族几丁质酶。MaChi3的最适反应温度为55 ℃，最适反应pH为7.0。在体外的稳定性较好，在55 ℃下孵育1 h后仍能保持大部分的酶活力 (67%)，在pH 6.0~9.0范围内孵育1 h后残余酶活力超过80%。吐温40、吐温80、Ca2+和Mg2+对MaChi3的酶活力有轻微的促进作用，Triton X-100、EDTA、Zn2+和Ba2+对MaChi3无显著性影响，其他化学试剂和金属离子对MaChi3表现为不同程度的抑制作用。该酶对α-几丁质、胶体几丁质、虾壳粉、不同脱乙酰度 (50%~95%) 的壳聚糖、淀粉以及纤维素均有一定的水解活力，其中胶体几丁质是MaChi3的最适底物，酶活力为2.24 U·mg−1。α-几丁质被制备成胶体几丁质后，纤维结构变的更加松散，因而更有利于MaChi3的水解作用。以胶体几丁质为底物时，MaChi3的Km和Vmax值分别为5.93 mg·mL−1和8.58 μmol·(min·mg)−1。此外，胶体几丁质经MaChi3水解后主要产物为 (GlcNAc)2，反应24 h后产率为285.5 mg·g−1几丁质，同时生成了少量的GlcNAc。MaChi3优良的酶学特性使其在 (GlcNAc)2制备方面有一定的应用前景。
-
![]()
图 1 2016—2020 年中西太平洋金枪鱼围网捕捞自由鱼群的作业位置分布
Figure 1. Distribution of operating positions for tuna purse seine fishing of free fish schools in Western and Central Pacific Ocean from 2016 to 2020
![]()
图 2 2016—2020 年模型性能的综合评价
Figure 2. Comprehensive evaluation of model performance from 2016 to 2020
![]()
图 3 2016—2020 年使用 Jackknife 检验得到的测试增益
Figure 3. Test gain obtained using Jackknife test from 2016 to 2020
![]()
图 4 海表温度对中西太平洋鲣栖息地适宜性的响应曲线
Figure 4. Response curves of SST to habitat suitability of skipjack tuna in Western and Central Pacific Ocean
![]()
图 5 200 m 水层溶解氧浓度对中西太平洋鲣栖息地适宜性的响应曲线
Figure 5. Response curves of DO200 to habitat suitability of skipjack tuna in Western and Central Pacific Ocean
![]()
图 6 2016—2020 年中西太平洋鲣自由鱼群潜在栖息地适宜性指数分布
注：图中白色区域表示无数据。
Figure 6. Distribution of potential habitat suitability index for free swimming school of skipjack tuna in Western and Central Pacific Ocean from 2016 to 2020
Note: The white area in the figure indicates no data.
表 1 2016—2020 年环境因子年平均贡献率
Table 1 Average annual contribution rate of environmental factors from 2016 to 2020
环境因子
Environmental factor年平均贡献率 Average annual contribution rate/% 2016 2017 2018 2019 2020 均值 Mean 海表温度 SST 59.35 61.11 68.61 63.15 55.47 61.54 50 m 叶绿素 a 浓度 CHL50 11.41 10.59 10.94 12.07 13.86 11.77 200 m 溶解氧浓度 DO200 6.14 9.90 7.12 5.50 5.38 6.81 200 m 水温 Temp200 4.73 3.79 3.80 5.45 7.44 5.04 表层叶绿素 a 浓度 CHL0 2.90 3.81 1.18 4.25 6.48 3.72 海表盐度 SSS 4.77 3.75 3.06 2.22 3.61 3.48 海流速度 EV 2.68 1.13 1.01 1.03 1.69 1.51 表层溶解氧浓度 DO0 1.77 1.84 1.13 1.09 1.63 1.49 混合层深度 MLD 1.77 1.18 1.06 2.35 0.87 1.45 200 m 叶绿素 a 浓度 CHL200 1.24 1.79 0.92 1.41 1.69 1.41 100 m 叶绿素 a 浓度 CHL100 2.09 0.39 0.62 0.72 0.60 0.88 50 m 溶解氧浓度 DO50 0.84 0.39 0.43 0.46 0.96 0.62 南北向的海流速度 NV 0.32 0.33 0.11 0.32 0.34 0.28 
下载: 导出CSV
-
[1] Food and Agriculture Organization of the United Nations. The State of World Fisheries and Aquaculture, 2021[R]. Rome: FAO, 2021: 127-148.
[2] TRYGONIS V, GEORGAKARAKOS S, DAGORN L, et al. Spatiotemporal distribution of fish schools around drifting fish aggregating devices[J]. Fish Res, 2016, 177: 39-49. doi: 10.1016/j.fishres.2016.01.013
[3] HALLIER J P, GAERTNER D. Drifting fish aggregation devices could act as an ecological trap for tropical tuna species[J]. Mar Ecol Prog Ser, 2008, 353: 255-264. doi: 10.3354/meps07180
[4] 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.
[5] MUGO R, SAITOH S E I I, NIHIRA A, et al. Habitat characteristics of skipjack tuna (Katsuwonus pelamis) in the western North Pacific: a remote sensing perspective[J]. Fish Oceanogr, 2010, 19(5): 382-396. doi: 10.1111/j.1365-2419.2010.00552.x
[6] DUERI S, BOPP L, MAURY O. Projecting the impacts of climate change on skipjack tuna abundance and spatial distribution[J]. Glob Chang Biol, 2014, 20(3): 742-753. doi: 10.1111/gcb.12460
[7] ADAM M S, SIBERT J R. Population dynamics and movements of skipjack tuna (Katsuwonus pelamis) in the Maldivian fishery: analysis of tagging data from an advection-diffusion-reaction model[J]. Aquat Living Resour, 2002, 15(1): 13-23. doi: 10.1016/S0990-7440(02)01155-5
[8] SALAZAR J E, BENAVIDES I F, CABRERA C V P, et al. Generalized additive models with delayed effects and spatial autocorrelation patterns to improve the spatiotemporal prediction of the skipjack (Katsuwonus pelamis) distribution in the Colombian Pacific Ocean[J]. Reg Stud Mar Sci, 2021, 45: 101-115.
[9] MUGO R, SAITOH S I. Ensemble modelling of skipjack tuna (Katsuwonus pelamis) habitats in the western north Pacific using satellite remotely sensed data: a comparative analysis using machine-learning models[J]. Remote Sens, 2020, 12(16): 2591-2605. doi: 10.3390/rs12162591
[10] ALABIA I D, SAITOH S I, MUGO R, et al. Seasonal potential fishing ground prediction of neon flying squid (Ommastrephes pastrami) in the western and central North Pacific[J]. Fish Oceanogr, 2015, 24(2): 190-203. doi: 10.1111/fog.12102
[11] CHEN B Y, HONG Z, HAO X Q, et al. Environmental models for predicting habitat of the Indo-Pacific humpback dolphins in Fujian, China[J]. Aquat Conserv, 2020, 30(4): 787-793. doi: 10.1002/aqc.3279
[12] YOSHINAGA D H, FRANK H A. Histamine-producing bacteria in decomposing skipjack tuna (Katsuwonus pelamis)[J]. Appl Environ Microb, 1982, 44(2): 447-452. doi: 10.1128/aem.44.2.447-452.1982
[13] DAI L, WANG X, STAPLES K W, et al. Factors influencing successful fishing of tuna free-swimming schools in the equatorial western Pacific Ocean[J]. Turk J Fish Aquat Sc, 2019, 20(5): 341-350.
[14] SIBERT J, SENINA I, LEHODEY P, et al. Shifting from marine reserves to maritime zoning for conservation of Pacific bigeye tuna (Thunnus obesus)[J]. P Natl Acad Sci USA, 2012, 109(44): 18221-18225. doi: 10.1073/pnas.1209468109
[15] PRINCE E D, GOODYEAR C P. Hypoxia-based habitat compression of tropical pelagic fishes[J]. Fish Oceanogr, 2006, 15(6): 451-464. doi: 10.1111/j.1365-2419.2005.00393.x
[16] LEHODEY P, BERTIGNAC M, HAMPTON J, et al. El Niño Southern Oscillation and tuna in the western Pacific[J]. Nature, 1997, 389(6652): 715-718. doi: 10.1038/39575
[17] DI Y U, CHANG F C, BIN W, et al. Characterization of acid-and pepsin-soluble collagens from spines and skulls of skipjack tuna (Katsuwonus pelamis)[J]. Chin J Nat Medicines, 2014, 12(9): 712-720. doi: 10.1016/S1875-5364(14)60110-2
[18] ROGER C. The plankton of the tropical western Indian Ocean as a biomass indirectly supporting surface tunas (yellowfin, Thunnus albacares and skipjack, Katsuwonus pelamis)[J]. Environ Biol Fish, 1994, 39: 161-172. doi: 10.1007/BF00004934
[19] LEROY B, ITANO D G, USU T, et al. Vertical behavior and the observation of FAD effects on tropical tuna in the warm-pool of the western Pacific Ocean[M]//NIELSEN J L, ARRIZABALAGA H, FRAGOSO N, et al. Tagging and tracking of marine animals with electronic devices. Reviews: methods and technologies in fish biology and fisheries, vol 9. Dordrecht: Springer, 2009: 161-179.
[20] KIYOFUJI H, AOKI Y, KINOSHITA J, et al. Northward migration dynamics of skipjack tuna (Katsuwonus pelamis) associated with the lower thermal limit in the western Pacific Ocean[J]. Prog Oceanogr, 2019, 175: 55-67. doi: 10.1016/j.pocean.2019.03.006
[21] SCHAEFER K M. Vertical movement pattern of skipjack tuna (Katsuwonus pelamis) in the eastern equatorial Pacific Ocean, as revealed with archival tags[J]. Fish Bull, 2007, 105: 379-389.
[22] 杨胜龙, 周甦芳, 周为峰, 等. 基于Argo数据的中西太平洋鲣渔获量与水温、表层盐度关系的初步研究[J]. 大连水产学院学报, 2010, 25(1): 34-40. [23] GARCIA C B, GARCIA J, LOPEZ MARTIN M M, et al. Collinearity: revisiting the variance inflation factor in ridge regression[J]. J Appl Stat, 2015, 42(3): 648-661. doi: 10.1080/02664763.2014.980789
[24] IGARASHI H, SAITOH S I, ISHIKAWA Y, et al. Identifying potential habitat distribution of the neon flying squid (Ommastrephes bartramii) off the eastern coast of Japan in winter[J]. Fish Oceanogr, 2018, 27(1): 16-27. doi: 10.1111/fog.12230
[25] 赵静, 柳晓雪, 吴建辉, 等. 零膨胀模型在珍稀鱼类资源时空分布预测中的应用: 以长江口刀鲚为例[J]. 生态学杂志, 2020, 39(9): 3155-3163. doi: 10.13292/j.1000-4890.202009.028 [26] Sagarese S R, Frisk M G, Cerrato R M, et al. Application of generalized additive models to examine ontogenetic and seasonal distributions of spiny dogfish (Squalus acanthias) in the Northeast (US) shelf large marine ecosystem[J]. Can J Fish Aquat Sci, 2014, 71(6): 847-877. doi: 10.1139/cjfas-2013-0342
[27] 马金, 黄金玲, 陈锦辉, 等. 基于GAM的长江口鱼类资源时空分布及影响因素[J]. 水产学报, 2020, 44(6): 936-946. [28] PHILLIPS S J, ANDERSON R P, SCHAPIRE R E. Maximum entropy modeling of species geographic distributions[J]. Ecol Model, 2006, 190(3/4): 231-259.
[29] MANEL S, WILLIAMS H C, ORMEROD S J. Evaluating presence absence models in ecology: the need to account for prevalence[J]. J Appl Ecol, 2001, 38(5): 921-931. doi: 10.1046/j.1365-2664.2001.00647.x
[30] ZU T P, KANG R, WEN M L, et al. Belief reliability distribution based on maximum entropy principle[J]. IEEE Access, 2017, 6: 1577-1582.
[31] ALLOUCHE O, TSOAR A, KADMON R. Assessing the accuracy of species distribution models: prevalence, kappa and the true skill statistic (TSS)[J]. J Appl Ecol, 2006, 43(6): 1223-1232. doi: 10.1111/j.1365-2664.2006.01214.x
[32] MEYSMAN F J R, BRUERS S. Ecosystem functioning and maximum entropy production: a quantitative test of hypotheses[J]. Philosophical Transactions of the Royal Society B:Biol Sci, 2010, 365(1545): 1405-1416. doi: 10.1098/rstb.2009.0300
[33] SHCHEGLOVITOVA M, ANDERSON R P. Estimating optimal complexity for ecological niche models: a jackknife approach for species with small sample sizes[J]. Ecol Model, 2013, 269: 9-17. doi: 10.1016/j.ecolmodel.2013.08.011
[34] ELITH J H, GRAHAM C P, ANDERSON R, et al. Novel methods improve prediction of species distributions from occurrence data[J]. Ecography, 2006, 29(2): 129-151. doi: 10.1111/j.2006.0906-7590.04596.x
[35] ZU T, KANG R, WEN M, et al. Belief reliability distribution based on maximum entropy principle[J]. IEEE Access, 2017, 6: 1577-1582.
[36] WEST A M, KUMAR S, BROWN C S, et al. Field validation of an invasive species MaxEnt model[J]. Ecol Inform, 2016, 36: 126-134. doi: 10.1016/j.ecoinf.2016.11.001
[37] KUMAR P. Assessment of impact of climate change on Rhododendrons in Sikkim Himalayas using Maxent modelling: limitations and challenges[J]. Biodivers Conserv, 2012, 21(5): 1251-1266. doi: 10.1007/s10531-012-0279-1
[38] MORENO R, ZAMORA R, MOLINA J R, et al. Predictive modeling of microhabitats for endemic birds in South Chilean temperate forests using Maximum Entropy (MaxEnt)[J]. Ecol Inform, 2011, 6(6): 364-370. doi: 10.1016/j.ecoinf.2011.07.003
[39] LIU S J, YANG J. Modeling spatial patterns of forest fire in Heilongjiang Province using Generalized Linear Model and Maximum Entropy Model[J]. Chin J Ecol, 2013, 32(6): 1620-1628.
[40] 梁阳阳, 陈康, 崔凯, 等. 气候变化情景下须鳗鰕虎鱼在中国的潜在地理分布[J]. 大连海洋大学学报, 2022, 37(5): 739-746. [41] NURSAN M, YONVITNER Y, AGUS S B. Distribution of skipjack (Katsuwonus pelamis) fishing areas using purse seine fishing equipment in WPP 573[J]. J Trop Fish Manag, 2022, 6(1): 11-20.
[42] ZHANG K L, YAO L J, MENG J S, et al. Maxent modeling for predicting the potential geographical distribution of two peony species under climate change[J]. Sci Total Environ, 2018, 634: 1326-1334. doi: 10.1016/j.scitotenv.2018.04.112
[43] 陈世泳. 中西太平洋正鲣渔场推移与表明水温变化之关系[D]. 基隆: 台湾海洋大学, 2006: 24-26. [44] ZHOU C, HU Y Y, CAO J, et al. Comparison of nominal and standardized catch per unit effort data in quantifying habitat suitability of skipjack tuna in the equatorial Pacific Ocean[J]. Acta Oceanol Sin, 2022, 41(3): 1-10. doi: 10.1007/s13131-021-1922-z
[45] PICAUT J, IOUALALEN M, MENKÈS C, et al. Mechanism of the zonal displacements of the Pacific warm pool: implications for ENSO[J]. Science, 1996, 274(5292): 1486-1489. doi: 10.1126/science.274.5292.1486
[46] ELY B, VIÑAS J, ALVARADO BREMER J R, et al. Consequences of the historical demography on the global population structure of two highly migratory cosmopolitan marine fishes: the yellowfin tuna (Thunnus albacares) and the skipjack tuna (Katsuwonus pelamis)[J]. BMC Evol Biol, 2005, 5(1): 1-9. doi: 10.1186/1471-2148-5-1
[47] LEHODEY P. The pelagic ecosystem of the tropical Pacific Ocean: dynamic spatial modelling and biological consequences of ENSO[J]. Prog Oceanogr, 2001, 49(1/2/3/4): 439-468.
[48] BOYCE D G, TITTENSOR D P, WORM B. Effects of temperature on global patterns of tuna and billfish richness[J]. Mar Ecol Prog Ser, 2008, 355: 267-276. doi: 10.3354/meps07237
[49] MATHIEU-COSTELLO O, BRILL R W, HOCHACHKA P W. Structural basis for oxygen delivery: muscle capillaries and manifolds in tuna red muscle[J]. Comp Biochem Physiol A, 1996, 113(1): 25-31. doi: 10.1016/0300-9629(95)02059-4
[50] BUSHNELL P G, BRILL R W. Responses of swimming skipjack (Katsuwonus pelamis) and yellowfin (Thunnus albacares) tunas to acute hypoxia, and a model of their cardiorespiratory function[J]. Physiol Zool, 1991, 64(3): 787-811. doi: 10.1086/physzool.64.3.30158207
[51] SCHAEFER K M, FULLER D W. Vertical movement pattern of skipjack tuna (Katsuwonus pelamis) in the eastern equatorial Pacific Ocean, as revealed with archival tags[J]. Fish Bull, 2007, 105: 379-989.
[52] ZAINUDDIN M, KIYOFUJI H, SAITOH K, et al. Using multi-sensor satellite remote sensing and catch data to detect ocean hot spots for albacore (Thunnus alalunga) in the northwestern North Pacific[J]. Deep-Sea Res II, 2006, 53(3/4): 419-431. doi: 10.1016/j.dsr2.2006.01.007
[53] BERNAL D, BRILL R W, DICKSON K A, et al. Sharing the water column: physiological mechanisms underlying species-specific habitat use in tunas[J]. Rev Fish Biol Fisher, 2017, 27(4): 843-880. doi: 10.1007/s11160-017-9497-7
[54] DEARY A L, MORET-FERGUSON S, ENGELS M, et al. Influence of central Pacific oceanographic conditions on the potential vertical habitat of four tropical tuna species1[J]. Pac Sci, 2015, 69(4): 461-475. doi: 10.2984/69.4.3
[55] ANDRADE H A. The relationship between the skipjack tuna (Katsuwonus pelamis) fishery and seasonal temperature variability in the south western Atlantic[J]. Fish Oceanogr, 2003, 12(1): 10-18. doi: 10.1046/j.1365-2419.2003.00220.x
[56] LLOPIZ J K, HOBDAY A J. A global comparative analysis of the feeding dynamics and environmental conditions of larval tunas, mackerels, and billfishes[J]. Deep-Sea Res II, 2015, 113: 113-124. doi: 10.1016/j.dsr2.2014.05.014
-
期刊类型引用(1)
1. 吴海星，黄安妮，高霞，申铉日，夏光华，张雪莹. Lysinibacillus sphaericus胶原酶的异源表达、鉴定及其在抗氧化活性肽制备中的应用. 食品工业科技. 2025(06): 206-216 . 
百度学术
其他类型引用(0)
计量
- 文章访问数:
- HTML全文浏览量:
- PDF下载量:
- 被引次数: 1
粤公网安备 44010502001741号
