Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Assessment of dietary supplementation with galactomannan oligosaccharides and phytogenics on gut microbiota of European sea bass (Dicentrarchus Labrax) fed low fishmeal and fish oil based diet

Abstract

There is an increasing interest from the aquafeed industry in functional feeds containing selected additives that improve fish growth performance and health status. Functional feed additives include probiotics, prebiotics, organic acids, and phytogenics (substances derived from plants and their extracts). This study evaluated the effects of dietary inclusion of a mucilage extract rich in galactomannan oligosaccharides (GMOS), a mixture of garlic and labiatae-plants oils (PHYTO), and a combination of them (GMOSPHYTO), on gut microbiota composition of European sea bass (Dicentrarchus labrax) fed with a low fishmeal (FM) and fish oil (FO) diet. Three experimental diets and a control diet (plant-based formulation with 10% FM and 6% FO) were tested in a 63-days feeding trial. To analyze the microbiota associated to feeds and the intestinal autochthonous (mucosa-adhered) and allochthonous (transient) microbial communities, the Illumina MiSeq platform for sequencing of 16S rRNA gene and QIIME2 pipeline were used. Metabarcoding analysis of feed-associated bacteria showed that the microbial communities of control (CTRL) feed deeply differed from those of experimental diets. The number of reads was significantly lower in CTRL feed than in other feeds. The OTU (operational taxonomic unit) number was instead similar between the feeds, ranging from 42 to 50 OTUs. The variation of resident gut microbiota induced by diet was lower than the variation of transient intestinal microbiota, because feedstuffs are a major source of allochthonous bacteria, which can temporarily integrate into the gut transient microbiome. However, the composition of transient bacterial communities was not simply a mirror of feed-borne bacteria. Indeed, the microbial profile of feeds was different from both faecal and mucosa profiles. Our findings suggest that the dietary inclusion of GMOS (0.5%) and PHYTO (0.02%) in a low FM and FO diet induces changes in gut microbiota composition of European sea bass. However, if on allochthonous microbiota the combined inclusion of GMOS and PHYTO showed an antagonistic effect on bactericidal activity against Vibrionales, at mucosa level, only GMOSPHYTO diet increased the relative abundance of Bacteroidales, Lactobacillales, and Clostridiales resident bacterial orders. The main beneficial effects of GMOS and PHYTO on gut microbiota are the reduction of coliforms and Vibrionales bacteria, which include several potentially pathogenic species for fish, and the enrichment of gut microbiota composition with butyrate producer taxa. Therefore, these functional ingredients have a great potential to be used as health-promoting agents in the farming of European sea bass and other marine fish.

Introduction

Farming of carnivorous fish species still relies on fishmeal (FM) and fish oil (FO) that represent the optimal sources of protein and lipids in aquafeeds. However, this dependence has to be overcome because FM and FO are no longer cost-effective and environmentally sustainable resources. Therefore, in the last decades, fish farmers and commercial feed producers have made substantial efforts to progressively reduce the proportion of these products in aquaculture feed, by replacing ground-up forage fish with terrestrial plant-derived protein and lipid sources [1, 2].

However, only a partial replacement of FM and FO is possible without negatively affecting fish growth performance and health status. The main drawbacks of using high levels of plant proteins in carnivorous fish feed are related to several nutritional imbalances and to the presence of a wide variety of anti-nutritional factors that challenge fish intestine functionality and health, thus reducing nutrient absorption and animal growth [3, 4]. Indeed, various studies have described a range of inflammatory gut reactions in carnivorous fish species, such as Atlantic salmon (Salmo salar), and European sea bass (Dicentrarchus labrax) in response to high replacement levels of marine-derived raw materials with plant-derived raw materials in the diet [4, 5, 6]. The shortening of the primary and secondary intestinal mucosal folds, higher level of infiltrated leucocytes in the lamina propria and submucosa, damaged microvilli, altered gut microbiota profiles, and increased gut permeability are the most typical signs of such inflammatory response [5, 6, 7]. These negative-side effects reduce the capacity of enterocytes to absorb nutrients and can promote the translocation of indigenous or opportunistic bacteria, potentially pathogenic for the host [8, 7]. In this context, there is an increasing interest in functional feeds containing selected additives that can help to prevent or mitigate the intestinal epithelium damage caused by extreme plant-based diet formulations [913]. Feed additives have quite diverse chemical natures and characteristics; they may be both, nutritive and non-nutritive ingredients and may interact with fish physiology by either direct or indirect mechanisms [14]. Functional feed additives include probiotics, prebiotics, immune-stimulants, organic acids, nucleotides, exogenous enzymes, and phytogenics (substances derived from herbs, spices, other plants and their extracts). Their application into diet formulations targets a specific purpose: probiotics, prebiotics, phytogenics, and immune-stimulants target the improvement of intestinal health, stress, and disease resistance, whereas acidifiers, exogenous enzymes and indirectly probiotics are used to improve fish performance by enhancing feed digestibility, or counteracting the negative effects of antinutrients [14].

Among prebiotics, galactomannans from different sources, have been related to enhanced fish growth performance and health status [1520]. Galactomannans are heteropolysaccharides structurally composed of D-mannose, which makes up the backbone, and D-galactose that forms single branches along the mannan chain. Galactomannans originate from two sources; the first and main source is represented by plants, in particular, the endosperm of Leguminosae family's dicotyledonous plants seeds; the second is represented by microorganisms, such as yeast and fungi, in which galactomannans constitute an essential component of the cell walls.

Galactomannans from plants are in the form of soluble dietary fibers, thereby acting as prebiotics. The basic criteria to consider an ingredient as a prebiotic have been firstly established by Gibson and Roberfroid [21]. They defined a prebiotic as a nondigestible food ingredient that beneficially affects the host by selectively stimulating the growth and/or activity of one or a limited number of bacterial species already resident in the intestine.

Unlike yeasts’ mannans that have been widely investigated in several fish species leading to contrasting results on their effects [2224, 9, 2529], galactomannans of plant origin are still scarcely explored as functional ingredients in aquafeeds. Nevertheless, it is known that the improvement of growth performance, health, and disease resistance in fish fed with prebiotics are frequently connected with changes in their gut microbial communities [18, 3032].

In addition to changes in gut microbiota, other indirect mechanisms of action of prebiotics include bacterial end-products of fermentation. By definition, prebiotics are not digested by the host itself, but are fermented by bacteria present in the host’s gut. The end products of prebiotic fermentation are short chain fatty acids, predominantly acetic, propionic and butyric acids. Among them butyrate is considered to be a preferred energy source for colonocytes and in vertebrates, including fish, plays a pivotal role in the maintenance of overall gut health, intestinal morphology, and function [3334, 11, 31].

By comparison with prebiotics, only limited information is available on the potential benefits of phytogenics on fish health and, in general, on their application in animal nutrition [35, 13]. These products have great potential to be used in aquaculture due to a number of beneficial biological activities that include growth promotion, appetite stimulation, modulation of immune and antioxidant response in addition to antiparasitic, antibacterial, anaesthetic, and antistress activities [3643]. These activities are positively reflected on feed palatability, fish digestive functions, and intestinal microbiota structure.

Phytogenic feed additives using plant extracts show greater modes of action in animal nutrition compared to synthetic, nature-identical substances [13, 38]. This advantage is based on the synergistic effects of all agents within a plant, which have not been reduced to the effects of a single lead substance. The challenge is finding the right combinations of natural substances while fully exploiting the synergy among their active ingredients.

Essential oils (EOs) represent the major group of phytogenic feed additives. EOs are a natural mixture of various organic substances synthesized by aromatic plants during secondary metabolism. EOs usually contain high amounts of different compounds, such as terpenes, alcohols, acetones, phenols, acids, aldehydes, and esters [44]. EOs are potent antimicrobial agents against different strains of pathogenic and food borne bacteria [45, 43]. Therefore, due to their antibacterial activity, phytochemicals can exert a prebiotic-like effect modulating fish gut bacterial composition [46, 12].

Accordingly, based on the aforementioned properties of galactomannans and phytogenics, the aim of this study was to investigate the effects of dietary inclusion of galactomannan oligosaccharides (GMOS), a mixture of garlic and labiatae-plants oils (PHYTO), and a combination of them (GMOSPHYTO), on gut microbiota composition of European sea bass fed with a plant-based diet. A 16S metabarcoding approach” was used to characterize gut bacterial communities of European sea bass.

Material and methods

Ethics statement

All procedures involving fish complied with the guidelines of the European Union Council (86/609/EU) and Spanish legislation (RD 53/2013) and were approved by Bioethical Committee of the University of Las Palmas de Gran Canaria (Ref. 007/2012 CEBA ULPG).

Feeding trial, diets, and sampling

Details of the feeding trial, set at Parque Científico-Tecnológico Marino (PCTM) of the University of Las Palmas de Gran Canaria (Telde, Canary Island, Spain), have been described by Torrecillas et al. [20]. Briefly, four isoenergetic and isonitrogenous experimental diets were manufactured by an extrusion process in the BioMar Tech-Centre (Brande, Denmark) starting from a plant-based formulation (10% FM, 6% FO) and containing different additives: 0.5% of galactomannan oligosaccharides from mucilage (GMOS; Delacon, Austria), 0.02% of a mixture of garlic and labiate-plants oils (PHYTO; Delacon, Austria) and a combination of both additives (GMOSPHYTO, 0.52%). Principal ingredients, and proximate composition of diets is reported in Table 1, modified from Torrecillas et al. [20]. GMOS was included in the diet in the mix pre-extrusion process, PHYTO was included in post extrusion process by vacuum coating and homogenized with the dietary fish oil. The GMOS and PHYTO dosages were included in the diets accordingly to commercial recommendations. For GMOS, these recommendations were based on internal trials and other previous studies with similar products for the same fish species. For the blend of herbal extracts, the dietary concentration tested was also based on commercial recommendations, which mainly depended on MIC internal assays against Vibrio sp. The stability of the supplements was evaluated before and after feed production, as well as at the beginning of the feeding trial.

thumbnail
Table 1. Main ingredients and proximate composition of the diets (modified from Torrecillas et al., 2019).

https://doi.org/10.1371/journal.pone.0231494.t001

Nine hundred European sea bass juveniles were randomly distributed in 12 fiberglass tanks of 500 L in an open water system and fed ad libitum with four different diets in triplicate for 63 days. At the end of the trial, two fish per replicate (6 fish/diet) were sampled and the whole intestine was aseptically removed. The animals used for sampling were sacrificed by an overdose of anaesthetic (clove oil) using water bath immersion, and all efforts were made to minimize pain, stress, and discomfort in the animals. The faecal matter containing allochthonous (transient) intestinal bacteria was removed from each intestine by squeezing and collected in a sterile tube with 800 μl of Xpedition Lysis/Stabilization Solution (Zymo Research). The autochthonous (adhered) microbiota was obtained by scraping the whole intestinal mucosa (excluding pyloric caeca) with a sterile cotton swab. The tip of swap was immediately immersed in 200 μl of Xpedition Lysis/Stabilization Solution and vortexed to facilitate bacteria releasing as described in Rimoldi et al. [47]. Both faecal and mucosa samples were stored at room temperature for up to 48 hrs until bacterial DNA extraction.

The remained fish following this research were housed at the Las Palmas University facilities to be used in practical lessons for undergraduate students of veterinary and marine sciences faculties, focused to the learning of aquaculture husbandry practices. Fish were not used in further research, due to the influence of experimental diets in the intestine of different groups, giving an erroneous starting point for other feeding trials. Besides, animal’s final weight was high for new trials.

DNA extraction

The bacterial DNA was extracted from 4 samples from each feed, 6 samples of faeces and 6 of intestinal mucosa per each dietary fish group. For each extraction, 200 mg of feed, 200 mg of faeces, and 200 microliters of mucosal bacteria suspension were used. The extraction of bacterial DNA from the feeds was done in parallel to biological samples, right after the feeding trail.

DNeasy PowerSoil Kit (Qiagen, Italy) was used for extraction, following the manufacturer’s instructions with few modifications at the lysis step. Specifically, lysis was performed in PowerBead Tubes by means of a TissueLyser II (Qiagen, Italy) set at 25 Hz for 2 min. As a negative control for the extraction procedure, a sample with only lysis buffer was processed in parallel with all samples. The concentration of extracted DNA was assessed using a NanoDropTM 2000 Spectrophotometer (Thermo Scientific, Italy). DNA was stored at– 20°C until the PCR reaction was performed.

16S Illumina library construction and high-throughput sequencing

Methodology applied for 16S rRNA gene library preparation and sequencing have been described in Rimoldi et al. [47] and Terova et al. [48]. The Illumina protocol “16S Metagenomic Sequencing Library Preparation for Illumina MiSeq System” (#15044223 rev. B) was applied for library preparation. The V3-V4 region was amplified from 50 ng of microbial genomic DNA using Platinum® Taq DNA Polymerase High Fidelity kit (Thermo Fisher Scientific, Italy) and tailed using forward and reverse primers Pro341F (5′-CCTACGGGNBGCASCAG -3′) and Pro805R (5′-GACTACNVGGGTATCTAATCC -3′) selected by Takahashi et al. [49]. The amplicon length was checked on Agilent 2100 Bioanalyzer trace and the expected size was ~550 bp.

Illumina paired-end adapters with unique Nextera XT indexes were ligated to 16S amplicons using Nextera XT Index Kit (Illumina, San Diego, CA). Libraries were purified and normalized using the SequalPrep ™ Normalization Plate kit (Thermo Fisher), pooled at equimolar concentrations and diluted at 6 pM prior to sequencing on the Illumina MiSeq platform (Illumina). Libraries were sequenced with v3 chemistry on 300PE MiSeq runs.

16S amplicon sequencing data analysis

Raw FASTQ sequencing data were processed using the open-source bioinformatics pipeline QIIME 2 (v. 2018.4) at the default setting [50]. To reconstruct the original amplicons, the reads generated by three different runs (one with all the samples, and two extra runs with mucosal and faecal samples that gave suboptimal Good’s coverage value (< 99.5%) and low number of reads in the first run) were separately processed with DADA2 [51]. R1 and R2 paired reads were trimmed at both 3’ and 5’ ends using Cutadapt v.2018.4.0 software, filtered for base quality (Q>30) and merged. The remaining high-quality reads were de-replicated to obtain unique sequence (uniques) and chimeric sequences removed using QIIME DADA2 denoise-paired command. Denoised sequences with 99% or higher identity against uniques were de novo clustered into Operational Taxonomic Units (OTUs). Only the OTUs that represented at least 0.005% of total reads were kept. Taxonomy was assigned, down to genus level, using FeatureData [Sequence] artefact against the reference database GreenGenes v.13-8. Reads of mitochondrial or eukaryotic origin were excluded.

Alpha and beta diversity statistics were performed by alpha-phylogenetic and beta-phylogenetic command, respectively. To evaluate the compositional alpha- and beta-diversity, the samples were divided into two macro-groups (faeces+feeds, and mucosae+feeds) and separately analysed. Alpha diversity was calculated based on rarefied OTU table (rarefied at the lowest sample size) using observed OTUs, Shannon, Pielou’s evenness, and Faith PD. Difference in taxonomic profiles among samples (beta diversity) was calculated using weighted (presence/absence/abundance matrix) and unweighted (presence/absence matrix) UniFrac distances [5253]. The dissimilarity matrices were graphically represented by three-dimensional PCoA plots.

Statistics

All data were verified for normality and homoscedasticity of variance by Shapiro-Wilk’s and Levene’s test, respectively. Depending on whether the assumptions were satisfied or not, differences between groups were analysed by one-way ANOVA with Tukey’s pairwise test or by nonparametric Kruskal-Wallis followed by Dunn’s post hoc test. Statistical significance was set at p<0.05. All analyses were performed using Past3 software [54].

The number of reads across samples was normalized by sample size and the relative abundance (%) of each taxon was calculated. OTUs assigned to the phylum Cyanobacteria (class Chloroplast) and to mitochondria, were considered potential plant contaminants and were removed from the analysis. For each dietary group, only taxa present in at least half samples and with a mean overall abundance of more than 1% (up to order level) and 0.5% at family and genus level were considered for statistical analysis. Less abundant taxa were indicated as “Others” in the histograms of microbial community profiles of fish. They were not considered for statistical analysis, since we decided to analyse the effects of diet only on the most representative OTUs.

Before being statistically analysed, the resulting microbial profiles relative abundances of each OTU were calculated as the angular transformation (arcsine of the square root). The angular transformation of data is commonly used to make normal a binomial distribution like in our case in which the data were expressed as decimal fractions and percentages.

Multivariate analysis of beta diversity was tested using non-parametric analysis of similarities (ANOSIM) and Adonis tests with 999 permutations.

Results

Fish growth parameters and biometry

After 63 days of feeding, the experimental diets did not induce differences in fish growth or diet utilization. All growth performance data concerning the present feeding trial have been recently reported by Torrecillas et al. [20].

Illumina MiSeq sequencing efficiency

The obtained raw sequencing data were divided into two macro-groups (i.e. faeces+feeds and mucosa+feeds) that were analysed separately, i.e each macro-group indicates the set of samples statistically analysed together by the QIIME2 pipeline. Feeds were sequenced only one time, but they were statistically analysed twice, once per each macro-group.

We obtained a total number of quality-filtered reads of 1,338,727, which corresponded to 34,326 ± 17,538 reads per sample, for faeces+feeds macro-group and 1,114,122 reads, corresponding to 29,319 ± 18,581 reads per sample, for mucosa+feeds macro-group.

An overall of 125 and 210 unique sequences were found in faeces+feeds and mucosa+feeds macro-groups, respectively (S1 Data File). Only two samples (one of GMOS feed and one of mucosa from GMOSPHYTO group) were excluded after rarefaction analysis, since their rarefaction curve did not reach the plateau. On the basis of rarefaction analysis results, faeces+feeds and mucosae+feeds samples were rarefied at 4500 and 5000 reads, respectively. Rarefaction curves of Observed OTUs numbers of two macro-groups of samples have been reported in S1 and S2 Figs. All fastq sequencing files were deposited in the European Nucleotide Archive (EBI ENA) public database under the accession project code: PRJEB32279.

Microbial profile of experimental feeds

Metabarcoding analysis of feed-associated bacteria showed that the microbial communities of CTRL feed deeply differed from those of experimental functional diets GMOSPHYTO, PHYTO, and GMOS. The number of reads taxonomically classified according to Greengenes database was 6,315 ± 760; 30,613 ± 5,521; 42,689 ± 8,558; and 30,953 ± 11,762 for CTRL, GMOS, GMOSPHYTO, and PHYTO diets, respectively. This number was significantly lower in CTRL feed (p<0.05) in comparison to other feed samples. In addition, Pielou’s evenness and Shannon alpha diversity indices of CTRL feed differed from those of functional feeds, being both indices significantly higher in CTRL samples. The OTU number was instead comparable among the diets, ranging from 42 to 50 OTUs (S1 Table).

After removing the OTUs assigned to eukaryotic sequences, the most abundant bacterial taxa were mainly comprised of 4 phyla, 7 classes, 10 orders, 14 families, and 11 genera. We have presented the profiles of microbial communities for each feed at phylum (Fig 1A) and genus (Fig 1B) taxonomic level and reported the relative abundance (%) of the most abundant taxa found in feed samples (S2 Table).

thumbnail
Fig 1.

Relative abundance (%) of the most prevalent bacteria in CTRL, GMOS, PHYTO, and GMOSPHYTO feeds at phylum (A), and genus (B) taxonomic level. Only bacteria with an overall abundance of ≥1% were reported. Bacteria with lower abundance were pooled and indicated as “Others”.

https://doi.org/10.1371/journal.pone.0231494.g001

At phylum level, functional feeds GMOS, GMOSPHYTO, and PHYTO were characterized by higher percentage of Firmicutes (55–61%) than CTRL feed (49%). Conversely, microbiota associated to CTRL feed showed a higher relative abundance of Proteobacteria (40%), mainly represented by the Gammaproteobacteria class, as compared to other feeds (29–34%) (Fig 1A, S2 Table). Accordingly, a high amount of the Shewanellaceae (9.9%), Enterobacteriaceae (4.6%), and Vibrionaceae (11%) families was found in the CTRL diet. Feeds containing prebiotic and/or phytogenic were instead enriched in Lactobacillaceae (66–73%).

At genus level, CTRL feed had higher relative abundance of Corynebacterium (3.8%), Enterococcus (5.5%), Cetobacterium (5.24%), Shewanella (9.9), Erwinia (4.6%), and Photobacterium (11.2%), whereas the Lactobacillus genus was dominant, with percentages ranging between 66 and 73%, in GMOS, GMOSPHYTO, and PHYTO diets (Fig 1B, S2 Table).

The permutational multivariate analysis made with Adonis and ANOSIM test on both, unweighted and weighted UniFrac distance data, statistically confirmed the differences between feed-associated microbial communities (beta-diversity). Adonis and ANOSIM tests revealed a significant difference (p<0.05) between CTRL and functional feeds (GMOS, GMOSPHYTO, and PHYTO) in both, type (R2>0.51, R>0.50) and abundance (R2>0.79, R = 1) of taxa found. The results of multivariate analysis are summarized in Table 2.

thumbnail
Table 2. Results of permutational multivariate analysis of variance (Adonis) and analysis of similarity (ANOSIM) based on unweighted and weighted UniFrac distance matrices, using abundance data of feed-associated bacterial communities.

https://doi.org/10.1371/journal.pone.0231494.t002

Microbial profile and dietary modulation of allochthonous gut communities

The whole microbial community profile of 24 faecal samples was mainly composed of 11 phyla, 15 classes, 27 orders, 52 families, and 59 genera. By considering only the most representative taxa, the overall allochthonous microbiota consisted of 6 phyla (Actinobacteria, Bacteriodetes, Firmicutes, Proteobacteria, Spirochaetes, and Tenericutes), 9 classes (Actinobacteria, Bacteroidia, Bacilli, Clostridia, Alphaproteobacteria, Betaproteobacteria, Gammaproteobacteria, [Brevinematae], and Mollicutes), 17 orders (Actinomycetales, Bacteroidales, Bacillales, Lactobacillales, Clostridiales, Rhizobiales, Rickettsiales, Burkholderiales, Neisseriales, Alteromonadales, Enterobacteriales, Oceanospirillales, Pasteurellales, Pseudomonadales, Vibrionales, [Brevimatales], and Mycoplasmatales), 23 families (Corynebacteriaceae, Propionibacteriaceae, Bacteroidaceae, Bacillaceae, Staphylococcaceae, Lactobacillaceae, Streptococcaceae, Lachnospiraceae, Ruminococcaceae, Aurantimonadaceae, Methylobacteriaceae, Rhodobacteraceae, Comamonadaceae, Oxalobacteraceae, Neisseriaceae, Shewanellaceae, Enterobacteriaceae, Halomonadaceae, Pasteurellaceae, Moraxellaceae, Vibrionaceae, Brevinemataceae, and Mycoplasmataceae), and 24 genera (Corynebacterium, Propionibacterium, Bacteroides, Staphylococcus, Lactobacillus, Streptococcus, Faecalibacterium, Ruminococcus, Methylobacterium, Paracoccus, Aquabacterium, Janthinobacterium, Glaciecola, Shewanella, Escherichia, Cobetia, Aggregatibacter, Acinetobacter, Enhydrobacter, Psychrobacter, Enterovibrio, Photobacterium, Vibrio, and Mycoplasma).

The profiles of faecal microbial communities for each feeding group and individual fish are presented at phylum and family taxonomic level in Figs 2 and 3, respectively. To elaborate alpha rarefaction analysis (alpha diversity), samples were normalized at a sequencing depth of 4,500 reads. Administration of combined functional diet GMOSPHYTO, but not GMOS and PHYTO diets, significantly decreased (p<0.05) the number of observed OTUs with respect to CTRL diet. Conversely, diet type did not affect either phylogenetic diversity (Faith PD) or entropy (Shannon and Pielou’s evenness) (Table 3).

thumbnail
Fig 2.

Relative abundance (%) of the most prevalent allochthonous bacterial phyla in each dietary group (A) and in individual fish (B). In the figure, all taxa with an overall abundance of ≥1% were reported. * indicates outlier samples excluded from relative abundance analysis.

https://doi.org/10.1371/journal.pone.0231494.g002

thumbnail
Fig 3.

Relative abundance (%) of the most prevalent allochthonous bacterial families in each dietary group (A) and in individual fish (B). In the figure, all taxa with an overall abundance of ≥1% were reported.

https://doi.org/10.1371/journal.pone.0231494.g003

thumbnail
Table 3. Original number of reads per group-treatment assigned to OTUs, and alpha diversity metrics values (rarefied at 4500 reads) of faecal microbial community in sea bass fed CTRL, GMOS, GMOSPHYTO, and PHYTO diets.

https://doi.org/10.1371/journal.pone.0231494.t003

Analysis of beta-diversity revealed an overall effect of diet on microbial communities both in presence/absence (unweighted UniFrac), and in relative abundance (weighted UniFrac) of OTUs (Fig 4A and 4B). However, the major effect of diet was observed in terms of relative abundance of taxa. The first principal coordinate PC1 of weighted UniFrac PCoA plot explained, indeed, 77% of the variation between individuals (Fig 4B). Interestingly, CTRL and GMOS fish groups clustered together and distinctly from PHYTO and GMOSPHYTO samples, which in turn grouped together (Fig 4B). Additionally, both in unweighted (Fig 4A) and weighted (Fig 4B) UniFrac PCoA, faecal samples appeared clearly separated from feed samples, thus indicating that observed differences between transient intestinal bacterial communities were not simply a consequence of undigested feed that might have been present in the intestinal lumen.

thumbnail
Fig 4. Beta diversity metrics.

Principal coordinate analysis (PCoA) of Unweighted (A) and Weighted (B) Unifrac distances of gut allochthonous microbial communities associated to different diet. The figures show the 3D plot of individual fish according to their microbial profile at genus level.

https://doi.org/10.1371/journal.pone.0231494.g004

The permutational multivariate analysis applying Adonis and ANOSIM tests on UniFrac distance data, fully confirmed the PCoA results. Multivariate analysis on weighted data revealed significant differences between CTRL and GMOS versus PHYTO and GMOSPHYTO dietary groups (R2 > 0.41, R > 0.38, p<0.05). Results of pairwise comparisons on phylogenetic distances are summarized in Table 4.

thumbnail
Table 4. Results of permutational multivariate analysis of variance (Adonis) and analysis of similarity (ANOSIM) based on unweighted and weighted Unifrac distance matrices, using abundance data of faecal bacterial communities at genus level.

https://doi.org/10.1371/journal.pone.0231494.t004

In line with beta-diversity analysis, the major differences in terms of transient microbiota composition and taxa abundance were found between CTRL and fish feed diets containing phytogenics (i.e. PHYTO and GMOSPHYTO).

At phylum level, there were no relevant changes in transient gut microbiota profiles, which were mainly constituted by Proteobacteria, Firmicutes, and Actinobacteria phyla (Fig 2). A significant decrease (p<0.05) of bacteria belonging to the Alpha- and Betaproteobacteria classes was found in fish fed with PHYTO and GMOSPHYTO diets in comparison to control and GMOS groups that were characterized by a high number of bacteria assigned to the Rickettsiales and Neisseriales orders (Table 5). PHYTO and GMOSPHYTO dietary groups were enriched in Gammaproteobacteria essentially represented by the Vibrionales order (85–88%). Other bacterial orders belonging to Gammaproteobacteria, such as Alteromonadales, Enterobacteriales, and Pseudomonadales were practically detected only in fish fed with CTRL and GMOS diets, whereas Oceanospirillales and Pasteurellales were exclusively found in GMOS faecal samples. PHYTO were also enriched in the Clostridia class.

thumbnail
Table 5. Mean relative abundance (%) ± SD (n = 6) of the most prevalent phyla, orders, classes, families, and genera found in faecal samples of sea bass fed with four experimental diets.

https://doi.org/10.1371/journal.pone.0231494.t005

At family level changes in taxa abundance practically reflected what was observed at higher taxonomic levels. Interestingly, Propionibacteriaceae, represented by the Propionibacterium genus, were negatively affected by phytogenics containing diets. Similarly, the number of lactic acid bacteria assigned to the Streptococcaceae family was significantly reduced by all functional diets (Table 5).

When the analysis was performed at genus taxonomic level, only few genera resulted significantly influenced by diet. Among them Streptococcus and Janthinobacterium were abundant in control samples, whereas Propionibacterium and Methylobacterium in CTRL and GMOS samples. The Bacteroides, Faecalibacterium and Ruminococcus genera were present only in PHYTO sample. Lastly, PHYTO and GMOSPHYTO diets seemed to have a bactericidal activity against the genus Escherichia.

Microbial profile and dietary modulation of autochtonous gut communities

After removing the reads corresponding to eukaryotic sequences, the entire microbial community profile of 23 intestinal mucosal samples consisted of 10 phyla, 21 classes, 40 orders, 75 families, and 112 genera. However, if we considered only the most representative taxa, the overall autochthonous microbiota was composed of 7 phyla (Actinobacteria, Bacteroidetes, Firmicutes, Fusobacteria, Proteobacteria, Spirochaetes, and Tenericutes), 11 classes (Actinobacteria, Bacteroidia, Flavobacteriia, Bacilli, Clostridia, Fusobacteriia, Alphaproteobacteria, Betaproteobacteria, Gammaproteobacteria, [Brevinematae], and Mollicutes), 21 orders (Actinomycetales, Bacteroidales, Flavobacteriales, Bacillales, Lactobacillales, Clostridiales, Fusobacteriales, RF32, Rhizobiales, Rhodobacterales, Sphingomonadales, Burkholderiales, Neisseriales, Alteromonadales, Enterobacteriales, Oceanospirillales, Pseudomonadales, Salinisphaerales, Vibrionales, [Brevinematales], and Mycoplasmatales), 42 families, and 51 genera (Tables 6 and 8). Profiles of gut mucosa communities for each feeding group and individual fish are graphically presented at phylum (Fig 5) and order (Fig 6) taxonomic level.

thumbnail
Fig 5.

Relative abundance (%) of the most prevalent autochthonous bacterial phyla in each dietary groups (A) and in individual fish (B). In the figure, all taxa with an overall abundance of ≥1% were reported.

https://doi.org/10.1371/journal.pone.0231494.g005

thumbnail
Fig 6.

Relative abundance (%) of the most prevalent autochthonous bacterial orders in each dietary groups (A) and in individual fish (B). In the figure, all taxa with an overall abundance of ≥1% were reported.

https://doi.org/10.1371/journal.pone.0231494.g006

thumbnail
Table 6. Original number of reads per group-treatment assigned to OTUs, and alpha diversity metrics values (rarefied at 5000 reads) of mucosa microbial community in sea bass fed CTRL, GMOS, GMOSPHYTO, and PHYTO diets.

https://doi.org/10.1371/journal.pone.0231494.t006

To elaborate alpha rarefaction analysis (alpha diversity), samples were rarefied at 5,000 reads. Albeit, PHYTO samples had the lowest number of reads, the autochthonous gut microbial communities of fish fed with this diet were characterized by higher biodiversity than control and GMOS dietary groups, as indicated by high value of Shannon diversity index. Contrariwise, functional diets did not affected species richness (Observed OTUs) (Tables 6 and 7).

thumbnail
Table 7. Results of permutational multivariate analysis of variance (Adonis), and analysis of similarity (ANOSIM) based on unweighted and weighted UniFrac distance matrices using abundance data of intestinal mucosa-associated bacterial communities.

https://doi.org/10.1371/journal.pone.0231494.t007

Both weighted and unweighted UniFrac analyses were performed to measure microbial community diversity (beta-diversity). Data of UniFrac matrices were projected onto three-dimensional plots using principal coordinates analysis (PCoA) (Fig 7). There was a diet effect on unweighted UniFrac distances, i.e. on the presence/absence of specific taxa. Indeed, unweighted UniFrac PCoA revealed a clustering of PHYTO samples, which grouped separately from CTRL and GMOSPHYTO fish along the first principal coordinate PC1 (47% of the variation) (Fig 7A). On the contrary, in weighted PCoA, most samples were broadly indistinguishable (Fig 7B).

thumbnail
Fig 7. Beta diversity metrics.

Principal coordinate analysis (PCoA) of Unweighted (A) and Weighted (B) Unifrac distances of gut autochthonous microbial communities associated to different diet. The figures show the 3D plot of individual fish according to their microbial profile at genus level.

https://doi.org/10.1371/journal.pone.0231494.g007

Likewise, for faecal bacterial communities, microbial profiles of feeds and mucosa differed to each other, both qualitatively and qualitatively. Multivariate analysis on UniFrac distance data confirmed PCoA results, indicating a significant divergence between feeding groups for only unweighted UniFrac distance matrix. Pairwise test Adonis and ANOSIM on the unweighted UniFrac data showed that fish fed diet PHYTO significantly differed (p<0.05) from CTRL (R2 = 0.19, R = 0.28) and GMOSPHYTO (R2 = 0.21, R = 0.30) (Table 6).

The variation of resident gut microbiota induced by diet was definitely lower than the variation of transient intestinal microbiota. The mucosa adhered microbial community of sea bass was mainly dominated, regardless of the diet, by three phyla: Proteobacteria, Firmicutes and Bacteroidetes (Fig 5, Table 8).

thumbnail
Table 8. Mean relative abundance (%) ± SD (n = 6) of the most prevalent phyla, orders, classes, families, and genera found in gut mucosa samples of sea bass fed with four experimental diets.

https://doi.org/10.1371/journal.pone.0231494.t008

Among them, the amount of Firmicutes was positively influenced (p<0.05) by administration of GMOSPHYTO diet (Table 8). This was essentially due to the enrichment in bacteria belonging to the Clostridiales order of the Clostridia class. No differences in relative abundance were found at lower taxonomic level, however several families and genera seemed to have a diet-specific association. The Prevotellaceae, Ruminococcaceae, Bacteroidaceae, and Veillonellaceae families were identified only in fish fed diet GMOSPHYTO, whereas Fusobacteriaceae were found solely in intestine of control feeding group (Table 8).

Accordingly, the orders Clostridiales, Lactobacillales and Pseudomonadales, as well as the Bacteroides and Prevotella genera were particularly abundant in resident gut microbiota of fish fed with GMOSPHYTO diet. Similarly, in the same fish group, several genera of Clostridiales, such as Blautia, Lachnospira, Roseburia, Faecalibacterium, Dialister, and Phascolarctobacterium, were identified. On the contrary, Cetobacterium, belonging to the Fusobacteriaceae family, was found only in the autochthonous gut microbiota of control fish (Table 8).

Discussion

Improving fish health is a major concern for fish farmers. Since 2006 when a European Union-wide ban on the use of antibiotics as growth promoters in animal feed entered into effect (EU Regulation No. 1831/2003), functional feeds with selected additives such as probiotics, prebiotics and phytogenics started to be used in aquaculture and other animal production industries to promote fish health [55].

In our study, the choice to use functional feed additives, such as galactomannan oligosaccharides (GMOS) from mucilage and phytogenics in low FM/FO diets, meets the consumer’s demand of eco-friendly production practices. European sea bass (Dicentrarchus labrax) promptly accepted all the experimental diets for all the duration of the feeding trial, but growth performance, feed intake and feed efficiency were unaffected by prebiotic and/or phytogenic incorporation [20].

Marker gene analysis of 16S rRNA was successfully applied to the study of both gut mucosa-associated (autochthonous) and lumen content (allochthonous) microbial communities. The analysis revealed that autochthonous microbiota was characterized by higher species richness than allochthonous microbiota. Number of bacterial genera observed in mucosa samples was indeed, higher than in faecal samples, however autochthonous microbiota was rarefied to a higher read count. Consistent with our results, similar differences in bacterial abundance between transient and mucosa-associated gut microbiota have been previously described in sea bass using 16S rRNA gene sequencing [56], as well as in other species [5759, 48].

As reported in literature, Proteobacteria resulted the most common bacteria phylum in both transient and resident microbial communities of sea bass intestine irrespective to the diet [56, 60]. No significant differences on allochthonous bacterial alpha-diversity and species richness were observed in response to diets containing functional ingredients. In contrast, the mucosa-associated microbiota of fish fed with PHYTO diet, supplemented with a mixture of garlic and labiate-plants oils, showed an increase of Shannon diversity index in comparison to the control feeding group. Carda-Diéguez and colleagues [56] did not find any differences in terms of autochthonous microbiota biodiversity in sea bass fed a functional diet containing β-glucans and essential oils, whereas in rainbow trout (Oncorhynchus mykiss), intestinal resident bacterial diversity and richness were negatively affected by garlic extract in the diet [61]. Nevertheless, in general, the lack of a negative effect on microbial diversity and species richness, like in our study, are considered desired features. High intestinal microbial biodiversity usually reflects a healthy status of the host, because a reduction in commensal bacterial diversity and/or species richness may result in diminished colonization resistance against incoming opportunistic pathogens, which have the potential to induce infection of the intestinal tract [6263].

Although transient and resident gut microbiota of sea bass was dominated by the same phyla (Proteobacteria, Firmicutes, Actinobacteria) irrespective to the diet, the present study clearly showed that gut microbiota can be modulated by dietary supplementation with functional ingredients.

Prebiotics, such as mannan oligosaccharides (MOS) and phytogenics have proved to be effective at modulating fish microbiota in several studies [64, 18, 46, 56, 65, 12, 61]. In the present study, allochthonous intestinal microbiota of sea bass fed GMOS diet, but not GMOSPHYTO, showed reduced number of Gammaproteobacteria, mainly represented by the Vibrionales order. This could be considered a positive effect of dietary GMOS administration, since the Vibrionales order includes several potentially pathogenic species for fish, such as Vibrio anguillarum, and Photobacterium damselae. Similarly, Guardiola et al. [66] showed in sea bream that dietary administration of fenugreek (Trigonella feonum-graecum L.) seeds, combined with probiotics strains, for 3 weeks, increased skin mucus bactericidal activity against the pathogenic bacterium P. damselae in comparison to fish of the control and other experimental groups. Seeds of fenugreek are a rich source of a wide variety of nutritive and bioactive compounds, including galactomannan, having several potent activities including antimicrobial and anti-inflammatory effects.

Accordingly, our tested diets with low levels of ingredients of marine origin, and supplemented with GMOS and PHYTO improved significantly European sea bass resistance to V. anguillarum after intestinal infection and stress challenge [20]. However, although the three functional diets reduced gut pathogen translocation rate, there was an evident greater protective effect when GMOS or PHYTO were supplemented to the diet individually. Indeed, the relative percentage of survival was 33%, 47%, 20%, and 40% for fish fed GMOS, PHYTO, GMOSPHYTO, and CTRL diets, respectively [20]. Furthermore, at morphology level, the combination of two additives, GMOS and PHYTO, invalidated their site-specific anti-inflammatory properties on fish intestine [20]. This observed effectiveness reduction could indicate a possible antagonistic effect between products, effect that appeared to be confirmed by our 16S analysis of gut microbial communities, which showed an increased amount of Vibrionales when GMOS and PHYTO were added together to the feed, but not when only GMOS was supplemented to the diet.

Similarly, previous evidences in European sea bass revealed a positive effect of prebiotic mannan-oligosaccharides (MOS) in terms of fish survival to V. anguillarum and V. alginolyticus infection [22, 67]. Also, juvenile trout fed with a diet supplemented with 0.2% MOS showed a significantly decrease of Aeromonas/Vibrio spp. in comparison to the control group [64]. MOS is able to bind the threadlike fimbriae of gram-negative pathogenic bacteria (i.e. Aeromonas, Vibrio, and Pseudomonas) preventing them from attaching to the gut wall thus avoiding their intestinal colonization [68]. Dietary inclusion of MOS caused also an enhancement in the number of cells secreting acid mucins in posterior gut of fish [25, 26]. Mucins are the major antiadhesive components of mucus and the improvement in mucus secretion represents, therefore, the primary defence mechanism against colonizing organisms of gut epithelial mucosal surface.

Bacteria assigned to the Clostridiales order were more abundant in intestinal content of fish fed PHYTO diet, as well as, in mucosa-associated microbiota of fish receiving GMOSPHYTO diet. In line with our results, a significant and positive correlation was found in resident intestinal microbiota of rainbow trout between the abundance of the genus Clostridium and the levels of garlic oil in the diet [61]. Previous studies have highlighted the importance of some members of the Clostridiales order, which includes several butyrate producers. It is known, indeed, that members of the genus Clostridium, besides being butyrogenic, can improve fish nutrition by providing essential fatty acids and vitamins [69, 70, 71]. Accordingly, among Clostridiales, the Faecalibacterium and Ruminococcus genera were found only in the intestinal lumen of fish fed PHYTO diet. In particular, the Ruminococcus genus plays an important role in the degradation of indigestible carbohydrate, such as resistant starch and dietary fibers, thus contributing to the more efficient energy utilization of feed and to intestinal health of host by producing butyrate as end product of dietary fiber fermentation [72, 73, 74]. Similarly, Faecalibacterium is the most important commensal butyrate-producing bacterium in the human colon and it is considered as a bio-indicator of human health [75].

Interestingly, dietary supplementation with phytogenics inhibited coliform bacteria; specifically Escherichia genus was practically undetectable in allochthonous microbiota of fish fed PHYTO and GMOSPHYTO diets. Accordingly, a reduction in Coliforms and in Escherichia coli counts were reported in the cecum, small and large intestine of broilers fed with diets containing phytogenic feed additives [7678].

In the last two decades EOs have been extensively investigated for their use as feed additives for terrestrial animals [79, 34]; only limited information is instead available about their potential benefits on fish health [41]. Studies reporting the effect of dietary EOs on gut bacterial composition of fish are indeed, still scarce, therefore our study represents a contribution. Consistent with our results, Ran et al. [65] tested the effects of a commercial blend of thymol and carvacrol essential oils in hybrid tilapia finding a reduction of E. coli/Coliform in the gut microbiota. Similarly, in another study in tilapia (Oreochromis niloticus), supplementation of diet with EOs from lemongrass, and geranium significantly decreased intestinal total bacteria, Coliforms, E. coli, and Aeromonas counts [12].

Numerous studies have demonstrated the antimicrobial effects of EOs in vitro, too. Among these, Helander et al. [80] demonstrated bactericidal properties of EOs from oregano, carvacrol, and thymol versus E. coli and Salmonella, whereas Ankri and Mirelman [81] described the wide spectrum antibacterial activity of garlic against Gram-negative and Gram-positive bacteria, including E. coli. The antimicrobial mechanism of essential oils is related to their lipophilic properties and chemical composition. EOs can damage bacteria cells affecting both, the membrane and cytoplasm. Indeed, the lipophilic nature of EO compounds allows them to penetrate cell membrane and remain between the phospholipid bilayer thus changing bacterial membrane structure and, consequently, its permeability. EOs can also affect quorum sensing systems [82, 43]. However, although, the antimicrobial properties of EOs have been clearly evidenced either in vitro or in vivo, it is also true that the commonly used microbiome profiling methods focus on the relative abundance or proportions of OTUs, meaning that that they cannot discriminate whether an enrichment in certain species is due to an increase in absolute abundance or a decrease in the abundances of other dominant taxa.

In the present study, changes in gut mucosa-adhered (resident) microbial communities were definitely less pronounced in comparison to dietary modulation of transient gut microbiota. This is a direct consequence of the fact that feedstuffs are a major source of allochthonous bacteria, which can temporarily integrate into the gut transient microbiome. However, it does not mean that the composition of transient bacterial communities is simply a mirror of feed-borne bacteria. Indeed, microbial profile of feeds resulted different from faecal microbial profiles. Otherwise, the data shows that even though there was a very large relative abundance of LAB in the diets, they failed to colonize/establish in the gut. They were not observed in the allochthonous bacteria, whereas in the autochthonous, they only represented an important order in the GMOSPTHYTO group and to a lesser extent in the GMOS group.

At resident microbiota level, only the combination of both functional ingredients (GMOSPHYTO) caused a significant increase of Firmicutes, essentially due to the enrichment in bacteria belonging to the Clostridiales order. In addition to Clostridium, Firmicutes phylum includes different genera of lactic acid bacteria (LAB) such as Streptococcus, Lactobacillus, and Leuconostoc. LAB are generally considered beneficial microorganisms and used as probiotics for fish and for other vertebrates [83, 84]. Therefore, an increase in their number is usually considered desirable, because they are associated with a healthy intestine. However, in sea bass, the number of LAB did not increase in response to dietary administration of prebiotic GMOS. This was an unexpected result since galactomannan-oligosaccharides are complex and indigestible carbohydrates (dietary fiber) that can be used as a substrate by LAB. Previous studies using similar feed additives reported an increase in LAB number in sea bream (Sparus aurata), Nile tilapia (O. niloticus), and rainbow trout (O. mykiss) gut microbiota [18, 8586]. Similarly, an enrichment in Lactobacillus genus was found in the intestine of trout fed with a diet containing insect meal, which is rich in chitin, a mucopolysaccharide polymer, which acts as a prebiotic being hardly digested by many fish species [48]. The combined administration of GMOS and phytogenic (GMOSPHYTO) also caused an enrichment in the Bacteroides and Prevotella genera in mucosa-associated microbial communities of sea bass. In humans, Bacteroides is clinically important genus, whereas the Prevotella species are commensal colonizers at mucosal sites, generally characterized by low pathogenicity and only few strains belonging to this genus exhibit opportunistic properties by promoting inflammatory diseases [87]. Interestingly, in line with our results, the number of counts of Prevotella was greater among mice fed galactomannan derived from citrus and fenugreek than among normal diet fed animals [88]. In the same study the increased amount of Prevotella was associated to an improved glucose metabolism.

The aforementioned differences in bacterial communities were also confirmed by the unweighted and weighted UniFrac analyses. PCoA plots showed that the major effect of diet was observed in faecal samples in terms of both, relative abundance and presence of specific taxa.

Conclusions

In summary, CTRL feed-associated bacteria showed a higher diversity in comparison to the bacteria found in functional feeds. We hypothesize that the additives GMOS and PHYTO, having an antimicrobial effect, have favoured the dominance of specific taxa, such as LAB thus lowering the biodiversity of functional feeds. However, the high amount of Lactobacillus found in functional feeds is in terms of relative abundance and it is not an absolute quantification. It means that we cannot discriminate whether an enrichment in certain species is due to an increase in absolute abundance or a decrease in the abundances of other dominant taxa.

As for the gut-associated microbiota, our findings suggest that the dietary inclusion of GMOS (0.5%) and a phytogenic composed by garlic and labiatae plants oils (0.02%) in a low FM and FO diet induces changes in bacterial community composition of European sea bass. However, if on allochthonous microbiota the combined inclusion of GMOS and PHYTO showed an antagonistic effect, at mucosa level only GMOSPHYTO diet increased the relative abundance of the Bacteroidales, Lactobacillales and Clostridiales resident microbiota orders. The main beneficial effects of GMOS and PHYTO on gut microbiota are on one side, the reduction of coliforms and Vibrionales bacteria, which include several potentially pathogenic species for fish, and the other, the enrichment of gut microbiota composition with butyrate producer taxa. Therefore, it is evident that these functional ingredients have a great potential to be used as health-promoting agents in in the farming of European sea bass and other marine fish.

Supporting information

S1 Dataset. List of OTUs found in two macro-groups of samples.

https://doi.org/10.1371/journal.pone.0231494.s001

(XLSX)

S1 Fig. Alpha diversity rarefaction curves of observed OTUs index of faeces and feed samples.

https://doi.org/10.1371/journal.pone.0231494.s002

(TIFF)

S2 Fig. Alpha diversity rarefaction curves of observed OTUs index of mucosa and feed samples.

https://doi.org/10.1371/journal.pone.0231494.s003

(TIF)

S1 Table. Original number of reads assigned to OTUs and alpha diversity metrics values (rarefied at 4500 reads) of feed-associated microbial community.

Reported data are expressed as means ± SD (n = 4). The means were compared by Kruskal-Wallis test (p<0.05). Different superscript letters on the same column indicate significant differences.

https://doi.org/10.1371/journal.pone.0231494.s004

(DOCX)

S2 Table. Mean relative abundance (%) ± SD of the most prevalent phyla, orders, classes, families, and genera found in feed samples.

Means in the same row with different letters indicate statistical significance between taxonomic groups’ abundances (p<0.05).

https://doi.org/10.1371/journal.pone.0231494.s005

(DOCX)

Acknowledgments

Special thanks to Delacon Biotechnik GmbH and to BioMar A/S, which provided the functional products and diets, respectively.

References

  1. 1. Shepherd CJ, Jackson AJ. Global fishmeal and fish-oil supply: inputs, outputs and markets. J. Fish. Biol. 2013; 83:1046–1066. pmid:24090562
  2. 2. FAO. World Fisheries and Aquaculture, 2018.
  3. 3. Francis G, Makkar HPS, Becker K. Anti-nutritional factors present in plant derived alternate fish feed ingredients and their effects in fish. Aquaculture 2001; 199:197–227. https://doi.org/10.1016/S0044-8486(01)00526-9
  4. 4. Krogdahl Å, Penn M, Thorsen J, Refstie S, Bakke AM. Important antinutrients in plant feedstuffs for aquaculture: an update on recent findings regarding responses in salmonids. Aquac. Res. 2010; 41:333–344. https://doi.org/10.1111/j.1365-2109.2009.02426.x
  5. 5. Øverland M, Sørensen M, Storebakken T, Penn M, Krogdahl Å, Skrede A. Pea protein concentrate substituting fish meal or soybean meal in diets for Atlantic salmon (Salmo salar). Effect on growth performance, nutrient digestibility, carcass composition, gut health, and physical feed quality. Aquaculture 2009; 288:305–311. https://doi.org/10.1016/j.aquaculture.2008.12.012
  6. 6. Torrecillas S, Monpel D, Caballero MJ, Montero D, Merrifield D, Rodies A, et al. Effect of fishmeal and fish oil replacement by vegetable meals and oils on gut health of European sea bass (Dicentrarchus labrax). Aquaculture 2017; 468:386–398. https://doi.org/10.1016/j.aquaculture.2016.11.005
  7. 7. Benedito-Palos L, Navarro JC, Sitjà-Bobadilla A, Bell JG, Kaushik S, Pérez-Sánchez J. High levels of vegetable oils in plant protein-rich diets fed to gilthead sea bream (Sparus aurata L.): growth performance, muscle fatty acid profiles and histological alterations of target tissues. Br. J. Nutr. 2008; 100:992–1003. pmid:18377678
  8. 8. Ringø E, Myklebust R, Mayhew TM, Olsen RE. Bacterial translocation and pathogenesis in the digestive tract of larvae and fry. Aquaculture 2007; 268:251–264. https://doi.org/10.1016/j.aquaculture.2007.04.047
  9. 9. Torrecillas S, Montero D, Izquierdo M. Improved health and growth of fish fed mannan oligosaccharides: Potential mode of action. Fish Shellfish Immunol. 2014; 36:525–544. pmid:24412165
  10. 10. Song SK, Beck BR, Kim D, Park J, Kim J, Kim HD, et al. Prebiotics as immunostimulants in aquaculture: A review. Fish Shellfish Immunol. 2014; 40:40–48. pmid:24973515
  11. 11. Rimoldi S, Finzi G, Ceccotti C, Girardello R, Grimaldi A, Ascione C, et al. Butyrate and taurine exert a mitigating effect on the inflamed distal intestine of European sea bass fed with a high percentage of soybean meal. Fish. Aquat. Sci. 2016; 19:40. https://doi.org/10.1186/s41240-016-0041-9
  12. 12. Al-Sagheer AA, Mahmoud HK, Reda FM, Mahgoub SA, Ayyat MS. Supplementation of diets for Oreochromis niloticus with essential oil extracts from lemongrass (Cymbopogon citratus) and geranium (Pelargonium graveolens) and effects on growth, intestinal microbiota, antioxidant and immune activities. Aquacult. Nutr. 2018; 24:1006–1014. https://doi.org/10.1111/anu.12637
  13. 13. Stevanović ZD, Bošnjak-Neumüller J, Pajić-Lijaković I, Raj J, Vasiljević M. Essential oils as feed additives-future perspectives. Molecules 2018; 23:1717. pmid:30011894
  14. 14. Encarnação P. Functional feed additives in aquaculture feeds. In: Aquafeed formulation (ed. by Nates S. F.) Elsevier, Amsterdam 2016; 217–236. https://doi.org/10.1016/B978-0-12-800873-7.00005-1
  15. 15. Grisdale-Helland B, Helland SJ, Gatlin DM. The effects of dietary supplementation with mannanoligosaccharide, fructooligosaccharide or galactooligosaccharide on the growth and feed utilization of Atlantic salmon (Salmo salar). Aquaculture 2008; 283:163–167. https://doi.org/10.1016/j.aquaculture.2008.07.012
  16. 16. Ringø E, Olsen RE, Gifstad TØ, Dalmo RA, Amlund H, Hemre GI, et al. Prebiotics in aquaculture: a review. Aquacult. Nutr. 2010; 16:117–36. https://doi.org/10.1111/j.1365-2095.2009.00731.x
  17. 17. Zhou Q-C, Buentello JA, Gatlin DM. Effects of dietary prebiotics on growth performance, immune response and intestinal morphology of red drum (Sciaenops ocellatus). Aquaculture 2010; 309:253–257. https://doi.org/10.1016/j.aquaculture.2010.09.003
  18. 18. Dimitroglou A, Merrifield DL, Spring P, Sweetman J, Moate R, Davies SJ. Effects of mannan oligosaccharide (MOS) supplementation on growth performance, feed utilisation, intestinal histology and gut microbiota of gilthead sea bream (Sparus aurata). Aquaculture 2010; 300(1): 182–188. https://doi.org/10.1016/j.aquaculture.2010.01.015
  19. 19. Guerreiro I, Serra CR, Pousão-Ferreira P, Oliva-Teles A, Enes P. Prebiotics effect on growth performance, hepatic intermediary metabolism, gut microbiota and digestive enzymes of white sea bream (Diplodus sargus). Aquacult. Nutr. 2018a; 24:153–163. https://doi.org/10.1111/anu.12543
  20. 20. Torrecillas S, Terova G, Makol A, Serradell A, Valdenegro V, Gini E, et al. Dietary phytogenics and galactomannan oligosaccharides in low fish meal and fish oil-based diets for European sea bass (Dicentrarchus labrax) juveniles: Effects on gut health and implications on in vivo gut bacterial translocation. 2019; PlosOne. 14(9): e0222063. https://doi.org/10.1371/journal.pone.0222063 pmid:31532807
  21. 21. Gibson GR, Roberfroid MB. Dietary modulation of the human colonic microflora introducing the concept of prebiotics. J. Nutr. 1995; 125:1401–1412. pmid:7782892
  22. 22. Torrecillas S, Makol A, Caballero MJ, Montero D, Robaina L, Real F, et al. Immune stimulation and improved infection resistance in European sea bass (Dicentrarchus labrax) fed mannan oligosaccharides. Fish Shellfish Immunol. 2007; 23:969–981. pmid:17766145
  23. 23. Torrecillas S, Makol A, Caballero MJ, Montero D, Ginés R, Sweetman J, et al. Improved feed utilization, intestinal mucus production and immune parameters in sea bass (Dicentrarchus labrax) fed mannan oligosaccharides (MOS). Aquacult. Nutr. 2011; 17:223–233. https://doi.org/10.1111/j.1365-2095.2009.00730.x
  24. 24. Torrecillas S, Makol A, Betancor MB, Montero D, Caballero MJ, Sweetman J, et al. Enhanced intestinal epithelial barrier health status on European sea bass (Dicentrarchus labrax) fed mannan oligosaccharides. Fish Shellfish Immunol. 2013; 34:1485–1495. pmid:23528875
  25. 25. Torrecillas S, Montero D, Caballero MJ, Robaina L, Zamorano MJ, Sweetman J, et al. Effects of dietary concentrated mannan oligosaccharides supplementation on growth, gut mucosal immune system and liver lipid metabolism of European sea bass (Dicentrarchus labrax) juveniles. Fish Shellfish Immunol. 2015a; 42:508–516. https://doi.org/10.1016/j.fsi.2014.11.033
  26. 26. Torrecillas S, Montero D, Caballero MJ, Pittman KA, Custódio M, Campo A, et al. Dietary mannan oligosaccharides: counteracting the side effects of soybean meal oil inclusion on European sea bass (Dicentrarchus labrax) gut health and skin mucosa mucus production? Front. Immunol. 2015b; 6:397 pmid:26300883
  27. 27. Torrecillas S, Caballero MJ, Montero D, Sweetman J, Izquierdo M. Combined effects of dietary mannan oligosaccharides and total fish oil substitution by soybean oil on European sea bass (Dicentrarchus labrax) juvenile diets. Aquacult. Nutr. 2016; 22:1079–1090. https://doi.org/10.1111/anu.12322
  28. 28. Torrecillas S, Rivero-Ramírez F, Izquierdo MS, Caballero MJ, Makol A, Suarez-Bregua P, et al. Feeding European sea bass (Dicentrarchus labrax) juveniles with a functional synbiotic additive (mannan oligosaccharides and Pediococcus acidilactici): An effective tool to reduce low fishmeal and fish oil gut health effects? Fish Shellfish Immunol. 2018; 81:10–20. pmid:29981880
  29. 29. Terova G, Forchino A, Rimoldi S, Brambilla F, Antonini M, Saroglia M. Bio-Mos®: An effective inducer of dicentracin gene expression in European sea bass (Dicentrarchus labrax). Comp. Biochem. Physiol. 2009; Part B, 153:372–377. pmid:19393760
  30. 30. Dimitroglou A, Merrifield DL, Carnevali O, Picchietti S, Avella M, Daniels C, et al. Microbial manipulations to improve fish health and production–A Mediterranean perspective. Fish Shellfish Immunol. 2011; 30:1–16. pmid:20801223
  31. 31. Guerreiro I, Oliva-Teles A, Enes P. Prebiotics as functional ingredients: focus on Mediterranean fish aquaculture. Rev. Aquac. 2018b; 10:800–832. https://doi.org/10.1111/raq.12201
  32. 32. Butt RL, Volkoff H. Gut microbiota and energy homeostasis in fish. Front. Endocrinol. 2019; (Lausanne) 10. https://doi.org/10.3389/fendo.2019.00009 30733706
  33. 33. Wong JM, de Souza R, Kendall CW, Emam A, Jenkins DJ. Colonic health: fermentation and short chain fatty acids. J. Clin. Gastroenterol. 2006; 40(3):235–43. pmid:16633129
  34. 34. Roberfroid M, Gibson GR, Hoyles L, McCartney AL, Rastall R, Rowland I, et al. Prebiotic effects: metabolic and health benefits. Br. J. Nutr. 2010; 104:S1–S63. pmid:20920376
  35. 35. Zeng Z, Zhang S, Wang H, Piao X. Essential oil and aromatic plants as feed additives in non-ruminant nutrition: A review. J. Anim. Sci. Biotechnol. 2015; 6(1):7. pmid:25774291
  36. 36. Harikrishnan R, Balasundaram C, Heo M. Impact of plant products on innate and adaptive immune system of cultured finfish and shellfish. Aquaculture 2011; 317:1–15. https://doi.org/10.1016/j.aquaculture.2011.03.039
  37. 37. Chakraborty SB, Horn P, Hancz C. Application of phytochemicals as growth-promoters and endocrine modulators in fish culture. Rev. Aquacult. 2014; 6:1–19. https://doi.org/10.1111/raq.12021
  38. 38. Reverter M, Bontemps N, Lecchini D, Banaigs B, Sasal P. Use of plant extracts in fish aquaculture as an alternative to chemotherapy: current status and future perspectives. Aquaculture 2014; 433:50–61. https://doi.org/10.1016/j.aquaculture.2014.05.048
  39. 39. Zaki MA, Labib EM, Nour AM, Tonsy HD, Mahmoud SH. Effect of some medicinal plants diets on mono sex Nile tilapia (Oreochromis niloticus), growth performance, feed utilization and physiological parameters. APCBEE Procedia 2012; 4:220–227. https://doi.org/10.1016/j.apcbee.2012.11.037
  40. 40. Awad E, Cerezuela R, Esteban MÁ. Effects of fenugreek (Trigonella foenum graecum) on gilthead seabream (Sparus aurata L.) immune status and growth performance. Fish Shellfish Immunol. 2015; 45:454–464. pmid:25956720
  41. 41. Sutili FJ, de Lima Silva L, Gressler LT, Gressler LT, Battisti EK, Heinzmann BM, et al. Plant essential oils against Aeromonas hydrophila: in vitro activity and their use in experimentally infected fish. J. Appl. Microbiol. 2015; 119(1):47–54. pmid:25810355
  42. 42. Sutili FJ, Gatlin DM, Heinzmann BM, Baldisserotto B. Plant essential oils as fish diet additives: benefits on fish health and stability in feed. Rev. Aquac. 2018; 10:716–726. https://doi.org/10.1111/raq.12197
  43. 43. da Cunha JA, Heinzmann BM, Baldisserotto B. The effects of essential oils and their major compounds on fish bacterial pathogens—a review. J. Appl. Microbiol. 2018; 125:328–344. pmid:29742307
  44. 44. Hüsnü K, Baser C, Demirci F. Chemistry of essential oils. In: Flavours and fragrances: chemistry, bioprocessing and sustainability. Edited by Berger RG. New York: Springer 2007; 43–86. http://dx.doi.org/10.1007/978-3-540-49339-6_4
  45. 45. Burt SA. Essential oils: their antibacterial properties and potential applications in foods–a review. Int. J. Food Microbiol. 2004; 94:223–253. pmid:15246235
  46. 46. Giannenas I, Triantafillou E, Stavrakakis S, Margaroni M, Mavridis S, Steiner T, et al. Assessment of dietary supplementation with carvacrol or thymol containing feed additives on performance, intestinal microbiota and antioxidant status of rainbow trout (Oncorhynchus mykiss). Aquaculture 2012; 350-353-26–32. https://doi.org/10.1016/j.aquaculture.2012.04.027
  47. 47. Rimoldi S, Gini E, Iannini F, Gasco L, Terova G. The effects of dietary insect meal from Hermetia illucens prepupae on autochthonous gut microbiota of rainbow trout (Oncorhynchus mykiss). Animals 2019; 9(4):143. pmid:30987067
  48. 48. Terova G, Rimoldi S, Ascione C, Gini E, Ceccotti C, Gasco L. Rainbow trout (Oncorhynchus mykiss) gut microbiota is modulated by insect meal from Hermetia illucens prepupae in the diet. Rev. Fish Biol. Fisheries 2019; 29:465–486. https://doi.org/10.1007/s11160-019-09558-y
  49. 49. Takahashi S, Tomita J, Nishioka K, Hisada T, Nishijima M. Development of a prokaryotic universal primer for simultaneous analysis of bacteria and archaea using next-generation sequencing. 2014; PLoS ONE 9(8), e105592. pmid:25144201
  50. 50. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet C, Al-Ghalith GA, et al. QIIME 2: Reproducible, interactive, scalable, and extensible microbiome data science. 2018; PeerJ Preprints 6, e27295v2. https://doi.org/10.7287/peerj.preprints.27295v2
  51. 51. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: High-resolution sample inference from Illumina amplicon data. 2016; Nat. Methods 13:581–583. pmid:27214047
  52. 52. Lozupone C, Knight R. UniFrac: a New Phylogenetic Method for Comparing Microbial Communities. Appl. Environ. Microbiol. 2005; 71:8228–8235. pmid:16332807
  53. 53. Lozupone CA, Hamady M, Kelley ST, Knight R. Quantitative and qualitative diversity measures lead to different insights into factors that structure microbial communities. Appl. Environ. Microbiol. 2007; 73:1576–1585. pmid:17220268
  54. 54. Hammer DAT, Ryan PD, Hammer Ø, Harper DAT. Past: Paleontological Statistics Software Package for Education and Data Analysis. Palaeontol. Electron. 2001; 4(1): 4–9. http://palaeo-electronica.org
  55. 55. Dawood MAO, Koshio S, Esteban MA. Beneficial roles of feed additives as immunostimulants in aquaculture: a review. Rev. Aquacult. 2018; 10: 950–974. https://doi.org/10.1111/raq.12209
  56. 56. Carda-Diéguez M, Mira A, Fouz B. Pyrosequencing survey of intestinal microbiota diversity in cultured sea bass (Dicentrarchus labrax) fed functional diets. FEMS Microbiol. Ecol. 2014; 87:451–459. pmid:24730648
  57. 57. Kim D-H, Brunt J, Austin B. Microbial diversity of intestinal contents and mucus in rainbow trout (Oncorhynchus mykiss). J. Appl. Microbiol. 2007; 102:1654–1664. pmid:17578431
  58. 58. Wu S, Gao T, Zheng Y, Wang W, Cheng Y, Wang G. Microbial diversity of intestinal contents and mucus in yellow catfish (Pelteobagrus fulvidraco). Aquaculture 2010; 303:1–7. https://doi.org/10.1016/j.aquaculture.2009.12.025
  59. 59. Rimoldi S, Terova G, Ascione C, Giannico R, Brambilla F. Next generation sequencing for gut microbiome characterization in rainbow trout (Oncorhynchus mykiss) fed animal by-product meals as an alternative to fishmeal protein sources. 2018; PLoS ONE 13(3):e0193652. pmid:29509788
  60. 60. Gatesoupe F-J, Huelvan C, Le Bayon N, Le Delliou H, Madec L, Mouchel O, et al. The highly variable microbiota associate to intestinal mucosa correlates with growth and hypoxia resistance of sea bass, Dicentrarchus labrax, submitted to different nutritional histories. BMC Microbiology 2016; 16:266. pmid:27821062
  61. 61. Etyemez M, Balcázarb JL, Demirkalea I, Dikela S. Effects of garlic-supplemented diet on growth performance and intestinal microbiota of rainbow trout (Oncorhynchus mykiss). Aquaculture 2018; 486: 170–174. https://doi.org/10.1016/j.aquaculture.2017.12.022
  62. 62. Sekirov I, Russell SL, Caetano L, Antunes M, Finlay BB. Gut microbiota in health and disease. Physiol. Rev. 2010; 90: 859–904. pmid:20664075
  63. 63. Apper E, Weissman D, Respondek F, Guyonvarch A, Baron F, Boisot P, et al. Hydrolysed wheat gluten as part of a diet based on animal and plant proteins supports good growth performance of Asian seabass (Lates calcarifer), without impairing intestinal morphology or microbiota. Aquaculture 2016; 453: 40–48. https://doi.org/10.1016/j.aquaculture.2015.11.018
  64. 64. Dimitroglou A, Merrifield DL, Moate R, Davies SJ, Spring P, Sweetman J, et al. Dietary mannan oligosaccharide supplementation modulates intestinal microbial ecology and improves gut morphology of rainbow trout, Oncorhynchus mykiss (Walbaum). J. Anim. Sci. 2009; 87: 3226–3234. pmid:19617514
  65. 65. Ran C, Hu J, Liu W, Liu Z, He S, Dan BC, et al. Thymol and carvacrol affect hybrid tilapia through the combination of direct stimulation and an intestinal microbiota-mediated effect: insights from a germ-free zebrafish model. J. Nutr. 2016; 146(5): 1132–40. pmid:27075912
  66. 66. Guardiola FA, Bahi A,Bakhrouf A, Esteban MA. Effects of dietary supplementation with fenugreek seeds, alone or in combination with probiotics, on gilthead seabream (Sparus aurata L.) skin mucosal immunity. Fish Shellfish Immunol. 2017; 65: 169–178. pmid:28433714
  67. 67. Torrecillas S, Makol A, Caballero MJ, Montero D, Dhanasiri AKS, Sweetman J, et al. Effect of mortality and stress response in European sea bass, Dicentrarchus labrax (L.), fed mannan oligosaccharides (MOS) after Vibrio anguillarum exposure. J. Fish Dis. 2012; 35: 591–602. pmid:22690841
  68. 68. Spring P, Wenk C, Dawson KA, Newman KE. The effects of dietary mannan oligosaccharides on cecal parameters and the concentrations of enteric bacteria in the cecae of Salmonella-challenged broiler chicks. Poult. Sci. 2000; 79:205–211. pmid:10735748
  69. 69. Ringø E, Strøm E, Tabachek JA. Intestinal microflora of salmonids: a review. Aquac. Nutr. 1995; 26: 773–789. https://doi.org/10.1111/j.1365-2109.1995.tb00870.x
  70. 70. Pryde SE, Duncan SH, Hold GL, Stewart CS, Harry JF. The microbiology of butyrate formation in the human colon. FEMS Microbiol. Lett. 2002; 217: 133–139. pmid:12480096
  71. 71. Esquivel-Elizondo S, Ilhan ZE, Garcia-Penã EI, Krajmalnik-Brown R. Insights into butyrate production in a controlled fermentation system via gene predictions. 2017; mSystems 2: e00051–17. pmid:28761933
  72. 72. Walker AW, Ince J, Duncan SH, Webster LM, Holtrop G, Ze X, et al. Dominant and diet-responsive groups of bacteria within the human colonic microbiota. ISME J. 2011; 5(2): 220–30. pmid:20686513
  73. 73. Ze X, Duncan SH, Louis P, Flint HJ. Ruminococcus bromii is a keystone species for the degradation of resistant starch in the human colon. ISME J. 2012; 6(8): 1535–43. pmid:22343308
  74. 74. Bonder MJ, Tigchelaa EF, Cai X, Trynka G, Cenit MC, Hrdlickova B, et al. The influence of a short-term gluten-free diet on the human gut microbiome. Genome Med. 2016; 8: 45. pmid:27102333
  75. 75. Ferreira-Halder CV, de Sousa Faria AV, Andrade SS. Action and function of Faecalibacterium prausnitzii in health and disease. Best Pract. Res. Cl. Ga. 2017; 31: 643–648. pmid:29566907
  76. 76. Cho JH, Kim HJ, Kim IH. Effects of phytogenic feed additive on growth performance, digestibility, blood metabolites, intestinal microbiota, meat color and relative organ weight after oral challenge with Clostridium perfringens in broilers. Livest. Sci. 2014; 160: 82–88. https://doi.org/10.1016/j.livsci.2013.11.006
  77. 77. Mountzouris KC, Paraskevas V, Tsirtsikos P, Palamidi I, Steiner T, Schatzmayr G, et al. Assessment of a phytogenic feed additive effect on broiler growth performance, nutrient digestibility and caecal microflora composition. Anim. Feed Sci. Tech. 2011; 168: 223–231. https://doi.org/10.1016/j.anifeedsci.2011.03.020
  78. 78. Franciosini MP, Casagrande-Proietti P, Forte C, Beghelli D, Acuti G, Zanichelli D, et al. Effects of oregano (Origanum vulgare L.) and rosemary (Rosmarinus officinalis L.) aqueous extracts on broiler performance, immune function and intestinal microbial population. J. Appl. Anim. Res. 2016; 44(1): 474–479. https://doi.org/10.1080/09712119.2015.1091322
  79. 79. Giannenas I, Bonos E, Christaki E, Florou-Paneri P. Essential oils and their applications in animal nutrition. Med. Aromat. Plants. 2013; 2: 1–12. https://doi.org/10.4172/2167-0412.1000140
  80. 80. Helander IM, Alakomi HL, Latva-Kala K, Mattıla-Sandholm T, Pol I, Smıd EJ, et al. Characterization of the action of selected essential oil components on gram-negative bacteria. J. Agric. Food Chem. 1998; 46: 3590–3595. https://doi.org/10.1021/jf980154m
  81. 81. Ankri S, Mirelman D. Antimicrobial properties of allicin from garlic. Microbes Infect. 1999; 1(2): 125–129. pmid:10594976
  82. 82. Nazzaro F, Fratianni F, De Martino L, Coppola R, De Feo V. Effect of essential oils on pathogenic bacteria. Pharmaceuticals (Basel) 2013; 6: 1451–74. pmid:24287491
  83. 83. Ringø E, Gatesoupe F-J. Lactic acid bacteria in fish: a review. Aquaculture 1998; 160: 177–203. https://doi.org/10.1016/S0044-8486(97)00299-8
  84. 84. Kim SK, Bhatnagar I, Kang KH. Development of marine probiotics: prospects and approach. Adv. Food Nutr. Res. 2012; 65: 353–362. pmid:22361199
  85. 85. Gonçalves AT, Gallardo-Escárate C. Microbiome dynamic modulation through functional diets based on pre- and probiotics (mannan-oligosaccharides and Saccharomyces cerevisiae) in juvenile rainbow trout (Oncorhynchus mykiss). J. Appl. Microbiol. 2017; 122: 1333–1347. pmid:28256031
  86. 86. Levy-Pereira N, Yasui GS, Vedovelli Cardozo M, Neto JD, Vaz Farias TH, Sakabe R, et al. Immunostimulation and increase of intestinal lactic acid bacteria with dietary mannan-oligosaccharide in Nile tilapia juveniles. R. Bras. Zootec. 2018; 47: e20170006. http://dx.doi.org/10.1590/rbz4720170006
  87. 87. Larsen JM. The immune response to Prevotella bacteria in chronic inflammatory disease. Immunology. 2017; 151(4): 363–374. pmid:28542929
  88. 88. Shtriker MG, Hahn M, Taieb E, Nyska A, Moallem U, Tirosh O, et al. Fenugreek galactomannan and citrus pectin improve several parameters associated with glucose metabolism and modulate gut microbiota in mice. Nutrition. 2018; 46: 134–142. pmid:28993009