Research article | Open | Published:
The highly variable microbiota associated to intestinal mucosa correlates with growth and hypoxia resistance of sea bass, Dicentrarchus labrax, submitted to different nutritional histories
BMC Microbiologyvolume 16, Article number: 266 (2016)
The better understanding of how intestinal microbiota interacts with fish health is one of the key to sustainable aquaculture development. The present experiment aimed at correlating active microbiota associated to intestinal mucosa with Specific Growth Rate (SGR) and Hypoxia Resistance Time (HRT) in European sea bass individuals submitted to different nutritional histories: the fish were fed either standard or unbalanced diets at first feeding, and then mixed before repeating the dietary challenge in a common garden approach at the juvenile stage.
A diet deficient in essential fatty acids (LH) lowered both SGR and HRT in sea bass, especially when the deficiency was already applied at first feeding. A protein-deficient diet with high starch supply (HG) reduced SGR to a lesser extent than LH, but it did not affect HRT. In overall average, 94 % of pyrosequencing reads corresponded to Proteobacteria, and the differences in Operational Taxonomy Units (OTUs) composition were mildly significant between experimental groups, mainly due to high individual variability. The highest and the lowest Bray-Curtis indices of intra-group similarity were observed in the two groups fed standard starter diet, and then mixed before the final dietary challenge with fish already exposed to the nutritional deficiency at first feeding (0.60 and 0.42 with diets HG and LH, respectively). Most noticeably, the median percentage of Escherichia-Shigella OTU_1 was less in the group LH with standard starter diet. Disregarding the nutritional history of each individual, strong correlation appeared between (1) OTU richness and SGR, and (2) dominance index and HRT. The two physiological traits correlated also with the relative abundance of distinct OTUs (positive correlations: Pseudomonas sp. OTU_3 and Herbaspirillum sp. OTU_10 with SGR, Paracoccus sp. OTU_4 and Vibrio sp. OTU_7 with HRT; negative correlation: Rhizobium sp. OTU_9 with HRT).
In sea bass, gut microbiota characteristics and physiological traits of individuals are linked together, interfering with nutritional history, and resulting in high variability among individual microbiota. Many samples and tank replicates seem necessary to further investigate the effect of experimental treatments on gut microbiota composition, and to test the hypothesis whether microbiotypes may be delineated in fish.
There is growing evidence that intestinal microbes play functional roles that are essential to health and nutrition, and a better understanding of the relationship between fish and their gut microbiota is crucial for sustainable aquaculture development . Many factors have been suggested as influencing the origin and composition of gut microbiota of fish, including genetic background [2, 3], diet , stress , and many environmental factors (e.g., temperature ). Regarding European sea bass, Dicentrarchus labrax, one of the two main marine fish species produced by aquaculture in southern Europe , high similarity was noticed among the faecal microbial communities collected from individuals grouped in the same tank . However, dominant and subordinate individuals of Arctic charr harboured distinct aerobic microbiota associated with intestinal mucosa , suggesting that social interaction may affect autochthonous gut microbiota in individuals cohabiting in the same tank. High-throughput sequencing methods have been recently applied to intestinal microbiota in fish. The pyrosequencing of 16S rRNA gene fragments allowed Roeselers et al.  to detect bacterial taxa that were shared by the intestinal communities in zebrafish of different origins, including specimens caught in the wild. The authors concluded that fish have a specific core intestinal microbiota, as is the case with higher vertebrates . However, the life history and diet may deeply influence gut microbiome in fish .
During the last decade, the composition of fish feeds has considerably changed to find the way to expand aquaculture without depleting natural fish stocks. Fish meal can be almost completely replaced by plant protein sources in the diet of sea bass . The replacement of fish oil by vegetable oils has been also achieved at least partially , or in short-term experiments , but the high requirement of sea bass for highly unsaturated fatty acids (HUFA) has hindered the simultaneous substitution of fish oil and fish meal by plant ingredients. Scarce information is available about the effect of dietary fatty acids on gut microbiota in fish, and most studies concerned salmonids . Starch and other polysaccharides from vegetal protein sources are known to influence gut microbiota in fish, including sea bass .
The composition of starter diets impacted the bacterial community associated with sea bass larvae , but the long-term effects of early nutritional treatments on intestinal microbiota have not yet been tested in this species. In rainbow trout fry, a short hyperglucidic hypoproteic dietary stress had short-term and long-term effects on intestinal fungi in juveniles, while intestinal bacteria were not significantly affected .
The present experiment was aimed at investigating the short-term and long-lasting effects of dietary stress on gut microbiota in sea bass fed diets either deficient in HUFA, or hypoproteic with starch as substitute for energy supply. To this end, sea bass were challenged with the unbalanced diets, or fed a standard diet, at first feeding (phase 1). After a transition period of 5 months, juvenile fish originating from the different dietary groups were transferred into the challenge tanks for phase 2. This second phase corresponded to a common garden experiment, which mixed in the same tanks animals that were already challenged or not during phase 1. These fish were fed the unbalanced diets for two months, while other individuals fed in standard conditions from start feeding onwards were maintained in a control tank. The physiological status of the individuals was addressed through growth and resistance to hypoxia. These health criteria were compared between dietary groups, and correlated with the composition of the intestinal community of autochthonous bacteria in each individual.
Individual variability and dissimilarity of intestinal microbiota
The microbial community was analysed in the intestinal mucosa of 54 fish after two days of fasting at the end of the experiment. The bacterial composition was computed from pyrosequencing data after RNA extraction and reverse transcription, with a view to compare the relative activity of every taxon. The number of valid reads per intestinal sample was highly variable, with significant differences between groups (Table 1). Normalization was thus critical before comparing the rarefaction curves and alpha-diversity. Within each experimental group, a wide variability appeared among samples. The rarefaction curves showed that the number of OTUs was still exponentially increasing even above 11,000 reads in some samples, whereas a rarefaction plateau appeared after much less reads in other ones [see Additional file 1]. Due to individual variability, there was no significant difference between the mean alpha-diversity indices of the experimental groups (Table 1).
The high variability could be further illustrated by comparing the numbers of OTUs shared among individuals and experimental groups. Over a total of 1111 OTUs, Escherichia-Shigella OTU_1 was the only taxon detected in every sample, with a relative abundance varying between 82.5 and 9.6 %. When the experimental groups were compared, 42 OTUs were shared by every group, while 135–159 OTUs were detected only in one group [see Additional file 2].
The intra-group similarity appeared relatively low in the three groups without transfer before the final dietary challenge, with a mean Bray-Curtis index of 0.52–0.53 (Table 1). Interestingly, the same index was differentially affected when the fish were transferred to the challenge tanks, depending on the diet. The individuals submitted to HUFA restriction only during phase 2 presented the lowest intra-group similarity in gut microbiota composition (group C1-LH2; mean Bray-Curtis index: 0.42; difference significant with every other group), while the highest similarity was observed among the fish transferred to the tank challenged with the starchy hypoproteic diet (group C1-HG2; mean Bray-Curtis: 0.60; difference significant only with groups HG1-HG2 and C1-LH2).
Despite individual variability, the non-parametric multivariate analysis of variance (PERMANOVA) revealed an overall significant difference between the mean bacterial profiles of the experimental groups, but none of the Bonferroni-corrected p values was less than 0.05 in the post-hoc pairwise comparison (Table 2). The overall significant difference was mainly due to the dissimilarity between the two groups submitted to HUFA deficiency, and between the control group reared in standard conditions (C1-C2) compared to each of the two groups transferred for phase 2.
Phylogenetic analysis of intestinal microbiota
The dominant phylum was Proteobacteria in every sample (94.4 ± 1.0 %), with mainly Gammaproteobacteria (56.5 ± 2.6 %) and Alphaproteobacteria (35.0 ± 2.9 %). Three other phyla were significantly represented: Bacteroidetes (2.3 ± 0.5 %), Actinobacteria (1.4 ± 0.3 %), and Firmicutes (1.1 ± 0.3 %). At the phylum level, the only significant difference after ANOVA corresponded to Actinobacteria, which were more abundant in group C1-C2, compared to C1-HG2, but some other differences were observed between groups after Linear discriminant Effective Size (LEfSe) pairwise analysis [see Additional file 3].
The proportion of 41 OTUs exceeded 5 % in at least one sample, with great variability among individuals [see Additional file 4]. At the OTU level, the most noticeable difference between groups concerned Escherichia-Shigella OTU_1, which was generally dominant, but significantly less prevalent in group C1-LH2 compared to group LH1-LH2 or C1-C2 by LEfSe, though there was no significant difference after ANOVA on every group, due to high individual variability (Fig. 1).
Pseudomonas OTU_3 was the only major OTU with a significant difference detected by ANOVA, with a mean relative prevalence of 14.7 % in the control group C1-C2, which represented more than the double of the mean proportions in the other groups. The hypoproteic diet with high starch supply had little effect on gut microbiota composition, which appeared rather similar between the two groups, with or without dietary challenge at first feeding. The most noticeable difference concerned Alkanindiges sp. OTU_14, which was not detected in group C1-HG2 [see Additional file 5].
There were other significant differences between groups at every phylogenetic level, as revealed by ANOVA, Kruskal-Wallis or LEfSe. It concerned some Gammaproteobacteria [see Additional file 5], Alpha- and Beta-Proteobacteria [see Additional file 6], Actinobacteria and Firmicutes [see Additional file 7], Bacteroidetes and Spirochaetae [see Additional file 8].
However, some of these differences appeared more visible between individuals than between groups. For example, three distinct OTUs of Bacillus sp. were relatively active, each one in three different samples from group C1-LH2 (OTU_178, OTU_196, and OTU_202, accounting for 1.2, 0.8, and 0.6 % of total reads in each of the three samples, respectively), while this genus was not detected elsewhere, except at a low level in one sample from group C1-HG2 (OTU_202, 0.05 % total reads).
Growth and resistance to hypoxia of sea bass
The highest mean weight corresponded to group HG1-HG2 at the three sampling dates, even at the beginning of phase 2 at 225 dph [see Additional file 9]. At this date, the mean weights of groups C1-C2, C1-LH2 and C1-HG2 were significantly lower than that of HG1-HG2. At 266 and 287 dph the mean weights of LH1-LH2 and C1-LH2 were significantly lower than that of HG1-HG2. The Specific Growth Rate (SGR) during this last three weeks was chosen as the best growth indicator, with a view to reduce the incidence of the initial differences in mean weights at the beginning of phase 2. During the final three weeks, the highest mean SGR was observed in the control group C1-C2 (Table 3). The lowest SGR were observed in the groups fed the HUFA-deficient diet, especially when the dietary challenge was applied at both phases (group LH1-LH2). The hypoproteic diet with high starch supply resulted in mild SGR.
At the end of the experiment, the group fed HUFA-deficient diet at both periods was significantly less resistant to hypoxia than those with normal HUFA dietary supply. The two groups that were fed unbalanced diets only during the last two months presented also a significant difference in mean Hypoxia Resistance Time (HRT), while the other pairwise comparisons did not indicate further significant differences (Table 3).
Correlation between gut microbiota composition and host’s growth or resistance to hypoxia
Due to some missing data, the complete correlation analysis could be performed only with 51 fish [see Additional file 10]. The individual data of SGR and HRT were put in front of the two indices used to describe alpha-diversity in gut microbiota. The means of OTU richness and dominance were not significantly different between experimental groups (Table 1), but disregarding the groups, the individual scores appeared dependent on the physiological traits. The Partial Least Squares (PLS) analysis showed that SGR correlated strongly with OTU richness, whereas HRT correlated with dominance (Fig. 2a). The two physiological traits correlated also with the relative prevalence of two distinct sets of OTUs (Fig. 2b). SGR correlated with the prevalence of Herbaspirillum sp. OTU_10 and Pseudomonas sp. OTU_3 (correlation scores around 0.8 and 0.6, respectively), while HRT correlated positively with Vibrio sp. OTU_7 and Paracoccus sp. OTU_4, but negatively with Rhizobium sp. OTU_9 (correlation scores around 0.8 and 0.6, and −0.8, respectively). The correlation between SGR and Pseudomonas sp. OTU_3 could be partly explained by the higher prevalence of the bacterium in the fast-growing control group, but the four other correlations seemed independent from the experimental conditions.
There is growing evidence indicating the impact of dietary components on fish gut microbiota, which seems essential for host health and well-being . However, apparent contradiction may arise between some observations obtained under different conditions. The present experiment attempted to evaluate the short and long-term influence of two kinds of nutritional deficiencies on the bacterial community associated with intestinal mucosa in sea bass individuals. When animals are subjected to nutritional stress, not all individuals may react in a similar way, partly depending on intestinal microbiota . A common garden approach was used to test the effects of possible social interaction and inter-individual contamination.
High-throughput sequencing has notably changed the insight on intestinal microbiota in fish, including European sea bass, which has been recently studied with such methods. Carda-Diéguez et al.  identified around 78 bacterial families in autochthonous intestinal microbiota of sea bass from genomic DNA, which was in the same range as the 90 families detected in the present dataset obtained from cDNA. However, the phylogenetic profiles were quite different with large proportions of Bacteroidetes, Tenericutes, and Firmicutes in the fresh samples analysed in the previous study. Surprisingly, after storage at −80 °C, Carda-Diéguez et al.  reported also a shift of the bacterial profile towards an overwhelming dominance of Proteobacteria. The main genera were not the same as in the present dataset, also dominated by Proteobacteria, but where Ralstonia sp. and Methylobacterium sp. accounted only for 0.05 and 0.02 % of total reads, respectively, and where Bradyrhizobium sp. was not even detected. A large drift of bacterial profiles during cold storage was unlikely in the present study, as the samples were immediately soaked in RNAlater before deep-freezing. Another potential source of bias may be due to resorting to nested PCR for mucosal samples with low concentration of bacterial 16S rRNA. Yu et al.  expressed reservation about the detection of rare OTUs after using nested PCR. The present analysis focused on the most active OTUs, which were likely the most susceptible to dietary influence. For this reason, the dataset was limited to the 1111 OTUs that represented at least 0.0002 % of total reads. The most prevalent OTU was the only taxon shared by every individual, but this did not necessarily contradict the core theory applied to fish intestinal microbiota , as the rare or weakly active bacteria could not be detected with the present method.
A striking feature of the present data was the high variability of the bacterial profiles among individuals reared in the same conditions. This may also happen in wild fish living in the same environment. Star et al.  displayed heterogeneous pyrosequencing profiles among the intestinal contents of 11 specimens of Atlantic cod, which were caught in one location, and then kept in a common tank for seven to twelve days of fasting. These results confirmed that gut microbiota composition may keep original features in fish individuals, even after mix in the same environment, in the absence of feed supply, which could blur the analysis done from intestinal contents. Some individual variation was already observed in the faecal microbiota from European sea bass, but the intestinal bacterial profiles of fish confined in the same tank were much more similar than those reared in replicate tanks . The study was based on the analysis of faeces that were emitted during the overnight isolation of individuals in separate aquaria, which led to conclude to a close intra-tank similarity between individual samples. After release, faecal microbiota evolved under the influence of nutrients and aquatic environment. Faeces collected after some delay cannot accurately reflect the intimate interaction between host intestinal mucosa and associated microbiota. Before comparing different microbial datasets, it is essential to distinguish the nature of the samples, which may come from faeces, digestive tract contents, or intestinal mucosa. The transient bacterial populations in the intestinal content (“allochthonous”) are known to differ from the community adhering to mucus in fish intestine (“autochthonous”) [23–25]. In the autochthonous community of sea bass intestine, Carda-Diéguez et al.  noted a higher diversity than in the samples of intestinal content, which appeared less reliable and less representative of the complexity of the bacterial consortium living in fish intestine. This later result confirmed the interest to focus on mucosal samples when studying the interaction between gut microbiota and the host, whereas the analysis of intestinal contents would be more dependent on the interaction with the diet.
Some relationship between individual characteristics and intestinal microbiota were reported in the literature. For example, the distribution of cultivable bacteria differed between the intestinal samples of dominant and subordinate Arctic charr . Two classes of orange-spotted with slow and fast growth rates presented two distinct bacterial profiles . Similar observations were done on Atlantic cod larvae, and less clearly on mangrove killifish larvae . Sea bass can also exhibit dominant behaviour, leading to different individual growth rates, which may account partly for the variability in intestinal microbiota, besides genetic factors. In flow-through water systems, there may be cross contamination between the bacterial communities in fish intestine, fish skin, and those adhering to tank wall. Microbiota associated with cutaneous mucus was influenced by the diet in Atlantic salmon , and by genetic background in brook charr . The interactions between fish individuals and bacterial communities within each tank may therefore cause partly the dissimilarity that was observed between intestinal microbiota collected from different tanks. The lack of replication of the experimental treatments in several tanks limited thus the conclusiveness about the effects of the diets in the present study. Some mean proportions of OTUs in the control group were significantly different from others in groups treated with a deficient diet, but it cannot be excluded that these differences might be due to tank cohabitation, rather than actually due to the diet (e.g., Pseudomonas sp. OTU_3).
The maintenance of discernible characteristics in gut microbiota composition after three months of common garden in fish individuals with two different nutritional histories was an innovative finding. The evidence came from the tank subjected to HUFA restriction during the final two months of dietary challenge, which was applied to fish either reared previously in standard conditions, or already deprived from normal HUFA supply at first feeding. Ringø et al.  “suggested that dietary fatty acids affect the attachment sites for the gastrointestinal microbiota, possibly by modifying the fatty acid composition of the intestine wall”. Dietary HUFA are known to influence the immune status of sea bass , and that might also impact the association of bacteria to intestinal mucosa. Bacterial colonization of marine fish larvae may be affected by dietary fatty acids , but the hypothetical long-term effect of such initial colonization remains to be investigated in fish. However, the restriction of dietary HUFA at start feeding compromised durably the growth potential of sea bass, as indicated by the lower SGR observed in the individuals already challenged during the larval stages, compared to those submitted to similar deficiency only during the final period. The long-term effect of the initial deficiency might also undermine the development of the immune system and the cell-wall defences in the intestine, which could explain why gut microbiota remained dissimilar between the two groups in the same tank. The most visible differences between these two groups lied in the highest individual variability and in the lowest average proportion of Escherichia-Shigella OTU_1 that were observed in the group transferred from the control tank for feeding the HUFA-deficient diet in phase 2. The high individual variability can be perceived through the diversity of the prevalent OTUs that competed with OTU_1 in the latter group, though their average proportions did not result in significant differences between groups. In the other groups, the apparent similarity of the relative proportion of OTU_1 might mask possible differences in bacterial load, especially when comparing the control group to those submitted to HUFA deficiency.
Besides HUFA, other feed components could affect the host-microbe interaction. In particular, the proportion of lupin meal was much higher in the HUFA-deficient diet than in the two others diets of the final challenge (Table 4). This alternative protein source is rich in non-starch polysaccharides, and it influenced gut microbiota composition in sea bass . It might induce the relatively high activity of Bacillus sp. in some samples, as these OTUs were close to strain DFEL3.4, previously isolated from the stomach of gilthead sea bream fed lupin meal in the same laboratory  [see Additional file 11].
The hypoproteic diet with high starch supply had little effect on gut microbiota composition. This was different from the results obtained previously in sea bass fed isoproteic diets with either lupin meal, or starch, or cellulose as carbohydrate source . The impact of the diet on gut microbiota was likely moderated in the present experiment by introducing lupin meal in every diet during the final challenge.
Beyond the experimental treatments, the correlative study allowed to highlight the interaction between gut microbiota composition and individual characteristics of the host, probably depending on genetic background and/or traits of behaviour. Such relationship should be carefully interpreted, as the correlation network is likely intricate. However it seems possible to associate physiological traits such as growth potential and hypoxia resistance to some characteristics of gut microbiota in fish, as already attempted between quantitative trait loci and specific bacterial strains associated with the skin of brook charr . Using DGGE, Forberg et al.  noted higher band richness in large killifish larvae compared to small individuals, but Shannon index and evenness were not significantly different, whereas the reverse observations were done with Atlantic cod larvae. Pyrosequencing allowed much deeper insight into bacterial diversity than DGGE, and the dominance index appeared clearly independent from OTU richness, which strongly correlated with SGR in sea bass. That suggested that the gut microbiota of fast-growing individuals might be more flexible, due to the presence of numerous taxa, which were not necessarily prevalent, but which could be activated in response to environmental changes, possibly benefiting to the host. The relationship between fish growth and gut microbiota may however depend on rearing conditions, especially on those aimed at managing the intestinal community. A synbiotic treatment with a probiotic strain of Lactococcus lactis and oligosaccharides increased growth in Siberian sturgeon, while decreasing gut microbial richness and Shannon index . More surprising was the strong correlation between the dominance index in gut microbiota and the host’s capacity to resist hypoxia. It is admitted in fish as in other vertebrates that the resistance to hypoxia depends on several physiological traits including strong capacity for metabolic depression and high energy reserves to fuel anaerobic metabolism [34, 35]. A longer exposure to hypoxia might favour the relative prevalence of Vibrio sp. OTU_7 in some resistant individuals, possibly by stimulating specific metabolic reactions as in Vibrio cholerae [36–38]. However, the fish were exposed to hypoxia for 6–8 h, and this range of variation seemed rather short to affect differentially intestinal microbiota. The diet might interfere and, as contrary to lupin meal, dietary starch stimulated Vibrio spp. in autochthonous gut microbiota of sea bass . In the present experiment, the HUFA-deficient diet, which depressed HRT, was rich in lupin meal with limited starch supply. Circumspection was nevertheless required to interpret partial least square correlation between HRT or SGR and the most prevalent OTUs, which were highly variable among individuals. For example, the negative correlation between HRT and Rhizobium sp. OTU_9 seemed mainly due to one individual, which showed the lowest HRT (6.04 h) and the highest percentage of OTU_9 (25.4 % total reads), whereas this OTU was detected only in seven individuals.
The main lesson is that active microbiota associated with intestinal mucosa may considerably vary among sea bass individuals, and large samples collected in several replicate tanks will be necessary to attempt at understanding the possible roles of gut microbiome. The severe HUFA restriction highlighted the interference between fish phenotype and gut microbiota, showing that its variability is not merely stochastic, but linked to life history or genetic background. This may suggest to investigate whether microbiotypes could be delineated in fish as in mammals , likely in terms of bacterial functions, rather than phylotypes.
Fish and rearing protocols
Newly-hatched sea bass larvae were provided by Aquastream (Ploemeur, France), allotted at 2 dph (day post hatch) in 15 tanks of 35 L, and then reared in the general conditions described elsewhere . The 15 tanks were divided in three groups of five tanks, which were fed the compound diets from 7 dph onwards (Fig. 3). Two unconventional diets, LH1 and HG1, were tested in comparison with a control diet C1 (Table 4). Diet LH1 consisted in low HUFA content by using defatted fish meal and soybean oil (less than 0.3 % EPA + DHA, dry matter basis). This diet was administered for the 22 first days of feeding, till 28 dph. Starch replaced 40 % of fish meal in the hypoproteic and hyperglucidic diet HG1, which was used from 7 to 22 dph, and then the larvae of this group were fed diet C1 till 28 dph. At 29dph, all the larvae were fed Artemia nauplii, with a view to compensate for the growth deficit caused by the deficient diets. At 35 dph, the larvae of the control group, fed diet C1 from 7 to 28 dph, were big enough to be grouped together in one 450-L tank (10,014 fish in total from the five replicated tanks), and then progressively weaned onto standard diet from 35 to 52 dph. Totals of 3662 and 6269 larvae previously fed LH1 and HG1, respectively, were transferred to two other 450-L tanks at 44 dph, and then progressively weaned onto standard diet till 52 dph. The water temperature was maintained at 20 °C during all the rearing period, in an open system without recirculation. At 191 dph, as the juveniles grew up, 500 individuals were randomly selected from each of the three dietary groups C1, LH1 and HG1 (individual mean weight of 16.6, 16.7, and 15.4 g, respectively). At 202 dph, a PIT-tag (PIT: passive integrated transponder) was subcutaneously implanted in every fish that was selected for the second phase of the experiment. Five groups of 100 tagged individuals were named after the diets fed during the two experimental phases (Table 4). The groups were formed as follows: three lots were randomly selected from the initial group fed diet C1, one still fed control diet (final group C1-C2), while the two other groups were mixed either with fish from group LH1 or HG1 for the common garden test. At 227 dph, the second phase of nutritional challenge started by feeding either low HUFA diet LH2 or high-starch/low-protein diet HG2 in each of the two tanks where two groups cohabited, namely, groups C1-LH2 and LH1-LH2 on the one hand, and groups C1-HG2 and HG1-HG2 on the other hand. HUFA were restricted in diet LH2 by replacing 65 % of the protein sources by lupin meal (c. 0.5 % HUFA, dry matter basis), while 30 % of the protein sources were replaced by starch in diet HG2. The fish were fed these experimental diets for two months, until the end of the experiment. They were individually weighed at 225, 266 and 287 dph after light anaesthesia (2-phenoxyethanol, 200 μL L−1). The individual specific growth rate (SGR) was computed between the last two weighing times (266 and 287 dph).
Final hypoxic challenge test
At 292 dph, the fish were grouped in one tank, and fasted for one day before being challenged for resistance to hypoxia. The method of this challenge test was previously described [41, 42]. Briefly, the rearing water was deprived of oxygen by bubbling nitrogen gas in the tank, under constant monitoring of dissolved O2. The oxygenation level was first dropped by 90 % from saturation within 1 h, and then further decreased by 1.2 % per hour, until reaching a minimum of c. 4 % air saturation. In the meantime, the fish were constantly observed, and each individual losing its maintenance of equilibrium was identified by reading its PIT-tag, and immediately placed in a fully aerated tank for recovery. The challenge lasted c. 8 h in total, and the interval of time from start to equilibrium lost was noted for every fish as its hypoxia resistance time (HRT).
Sampling for microbiological analysis
In each group, the fish were sorted based on HRT, and 12 individuals were selected per group for microbiological sampling (60 fish in total). Four individuals were selected among the most sensitive to hypoxia in each group, four others were among the mildly sensitive, and the four last were among the most resistant. As there were significant differences in hypoxia sensitivity between groups, the selection grid was not the same in every group [see Additional file 10].
At 294 dph, after two days of fasting, the 60 selected fish were euthanized with an overdose of anaesthetic (2-phenoxyethanol, 1 mL L−1). The intestines were empty, and each one was dissected under sterile conditions, and separated from the perivisceral fat in a Petri dish on ice. The intestine was immediately plunged into a microtube with 1.5 mL RNAlater (Qiagen). After soaking for 24 h at c. 4 °C, the tubes were stored at −80 °C.
RNA extraction, RT-PCR, and pyrosequencing
The microbial profiles compared in the present experiment were based on the analysis of 16S rRNA after reverse transcription, which was preferred to the method based on genomic DNA. As explained previously , this is a way to focus on the relative ribosomal activity among bacteria, whereas the relative abundance of rDNA cannot provide information about bacterial activity, and the result is biased by the variable number of gene copies among species.
After thawing, the intestine was removed from RNAlater with sterile tweezers, cut open longitudinally along the entire length in a Petri dish, and plunged into Extract-All (Eurobio) chilled on ice. RNA was extracted according to the instructions of the manufacturer for biological tissues, with an additional step of bead-beating for 10 min after the initial step of homogenization with a dispersing aggregate unit. After purification, the RNA concentration was estimated by NanoDrop (Thermo Scientific), and aliquoted for reverse transcription (RT). The surplus was precipitated and stored at −80 °C. cDNA was transcribed with the QuantiTect® RT kit (Qiagen), and stored at −20 °C.
Due to the small proportion of 16S rRNA in the samples, nested PCR was required before pyrosequencing. The PCR mix contained Taq DNA polymerase (0.025 U μl−1; MP Biomedicals), 0.2 mM of each dNTP (deoxyribonucleotide triphosphate; Eurogentec premix), and 0.4 μM of each primer (first round: EUB-8-f, 907-r; second round, V3-V4 region: PCR1F_460 and PCR1R_460 [see Additional file 12]). After initial denaturation at 94 °C for 2 min in T100 Thermal Cycler (Bio-Rad), 20 and 25 cycles were run in the first and second round, respectively, with 30 s denaturation at 94 °C, 30 s annealing at 55 °C (first round) or 62 °C (second round), and 1 min elongation at 72 °C. Both rounds ended with 1 min extension at 72 °C. Six of the sixty samples were ruled out because of insufficient PCR yield [see Additional file 10]. The 54 other PCR products were purified with GenElute PCR Clean-Up Kit (Sigma), further prepared by GenPhySE (INRA, UMR1388, Toulouse, France), and sequenced with Illumina MiSeq at GeT-PlaGe . The pyrosequencing data were deposited in the National Center for Biotechnology Information-Short Reads Archive (NCBI-SRA) under the BioProject accession number PRJNA294963 .
Bioinformatics data processing
The raw sequence dataset was first treated with FROGS (Find Rapidly OTU with Galaxy Solution) . Briefly, after merging the paired 250 bp reads, the software denoised the dataset, which was clustered with Swarm . A first range of chimera was removed with vsearch , and then the dataset was further filtered using PhiX and removing the singletons. A second filtration level was obtained by keeping only the clusters that represented at least 0.0002 % of total reads. After double identification with RDP and Blast + the dataset was restricted to the bacterial kingdom. A total of 1160 different sequences were thus detected among the 1,017,693 remaining reads. Between 11,599 and 28,067 valid reads were counted in each sample, and the data were normalized on the basis of 11,599 reads per sample, before computing alpha-diversity and the rarefaction curves. The normalization resulted in 1111 operational taxonomic units (OTUs). The final identification was assigned at the lowest phylogenetic level of RDP-Blast concordance, after correcting some misleading affiliations.
Fish growth and resistance to hypoxia, the diversity indices of intestinal microbiota, and the relative abundance of OTUs were compared between experimental groups by ANOVA or Kruskal-Wallis test, depending on normality and homoscedasticity. Post-hoc tests were used for multiple comparisons between dietary groups (Tukey’s and Dunn’s tests after ANOVA and Kruskal-Wallis, respectively). It must be noticed that the most relevant pairwise comparisons in the present experiment were those between the groups (1) with or without the deficient diet at first feeding and reared in the same tank during the final phase (LH1-LH2 vs. C1-LH2; HG1-HG2 vs. C1-HG2) and (2) with standard diet at first feeding, but either transferred or not from the control tank before the final phase (C1-LH2, C1-HG2, and C1-C2). The Bray-Curtis index was used for comparing the similarity between bacterial profiles by PERMANOVA with PAST . The bacterial profiles were further compared between two groups by Linear Discriminant Analysis (LDA) Effective Size (LEfSe) pairwise analysis under Galaxy environment [49, 50]. The canonical correlation between intestinal microbiota and fish growth or resistance to hypoxia was analysed after (sparse) Partial Least Squares, (s)PLS, classification with mixOmics [51, 52].
Analysis Of VAriance
Basic local alignment search tool
Control diet (C1 and C2, also used to name experimental groups)
Denaturing Gradient Gel Electrophoresis
Day post hatch
Find Rapidly OTU with Galaxy Solution
Hyper-Glucidic diet (HG1 and HG2, also used to name experimental groups)
Hypoxia Resistance Time
Highly poly-Unsaturated Fatty Acids
Linear Discriminant Analysis
Linear discriminant Effective Size
Low-HUFA diet (LH1 and LH2, also used to name experimental groups)
Operational Taxonomy Unit
PAleontological STatistics software package
Polymerase Chain Reaction
PERmutational Multivariate Analysis Of Variance
Partial Least Squares
Ribosomal Database Project
ribosomal Ribo-Nucleic Acid
Reverse Transcription Polymerase Chain Reaction
Specific Growth Rate
sparse Partial Least Squares
Nayak SK. Role of gastrointestinal microbiota in fish. Aquac Res. 2010;41:1553–73.
Navarrete P, Magne F, Araneda C, Fuentes P, Barros L, Opazo R, et al. PCR-TTGE analysis of 16S rRNA from rainbow trout (Oncorhynchus mykiss) gut microbiota reveals host-specific communities of active bacteria. PLoS One. 2012;7:e31335.
Forberg T, Sjulstad EB, Bakke I, Olsen Y, Hagiwara A, Sakakura Y, et al. Correlation between microbiota and growth in Mangrove Killifish (Kryptolebias marmoratus) and Atlantic cod (Gadus morhua). Sci Rep. 2016;6:21192.
Ringø E, Zhou Z, Vecino JLG, Wadsworth S, Romero J, Krogdahl Å, et al. Effect of dietary components on the gut microbiota of aquatic animals. A never-ending story? Aquacult Nutr. 2016;22:219–82.
Olsen R, Sundell K, Hansen T, Hemre GI, Myklebust R, Mayhew T, et al. Acute stress alters the intestinal lining of Atlantic salmon, Salmo salar L.: an electron microscopical study. Fish Physiol Biochem. 2002;26:211–21.
Leamaster BR, Walsh WA, Brock JA, Fujioka RS. Cold stress-induced changes in the aerobic heterotrophic gastrointestinal tract bacterial flora of red hybrid tilapia. J Fish Biol. 1997;50:770–80.
Federation of European Aquaculture Producers. Annual Report. 2015. p. 20. http://www.feap.info/shortcut.asp?FILE=1361. Accessed 28 Oct 2016.
De Schryver P, Dierckens K, Quyen QBT, Amalia R, Marzorati M, Bossier P, et al. Convergent dynamics of the juvenile European sea bass gut microbiota induced by poly-beta-hydroxybutyrate. Environ Microbiol. 2011;13:1042–51.
Ringø E, Olsen RE, Øverli Ø, Løvik F. Effect of dominance hierarchy formation on aerobic microbiota associated with epithelial mucosa of subordinate and dominant individuals of Arctic charr, Salvelinus alpinus (L.). Aquac Res. 1997;28:901–4.
Roeselers G, Mittge EK, Stephens WZ, Parichy DM, Cavanaugh CM, Guillemin K, et al. Evidence for a core gut microbiota in the zebrafish. ISME J. 2011;5:1595–608.
Turnbaugh PJ, Gordon JI. The core gut microbiome, energy balance and obesity. J Physiol (London). 2009;587:4153–8.
Sullam KE, Rubin BER, Dalton CM, Kilham SS, Flecker AS, Russell JA. Divergence across diet, time and populations rules out parallel evolution in the gut microbiomes of Trinidadian guppies. ISME J. 2015;9:1508–22.
Kaushik SJ, Coves D, Dutto G, Blanc D. Almost total replacement of fish meal by plant protein sources in the diet of a marine teleost, the European seabass, Dicentrarchus labrax. Aquaculture. 2004;230:391–404.
Ozsahinoglu I, Eroldogan T, Mumogullarinda P, Dikel S, Engin K, Yilmaz AH, et al. Partial replacement of fish oil with vegetable oils in diets for European seabass (Dicentrarchus labrax): effects on growth performance and fatty acids profile. Turk J Fish Quat Sci. 2013;13:819–25.
Castro C, Corraze G, Panserat S, Oliva-Teles A. Effects of fish oil replacement by a vegetable oil blend on digestibility, postprandial serum metabolite profile, lipid and glucose metabolism of European sea bass (Dicentrarchus labrax) juveniles. Aquacult Nutr. 2015;21:592–603.
Gatesoupe FJ, Huelvan C, Le Bayon N, Severe A, Aasen IM, Degnes KF, et al. The effects of dietary carbohydrate sources and forms on metabolic response and intestinal microbiota in sea bass juveniles, Dicentrarchus labrax. Aquaculture. 2014;422:47–53.
Delcroix J, Gatesoupe FJ, Desbruyères E, Huelvan C, Le Delliou H, Le Gall MM, et al. The effects of dietary marine protein hydrolysates on the development of sea bass larvae, Dicentrarchus labrax, and associated microbiota. Aquacult Nutr. 2015;21:98–104.
Geurden I, Mennigen J, Plagnes-Juan E, Veron V, Cerezo T, Mazurais D, et al. High or low dietary carbohydrate: protein ratios during first-feeding affect glucose metabolism and intestinal microbiota in juvenile rainbow trout. J Exp Biol. 2014;217:3396–406.
Tsukumo DM, Carvalho BM, Carvalho Filho MA, Saad MJA. Translational research into gut microbiota: new horizons on obesity treatment: updated 2014. Arch Endocrinol Metab. 2015;59:154–60.
Carda-Dieguez 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–9.
Yu G, Fadrosh D, Goedert JJ, Ravel J, Goldstein AM. Nested PCR biases in interpreting microbial community structure in 16S rRNA gene sequence datasets. PLoS One. 2015;10:e0132253.
Star B, Haverkamp THA, Jentoft S, Jakobsen KS. Next generation sequencing shows high variation of the intestinal microbial species composition in Atlantic cod caught at a single location. BMC Microbiol. 2013;13:248.
Bakke McKellep AM, Penn MH, Salas PM, Refstie S, Sperstad S, Landsverk T, et al. Effects of dietary soyabean meal, inulin and oxytetracycline on intestinal microbiota and epithelial cell stress, apoptosis and proliferation in the teleost Atlantic salmon (Salmo salar L.). Brit J Nutr. 2007;97:699–713.
Kim DH, Brunt J, Austin B. Microbial diversity of intestinal contents and mucus in rainbow trout (Oncorhynchus mykiss). J Appl Microbiol. 2007;102:1654–64.
Feng JB, Luo P, Dong JD, Hu CQ. Intestinal microbiota of mangrove red snapper (Lutjanus argentimaculatus Forsskål, 1775) reared in sea cages. Aquac Res. 2011;42:1703–13.
Sun YZ, Yang HL, Ling ZC, Chang JB, Ye JD. Gut microbiota of fast and slow growing grouper Epinephelus coioides. Afr J Microbiol Res. 2009;3:713–20.
Landeira-Dabarca A, Sieiro C, Álvarez M. Change in food ingestion induces rapid shifts in the diversity of microbiota associated with cutaneous mucus of Atlantic salmon Salmo salar. J Fish Biol. 2013;82:893–906.
Boutin S, Sauvage C, Bernatchez L, Audet C, Derome N. Inter individual variations of the fish skin microbiota: host genetics basis of mutualism? PLoS One. 2014;9:e102649.
Ringø E, Bendiksen HR, Gausen SJ, Sundsfjord A, Olsen RE. The effect of dietary fatty acids on lactic acid bacteria associated with the epithelial mucosa and from faecalia of Arctic charr, Salvelinus alpinus (L.). J Appl Microbiol. 1998;85:855–64.
Mourente G, Good JE, Thompson KD, Bell JG. Effects of partial substitution of dietary fish oil with blends of vegetable oils, on blood leucocyte fatty acid compositions, immune function and histology in European sea bass (Dicentrarchus labrax L.). Brit J Nutr. 2007;98:770–9.
Seychelles LH, Audet C, Tremblay R, Lemarchand K, Pernet F. Bacterial colonization of winter flounder Pseudopleuronectes americanus fed live feed enriched with three different commercial diets. Aquacult Nutr. 2011;17:E196–206.
Silva FCP, Nicoli JR, Zambonino-Infante JL, Kaushik S, Gatesoupe FJ. Influence of the diet on the microbial diversity of faecal and gastrointestinal contents in gilthead sea bream (Sparus aurata) and intestinal contents in goldfish (Carassius auratus). FEMS Microbiol Ecol. 2011;78:285–96.
Geraylou Z, Souffreau C, Rurangwa E, De Meester L, Courtin CM, Delcour JA, et al. Effects of dietary arabinoxylan-oligosaccharides (AXOS) and endogenous probiotics on the growth performance, non-specific immunity and gut microbiota of juvenile Siberian sturgeon (Acipenser baerii). Fish Shellfish Immunol. 2013;35:766–75.
Guppy M, Withers P. Metabolic depression in animals: physiological perspectives and biochemical generalizations. Biol Rev. 1999;74:1–40.
Bickler PE, Buck LT. Hypoxia tolerance in reptiles, amphibians, and fishes: life with variable oxygen availability. Annu Rev Physiol. 2007;69:145–70.
Marrero K, Sanchez A, Rodriguez-Ulloa A, Gonzalez LJ, Castellanos-Serra L, Paz-Lago D, et al. Anaerobic growth promotes synthesis of colonization factors encoded at the Vibrio pathogenicity island in Vibrio cholerae El Tor. Res Microbiol. 2009;160:48–56.
Fan FX, Liu Z, Jabeen N, Birdwell LD, Zhu J, Kan B. Enhanced interaction of Vibrio cholerae virulence regulators Tcpp and Toxr under oxygen-limiting conditions. Infect Immun. 2014;82:1676–82.
Hiremath G, Hyakutake A, Yamamoto K, Ebisawa T, Nakamura T, Nishiyama S, et al. Hypoxia-induced localization of chemotaxis-related signaling proteins in Vibrio cholerae. Mol Microbiol. 2015;95:780–90.
Waldram A, Holmes E, Wang Y, Rantalainen M, Wilson ID, Tuohy KM, et al. Top-down systems biology modeling of host metabotype-microbiome associations in obese rodents. J Proteome Res. 2009;8:2361–75.
Lamari F, Castex M, Larcher T, Ledevin M, Mazurais D, Bakhrouf A, et al. Comparison of the effects of the dietary addition of two lactic acid bacteria on the development and conformation of sea bass larvae, Dicentrarchus labrax, and the influence on associated microbiota. Aquaculture. 2013;376:137–45.
Roze T, Christen F, Amerand A, Claireaux G. Trade-off between thermal sensitivity, hypoxia tolerance and growth in fish. J Therm Biol. 2013;38:98–106.
Vanderplancke G, Claireaux G, Quazuguel P, Madec L, Ferraresso S, Sévère A, et al. Hypoxic episode during the larval period has long-term effects on European sea bass juveniles (Dicentrarchus labrax). Mar Biol. 2015;162:367–76.
Genotoul (Toulouse, France). Genome & transcriptome core facilities. http://get.genotoul.fr/fileadmin/user_upload/Presentation/poster_GenoToul_2012_english_combined.pdf. Accessed 28 Oct 2016.
fish gut metagenome. Intestinal microbiota of sea bass juveniles Metagenome. Long-term effects of early nutritional stimuli on growth, resistance to hypoxia, and intestinal microbiota in European sea bass, Dicentrarchus labrax, juveniles fed unconventional diets. Accession: PRJNA294963, ID: 294963. https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA294963. Accessed 28 Oct 2016.
Escudie F, Auer L, Cauquil L, Vidal K, Maman S, Mariadassou M, et al. FROGS: Find Rapidly OTU with Galaxy Solution. JOBIM Conf, 6–9 July 2015, Clermont-Ferrand, France. http://bioinfo.genotoul.fr/fileadmin/user_upload/FROGS_poster_Jobim_2015.pdf. Accessed 28 Oct 2016.
Mahé F, Rognes T, Quince C, de Vargas C, Dunthorn M. Swarm: robust and fast clustering method for amplicon-based studies. PeerJ. 2014;2:e593.
Rognes T, Mahé F, xflouris. vsearch: VSEARCH version 1.0.16. http://zenodo.org/record/15524. Accessed 28 Oct 2016.
Hammer Ø, Harper DAT, Ryan PD. PAST: paleontological statistics software package for education and data analysis. Palaeontol Electron. 2001;4. http://palaeo-electronica.org/2001_1/past/issue1_01.htm. Accessed 28 Oct 2016.
Segata N, Izard J, Walron L, Gevers D, Miropolsky L, Garrett W, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12:R60.
Galaxy/Hutlab. http://huttenhower.sph.harvard.edu/galaxy/. Accessed 28 Oct 2016.
Lê Cao KA, Martin PGP, Robert-Granié C, Besse P. Sparse canonical methods for biological data integration: application to a cross-platform study. BMC Bioinformatics. 2009;10:34.
Omics Data Integration Project. mixOmics. http://mixomics.org/. Accessed 28 Oct 2016.
This work was carried out with financial support from the Commission of the European Communities, specific RTD programme of Framework Programme 7 (project FP7-KBBE-2011-5, ARRAINA project no. 288925, Advanced Research Initiatives for Nutrition and Aquaculture). We are grateful to the Genotoul bioinformatics platform Toulouse Midi-Pyrénées, France, and to Sigenae group for providing help in pyrosequencing, computing and storage resources thanks to Galaxy instance (http://sigenae-workbench.toulouse.inra.fr). Special thanks are also due to Drs G. Pascal and O. Zemb, INRA, UMR1388 GenPhySE, for their precious help about pyrosequencing and bioinformatics.
This work was funded by the Commission of the European Communities, specific RTD programme of Framework Programme 7 (project FP7-KBBE-2011-5, ARRAINA project no. 288925, Advanced Research Initiatives for Nutrition and Aquaculture).
Availability of data and materials
The pyrosequencing data were deposited in the National Center for Biotechnology Information-Short Reads Archive (NCBI-SRA) under the accession number PRJNA294963 .
JLZ and DM conceived and designed the experiments on sea bass, and revised critically the manuscript. PQ and NLB were particularly involved in the rearing management of the larval and juvenile stages, respectively. HLD, LM, and OM were particularly involved in the biochemical analyses concerning diets and physiological data. CH was particularly involved in the statistical analyses regarding rearing and physiological data. FJG conceived and conducted the microbiological tasks, from intestinal sampling to bioinformatics and statistics, before writing the manuscript. All authors participated in data acquisition processes, and read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The present work followed the French and European guidelines on animal ethics (agreement number of the laboratory: 01964.01).
Rarefaction curves of OTU richness of each sample from the experimental groups, computed after normalization to 11599 reads per sample. (PPTX 916 kb)
Venn diagram showing the distribution of the 1111 OTUs among the intestinal samples of the five experimental groups. (PPTX 233 kb)
Histogram showing the overwhelming dominance of Alpha- and Gamma-Proteobacteria among the OTUs detected in every experimental group. (PPTX 65 kb)
Histogram showing the distribution of the most dominant OTUs (more than 5 % total reads in at least one sample). (PPTX 2597 kb)
Mean relative abundance of phylogenetic clusters among Gammaproteobacteria with significant differences between experimental groups. (DOCX 21 kb)
Mean relative abundance of phylogenetic clusters among Alpha- and Beta-Proteobacteria with significant differences between experimental groups. (DOCX 19 kb)
Mean relative abundance of phylogenetic clusters among Actinobacteria and Firmicutes with significant differences between experimental groups. (DOCX 19 kb)
Mean relative abundance of phylogenetic clusters among Bacteroidetes and Spirochaetae with significant differences between experimental groups. (DOCX 18 kb)
Mean weights (± SE) of sea bass before, at midterm, and after the final dietary challenge of sea bass. (DOCX 15 kb)
Hypoxia resistance time (h) of the individuals selected for the microbiological analysis. (DOCX 18 kb)
Phylogenetic tree derived by neighbour joining of the three Bacillus sp. OTUs from the present intestinal samples, aligned with GenBank sequences. (PPTX 64 kb)
Sequences of the primers used in the nested PCR for pyrosequencing. (DOCX 15 kb)