Skip to main content

Enterococcus faecium are associated with the modification of gut microbiota and shrimp post-larvae survival



Probiotics are widely used to promote host health. Compared to mammals and terrestrial invertebrates, little is known the role of probiotics in aquatic invertebrates. In this study, eighteen tanks with eight hundred of shrimp post-larvae individuals each were randomly grouped into three groups, one is shrimps administered with E. faecium as probiotic (Tre) and others are shrimps without probiotic-treatment (CK1: blank control, CK2: medium control). We investigated the correlations between a kind of commercial Enterococcus faecium (E. faecium) powder and microbiota composition with function potentials in shrimp post-larvae gut.


We sequenced the 16S rRNA gene (V4) of gut samples to assess diversity and composition of the shrimp gut microbiome and used differential abundance and Tax4Fun2 analyses to identify the differences of taxonomy and predicted function between different treatment groups. The ingested probiotic bacteria (E. faecium) were tracked in gut microbiota of Tre and the shrimps here showed the best growth performance especially in survival ratio (SR). The distribution of SR across samples was similar to that in PCoA plot based on Bray-Curits and two subgroups generated (SL: SR < 70%, SH: SR ≥ 70%). The gut microbiota structure and predicted function were correlated with both treatment and SR, and SR was a far more important factor driving taxonomic and functional differences than treatment. Both Tre and SH showed a low and uneven community species and shorted phylogenetic distance. We detected a shift in composition profile at phylum and genus level and further identified ten OTUs as relevant taxa that both closely associated with treatment and SR. The partial least squares path model further supported the important role of relevant taxa related to shrimp survival ratio.


Overall, we found gut microbiota correlated to both shrimp survival and ingested probiotic bacteria (E. faecium). These correlations should not be dismissed without merit and will uncover a promising strategy for developing novel probiotics through certain consortium of gut microbiota.


The gut microbiota has emerged as a key regulator of host health, affecting digestion, immune activities, and homeostasis maintenance [1,2,3,4]. As such, the characterization and subsequent manipulation of this microscopic community is a supportive proposition for aquaculture research. Litopenaeus vannamei, a valuable and favored specie for many Asian countries [5], has subsequently dominated the field of aquatic invertebrate gut microbiomes [6]. With increasingly threatened by the onset of various diseases in shrimp aquaculture [7], the most common interventions, such as antibiotics [8, 9] and probiotics [10, 11], have been conducted to produce numbers of possible positive effects on the shrimp host. Considering the pernicious side effects of antibiotics and other chemotherapeutic agents, such as antibiotic resistance and the production of acute toxic to aquatic organisms [8], probiotics have become promising candidates for maintaining host health for its “green” and long-term research [6]. However, in comparison with mammals and terrestrial invertebrates, relatively very little is known about how probiotics impact gut microbiota of aquatic invertebrates.

It has been recognized that probiotics can increase competition, potentially supporting colonization resistance against pathogens in the gut, and thereby to promote shrimp survival and other growth performance [12]. For instance, shrimps treated with Lactobacillus. plantarum in their diet showed a significant increase in survival rate and growth, along with a reduction of Vibrio. harvey that is a common cause of shrimp mortalities [10]. The probiotic consortium has demonstrated an effective way to reduce the abundance of pathogenic Vibrio species and to prevent the mortality during Vibrio challenges [13]. Such studies on shrimp have mainly focused on the physiological and immunological responses to the applied probiotics. Even there were a range of probiotics showed unstable effects or they were not beneficial to the host [14, 15]. Most importantly, probiotics must be able to pass through the gut and information about whether the ingested probiotics can be detected in the gut microbiota remains scarce and contradictory [16, 17]. To increase resolution and enhance the likelihood of identifying ecotypes of ingested probiotic bacteria, a novel supervised computational pipeline, named oligotyping pipeline, has been further conducted [18].

Previous studies in human model focusing on the microecological effects of probiotics on gut microbiota were conflicting. Some studies disclosed that probiotics may perturb the function of indigenous microbiota [19, 20], whereas Kristensen et al. [21] argued that probiotics may promote the homeostasis of the gut microbiota without altering its composition. Regarding to shrimp, it is often suggested that a reduction in bacterial diversity (Shannon and Richness) within the gut or the differential abundance of certain bacteria may be responsible for the onset of pathogenesis [22, 23]. However, lower species diversity was also in healthy one [24]. Since many debates over these still exist and it is often impossible to discern between cause and effect, related research is still going on. Furthermore, the significance of shaping gut microbiome in larvae stage or early has been recognized recently [25]. The initial colonizers of the gut microbiota may act as a natural barrier to subsequent bacterial colonization and thus have long-term effects on host development [22, 26]. However, several probiotic species that are not indigenous to the marine environment with limited proliferation potential are known to be bad colonizers [27]. Given the host-specific gut microbiota, exploration of its composition and function potentials coupled with administrated probiotic bacteria will extend our knowledge of their relationships in how to regulate host health.

Modular co-occurrence network has been used to provide novel insights into potential interactions by revealing the niche spaces shared by modular members [28]. The clustering of microbial species into distinct modules has been used to infer gathering mechanism [29, 30]. In this study, we used 16S rRNA gene (V4 region) amplicon sequencing and Tax4Fun2 with six biological replicate tanks to measure gut microbiota and predicted function of post-larval shrimp. We compared the difference of gut microbiota composition and function potentials between shrimps administered with Enterococcus faecium as probiotic and those without probiotic-treatment, and thereby to address three key questions: (i) Can administrated probiotics be identified in the shrimp gut after 51-day continuous feeding? (ii) How does this probiotic impact the gut microbiota composition and predicted function? (iii) What are the relationships between gut microbiota and shrimp survival?


Shrimp performance and environmental factors

During shrimp culture, there were no significant variation in environmental factors, such as water temperature, dissolved oxygen (DO), pH, ammonium, and nitrite across all groups (Additional file 1: Fig. S1). At the end of cultivation, we calculated growth performance of all survivors each tank (Additional file 1: Table S1). The survival ratio (SR) is recognized as a crucial indicator for shrimp in post-larvae or before. We observed that the survival ratio of shrimp with probiotic administration (Tre) was significantly higher than that of the shrimp without probiotic administration (CK1, blank control with additive-free feed; CK2, medium control with medium-treated feed) (Table 1, P < 0.05). Notably, a larger difference of survival ratio was found within shrimps of the controls, while shrimps of Tre showed a more consistent survival pattern (Additional file 1: Fig. S2). The other indexes characterizing shrimp growth performance, such as unit yield (UY) and specific growth ratio (SGR) of shrimp was also higher in Tre. Among these, UY of Tre was significantly higher than CK2 (P < 0.05), and FCR of Tre was significantly lower than CK2 (P < 0.05), but both were non-significant compared to CK1 (Table 1).

Table 1 Survival ratio (SR), unit yield (UY), and specific growth ratio (SGR) of shrimp and the feed conversion ratio (FCR) after 51 days

Tracking probiotic bacteria in gut microbiota

In order to increase the amount of effective probiotic bacteria, activation is the regular step before its application into aquaculture [14]. After activation, the relative abundance of Enterococcus significantly increased from 24.29% in E. faecium powder (EF-P) to 87.42% in activated E. faecium (A-EF), and then sharply dropped to 10.85% in E. faecium-fermented feeds (EF-F) (Fig. 1, Additional file 1: Table S2). After 51 days, Enterococcus, which accounted for 0.5% relative abundance of all samples, was the only genus from EF-P that could be traced in the Tre gut. Oligotyping pipeline [18] was conducted to finely track the source of Enterococcus in the gut microbiota. The OTUs were limited to those annotated as Enterococcus in our dataset. Two information-rich positions of the aligned 16S rRNA V4 sequences were identified by Shannon entropy analysis. Overall, Enterococcus comprised 12 major oligotypes that were almost exclusively present in gut samples of Tre and EF-P samples (Additional file 1: Fig. S3). Up to 92% of these sequences were represented by OTU 1462 (Enterococcus), which was identified as E. faecium except oligotype “CT” from the EzBioCloud database [31]. Moreover, OTU1462 was the most abundant taxa affiliated to Enterococcus in EF-P (71.19%), A-EF (86.36%), and EF-F (82.85%), and it was mostly found in gut microbiota of Tre (Additional file 1: Fig. S4).

Fig. 1
figure 1

Experimental design and the dynamic compositional variation of Enterococcus faecium powder (EF-P) (bar chart, abundant (average relative abundance > 1%) genera in EFP)

The diversity analysis of Gut microbiota composition

Principal coordinate analysis (PCoA) based on Bray–Curtis distance was carried out to study the effect of treatment on the gut microbiota structure (Fig. 2A), and objects driven by treatment didn’t exhibit evident clusters. However, there was a shorter distance from objects to centroid in Tre compared with the controls (Fig. 2B), which was supported by ANOSIM analysis that showed samples outside Tre group are more different than those within Tre group (Table 2). Unexpectedly, a clearly more-consistent cluster with 32.75% explained variance at the fist dimension was observed (Fig. 2A). We tried to map these objects to SR for their similar separation pattern of samples (Additional file 1: Fig. S2), and thus two subgroups bounded by 70% SR were generated; one is tanks with less survivors (SL, SR < 70% with total 24 samples), another is tanks with more survivors (SH, SR ≥ 70% with total 30 samples). Changes in the gut microbiota structure occurred between SL and SH were also supported by ANOSIM (Table 2). Furthermore, a PERMANOVA test was applied to identify the most important factors driving the structure differences between gut microbiomes. The results identified that treatment combined with SR (F test = 7.28, R2 = 0.43, P < 0.001) explained 43% of the between-sample variance. Specifically, SR (F test = 3.40, R2 = 0.12, P < 0.001) was a more important factor driving the structure differences than treatment (F test = 19.33, R2 = 0.27, P < 0.001). Meantime, vectors were loaded to estimate the contribution of two factors, and the result also emphasized that SR contributes far more to the gut microbiota structure than treatment (Fig. 2A).

Fig. 2
figure 2

Diversity of gut microbiota. A Principal coordinates analysis plots of gut microbiota Bray–Curtis beta diversity for treatment and survival ratio (SR). The dotted ellipses with corresponding color for treatment represent the 95% confidence interval. The hollow circle represents shrimp survival ratios of ≤ 70% (SL), while the solid circle represents shrimp survival ratios of > 70% (SH). B Boxplot of Bray–Curtis distance between samples and group centroids (**P < 0.01, Tukey's 'Honest Significant Difference'). C Alpha-diversity indexes of gut microbiota for treatment (*P < 0.05, **P < 0.01, Mann–Whitney U test). D Alpha-diversity indexes of gut microbiota for SR in each treatment group. The hollow represents SL, while the solid represents SH (*P < 0.05; **P < 0.01; Mann–Whitney U test)

Table 2 Analysis of similarity (ANOSIM) of the difference in gut microbiota grouped by treatment and the survival ratio for each group, based on Bray–Curtis distance

Regarding to bacterial α-diversity, the community of Tre showed lower values than that observed in the controls (Fig. 2C), indicating a less species, more similar taxonomy, albeit uneven species in the gut microbiota of Tre. However, significant differences (P < 0.05) were only observed between Tre and CK2, not between Tre and CK1 (Fig. 2C). For grouping by SR in each group, only CK2 showed a slightly higher community species in SH compared with that observed in SL (Fig. 2D). Both Shannon and Evenness indexes showed significantly lower values in SL than that of SH (P < 0.01). Although not statistically, the values of PD also showed a certain difference between SH and SL (Fig. 2D). In addition, the community-level habitat niche breadths were further estimated to determine the range of resources a species uses. A relatively wider niche was observed in the gut microbiota of Tre (Additional file 1: Fig. S5A), but they were non-significant. The gut microbiota of SH exhibited significantly higher niche breadth values than that of SL (P < 0.001) (Additional file 1: Fig. S5B).

The difference in community composition among shrimp categories

The distribution of the main bacterial phyla/genera characterizing the gut microbiota were related to treatment and SR in the study (Additional file 1: Fig. S6). Compared with two control groups, Proteobacteria (CK1: 35.70%, CK2: 37.31, Tre: 27.09%) was subjected to a significant reduction in the gut microbiota of Tre where Bacteroides (CK1: 26.55%, CK2: 25.51%, Tre: 34.44%) accounted a major proportion (Additional file 1: Table S3 and Fig. S6A). Regarding to SR subgroups, the major phyla distribution in the gut microbiota of SH were the same as that of Tre (Additional file 1: Table S3 and Fig. S6A). Besides, the relative abundances of Planctomycetes, Chloroflexi, Tenericutes were significantly higher in SL compared with SH, whereas Verrucomicrobia and Saccharibacteria demonstrated an opposite trend (P < 0.05, Additional file 1: Table S3 and Fig. S6A).

At the genus level, the relative abundance of Algoriphagus was more differentially abundant in gut microbiota of Tre than that of two control groups while Vibrio was observed enriched in gut microbiota of CK1 and CK2 (Additional file 1: Table S3 and Fig. S6B). For grouping by SR, the relative abundances of Planctomyces, Haliea, Vibrio, Mycobacterium, Rhodopirellula, Ruegeria, Winogradskyella, Candidatus Bacilloplasma, and Nitrosomona were significantly higher in SL compared with that of SH. On the contrary, the gut microbiota of SH significantly enriched Haloferula, Formosa, Algoriphagus and Pirellula (P < 0.05, Additional file 1: Table S3 and Fig. S6B).

Differentially abundant OTUs were further identified by ANCOM analysis (Additional file 1: Fig. S7). OTU1462 was recognized as the applied probiotic bacteria according to oligotyping pipeline and was enriched in gut microbiota of Tre (Additional file 1: Fig. S4). Corresponding to this result, several OTUs including OTU1203 (Flavobacteriaceae), OTU996 (Vibrio), OTU1462 (Enterococcus), etc. were observed enriched in gut microbiota of Tre. However, it is worth mentioning that higher relative abundance of OTU996 (Vibrio) here was mainly attributed to gut samples from low-survival shrimp that contributed 95.13% (Additional file 1: Fig. S7). Significantly less abundant OTUs, for instance, OTU1189 (Vibrio), OTU879 (Phyllobactereriaceae) and OTU368 (Erythrobacter) were observed in gut microbiota of Tre (Additional file 1: Fig. S7). Further, ANCOM analysis identified 44 differentially abundant OTUs based on SR (Fig. 3). Among these, 23 OTUs with higher abundance in SH were mainly members of families Flavobacteriaceae (8/23), Halieaceae (5/23), and Verrucomicrobiaceae (4/23). However, the taxon distribution of 21 OTUs with higher abundances in SL was relatively more dispersed and was mainly attributed to families Planctomycetaceae (2/21), Rhodobacteraceae (2/21), and Microbacteriaceae (2/21). Among them, two potential pathogens, OTU996 (Vibrio) and OTU1469 (Chlamydiales Incertae Sedis) were found (Fig. 4).

Fig. 3
figure 3

Heatmap of differentially abundant OTUs for shrimp survival (SR) using ANCOM with custom value. The purple bar represents the OTUs were more prevalent in the gut of shrimp who survived less well, whereas the pink bar represents OTUs were more prevalent in the gut of shrimp who survived well

Fig. 4
figure 4

Diversity of predicted function (Kos) in gut microbiome. A Alpha diversity (Shannon) of predicted function in gut microbiota (* P < 0.05; ** P < 0.01; Mann–Whitney U test). B Principal coordinates analysis plots of predicted function Bray–Curtis beta diversity for treatment and survival ratio (SR). The dotted ellipse border for SL and the solid ellipse border for SH both represent the 95% confidence interval. The hollow circle represents shrimp survival ratios of ≤ 70% (SL), while the solid circle represents shrimp survival ratios of > 70% (SH). C Relationships between function dissimilarity and Shannon in gut microbiota. D Relationships between function dissimilarity and community distance in gut microbiota

Predicted functions of gut microbiota

Tax4Fun2 prediction was applied to assess how treatment and SR affect the potential function and pathway of gut microbiota and 7843 Kos generated based on 3553 OTUs. Unexpectedly, treatment didn’t shape the within-sample diversity (Shannon) of predicted function (P > 0.05). Otherwise, the values of Shannon index showed a significant difference between SL and SH (P < 0.05, Fig. 4A), which supported by PCoA where objects grouped by SR showed an almost clear separation based on Bray–Curtis distance (Fig. 4B). Consistently, PERMANOVA indicated that treatment and SR explained 9.04% and 28.63 variances in predicted function, respectively. In addition, function dissimilarity presented significant correlation with Shannon index (R2 = 0.6035, P < 0.001) and community distance (R2 = 0.4605, P < 0.001), indicating the more species, the more similar the predicted function composition. Furthermore, 110 significantly abundant pathways (organismal systems, cellular process, metabolism, etc.) were identified using Monte Carlo multiple randomized approach at 95% confidence interval (Fig. 5). The data was standardized by “range” and showed a clear cluster enriched in SH (Cluster A) (Additional file 2). Regarding to metabolism in cluster A, these were mainly involved amino acid metabolism, biosynthesis of other secondary metabolites, energy metabolism and metabolism of terpenoids and polyketides (Additional file 2). These results indicated that close interconnection between predicted function and gut microbiota.

Fig. 5
figure 5

Heatmap of differentially abundant Kos and the relative abundance of Ko was normalized by “range”

Ecological network analysis on the gut microbiota

A metacommunity network was constructed to estimate the bacterial interaction based on Spearman’s rank correlation (Additional file 1: Fig. S8). The entire network was modularized, and top eight modules (with 65.7% node coverage) were studied (Additional file 1: Fig. S8A). Out of the total significant co-occurrent features, 17.96% were allocated in module I, 15.29 in module II, etc. (Fig. S8B). Most modules were almost specific (relatively more abundant) to a single group, either to two control groups (Additional file 1: Fig. S8B). In addition, OTUs from module I were mainly affiliated to families Rhodobacteraceae, Vibrionaceae, and Microbacteriaceae, whereas OTUs from module II were primarily affiliated to families Halieaceae, Flavobacteriaceae, and Verrucomicrobiaceae.

The network for each subgroup was separately generated to explore their specific co-occurrence pattern (Fig. 6). The topological properties of each subgroup showed more diverse co-occurrence pattern in treatment group while SL and SH showed a more similar pattern (Fig. 6, Additional file 1: Fig. S9, Table S5). Specifically, the network of Tre showed the highest values of degree and eigenvector centrality compared to the two control groups. The network of CK2 showed the highest value of betweenness centrality and extremely low value of closeness centrality (Additional file 1: Fig. S9A). For SR subgroup, the nodes in network of SH with significantly higher closeness centrality and lower eigenvector centrality compared with that of SL (Additional file 1: Fig. S9B). Considering the high level of modularity in these network (Additional file 1: Table S5), we further analyzed the taxonomy of nodes in top three modules that account for 19.29%, 6.41%, 16.47%, 12.46% and 20.27% in the network of CK1, CK2, Tre, SL and SH, respectively (Additional file 3). OTUs from module 51, 143, 121 in CK1 were mainly affiliated to families Planctomycetaceae, Rhodobacteraceae, Flavobacteriaceae, and OTUs from module 77, 116, 76 in CK2 mainly affiliated to families Flavobacteriaceae, Halieaceae. As for Tre, OTUs from module 4, 2, 44 mainly affiliated to families Marinilabiaceae, Vibrionaceae, Rhodobacteraceae. In SL, OTUs from module 61, 35, 34 mainly affiliated to families Rhodobacteraceae, Planctomycetaceae, while most of the OTUs affiliated to family Flavobacteriaceae in SH (Additional file 3).

Fig. 6
figure 6

Co-occurrence patterns of gut microbiota for subgroup. There are 156, 159, 152, 149 and 82 modules in CK1, CK2, Tre, SL and SH, respectively. The top 3 modules with more OTUs are marked in each network. The size of each node is proportional to the number of relative abundances. The nodes are colored according to phylum

Relevant taxa related to shrimp survival

Using Boruta analysis, 10 OTUs that were possibly related to probiotics and shrimp survival were identified and defined as relevant OTUs (Additional file 4). To better analyze the relationship among applied probiotic bacteria (OTU 1462), relevant OTUs, gut microbiota, and shrimp survival, a partial least squares path model (PLS-PM) was constructed. This model was assessed using the Goodness of Fit (GoF) statistic, and the GoF value of which was 0.534. Unexpectedly, probiotic bacteria exerted a direct and weak association with gut microbiota (0.033) and shrimp survival (− 0.202) but showed a strong association with relevant OTUs (− 0.409) (Fig. 7). Relevant OTUs showed a strong and direct association with gut microbiota (0.965) and exhibited a non-significant association with shrimp survival (0.729) with a combination of direct (− 0.130) and indirect (0.859) associations. As expected, gut microbiota exerted a significant and direct association with shrimp survival (0.890) (Fig. 7). In this model, the average redundancy for shrimp survival represents that probiotic bacteria, relevant OTUs and gut microbiota predict 74% of the variability of survival indicator (Fig. 7).

Fig. 7
figure 7

The final path of partial least squares path model (PLS-PM) quantification of the direct and indirect effects of probiotic bacteria and relevant OTUs on gut microbiota, based on Bray–Curtis distance and shrimp survival. Blue and red lines represent positive and negative effects, respectively. Dashed lines indicate the coefficients that did not significantly (P > 0.05) differ from 0. The wider the line, the higher the absolute value of the coefficient obtained with PLS-PM. (** P < 0.01; *** P < 0.001)


Enterococcus faecium powder, a popular probiotic agent in aquaculture, was used to study the correlations between probiotics and gut microbiota. In practice, probiotics are activated to increase the number of effective bacteria before application. In this study, the relative abundance of Enterococcus greatly varied after activation and a 9-h fermentation with feeds (Fig. 1). The effects of feed fermentation on the gut microbial community varied with substrate, and inappropriate methods may counteract the beneficial effects of probiotics [14, 32]. Although not all of them require such operation, more attention should be paid to probiotic fermentation method especially the specific requirements and appropriate culturing conditions of target probiotics [33].

Although we could not determine whether E. faecium would be a good gut colonizer without supported result from Fluorescence in situ Hybridization (FISH), it was indeed tracked in gut microbiota with a total 0.5% abundance. Previous studies suggested that probiotic E. faecium could be isolated from healthy shrimps with antibacterial and adhesive activities [34, 35]. Although stable colonization is not necessarily required for probiotics to be beneficial to host health, a substantial population must, at least transiently, establish for them to have a metabolic impact on the host and the indigenous gut microbiota [36]. However, the strong correlation between probiotic and shrimp survival along with a small within-group difference in Tre was observed (Tables 1 and 2, Additional file 1: Fig. S2). One possible supported explanation for this is the critical role of rare species in regulating the bacterial community, which has been widely demonstrated in various natural ecosystems [37, 38] and host-associated microbiomes [39]. Notably, it may vary greatly depending on the niche and the species that is introduced, the time it is analyzed and the immune response in case of animals. In addition, the shrimp post-larvae used in this study may contribute to the prominent effect of such few proportions of Enterococcus. The development of microbiota within the gut is reported to be influenced by priority effects, in which the order and timing of species arrival determine how the species affect each other [40]. Previous studies proposed that probiotics may show a weaker effect on older than younger individuals [41,42,43]. Hence, we inferred that probiotic administration should be adopted sufficiently early to achieve the optimal effect and further study are required to support it.

Changes to the gut microbiota have been implicated in a wide range of host health in humans [44], as well as in shrimps [2, 22, 23, 45]. Specifically, the “healthy” gut microbiota was often characterized by high species diversity and function redundancy, less potential pathogens. The correlations between “healthy” gut microbiota and these “healthy” features should not be dismissed without merit and it revealed a general pattern in kinds of animals. However, diversity is the question, not the answer [46]. In our study, gut microbiota of Tre and SH indicated a simplified species community (lower values of Shannon and Richness) with a cohort of dominant taxa (lower value of Evenness) (Fig. 2C, D), which was also observed in gut microbiota of probiotic-treated fish [47]. The predicted function of a metacommunity appeared to be more unitary relative to the diverse community species (Fig. 4C, D), indicating the whole community may exist stronger function redundancy. In that case, we can infer that even if the species diversity of the gut microbiota of Tre and SH is low, their function potentials may be sufficient to maintain the whole ecosystem. Some predicted metabolic pathways, such as amino acid metabolism and energy metabolism, were enriched in SH group (Fig. 5), playing a key role in energy homeostasis [48]. In addition, it also can be well supported by a homogeneous composition structure (Fig. 2B) and a shorter phylogenetic distance (Fig. 2C, D) that were observed in Tre and SH group, suggesting that close microbial lineages may display ecological similarities, potentially accounting for the community cluster [49, 50].

Most pathogens that impact aquaculture are deemed opportunistic, including Vibrio. Previous studies have reported that was closely related to the onset of various diseases and mass mortality of shrimp [10, 11, 13]. As expected, we observed significantly abundant Vibrio in gut microbiota of shrimps who did not survive as well (Additional file 1: Table S3). However, those shrimps in SL with more abundant Vibrio were still survivors, that is to say, Vibrio strains can co-exist with other bacteria even with those ingested probiotic bacteria (Additional file 1: Fig. S7). Once dysbiosis in aquaculture, they will apparently exploit chance to emerge as virulent and damaging [51]. In that case, the antibacterial activity of Enterococcus strains against Vibrio may be shown [35, 52], which was supported by the result of this study that Vibrio was less abundant in gut microbiota of probiotic-treated shrimp (Additional file 1: Table S3). Otherwise, the modularity of ecological networks is thus essential to uncover how microbiome assembly influences microbiome structure [53]. The module 4 with highly modulization in gut microbiota of Tre showed members of genus Vibrio were more co-occurred with other taxa (Fig. 6, Additional file 3), which could infer that probiotic was associated with the more interactions occurred between certain taxa and Vibrio species. Besides, those centrality metrics of network grouped by treatment seemed more closely related to total graph sizes, for example, closeness is essentially a global centrality version of degree [30]. Gut microbiota of Tre showed a more co-occurrence within the community, which can be inferred that there are more underlying interactions that associated with probiotics.

Finally, we identified 10 relevant OTUs that were differentially abundant and associated with surviving individuals (Additional file 4). Among them, some members of Flavobacteriaceae are known for their ability to degrade and metabolize refractory organic compounds [54]. Certain planktonic Verrucomicrobiaceae phylotypes have the potential to be active polysaccharide degraders [55, 56], promoting food utilization to achieve enhanced growth performance. Additionally, several genera of Verrucomicrobiaceae have been isolated from marine animals, which might be closely related to the host health. Halieaceae, a newly found family, has been reported its metabolic capabilities, such as anoxygenic photosynthesis [57]. However, the PLS-PM model showed that these relevant OTUs strongly correlated with applied probiotic and the modification of gut microbiota and shrimp health (Fig. 7). Large-scale studies have demonstrated that probiotics indeed cause shifts in the compositions of gut microbiota [14, 58, 59]. However, few reports have proposed that probiotics play growth-promoting roles in shrimps in collaboration with gut “partners” in shrimps. This relationship may explain why multi-strain probiotics appear to show greater efficacy than single‐strain probiotics [60]. Notably, uncultured bacteria accounted for a considerable proportion of the indicators in this work (Additional file 4), which may indicate that they play a vital role in host health, with unexplored function. In the future, it will be more conducive to accurately screen host-specific taxa at the strain level. Deeper verification tests based on metagenomics and culturomics are necessary, which will be further applied in the probiotic field and provide insights into the relationship between interspecies interaction and shrimp survival.

Materials and methods

Experimental design and sample collection

The post-larval shrimp (Litopenaeus vannamei) tanks in this study are located in Zhejiang Mariculture Research Institute, Wenzhou, China (27° 51′ N, 120° 50′ E). All test tanks were equal in size (600 L) with 400 L aquaculture water. Both the daily water exchange rates and feed amount (Xiamiaobao, Chia Tai Group, China) were uniformly managed. Shrimps (length, 5 mm) were cultured in 18 tanks with 800 individuals each and were randomly divided into three groups (two control groups and one treatment group) with six replicates. Commercial Enterococcus faecium powder (EF-P, number of live bacteria ≈ 100 × 108 CFU/g, Yichun Strong Microbial Technology Co., Ltd., Jiangxi, China) was selected for the treatment group (Tre). The blank control (CK1) was provided normal feeds, while the medium control (CK2) was provided feeds mixed with sterilized de-Man Rogosa Sharpe (MRS) medium (MM-F, 5 mL of MRS medium per gram of feed). Tre was provided feeds mixed with activated commercial E. faecium powder (A-EF) at the same ratio three times a day. Detailed information on the activation, fermentation, and experimental setup are shown in Fig. 1.

The shrimps were starved for one day before sampling. After 51 days, we collected all shrimps and survival ratio (SR), unit yield (UY), specific growth rate (SGR), and feed conversion ratio (FCR) were estimated. Ten gut samples of large, medium, small shrimps including digesta were randomly collected as one biological replicate respectively to eliminate bias and a total of three biological replicates each tank was performed for DNA extraction. Importantly, there was no obvious disease or substantially abnormal death rates observed during cultivation management. Water temperature, dissolved oxygen (DO), oxidation–reduction potential (ORP), salinity and pH were measured in each tank every day with a probe (YSI 550A, USA). The concentrations of nitrite (NO2), nitrate (NO3), phosphate (PO43−), and ammonium (NH4+) were analyzed according to standard methods (GB 17378.4-2007) at day 11, 41, and 51.

DNA extraction and bacterial 16S rRNA gene amplification

The microbial DNA of shrimp gut, EF-P, and products of activation (A-EF) and fermentation (EF-F) were extracted according to the user’s manual of the QIAamp DNA Stool mini kit (Qiagen GmbH, Hilden, Germany). The concentrations of extracted DNA were measured on a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA). The bacterial primers 515F (5′–GTGCCAGCMGCCGCGGTAA–3′) and 806R (5′–GGACTACNNGGGTATCTAAT–3′) were used to amplify the V4 region of the 16S rRNA gene. Each sample was amplified in triplicate (50-μL reaction system) to reduce bias under the listed conditions: 25 cycles of denaturation at 95 °C for 30 s, annealing at 55 °C for 30 s, and extension at 72 °C for 45 s, with a final extension for 10 min. To check the expected amplification effect, PCR products were visualized on a 1.5% agarose gel. Triplicate amplicons were combined and purified by a PCR fragment purification kit (TaKaRa Biotech, Japan). The purified amplicons were quantified with a PicoGreen-iT dsDNA Assay Kit (Invitrogen, Carlsbad, CA, USA). Equal concentrations of amplicons were combined into one pooled sample and sequenced on a MiSeq PE 2 × 250 platform (Illumina, Majorbio, Shanghai, China).

Processing of sequencing data

Raw data were processed with QIIME v1.9.1 pipeline [61]. Using the UCHIME algorithm, chimeric sequences were identified [62] and removed from subsequent analysis. Bacterial phylotypes were identified using UCLUST [63] and were classified into operational taxonomic units (OTUs) at a 97% similarity level. For each OTU, the most abundant sequence was regarded as the representative sequence and was then taxonomically assigned using the SILVA_128 database ( The matching results yielded 2,702,120 reads, and 26,002 reads were unclassified due to chimeras. Archaea and Chloroplast sequences were removed using “”, as well as other sequences that could not be assigned to bacteria. Afterwards, singletons (n = 2) were also discarded. To correct for varying sampling efforts, the data were randomly rarefied at the minimum sequencing depth (24,000 sequences). After that, the sequencing of 16S rRNA genes generated 1,296,000 high-quality sequences and 7656 OTUs based on 97% similarity. Finally, OTUs in no less than two samples were selected for β-diversity and function analysis and 3553 OTUs were generated.

Bioinformatic analysis

QIIME was used to calculate the α-diversity indexes, including observed species (Richness), Shannon index, and phylogenetic diversity (PD). Pielou’s evenness (Evenness) was estimated with the Vegan R package. Nonparametric test (Mann–Whitney U test for two groups and Krustkal-Wallis H test for more than three groups) were performed on all data using SPSS v.16.0. Significance levels were considered at α < 0.05 with 95% confidence interval. Based on Bray–Curtis distances, principal coordinates analysis (PCoA) was applied to analyze the difference in gut microbiota and two factors (treatment and survival ratio) as dependent variables that are explained by the ordination scores were added using “envfit” in Vegan R package. The “betadisper” in Vegan R package was used to check the assumption of heterogeneity in dispersions according to Bray–Curtis distance [64]. Pairwise analyses of similarity (ANOSIM) and permutational multivariate analysis of variance (PERMAVONA) were applied to separate and quantitatively estimate the effects of treatment and survival ratio on the composition of gut microbiota by using “anosim” and “adonis” in Vegan R package [65]. We expected bacterial population with wider niche breadth at the community level to be more metabolically flexible [66], and the analysis was conducted using the “niche.width” in spaa R package [67]. Functional profiling of bacterial taxa was carried out using the Tax4Fun2 R package [68]. Shannon index of by KEGG orthology (KO) was measured. Additionally, PCoA and PERMAVONA were used to estimate the effects of treatment and on bacterial function represented by KO matrix based on Bray-Curits distance. Linear regression models were applied to estimate for significant correlations between Shannon index or community distance (PCoA1) or and function dissimilarity (PCoA1). Analysis of composition of microbiomes (ANCOM) was applied to detect the differentially abundant OTUs in the groups based on treatment and survival ratio [69], and 44 OTUs were further identified by custom values, with significant changes in their relative abundances between SH and SL in each group. The Boruta R package was used to find the relevant OTUs that caused differences in the gut microbiota of six groups (treatment × survival ratio) [70]. To account for the total effects of probiotic bacteria (E. faecium) and relevant OTUs on both gut microbiota and shrimp survival, partial least squares path modeling (PLS-PM) was run using functions available in the plspm R package.

Oligotyping analysis

The python-pipeline v.2.7 (available from was used for oligotyping that uses highly variable single-nucleotide polymorphisms to achieve strain-level resolutions. After the initial calculation of Shannon entropy (variability) across each nucleotide position with the analyze-entropy script, members of genus Enterococcus were done, and the process were supervised until each oligotype had converged [18]. Considering the small number of sequences, we kept those almost resulting in oligotypes with regards to the detailed parameters: c (2), s (1), A (0), M(0), which specified that a unique oligotype represented at least 1% of sequence reads, was detected in at least 1 samples, and showed all unique read.

Network analysis

To reduce the complexity and increase the credibility of the data sets, OTUs in no less than one-third of the samples were selected to construct the whole network and subgroups. All possible pairwise Spearman’s rank correlations (r) between the selected OTUs were estimated using the Picante R package. Only extraordinarily robust (r > 0.8 or r <  − 0.8) and statistically significant (P < 0.001) correlations were incorporated into the network analyses. All the networks exhibited scale-free characteristics, indicating non-random network structures. Network visualization, node-level topological features (i.e., degree, betweenness, closeness centrality, and eigenvector centrality), and modular analysis were performed in Gephi v.0.8.2. Afterward, 1000 Erdös–Réyni random networks, assigned the identical number of nodes and edges as the real metacommunity network, were generated within the igraph R package [71], where the topology characteristics, including modularity, clustering coefficient, and average path length, of both real and random networks were calculated and compared.

Availability of data and materials

The raw sequencing files from this study are available from the NCBI Sequence Read Archive (SRA) under accession number PRJNA591260.



De-Man Rogosa Sharpe


Unit yield


Survival ratio


Feed conversion ratio


Specific growth rate


Dissolved oxygen


Deoxyribonucleic acid


Ribosomal ribonucleic acid


Polymerase chain reaction


Operational taxonomic unit


Quantitative insights into microbial ecology


Principal coordinate analysis


Analysis of similarity


Analysis of composition of microbiomes


Partial least squares path modeling


  1. 1.

    Yang H, Yang M, Sun J, Shi X, Zhao X, Wang J. Dual oxidases participate in the regulation of intestinal microbiotic homeostasis in the kuruma shrimp Marsupenaeus japonicus. Dev Comp Immunol. 2016;59:153–63.

    CAS  PubMed  Google Scholar 

  2. 2.

    Xiong J, Dai W, Zhu J, Liu K, Dong C, Qiu Q. The underlying ecological processes of gut microbiota among cohabitating retarded, overgrown and normal shrimp. Microb Ecol. 2017;73:988–99.

    PubMed  Google Scholar 

  3. 3.

    Dai W, Yu W, Zhang J, Zhu J, Tao Z, Xiong J. The gut eukaryotic microbiota influences the growth performance among cohabitating shrimp. Appl Microbiol Biot. 2017;101:6447–57.

    CAS  Google Scholar 

  4. 4.

    Dai W, Zhang J, Qiu Q, Chen J, Yang W, Ni S, Xiong J. Starvation stress affects the interplay among shrimp gut microbiota, digestion and immune activities. Fish Shellfish Immun. 2018;80:191–9.

    Google Scholar 

  5. 5.

    Liao I, Chien Y. The pacific white shrimp, Litopenaeus vannamei, in Asia: the world’s most widely cultured alien crustacean. In: The wrong place-alien marine crustaceans: distribution, biology and impacts. Springer, Dordrecht. 2011; 489–519.

  6. 6.

    Holt C, Bass D, Stentiford G, van der Giezen M. Understanding the role of the shrimp gut microbiome in health and disease. J Invertebr Pathol. 2020;186:107387.

    PubMed  Google Scholar 

  7. 7.

    Thitamadee S, Prachumwat A, Srisala J, Jaroenlak P, Salachan PV, Sritunyalucksana K, et al. Review of current disease threats for cultivated penaeid shrimp in Asia. Aquaculture. 2016;452:69–87.

    Google Scholar 

  8. 8.

    Holmström K, Gräslund S, Wahlström A, Poungshompoo S, Bengtsson B, Kautsky N. Antibiotic use in shrimp farming and implications for environmental impacts and human health. Int J Food Sci Technol. 2003;38(3):255–66.

    Google Scholar 

  9. 9.

    Tu H, Silvestre F, Phuong N, Kestemont P. Effects of pesticides and antibiotics on penaeid shrimp with special emphases on behavioral and biomarker responses. Environ Toxicol Chem. 2010;29(4):929–38.

    CAS  PubMed  Google Scholar 

  10. 10.

    Kongnum K, Hongpattarakere T. Effect of Lactobacillus plantarum isolated from digestive tract of wild shrimp on growth and survival of white shrimp (Litopenaeus vannamei) challenged with Vibrio harveyi. Fish Shellfish Immunol. 2012;32(1):170–7.

    PubMed  Google Scholar 

  11. 11.

    Mazón-Suástegui J, Salas-Leiva J, Medina-Marrero R, Medina-García R, García-Bernal M. Effect of Streptomyces probiotics on the gut microbiota of Litopenaeus vannamei challenged with Vibrio parahaemolyticus. Microbiologyopen. 2020;9(2):e967.

    PubMed  Google Scholar 

  12. 12.

    Knipe H, Temperton B, Lange A, Bass D, Tyler C. Probiotics and competitive exclusion of pathogens in shrimp aquaculture. Rev Aquac. 2021;13(1):324–52.

    Google Scholar 

  13. 13.

    Quiroz-Guzmán E, Vázquez-Juárez R, Luna-González A, Balcázar JL, Barajas-Sandoval DR, Martínez-Díaz SF. Administration of probiotics improves the brine shrimp production and prevents detrimental effects of pathogenic vibrio species. Mar Biotechnol. 2018;20:512–9.

    Google Scholar 

  14. 14.

    Liu Z, Qiuqian L, Yao Z, Wang X, Huang L, Zheng J, et al. Effects of a commercial microbial agent on the bacterial communities in shrimp culture system. Front Microbiol. 2018;9:2430.

    PubMed  PubMed Central  Google Scholar 

  15. 15.

    Rubia A, Rahman M, Rahman BM, Ahmed KK. The effect of commercial probiotics on health and production of shrimp (Penaeus monodon). Asian J Med Biol Res. 2017;3:88–93.

    Google Scholar 

  16. 16.

    Suez J, Zmora N, Zilberman-Schapira G, Mor U, Dori-Bachash M, Bashiardes S, et al. Post-antibiotic gut mucosal microbiome reconstitution is impaired by probiotics and improved by autologous FMT. Cell. 2018;174(6):1406–23.

    CAS  PubMed  Google Scholar 

  17. 17.

    Zmora N, Zilberman-Schapira G, Suez J, Mor U, Dori-Bachash M, Bashiardes S, et al. Personalized gut mucosal colonization resistance to empiric probiotics is associated with unique host and microbiome features. Cell. 2018;174(6):1388–405.

    CAS  PubMed  Google Scholar 

  18. 18.

    Eren AM, Maignien L, Sul WJ, Murphy LG, Grim SL, Morrison HG, et al. Oligotyping: differentiating between closely related microbial taxa using 16S rRNA gene data. Methods Ecol Evol. 2013;4:1111–9.

    PubMed Central  Google Scholar 

  19. 19.

    Eloe-Fadrosh EA, Brady A, Crabtree J, Drabek EF, Ma B, Mahurkar A, et al. Functional dynamics of the gut microbiome in elderly people during probiotic consumption. mBio. 2015;6:e00231-15.

    PubMed  PubMed Central  Google Scholar 

  20. 20.

    Maldonado-Gómez MX, Martínez I, Bottacini F, O’Callaghan A, Ventura M, van Sinderen D, et al. Stable engraftment of Bifidobacterium longum AH1206 in the human gut depends on individualized features of the resident microbiome. Cell Host Microbe. 2016;20:515–26.

    PubMed  Google Scholar 

  21. 21.

    Kristensen NB, Bryrup T, Allin KH, Nielsen T, Hansen TH, Pedersen O. Alterations in fecal microbiota composition by probiotic supplementation in healthy adults: a systematic review of randomized controlled trials. Genome Med. 2016;8:52.

    PubMed  PubMed Central  Google Scholar 

  22. 22.

    Xiong J. Progress in the gut microbiota in exploring shrimp disease pathogenesis and incidence. Appl Microbiol Biot. 2018;102:7343–50.

    CAS  Google Scholar 

  23. 23.

    Yao Z, Yang K, Huang L, Huang X, Qiuqian L, Wang K, et al. Disease outbreak accompanies the dispersive structure of shrimp gut bacterial community with a simple core microbiota. AMB Expr. 2018;8:120.

    Google Scholar 

  24. 24.

    Zakrzewski M, Simms L, Brown A, Appleyard M, Irwin J, Waddell N, et al. IL23R-protective coding variant promotes beneficial bacteria and diversity in the ileal microbiome in healthy individuals without inflammatory bowel disease. J Crohns Colitis. 2019;13(4):451–61.

    PubMed  Google Scholar 

  25. 25.

    Wang Y, Wang K, Huang L, Dong P, Wang S, Chen H, et al. Fine-scale succession patterns and assembly mechanisms of bacterial community of Litopenaeus vannamei larvae across the developmental cycle. Microbiome. 2020;8:106.

    CAS  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Xiong J, Dai W, Li C. Advances, challenges, and directions in shrimp disease control: the guidelines from an ecological perspective. Appl Microbiol Biot. 2016;100:6947–54.

    CAS  Google Scholar 

  27. 27.

    Vargas-Albores F, Porchas-Cornejo M, Martínez-Porchas M, Villalpando-Canchola E, Gollas-Galván T, Martínez-Córdova L, et al. Bacterial biota of shrimp intestine is significantly modified by the use of a probiotic mixture: a high throughput sequencing approach. Helgol Mar Res. 2017;71(1):1–10.

    Google Scholar 

  28. 28.

    Faust K, Raes J. Microbial interactions: from networks to models. Nat Rev Microbiol. 2010;10:538–50.

    Google Scholar 

  29. 29.

    Jiao S, Yang Y, Xu Y, Zhang J, Lu Y. Balance between community assembly processes mediates species coexistence in agricultural soil microbiomes across eastern China. ISME J. 2020;14:202–16.

    PubMed  Google Scholar 

  30. 30.

    Mamet SD, Redlick E, Brabant M, Lamb EG, Helgason BL, Stanley K, et al. Structural equation modeling of a winnowed soil microbiome identifies how invasive plants re-structure microbial networks. ISME J. 2019;13:1988–96.

    PubMed  PubMed Central  Google Scholar 

  31. 31.

    Yoon S, Ha S, Kwon S, Lim J, Kim Y, Seo H, et al. Introducing EzBioCloud: a taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies. Int J Syst Evol Microbiol. 2017;67(5):1613.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Choi J, Shinde P, Ingale S, Kim J, Kim Y, Kim K, et al. Evaluation of multi-microbe probiotics prepared by submerged liquid or solid substrate fermentation and antibiotics in weaning pigs. Livest Sci. 2011;138:144–51.

    Google Scholar 

  33. 33.

    Wang C, Shi C, Zhang Y, Song D, Lu Z, Wang Y. Microbiota in fermented feed and swine gut. Appl Microbiol Biot. 2008;102:2941–8.

    Google Scholar 

  34. 34.

    Maeda M, Shibata A, Biswas G, Korenaga H, Kono T, Itami T, et al. Isolation of lactic acid bacteria from kuruma shrimp (Marsupenaeus japonicus) intestine and assessment of immunomodulatory role of a selected strain as probiotic. Mar Biotechnol. 2014;16(2):181–92.

    CAS  Google Scholar 

  35. 35.

    Sha Y, Wang L, Liu M, Jiang K, Xin F, Wang B. Effects of lactic acid bacteria and the corresponding supernatant on the survival, growth performance, immune response and disease resistance of Litopenaeus vanname. Aquaculture. 2016;452:28–36.

    CAS  Google Scholar 

  36. 36.

    Ojima M, Yoshida K, Sakanaka M, Jiang L, Odamaki T, Katayama T. Ecological and molecular perspectives on responders and non-responders to probiotics and prebiotics. Curr Opin Biotechnol. 2022;73:108–20.

    CAS  Google Scholar 

  37. 37.

    Zhang H, Hou F, Xie W, Wang K, Zhou X, Zhang D, et al. Interaction and assembly processes of abundant and rare microbial communities during a diatom bloom process. Environ Microbiol. 2020;22:1707–19.

    CAS  PubMed  Google Scholar 

  38. 38.

    Du S, Dini-Andreote F, Zhang N, Liang C, Yao Z, Zhang H, Zhang D. Divergent co-occurrence patterns and assembly processes structure the abundant and rare bacterial communities in a salt marsh ecosystem. Appl Environ Microb. 2020;86:e00322-e420.

    CAS  Google Scholar 

  39. 39.

    Jousset A, Bienhold C, Chatzinotas A, Gallien L, Gobet A, Kurm V, et al. Where less may be more: how the rare biosphere pulls ecosystems strings. ISME J. 2017;11:853–62.

    PubMed  PubMed Central  Google Scholar 

  40. 40.

    Sprockett D, Fukami T, Relman DA. Role of priority effects in the early-life assembly of the gut microbiota. Nat Rev Gastro Hepat. 2018;15:197–205.

    Google Scholar 

  41. 41.

    Gebru E, Lee J, Son J, Yang S, Shin S, Kim B, et al. Effect of probiotic-, bacteriophage-, or organic acid-supplemented feeds or fermented soybean meal on the growth performance, acute-phase response, and bacterial shedding of grower pigs challenged with Salmonella enterica serotype typhimurium. J Anim Sci. 2010;88:3880–6.

    CAS  PubMed  Google Scholar 

  42. 42.

    Giatsis C, Sipkema D, Smidt H, Heilig H, Benvenuti G, Verreth J, et al. The impact of rearing environment on the development of gut microbiota in tilapia larvae. Sci Rep. 2015;5:18206.

    CAS  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Yan Q, Li J, Yu Y, Wang J, He Z, Van Nostrand JD, et al. Environmental filtering decreases with fish development for the assembly of gut microbiota. Environ Microbiol. 2016;18:4739–54.

    CAS  PubMed  Google Scholar 

  44. 44.

    Claesson M, Clooney A, O’toole P. A clinician’s guide to microbiome analysis. Nat Rev Gastroenterol Hepatol. 2017;14(10):585–95.

    PubMed  Google Scholar 

  45. 45.

    Xiong J, Zhu J, Dai W, Dong C, Qiu Q, Li C. Integrating gut microbiota immaturity and disease-discriminatory taxa to diagnose the initiation and severity of shrimp disease. Environ Microbiol. 2017;19:1490–501.

    PubMed  Google Scholar 

  46. 46.

    Shade A. Diversity is the question, not the answer. ISME J. 2017;11:1–6.

    PubMed  Google Scholar 

  47. 47.

    Gupta S, Fečkaninová A, Lokesh J, Koščová J, Sørensen M, Fernandes J, et al. Lactobacillus dominate in the intestine of Atlantic salmon fed dietary probiotics. Front Microbiol. 2019;9:3247.

    PubMed  PubMed Central  Google Scholar 

  48. 48.

    Baniel A, Amato K, Beehner J, Bergman T, Mercer A, Perlman R, et al. Seasonal shifts in the gut microbiome indicate plastic responses to diet in wild geladas. Microbiome. 2021;9(1):1–20.

    Google Scholar 

  49. 49.

    Tapia-Paniagua ST, Chabrillón M, Patricia Díaz-Rosales P, de la Banda IG, Lobo C, Balebona MC, et al. Intestinal microbiota diversity of the flat fish Solea senegalensis (Kaup, 1858) following probiotic administration. Microb Ecol. 2010;60:310–9.

    PubMed  Google Scholar 

  50. 50.

    Philippot L, Andersson S, Battin T, Prosser JI, Schimel JP, Whitman WB, et al. The ecological coherence of high bacterial taxonomic ranks. Nat Rev Microbiol. 2010;8:523–9.

    CAS  PubMed  Google Scholar 

  51. 51.

    Derome N, Gauthier J, Boutin S, Llewellyn M. Bacterial opportunistic pathogens of fish. In: The Rasputin effect: when commensals and symbionts become parasitic. Cham: Springer; 2016. p. 81–108.

  52. 52.

    Swain S, Singh C, Arul V. Inhibitory activity of probiotics Streptococcus phocae PI80 and Enterococcus faecium MC13 against vibriosis in shrimp Penaeus monodon. World J Microbiol Biotechnol. 2009;25(4):697–703.

    Google Scholar 

  53. 53.

    Lurgi M, Thomas T, Wemheuer B, Webster N, Montoya J. Modularity and predicted functions of the global sponge-microbiome network. Nat Commun. 2019;10(1):1–12.

    CAS  Google Scholar 

  54. 54.

    Oberbeckmann S, Osborn AM, Duhaime MB. Microbes on a bottle: Substrate, season and geography influence community composition of microbes colonizing marine plastic debris. PLOS ONE. 2016;11:e0159289.

    PubMed  PubMed Central  Google Scholar 

  55. 55.

    Freitas S, Hatosy S, Fuhrman JA, Huse SM, Welch DB, Sogin ML, et al. Global distribution and diversity of marine Verrucomicrobia. ISME J. 2012;6:1499–505.

    CAS  PubMed  PubMed Central  Google Scholar 

  56. 56.

    Martinez-Garcia M, Brazel DM, Swan BK, Arnosti C, Chain PSG, Reitenga KG, et al. Capturing single cell genomes of active polysaccharide degraders: an unexpected contribution of Verrucomicrobia. PLOS ONE. 2012;7:e35314.

    CAS  PubMed  PubMed Central  Google Scholar 

  57. 57.

    Li SH, Seo JH, Kang I, Cho JC. Inferred metabolic capabilities of novel marine bacterial clades of the family Halieaceae, abstr, p-49. 2019 International meeting of the microbiological society of Korea.

  58. 58.

    Vargas-Albores F, Porchas-Cornejo MA, Martínez-Porchas M, Villalpando-Canchola E, Gollas-Galván T, Martínez-Córdova LR. Bacterial biota of shrimp intestine is significantly modified by the use of a probiotic mixture: a high throughput sequencing approach. Helgol Mar Res. 2017;71:5.

    Google Scholar 

  59. 59.

    Xie JJ, Liu Q, Liao S, Fang HH, Yin P, Xie SW, et al. Effects of dietary mixed probiotics on growth, non-specific immunity, intestinal morphology and microbiota of juvenile pacific white shrimp Litopenaeus vannamei. Fish Shellfish Immun. 2019;90:456–65.

    CAS  Google Scholar 

  60. 60.

    Chapman CMC, Gibson GR, Rowland I. Health benefits of probiotics: are mixtures more effective than single strains? Eur J Nutr. 2011;50:1–17.

    CAS  Google Scholar 

  61. 61.

    Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  62. 62.

    Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27:2194–200.

    CAS  PubMed  PubMed Central  Google Scholar 

  63. 63.

    Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.

    CAS  PubMed  Google Scholar 

  64. 64.

    Anderson M. Distance-based tests for homogeneity of multivariate dispersions. Biometrics. 2006;62(1):245–53.

    PubMed  Google Scholar 

  65. 65.

    Anderson M. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001;26(1):32–46.

    Google Scholar 

  66. 66.

    Wu W, Lu HP, Sastri A, Yeh YC, Gong GC, Chou WC, et al. Contrasting the relative importance of species sorting and dispersal limitation in shaping marine bacterial versus protist communities. ISME J. 2017;12:485–94.

    PubMed  PubMed Central  Google Scholar 

  67. 67.

    Zhang J. Spaa: Species association analysis. R package version 0.2.1,; 2013.

  68. 68.

    Wemheuer F, Taylor J, Daniel R, Johnston E, Meinicke P, Thomas T, et al. Tax4Fun2: prediction of habitat-specific functional profiles and functional redundancy based on 16S rRNA gene sequences. Environ Microbiome. 2020;15:1–12.

    Google Scholar 

  69. 69.

    Mandal S, Van Treuren W, White RA, Eggesbø M, Knight R, Peddada SD. Analysis of composition of microbiomes: a novel method for studying microbial composition. Microb Ecol Health Dis. 2015;26:27663.

    PubMed  Google Scholar 

  70. 70.

    Kursa MB, Jankowski A, Rudnicki WR. Boruta-A system for feature selection. Fund Inform. 2010;101:271–86.

    Google Scholar 

  71. 71.

    Erdös P, Réyni A. On the evolution of random graphs. Publ Math Inst Hung Acad Sci. 1960;5:17–60.

    Google Scholar 

Download references


The authors would like to thank the entire research group in Zhang lab of the Ningbo University.


The work was supported by National Natural Science Foundation of China (31672658), the Natural Science Foundation of Ningbo (2019A610445) and K.C. Wong Magna Fund in Ningbo University.

Author information




Designed the experiment: DZ and SD. Performed the experiments and sequencing: SD, WC, XH and CC. Analyzed the microbiome data: SD. Drafting the manuscript: SD. Contributed to writing the manuscript: ZY, HG, and DZ. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Zhiyuan Yao or Demin Zhang.

Ethics declarations

Ethics approval and consent to participate

All animal experiments were carried out in accordance with the UK Animals (Scientific Procedures) Act, 1986, and associated guidelines and the EU Directive 2010/63/EU for animal experiments.

Consent for publication

Not applicable.

Competing interests

The authors have no conflict of interest to declare.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

 Supplementary tables and figures.

Additional file 2:

 Ko pathways of cluster A.

Additional file 3:

 OTUs from major modules.

Additional file 4:

 Relevant OTUs.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Du, S., Chen, W., Yao, Z. et al. Enterococcus faecium are associated with the modification of gut microbiota and shrimp post-larvae survival. anim microbiome 3, 88 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Shrimp post-larvae
  • Survival
  • Enterococcus faecium
  • Gut microbiota
  • Predicted function