In total, 3,110,414 raw reads were detected in 32 sequenced chyme samples. After splicing, filtering, and removing chimeras, non-chimeric reads were clustered into ASVs using Silva (v138) with standard > 70% (default). Sequencing depth, like rarefaction curves (Fig. S1), reflected the total microbial species richness. In Figs. 3a-e, the HP diet boosted higher richness in α-diversity, as revealed by higher observed-ASVs (P = 0.013) and Chao1 (P = 0.024) than the MP diet. While in Figs. 3c-3e, the interaction indicated that the addition of ENZs to the HP diet induced higher Shannon index and Pielou's evenness whereas decreased Shannon index in the MP diet (P < 0.05). Moreover, the interactive effect of MP and ENZs tended to reduce the Simpson index (P = 0.077).
A PERMANOVA analysis (Table S2) including Bray-Curtis and unweighted_ UniFrac indicated that the microbial community was separated among four groups, as visualized in Figs. 4a-d.
Figures 5a-j display the top 10 phyla. The MP diet promoted the growth of Firmicutes in the colon (P = 0.014) while the HP diet enhanced the abundance of Spirochaetota, Verrucomicrobiota, Desulfobacterota, and Fibrobacterota (P < 0.05). By supplementing ENZs in the diet, colonic Actinobacteriota and Desulfobacterota abundance were elevated while Verrucomicrobiota abundance was lowered (P < 0.05) compared with the no-ENZs groups. The interaction between HP and ENZs reduced the abundance of Firmicutes but increased that of Bacteroidetes and Desulfobacterota abundance compared with the interactive effect of MP and ENZs (P < 0.05).
Figures 6a-j show the top 10 bacteria at the family level. Compared to the MP diet, the HP diet increased the abundance of Lactobacillaceae, Oscillospiraceae, Muribaculaceae, and Spirochaetaceae (P < 0.01). Meanwhile, the MP diet prompted higher abundance of Clostridiaceae and Selenomonadaceae (P < 0.01). The ENZs intervention impaired the growth of Clostridiaceae and Peptostreptococcaceae but enriched Lactobacillaceae abundance (P < 0.01). The interactions revealed that when ENZs was added to HP and MP diets, Lactobacillaceae and Muribaculaceae abundance were increased, while Peptostreptococcaceae and Selenomonadaceae abundance were just lowered in the HP diet (P < 0.001).
The top 30 genera were listed in Table 8. Dietary HP boosted the abundance of Lactobacillus, Muribaculaceae, Romboutsia, Treponema, Rikenellaceae_RC9_gut_group, Oscillospiraceae-NK4A214_group, Streptococcus, Turicibacter, [Eubacterium]_coprostanoligenes_group, Ruminococcus, and UCG-010 (P < 0.05), compared with the MP diet. On the other hand, the MP diet increased the abundance of Clostridium_sensu_stricto_1, Terrisporobacter, Prevotellaceae_NK3B31_group, Ruminococcaceae-undifitened, Selenomonadaceae-unculture, and Shuttleworthia. Dietary ENZs enriched the abundance of Lactobacillus, Rikenellaceae_RC9_gut_group, Oscillospiraceae-UCG-002, Megasphaera, Coprococcus, and Succinivibrio abundance whereas lowered that of Clostridium_sensu_stricto_1, Terrisporobacter, Oscillospiraceae-UCG-005, Prevotellaceae_NK3B31_group, Romboutsia, and Turicibacter, compared with groups without ENZs. There were interactions between CP levels and ENZs at the genus level. The ENZs supplement to the HP diet lowered Clostridium_sensu_stricto_1, Terrisporobacter, Romboutsia, and Turicibacter abundance but enriched Lactobacillus, Muribaculaceae, Ruminococcaceae-undifitened, Rikenellaceae_RC9_gut_group, and Blautia abundance. Differently, ENZs addition to the MP diet increased Lactobacillus but lowered Muribaculaceae abundance.
Table 8
Effect of dietary ENZs on genera in colonic digesta of finishing pigs (%)1
CP levels | HP | MP | | | CP × ENZs | Pooled SEM | P-value |
ENZs | | | + | - | HP- | HP+ | MP- | MP+ | CP | ENZs | Interaction |
Clostridium-sensu-stricto-1 | 27.51b | 42.27a | 26.27b | 43.51a | 45.83a | 9.19b | 41.19a | 43.35a | 2.24 | < 0.001 | < 0.001 | < 0.001 |
Lactobacillus | 17.70a | 6.46b | 21.72a | 2.43b | 2.13c | 33.26a | 2.73c | 10.19b | 1.35 | < 0.001 | < 0.001 | < 0.001 |
Terrisporobacter | 4.30b | 5.93a | 3.52b | 6.72a | 7.36a | 1.25b | 6.08a | 5.79a | 0.57 | 0.008 | < 0.001 | < 0.001 |
Oscillospiraceae-UCG-005 | 3.83 | 2.93 | 2.52b | 4.24a | 4.74 | 2.92 | 3.73 | 2.12 | 0.45 | 0.052 | 0.001 | 0.812 |
Prevotella | 2.86 | 3.59 | 4.23 | 2.19 | 1.64 | 4.07 | 2.74 | 4.44 | 1.14 | 0.527 | 0.081 | 0.753 |
Muribaculaceae | 2.75a | 1.32b | 2.24 | 1.83 | 2.24b | 3.27a | 1.43bc | 1.21c | 0.25 | < 0.001 | 0.117 | 0.020 |
Prevotellaceae-NK3B31-group | 1.09b | 2.62a | 1.03b | 2.68a | 0.85b | 1.33b | 4.51a | 0.72b | 0.65 | 0.026 | 0.017 | 0.003 |
Ruminococcaceae- undifitened | 0.51a | 0.37b | 0.40b | 0.47a | 0.37b | 0.65a | 0.44b | 0.30b | 0.05 | 0.008 | 0.170 | < 0.001 |
Romboutsia | 2.05a | 1.29b | 0.83b | 2.51a | 3.83a | 0.26c | 1.18b | 1.40b | 0.09 | < 0.001 | < 0.001 | < 0.001 |
Treponema | 2.75a | 0.44b | 1.27 | 1.92 | 2.31 | 3.20 | 0.24 | 0.64 | 0.63 | 0.001 | 0.317 | 0.705 |
Unidentified | 0.21 | 0.10 | 0.08b | 0.24a | 0.05c | 0.38a | 0.11b | 0.10b | 0.07 | 0.100 | 0.031 | 0.019 |
Roseburia | 1.12 | 1.59 | 1.27 | 1.45 | 0.40b | 1.85a | 2.51a | 0.69b | 0.24 | 0.064 | 0.464 | < 0.001 |
Christensenellaceae-R-7-group | 1.24 | 0.82 | 1.27 | 0.79 | 1.27 | 1.21 | 0.31 | 1.34 | 0.37 | 0.276 | 0.212 | 0.156 |
Rikenellaceae-RC9-gut-group | 1.52a | 0.49b | 1.31a | 0.73b | 0.89b | 2.16a | 0.57b | 0.47b | 0.20 | < 0.001 | 0.006 | 0.002 |
Oscillospiraceae-UCG-002 | 1.24 | 0.79 | 1.32a | 0.71b | 0.87 | 1.60 | 0.55 | 1.03 | 0.22 | 0.053 | 0.011 | 0.582 |
Megasphaera | 0.99 | 1.01 | 1.66a | 0.34b | 0.87 | 1.60 | 0.55 | 1.03 | 0.22 | 0.948 | 0.001 | 0.077 |
Coprococcus | 0.73 | 0.96 | 1.06a | 0.64b | 0.405 | 1.06 | 0.87 | 1.05 | 0.131 | 0.097 | 0.004 | 0.077 |
Oscillospiraceae-NK4A214-group | 1.06a | 0.58b | 0.94 | 0.70 | 0.97 | 1.15 | 0.44 | 0.72 | 0.14 | 0.001 | 0.194 | 0.904 |
Blautia | 0.79 | 0.84 | 0.94 | 0.69 | 0.38c | 1.20a | 1.00ab | 0.67bc | 0.17 | 0.796 | 0.149 | 0.002 |
Turicibacter | 1.00a | 0.57b | 0.44b | 1.13a | 1.86a | 0.14c | 0.40bc | 0.75b | 0.13 | 0.002 | < 0.001 | < 0.001 |
Phascolarctobacterium | 0.78 | 0.74 | 0.87 | 0.65 | 0.60 | 0.97 | 0.70 | 0.78 | 0.12 | 0.709 | 0.072 | 0.231 |
Streptococcus | 1.21a | 0.28b | 0.65 | 0.85 | 1.66a | 0.76ab | 0.04b | 0.53ab | 0.31 | 0.005 | 0.524 | 0.032 |
[Eubacterium]-coprostanoligenes-group | 0.99a | 0.45b | 0.80 | 0.64 | 0.93 | 1.04 | 0.34 | 0.55 | 0.16 | 0.002 | 0.336 | 0.746 |
Succinivibrio | 0.77 | 0.62 | 1.02a | 0.37b | 0.22 | 1.32 | 0.51 | 0.72 | 0.23 | 0.493 | 0.008 | 0.059 |
Ruminococcus | 0.85a | 0.44b | 0.73 | 0.56 | 0.70 | 1.00 | 0.42 | 0.45 | 0.08 | < 0.001 | 0.052 | 0.110 |
Anaerovibrio | 0.41 | 0.83 | 0.66 | 0.59 | 0.35 | 0.47 | 0.82 | 0.85 | 0.23 | 0.077 | 0.757 | 0.844 |
UCG-010 | 0.81a | 0.43b | 0.67 | 0.57 | 0.67 | 0.95 | 0.46 | 0.39 | 0.13 | 0.007 | 0.429 | 0.205 |
Selenomonadaceae-uncultured | 0.53b | 2.88a | 2.02 | 1.38 | 0.13c | 0.93bc | 1.84b | 3.92a | 0.33 | < 0.001 | 0.063 | < 0.001 |
Shuttleworthia | 0.06b | 0.98a | 0.67 | 1.37 | 0 | 0.12 | 0.75 | 1.21 | 0.24 | 0.001 | 0.228 | 0.484 |
Acetitomaculum | 0.51a | 0.13b | 0.47a | 0.16b | 0.73a | 0.28b | 0.21b | 0.04c | 0.04 | < 0.001 | < 0.001 | 0.001 |
1Each value represents the means of 8 replicate pens a, b Represents a row that differs significantly (P < 0.05). Abbreviation: HP, 16.73% protein level, MP, 14.73% protein level; ENZs: the combination of 2,500 U/kg α-amylase,100 U/kg protease, and 10,000 U/kg β-glucanase. |
Short chain fatty acid concentrations in the colonic digesta.
As shown in Table 9, HP contributed to higher SCFA and acetate concentrations and tended to increase isobutyric acid concentration, compared to the MP diet (P < 0.001). The interactive effect of HP and ENZs resulted in higher SCFA and acetate concentrations (P ≤ 0.05) and had the tendency to increase butyrate concentration (P = 0.087).
Table 9
The effects of dietary CP levels and multi-enzyme complex on the VFA of finishing pigs (umoL/ g dry feces)1
CP levels | HP | MP | | | CP × ENZs | Pooled SEM | P-value |
ENZs | | | + | - | HP- | HP+ | MP- | MP+ | CP | ENZs | Interaction |
SCFA | 91.08a | 83.52b | 87.75 | 86.84 | 89.44ab | 92.71a | 84.25bc | 82.79c | 1.15 | < 0.001 | 0.439 | 0.050 |
Acetate | 36.14a | 29.67b | 32.75 | 33.06 | 35.07a | 37.21a | 30.44b | 28.91b | 0.76 | < 0.001 | 0.694 | 0.023 |
Propionate | 22.50 | 21.94 | 21.95 | 22.49 | 22.16 | 22.85 | 21.75 | 22.13 | 0.56 | 0.325 | 0.347 | 0.778 |
Butyrate | 11.16 | 10.81 | 10.84 | 11.13 | 10.82 | 11.51 | 10.86 | 10.76 | 0.22 | 0.129 | 0.204 | 0.087 |
Isovaleric acid | 15.40 | 15.01 | 15.49 | 4.96 | 15.78 | 15.01 | 15.20 | 14.90 | 0.37 | 0.358 | 0.160 | 0.539 |
Isovaleric acid | 5.91 | 6.23 | 5.94 | 6.20 | 5.68 | 6.13 | 6.19 | 6.27 | 0.17 | 0.070 | 0.133 | 0.279 |
BCFA | 21.30 | 21.28 | 21.15 | 21.43 | 21.46 | 21.39 | 21.39 | 21.16 | 0.44 | 0.952 | 0.541 | 0.922 |
1Each value represents the means of 8 replicate pens a b Represents a row that differs significantly (P < 0.05) Abbreviation: HP, 16.73% protein level, MP, 14.73% protein level; ENZs: the combination of 2,500 U/kg α-amylase, 100 U/kg protease, and 10,000 U/kg β-glucanase. |
The spearman’s correlation analysis for short-chain fatty acids and colonic microbiota at the genus level.
Figure 7 shows that acetate concentration was positively correlated with the abundance of Lactobacillus, Muribaculaceae, Rikenellaceae_RC9_gut_group, Oscillospiraceae-NK4A214_group, Ruminococcus, UCG-010, and [Eubacterium]_coprostanoligenes_group, but negatively correlated with Clostridium_sensu_stricto_1, Shuttleworthia, Terrisporobacter, and Anaerovibrio abundance. Christensenellaceae_R-7_group, Rikenellaceae_RC9_gut_group, Oscillospiraceae-NK4A214_group, Oscillospiraceae-UCG-002, and UCG-010. Butyrate production was positively related to Lactobacillus, Muribaculaceae, Rikenellaceae_RC9_gut_group, and Oscillospiraceae-NK4A214_group abundance, but Clostridium_sensu_stricto_1 played the contrary function to produce butyrate. Isovaleric acid was associated with the abundance of Muribaculaceae.