1 材料与方法
1.1 试验材料
1.2 试验设计
1.2.1 育肥试验
表1 基础饲粮组成及营养水平(干物质基础)Table 1 Composition and nutrient levels of experimental diets (DM basis)% |
| 项目 Items | 对照组 Control group | 试验Ⅰ组 Experiment group Ⅰ | 试验Ⅱ组 Experiment group Ⅱ | 试验Ⅲ组 Experiment group Ⅲ |
|---|---|---|---|---|
| 原料 Ingredients | ||||
| 豆粕 Soybean meal | 14.0 | 14.0 | 14.0 | 14.0 |
| 玉米 Corn | 21.0 | 21.0 | 21.0 | 21.0 |
| 花生秧 Peanut seedling | 26.0 | 17.0 | 9.0 | |
| 杂交狼尾草 Hybrid Penisetum | 35.0 | 35.0 | 35.0 | 35.0 |
| 大蒜皮Garlic skin | 8.0 | 16.0 | 24.0 | |
| 食盐 NaCl | 0.5 | 0.5 | 0.5 | 0.5 |
| 石粉 Limestone | 0.4 | 0.6 | 1.0 | 1.5 |
| 磷酸氢钙 CaHPO4 | 0.7 | 1.0 | 0.2 | 0.2 |
| 磷酸脲 Urea phosphate | 0.4 | 0.9 | 1.3 | 1.8 |
| 预混料 Premix1) | 2.0 | 2.0 | 2.0 | 2.0 |
| 合计 Total | 100.0 | 100.0 | 100.0 | 100.0 |
| 营养水平 Nutrient levels2) | ||||
| 干物质(鲜样) DM (fresh) | 34.17 | 34.24 | 34.28 | 34.35 |
| 有机物 OM | 88.15 | 87.84 | 87.89 | 87.63 |
| 总能 GE/(MJ/kg) | 15.35 | 15.35 | 15.45 | 15.37 |
| 代谢能 ME/(MJ/kg) | 9.24 | 9.22 | 9.24 | 9.24 |
| 粗蛋白质 CP | 14.26 | 14.21 | 14.26 | 14.28 |
| 中性洗涤纤维 NDF | 44.42 | 44.64 | 44.47 | 44.47 |
| 酸性洗涤纤维 ADF | 22.11 | 22.14 | 22.12 | 22.10 |
| 钙 Ca | 0.98 | 0.97 | 0.99 | 0.99 |
| 总磷 TP | 0.54 | 0.54 | 0.56 | 0.55 |
1)预混料为每千克饲粮提供 The premix provided the following per kg of diets:Cu 20.0 mg,Fe 80.0 mg,Mn 30.0 mg,Zn 80.0 mg,I 1 mg,Se 0.30 mg,VA 20 000 IU,VD 5 000 IU,VE 50.0 mg。 2)代谢能为计算值,其他均为实测值。ME was a calculated value, while the others were measured values. |
1.2.2 消化代谢试验
1.3 饲养管理
1.4 样品采集
1.5 指标测定
1.5.1 生长性能测定
1.5.2 营养成分和营养物质表观消化率测定
1.5.3 瘤胃发酵参数测定
1.5.4 瘤胃微生物结构测定
1.6 数据统计方法
2 结果
2.1 饲粮中添加不同水平大蒜皮对福清山羊生长性能和营养物质表观消化率的影响
表2 饲粮中添加不同水平大蒜皮对福清山羊生长性能和营养物质表观消化率的影响Table 2 Effects of dietary supplemented with different levels of garlic skin on growth performance and nutrient apparent digestibility of Fuqing goats |
| 项目 Items | 对照组 Control group | 试验Ⅰ组 Experiment group Ⅰ | 试验Ⅱ组 Experiment group Ⅱ | 试验Ⅲ组 Experiment group Ⅲ | SEM | P值 P-value |
|---|---|---|---|---|---|---|
| 生长性能 Growth performance | ||||||
| 干物质采食量 DMI/(g/d) | 640.07 | 643.01 | 642.04 | 645.52 | 1.250 | 0.484 |
| 初始体重 IBW/kg | 15.33 | 15.40 | 15.52 | 15.43 | 0.523 | 0.818 |
| 终末体重 FBW/kg | 23.33 | 23.45 | 23.16 | 23.18 | 0.984 | 0.842 |
| 平均日增重 ADG/(g/d) | 88.90 | 89.41 | 84.95 | 86.19 | 1.432 | 0.652 |
| 营养物质表观消化率 Nutrient apparent digestibility/% | ||||||
| 干物质 DM | 56.25 | 58.64 | 59.14 | 58.94 | 0.450 | 0.054 |
| 有机物 OM | 62.34 | 61.82 | 61.75 | 61.96 | 0.363 | 0.945 |
| 粗蛋白质 CP | 67.85 | 72.69 | 71.48 | 73.03 | 0.809 | 0.075 |
| 总能 GE | 57.18b | 58.67b | 57.19b | 61.98a | 0.596 | 0.017 |
| 中性洗涤纤维 NDF | 51.77b | 54.40ab | 55.09ab | 56.35a | 0.557 | 0.021 |
| 酸性洗涤纤维 ADF | 35.06Bb | 35.07Bb | 37.37Bb | 40.35Aa | 0.628 | 0.007 |
同行数据肩标不同小写字母表示差异显著(P<0.05),不同大写字母表示差异极显著(P<0.01),相同小写字母或无字母表示差异不显著(P>0.05)。下表同。 | |
In the same row, values with different small letter superscripts mean significant difference (P<0.05), and with different capital letter superscripts mean significant difference (P<0.01), while with the same small letter or no letter superscripts mean no significant difference (P>0.05). The same as below. |
2.2 饲粮中添加不同水平大蒜皮对福清山羊瘤胃发酵参数的影响
表3 饲粮中添加不同水平大蒜皮对福清山羊瘤胃发酵参数的影响Table 3 Effects of dietary supplemented with different levels of garlic skin on rumen fermentation parameters of Fuqing goats |
| 项目 Items | 对照组 Control group | 试验Ⅰ组 Experiment group Ⅰ | 试验Ⅱ组 Experiment group Ⅱ | 试验Ⅲ组 Experiment group Ⅲ | SEM | P值 P-value |
|---|---|---|---|---|---|---|
| pH | 6.63 | 6.61 | 6.68 | 6.63 | 0.034 | 0.900 |
| 氨态氮 NH3-N/(mg/dL) | 16.35a | 14.75ab | 13.91b | 13.46b | 0.350 | 0.014 |
| 乙酸 Acetate/(mmoL/L) | 27.38 | 29.68 | 28.08 | 28.92 | 0.612 | 0.588 |
| 丙酸 Propionate/(mmoL/L) | 9.76b | 10.35ab | 11.09a | 11.49a | 0.230 | 0.039 |
| 异丁酸 Isobutyrate/(mmoL/L) | 4.00 | 4.53 | 4.44 | 4.31 | 0.277 | 0.926 |
| 丁酸 Butyrate/(mmoL/L) | 8.72 | 8.29 | 9.04 | 9.35 | 0.236 | 0.470 |
| 异戊酸 Isovalerate/(mmoL/L) | 1.81 | 1.58 | 1.57 | 1.82 | 0.098 | 0.721 |
| 戊酸 Valerate/(mmoL/L) | 2.81 | 2.77 | 2.96 | 2.76 | 0.109 | 0.914 |
| 乙酸/丙酸 Acetate/propionate | 2.82Aa | 2.87Aa | 2.53Bb | 2.52Bb | 0.050 | 0.005 |
| 总挥发性脂肪酸 TVFA/(mmoL/L) | 54.48 | 57.20 | 57.17 | 58.64 | 0.826 | 0.426 |
2.3 饲粮中添加不同水平大蒜皮对福清山羊瘤胃微生物结构的影响
2.3.1 瘤胃微生物Alpha多样性及Beta多样性
表4 饲粮中添加不同水平大蒜皮对福清山羊瘤胃微生物Alpha多样性指数的影响Table 4 Effects of dietary supplemented with different levels of garlic skin on rumen microbiota Alpha diversity indexes of Fuqing goats |
| 项目 Items | 对照组 Control group | 试验Ⅰ组 Experiment group Ⅰ | 试验Ⅱ组 Experiment group Ⅱ | 试验Ⅲ组 Experiment group Ⅲ | SEM | P值 P-value |
|---|---|---|---|---|---|---|
| Ace指数 Ace index | 1 099.425b | 1 329.230a | 1 419.717a | 1 259.012a | 41.981 | 0.031 |
| Chao1指数 Chao1 index | 1 088.841b | 1 316.940a | 1 409.307a | 1 247.844a | 41.681 | 0.029 |
| Simpson指数 Simpson index | 0.989 | 0.993 | 0.994 | 0.991 | 0.001 | 0.150 |
| Shannon指数 Shannon index | 8.190b | 8.635a | 8.695a | 8.356ab | 0.071 | 0.022 |
| PD_whole_tree指数 PD_whole_tree index | 116.352b | 145.142a | 148.702a | 142.773a | 4.899 | 0.045 |
| 覆盖度 Coverage | 0.999 | 0.999 | 0.999 | 0.999 | <0.001 | 0.135 |
2.3.2 瘤胃微生物菌群相对丰度
表5 饲粮中添加不同水平大蒜皮对福清山羊瘤胃微生物门水平和属水平微生物菌群相对丰度的影响Table 5 Effects of dietary supplemented with different levels of garlic skin on relative abundance of rumen microbial flora at phylum level and genus level of Fuqing goats% |
| 项目 Items | 对照组 Control group | 试验Ⅰ组 Experiment group Ⅰ | 试验Ⅱ组 Experiment group Ⅱ | 试验Ⅲ组 Experiment group Ⅲ | SEM | P值 P-value |
|---|---|---|---|---|---|---|
| 门水平 Phylum level | ||||||
| 厚壁菌门 Firmicutes | 49.43b | 50.83b | 44.01b | 61.65a | 2.286 | 0.049 |
| 拟杆菌门 Bacteroidetes | 33.81a | 33.56a | 40.45a | 19.95b | 2.559 | 0.034 |
| 髌骨菌门 Patescibacteria | 2.96 | 4.04 | 3.45 | 4.43 | 0.518 | 0.826 |
| 疣微菌门 Verrucomicrobiota | 3.59 | 2.08 | 3.17 | 4.10 | 0.300 | 0.089 |
| 变形菌门 Proteobacteria | 3.53 | 2.24 | 3.69 | 2.63 | 0.460 | 0.632 |
| 脱硫菌门 Desulfobacterota | 2.64 | 2.09 | 1.20 | 2.47 | 0.326 | 0.414 |
| 互养菌门 Synergistetes | 2.19 | 1.92 | 1.93 | 1.14 | 0.235 | 0.550 |
| 放线菌门 Actinobacteria | 0.32 | 1.19 | 0.20 | 2.30 | 0.353 | 0.164 |
| 蓝藻菌门 Cyanobacteria | 0.70 | 0.91 | 1.02 | 0.24 | 0.147 | 0.299 |
| 未分类细菌 Unclassified_Bacteria | 0.29 | 0.22 | 0.44 | 0.37 | 0.065 | 0.640 |
| 属水平 Genus level | ||||||
| 普雷沃氏菌属 Prevotella | 9.15ab | 13.32a | 14.66a | 4.49b | 1.399 | 0.048 |
| 未培养的瘤胃菌属 Uncultured_rumen_bacterium | 14.79 | 10.21 | 9.42 | 9.26 | 1.209 | 0.363 |
| 未分类毛螺菌科 Unclassified_Lachnospiraceae | 6.26 | 4.75 | 3.72 | 5.62 | 0.453 | 0.209 |
| 毛螺菌科XPB1014群 Lachnospiraceae_XPB1014_group | 3.14b | 3.66b | 4.37b | 9.53a | 0.714 | 0.005 |
| 奎氏菌属 Quinella | 8.98 | 3.95 | 4.66 | 0.19 | 1.091 | 0.055 |
| 理研菌科RC9肠道群 Rikenellaceae_RC9_gut_group | 2.77b | 4.15ab | 7.21a | 1.67b | 0.766 | 0.047 |
| NK4A214群 NK4A214_group | 3.65b | 3.40b | 2.28b | 8.21a | 0.587 | 0.001 |
| 克里斯滕森菌科R_7群 Christensenellaceae_R_7_group | 3.53b | 2.79b | 1.84b | 5.86a | 0.415 | 0.002 |
| 候选单胞生糖菌属 Candidatus_Saccharimonas | 1.92 | 3.12 | 2.86 | 3.99 | 0.423 | 0.487 |
| 未分类F082 Unclassified_F082 | 2.02 | 2.74 | 3.74 | 2.70 | 0.336 | 0.338 |
图2 基于福清山羊瘤胃微生物菌群的随机森林模型分析Lachnospiraceae_XPB1014_group:毛螺菌科XPB1014群;Quinella:奎氏菌属;Sphingobacterium:鞘脂杆菌属;unclassified_Cyanobacteriales:未分类蓝藻菌目;[Anaerorhabdus]_furcosa_group:棍状厌氧菌属;Candidatus_Saccharimonas:糖类假丝酵母菌属;unclassified_Oxalobacteraceae:未分类草酸杆菌科;Lachnoclostridium:乳酸菌属;Lactococcus:乳球菌属;Pseudomonas:假单胞菌属;uncultured_rumen_bacterium_5C0d_16:未培养瘤胃细菌5C0d_16;unclassified_Eggerthellaceae:未分类的伊格尔兹氏菌科;[Eubacterium]_ventriosum_group:凸腹真杆菌属;unclassified_Coriobacteriales:未分类红蝽菌目;Atopobium:奇异菌属;Incertae Sedis:分类未定;Pyramidobacter:锥形杆菌属;Papillibacter:乳头杆菌属;Lachnospiraceae_ND3007_group:毛螺菌科ND3007群;unclassified_Coriobacteriales_Incertae_Sedis:未分类的赤毛杆菌目分类未定;uncultured_Erysipelotrichaceae_bacterium:未培养丹毒丝菌科细菌;Syntrophococcus:互养球菌属;Lachnospiraceae_NK4A136_group:毛螺菌科NK4A136群;unclassified_Diplorickettsiaceae:未分类梁立克星科;unclassified_Erysipelotrichaceae:未分类丹毒丝菌科;Chryseobacterium:金黄杆菌属;Family_ⅩⅢ_AD3011 group:Family_ⅩⅢ_AD3011群;Levilactobacillus:促生乳杆菌属;Mean Decrease Gini:基尼指数。 Fig.2 A random-forest model based on rumen microbial flora of Fuqing goats |
2.4 营养物质表观消化率与瘤胃微生物菌群相关性分析
图4 营养物质表观消化率与菌属和Alpha指数相关性分析网络组合图Mantel’s p:环境因子与菌属和Alpha指数Mantel分析的P值 P-value of Mantel analysis between environmental factors and bacterial genera and Alpha index;Mantel’s r:环境因子与菌属和Alpha指数Mantel分析的r值 r value of Mantel analysis between environmental factors and bacterial genera and Alpha index;Pearson’s r:环境因子与菌属和Alpha指数相关性P值 P-value of correlation between environmental factors and bacterial genera and Alpha index;DM:干物质表观消化率 dry matter apparent digestibility;OM:有机物表观消化率 organic matter apparent digestibility;CP:粗蛋白质表观消化率 crude protein apparent digestibility;GE:总能表观消化率 gross energy apparent digestibility;NDF:中性洗涤纤维表观消化率 neutral detergent fiber apparent digestibility;ADF:酸性洗涤纤维表观消化率 acid detergent fiber apparent digestibility;Genus:菌属相对丰度 bacterial genera relative abundance;Shannon:Shannon指数 Shannon index;Simpson:Simpson指数 Simpson index;chao1:Chao1指数 Chao1 index;ACE:Ace指数 Ace index。 Fig.4 Network combination diagram of correlation analysis between nutrient apparent digestibility and bacterial genera and Alpha index |
图5 营养物质表观消化率与相对丰度前10菌属相关性分析热图GE:总能表观消化率 gross energy apparent digestibility;ADF:酸性洗涤纤维表观消化率 acid detergent fiber apparent digestibility;DM:干物质表观消化率 dry matter apparent digestibility;CP:粗蛋白质表观消化率 crude protein apparent digestibility;OM:有机物表观消化率 organic matter apparent digestibility;NDF:中性洗涤纤维表观消化率 neutral detergent fiber apparent digestibility;Quinella:奎氏菌属;uncultured_rumen_bacterium:未培养瘤胃细菌;unclassified_Lachnospiraceae:未分类毛螺菌科;Christensenellaceae_R_7_group:克里斯滕森菌科R_7群;NK4A214_group:NK4A214群;Prevotella:普雷沃氏菌属;Rikenellaceae_RC9_gut_group:理研菌科RC9肠道群;Candidatus_Saccharimonas:候选单胞生糖菌属;Lachnospiraceae_XPB1014_group:毛螺菌科XPB1014群。 *:显著相关(P<0.05)。*: significant correlation (P<0.05). Fig.5 Heat map of correlation analysis between nutrient apparent digestibility and relative abundance of top 10 bacterial genera |
