水体投放枯草芽孢杆菌对不同习性鱼类肠道菌群和水质影响的研究.pdf

渔业研究，（）：收稿日期：基金项目：国家现代农业产业技术体系建设专项（）作者简介：张晶晶（），女，助理农艺师，硕士，从事水产养殖病害研究。：通信作者：樊海平（），男，研究员，硕士，从事水产养殖病害研究。：张晶晶，林 煜，陈度煌，等 水体投放枯草芽孢杆菌对不同习性鱼类肠道菌群和水质影响的研究 渔业研究，（）：水体投放枯草芽孢杆菌对不同习性鱼类肠道菌群和水质影响的研究张晶晶，林 煜，陈度煌，陈 斌，薛凌展，樊海平（福建省淡水水产研究所，福建 福州）摘要：为了探讨水体投放枯草芽孢杆菌（）对黑脊倒刺鲃（）、半刺厚唇鱼（）和鲢（）等不同习性鱼类肠道微生物菌群的影响以及对养殖水体的净化效果，实验采用 测序分析了该 种实验鱼前、中、后肠的菌群结构，并测定分析养殖水体、氨氮和亚硝酸盐含量变化。结果表明：投放枯草芽孢杆菌 起能够稳定降低水体中 氨氮含量，自 起能够稳定降低水体中 亚硝酸盐含量。肠道菌群分析结果显示：枯草芽孢杆菌实验组和对照组各肠道部位菌群整体区分明显，实验组肠道菌群丰度降低。整体而言，实验组放线菌门呈上升趋势，而浮霉菌门呈下降趋势，芽孢杆菌属在各肠道丰度提高，常见水产致病菌气单胞属和弧菌属丰度均降低。信号通路富集结果显示，鲢投放枯草芽孢杆菌组肠道菌群代谢功能加强。水体投放枯草芽孢杆菌实验证明，枯草芽孢杆菌可以降低水体氨氮、亚硝酸盐含量，起到净化水质的作用，同时可抑制鱼体肠道有害菌群的繁殖，提高鲢肠道菌群代谢水平。关键词：枯草芽孢杆菌；水体净化；肠道菌群；黑脊倒刺鲃；半刺厚唇鱼；鲢中图分类号：文献标识码：文章编号：（）枯草芽孢杆菌（），隶属于厚壁菌门（）、芽孢杆菌科（）、芽孢杆菌属（），具有耐高温、耐酸、耐碱、稳定性高的特性。养殖水体中投放适量的枯草芽孢杆菌可以降低水体亚硝酸盐含量，维持水体 及溶解氧含量，起到净化水质的作用；饲料中添加枯草芽孢杆菌可提高罗非鱼（）的生长速度、消化酶活性以及先天性免疫反应和抗病力；口服芽孢杆菌可以提高南亚野鲮（）的呼吸暴发活性以及先天性免疫反应。作为水产养殖上应用较广泛、极具潜力的有益菌种之一，枯草芽孢杆菌可以抑制有害菌生长，调整肠道微生态菌群，提高消化酶，促进脂肪代谢，降低饵料系数，提高鱼体对营养物质的吸收。国内外关于枯草芽孢杆菌改善鱼类肠道微生物的报道常见于日粮添加，而水体中投放该菌常被应用于对虾养殖研究，其他鱼类则较为少见。本研究选取淡水中下层偏草食性杂食鱼类黑脊倒刺鲃（）、底栖杂食性鱼类半刺厚唇鱼（）、中上层滤食性鱼类鲢（）为研 第 期张晶晶等：水水体体投投放放枯枯草草芽芽孢孢杆杆菌菌对对不不同同习习性性鱼鱼类类肠肠道道菌菌群群和和水水质质影影响响的的研研究究究对象，于其养殖水体中投放枯草芽孢杆菌，探讨该菌对上述 种生活于不同水层、不同食性鱼类的肠道菌群的影响，旨在为后续相关水体投放枯草芽孢杆菌对鱼类肠道的研究提供参考。材料与方法 实验材料于 年 月 日 月 月在福建省淡水水产研究所榕桥中试基地开展实验，水源为水库水。实验用鱼均由福建省淡水水产研究所榕桥中试基地提供，健康黑脊倒刺鲃（记为）平均体质量和平均体长分别为（）和（）；半刺厚唇鱼（记为）平均体质量和平均体长分别为（）和（）；鲢（记为）平均体质量和平均体长分别为（）和（）。枯草芽孢杆菌由北京渔经生物技术有限公司提供，芽孢含量达 。投喂粗蛋白含量为的淡水鱼膨化颗粒配合饲料，配合饲料购自天马科技集团股份有限公司。鱼种放养与日常管理实验容器为 玻璃纤维桶，每桶装水，混养上述 种实验鱼各 尾，连续充气，暂养，待稳定后开始实验，实验设置投放芽孢杆菌的实验组（，）和不投放芽孢杆菌的对照组（，）。实验组按照 枯草芽孢杆菌粉剂的水体浓度进行泼洒，使用前以枯草芽孢杆菌 红糖 水 比例进行洁净水活化。实验期间上午、下午各投喂饲料 次，日投饵量约为 种鱼体总体质量的，日换水 次，换水量为，实验组换水后，仍按照 枯草芽孢杆菌剂量并活化 后再泼洒，将水体中枯草芽孢杆菌含量维持在 ，对照组不投放芽孢杆菌。每天早上：，在投喂前检测养殖水体水温、值、氨氮和亚硝酸盐含量，检测仪器为格蓝普智能多参数水质检测仪（型号 ）。肠道菌群 测序实验后，随机采集实验组和对照组养殖桶黑脊倒刺鲃、半刺厚唇鱼和鲢各 尾，合计 尾，分别取鱼体的前肠（，）、中肠（，）和后肠（，）内容物，黑脊倒刺鲃实验组和对照组前、中、后肠分别命名为、和 、；半刺厚唇鱼实验组和对照组前、中、后肠分别命名 、和、；鲢实验组和对照组 前、中、后 肠 分 别 命 名 为 、和 、（以下同此）。液氮速冻后 冰箱保存，送至杭州联川生物技术股份有限公司进行 测 序 分 析。扩 增 引 物 为：；：，测 序 平 台 为 测序仪，测序策略为 双端测序。生物信息学分析 数据处理对测序 进行质控，获得高质量的。通过 （）进行长度过滤和去噪，以 相似度进行单元聚类获得 （）特征序列和丰度表格。多样性分析基于得到的 特征序列和丰度表格进行 多样性分析。用 指数和 指数评估群落中所包含物种的数目；用 指数和 指数来评估群落的多样性。评判标准为 指数和 指数越大，物种数目越多；指数越大，多样性越高；指数越小，群落多样性越高。分析主坐标分析 是基于距离矩阵，通过一系列的特征值和特征向量排序后，以可视化二维空间重新排列样品，最大化展示样品之间的关系信息。分析中样品点之间的距离越近，说明样品之间物种组成结构越相似。功能预测分析运用 软件与 数据库进行物种组成信息比对，计算相应 的信息及丰度，得到相应代谢通路，在 置信度区间内对各实验组和对照组进行 和 检验，分析其在功能上的差异。渔业研究第 卷 结果与分析 水质指标变化实验期间水温为 ，其中实验组水体芽孢杆菌含量维持在（）之间、值为 ，而对照组 值为 。实验开始后，实验组水体氨氮含量持续上升，第 天达到 ，而后持续下降，至第 天后稳定在 左右；对照组氨氮含量较为稳定，维持在 之间。自第 天起，实验组氨氮含量比对照组稳定下降了 左右（图）。实验组水体亚硝酸盐含量一开始持续上升，第 天达到峰值 ，而后持续下降，至第 天后稳定在 左右；对照组亚硝酸盐含量持续上升，至 后稳定在 左右。自第 天起，实验组亚硝酸盐含量比对照组稳定下降了 左右（图）。图1实验期间水体氨氮（a）、亚硝酸盐（b）含量变化Fig.1The content variation of ammonia nitrogen(a)and nitrite(b)in the test duration 测序样本序列信息测序 个样本下机数据量为 ，平均测序数据量为 。数据共产生 条 、，每个样本平均有 条 、，数据有效性为，序列分布在 之间，平均 含量为 。测序样本覆盖率均大于，说明测序数据量足以代表样品的细菌群落。微生物群落 多样性分析实验组黑脊倒刺鲃前肠、半刺厚唇鱼后肠 指数上升，上升比例分别为 和；指数上升，上升比例分别为 和。其他鱼体肠道部位 和 指数均下降，下降比例分别为 和。黑脊倒刺鲃中肠、半刺厚唇鱼后肠 指数上升，上升比例分别为 和；指数上升，上升比例分别为 和。其他鱼体肠道部位 和 指数均下降，下降比例分别为 和 （图）。第 期张晶晶等：水水体体投投放放枯枯草草芽芽孢孢杆杆菌菌对对不不同同习习性性鱼鱼类类肠肠道道菌菌群群和和水水质质影影响响的的研研究究整体而言，种鱼类各肠道部位实验组 和 指 数 下 降，和 指数下降，细菌丰富度降低。分析黑脊倒刺鲃、半刺厚唇鱼和鲢 分析最主要特征值 和 对样本的差异的贡献值分别合计达到、和，可以代表大部分变量信息。除黑脊倒刺鲃前肠和鲢后肠外，种鱼类其他不同肠道部位中，实验组与对照组各自区分聚类，可用于后续信息分析（图）。基于门水平的微生物群落结构分析如图 所示，黑脊倒刺鲃前肠对照组和实验组的优势菌门为变形菌门、浮霉菌门和放线菌门，梭杆菌门和厚壁菌门分别是两组区别优势菌门，优 势 菌 门 合 计 占 比 分 别 为 和。中肠对照组和实验组的优势菌门为变形菌门和厚壁菌门，浮霉菌门、梭杆菌门和放线菌门、拟杆菌门分别是两组区别优势菌门，优势菌门合计占比分别为 和 。后肠实验 组 的 优 势 菌 门 为 变 形 菌 门，占 比 高 达，对照组共有优势菌门为变形菌门，区别优势菌门为浮霉菌门、放线菌门和拟杆菌门，优势菌门合计占比为 。与对照组相比，浮霉菌门和梭杆菌门在实验组前、中、后肠的丰度均下降；变形菌门在前、后肠的丰度上升，在中肠下降；放线菌门和厚壁菌门在前、中肠的丰度上升，在后肠下降。半刺厚唇鱼前肠实验组、对照组和中肠实验组、对照组的优势菌门为变形菌门、浮霉菌门、放线菌门和蓝细菌门，合计占比分别为 、和 、。后肠对照组优势图2三种鱼体不同肠道部位微生物群落Alpha多样性变化Fig.2Changes of Alpha diversity of microbial community in different intestinal sites of three fish species组别Groups组别Groups组别Groups组别Groups注：黑脊倒刺鲃前肠；黑脊倒刺鲃中肠；黑脊倒刺鲃后肠；半刺厚唇鱼前肠；半刺厚唇鱼中肠；半刺厚唇鱼后肠；鲢前肠；鲢中肠；鲢后肠。以下同此。：；渔业研究第 卷图3三种鱼类前、中、后肠道菌群PCoA分析Fig.3PCoA analysis of gut microbiota in the foregut,midgut and hindgut of the three fishes组别Group组别Group组别GroupPCoA plot(P=0.008)PCoA1(66.24%)PCoA2(16.11%)PCoA2(26.01%)PCoA2(28.23%)PCoA plot(P=0.001)PCoA1(54.79%)PCoA plot(P=0.014)PCoA1(53.93%)0.20.10.0-0.1-0.20.100.050.00-0.05-0.100.2-0.2-0.20.10.0-0.10.10.0-0.10.20.00.00.1-0.1注：组别编号最后一个字母 表示实验组，表示对照组。以下同此。：，图4三种鱼类肠道菌群门水平差异分析Fig.4Variation analysis of gut microbiota at phylum level among three fish species相对丰度/%Relative abundance0255075100Sc_F31CAh_H31EAh_H31CAh_M31EAh_M31CAh_F31EAh_F31CHm_H31EHm_H31CHlm_M31EHm_M31CHm_F31EHm_F31CSc_H31ESc_H31CSc_M31ESc_M31CSc_F31E组别GroupsPhylum变形菌门Protcobactcria浮霉菌门Planctomycetes放线菌门Actinobactcria蓝细菌门Cyanobactcria厚壁菌门Firmicutcs拟杆菌门Bactcroidctes梭杆菌门Fusobacteria疣微菌门Vemcomicrobia绿弯菌门Chloroflex未分类门UnclassifiedCandidatus_Saccharibactcria酸杆菌门Acidobactcria衣原体门Chlamydiae芽单胞菌门GemmatimonadctesCandidatus_Melainabactcria硝化螺旋菌门Nitrospirace绿菌门Chlorobi装甲菌门ArmatimonadetcsCandidatus_Latescibacteria栖热菌门Deinococcus-Thermus其他Others菌门为变形菌门和梭杆菌门，合计占比 ，实验组共有优势菌门为变形菌门，区别优势菌门为梭杆菌门、放线菌门和厚壁菌门，合计占比为。实验组放线菌门、蓝细菌门、厚壁菌门、梭杆菌门在前、中、后肠的丰度均上升，变形菌门在前肠的丰度上升，在中、后肠的丰度下降；浮霉菌门在前、中肠的丰度上升，在后肠的丰度下降。鲢前、中、后肠道实验组和对照组的优势菌门一致。前肠实验组、对照组和中肠实验组、对照组的优势菌群均为变形菌门、浮霉菌门、放线菌门和蓝细菌门，合计占比分别为 、和 、。鲢后肠实验组和对照组的优势菌群为变形菌门、浮霉菌门、放线菌门和拟杆菌门，合计占比分别为 和。放线菌门在前、中、后肠道的丰度上升，蓝细菌门、浮霉菌门在前、中、后肠道的丰度下降；变形菌门在前、中肠的丰度上升，在后肠的丰度下降；拟杆菌门在中、后肠的丰度上升，在前肠的丰度下降。第 期张晶晶等：水水体体投投放放枯枯草草芽芽孢孢杆杆菌菌对对不不同同习习性性鱼鱼类类肠肠道道菌菌群群和和水水质质影影响响的的研研究究 种鱼类共有的优势菌门均为变形菌门，其他优势菌门在各组中均有一定的差异。整体而言，实验组放线菌门在 种鱼类各个肠道部位的丰度均上升，浮霉菌门在 种鱼类各个肠道部位的丰度均下降。种鱼的前肠变形菌门和放线菌门丰度均上升，中肠道放线菌门丰度均上升，后肠道浮霉菌门丰度均下降。基于属水平的微生物群落结构分析如图 所示，黑脊倒刺鲃前肠对照组和实验组共有优势菌属为浮霉菌门未分类属和蓝细菌门未分类属，鲸杆菌属和 属是各自区别优势菌属，各优势菌属合计占比分别为 和 。气单胞属为中肠对照组第一优势菌属，占比为 ，其他优势菌属为鲸杆菌属、厚壁菌门未分类属；实验组优势菌属为厚壁菌门未分类属、浮霉菌门未分类属和芽孢杆菌属，合计占比为 。后肠对照组优势菌属为 变形菌纲未分类属、鲸杆菌属、浮霉菌门未分类属，合计占比为 ；实验组优势菌属主要为气单胞属，占比为 。与对照组相比，实验组浮霉菌门未分类属、鲸杆菌属、弧菌属、出芽菌属在前、中、后肠的丰度均下降，芽孢杆菌属在前、中、后肠的丰度均上升。蓝细菌门未分类属和 变形菌纲未分类属是半刺厚唇鱼前肠对照组和实验组的优势菌属，浮霉菌门未分类属和放线菌门未分类属分别是两组的差别优势菌属，各优势菌属合计占比分别为 和 。蓝细菌门未分类属是中肠对照组和实验组的共有优势菌属，浮霉菌门未分类属、根瘤菌目未分类属和 变形菌纲未分类属、属分别是两组的差别优势菌属，合计占比分别为 和 。属是后肠对照组第一优势菌属，占比为 ，第二、三优势菌属为弧菌属、气 单 胞 属，占 比 分 别 为 和 ；后肠实验组优势菌属为鲸杆菌属、厚壁杆菌门未分类属和气单胞属，合计占比为 。蓝细菌门未分类属、属、芽孢杆菌属在前、中、后肠道的丰度均上升，气单胞属、弧菌属和 属在前、中、后肠道的丰度均下降。图5三种鱼类肠道菌群属水平差异分析Fig.5Variation analysis of gut microbiota at genus level among three fish species气单胞属Aeromonas蓝细菌门未分类属Cyanobacteria_unclassified浮霉菌门未分类属Planctomycetales_unclassified-变形菌纲未分类属Alphaproteobacteria_unclassified放线菌门未分类属Actinobacteria_unclassifiedDemequina-变形菌纲未分类属Betaproteobacteria_unclassified鲸杆菌属Cetobacterium厚壁菌门未分类属Firmicutes_unclassifiedLawsonia根瘤菌目未分类属Rhizobiales_unclassified拟杆菌门未分类属Bacteroidetes_unclassified放线菌目未分类属Actinomycetales_unclassified未分类属Unclassified疣微菌门未分类属Verrucomicrobia_unclassified生丝微菌属Hyphomicrobium浮霉菌科未分类属Planctomycetaceae_unclassified出芽菌属Gemmata芽孢杆菌属Bacillus弧菌属Vibrio其他Others相对丰度/%Relative abundance0255075100Sc_F31CAh_H31EAh_H31CAh_M31EAh_M31CAh_F31EAh_F31CHm_H31EHm_H31CHlm_M31EHm_M31CHm_F31EHm_F31CSc_H31ESc_H31CSc_M31ESc_M31CSc_F31E属水平Genus组别Groups 鲢前肠对照组和实验组共有优势菌属为蓝细菌门未分类属，浮霉菌门未分类属、放线菌门未分类属和 变形菌纲未分类属、属分别是两组的区别优势菌属，合计占比分别为 和 。鲢中肠对照组和实验组共有优势菌属为蓝细菌门未分类属，浮霉菌门未分类属、变形菌纲未分类属和放线菌门未分类属和 属是两组的区别优势菌属，合计渔业研究第 卷占比分别为 和 。鲢后肠对照组和实验组共有优势菌属为 变形菌纲未分类属，蓝细菌门未分类属、浮霉菌门未分类属和放线菌门未分类属、属是两组的区别优势菌属，合计占比分别为 和 。属和芽孢杆菌属在前、中、后肠的丰度均上升，气单胞属、蓝细菌门未分类属、浮霉菌门未分类属和弧菌属在前、中、后肠的丰度均下降。种鱼类肠道菌群以及同种鱼类不同肠道在属水平优势菌群差异较大，各肠道样本菌属比较丰富，且丰度较为均一。整体而言，气单胞属、弧菌属在 种鱼类各肠道部位的丰度均下降，芽孢杆菌属在 种鱼类各肠道部位的丰度均上升。三种鱼类微生物菌落功能差异分析黑脊倒刺鲃中肠外源性物质生物降解与代谢和后肠酶家族、代谢、碳水化合物代谢以及氨基酸代谢提高，前肠多糖生物合成与代谢、后肠外源性物质生物降解与代谢、萜类和多酮类化合物代谢、酯类代谢降低。半刺厚唇鱼前肠和中肠酶家族降低，中肠的其他氨基酸代谢，后肠脂类代谢、萜类和多酮类化合物代谢提高。鲢前肠碳水化合物代谢，中肠其他氨基酸代谢，后肠其他次生代谢物生物合成、萜类和多酮类化合物代谢提高（表）。种鱼类实验组非代谢信号途径中，黑脊倒刺鲃后肠信号转导、细胞过程和信号转导以及细胞运动提高；半刺厚唇鱼前肠膜运输信号通路提高，后肠细胞运动和遗传信息处理水平降低；鲢前肠遗传信息处理和转录水平降低。讨论枯草芽孢杆菌影响水体环境和鱼体肠道需要一定的时间，一般为 周，。因本研究中使用的 种鱼类习性不同，为使枯草芽孢杆菌在 种鱼体中能够稳定定植，研究周期延长至，结果显示投放枯草芽孢杆菌株后，水体 变化不大，这与杭小英等、刘慧玲等的研究报道一致。国内研究大量报道了枯草芽孢杆菌能够降低水体氨氮含量、亚硝酸盐含量，起到净化水质的作用，主要是因为枯草芽孢杆菌可通过自养反硝化和自养亚硝化而降低水中氨氮、亚硝酸态氮的含量。本研究初期投放枯草芽孢杆菌时养殖水体氨氮含量上升，可能是由于枯草芽孢杆菌进入水体，消耗了大量的有机氮，从而造成氨氮短时间的积累。水体中保持一定浓度的枯草芽孢杆菌后，种鱼类实验组各肠道部位菌群整体种类减少，丰富度下降。研究证明，水产动物成鱼的肠道菌群 多样性一般低于幼鱼，低水平的 多样性有利于肠道菌群的稳定。枯草芽孢杆菌通过与致病菌竞争生存空间、定居部位以及营养物质等，可以有效地抑制有害菌的生长和繁殖，改善鱼类肠道内微生物群落的组成。同时，枯草芽孢杆菌可以分泌产生酸性物质、细菌素等抑制或杀灭有害菌，调整肠道微生态平衡。有研究报道，日粮中添加该菌可以减少草鱼（）肠道菌群中假单胞菌和和黄杆菌等条件致病菌的增殖；水体中投放该菌可以抑制克氏原鳌虾（）大肠菌群和金黄色葡萄球菌等有害致病菌的增殖。本研究水体中的枯草芽孢杆菌抑制了 种鱼类中肠道菌群气单胞属、弧菌属的增殖。由此可见，枯草芽孢杆菌作为水体或饲料添加剂，均能改善鱼体肠道菌群，抑制有害病原菌的增殖。本研究中气单胞属在黑脊倒刺鲃实验组大量繁殖的原因有待进一步研究。有报道称气单胞属也是黄颡鱼（）的优势菌属，其可能对黄颡鱼的免疫系统发育起到一定的促进作用。芽孢杆菌作为水产常用的微生态制剂，可以合成益菌物质，抑制细菌、真菌和病毒生长，同时形成生物大分子降解酶，促进鱼体消化。本研究中芽孢杆菌属在 种鱼类各肠道部位的丰度均上升，说明枯草芽孢杆菌可以促进肠道有益菌的生长与定植。大量文献 报道枯草芽孢杆菌可以增强鱼类代谢功能、提高脂肪利用率、促进鱼类生长。使用枯草芽孢杆菌投喂杂交鲟 西伯利亚鲟 施氏鲟 ，可显著降低饵料系数，增强鱼体对血糖和脂肪的利用和转化能力；投喂草鱼枯草芽孢杆菌，其肝脏油及脂肪含量显著下降，脂蛋白脂酶及脂肪甘油三酯脂肪酶的表达水平上调，跟脂质代谢通路相关的基因表达几乎上调。本研究 第 期张晶晶等：水水体体投投放放枯枯草草芽芽孢孢杆杆菌菌对对不不同同习习性性鱼鱼类类肠肠道道菌菌群群和和水水质质影影响响的的研研究究中枯草芽孢杆菌加强了鲢的肠道菌群代谢功能，对其他 种鱼类则效果不佳，产生这一差异可能是因为鲢是滤食性鱼类，具海绵状的鳃耙和较长的肠道，能够有效过滤吸收细菌。胡保同在鲢肠道前端找到被黏液包裹的细菌凝聚体；朱学宝以细菌絮凝体作为直接饵料饲养鲢，后鲢鱼增重 。在自然湖泊中，鲢摄取的细菌数量可占食物总量的。因此相较于半刺厚唇鱼和黑脊倒刺鲃，鲢更容易摄食和吸收水体中的枯草芽孢杆菌。表 三种鱼类不同肠道部位菌群 信号通路相对丰度差异 肠道部位 途径相对丰度 实验组 对照组 黑脊倒刺鲃 前肠中肠后肠多糖生物合成与代谢 外源性物质生物降解与代谢 外源性物质生物降解与代谢 信号转导 萜类和多酮类化合物代谢 代谢 脂类代谢 酶家族 细胞过程和信号转导 细胞运动 碳水化合物代谢 氨基酸代谢 半刺厚唇鱼 前肠中肠后肠酶家族 折叠、分类和降解 膜运输 酶家族 其他氨基酸代谢 细胞运动 遗传信息处理 脂类代谢 萜类和多酮类化合物代谢 鲢 前肠中肠后肠碳水化合物代谢 遗传信息处理 转录 其他氨基酸代谢 其他次生代谢物生物合成 萜类和多酮类化合物代谢 渔业研究第 卷 国内外关于水体投放枯草芽孢杆菌对鱼类肠道菌群影响的报道较少，本研究选取生活于不同水层、不同食性的鱼类作为研究对象，探讨该菌对 种鱼类不同肠道部位菌群和代谢的变化。然而由于样本取样数目的限制导致相关统计学分析无法展开，因此本研究只能基于已有的数据，分析不同肠道门、属和功能的差异变化。由于不同鱼类以及同一鱼类不同肠道菌群差异较大，本文只能简单整合、比较同一鱼类和同一肠道部位共有的菌群变化，以期为后续相关水体投放枯草芽孢杆菌对鱼类肠道菌群影响的研究提供一定的数据支撑和借鉴价值。参考文献：杨明容，陈再炜，曹岩，等 枯草芽孢杆菌对泥鳅养殖水体水质的影响 江西水产科技，（）：李松 外源性枯草芽孢杆菌在水体和斑马鱼肠道的定植情况及其对所处环境菌群的影响 杭州：浙江大学，：，（），（）：邱月 枯草芽孢杆菌在现代农业中的应用 园艺与种苗，（）：，（）：，：宋雪飞，周庆，张迎颖，等 四种非营养性饲料添加剂对澳洲淡水龙虾消化酶活力、免疫功能及离体消化率的影响 渔业研究，（）：王成强，王际英，黄炳山，等 珍珠龙胆石斑鱼肠道枯草芽孢杆菌的分离鉴定及产酶能力分析 渔业研究，（）：，：陆家昌，李活，黄翔鹄 枯草芽孢杆菌对水质及凡纳滨对虾幼体免疫指标影响的研究 南方水产，（）：，（），：，：，：，：，（）：杭小英，叶雪平，施伟达，等 枯草芽孢杆菌制剂对罗氏沼虾养殖池塘水质的影响 浙江海洋学院学报（自然科学版），（）：刘慧玲，黄翔鹄，李长玲，等不同浓度的枯草芽孢杆菌对罗非鱼鱼苗的养殖水体水质及其抗病力的影响 水产养殖，（）：王华，王京伟，杜刚养殖水体中微生物全程自养脱氮初步研究 水产科学，（）：，第 期张晶晶等：水水体体投投放放枯枯草草芽芽孢孢杆杆菌菌对对不不同同习习性性鱼鱼类类肠肠道道菌菌群群和和水水质质影影响响的的研研究究 ，（）：，：，：曹树威，余昌花，罗鲜青，等 芽孢杆菌特性及其在畜牧生态养殖中应用研究进展 中国兽药志，（）：，（），（）：彭春雷 枯草芽孢杆菌抑制水体微生物及对克氏原螯虾生存状态的影响 荆楚理工学院学报，（）：，：张书环，吴金平，褚志鹏，等 饲喂枯草芽孢杆菌对杂交鲟生长和肠道菌群结构的影响 海洋渔业，（）：罗燕儿，赵慧，郭道远，等 枯草芽孢杆菌对草鱼肝脏脂质代谢的调节作用 水生生物学报，（）：，（），（）：肖霞，赵娟娟 养殖滤食性鱼类对大水面水体净化作用的探讨 河南水产，（）：胡保同 滤食性鱼类摄食机制问题 水产科技情报，（）：朱学宝 细菌絮凝体对滤食性鱼类饵料效果的研究 水产学报，（）：渔业研究第 卷 ，（，）：，：；