1 材料与方法
1.1 试验设计
表1 基础饲粮组成及营养水平(干物质基础)Table 1 Composition and nutrient levels of the basal diet (DM basis) % |
| 项目Items | 含量Content | ||
|---|---|---|---|
| 原料Ingredients | |||
| 玉米Corn | 12.00 | ||
| 小麦麸Wheat bran | 13.90 | ||
| 豆粕Soybean meal | 13.00 | ||
| 玉米胚芽粕Corn germ meal | 10.60 | ||
| 次粉Wheat middling | 6.00 | ||
| 艾叶粉Mugwort powder | 5.00 | ||
| 花生秧粉Peanut powder | 20.00 | ||
| 花生壳粉Peanut shell powder | 10.00 | ||
| 谷糠Chaff | 5.00 | ||
| 氯化钠NaCl | 0.30 | ||
| DL-蛋氨酸DL-Met | 0.20 | ||
| L-赖氨酸L-Lys | 0.30 | ||
| 预混料Premix1) | 4.00 | ||
| 合计Total | 100.00 | ||
| 营养水平Nutrient levels2) | |||
| 消化能DE/(MJ/kg) | 10.61 | ||
| 粗蛋白质CP | 15.39 | ||
| 粗纤维CF | 18.50 | ||
| 粗脂肪EE | 2.25 | ||
| 钙Ca | 0.98 | ||
| 有效磷AP | 0.56 | ||
| 赖氨酸Lys | 0.81 | ||
| 蛋氨酸+半胱氨酸Met+Cys | 0.61 | ||
| 中性洗涤纤维NDF | 36.44 | ||
| 酸性洗涤纤维ADF | 22.09 | ||
| 酸性洗涤木质素ADL | 5.60 |
1)预混料为每千克饲粮提供 The premix provided the following per kg of the diet:Fe 70 mg,Cu 20 mg,Zn 70 mg,Mn 10 mg,Co 0.15 mg,I 0.2 mg,Se 0.25 mg,VA 10 000 IU,VD 900 IU,VE 50 mg,VK2 2 mg,硫胺素 thiamine 2 mg,核黄素 riboflavin 6 mg,泛酸 pantothenic acid 50 mg,吡哆醇 pyridoxine 2 mg,VB12 0.02 mg,烟酸 nicotinic acid 50 mg,叶酸 folic acid 44 mg,胆碱 choline 1 000 mg,生物素 biotin 0.2 mg。 2)计算值 Calculated values。 |
1.2 饲养管理
1.3 样品采集与指标测定
1.3.1 样品采集
1.3.2 盲肠菌群高通量测序
1.3.3 粪便SCFAs含量测定
1.4 数据分析
2 结果
2.1 饲粮中添加丁酸梭菌对肉兔盲肠菌群结构的影响
2.1.1 Alpha多样性
表2 饲粮中添加丁酸梭菌对肉兔盲肠菌群Alpha多样性的影响Table 2 Effects of dietary C. butyricum on cecal flora Alpha diversity of meat rabbits |
| 项目 Items | 组别Groups | P值 P-value | |||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| CG | LC | MC | HC | ||||||||
| 32日龄32 days of age | |||||||||||
| Chao1指数Chao1 index | 747.18±14.05c | 979.72±26.90a | 881.43±22.23b | 858.92±75.60b | <0.001 | ||||||
| Ace指数Ace index | 728.86±21.60c | 985.76±10.06a | 844.77±3.12b | 829.22±17.35b | <0.001 | ||||||
| Shannon指数Shannon index | 4.37±0.28b | 5.32±0.03a | 4.56±0.11b | 2.75±0.68c | <0.001 | ||||||
| 39日龄39 days of age | |||||||||||
| Chao1指数Chao1 index | 1 020.38±23.66 | 1 039.74±7.07 | 983.04±15.94 | 986.91±39.73 | 0.052 | ||||||
| Ace指数Ace index | 973.37±14.95b | 1 032.04±9.31a | 987.12±13.37ab | 936.84±48.69b | <0.001 | ||||||
| Shannon指数Shannon index | 5.11±0.10ab | 6.00±0.16a | 4.92±0.24c | 5.22±0.21b | <0.001 | ||||||
| 46日龄46 days of age | |||||||||||
| Chao1指数Chao1 index | 1 155.45±29.38 | 1 108.54±76.29 | 1 132.39±15.05 | 1 073.95±28.47 | 0.202 | ||||||
| Ace指数Ace index | 1 132.09±13.08b | 1 078.33±12.44c | 1 194.38±26.51a | 1 074.92±48.40c | <0.001 | ||||||
| Shannon指数Shannon index | 6.21±0.23a | 6.11±0.12a | 6.14±0.04a | 5.33±0.23b | <0.001 | ||||||
| 53日龄53 days of age | |||||||||||
| Chao1指数Chao1 index | 1 063.97±43.26bc | 1 104.43±79.83b | 1 270.89±17.67a | 1 005.48±18.67c | <0.001 | ||||||
| Ace指数Ace index | 1 100.74±53.61b | 1 023.50±11.41c | 1 259.40±7.69a | 972.00±21.97d | <0.001 | ||||||
| Shannon指数Shannon index | 5.80±0.04a | 5.06±0.16b | 5.85±0.22a | 5.58±0.19a | <0.001 | ||||||
| 60日龄60 days of age | |||||||||||
| Chao1指数Chao1 index | 1 053.57±35.97b | 1 138.02±27.50a | 1 071.79±8.70b | 1 090.20±13.08b | 0.007 | ||||||
| Ace指数Ace index | 1 111.45±76.05ab | 1 131.00±13.49a | 1 046.96±43.97b | 1 113.27±7.71ab | <0.001 | ||||||
| Shannon指数Shannon index | 5.67±0.27b | 5.71±0.09b | 5.85±0.21b | 6.34±0.16a | 0.003 | ||||||
| 67日龄67 days of age | |||||||||||
| Chao1指数Chao1 index | 1 062.78±31.48b | 1 163.47±17.26a | 1 134.32±24.85a | 1 079.31±21.77b | 0.002 | ||||||
| Ace指数Ace index | 1 101.47±22.51c | 1 257.00±25.45a | 1 141.40±8.68b | 1 091.22±19.49c | <0.001 | ||||||
| Shannon指数Shannon index | 5.46±0.12ab | 5.16±0.15c | 5.63±0.09a | 5.30±0.29bc | 0.019 | ||||||
同行数据肩标不同小写字母表示差异显著(P<0.05),相同或无字母表示差异不显著(P>0.05)。下表同。 | |
In the same row, values with different small letter superscripts mean significant difference (P<0.05), while with the same or no letter superscripts mean no significant difference (P>0.05). The same as below. |
2.1.2 Beta多样性
2.1.3 门水平上相对丰度
表3 饲粮中添加丁酸梭菌对肉兔盲肠菌群门水平上相对丰度的影响Table 3 Effects of dietary C. butyricum on cecal flora relative abundance of meat rabbits at phylum level % |
| 项目 Items | 组别Groups | P值 P-value | |||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| CG | LC | MC | HC | ||||||||
| 32日龄32 days of age | |||||||||||
| 厚壁菌门Firmicutes | 34.76±8.54 | 43.47±6.16 | 44.46±5.46 | 44.37±4.14 | 0.095 | ||||||
| 拟杆菌门Bacteroidetes | 53.67±8.00 | 45.00±6.87 | 46.26±5.64 | 44.36±4.77 | 0.152 | ||||||
| 疣微菌门Verrucomicrobia | 6.70±1.43 | 6.80±0.13 | 6.37±0.34 | 6.47±0.35 | 0.815 | ||||||
| 放线菌门Actinobacteria | 2.86±0.48a | 2.51±0.54a | 1.77±0.49b | 2.61±0.23a | 0.014 | ||||||
| 其他Other | 1.59±0.69 | 1.74±0.95 | 0.72±0.43 | 1.67±0.94 | 0.176 | ||||||
| 未归类Unclassified | 0.39±0.24 | 0.44±0.36 | 0.39±0.37 | 0.49±0.41 | 0.972 | ||||||
| 39日龄39 days of age | |||||||||||
| 厚壁菌门Firmicutes | 42.20±11.84 | 47.47±10.32 | 51.88±11.06 | 55.21±17.64 | 0.483 | ||||||
| 拟杆菌门Bacteroidetes | 47.54±10.69 | 42.59±10.27 | 38.85±10.80 | 34.01±18.16 | 0.467 | ||||||
| 疣微菌门Verrucomicrobia | 5.36±0.93 | 6.03±0.76 | 5.95±0.27 | 5.76±0.22 | 0.395 | ||||||
| 放线菌门Actinobacteria | 2.70±0.37ab | 2.42±0.33bc | 2.12±0.45c | 2.98±0.16a | 0.012 | ||||||
| 其他Other | 1.45±0.14 | 1.12±0.46 | 0.95±0.49 | 1.58±0.43 | 0.128 | ||||||
| 未归类Unclassified | 0.71±0.28 | 0.34±0.28 | 0.22±0.12 | 0.43±0.34 | 0.059 | ||||||
| 46日龄46 days of age | |||||||||||
| 厚壁菌门Firmicutes | 55.33±8.92 | 53.31±11.34 | 65.83±8.28 | 65.18±12.74 | 0.140 | ||||||
| 拟杆菌门Bacteroidetes | 34.57±8.76 | 35.82±11.31 | 23.12±8.22 | 24.38±12.74 | 0.129 | ||||||
| 疣微菌门Verrucomicrobia | 5.52±0.35 | 6.07±0.59 | 6.05±0.42 | 6.07±0.23 | 0.144 | ||||||
| 放线菌门Actinobacteria | 2.53±0.26 | 2.76±0.35 | 2.91±0.31 | 2.51±0.36 | 0.181 | ||||||
| 其他Other | 1.51±0.41 | 1.49±0.41 | 1.53±0.36 | 1.32±0.31 | 0.781 | ||||||
| 未归类Unclassified | 0.51±0.27 | 0.51±0.10 | 0.53±0.34 | 0.52±0.43 | 1.000 | ||||||
| 53日龄53 days of age | |||||||||||
| 厚壁菌门Firmicutes | 54.10±16.83 | 56.43±15.21 | 58.47±14.00 | 61.43±7.56 | 0.834 | ||||||
| 拟杆菌门Bacteroidetes | 36.34±16.39 | 32.88±15.07 | 31.11±14.57 | 28.50±7.34 | 0.812 | ||||||
| 疣微菌门Verrucomicrobia | 5.66±0.77 | 5.81±0.47 | 5.98±0.35 | 5.85±0.66 | 0.890 | ||||||
| 放线菌门Actinobacteria | 2.43±0.35 | 2.63±0.35 | 2.70±0.21 | 2.46±0.19 | 0.416 | ||||||
| 其他Other | 0.92±0.56 | 1.75±1.10 | 1.14±0.61 | 1.26±0.50 | 0.351 | ||||||
| 未归类Unclassified | 0.52±0.42 | 0.47±0.27 | 0.58±0.40 | 0.48±0.23 | 0.966 | ||||||
| 60日龄60 days of age | |||||||||||
| 厚壁菌门Firmicutes | 61.44±6.18 | 66.77±9.33 | 65.79±6.51 | 65.63±10.05 | 0.794 | ||||||
| 拟杆菌门Bacteroidetes | 28.65±5.92 | 24.69±9.74 | 25.30±6.39 | 26.34±9.25 | 0.895 | ||||||
| 疣微菌门Verrucomicrobia | 5.18±0.21a | 4.14±0.59b | 4.14±0.46b | 3.56±0.51b | 0.001 | ||||||
| 放线菌门Actinobacteria | 2.82±0.14 | 2.42±0.36 | 2.72±0.34 | 2.46±0.32 | 0.183 | ||||||
| 其他Other | 1.37±0.40 | 1.46±0.34 | 1.54±0.51 | 1.39±0.65 | 0.952 | ||||||
| 未归类Unclassified | 0.51±0.22 | 0.49±0.27 | 0.48±0.30 | 1.75±2.96 | 0.443 | ||||||
| 67日龄67 days of age | |||||||||||
| 厚壁菌门Firmicutes | 68.33±8.62 | 67.85±12.13 | 72.53±14.61 | 72.71±11.36 | 0.854 | ||||||
| 拟杆菌门Bacteroidetes | 21.71±8.71 | 23.54±13.18 | 19.16±14.77 | 19.56±11.50 | 0.925 | ||||||
| 疣微菌门Verrucomicrobia | 5.34±0.17a | 4.11±0.38bc | 4.19±0.50b | 3.56±0.54c | <0.001 | ||||||
| 放线菌门Actinobacteria | 3.06±0.29 | 2.55±0.39 | 2.48±0.30 | 2.66±0.31 | 0.073 | ||||||
| 其他Other | 0.88±0.36 | 1.30±0.78 | 0.95±0.34 | 1.00±0.33 | 0.522 | ||||||
| 未归类Unclassified | 0.65±0.33 | 0.62±0.37 | 0.67±0.26 | 0.48±0.28 | 0.755 | ||||||
2.1.4 LEfSe分析结果
图2 32日龄组间差异LEfSe图Clostridium_sensu_stricto_1:严格梭菌属;Clostridiaceae_1:梭菌科1;Tannerellaceae:坦纳菌科;Parabacteroides:类芽孢杆菌属;Bacteroides:拟杆菌属;Acidobacteriia:酸杆菌门;Parabacteroides_distasonis:青春双歧杆菌;Bacteroides_fragilis:脆弱拟杆菌;Streptomycetaceae:链霉菌科;Streptomycetales:链霉菌目;Bacteroides_eggerthi:埃氏拟杆菌;Maihella_massiliensis:马赛耶尔森氏菌属;Sphingobacterium:鞘氨醇杆菌属;Sphingobacteriales:鞘脂杆菌目;Sphingobacteriaceae:鞘脂杆菌科;Sphingobacterium:鞘氨醇杆菌。 Fig.2 Difference LEfSe diagram at 32 days of age among groups |
图3 39日龄组间差异LEfSe图Bacteroides_nordi:诺迪拟杆菌属;Rhizobiales:根瘤菌属;Thermileophila:嗜热油菌纲;Pantoea_dispersa:分散泛菌;Pantoea:泛菌属;Acidobacteriales:酸杆菌目; Rhodobacteraceae:红杆菌科;Xanthomonadales:黄色单胞菌目;Phiascolarctobacterium_faecium:粪考拉杆菌;Enterobacter:肠杆菌属;Sphingomohadales:鞘氨醇单胞菌目;Sphingomonadaceae:鞘氨醇单胞菌科;Bacillaceae:芽孢杆菌科;Sphingomonas:鞘脂单胞菌属;Actinobacteria:放线菌门;Anaerovorax:优杆菌属;Oscilibacter:颤杆菌克属;Burkholderiaceae:伯克氏菌科;Betaprotpobacteriales:变形菌目;Clostridium_methylpentosum_DSM_5476:甲基戊糖梭菌DSM_5476;Ruminiclostridium:瘤胃梭菌属;Ruminiclostridium_6:瘤胃梭菌属6;Ruminococcaceae_V9D2013_group:瘤胃球菌科V9D2013_组;Akkermainsia:阿克曼菌属;Akkermansiaceae:阿克曼斯亚科;Verrucormicrobiales:疣微菌目;Verrucomicrobiae:疣微菌纲;Verrucormicrobia:疣微菌门;Burkholderiaceae:伯克氏菌科;Bacteroides:拟杆菌属;Lachnoclostridium:乳酸杆菌属;Butyricicoccus:丁酸球菌属;Ruminiclostridium:瘤胃梭菌属;Anaerovorax:优杆菌属;Xanthomonadales:黄色单胞菌目;Rhodobacterales:红杆菌目;Rhodanobacteraceae:罗河杆菌科;Rhodobacteraceae:红杆菌科;Bacillus:芽孢杆菌属;Shimia:沈氏菌属;Acidobacteriales:酸杆菌目;Rhodanobacter:罗河杆菌属;Acidobacteriia:酸杆菌科;Bacillaceae:芽孢杆菌科;Rhodanobacter_sp_C06:罗河杆菌属sp_C06。 Fig.3 Difference LEfSe diagram at 39 days of age among groups |
图4 46日龄组间差异LEfSe图Escherichia_Shigella:大肠埃希菌-志贺氏菌属;Enterobacteriaceae:肠杆菌科;Enterobacteriales:肠杆菌目;Coriobacteriales_bacterium_DNF00809:红蝽菌目杆菌bacterium_DNF00809;Caulobacteraceae:柄杆菌科;Caulobacterales:柄杆菌目;Brevundimonas:短波单细胞菌属;Roseburia:罗氏菌属。 Fig.4 Difference LEfSe diagram at 46 days of age among groups |
图5 53日龄组间差异LEfSe图Lachnoclostridium:乳酸杆菌属。 Fig.5 Difference LEfSe diagram at 53 days of age among groups |
图6 60日龄组间差异LEfSe图Ruminiclostridium:瘤胃梭菌属。 Fig.6 Difference LEfSe diagram at 60 days of age among groups |
图7 67日龄组间差异LEfSe图Thauera:陶厄氏菌属;Sphingomohadales:鞘氨醇单胞菌目;Sphingomonadaceae:鞘氨醇单胞菌科;Streptomycetaceae:链霉菌科;Streptomycetales:链霉菌目;Sphingomonas:鞘脂单胞菌属;Actinobacteria:放线菌门;Bacillaceae:芽孢杆菌科;Enterobacter:肠杆菌属;Pseudomoradales:假单胞菌目;Eisenbergiella:塔伊艾森伯格菌属;Enterobacteriales:肠杆菌目;Escherichia_Shigella:大肠埃希菌-志贺氏菌属;Enterobacteriaceae:肠杆菌科;Rhodocyclaceae:红环菌科;Thauera:陶厄氏菌属。 Fig.7 Difference LEfSe diagram at 67 days of age among groups |
2.1.5 肠道菌群动态变化
图8 不同分组内相对丰度前15菌属动态变化气泡图NA:未归类 unclassified;Ruminococcaceae_NK4A214_group:瘤胃球菌科NK4A214群;Ruminococcaceae_UCG-014:瘤胃球菌科UCG-014;Bacteroides:拟杆菌属;Lachnospiraceae_NK4A136_group:毛螺菌科NK4A136群;Akkermansia:阿克曼氏菌属;dgA-11_gut_group:dgA-11肠道群;Ruminococcaceae_UCG-005:瘤胃球菌科UCG-005;Alistipes:另枝菌属;Ruminococcaceae_V9D2013_group:瘤胃球菌科V9D2013群;Ruminococcus_1:瘤胃球菌属1;Ruminiclostridium_6:罗马尼亚梭菌6;Ruminococcaceae_UCG-013:瘤胃球菌科UCG-013;Other:其他。 Fig.8 Bubble plot of dynamic changes of relative abundance top10 bacteria in different groups |
2.2 饲粮中添加丁酸梭菌对肉兔粪便SCFAs含量的影响
表4 各试验组与CG组53日龄粪便SCFAs含量差异Table 4 Differences of fecal SCFAs contents between each test group and CG group on 53 days of age |
| 项目 Items | LC组较CG组 LC group vs CG group | MC组较CG组 MC group vs CG group | HC组较CG组 HC group vs CG group | ||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| P值 P-value | 差异倍数 Fold change | 趋势 Trend | P值 P-value | 差异倍数 Fold change | 趋势 Trend | P值 P-value | 差异倍数 Fold change | 趋势 Trend | |||
| 乙酸AA | 0.007 | 0.542 | ↑ | 0.800 | 1.036 | ↑ | 0.140 | 1.326 | ↑ | ||
| 丙酸PA | 0.065 | 0.701 | ↑ | 0.281 | 1.194 | ↑ | 0.171 | 1.247 | ↑ | ||
| 丁酸BA | 0.014 | 0.371 | ↓ | 0.520 | 0.790 | ↑ | 0.505 | 0.781 | ↑ | ||
| 异丁酸IBA | 0.073 | 0.458 | ↓ | 0.555 | 0.833 | ↑ | 0.725 | 0.884 | ↑ | ||
| 戊酸VA | 0.001 | 0.535 | ↑ | 0.657 | 0.878 | ↑ | 0.510 | 0.921 | ↑ | ||
| 异戊酸IVA | 0.088 | 0.515 | ↑ | 0.602 | 0.840 | ↑ | 0.868 | 1.065 | ↑ | ||
| 己酸HA | 0.560 | 0.703 | ↑ | 0.868 | 0.866 | ↑ | 0.412 | 0.575 | ↑ | ||
表5 各试验组与CG组67日龄粪便SCFAs含量差异Table 5 Differences of fecal SCFAs contents between each test group and CG group on 67 days of age |
| 项目 Items | LC组较CG组 LC group vs CG group | MC组较CG组 MC group vs CG group | HC组较CG组 HC group vs CG group | ||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| P值 P-value | 差异倍数 Fold change | 趋势 Trend | P值 P-value | 差异倍数 Fold change | 趋势 Trend | P值 P-value | 差异倍数 Fold change | 趋势 Trend | |||||||||||||
| 乙酸AA | 0.242 | 0.674 | ↑ | 0.767 | 1.078 | ↑ | 0.540 | 0.874 | ↑ | ||||||||||||
| 丙酸PA | 0.363 | 0.760 | ↑ | 0.998 | 1.001 | ↑ | 0.104 | 0.767 | ↑ | ||||||||||||
| 丁酸BA | 0.478 | 0.691 | ↑ | 0.554 | 1.319 | ↑ | 0.685 | 0.815 | ↑ | ||||||||||||
| 异丁酸IBA | 0.086 | 0.577 | ↑ | 0.375 | 0.770 | ↑ | 0.307 | 0.747 | ↑ | ||||||||||||
| 戊酸VA | 0.087 | 0.438 | ↓ | 0.947 | 1.028 | ↑ | 0.353 | 0.698 | ↑ | ||||||||||||
| 异戊酸IVA | 0.034 | 0.465 | ↓ | 0.411 | 0.751 | ↑ | 0.159 | 0.647 | ↑ | ||||||||||||
| 己酸HA | 0.237 | 0.318 | ↓ | 0.367 | 0.460 | ↓ | 0.664 | 0.741 | ↑ | ||||||||||||
2.3 SCFAs与盲肠菌群关联性分析
图9 肉兔盲肠相对丰度排名前10菌属与主要粪便SCFAs含量关联分析热图A:53日龄关联分析热图 heat map of correlation analysis at 53 days of age;B:67日龄关联分析热图 heat map of correlation analysis at 67 days of age。 图中颜色深浅表示菌属相对丰度与SCFAs含量相关性大小,其中红色表示正相关,蓝色表示负相关。*表示相关性显著(P<0.05),**表示相关性极显著(P<0.01)。 Fig.9 Heat map of correlation analysis between relative abundance of top 10 bacteria in cecum and content of SCFAs in feces of meat rabbits The color depth in the figure indicated the correlation between bacteria relative abundance and SCFAs content, red indicated positive correlation, and blue indicated negative correlation. * indicated significant correlation (P<0.05), ** indicated extremely significant correlation (P<0.01). |

京公网安备 11010802041926号