Dietary aquaculture by-product hydrolysates: impact on the transcriptomic response of the intestinal mucosa of European seabass (Dicentrarchus labrax) fed low fish meal diets

Aquaculture production is expected to double by 2030, and demands for aquafeeds and raw materials are expected to increase accordingly. Sustainable growth of aquaculture will require the development of highly nutritive and functional raw materials to efficiently replace fish meal. Enzymatic hydrolysis of marine and aquaculture raw materials could bring new functionalities to finished products. The aim of this study was to determine the zootechnical and transcriptomic performances of protein hydrolysates of different origins (tilapia, shrimp, and a combination of the two) in European seabass (Dicentrarchux labrax) fed a low fish meal diet (5%), for 65 days. Results were compared to a positive control fed with 20% of fish meal. Growth performances, anterior intestine histological organization and transcriptomic responses were monitored and analyzed. Dietary inclusion of protein hydrolysates in the low fish meal diet restored similar growth performances to those of the positive control. Inclusion of dietary shrimp hydrolysate resulted in larger villi and more goblet cells, even better than the positive control. Transcriptomic analysis of the anterior intestine showed that dietary hydrolysate inclusion restored a pattern of intestinal gene expression very close to the pattern of the positive control. However, as compared to the low fish meal diet and depending on their origin, the different hydrolysates did not modulate metabolic pathways in the same way. Dietary shrimp hydrolysate inclusion modulated more metabolic pathways related to immunity, while nutritional metabolism was more impacted by dietary tilapia hydrolysate. Interestingly, the combination of the two hydrolysates enhanced the benefits of hydrolysate inclusion in diets: more genes and metabolic pathways were regulated by the combined hydrolysates than by each hydrolysate tested independently. Protein hydrolysates manufactured from aquaculture by-products are promising candidates to help replace fish meal in aquaculture feeds without disrupting animal metabolism and performances.


Background
Fish captures have been stabilizing since the 1980's, and in most parts of the world seafood production is now supplied by the aquaculture industry [1]. Aquaculture production is expected to double by 2030 to meet customer demand, so demands for aquafeeds and raw materials will increase accordingly [2]. Originally, fish meals (FM) manufactured from wild fish were primarily used for growing carnivorous fish species. This kind of raw material is an excellent source of highly digestible protein and fat, is well balanced in essential amino acids, and has good palatability properties [3]. However, FM availability from wild fish is limited, its price is volatile, and its inclusion at high levels in aquafeeds is not sustainable; according to Duarte et al., FM and fish oil will be depleted by 2040 [4]. Consequently, sustainable growth of aquaculture will require the development of highly nutritive and functional raw materials to replace FM efficiently. Many studies have been conducted on carnivorous species to evaluate the performance of plant-based meal (PBM) diets to efficiently replace dietary FM; many of them were remarkably efficient, with partial or total substitution of FM by PBM [5][6][7][8]. However, replacing FM by high levels of PBM may also reduce feed palatability and fish growth [9,10]. Imbalanced amino acid composition of PBM leads to nutritional deficiency, and feed formulae incorporating high levels of such raw materials have to be supplemented in essential free amino acids [11]. It is also now well-documented that plants contain endogenous anti-nutritive compounds and complex carbohydrates that could affect nutrient digestibility and thereby negatively impact nutritional performances [12,13]. Moreover, high dietary inclusion of PBM in feeds for carnivorous species leads to enteritis [14,15], as well as depressed immunity [16,17]. Trials had been conducted about FM replacement by PBM on European seabass and it had been shown in different trials that it was possible to replace up to 50% of dietary FM by plant based meal without modifying fish growth performances (for review Kousoulaki et al., 2015) [18].
The large availability of marine and aquaculture by-products could allow for partial replacement of wild fish-derived FM [19,20]. By-products from the processing of fishery and aquaculture animals are more and more considered as a potential source of raw materials for sustainable FM production [21][22][23]. However, the quality of the processes applied to raw materials is critical to enhance their nutritional value. In this sense, enzymatic hydrolysis could improve the palatability, nutritional quality and functional properties of the finished product [24][25][26]. Enzymatic hydrolysis of fish proteins results in the formation of a mixture of free amino acids, di-, tri-and oligo-peptides, and enhances the occurrence of polar groups and the solubility of hydrolysate compounds [27]. Because low-molecular-weight nitrogenous compounds are important for the feeding behavior [28], nutrition [29][30][31] and health [32,33] of aquaculture species, protein hydrolysates could be good candidates for high FM substitution in aquafeeds. However, protein hydrolysate performances could be highly dependent on the methods used for their production: their nutritional and functional properties are closely related to their characteristics and composition, including the abundance and diversity of different oligo-peptides [34,35].
Dietary shrimp hydrolysates stimulate growth performances in fish, and also possess antimicrobial properties against aquaculture pathogens [36][37][38][39]. Similarly, tilapia hydrolysates as well as other fish hydrolysates have been evaluated; functional properties have been evidenced, such as antioxidant [24] and antimicrobial [40] activities. In particular, Khosravi et al. showed that high levels of FM could be replaced by low levels of protein hydrolysate combined with PBM in aquafeeds [38]. In addition, health benefits (i.e. immunity and gut cellular organization) of dietary protein hydrolysates have been reported in different fish species under challenging conditions [32,38,39,41].
The European seabass (Dicentrarchus labrax, Linnaeus, 1758) is a marine fish widely reared in the Mediterranean sea, with more than 156,449 tons produced in 2014 [1]. It is a strictly carnivorous species [42] that requires a high level of animal proteins in its diet. For instance, the natural diet of wild European seabass contains ca. 43-50% of animal proteins [43], whereas aquafeeds for this species generally contain at least 20% of FM to support good fish performances [44]. Thus, replacing FM in European seabass aquafeed is still a major objective for the aquaculture industry. Although the complete replacement of FM by PBM has been reported as possible for European seabass, feed formulation required palatability enhancers to improve feed intake, as well as amino acid supplementation to avoid nutritional deficiency linked to PBM [6]. Trials recently conducted in European seabass showed that 5% inclusion of shrimp-and tilapia-based protein hydrolysates associated with a combination of PBM successfully replaced 15% FM without affecting growth or health performances [45].
The intestine is involved not only in digestion and feed absorption, but also in water and electrolyte balance, nutrient sensing, and immunity. This functional diversity is gradually being elucidated in fish, as different histological and molecular approaches provide new items of knowledge regarding the many vital functions conducted along the gastrointestinal tract [46,47]. In this view, the numbers of transcriptomic studies in aquaculture have increased, mostly in the field of nutrition and immunity.
Many of them have focused on understanding how specific diets and functional ingredients could modulate metabolic pathways and regulate specific tissue expression [48][49][50]. But the effect of dietary protein hydrolysates on the regulation of fish metabolism has never been investigated so far.
We evaluated the effects of protein hydrolysates of different origins on European seabass fed a low FM diet. In addition to growth performances, we also studied the cellular organization and gene responses of the intestinal mucosa to investigate the effects of dietary protein hydrolysates on fish metabolism.

Diets
Five diets were formulated as follows: 2 diets containing FM at 5 and 20% of dry matter (diets FM5 and FM20) and 3 more diets containing 3 protein hydrolysates, shrimp-based hydrolysate (SH), tilapia hydrolysate (TH) and a 50/50 mixture of the two (MH), included at 5% of dry matter in the FM5 diet (Additional file 1). Both protein hydrolysates were provided by Aquativ (Diana Aqua, Symrise group, Elven, France). They were produced from the cephalothorax of white shrimp (Litopenaeus vannamei) and from Nile tilapia (Oreochromis niloticus) carcasses obtained from commercial food processing plants. These two protein hydrolysates have very different peptide profiles [37,40]. All the diets were balanced for deficient amino acids according to the requirements determined for European seabass [51]. Diets manufactured by BIOMAR (Tech Centre, Brande, Denmark) were extruded with 2 different pellet diameters ̶ 1.5 and 2.5 mm̶ for them to be adapted to the size of the fish during the trial. Diets were isoproteic (42.7 ± 1.1% of crude protein), isolipidic (19.3 ± 0.5% of crude fat) and isoenergetic (5.2 ± 0.1 kJ/kg).

Animals and feeding trial
The feeding trial was conducted at the Institut de Recerca i Tecnologia Agroalimentàries (IRTA, Sant Carles de la Rapita, Spain). European seabass (body weight, BW = 2.0 ± 0.2 g) were obtained from Piscicultura Marina Mediterránea SL (Andromeda Group, Burriana, Valencia, Spain). They were fed a commercial feed (OptibassL-2, Skretting, Spain; 48.5% of proteins, 16% of lipids, 3.7% of fibers, 6.4% of ashes) for 2 weeks for them to acclimate to the experimental facilities. Then they were randomly distributed into twenty 500-L fiberglass circular tanks (5 diets, 4 replicates per diet) at an initial density of 0.4 kg/m 3 (100 fish per tank). Before the feeding trial, sea bass were individually weighed for BW and measured for standard length (SL) (BW = 2.2 ± 0.01 g; SL = 5.1 ± 0.04 cm). During the study, average water temperature and pH (pH meter 507; Crison Instruments), salinity (MASTER-20 T; ATAGO Co. Ltd) and dissolved O 2 (OXI330; Crison Instruments) were 23.2 ± 0.5°C, 7.5 ± 0.2, 35.8 ± 0.3 ppm and 6.2 ± 1.2 mg/L, respectively. The water flow rate in the experimental tanks was maintained at approximately 9.0 L/min by means of a recirculation system (IRTAmar®) that maintained adequate water quality (ammonia: 0.08 ± 0.04 ppm, nitrites: 0.032 ± 0.02 ppm) through UV, biological and mechanical filtration. The photoperiod followed seasonal changes (February-April; latitude 40°37'41'N). Sea bass were fed in excess 6 times a day with automatic feeders (ARVO-TEC T Drum 2000™, Arvotec, Huutokoski, Finland) at a ration rate of 4.5%/day for 65 days. On day 31, the BW of 50 fish per tank was recorded to adjust the daily feed ration.
At the end of the trial, we measured the BW and SL of all sea bass from each tank (the fish were fasted for 24 h prior to sampling). Specific Growth Rates (SGR), Fulton's condition factors (K), feed conversion ratio (FCR) and Survival Rates (SR) were calculated as follows: SGR (% BW/day) = 100 x [ln final BW -ln initial BW]/duration of the trial (days); K = [BW/SL 3 ] × 100; FCR = [kg diet consumed]/[kg final biomass − kg initial biomass + kg sampled fish + mortalities]; SR (%) = [number of fish at the end of the trial/initial number of fish] × 100.

Ethics statement
All experimental procedures involving sea bass were conducted in compliance with the experimental research protocol approved by the Committee of Ethics and Animal Experimentation of the IRTA, the Departament Agricultura, Ramaderia, Pesca, Alimentació i Medi Natural (permit number 7962) and in accordance with the Guidelines of the European Union Council (86/609/ EU) for the use of laboratory animals.

Sample collection
During handling and weighing, sea bass were anesthetized with 50 mg/L of MS-222 (Sigma Aldrich, Saint-Louis, MO, USA). At the end of the trial, they were sacrificed with an overdose of MS-222 (100 mg/mL) to collect the different tissues for analytical purposes. Sea bass were sampled on a cold plate (0-4°C), and the intestines from 15 individuals per diet were dissected and preserved in a 10% phosphate-formaldehyde buffer (pH = 7.2) for histological purposes. In addition, the anterior intestine from nine fish per dietary treatment was sampled and immediately flash-frozen in liquid nitrogen for transcriptomic analysis. These samples were kept at − 80°C until RNA extraction.

Histological organization of the intestine
After fixation, samples were dehydrated in a graded series of ethanol, cleared with xylene, embedded in paraffin (Histolab ZX-60 Myr, Medical Specialties MYR SL, Tarragona, Spain) and cut into serial sections (2-3 μm thick) (HM Microm, Leica Microsystems, Nussloch, Germany). Sections were stained with Hematoxin-Eosin for general histological descriptions, whereas slides were stained with Periodic Acid Schiff (PAS) and Alcian blue (AB) at two different pH values (1.0 and 2.5) to stain different types of mucins produced by goblet cells [52]. PAS stains neutral mucins produced by intestinal goblet cells in magenta, whereas AB weakly stains ionised sulphated glycoconjugates at pH = 1.0 and sialic acid at pH = 2.5. All sections were observed under a light microscope (Leica DM LB; Leica Microsystems) and photographed (Olympus DP70 Digital Camera; Olympus Imaging Europa GmbH, Hamburg, Germany). Digital images were processed and analyzed using ANALYSIS image analysis software package (Soft Imaging Systems GmbH). Total numbers of goblet cells (full and empty) and villi height were measured based on the analysis of eight to ten randomly chosen fields from the intestinal mucosa of 15 sea bass per diet. Goblet cell counts in intestinal villi were expressed over a contour length of 100 μm, whereas villi height and width were calculated according to the method of Escaffre et al. [53].

Illumina sequencing
Total RNAs from a pool of the proximal intestines of three sea bass from the same replicate tank were extracted separately. Tissues were homogenized (Mini-Beadbeater, Biospec Products Inc., USA) in 1 mL of TRIzol (Ambion, Life Technologies, Carlsbad, CA, USA), and solvent extraction was performed following the manufacturer's instructions. Final RNA concentrations were determined by spectrophotometry (NanoDrop 2000; Thermo Fisher Scientific, Waltham, Massachusetts, USA). RNA quality was assessed from A 260 /A 280 ratios, and RNA integrity was assessed by denaturing gel electrophoresis. The three RNA samples from same replicate tank were then pooled before generating cDNA libraries. The total RNA concentration of each sample was quantified using a NanoDrop spectrophotometer (Thermo Fisher Scientific). cDNA libraries were built, and sequencing was performed as described in [54], with slight modifications. One μg of total RNA from each sample was initially used. dsDNAs were cleaved into 300-bp fragments using a Covaris S220 sonicator (Covaris Inc., Woburn, MA, USA) (duty cycle: 5%, intensity 3, 200 bursts per second, duration: 50 s). Eight pM of cDNA libraries per lane were loaded onto flow cells (Illumina Inc., San Diego, California, USA). The sequencing of 100-pb paired-end reads was performed on an Illumina Miseq Sequencer at the SéSAME Platform (Centre de Lutte Contre le Cancer François Baclesse, Caen, France).

Bioinformatic analysis
Bioinformatic analysis was performed on RNA sequencing data from European seabass anterior intestine to study their physiological and molecular responses to dietary hydrolysate inclusion in a low fish meal diet. We tested 5 dietary conditions on the transcriptomic response of sea bass intestine: FM5 and FM20 diets as negative and positive control groups, respectively, and the SH, TH and MH hydrolysate-based diets as experimental groups. The whole raw dataset was filtered and trimmed using Trimmomatic v(0.30) [55], using the following parameters: ILLUMINACLIP:adapter.fa:2:30:10, LEADING:5, TRAILING:5, SLIDING WINDOW:4:5, MINLEN:25. Global sequence quality was checked using FastQC (v 0.11.3) (https://www.bioinformatics.babraham.ac.uk/projects/fastqc). Global assembly was conducted using Trinity 2.1.1. [56], a dedicated package for de novo transcriptomics. A normalization step was conducted according to kmer coverage (kmer of 25 nt, maximum coverage of 30) proposed by the Trinity package. The inconsistent contigs generated by Trinity were removed after a remapping of reads using Bowtie (v 1.1.2) [57], and estimating relative abundance using RSEM [58] to get the FPKM values (v 1.2.22) (the two software programs were launched through perl wrappers provided by the Trinity package). Finally, only transcripts with at least an FPKM value above 1 and isoforms corresponding to more than 1% of the total gene count were kept. Annotation of contig sequences was performed using both Blast2Go software [59] and the Trinotate pipeline (http://trinotate.github.io), as described in [60]. Sequences were blasted against the NCBI nr database (release 193) with the following set-up parameters: max BLAST hits 20, min Expect Value 10 − 3 , and against the human proteome Ensembl database (release 82) (BLASTX). Only hits with E-values < 0.001 were kept. Peptide prediction was performed using Transdecoder [60]. Similarity search (blastp of the Transdecoder-predicted peptides) was performed against the uniprot-swissprot database (release 2015-05). Peptide signal prediction was performed using signalP v4.1 [61]. Transmembrane peptide detection was performed using TMHMM v2.0c [62]. Protein domain search was performed using hmmscan from the hmmer v.3.1b1 suite against the Pfam-A database [63]. Finally, transcriptome functional annotation was performed using the Trinotate pipeline. A Gene Ontology (GO) classification was assigned to each predicted protein in BLASTX (E-Value hit: 10 − 6 , annotation cut-off: 55, GO weight: 5). Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation was based on best BLASTX and PFAM search results.

Differential expression analysis
Differential gene expression between dietary treatments was identified using R software [64] and the DESeq2 package [65]. Multiple testing was accounted for by controlling the false discovery rate (FDR) at 5% using the Benjamini-Hochberg procedure. The RNA targets with adjusted P-values < 0.2 and absolute fold changes ≥1.4 were considered as differentially expressed. The relatively low stringency of the cut-off criteria is consistent with other nutrigenomic studies [48]. The impacts of hydrolysate diets on specific pathways were monitored using Ingenuity Pathway Analysis software (Qiagen, Hilden, Germany) with standard parameters. All annotated regulated genes were used.

Statistics
Final BW, SL and K results were expressed as mean ± standard error of the mean (SEM) from all fish from each treatment. SGR, intestinal villi height, number of goblet cells, proximal composition and survival rates were expressed as mean ± standard error of the mean (SEM) from the value of each replicate tank. All data were analyzed by one-way analysis of variance (ANOVA) followed by Tukey's test (normality and homogeneity of variances were previously checked). Differences were considered significant at P < 0.05. For GO term comparison with Wego tool, Pearson Chi-Square tests were performed to show significant differences between diets (P < 0.05) [66].

Zootechnical performances
Sea bass and feed performances are summarized in Table 1. At the end of the trial, dietary FM reduction (FM5 diet) had significantly impaired growth performances as compared to the positive control group (FM20) (P < 0.05). However, diets including protein hydrolysates (FM5 + 5% TH, FM5 + 5% SH and FM5 + 5% MH) significantly improved growth performances as compared to FM5 (P < 0.05) and allowed for similar growth performances as in sea bass fed the FM20 diet. The condition factor (K) and the feed intake were not affected by dietary FM reduction or dietary hydrolysate inclusion (P > 0.05). Even if FCR values were slightly improved in fish fed dietary hydrolysates compared to negative control (FM5), the only statistical difference was between the negative and the positive controls (P < 0.05). Dietary treatments did not impact fish survival (P > 0.05). Figure 1 shows the mean intestinal villi height of European seabass fed each diet. A significant decrease of villi height was observed in sea bass fed the FM5 diet as compared to the FM20 diet (P < 0.05). Inclusion of protein hydrolysates in the FM5 diet significantly increased intestinal villi height to reach values close to or even higher than in sea bass fed the FM20 diet (P < 0.05). The response of intestinal globlet cell density to dietary FM reduction and protein hydrolysate supplementation showed the same trend as intestinal villi height.

Global overview of the RNA-seq
The whole European seabass transcriptome project includes 108 Illumina libraries corresponding to the anterior intestine, but also the liver and kidney for each diet, for a total of 832,132,824 paired-end reads of 100 bp. Table 2 presents an overview of the sequencing project. Mapping was conducted both on the reference seabass genome (http://seabass.mpipz.de/DOWNLOADS/) and on the de novo assembled transcriptome from all samples, corresponding to 36 experimental conditions (6 diets * 3 tissues * 2 fish conditions) realized in triplicate.
Mapping scores were as follows: 83.55% for the mRNA extracted from the reference seabass genome annotation, and 92.77% for the de novo assembled transcriptome of the anterior intestine. In addition, the numbers of DE genes and annotations were higher with the de novo assembled transcriptome of the anterior intestine. Further differential expression analyses were then conducted on the global de novo transcriptome, composed of 56,246 transcripts after filtering with FPKM < 1 and isoform < 1%. The transcriptome corresponding to the total expression in

Quantitative analysis of differential gene expression
Differential gene expression analysis of RNA-seq data was performed to compare diets containing hydrolysates included to low (FM5) and high (FM20) fish meal diets. A total of 383 unique genes were differentially expressed ( Table 3). Compared to the FM5 diet, 46 and 50 genes were differentially regulated by diets containing TH and SH, respectively, and this number was higher in the FM20 (197) and MH (270) groups. Compared to the FM5 diet again, the distribution of differentially regulated genes between the FM20 diet and diets including hydrolysates indicated that 30.4% (14/46), 66.0% (33/50), and 39.3% (106/270) of regulated genes were shared with the TH, SH and MH diets, respectively. A total of 301 genes were only regulated by diets containing hydrolysates, but not by the FM20 diet, whereas some genes were regulated only by a specific hydrolysate diet. In particular, 19, 8, and 222 genes were specifically regulated by the TH, SH and MH diets, respectively. Finally, the TH, SH and MH diets induced gene expression patterns very similar to the FM20 diet pattern, with only 8, 6, and 6 genes differentially expressed among these dietary groups, respectively. The complete list and corresponding fold changes of regulated genes is provided in Additional files 2 and 3, for a comparison with the low-FM (FM5) and high-FM (FM20) diets, respectively.

Identification of regulated genes
GO term analysis revealed that dietary FM reduction, from 20 to 5%, led to the regulation of a wide variety of genes involved in the metabolism, immunity, and tissue development of European seabass intestinal mucosa (Fig. 2). GO term analysis of the genes regulated by hydrolysate-supplemented diets is presented in Fig. 3. The specific hydrolysates, which were influenced by by-product origin and the hydrolysis process, differentially impacted gene regulation (P < 0.05). In particular, under the TH diet, the level of GO term annotation was higher for the following categories: macromolecular complex, non-membrane-bounded organelle, intracellular non-membrane-bounded organelle, hydrolase activity, and ribonucleotide binding. In contrast, under the SH diet, the differentially regulated genes belonged mostly to the immune response and response to stimulus. Finally, under the MH diet, lipid metabolism-related GO terms were over-represented.
Based on GO term annotations, 74 differentially expressed transcripts were linked to nutrition as compared   (5). Under the FM20 diet, carrier genes were over-expressed, except for a few related to lipid transport (apolipoproteins). Diet formulations induced specific regulations with 12, 4, 4, and 32 transcripts differentially expressed by the FM20, SH, TH, and MH diets, respectively, as compared to FM5. Table 5 presents regulated transcripts related to the cholesterol and long-chain fatty acid pathways. Lipid metabolism was modulated by the MH diet even more than by the FM20 diet. As compared to the FM5 diet, 15, 0, 3, and 28 transcripts were differentially expressed by the FM20, TH, SH and MH diets, respectively. The TH diet did not impact the expression of the transcripts related to the cholesterol and long-chain fatty acid pathways.
Transcripts related to fish immunity and/or the stress response are presented in Table 6. As compared to the FM5 diet, 20, 4, 12, and 36 transcripts were differentially expressed by the FM20, TH, SH and, MH diets, respectively. In contrast to the FM20 diet, all the diets containing hydrolysate down-regulated the expression of transcripts related to interferon. The diets containing one kind of hydrolysate (TH or SH) did not impact the regulation of histo-compatibility antigen-related genes, and only slightly impacted the serum complement-and cellular damage-related genes. Lectin-related gene expression was not affected by dietary TH inclusion.

Impacted metabolic pathways
We conducted various analyses to study the metabolic pathways impacted by dietary hydrolysates. Firstly, we performed a string analysis (https://string-db.org/cgi/ input.pl?) to reveal relationships between regulated genes and a defined enrichment (Additional file 4). Regarding the TH diet, 11 genes, distributed into 3 groups, were interconnected, but no enrichment was observed. In the case of the SH diet, 11 regulated genes formed 2 groups, but without inducing overexpression of any specific metabolic pathway. Six regulated genes were common to the TH and SH diets (cebpd, cebpb, rsad2, herc6, mx1 and cmpk2). The g6pc, aldob, eif2ak4, ace, apob and hba1, scn2a, junb, fdp5, dhx58 genes were connected to the TH and SH diets, respectively. Under the MH diet, many more regulated genes were connected (87), with the ACLY gene at the center of the gene network. The enrichment scores calculated for biological processes, molecular functions, cellular component and KEGG pathways are shown in Table 7.
The main enrichments for biological processes were the following: small molecule metabolic process, cellular lipid metabolic process, lipid metabolic process, organophosphate metabolic process, single-organism catabolic process; cellular components were cytosol, cytoplasmic part, membrane-bounded vesicle, extracellular exosome, and extracellular region part; whereas molecular functions were catalytic activity, ethanolamine kinase activity, oxidoreductase activity, cofactor binding, coenzyme binding, and KEGG metabolic pathways were steroid biosynthesis, fat digestion and absorption, glycerophospholipid metabolism, and glycolysis/gluconeogenesis.
A KEGG analysis based on KEGG orthology annotation showed that the main impacted pathways differed depending on the origin of the dietary hydrolysate (Table 8). In particular, 9, 9, and 55 differentially regulated genes participated to metabolic pathways (KO: 01100) in response to the TH, SH, and MH diets, respectively. To be more precise, the metabolic pathways impacted by the TH diet mainly influenced carbohydrate metabolism, whereas the SH diet impacted the nucleotide metabolism and terpenoid backbone biosynthesis. In addition to these metabolic pathways, defense and immunity pathways were also affected by hydrolysate inclusion to the diet. The complete list of impacted pathways and corresponding genes is available in Additional file 5.
We also investigated the pathways regulated by dietary treatments using Ingenuity Pathway Analysis Software. Results are summarized in Table 9, focusing on the top canonical pathways as well as on molecular and cellular functions. The main metabolic pathways impacted by dietary hydrolysates were the ones related to the

Impact of hydrolysate inclusion in a low FM diet on European seabass growth performances and histological organization of the intestine mucosa
Our results showed that a 15% decrease of dietary FM negatively affected the growth performances of European seabass as compared to the 20% FM control diet. These findings are consistent with other studies about other carnivorous fish species [12,[67][68][69]. However, dietary inclusion of 5% of protein hydrolysate restored growth performances to the same level as the FM20 control diet. These positive effects of dietary hydrolysates are in agreement with other studies about fish fed diets incorporating such ingredients [29,32,38]. Feed intake was not modified by diets as fish were fed in excess to guarantee a high consumption. Feed conversion ratio improvement of fish fed the FM20 diet compared to the FM5 diet could be the result of a better nutritional balance of the diet combined to a better feed consumption. Dietary inclusion of hydrolysate allowed to recover a good feed intake as well as to improve enough feed efficiency to reach fish and feed performance very closed to the positive control. The high degree of protein hydrolysis of our hydrolysates [37,40], resulting in enrichment of the feed with highly palatable and digestible low-molecular-weight peptides and free amino acids, could explain why growth performances were improved. However, fish growth performances should not be the single criterion to evaluate the performance of a new aquafeed formula. It is also necessary to study the response of fish metabolism to ensure that modifications in the diet will not induce significant metabolic disturbances that could affect fish resistance to their environment. Maintaining the integrity of fish intestines is of critical concern when considering the performance of aquaculture feeds. Many studies have shown that enteritis events can occur when carnivorous fish are fed feeds formulated with high levels of plant-based proteins [17,[70][71][72][73][74]. A direct consequence of enteritis is deterioration of feed utilization and of the fish health status. The intestine is involved not only in digestion and feed absorption, but also in water and electrolyte balance, nutrient sensing, and immunity. This functional diversity is currently being elucidated in fish, and different histological and molecular approaches are helping to understand the many vital functions conducted along the gastrointestinal tract [46]. The general organization and morphometric parameters of the intestine, represented by villi size and the number of goblet cells, are good indicators of the health and condition of the fish intestinal mucosa: an increase of intestinal villi size reflects an improvement of the exchange surface, of the activity of the brush border enzymes and of the nutrient transport systems, with positive effects on digestion and absorption [75]. In addition, in rainbow trout (O. mykiss) goblet cells regulate proteins or peptides as well as ion and fluid transport, and also provide an effective immune barrier against potentially pathogenic gut bacteria [47,76]. Furthermore, enhanced mucin production by increased goblet cell populations physically displaces potentially pathogenic organisms, a more diverse microbiota leads to a thickening of the mucus layer, and this improves gut microniches inhabited by these beneficial bacteria [77]. In the present study, we demonstrate that European seabass fed a diet containing low fish meal levels (FM5) had an altered intestinal mucosa. These results are consistent with other studies that showed a decrease in villi height [39,72,[78][79][80] and/or in the number of goblet cells [38,39] in fish fed the low FM diet. Incorporation of hydrolysates at a relatively low inclusion rate (5%) in the low FM formula (FM5) improved the histological organization of the intestinal mucosa, with villi heights close to those observed in the control group (FM20). Moreover, the origin of the raw material used to manufacture the hydrolysates, as well as their hydrolysis level, seemed to be determining parameters for enhancing the morphological development of the intestine. The different specifications of the protein hydrolysates led to different responses of villi height, with a higher benefit from dietary SH than TH, although MH performance lay in-between. However, it is still unclear whether this improved response was related to a dietary effect of SH on the intestinal mucosa due to i) the higher intake of free amino acids, or ii) a protective effect of the intestinal epithelium due to the presence of bioactive peptides in the hydrolysate, or iii) a modulation of the intestinal microbiota and/or of antimicrobial peptides. These hypotheses are not mutually exclusive. Although the mode of action of hydrolysates remains to be elucidated, this study confirms that they are beneficial for the fish intestine, as already observed in olive flounder [39]. It also confirms the potential of hydrolysates for replacing FM in aquafeeds for carnivorous species.

Impact of hydrolysate inclusion in a low FM diet on intestine gene expression
The intestine is a complex organ because many metabolic interactions take place inside it [49]. We analyzed the differentially expressed genes from the intestinal mucosa of European seabass fed different diets. Low FM (FM5 diet) caused metabolic disturbances, and certain metabolic pathways were modulated in link with dietary inclusion of protein hydrolysates of different origins. Dietary FM reduction (FM5 diet) impaired several intestinal functions, i.e. nutrient transport, immune defense, gut morphogenesis (Fig. 2). These transcriptomic results are consistent with the lower performances recorded in sea bass fed the FM5 diet throughout the zootechnical trial. The different transcriptomic data analyses evidenced a few discrepancies in sea bass responses depending on dietary hydrolysate origin (Tables 8 and 9). The comparative analysis of the two SH and TH protein hydrolysates confirmed that the sea bass intestinal transcriptomic response seemed to be closely related to the biochemical properties of each hydrolysate. Thus, based on GO term analysis, a higher proportion of genes (34.1% vs. 12.7%) linked to hydrolases was differentially expressed in sea bass fed the TH diet as compared to the SH diet (Fig. 3). Interestingly, the regulated enzymes were not only implied in peptide/protein degradation, as expected from dietary protein hydrolysate inclusion, but also in carbohydrate degradation (malatase-gluco-amylase, aldolase B, glucose-6-phosphatase), and probably stimulate starch degradation in the PBM fraction of diets (Table 4). These differences in gene expression of intestinal hydrolases might be explained by the different peptide compositions of tilapia and shrimp hydrolysates. TH is mainly composed of higher-molecular-weight peptides than SH, and also exhibits higher peptide diversity than SH [40]. These higher-molecular-weight peptides could require additional hydrolysis by enterocyte enzymes, resulting in greater mRNA synthesis from genes related to protein digestion and absorption. On the contrary, GO terms of differentially expressed genes involved in the immune response (6.8% vs. 23.6%) and in the response to stimuli (34% vs. 56.4%) were less in sea bass fed the TH diet than in sea bass fed the SH diet (Fig. 3). In particular, the SH diet induced the expression of lectin-related genes (Table 6), which protect the intestinal mucosa from pathogenic bacterial invasion [81]. These findings support the immunostimulating effect of SH. SH is mainly composed of low-molecular-weight peptides derived from the enzymatic cleavage of haemocyanin [37], which has immunostimulatory properties used in the treatment of cancers such as melanomas or in bioadjuvants for vaccines [82][83][84]. In addition, peptides derived from haemocyanin degradation are also known to exhibit other bioactive properties such as responses to stressors, pathogens (bacteria and parasites) and antibacterial agents [85]. Marine protein hydrolysates have been reported to have in  vitro antibacterial [37,40], antioxidant and/or immunomodulatory [86] functional properties. Our results seem to confirm immuno-stimulatory effects (interferon, histocompatibility antigen, lectin, pentraxin or serum complement system expression) at a transcriptomic level, as already observed in several carnivorous fish fed diets containing the same hydrolysates [32,33,38,39,87]. Moreover, hydrolysates down-regulated the expression of interferon-related genes. In fish, genes linked to interferon are expressed during viral infection [81], and in vertebrates, interferons are also involved in intestinal homeostasis after infection [88]. These transcriptomic results confirm that dietary hydrolysate supplementation has beneficial effects on the intestinal mucosa.
String analysis did not evidence any GO term or metabolic enrichment in diets containing hydrolysates of a specific origin (shrimp or tilapia). This result is surprising given the marked difference in growth responses between sea bass fed the FM5 diet and sea bass fed diets containing hydrolysates, as well as the positive response of intestine histological organization in sea bass fed the SH diet. For the SH diet, this may be due to the fact that sea bass fed the FM5 diet did not show any histological signs of severe enteritis or inflammation of the intestinal mucosa; besides, there was no change at a molecular level. Enteritis induces the regulation of a large number of genes (apoptosis, proinflammatory, oxidative stress, endocytosis, and cell migration) as already shown in a study on the replacement of FM by soybean meal in Atlantic salmon [48]. In our study, the physiological state of the sea bass may not have been sufficiently altered to induce a transcriptomic response taking into account our filtering treatment (absolute fold change > 1.4 and P-adj < 0.2). Hence, the low number of regulated genes did not allow us to highlight any group of genes to explain the phenotypic differences observed between dietary treatments at a transcriptomic level.
On the contrary, in sea bass fed the diet containing the mixed hydrolysate (MH diet), the much larger number of differentially regulated genes (270 as compared to the FM5 diet) made it possible to carry out an enrichment analysis with String both in GO terms and in metabolic pathways. The analysis showed that the lipid metabolism was strongly impacted, and acetyl-CoA appeared to play a central role in this process. Interestingly, the MH diet decreased cholesterol and long-chain fatty acid metabolism (Table 5). Genes involved in cholesterol and long-chain fatty acid metabolism were overexpressed in Atlantic salmon fed protein diets containing a high level of soybean meal [48]. On the contrary, these two genes families were down-regulated in European seabass fed the MH diet. Moreover, many genes coding for apolipoprotein, which is involved in lipid transport, were down-regulated in sea bass fed with the mixed hydrolysate. These results provide evidence for a significant effect of the mixed hydrolysate diet on lipid metabolism at the transcriptomic level. Such an effect would deserve to be deeply investigated to understand the phenotypic consequences (lipid utilization and storage) of such gene responses to dietary protein hydrolysate. It is worth noting that modulation of lipid metabolism by dietary protein hydrolysates has already been observed in turbot and mice, with a beneficial effect on the reduction of visceral lipid accumulation [35,[89][90][91].
The MH diet also promoted European seabass health status by regulating the same immuno-stimulatory related genes as the TH and SH diets, but also by specifically modulating the expression of specific genes, such as H − 2 class I and class II histocompatibility antigen genes, which were over-expressed (Table 6). Moreover, the genes involved in the response to cellular damage (glutathione peroxidase, DNA damage repair) were under-expressed, which suggests an improvement of intestinal health and condition. Król et al. demonstrated that deregulation of the intestine (enteritis) induced the overexpression of genes linked to the oxidative defenses of epithelial cells [48]. A possible explanation for this larger number of differentially regulated genes may be that the combined hydrolysates provided a wider variety of functional peptides, an important characteristic linked to the performance of dietary protein hydrolysates [34].
Nevertheless, a focus could be performed further by qPCR in order to target the major actors of the main metabolic pathways highlighted by this study.
Restoring the intestinal transcriptomic profile of the high FM diet with a low FM diet supplemented with protein hydrolysates Hydrolysate diets induced a similar pattern of gene expression as the control diet (FM20), regardless of the included hydrolysate. Compared to the FM20 diet, the number of differentially regulated genes was very low for all hydrolysate diets, with 6, 8, and 6 genes for the TH, SH, and MH diets, respectively (Additional file 3). Once again, regulated genes were specific to hydrolysate origin, so that no gene was regulated by all three TH, SH and MH diets as compared to the FM20 group. Another mechanism involved in the differences in gene regulation between hydrolysate diets may be related to the glycosylation types in the resulting peptide fragments. The method used to generate these hydrolysates may also influence the degree of conservation of the glycosidic residues. Immune stimulation was more prominent in sea bass fed with SH; in parallel, glycosylation patterns conserved among plants and invertebrates have been found to be immunogenic in vertebrates [92]. Further work in the area of protein glycosylation may also be profitable and provide insights into this hypothesis. Although it is difficult to conclude about the relationship between regulated genes and the functional benefits of the TH diet, our results show that the SH diet down-regulated members of the GTPase IMAP family (members 6, 7 and 8) related to immunity. This highlights again the immunomodulatory effect of SH. Regarding the MH diet, Acyl-CoA desaturase and UDP-glucose 6-dehydrogenase were respectively overand under-expressed as compared to the FM20 diet. Therefore the lipid and carbohydrate pathways were modulated by this dietary treatment. Achieving the performance of a diet rich in FM during grow-out is a key point when dealing with low FM diets and alternative ingredients. In this study, the incorporation of 5% of hydrolysate compensated for the metabolic disturbances associated with the replacement of 15% of FM by a mixture of plant proteins, as observed in European seabass fed the FM5 diet. The intestinal transcriptomic response confirmed the results obtained at the zootechnical and gut histological levels, i.e. protein hydrolysates of aquaculture origin are valuable candidates to support FM replacement without deteriorating the zootechnical and functional performances of aquafeeds.

Conclusion
Within the global context of FM replacement in aquaculture feed formulation, evaluating the performance of new raw materials has become a crucial economic stake for aquafeed manufacturers. Although assessing the zootechnical performances of new formulations remains indispensable, the development of new protocols and tools for molecular and biochemical investigations can provide a more comprehensive view of the response of animals to novel diets. We combined zootechnical and transcriptomic approaches to evaluate the performances of protein hydrolysates as surrogates of dietary FM in European seabass. The reduction of the dietary FM level from 20 to 5% significantly impaired growth performances, intestinal histological organization, and induced significant changes in the transcriptomic profile of the intestine. By incorporating protein hydrolysate into the low FM diet (FM5) restored European seabass performances to a similar high level to those recorded with the high FM diet (FM20). In addition, this study demonstrates that the raw material used to manufacture protein hydrolysate is an important determinant in the transcriptomic response of fish.
Although the hydrolysates promoted the same zootechnical performances, the transcriptomic response of the intestine showed that genes and metabolic pathways were modulated in a hydrolysate-dependent manner. Moreover, including a combination of tilapia and shrimp hydrolysates into European seabass feed regulated a higher number of genes than the independent hydrolysates. It would be interesting to assess hydrolysate performances in fish grown under more challenging conditions (high density, low oxygen, low/ high temperature,…) to check how these differences in transcriptomic profiles are translated into phenotypic responses. These results offer an interesting scenario for the formulation of new high-performance feeds with a low level of FM inclusion. The analysis of the transcriptomic responses of liver and kidney sampled in the present study will allow us to complete this work and to better understand how metabolic pathways are modulated by the replacement of FM by protein hydrolysates in European seabass.

Availability of data and materials
The datasets supporting the conclusions of this article are included within the article (and its additional files). The raw data generated and/or analyzed during the current study is not publicly available due to the data mining still in progress. Nevertheless, some data could be available from the corresponding author on reasonable request.