Effects of P. monteilii JK-1 on grass carp growth performance
The effects of P. monteilii JK-1 on the growth performance of grass carp are shown in Table 1. After an 8-week feeding trial, fish fed with P. monteilii JK-1 supplemented diet had significantly higher final body weight, weight gain, and specific growth rates than that in control (P < 0.05).
Table 1
Effects of dietary P. monteilii JK-1 supplemented diet on growth performance of juvenile grass carp
Group | Initial body weight (g) | Final body weight (g) | Weight gain rate (%) | Specific growth rate (%/day) |
Control | 3 ± 0.05 | 10.76 ± 0.38a | 259.54 ± 14.02a | 2.28 ± 0.07a |
PM-JK1 | 3 ± 0.05 | 15.62 ± 0.74b | 416.88 ± 24.67b | 2.93 ± 0.08b |
Note: All data were expressed as mean ± SE. Different superscript letters in each column indicate significant differences between groups (P < 0.05). |
The protective effect of P. monteilii JK-1 on grass carp
We assessed the protected effects of P. monteilii JK-1 on the pathogen resistance of grass carp (Fig. 2). The results showed that grass carp in PM-JK1-T had considerably lower liver Aeromonas load than that in control-T (Fig. 2a). Notably, there was a significant decrease in the load of Aeromonas in the gut of grass carp after feeding P. monteilii JK-1. Also, grass carp in PM-JK1-T had a significant lower level of Aeromonas load in gut than that in the control-T after infection with Aeromonas (Fig. 2a). Meanwhile, the results of survival rate showed that grass carp in PM-JK1-T had higher survival rate than that in control-T (Fig. 2b), suggesting that P. monteilii JK-1 supplementation protect grass carp against the A. hydrophila infection.
Figure 2P. monteilii JK-1 protects grass carp against pathogen infection. a Aeromonas load (CFU/g) in liver and gut feeding with basal diet and P. monteilii JK-1 supplementation diet. The difference among the four group was constructed with Kruskal-Wallis test, significant differences (P < 0.05) among the four groups are denoted by different lowercase letters. b Kaplan-Meier graph of the grass carp survival after bath infection with A. hydrophila. One representative experiment of three independent biological repeats is shown. *P < 0.05, Gehan–Breslow–Wilcoxon test.
Effects of P. monteilii JK-1 on the immune response and antioxidant enzymes activity of grass carp
We evaluated the effects of fish-derived probiotics on the expression of immune-related genes in grass carp (Supplementary Fig. S1a, b). Compared to control, the expression of IL-1β, IL10, TNF-α, and TGF-β in head kidney was significantly improved in PM-JK1 (P < 0.05). Likewise, IL-1β, IL10, TNF-α, and TGF-β in gut were also recorded significantly higher in PM-JK1 than that in control.
Meanwhile, lower activities of SOD, CAT, and GPx in the liver were observed in control when compared to PM-JK1 (P < 0.05). Again, the PM-JK1 group displayed higher activities of SOD, CAT, and GPx in gut compared to control (Supplementary Fig. S1c, d).
Effects of P. monteilii JK-1 on gut microbiome
A total of 264,702 high-quality reads were obtained from 24 samples, after removing low quality reads, 61,785 ± 2916, 62,562 ± 2524, 71,744 ± 1272 and 66,123 ± 2259 were obtained for Control, group Control-T, group PM-JK1 and group PM-JK1-T, respectively. 4,041 OTUs with sequencing depths covering more than 99% were discovered in all samples, demonstrating significant sequencing coverage and the representativeness of the OTUs found in the samples. As shown in Fig. 5, compared to control, the Shannon index in PM-JK1 showed no significant difference, while in control-T showed a significant increase (Supplementary Fig. S2a). Meanwhile, there was no significant difference was observed between PM-JK1 and PM-JK1-T. The Richness also showed no significant change among the four groups (Supplementary Fig. S2b). PCoA based on Bray_curtis distance revealed there were significant differences among the groups in the distribution pattern of microbial communities (Supplementary Fig. S2c). There was no significant clear limit among the control, PM-JK1 and PM-JK1-T. However, the microbial communities in control-T were significantly different from the other three groups (control-T vs control: p = 0.006; control-T vs PM-JK1: p = 0.006; PM-JK1 vs PM-JK1-T: p = 0.006).
Analysis at the phylum level indicated that the most abundant phyla were Fusobacteria, Proteobacteria, and Actinobacteria, which accumulated abundance was higher than 90% (Fig. 3a). Fusobacteria was the most dominant phyla in control (81.42%), PM-JK1 (85.5%), and PM-JK1-T (69.37%), while Proteobacteria took dominance in control-T (48.7%) (Fig. 3a, b). At genus level, a similar bacterial composition was observed in control, PM-JK1, and PM-JK1-T, and the most abundance of genus was Cetobacterium (mean 81.2%, 85.4%, and 69.3%, respectively) (Fig. 3c). Of interest, the abundance of Bosea, Nocardia, Rhodobacteraceae_unclassified, Rhodococcus, and Phyllobacteriaceae_unclassified was significantly increased in control-T. The changes in bacterial composition among the four groups were then assessed using the linear discriminant analysis effect size (LEfSe) to pinpoint specific bacterial genera that were distinctive of the different treatments (Fig. 3d). The results revealed that Clostridiales_unclassified, Lachnospiraceae_unclassified, Propionibacterium, and Rhodobacter were enriched in control. A significant increase in Bosea, Nocardia, Rhodobacteraceae_unclassified, Rhodococcus, Phyllobacteriaceae_unclassified, Bradyrhizobiaceae_unclassified, Rhizobiaceae_unclassified, Alsobacter, Mycobacterium, and Microbacteriaceae_unclassified was observed in control-T. Furthermore, the abundance of Cetobacterium, Reyranells_unclassified, and Ruminococcaceae_unclassified was enriched in PM-JK1, while Romboutsia and Akkermansia were enriched in PM-JK1-T. Meanwhile, we used the STAMP to determine the significant differences in the abundance of microbiota between control and P. monteilii JK-1. As shown in Supplementary Fig. S3, dietary administration of P. monteilii JK-1 significantly improved the abundance of Akkermansia, Tessaracoccus, Ruminococcaceae_unclassified and Microbacterium, while decreasing the relative abundance of Firmicutes_unclassified and Thermus.
Figure 3 Bacterial compositions of the gut in different groups. a Relative abundance of the top 6 phyla in the four different groups. b Box plots showing the relative abundances of Proteobacteria and Fusobacteria in the four different groups. c Relative abundance of the top 20 genera in the gut samples from the four different groups. d Liner discriminant analysis effect size (LEfSe) analysis among the four groups, The LDA value threshold was set at 3.5.
Effects of P. monteilii JK-1 on gut ecological network
In this study, the phylogenetic molecular ecology networks (pMENs) were built using an RMT-based approach. There was significant difference between empirical networks and random networks (Table 2, Fig. 4). These four ecological networks may have exhibited typical small-word features because the average clustering coefficient and average path distance in empirical networks were much higher than those in random networks. In addition, the density, average clustering coefficient and average degree in PM-JK1 were significantly higher than that in control, while the average path distance was significantly lower, suggesting that dietary administration of P. monteilii JK-1 had significant effects on the gut microbiome ecological networks of grass carp. It's noteworthy that PM-JK1 has the greatest average degree, suggesting that the gut bacteria in PM-JK1 have the most complex network (Fig. 4).
Table 2
The topological properties of the gut microbiota co-occurrence networks in different groups and their respective identically sized random networks
Topological properties | Control | PM-JK1 | Control-T | PM-JK1-T |
Empirical Networks | | | | |
Nodes | 88 | 57 | 88 | 72 |
Total links | 189 | 281 | 256 | 259 |
Density | 0.049 | 0.176 | 0.067 | 0.101 |
Average clustering coefficient (ACC) | 0.303a | 0.418c | 0.379b | 0.440c |
Average degree (AD) | 4.295a | 9.860c | 5.818b | 7.194c |
Average path distance (GD) | 4.626c | 2.416a | 3.782b | 3.656b |
Modularity | 0.595b | 0.320a | 0.549b | 0.399a |
Random Networks | | | | |
Average path distance | 3.207 ± 0.058 | 2.163 ± 0.041 | 2.797 ± 0.042 | 2.472 ± 0.043 |
Average clustering coefficient | 0.052 ± 0.014 | 0.315 ± 0.026 | 0.091 ± 0.018 | 0.193 ± 0.022 |
Modularity | 0.426 ± 0.015 | 0.184 ± 0.010 | 0.334 ± 0.011 | 0.261 ± 0.011 |
Note: Random networks were generated by resetting all of the links of a matching empirical network with the same nodes and links. Data were generated from 100 random runs and SD indicates the standard deviation from the 100 runs. Different superscript letters in each row indicate significant differences between different groups (P < 0.05). |
Figure 4 The interspecific interactions among the gut bacterial community of grass carp are described by an ecological network. A link defined as a significant (P-value < 0.05) and strong (Spearman’s r > 0.6) relationship. Blue edges stand for negative correlations, while red edges stand for positive correlations. a The gut microbiota network in control. b The gut microbiota network in PM-JK1. c The gut microbiota network in control-T. d The gut microbiota network in PM-JK1-T.
The change in network complexity could cause changes in the role of individual nodes. To understand the ecological roles of the genera that were altered in the network, the within-module connectivity (Zi) and among-module connectivity (Pi) were calculated [20, 21] (Fig. 5). We noticed that the roles that the species play in the gut network in control were different from those in PM-JK1 (Fig. 5a). Meanwhile, the nodes in control were almost generalized to non-hubs, which means that the genera belonged to within-module connectivity (Fig. 5b). However, there were 20% of nodes in PM-JK1 were generalized to hubs, which means that the genera belonged to intra-module connectivity (Fig. 5b). Most importantly, the genus Cetobacterium was the only species that had the ability to connector hubs in PM-JK1, indicating that Cetobacterium played a key role in maintaining interactions between modules.
Figure 5 Node structural roles in networks within various groups were used to identify potential keystone taxa. a Nodes were partitioned into ultra-peripheral (Pi ≈ 0, Zi < 0), peripheral (Pi < 0.625, Zi < 0), non-hub connectors (0.625 < Pi < 0.8, Zi < 0), non-hub kinless nodes (Pi > 0.8, Zi < 0), provincial hubs (0 < Pi < 0.3, Zi > 0), connector hubs (0.3 < Pi < 0.75, Zi > 0), or kinless hubs (Pi > 0.75, Zi > 0). Red dots represent the nodes that in Control; Blue dots represent the nodes that in PM-JK1. b Relative abundance (%) of Ultra-peripheral node, Peripheral node, Connector hubs, Provincial hubs, Non-hub connector, and Non-hub kinless node within each group.
The closeness centrality, which is based on the distances from the node to all other nodes, is an importance of a vertex within a given complex network. A node with high closeness centrality plays an important role in ecological networks. We calculated the closeness centrality of each genus in different groups to further evaluate the effect of P. monteilii JK-1 on the gut ecological network (Supplementary Table S3). The results showed significant changes in the role of the bacteria in the ecological network of the gut before and after the infection. Microbacteriaceae_unclassified had the highest connectivity in control, while Cetobacterium had the highest connectivity in control-T, suggesting that Cetobacterium plays an important role within the ecological network of the gut after the infection (Supplementary Table S3). Moreover, dietary administration of P. monteilii JK-1 significantly increased the importance of Cetobacterium in the gut ecological network and suggested a role for Cetobacterium in the stability of gut ecological network.
Correlation Between The Gut Microbiota And Growth Performance And Disease Resistance
As shown in Supplementary Fig. S4, we used the Spearman’s correlation heatmap to demonstrate the relationship between the relative abundance of gut microbiota and growth performance and pathogen resistance. The results showed that Cetobacterium had significantly positive correlations with weight gain and specific growth rate, but Clostridiales_unclassified, Bradyrhizobiaceae_unclassified and Candidatus_Saccharibacteria_unclassified were opposite, suggesting that that Cetobacterium may contribute to the growth of grass carp. Meanwhile, Rhodobacteriaceae_unclassified and Akkemansia had positive relationships with the gut antioxidant response, and Akkemansia had a significant positive correlation with gut immune response, suggesting that Akkemansia may activate the intestinal immune response and antioxidant enzyme system.
To better understand the relationship among the P. monteilii JK-1 supplementation, keystone taxa, gut network, grass carp growth performance and infection level, a partial least squares path model (PLS-PM) was constructed (Fig. 6). The results proved that P. monteilii JK-1 supplementation had significant positive correlation with keystone taxa (0.73) and gut network (0.50), but had no significant correlation with growth performance and infection level, suggesting that P. monteilii JK-1 affect the growth performance and disease resistance of grass carp mainly through an indirect way (Fig. 6a, b). The keystone taxa (Cetobacterium and Akkemansia) positively correlated with gut work (0.52) and growth performance (0.40) but not significantly correlated with the infection level, indicating that Cetobacterium and Akkemansia had significant impacts on the gut network structure and growth performance of grass carp. Moreover, the gut network showed significant impacts on the growth performance (1.01) and the infection level (-0.80), suggesting that a stable gut network is conducive to improving growth performance and reducing infection levels in grass carp.
Figure 6 Effects of the major factors on the growth performance and pathogen resistance as determined by the partial least squares path model (PLS-PM) analysis. a PLS-PM showing the cascading relationships of different factors. The dashed rectangles represent the loading for the keystone taxa, the gut network, growth performance, and infection level, which generate the latent variables. After 1,000 bootstraps, path coefficients are calculated and represented by the width of the arrow (red stands for positive relationship, green stands for negative relationship). The dashed arrow indicates a coefficient that did not differ substantially from 0 (P > 0.05). The GoF statistic was used to evaluate the model, and the GoF value was 0.83. b Standardized effects of each factor on grass carp growth performance and pathogen resistance profiles calculated from the results of partial least squares path modeling.