Enterococcus faecium are associated with the modification of gut microbiota and shrimp post-larvae survival
Animal Microbiome volume 3, Article number: 88 (2021)
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 , has subsequently dominated the field of aquatic invertebrate gut microbiomes . With increasingly threatened by the onset of various diseases in shrimp aquaculture , 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 , probiotics have become promising candidates for maintaining host health for its “green” and long-term research . 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 . 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 . The probiotic consortium has demonstrated an effective way to reduce the abundance of pathogenic Vibrio species and to prevent the mortality during Vibrio challenges . 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 .
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.  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 . 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 . 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 . 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 . 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).
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 . 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  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 . 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).
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).
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).
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.
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).
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).
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 .
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 . 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 . 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 . 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 , 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 . 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 . 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 . 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 . 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 . 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 . 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 . 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 . 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 . 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 . Using the UCHIME algorithm, chimeric sequences were identified  and removed from subsequent analysis. Bacterial phylotypes were identified using UCLUST  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 (https://www.arb-silva.de/documentation/release-128/). The matching results yielded 2,702,120 reads, and 26,002 reads were unclassified due to chimeras. Archaea and Chloroplast sequences were removed using “filter_taxa_from_otu_table.py”, 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.
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 . 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 . We expected bacterial population with wider niche breadth at the community level to be more metabolically flexible , and the analysis was conducted using the “niche.width” in spaa R package . Functional profiling of bacterial taxa was carried out using the Tax4Fun2 R package . 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 , 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) . 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.
The python-pipeline v.2.7 (available from http://oligotyping.org) 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 . 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.
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 , 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
Feed conversion ratio
Specific growth rate
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
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Xiong J. Progress in the gut microbiota in exploring shrimp disease pathogenesis and incidence. Appl Microbiol Biot. 2018;102:7343–50.
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.
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.
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.
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.
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.
Faust K, Raes J. Microbial interactions: from networks to models. Nat Rev Microbiol. 2010;10:538–50.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Claesson M, Clooney A, O’toole P. A clinician’s guide to microbiome analysis. Nat Rev Gastroenterol Hepatol. 2017;14(10):585–95.
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.
Shade A. Diversity is the question, not the answer. ISME J. 2017;11:1–6.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Chapman CMC, Gibson GR, Rowland I. Health benefits of probiotics: are mixtures more effective than single strains? Eur J Nutr. 2011;50:1–17.
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.
Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27:2194–200.
Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.
Anderson M. Distance-based tests for homogeneity of multivariate dispersions. Biometrics. 2006;62(1):245–53.
Anderson M. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001;26(1):32–46.
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.
Zhang J. Spaa: Species association analysis. R package version 0.2.1, http://CRAN.R-project.org/package=spaa; 2013.
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.
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.
Kursa MB, Jankowski A, Rudnicki WR. Boruta-A system for feature selection. Fund Inform. 2010;101:271–86.
Erdös P, Réyni A. On the evolution of random graphs. Publ Math Inst Hung Acad Sci. 1960;5:17–60.
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.
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
The authors have no conflict of interest to declare.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
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). https://doi.org/10.1186/s42523-021-00152-x
- Shrimp post-larvae
- Enterococcus faecium
- Gut microbiota
- Predicted function