1 材料与方法
1.1 试验菌株
1.2 试验饲料
表1 基础饲料组成及营养水平(风干基础)Table 1 Composition and nutrient levels of the basal diet (air-dry basis) % |
| 项目Items | 含量Contents |
|---|---|
| 原料Ingredients | |
| 鱼粉Fish meal | 52.00 |
| 酪蛋白Casein | 11.98 |
| 明胶Gelatin | 3.00 |
| 玉米淀粉Corn starch | 17.73 |
| 鱼油Fish oil | 0.82 |
| 豆油Soy oil | 3.50 |
| 大豆磷脂Soybean lecithin | 2.00 |
| 磷酸二氢钙Ca(H2PO4)2 | 1.50 |
| 氯化胆碱Choline chloride | 0.40 |
| 稳定型维生素C Stable vitamin C | 0.02 |
| 维生素预混料Vitamin premix1) | 0.40 |
| 矿物质预混料Mineral premix2) | 0.50 |
| 微晶纤维素Microcrystalline cellulose | 6.15 |
| 合计Total | 100.00 |
| 营养水平Nutrient levels3) | |
| 干物质Dry matter | 91.86 |
| 粗蛋白质Crude protein | 49.90 |
| 粗脂肪Crude lipid | 13.36 |
| 粗灰分Ash | 11.63 |
1)维生素预混料为每千克饲料提供 Vitamin premix provided the following per kg of the diet:VA 10 mg,VD 10 mg,VE 100 mg,VK3 10 mg 10 mg,VB1 10 mg,VB2 20 mg,VB6 20 mg,VB12 0.05 mg,尼克酸 nicotinic acid 50 mg,泛酸钙 calcium pantothenate 100 mg,生物素 biotin 1 mg,肌醇 inositol 500 mg,叶酸 folic acid 4 mg。 2)矿物质预混料为每千克饲料提供 Mineral premix provided the following per kg of the diet:柠檬酸铁 ferric citrate 497 mg,CuSO4·5H2O 24 mg,ZnSO4·7H2O 176 mg, MnSO4·4H2O 122 mg,CoCl2·6H2O 0.18 mg, KIO3 0.51 mg,Na2SeO3 0.33 mg。 3)营养水平为实测值。Nutrient levels were measured values. |
1.3 试验设计及饲养管理
1.4 样品采集与制备
1.5 指标测定
1.5.1 生长性能计算
1.5.2 肠道抗氧化指标测定
1.5.3 肠道短链脂肪酸含量测定
1.6 肠道组织切片制作
1.7 肠道菌群测定
1.8 肠道代谢物测定
1.8.1 样品选择和代谢物提取
1.8.2 色谱-质谱分析
1.8.3 代谢物指纹图谱信息分析
1.9 数据统计分析
2 结果
2.1 生长性能
表2 丁酸梭菌对珍珠龙胆石斑鱼生长性能的影响Table 2 Effects of Clostridium butyricum on growth performance of hybrid grouper |
| 项目 Items | 组别Groups | ||
|---|---|---|---|
| CON | CB1 | CB2 | |
| 初体重IBW/(g/尾) | 24.89±0.05 | 24.93±0.04 | 24.85±0.03 |
| 末体重FBW/(g/尾) | 126.36±5.88 | 127.82±2.79 | 128.02±6.17 |
| 增重率WGR/% | 407.67±22.65 | 412.79±10.5 | 415.19±25.53 |
| 特定生长率SGR/(%/d) | 2.90±0.08 | 2.92±0.04 | 2.93±0.09 |
| 饲料效率FE/% | 110.82±3.04 | 112.10±0.71 | 109.16±1.46 |
| 日摄食率DFR/(%/d) | 2.15±0.08 | 2.13±0.01 | 2.18±0.01 |
| 肝体比HSI/% | 2.39±0.19 | 2.65±0.19 | 2.44±0.30 |
| 肥满度CF/(g/cm3) | 2.97±0.04 | 2.86±0.05 | 2.85±0.05 |
同行数据肩标无字母或相同字母表示差异不显著(P>0.05),不同小写字母表示差异显著(P<0.05)。下表同。 | |
In the same row, values with no letter or the same letter superscripts mean no significant difference (P>0.05), while with different small letter superscripts mean significant difference (P<0.05). The same as below. |
2.2 肠道生化指标
2.2.1 肠道抗氧化指标
表3 丁酸梭菌对珍珠龙胆石斑鱼肠道抗氧化指标的影响Table 3 Effects of Clostridium butyricum on intestinal antioxidant indexes of hybrid grouper |
| 项目 Items | 组别Groups | ||
|---|---|---|---|
| CON | CB1 | CB2 | |
| 总抗氧化能力 Total antioxidant capacity/(mmol/g prot) | 1.18±0.27 | 1.73±0.27 | 1.89±0.14 |
| 谷胱甘肽过氧化物酶 Glutathione peroxidase/(U/mg prot) | 502.29±24.34a | 943.57±45.56b | 851.37±23.90b |
| 过氧化氢酶 Catalase/(U/mg prot) | 21.83±0.08 | 20.53±0.67 | 19.66±0.32 |
| 超氧化物歧化酶 Superoxide dismutase/(U/mg prot) | 26.60±0.49a | 32.67±1.44b | 32.69±1.76b |
| 丙二醛 Malondialdehyde/(nmol/mg prot) | 3.22±0.73b | 3.93±0.01b | 1.15±0.36a |
2.2.2 肠道短链脂肪酸含量
表4 丁酸梭菌对珍珠龙胆石斑鱼肠道短链脂肪酸含量的影响Table 4 Effects of Clostridium butyricum on intestinal SCFA contents of hybrid grouper μmol/g |
| 项目 Items | 组别Groups | ||
|---|---|---|---|
| CON | CB1 | CB2 | |
| 乙酸Acetic acid | 49.03±0.81a | 59.82±2.89b | 59.39±3.11b |
| 丙酸Propionic acid | 2.21±0.18 | 1.68±0.19 | 2.60±0.89 |
| 丁酸Butyric acid | 0.62±0.02 | 0.70±0.05 | 0.72±0.001 |
| 异丁酸Isobutyric acid | 0.61±0.02 | 0.70±0.05 | 0.81±0.16 |
| 戊酸Valeric acid | 0.43±0.04 | 0.42±0.03 | 0.43±0.05 |
| 异戊酸Isovaleric acid | 1.67±0.25 | 1.93±0.05 | 2.54±0.80 |
| 己酸Caproic acid | 1.18±0.16 | 0.83±0.02 | 1.14±0.26 |
| 总短链脂肪酸Total SCFA | 55.79±0.64a | 65.88±2.99b | 67.58±1.02b |
2.3 肠道组织结构
表5 丁酸梭菌对珍珠龙胆石斑鱼肠道组织形态指标的影响Table 5 Effects of Clostridium butyricum on intestinal morphology indexes of hybrid grouper |
| 项目 Items | 组别Groups | |||
|---|---|---|---|---|
| CON | CB1 | CB2 | ||
| 前肠 Proximal intestine | 肌层厚度Muscle thickness/μm | 189.82±43.87 | 175.15±32.39 | 178.95±12.15 |
| 黏膜褶皱Number of fold mucosa/个 | 61.50±6.50 | 69.33±5.24 | 73.00±4.00 | |
| 黏膜褶皱高度Height of mucosal folds/μm | 459.73±36.73 | 447.43±14.36 | 521.06±7.69 | |
| 中肠 Middle intestine | 肌层厚度Muscle thickness/μm | 104.06±12.43a | 105.92±3.46a | 138.70±9.91b |
| 黏膜褶皱Number of fold mucosa/个 | 33.50±0.50 | 36.50±0.50 | 35.42±1.46 | |
| 黏膜褶皱高度Height of mucosal folds/μm | 351.19±29.99a | 391.89±14.32ab | 444.23±12.15b | |
| 后肠 Distal intestine | 肌层厚度Muscle thickness/μm | 224.40±94.44 | 187.28±36.39 | 214.97±40.56 |
| 黏膜褶皱Number of fold mucosa/个 | 51.50±6.50 | 51.67±6.49 | 52.00±6.00 | |
| 黏膜褶皱高度Height of mucosal folds/μm | 465.42±90.64 | 413.36±33.93 | 407.22±16.16 | |
图1 丁酸梭菌对珍珠龙胆石斑鱼肠道组织形态的影响PI:前肠 proximal intestine:MI:中肠middle intestine;DI:后肠 distal intestine;MT:肌层厚度 muscle thickness;LP:固有层 lamina propria;GC:杯状细胞 goblet cells。 A、B、C分别表示CON、CB1、CB2组;肌层厚度由黑色双向箭头表示;固有层由黑色单向箭头表示;杯状细胞由蓝色单向箭头表示。A, B, and C represent the CON, CB1, and CB2 groups, respectively. The black solid bidirectional arrows indicate the muscle thickness, the black solid unidirectional arrows indicate the lamina propria, and the blue solid unidirectional arrows indicate goblet cells. Fig.1 Effects of Clostridium butyricum on intestinal morphology of hybrid grouper (100×) |
2.4 肠道菌群变化
图2 丁酸梭菌对珍珠龙胆石斑鱼肠道菌群在门水平上的相对丰度的影响Proteobacteria:变形菌门;Firmicutes:厚壁菌门;Bacteroidota:拟杆菌门;Actinobacteriota:放线菌门;Chloroflexi:绿弯菌门;Campilobacterota:弯曲杆菌门;Spirochaetota:螺旋体门;Desulfobacterota:脱硫杆菌门;Other:其他。 Fig.2 Effects of Clostridium butyricum on relative abundance of gut microbiota at phylum level of hybrid grouper |
图3 丁酸梭菌对珍珠龙胆石斑鱼肠道菌群在属水平上的相对丰度的影响(前10)Pseudomonas:假单胞菌属;Lactobacillus:乳杆菌属;Massilia:马赛菌属;Brevundimonas:短波单胞菌属;Bacteroides:拟杆菌属;Subdoligranulum:罕见小球菌属;Eubacterium_coprostanoligenes_group:真杆菌属群;Other:其他。 Fig.3 Effects of Clostridium butyricum on relative abundance of gut microbiota at genus level of hybrid grouper (top 10) |
2.5 代谢组学分析
2.5.1 PCA
2.5.2 OPLS-DA模型
2.5.3 差异代谢物筛选
图6 差异代谢物筛选火山图Differential metabolites:差异代谢物;down regulated:下调;not significant:不显著;up regulated:上调;VIP:变量投影重要性 variable importance in projection;-log10 (P-value):-log10(P值);log2 (fold change):对数变化倍数 logarithm base 2 of fold change。 Fig.6 Volcano diagram for screening differential metabolites |
2.5.4 KEGG代谢通路富集分析
图7 肠道差异代谢物的KEGG代谢通路图(前20)Hits_number:命中数 number of hits;Raw_P:未校正P值 unadjusted P-values;Impact:影响;Insulin signaling pathway:胰岛素信号通路;Glycosylphosphatidylinostol (GPI)-anchor biosynthesis:糖基磷脂酰肌醇(GPI)锚生物合成;Apoptosis:细胞凋亡;Regulation of actin cytoskeleton:肌动蛋白骨架的调控;PPAR signaling pathway:PPAR信号通路;Necroptosis:坏死性凋亡;FoxO signaling pathway:FoxO信号通路;Linoleic acid metabolism:亚油酸代谢;Adipocytokine signaling pathway:脂肪细胞因子信号通路;Starch and sucrose metabolism:淀粉和蔗糖代谢;Vascular smooth muscle contraction:血管平滑肌收缩;Galactose metabolism:半乳糖代谢;Sphingolipid metabolism:鞘脂类代谢;ABC transporters:ABC转运蛋白;Biosynthesis of unsaturated fatty acids:不饱和脂肪酸生物合成;Pentose phosphate pathway:磷酸戊糖途径;Glycerophospholipid metabolism:甘油磷脂代谢;Nucleotide metabolism:核苷酸代谢;Fatty acid biosynthesis:脂肪酸合成;Arachidonic acid metabolism:花生四烯酸代谢。 Fig.7 KEGG metabolic pathway diagram of intestinal differential metabolites (top 20) |
图8 肠道差异代谢通路的差异丰度得分图(前20)Hits:命中;DA:差异丰度 differential abundance;DA score:差异丰度得分 differential abundance score;Starch and sucrose metabolism:淀粉和蔗糖代谢;Sphingolipid metabolism:鞘脂类代谢;Regulation of actin cytoskeleton:肌动蛋白骨架的调控;Pentose phosphate pathway:磷酸戊糖途径;Necroptosis:坏死性凋亡;Insulin signaling pathway:胰岛素信号通路;Glycosylphosphatidylinostol (GPI)-anchor biosynthesis:糖基磷脂酰肌醇(GPI)锚生物合成;Glycerophospholipid metabolism:甘油磷脂代谢;Galactose metabolism:半乳糖代谢;FoxO signaling pathway:FoxO信号通路;Apoptosis:细胞凋亡;Adipocytokine signaling pathway:脂肪细胞因子信号通路;ABC transporters:ABC转运蛋白;Nucleotide metabolism:核苷酸代谢;Vascular smooth muscle contraction:血管平滑肌收缩;Arachidonic acid metabolism:花生四烯酸代谢;Linoleic acid metabolism:亚油酸代谢;PPAR signaling pathway:PPAR信号通路;Fatty acid biosynthesis:脂肪酸合成;Biosynthesis of unsaturated fatty acids:不饱和脂肪酸生物合成。 Fig.8 Differential abundance score diagram of intestinal differential metabolic pathways (top 20) |
表6 显著富集通路中的差异代谢物的变化Table 6 Changes in differential metabolites within significantly enriched pathways |
| 项目 Items | 对照组 CON group/ (mg/g) | CB2组 CB2 group/ (mg/g) | Log2 | P值 P-value | 变量投影 重要性 VIP | 变化趋势 Trends of change | |||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 蜜二糖Melibiose | 0.38±0.07 | 0.73±0.13 | 0.92 | 0.05 | 1.24 | Up | |||||||
| 棉子糖Raffinose | 0.06±0.01 | 0.22±0.04 | 1.91 | 0.00 | 1.67 | Up | |||||||
| 麦芽三糖Maltotriose | 0.05±0.01 | 0.18±0.03 | 1.84 | 0.00 | 1.76 | Up | |||||||
| 海藻糖Trehalose | 0.53±0.010 | 1.47±0.14 | 1.47 | 0.00 | 1.83 | Up | |||||||
| 葡萄糖Glucose | 7.18±0.65 | 13.11±1.38 | 0.87 | 0.00 | 1.65 | Up | |||||||
| 蔗糖Sucrose | 0.53±0.10 | 1.47±0.14 | 1.47 | 0.00 | 1.83 | Up | |||||||
| 谷胱甘肽Glutathione | 0.15±0.02 | 0.26±0.03 | 0.75 | 0.02 | 1.45 | Up | |||||||
| 腺苷Adenosine | 19.31±1.62 | 32.46±5.27 | 0.75 | 0.04 | 1.27 | Up | |||||||
| 鸟嘌呤核苷Guanosine | 4.97±0.82 | 2.26±0.38 | -1.14 | 0.01 | 1.45 | Down | |||||||
| 异麦芽糖Isomaltose | 0.38±0.01 | 0.73±0.13 | 0.92 | 0.17 | 1.24 | Up | |||||||
| 麦芽糖Maltose | 0.53±0.10 | 1.47±0.14 | 1.47 | 0.00 | 1.83 | Up | |||||||
| 磷脂酰胆碱PC (44:11) | 0.22±0.06 | 0.47±0.02 | 1.09 | 0.00 | 1.71 | Up | |||||||
| γ-亚麻酸γ-linolenic acid | 30.48±5.45 | 13.63±1.14 | -1.16 | 0.01 | 1.46 | Down | |||||||
| 亚油酸酯Linoleate | 1.71±0.40 | 0.69±0.04 | -1.30 | 0.03 | 1.31 | Down | |||||||
| 油酸Oleic acid | 20.38±4.81 | 7.11±0.54 | -1.52 | 0.02 | 1.38 | Down | |||||||
| 十六烷酸Hexadecanoic acid | 0.67±0.17 | 0.11±0.02 | -2.64 | 0.01 | 1.54 | Down | |||||||
| α-亚麻酸α-linolenic acid | 30.48±5.45 | 13.63±1.14 | -1.16 | 0.01 | 1.46 | Down | |||||||
| 13-酮-9z,11e-十八碳二烯酸 13-OxoODE | 2.16±0.20 | 1.43±0.26 | -0.60 | 0.05 | 1.21 | Down | |||||||
| 13(s)-羟基-9z,11e-十八碳二烯酸 13(S)-HODE | 7.58±1.90 | 0.85±0.20 | -3.15 | 0.01 | 1.57 | Down |
Log2为该物质在该组对比的相对丰度比值取log2;P-value为该物质在该组对比的t检验得到的P值;VIP为该物质在OPLS-DA模型中的VIP值;Up为CB2组含量上升;Down为CB2组含量下降。 | |
Log2 is the relative abundance ratio of the substance in the comparison group, taken as log2; P-value is the P-value obtained from the t-test of the substance in the comparison group; VIP is the VIP value of the substance in the OPLS-DA model; Up indicates that the content of CB2 group increased; Down indicates that the content of CB2 group decreases. |
图9 饲料中添加丁酸梭菌对珍珠龙胆石斑鱼肠道代谢影响的代谢物网络图Clostridium btyricum:丁酸梭菌;Direct action:直接作用;Indirect action:间接作用;Up-regulation:上调;Down-regulation:下调;Glucose:葡萄糖;Melibiose:蜜二糖;Raffinose:棉子糖;Galactinol:肌醇半乳糖苷;Sucrose:蔗糖;Levulose:果糖;Galactose metabolism:半乳糖代谢;Hexadecanoic acid:十六烷酸;Palmitoyl-CoA:棕榈酰辅酶A palmitoyl CoA;Oleoyl-CoA:油酰辅酶A oleyl CoA;Oleic acid:油酸;Linoleoyl-CoA:亚油酰基辅酶A linoleyl-CoA;Linoleate:亚油酸脂;Alpha-linoleoyl-CoA:α-亚油酰基辅酶A;Alpha-linoleate:α-亚油酸脂;Gamma-linoleoyl-CoA:γ-亚油酰基辅酶A;Gamma-linolenic acid:γ-亚麻酸;Biosynthesis of unsaturated fatty acids:不饱和脂肪酸的生物合成;PC:磷脂酰胆碱 phosphatidylcholine;Arachidonic acid metabolism:花生四烯酸代谢;13(S)-HPODE:13(S)-羟基-9z,11e-二十碳四烯酸过氧化物 13(S) -hydroxy-9z, 11e-eicosatetraenoic acid peroxide;13(S)-HODE:13(s)-羟基-9z,11e-十八碳二烯酸 13(s) -hydroxy-9z, 11e-octadecadienoic acid;13-OxoODE:13-酮-9z,11e-十八碳二烯酸 13-keto-9z, 11e-octadecadienoic acid;Linoleic acid metabolism:亚油酸代谢;Starch and sucrose metabolism:淀粉和蔗糖代谢;Starch:淀粉;Maltose:麦芽糖;Trehalose:海藻糖;Dextran:葡聚糖;Isomaltose:异麦芽糖;Maltotriose:麦芽三糖;Adenosine:腺苷;Guanosine:鸟嘌呤核苷;Glutathione:谷胱甘肽;ABC transporters:ABC转运蛋白;Oligosaccharide:低聚糖;Polyol:多元醇;Lipid transporters:脂质转运蛋白;Monosaccharide transporters:单糖转运蛋白;Peptide and nickel transporters:肽和镍转运蛋白。 Fig.9 Metabolite network diagram of effects of dietary Clostridium butyricum on intestinal metabolism of hybrid grouper |
