Study of suitable habitats for Sepiella maindroni in Zhoushan sea areas based on MaxEnt model
-
摘要:
曼氏无针乌贼 (Sepiella maindroni) 具有较高的经济和营养价值，近年因过度捕捞、生境破坏等导致其资源衰退。为了定量评估多元生态因子对曼氏无针乌贼分布的影响，以曼氏无针乌贼的重要产卵地和栖息地—舟山海域为研究区域，基于2017—2021年浙江近海的底拖网调查数据、2021—2022年自主航次的调查数据和全球生物多样性信息机构 (Global Biodiversity Information Facility, GBIF) 数据，利用最大熵 (Maximum Entropy, MaxEnt) 模型分析曼氏无针乌贼的生境适宜性，探索其修复潜力和适生区的空间分布格局。结果表明，海表温度最大值、海表盐度最小值和pH是影响曼氏无针乌贼分布的主要生态因子。曼氏无针乌贼偏向于在海洋表面温度较高、盐度较大的海域活动，且能接受弱碱性环境。当舟山海域内海表温度最大值介于25.5~28.5 ℃、海表盐度最小值介于20‰~30‰、pH介于7.8~8.3时，曼氏无针乌贼存在的概率较大。曼氏无针乌贼在舟山海域的适生区面积分布较广 (30 213.40 km2)，占研究区域总面积的75.11%。其中高适生区总体分布在中街山列岛、嵊泗列岛和马鞍列岛，面积为4 600.78 km2。次高适生区、中适生区和低适生区的面积分别为8 205.52、9 723.31和7 683.78 km2。从修复潜力上看，曼氏无针乌贼在舟山海域有较大的恢复潜力。
Abstract:Cuttlefish (Sepiella maindroni) has high economic and nutritional value, but its resources have declined in recent years due to overfishing and habitat destruction. To quantitatively assess the impact of multiple ecological factors on the distribution of S. maindroni, we selected the Zhoushan sea area which is an important spawning ground and habitat for S. maindroni as the research area, and used the MaxEnt (Maximum Entropy) model to analyze the habitat suitability of S. maindroni as well as to explore its restoration potential and spatial distribution of suitable habitats, based on the bottom trawl survey data from the Zhejiang coastal area from 2017−2021, the survey data from autonomous cruises in 2021−2022, and data from the Global Biodiversity Information Facility (GBIF). The research results show that the maximum sea surface temperature (SST), minimum sea surface salinity (SSS) and pH were the main ecological factors affecting the S. maindroni distribution. S. maindroni tended to be active in areas with higher sea surface temperatures and higher salinity, and could tolerate a weakly alkaline environment. In the Zhoushan sea area, S. maindroni was more likely to exist with maximum SST of 25.5−28.5 ℃, minimum SSS of 20‰−30‰, and pH value of 7.8−8.3. The suitable habitat area in the Zhoushan sea area was quite extensive, covering 30 213.40 km2, which accounted for 75.11% of the total area of the research area. Among these, the highly suitable habitats were mainly distributed in the sea areas of Zhongjieshan Islands, Shengsi Islands and Ma'an Islands, with an area of 4 600.78 km2. The areas of the sub-highly suitable, moderately suitable and lowly suitable habitats were 8 205.52, 9 723.31 and 7 683.78 km2, respectively. In terms of restoration potential, S. maindroni has significant recovery potential in the Zhoushan sea area.
-
Keywords:
- Sepiella maindroni /
- Habitat suitability /
- MaxEnt model /
- Zhoushan sea area
-
细菌性疾病会对水产养殖业造成重大经济损失，目前水产动物细菌性疾病的治疗仍以抗生素为主，但由于不科学的药物使用，目前水产病原菌的耐药问题十分严重[1-2]。乳酸菌作为一类普遍认同的益生菌，其生长过程中能产生细菌素、过氧化氢和有机酸等代谢产物，这些代谢产物能抑制病原菌的繁殖和生长[3]。因此，作为一种抗生素替代物的乳酸菌受到越来越多的关注。
乳酸菌应用于水产行业以来，研究多集中于促进水产动物生长、提高动物免疫力等方面，而其益生特性的研究则多集中于水产动物肠道定植、体外抑菌[4]以及对消化道低pH[5]和胆盐[6]环境的适应性方面。但乳酸菌是一类能利用可发酵糖产生大量乳酸的细菌统称，现有研究显示这类细菌中包含能分解明胶、导致溶血的致病菌[7]；此外，细菌耐药基因转移的情况普遍存在[8]，且已有研究指出大肠杆菌(Escherichia coli)和乳链球菌(Streptococcus lactis)之间能转移卡那霉素抗性[9]。因此，以具有溶血、分解明胶能力或耐药特性未知的乳酸菌菌株作为口服益生菌对机体有严重的安全隐患。目前对于水产领域应用乳酸菌的安全性评估相对薄弱，尤其是相关菌株耐药性的研究鲜有报道。故本研究对实验室保存的2株乳酸菌[乳酸菌S60 (干酪乳杆菌Lactobacillus casei) 和S72 (植物乳杆菌L. plantarum)]开展益生特性和安全性研究，评价其作为鱼类口服益生菌的可行性，为其后续的实际应用提供理论依据。
1. 材料与方法
1.1 主要材料和仪器
实验中所用乳酸菌S60、S72和指示菌嗜水气单胞菌 (Aeromonas hydrophila)、维氏气单胞菌 (A. veronii)、无乳链球菌 (S. agalactiae)、迟缓爱德华氏菌 (Edwardsiella tarda)、创伤弧菌 (Vibrio vulnificus)、副溶血弧菌 (V. parahemolyticus)均由本实验室分离保存，且经16S rDNA和常规细菌学鉴定。
1.1.1 试剂
LB培养基、MRS培养基、BHI培养基均购于青岛高科技工业园海博生物科技有限公司，明胶培养基 (杭州滨和微生物试剂有限公司)，血平板 (广东环凯微生物科技有限公司)，牛胆盐 (胆酸质量分数≥60%，叶源生物)，5%脱纤维绵羊血(南京森贝伽生物科技有限公司)，氨苄青霉素 (USP级)、硫酸卡那霉素 (USP级)、链霉素 (720 IU·mg−1)、恩诺沙星 (质量分数98%)、氟苯尼考 (质量分数98%) 和盐酸多西环素 (质量分数≥98%) 均购于上海阿拉丁生化科技股份有限公司，磺胺甲唑 (质量分数98%，上海麦克林生化科技有限公司)。血液生化试剂盒[(谷丙转氨酶(ALT)、谷草转氨酶 (AST)、乳酸脱氢酶 (LDH)]均购于上海执诚生物科技有限公司。
1.1.2 仪器
全波长酶标仪 (Multiskan G，Thermo Scientific)，高速冷冻离心机 (赛默飞世尔科技有限公司)，细菌浊度仪SGZ-6AXJ (上海悦丰仪器仪表有限公司)，全自动生化分析仪BS-400 (深圳迈瑞生物医疗电子股份有限公司)，智能生化培养箱 (无锡华泽科技有限公司)，震荡培养箱ZQL-180S (上海知楚仪器有限公司)。
实验鱼为异育银鲫 (Carassius auratus)，由中国水产科学研究院淡水渔业研究中心南泉基地提供，鱼体健康无伤，体质量为 (150±12) g，体长为 (17±3) cm (n=100)。实验前，异育银鲫于循环养殖系统中暂养2周，水温为 (26±1) ℃，持续充氧，光周期同自然光，使用商品化颗粒饲料饲喂，日投饵量按体质量的3%计算，一天3次。
1.2 生长曲线与产酸能力
取乳酸菌菌种于MRS培养基中过夜培养活化 (30 ℃、180 r·min−1)，次日将活化好的培养物100 μL接入5 mL新鲜MRS培养基 (30 ℃、180 r·min−1) 培养18 h，然后按1∶100稀释作为菌种，再以20 μL·mL−1接种，设置3个平行。分别在接种后的第0、第3、第6、第9、第12、第15、第24、第27、第30、第33和第36小时取菌液，使用酶标仪测定光密度 (OD600)。以时间为横坐标，OD600为纵坐标，绘制生长曲线。同时在接种后的第0、第12、第24、第36、第48、第60和第72小时检测发酵液pH，评价其产酸能力。
1.3 压力耐受能力
参照杨媛媛等[10]方法检测耐酸能力；参照任世英等[11]方法制备含乳酸菌的MRS平板以检测菌株胆盐耐受能力。
1.4 抑菌活性
1.4.1 抑菌能力评价
选择常见的水产病原菌——嗜水气单胞菌、维氏气单胞菌、无乳链球菌、迟缓爱德华氏菌、创伤弧菌、副溶血弧菌作为指示菌，参照杨媛媛等[10]方法，检测实验菌株的抑菌能力。最终将数据统一调整至菌斑直径为5 mm时的抑菌直径。计算公式为抑菌直径=(抑菌圈直径1/菌斑直径1+抑菌圈直径2/菌斑直径2+抑菌圈直径3/菌斑直径3)/3×5。
1.4.2 抑菌物质检测
参照任世英等[11]方法制备含指示菌的LB琼脂平板，对发酵液进行预处理以分析其所含抑菌物质特性。
1.5 安全性评价
1.5.1 明胶实验
将含实验菌的菌液接种于明胶培养基，30 ℃培养48 h，取出置于4 ℃ 30 min。取出检查明胶培养基的状态；同时设置阴性对照组 (以无菌生理盐水代替实验菌) 和阳性对照组 (以嗜水气单胞菌代替实验菌)。接种实验菌的明胶培养基与阴性对照一样为完全凝固的状态，表示实验菌不分解明胶，若同阳性对照一样完全呈液态则表明实验菌能分解明胶，介于两者之间呈半固体则表明实验菌有较弱的分解明胶的能力。
1.5.2 溶血实验
将活化的菌株划线接种于含有体积分数为5%脱纤维绵羊血的血琼脂平板上，30 ℃倒置培养24~48 h，观察记录菌落周围是否有溶血的透明圈。为避免实验菌产酸引起的溶血对血平板结果造成干扰，参照Mao等[12]的方法同时开展体外溶血实验。
1.5.3 药敏实验
同1.4.2制备含实验菌的MRS平板，打孔备用 (每个平板打孔5个)。以生理盐水配置抗生素质量浓度为256 μg·mL−1、128 μg·mL−1、64 μg·mL−1、32 μg·mL−1和16 μg·mL−1的药液，按150 μL每孔加入到之前准备好的含实验菌的MRS平板，每个浓度重复3次。4 ℃扩散4 h后，转入培养箱30 ℃ 培养48 h后，测量记录抑菌圈直径，判读其抗生素耐受能力。
1.5.4 动物实验
将过夜培养的乳酸菌和嗜水气单胞菌，以5 000 r·min−1离心10 min收集菌体，PBS缓冲液洗涤3次。以PBS缓冲液调节乳酸菌浓度为1.0×109 CFU·mL−1，嗜水气单胞菌浓度为5.0×106 CFU·mL−1。设置4个处理组，每组3个平行，每个平行10尾鱼[空白对照组、S60组、S72组和嗜水气单胞菌(A.h)组]。对实验组每尾银鲫腹腔注射0.5 mL菌液，阴性对照组腹腔注射0.5 mL PBS缓冲液。每天定时检查实验鱼的死亡情况，持续观察96 h。96 h后将实验鱼捞出，麻醉 (MS-222 100 mg·mL−1) 检查鱼体，尾椎静脉取血，4 ℃、5 000 r·min−1离心10 min收集血浆供血液生化检测 (ALT、AST和LDH)，剖检成活的实验鱼，观察内脏器官有无病变。
1.6 数据统计
实验数据差异性分析使用SPSS 24.0统计软件中的单因素方差分析 (One-Way ANOVA) 和独立样本t-检验分析，显著性水平设置为P<0.05。
2. 结果
2.1 生长曲线
乳酸菌S60和S72的生长曲线呈典型的“S”形，可以观察到明显的迟缓期、对数生长期和稳定期 (图1)。从S72生长曲线上可以发现其迟缓期约为6 h，S60的迟缓期约为9 h，明显长于S72。表明在相同培养条件下，S72具有更快的生长速度，这与菌落形态方面的结果一致。
2.2 产酸能力
发酵液pH检测结果显示，乳酸菌S60和S72这2株的产酸最大值均出现在培养第24小时附近，2株乳酸菌的产酸能力基本一致，最低pH均约4.0 (图2)。
2.3 压力耐受能力
耐酸实验结果显示乳酸菌S60和S72经pH为4.5的PBS处理2 h后均能成活，由此显示2株乳酸菌能耐受pH为4.5的酸性环境。胆盐耐受实验结果显示，乳酸菌S60和S72能完全耐受0.1%的胆盐浓度，而在胆盐浓度>0.2%时，其生长受到了一定程度的抑制，可以观察到明显的抑菌圈 (表1)。
表 1 乳酸菌对酸、胆盐的耐受能力Table 1. Acid and bile tolerance of two strains of lactic acid bacterian=3 pH 乳酸菌 lactic acid bacteria 胆盐/%
bile salt抑菌圈直径/mm inhibition zone diameter S60 S72 S60 S72 3.5 − − 0 8 8 4.5 + + 0.1 8 8 5.5 + + 0.2 12 12 7.4 + + 0.3 15 16 0.5 19 19 注：−. 无菌生长；+. 有菌生长
Note: −. bacteria can not grow; +. bacteria can grow2.4 抑菌物质检测
用嗜水气单胞菌作为指示菌，以1 mol·L−1 氢氧化钠 (NaOH)加热、胰蛋白酶和蛋白酶K分别处理发酵液上清作为实验组，未处理的发酵上清液为对照组，测定抑菌圈大小。结果显示，对照组有明显的抑菌圈，而实验组无。说明2株乳酸菌主要通过发酵过程中产生的有机酸 (低pH)实现了对嗜水气单胞菌的抑制作用。
2.5 抑菌活性
抑菌实验结果显示，2株乳酸菌均能不同程度地抑制6株指示菌的生长 (表2)，其中对嗜水气单胞菌、维氏气单胞菌和创伤弧菌的抑菌能力最强，对副溶血弧菌的抑菌效果次之，对无乳链球菌和迟缓爱德华氏菌的抑菌能力最弱。2株乳酸菌在同一种病原菌中的抑菌效果差异不显著 (P>0.05)，而同一乳酸菌菌株对不同病原菌的抑菌效果差异显著 (P<0.05)。
表 2 乳酸菌对鱼源致病菌的抑菌活性Table 2. Antimicrobial activity of two strains of lactic acid bacteria against six fish pathogens$\overline{{\mathit{\boldsymbol{ X }}}} {\bf \pm} {\bf{SD,}}\;{\mathit{\boldsymbol{n}}}{\bf = 5} $ 致病菌
pathogens抑菌圈直径/mm inhibition zone diameter S60 S72 嗜水气单胞菌 A. hydrophila 26.88±2.43b 24.9±1.48C 维氏气单胞菌 A. veronii 28.85±2.45b 28.25±2.46B 创伤弧菌 V. vulnificus 33.87±1.96a 35.6±1.07A 无乳链球菌 S. agalactiae 8.33±2.38d 9.7±0.96E 副溶血弧菌 V. parahemolyticus 15.23±0.78c 13.1±1.82D* 迟缓爱德华氏菌 E. tarda 9.23±2.42d 13.05±1.92D* 注：不同大小写字母表示2株乳酸菌分别针对不同病原菌的抑菌能力的差异显著(P<0.05)；*. 2株乳酸菌对同一病原菌的抑菌能力差异显著(P<0.05)Note: Different letters indicate statistical difference in bacteriostatic ability of the same lactobacillus strain among different pathogenic bacteria (P<0.05); *. statistical difference in bacteriostatic ability of different lactobacillus strains on the same pathogenic bacteria (P<0.05) 2.6 安全性评价
划线接种血平板24 h后，2株乳酸菌形成的菌落周围无明显变化，阳性对照组 (嗜水气单胞菌) 菌落周围出现了明显的透明圈。第48小时菌落周围颜色变淡，无明显透明圈出现，阳性对照组的透明圈进一步增大 (图3)。体外溶血实验显示，2株乳酸菌发酵液的溶血率与对照组无显著差异，溶血率约5.0%，而嗜水气单胞菌发酵液引起的溶血率显著高于对照组 (P<0.05)，溶血率约72.0% (图4)。由此可见2株乳酸菌不存在溶血相关的毒力因子，也不存在溶血方面的隐患。
![]()
图 3 2株乳酸菌在血平板上的溶血表现Figure 3. Hemolysis activity of two strains of lactic acid bacteria on blood plate2株乳酸菌接入明胶培养基培养48 h后，经4 ℃处理，与阴性对照组一样仍然保持凝固状态，而接种指示菌 (嗜水气单胞菌) 的阳性对照组呈液态 (图5)。由此显示2株乳酸菌均无降解明胶的能力，不存在相关风险。
![]()
图 5 2株乳酸菌分解明胶的能力Figure 5. Gelatinolytic activity of two strains of lactic acid bacteria药敏实验结果显示，2株乳酸菌总体上对硫酸卡那霉素、链霉素和磺胺甲唑耐药，对氨苄青霉素、恩诺沙星、氟苯尼考和多西环素敏感。然而通过进一步的比较分析则发现2株乳酸菌对2种氨基糖胺类抗生素耐药；但S72对2种氨基糖胺类抗生素的高浓度剂量仍有一定的敏感性，对恩诺沙星低浓度剂量则表现出完全耐受 (表3)。
表 3 2株乳酸菌的耐药特性Table 3. Results of antibiotic susceptibility test of two strains of lactic acid bacteria菌株
bacterial strain药物质量浓度/μg·mL−1 drug concentration S60 S72 抑菌圈直径/mm
inhibition zone diameter256 128 64 32 16 256 128 64 32 16 氨苄青霉素 ampicillin 30 26 24 20 17 40 35 32 30 28 硫酸卡那霉素 kanamycin 8 8 8 8 8 11 8 8 8 8 链霉素 streptomycin 8 8 8 8 8 12 9 8 8 8 恩诺沙星 enoxacin 26 21 19 17 10 19 16 13 8 8 氟苯尼考 flufenicol 32 29 26 20 15 34 29 25 24 18 盐酸多西环素 doxycycline 37 34 28 24 21 34 28 24 21 16 磺胺甲唑 sulfamethoxazole 8 8 8 8 8 8 8 8 8 8 生理盐水 physiological saline 8 8 8 8 8 8 8 8 8 8 鱼体攻毒后，S72组仅在第72小时出现了1尾死鱼，死亡鱼类注射部位有炎症和出血症状，体表其他部分也有广泛的出血现象。其他时间点并未出现死鱼现象。对照组和S60组在整个实验期间均无死鱼出现。阳性对照组在攻毒24 h后陆续出现死鱼的情况，在实验结束时的累计死亡率达46.67% (图6)，且死亡鱼类在注射部位有炎症和出血现象，鱼体其他部位也有严重的出血现象。96 h后对存活的实验鱼进行检查，结果显示S60和S72组剩余的实验鱼体表在注射部位均有轻微的炎症迹象，而进一步解剖检查显示，鱼体腹腔无腹水存在，体内脏器无明显病变情况。
攻毒实验96 h后，鱼类血液生化检测结果显示S60和S72组鱼类血液ALT、AST和LDH 的水平和对照组无显著差异 (P>0.05，表4)。说明S60和S72攻毒并未造成鱼类组织损伤。而高达1.0×109 CFU·mL−1的细菌腹腔注射仅出现了1尾实验鱼的死亡，且第96小时血液生化检查结果显示剩余鱼类并未出现明显的组织损伤；与此同时嗜水气单胞菌以5.0×106 CFU·mL−1的浓度攻毒，却导致近50%的实验鱼死亡。由此可见，2株乳酸菌仍具有较高的安全性。
表 4 2株乳酸菌对银鲫血液生化指标的影响Table 4. Effects of two strains of lactic acid bacteria on blood biochemical index of C. auratus challenged with two lactic acid bacteriaU·L−1；n=5 项目
index对照组
control乳酸菌 lactic acid bacteria 嗜水气单胞菌组
A.hS60 S72 谷丙转氨酶 ALT 6.39±1.94 6.19±3.24 4.96±1.84 4.21±1.19 谷草转氨酶 AST 323.52±65.91 321.44±105.27 391±133.81 330.55±118.34 乳酸脱氢酶 LDH 116.53±38.97 88.54±43.39 65.14±37.26 73.35±31.90 3. 讨论
3.1 乳酸菌对鱼类肠道环境的适应性
消化道中的pH和胆盐对细菌繁殖生长具有抑制作用，因此益生菌在消化道中的益生作用与其对消化道环境的耐受能力密切相关。人类医学研究表明消化道pH受饮食组成影响，食物的存在能提高益生菌在胃内的成活率[13]。此外，将壳聚糖等生物高分子与菌剂混合应用也能显著提高益生菌在低pH条件下的成活率[14]。而现有资料显示有胃鱼类的胃内pH最低能达2[15]，因此本研究中的2株乳酸菌需要辅以饵料、生物高分子或经过一段时间的耐酸驯化[13]才能更好地应用于有胃鱼类。而无胃鱼类，如鲤科鱼类，其消化道内pH近中性[10]，因此本研究中的2株乳酸菌能适应其消化道pH。
胆汁能阻止细菌黏附到肠道黏膜的顶端，避免细胞损伤[6]，但也对乳酸菌的生存产生了胁迫。本研究中的2株乳酸菌能耐受0.1%的胆盐，而参考海水鱼类数据，胆囊内胆盐含量在0.2%~1.0%，进入肠道会被进一步稀释[15]。并且胆盐质量浓度在0.03%~0.3%范围内波动的肠道环境，尚能分离到不耐受0.1%胆盐浓度的乳酸菌[13]。因此认为2株乳酸菌能耐受鱼类肠道中的胆盐环境。
3.2 乳酸菌的生长特性、产酸能力和抑菌活性物质
乳酸菌是目前公认的肠道有益微生物，大量研究显示乳酸菌对于肠道微生物群落稳定[14]、与人类[16]和水产养殖动物[4]的健康密切相关。虹鳟 (Oncorhynchus mykiss)[17]、凡纳滨对虾(Litopenaeus vannamei)[18]和海参 (Apostichopus japonicus)[19]等养殖实验均显示在日粮中添加乳酸菌能促进养殖动物的生长，提高抗病力和肠道内有害菌产生竞争，发挥保护肠道的作用[20]。本研究结果显示乳酸菌S72比S60的繁殖速度快，更易取得竞争优势，和其他细菌竞争肠道上皮细胞的黏附位点[4,21]、营养物质[8]以及减少病原菌的定植[22]；此外，乳酸菌能通过产生乳酸[23-24]、细菌素、过氧化氢及小分子肽类等[20,25]代谢产物对肠道致病菌产生抑菌或杀菌作用，而快速繁殖的能力也意味着更多活性物质的产生，从而发挥更强的抑菌或杀菌效果。本研究结果显示2株乳酸菌产生的抑菌物质均为代谢过程中产生的有机酸类，与其他关于乳酸菌抑菌物质的研究结果[3,23]类似，并非每一株乳酸菌都能产生所有类型的抑菌物质，可见乳酸菌产生的抑菌物质种类和菌株、菌种有关。而Kiymaci等[23]发现乳酸菌所产生的有机酸对病原菌不仅有抑制效果，而且能削弱病原菌的毒力。此外，胞外的一些成分如多糖、蛋白等也具有改善机体免疫[26]和吸附肠道内有害物质[27]的能力，但这些潜在功能有待研究论证。
3.3 乳酸菌在水产动物中应用的安全性问题
溶血素是病原菌致病的重要毒力因子，研究表明溶血素不仅表现为溶解红细胞，其具备的穿孔能力、酯酶活性或表面活性作用能损伤多种细胞，诱发细胞凋亡甚至造成细胞裂解死亡[28]。而乳酸菌是一类能利用糖产生乳酸的细菌统称，其中某些种属具有溶血能力，如闫肃等[7]研究发现，所检视的48株乳酸菌中有9株具有溶血活性。因此，需要对候选益生菌的溶血能力进行检测，从而避免相关毒力的风险。本研究中2株乳酸菌溶血活性的检测结果显示，在血平板实验培养后期 (48 h) 乳酸菌菌落周围有颜色变淡的情况出现，与闫肃等[7]研究中将近68.75%的乳酸菌菌株 (检视菌株总数48株) 的溶血表型一致，因此认为这一现象与低pH相关，且在进一步的体外溶血实验中得到证实。由此可见，2株乳酸菌不产生溶血活性物质。
抗菌药物在保证人类和动物健康的同时[8]，也带来了严重的耐药问题，现有资料显示，水产养殖动物主要细菌性病原的耐药现象普遍存在[1,29]。恩诺沙星、氟苯尼考、多西环素和磺胺甲唑是目前用于水产养殖动物细菌性疾病防治的主要药物[30]。鉴于乳酸菌可能涉及的青霉素类[31]、氨基糖苷类[32]和四环素类[33]的耐药基因及其耐药基因可转移的风险。本文结合药敏实验结果，首先排除了本研究中2株乳酸菌在恩诺沙星、氟苯尼考和多西环素耐药性转移方面的风险。而磺胺类药物的耐药现象是否和获得外源性耐药基因 (sul1、sul2和sul3)[33-34]相关需进一步研究。尤其是确认是否携带sul1耐药基因，因为该耐药基因的扩散传播与整合子相关，不仅水平转移概率高，而且拥有更广泛的宿主范围[35]。本研究明确了候选益生菌的耐药特性不仅能有效的规避耐药转移风险，且为科学用药提供指导意见。
在肠道菌群稳态失衡时，常常会出现细菌位移，出现非致病菌感染致病的结果[36]。虽然鲜有哺乳动物食源性乳酸菌引起疾病的报道[8]，但应用到水产养殖动物中的一些菌株却暴露出一些安全隐患，如一株植物乳杆菌[37]在应用过程中对达氏鳇(Huso huso)[37]引起了组织损伤。因此，为评估2株乳酸菌潜在的条件致病风险，本研究采取腹腔注射的形式，以高达1.0×109 CFU·mL−1的浓度进行了动物攻毒实验。结果表明乳酸菌S60在整个养殖实验过程中未引起实验鱼的死亡；而乳酸菌S72仅导致了1尾实验鱼死亡 (死亡率3.33%)，相较于浓度仅为5.0×106 CFU·mL−1的嗜水气单胞菌 (典型条件致病菌) 所造成将近50%的实验鱼死亡，2株乳酸菌的安全性仍然很高。并且在攻毒处理后96 h的血液生化分析结果表明，实验组鱼类和对照组鱼类在组织损伤方面并未表现出显著差异，由此可见，在第96小时实验鱼基本已恢复到正常状态。综上，2株乳酸菌具备抑制常见鱼类病原菌的能力，对银鲫无明显致病性，可作为银鲫口服益生菌开发的资源储备。
-
![]()
图 3 刀切法检验生态因子对曼氏无针乌贼分布影响的重要程度
Figure 3. Importance of ecological factors on distribution of S. maindroni by Jackknife test
![]()
图 5 曼氏无针乌贼在舟山海域的适生区分布
Figure 5. Distribution of suitable habitat of S. maindroni in Zhoushan Islands
![]()
图 6 舟山海域内海洋特别保护区地理位置图
Figure 6. Geographical location map of marine special reserves in Zhoushan sea area
表 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 
下载: 导出CSV
表 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
下载: 导出CSV
表 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
下载: 导出CSV
-
[1] RICHARDSON L E, LENFANT P, CLARKE L J, et al. Examining current best-practices for the use of wild post-larvae capture, culture, and release for fisheries enhancement[J]. Front Mar Sci, 2023, 9: 1058497. doi: 10.3389/fmars.2022.1058497
[2] 钟明, 侍昊, 安树青, 等. 中国野生动物生境适宜性评价和生境破碎化研究[J]. 生态科学, 2016, 35(4): 205-209. doi: 10.14108/j.cnki.1008-8873.2016.04.028 [3] 吕国敏, 吴进锋, 陈利雄. 我国头足类增养殖研究现状及开发前景[J]. 南方水产, 2007, 3(3): 61-66. [4] HEIKKINEN R K, LUOTO M, ARAÚJO M B, et al. Methods and uncertainties in bioclimatic envelope modelling under climate change[J]. Prog Phys Geogr, 2006, 30(6): 751-777. doi: 10.1177/0309133306071957
[5] 徐晓萱, 谢玉, 刘姝含, 等. 基于GAM模型的浙江近海曼氏无针乌贼时空分布研究[J]. 浙江海洋大学学报 (自然科学版), 2022, 41(5): 400-407. [6] 魏裙倚, 崔国辰, 玄文丹, 等. 海表面温度及叶绿素a浓度对西北印度洋鸢乌贼时空分布的影响[J]. 中国水产科学, 2022, 29(3): 388-397. [7] 李杰, 张鹏, 晏磊, 等. 南海中南部海域鸢乌贼CPUE影响因素的GAM分析[J]. 中国水产科学, 2020, 27(8): 906-915. [8] 杨洋. 基于四种模型的肉苁蓉潜在地理分布预测及空间格局变化分析[D]. 西安: 陕西师范大学, 2017: 25-28. [9] 图雅. 基于物种分布模型的中国针茅属植物潜在分布区及其影响因子分析[D]. 北京: 北京林业大学, 2020: 72-73. [10] MEROW C, SILANDER J A. A comparison of Maxlike and Maxent for modelling species distributions[J]. Methods Ecol Evol, 2014, 5(3): 215-225. doi: 10.1111/2041-210X.12152
[11] 陈新美, 雷渊才, 张雄清, 等. 样本量对MaxEnt模型预测物种分布精度和稳定性的影响[J]. 林业科学, 2012, 48(1): 53-59. doi: 10.11707/j.1001-7488.20120110 [12] TENG S, SU N, LEE M, et al. Modeling the habitat distribution of Acanthopagrus schlegelii in the coastal waters of the Eastern Taiwan Strait using MaxEnt with fishery and remote sensing data[J]. J Mar Sci Eng, 2021, 9(12): 1442. doi: 10.3390/jmse9121442
[13] HARTE J, NEWMAN E A. Maximum information entropy: a foundation for ecological theory[J]. Trends Ecol Evol, 2014, 29(7): 384-389. doi: 10.1016/j.tree.2014.04.009
[14] 杨继超, 董民星, 种衍飞, 等. 最大熵模型在海洋生物适生区预测中的应用[J/OL]. 应用海洋学学报. https://kns.cnki.net/kcms/detail/35.1319.P.20230302.1347.002.html. [15] 周海涛, 那晓东, 臧淑英, 等. 最大熵 (MaxEnt) 模型在物种栖息地研究中的应用[J]. 环境科学与管理, 2016, 41(3): 149-151. doi: 10.3969/j.issn.1673-1212.2016.03.035 [16] 张嘉容, 杨晓明, 田思泉. 基于最大熵模型的南太平洋长鳍金枪鱼栖息地预测[J]. 中国水产科学, 2020, 27(10): 1222-1233. [17] 龚彩霞, 陈新军, 高峰, 等. 水温变暖对西北太平洋柔鱼潜在栖息地分布的影响[J]. 海洋学报 (中文版), 2022, 44(7): 95-102. [18] PIRTLE J L, SHOTWELL S K, ZIMMERMANN M, et al. Habitat suitability models for groundfish in the Gulf of Alaska[J]. Deep-Sea Res II, 2017, 165: 303-321.
[19] 曹子豪, 迟长凤, 刘慧慧, 等. 不同地理群体曼氏无针乌贼肌肉营养成分分析比较与评价[J]. 食品科学, 2015, 36(4): 101-105. doi: 10.7506/spkx1002-6630-201504019 [20] 倪正雅, 徐汉祥. 浙江近海乌贼资源评估及乌贼渔业管理[J]. 海洋渔业, 1986(2): 51-54. [21] 李继姬, 郭宝英, 吴常文. 浙江海域曼氏无针乌贼资源演变及修复路径探讨[J]. 浙江海洋学院学报 (自然科学版), 2011, 30(5): 381-385, 396. [22] 覃涛, 俞存根, 陈全震, 等. 舟山渔场及邻近海域头足类 (Cephalopod) 种类组成和数量分布[J]. 海洋与湖沼, 2011, 42(1): 124-130. doi: 10.11693/hyhz201101019019 [23] 陈伟峰, 叶深, 余玥, 等. 浙南近海头足类种类组成及生态位分析[J]. 水生生物学报, 2021, 45(2): 428-435. doi: 10.7541/2021.2018.262 [24] PANG Y M, TIAN Y J, FU C H, et al. Variability of coastal cephalopods in overexploited China seas under climate change with implications on fisheries management[J]. Fish Res, 2018, 208: 22-33. doi: 10.1016/j.fishres.2018.07.004
[25] 张义浩, 王志铮, 吴常文, 等. 舟山群岛定生海藻种类组成、生态分布及区系特征研究[J]. 浙江海洋学院学报 (自然科学版), 2002, 21(2): 98-105. [26] 刘连为, 隋宥珍, 徐开达, 等. 曼氏无针乌贼两个地理群体繁殖生物学比较研究[J]. 海洋湖沼通报, 2022, 44(4): 57-62. doi: 10.13984/j.cnki.cn37-1141.2022.04.008 [27] WARREN D L, MATZKE N J, CARDILLO M, et al. ENMTools 1.0: an R package for comparative ecological biogeography[J]. Ecography, 2021, 44(4): 504-511. doi: 10.1111/ecog.05485
[28] 尹飞, 王春琳, 宋微微. 曼氏无针乌贼幼体生态因子耐受性的研究[J]. 湛江海洋大学学报, 2005, 25(4): 39-43. [29] 黄伟卿, 陈宇光, 张艺, 等. 盐度胁迫对曼氏无针乌贼胚胎发育和生长性能的影响[J]. 渔业现代化, 2022, 49(6): 100-107. doi: 10.3969/j.issn.1007-9580.2022.06.013 [30] 郑美丽, 肖金华, 郑微云, 等. 曼氏无针乌贼的趋光特性[J]. 厦门大学学报 (自然科学版), 1980(3): 91-99. [31] GUISAN A, ZIMMERMANN N E. Predictive habitat distribution models in ecology[J]. Ecol Model, 2000, 135(2): 147-186.
[32] 谢平, 赵羽西, 桑燕芳, 等. 基于相关系数的水文周期变异分级方法及验证[J]. 水力发电学报, 2018, 37(12): 33-43. doi: 10.11660/slfdxb.20181204 [33] COBOS M E, PETERSON A T, BARVE N, et al. kuenm: an R package for detailed development of ecological niche models using MaxEnt[J]. PeerJ, 2019, 7: e6281. doi: 10.7717/peerj.6281
[34] 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
[35] 阮欧, 刘绥华, 陈芳, 等. 基于多源遥感的贵州草海国家级自然保护区黑颈鹤生境适宜性评价[J]. 生态学报, 2022, 42(5): 1947-1957. [36] MASTRANDREA M D, FIELD C B, STOCKER T F, et al. Guidance Note for Lead Authors of the IPCC Fifth Assessment Report on Consistent Treatment of Uncertainties[C]. The IPCC Panel, Intergovernmental Panel on Climate Change (IPCC), CA, USA, IPCC Cross-Working Group Meeting on Consistent Treatment of Uncertainties Jasper Ridge, 2010: 2-3.
[37] PETERSON A T, PAPE M, SOBERÓN J. Rethinking receiver operating characteristic analysis applications in ecological niche modeling[J]. Ecol Model, 2017, 213(1): 63-72.
[38] 王运生, 谢丙炎, 万方浩, 等. ROC曲线分析在评价入侵物种分布模型中的应用[J]. 生物多样性, 2007, 15(4): 365-372. doi: 10.3321/j.issn:1005-0094.2007.04.005 [39] 曹杰, 陈新军, 刘必林, 等. 鱿鱼类资源量变化与海洋环境关系的研究进展[J]. 上海海洋大学学报, 2010, 19(2): 232-239. [40] 肖戈, 徐博, 张衡, 等. 阿拉伯海外海鸢乌贼渔场时空分布与海洋环境要素的研究[J]. 南方水产科学, 2022, 18(4): 10-19. doi: 10.12131/20210217 [41] 陈峰, 李楠, 方舟, 等. 浙江近岸海域春夏季剑尖枪乌贼栖息地分布变化规律[J]. 上海海洋大学学报, 2021, 30(5): 847-855. doi: 10.12024/jsou.20201203245 [42] 吴天, 梁君, 周永东, 等. 中街山列岛海域曼氏无针乌贼栖息地偏向性研究[J]. 浙江海洋大学学报 (自然科学版), 2022, 41(5): 408-417. [43] 刘姝含, 刘连为, 徐开达, 等. 浙江中北部近海曼氏无针乌贼资源现状[J]. 浙江海洋大学学报 (自然科学版), 2022, 41(4): 286-293. [44] 袁健美, 张虎, 贲成恺, 等. 江苏近海曼氏无针乌贼增殖放流效果评估[J]. 水产养殖, 2023, 44(2): 12-16. doi: 10.3969/j.issn.1004-2091.2023.02.003 [45] 张建设. 曼氏无针乌贼养殖生物学特性和血细胞免疫功能研究[D]. 厦门: 厦门大学, 2007: 155. [46] YIN F, SUN P, PENG S M, et al. The respiration, excretion and biochemical response of the juvenile common Chinese cuttlefish, Sepiella maindroni at different temperatures[J]. Aquaculture, 2013, 402/403: 127-132. doi: 10.1016/j.aquaculture.2013.03.018
[47] PENG S M, YIN F, SHI Z H, et al. Optimum water temperature for the growth of juvenile common Chinese cuttlefish, Sepiella maindroni (de Rochebrune 1884)[J]. J Shellfish Res, 2011, 30(2): 205-209. doi: 10.2983/035.030.0202
[48] 雷舒涵, 吴常文, 高天翔, 等. 金乌贼和曼氏无针乌贼胚胎发育及其盐度耐受能力的比较研究[J]. 中国水产科学, 2011, 18(2): 350-359. [49] HALVORSEN R, MAZZONI S, BRYN A, et al. Opportunities for improved distribution modelling practice via a strict maximum likelihood interpretation of MaxEnt[J]. Ecography, 2015, 38(2): 172-183. doi: 10.1111/ecog.00565
[50] 朱海洋. 舟山渔场: 从“赶海”到“养海”[N]. 农民日报, 2023-02-16(6). [51] 宋旻鹏, 汪金海, 郑小东. 中国经济头足类增养殖现状及展望[J]. 海洋科学, 2018, 42(3): 149-156. doi: 10.11759/hykx20180128001 [52] 季乾昭, 王荣兴, 黄志旁, 等. 样本量与研究范围变化对MaxEnt模型准确度的影响: 以黑白仰鼻猴为例[J]. 兽类学报, 2019, 39(2): 126-133. -
期刊类型引用(14)
1. 马雪莲，刘丹，杨希，马睿，高强. 一株海氏肠球菌的分离鉴定及益生效果评价. 微生物学杂志. 2025(01): 51-58 . 
百度学术
2. 李秋月，杨雪娇，谭春艳，陈家裕，林连兵，邓先余. 罗非鱼源肠道益生菌的分离鉴定及其体外益生特性分析. 食品与生物技术学报. 2023(12): 1-9 . 百度学术
3. 张建飞. 一株产细菌素粪链球菌N9301的分离鉴定及生物学特性研究. 饲料研究. 2022(08): 78-82 . 百度学术
4. 王鸿蕉，蔡昀洁，岳明祥，俞佳，刘祥. 植物乳杆菌的抑菌性能及其影响因素研究进展. 食品研究与开发. 2022(14): 212-218 . 百度学术
5. 聂紫玉，吴艳阳，王增光，李子晗，康文丽，潘丽娜，汪家琦，戴智勇，赵玲艳，邓放明. 植物源益生乳酸菌的筛选及其特性. 食品科学. 2022(18): 143-151 . 百度学术
6. 张甜甜，鹿瑶，吕文亮，林显华，徐海燕，谷巍. 具有沙门氏菌拮抗能力的植物乳杆菌益生特性评价. 河南农业科学. 2021(02): 130-137 . 百度学术
7. 郭玉媛，张鑫，张瑞婷，金铁岩. 低温储存对2种朝鲜族传统米酒品质的影响. 食品工业. 2021(05): 76-79 . 百度学术
8. 万倩，李启明，吴华星，刘绒梅，曹珺，黄晓英，于基成，唐俊妮. 传统发酵食品中乳酸菌的安全性评估. 现代食品科技. 2021(06): 276-286 . 百度学术
9. 曲木，唐子鹏，赵子续，翟胜利，张宝龙. 乳酸菌在水产养殖中的应用. 生物化工. 2021(03): 131-134 . 百度学术
10. 王楠，尹纪元，王英英，李莹莹，吴斯宇，石存斌，李家豪，曹际振，王庆. 草鱼源乳酸菌的分离鉴定及其生物学特性研究. 南方水产科学. 2021(06): 74-84 . 本站查看
11. 杨晓萍，陈梓琦，钟先锋. 乳酸菌的分离鉴定方法及益生功能应用研究. 农产品加工. 2020(13): 50-54 . 百度学术
12. 赵永慧，杨杰，朱卿正，郭长明，柴金龙，王淑军，房耀维，盘赛昆. 一株植物乳杆菌MMB-03的筛选鉴定及抑菌特性研究. 安徽农业科学. 2020(17): 189-192+196 . 百度学术
13. 艾丹，常若葵，张斌，于玮玮，王远宏. 美极梅奇酵母T-2灌胃对小鼠的毒性作用. 食品研究与开发. 2020(22): 1-6 . 百度学术
14. 郭同旺，周书洪，郑守伟，吕志敏，张云泽，李炳乾，辛乃宏. 几种常见水产养殖益生菌的应用研究现状及展望. 盐科学与化工. 2020(11): 3-7 . 百度学术
其他类型引用(15)
计量
- 文章访问数: 395
- HTML全文浏览量: 136
- PDF下载量: 78
- 被引次数: 29
粤公网安备 44010502001741号
