# Integrative Metabolomics for Assessing the Effect of Insect (Hermetia illucens) Protein Extract on Rainbow Trout Metabolism

Simon Roques, Catherine Deborde, Laurence Guimas, Yann Marchand, Nadege Richard, Daniel Jacob, Sandrine Skiba-Cassy, Annick Moing, Benoit Fauconneau

► To cite this version:

Simon Roques, Catherine Deborde, Laurence Guimas, Yann Marchand, Nadege Richard, et al. Integrative Metabolomics for Assessing the Effect of Insect (Hermetia illucens) Protein Extract on Rainbow Trout Metabolism. Metabolites, 2020, 10 (3), pp.83. 10.3390/metabo10030083. hal-02570723

HAL Id: hal-02570723

https://hal.inrae.fr/hal-02570723

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Integrative Metabolomics for Assessing the Effect of Insect (*Hermetia illucens*) Protein Extract on Rainbow Trout Metabolism

Simon Roques<sup>1,2,3</sup>, Catherine Deborde<sup>3,4</sup>, Laurence Guimas<sup>5</sup>, Yann Marchand<sup>6</sup>, Nadège Richard<sup>2</sup>, Daniel Jacob<sup>3,4</sup>, Sandrine Skiba-Cassy<sup>1</sup>, Annick Moing<sup>3,4,*</sup> and Benoit Fauconneau<sup>1</sup>

Received: 3 February 2020; Accepted: 25 February 2020; Published: 27 February 2020

**Abstract:** Nutrition of high trophic species in aquaculture is faced with the development of sustainable plant-based diets. Insects seem particularly promising for supplementing plant-based diets. However, the complex effect of whole insect meal on fish metabolism is not well understood, and even less is known about insect meal extracts. The purpose of this work was to decipher the metabolic utilization of a plant-based diet supplemented with the gradual addition of an insect protein extract (insect hydrolysate at 0%, 5%, 10% and 15%). ¹H-NMR profiling was used to assess metabolites in experimental diets and in fish plasma, liver and muscle. A significant dose-dependent increase in growth and feed efficiency with increasing insect extract amounts was observed. The incremental incorporation of the insect extract in diet had a significant and progressive impact on the profile of dietary soluble compounds and trout metabolome. The metabolites modulated by dietary insect extracts in plasma and tissues were involved in protein and energy metabolism. This was associated with the efficient metabolic use of dietary free amino acids toward protein synthesis through the concomitant supply of balanced free amino acids and energy substrates in muscle. The findings provide new insights into how the dietary food metabolome affects fish metabolism.

**Keywords:** integrative metabolomics; fish nutrition; insect; proton NMR; profiling *Oncorhynchus mykiss*, *Hermetia illucens*

## 1. Introduction

Fish meal (FM) will no longer be sustainable for aquaculture in the foreseeable future because it is produced from limited stocks of pelagic fish, while the amount of feed used for farmed fish continues to increase [1]. Despite the widespread replacement of FM by plant-based ingredients, full plant-based feeds for high-trophic species still induce negative effects on fish growth and metabolism [2]. Thus, present-day commercial feeds always include FM. However, alternatives to marine resources exist to supplement plant-based feeds devoid of FM, and especially for protein ingredients

<footnote n="1">¹ INRAE, University Pau & Pays Adour, E2S UPPA, UMR NuMeA Nutrition, Metabolism and Aquaculture, F-64310 Saint Pée sur Nivelle, France; simon.roques@inrae.fr (S.R.); sandrine.skiba@inrae.fr (S.S.-C); benoit.fauconneau@inrae.fr (B.F.)</footnote>

<footnote n="2">² Phileo by Lesaffre, F-59700 Marcq-en-Baroeul, France; n.richard@phileo.lesaffre.com</footnote>

<footnote n="3">³ PMB-Metabolome, INRAE, 2018 Bordeaux Metabolome Facility, MetaboHUB, F-33140 Villenave d'Ornon, France; catherine.deborde@inrae.fr (C.D.); daniel.jacob@inrae.fr (D.J.)</footnote>

<footnote n="4">⁴ INRAE, University Bordeaux, UMR Fruit Biology and Pathology, F-33140 Villenave d'Ornon, France</footnote>

<footnote n="5">5 Copalis Industrie, F-62480 Le Portel, France; l.guimas@copalis.fr</footnote>

<footnote n="6">⁶ Le Gouessant, F-22402 Lamballe, France; yann.marchand@legouessant.fr</footnote>

<footnote>* Correspondence: annick.moing@inrae.fr; Tel.: +33-5-57-12-25-28</footnote>
with the use of single-cell proteins or land-animal proteins [3]. Amongst the sources of land-animal proteins, selected insect species can be incorporated into novel aquaculture diets, the most known being yellow mealworm beetle (Tenebrio molitor) and black soldier fly (Hermetia illucens). The addition of *H. illucens* larvae in fish feed composition already gives encouraging results on growth performance and is a good nutritional alternative to FM [4–7]. It has a high protein content, ranging from 40% to 63% of the total biomass depending on the rearing substrate [8], and an adequate amino acid profile quite similar to that of FM, especially regarding sulfur amino acids and lysine, which are limiting in most plant protein ingredients [9]. Protein extracts and hydrolysates are useful dietary ingredients for fish because they supply free compounds that can act as feed attractant, balance the absorption kinetic of amino acids in complete feed, and may contain bioactive peptides [10–13]. Furthermore, teleosts can specifically benefit of the direct absorption of small peptides in the larval stage and in the later stages, especially through the PepT1 transporter, which is activated during refeeding after a fasting state to ensure compensatory growth as compared to mammals and birds [14,15].

To assess a new ingredient in a feed, key indicators are the palatability of the whole feed, its digestibility and its utilization (growth, feed conversion, nutrient retention, etc.) and even somatic indexes [16]. These classic zootechnical indicators may be used to assess feed efficiency and growth performances but do not give any information about how feed influences these parameters. Omics approaches have led to deeper understanding of the effects of FM substitution in feeds on fish organisms. Nutrigenomics studies [17,18] have revealed several hepatic genes whose transcription is modified by a full plant-based diet in high-trophic species such as rainbow trout (*Oncorhynchus mykiss*) and sea bass (*Dicentrarchus labrax*). Similarly, metabolomic studies have been conducted on high-trophic species to assess the effect of full plant-based diets [19] or of diets based on alternative ingredients such as land proteins or single-cell proteins [20,21]. The combination of classical and omics approaches seems promising to decipher the metabolic mechanisms involved in phenotypic responses.

An integrative approach in nutritional studies should include multiple sources of data at different levels, e.g., feed characterization and fish phenotype, microbiome, transcriptome, proteome and metabolome, to cover the wide range of effects induced by a dietary treatment. By doing so, it could be possible to solve the complexity of effects induced by alternative ingredients used to complement plant-based diets. For example, multiple omic approaches in fish nutritional studies have been combined to assess the relationships between the gut microbiome and the metabolome in rainbow trout as well as in different coastal and estuarine species [22,23], or to link the transcriptome with the metabolome to assess the role of circadian rhythm on the nutritional state of leopard coral grouper (*Plectropomus leopardus*) [24]. However, studies on both feed compounds and fish metabolism by omic approaches are rare. Moreover, we have previously shown that numerous soluble compounds are present in fish feeds and that their characterization could be included in the nutritional integrative approach [25]. For example, Jasour et al. (2018) [26] demonstrated that biogenic amines in fish meal impact the rainbow trout metabolome.

Consequently, we investigated how a plant-based feed supplemented with graded amounts of insect protein extract (*H. illucens*) might affect the growth performance of rainbow trout and its metabolism. Metabolome was characterized in plasma, liver and muscle, using an integrative ¹H-NMR metabolomic approach. To assess whether metabolic effects and growth performance could be related with the soluble compounds present in the feeds, the main trends of the soluble compounds in the feeds were related to those of the metabolite profiles of fish depending on the amount (0%, 5%, 10% and 15%) of insect extract in the diets. We propose a synthetic vision of the metabolism modifications corresponding to the changes observed in the diets.
## 2. Results

### 2.1. Fish Growth

Initial and final body weights, feed conversion ratio (FCR), daily feed intake, hepatosomatic and viscerosomatic indexes and nitrogen retention are presented in Table 1. Supplementation of the plant-based (PB) diet with the insect extract had a positive and significant effect on fish final body weight (ANOVA, *P*-value = 4.10⁻⁴) with a progressive increase in final body weight depending on the amount of insect extract (INS). The final body weight of fish fed the diet with 5% of insect extract (INS05) was not significantly different from that of the PB diet and the same was observed between fish fed the INS10 and INS15 diets. The final body weight of fish fed the INS10 or INS15 diet was significantly different from that of fish fed the PB and INS05 diets (Tukey's test, *P* < 0.05). Daily feed intake was significantly modified by diet (ANOVA, *P* = 0.0321): with a lower daily feed intake in fish fed diets supplemented with the insect ingredient (Tukey's test, *P* < 0.05). Similar effects of the diets were observed on FCR (ANOVA, *P* = 0.0012) i.e., a progressive decrease from PB to INS15 diet with a significant difference between PB-INS05 diets, on one hand, and INS10-INS15, on the other (Tukey's test, *P* < 0.05). The nitrogen gain was similar between the four diets, but the nitrogen intake decreased when fish were fed insect-supplemented diets (Table S1). Thus, a significant effect (ANOVA, *P* = 0.003) of diet was observed on nitrogen efficiency ratio associated with a significant increase (Tukey's test, *P* < 0.05) for INS10 and INS15 diets compared to PB diet. Diet had no significant effect on hepatosomatic and viscerosomatic indexes (Kruskal-Wallis, *P* = 0.27 and *P* = 0.18, respectively).

<table><caption>**Table 1.** Growth and feed intake data for trout fed plant-based (PB) diet or diets with increasing amounts of insect ingredient (INS05, INS10 and INS15) for 84 days. Mean ± standard deviation. n = 3 tanks with 30 fish in each tank, or n = 9 fish for each diet.</caption><thead><tr><th>Zootechnical Parameters</th><th>PB</th><th>INS05</th><th>INS10</th><th>INS15</th></tr></thead><tbody><tr><td>Initial body weight (g) †</td><td>48.67 ± 1.35</td><td>48.90 ± 1.01</td><td>48.90 ± 1.01</td><td>48.90 ± 1.01</td></tr><tr><td>Final body weight (g) †</td><td>247.33 ± 13.99 <sup>a*</sup></td><td>258.67 ± 11.57 <sup>a</sup></td><td>293.56 ± 5.35 <sup>b</sup></td><td>303.33 ± 8.33 <sup>b</sup></td></tr><tr><td>Daily feed intake (% of body weight per day) †</td><td>1.57 ± 0.03 <sup>a</sup></td><td>1.50 ± 0.03 <sup>ab</sup></td><td>1.51 ± 0.03 <sup>ab</sup></td><td>1.48 ± 0.03 <sup>b</sup></td></tr><tr><td>Feed conversion ratio†</td><td>0.98 ± 0.03 <sup>a</sup></td><td>0.92 ± 0.03 <sup>ab</sup></td><td>0.89 ± 0.02 <sup>bc</sup></td><td>0.86 ± 0.02 <sup>c</sup></td></tr><tr><td>Nitrogen efficiency ratio (N retention % of N intake) †</td><td>40.32 ± 1.55 <sup>a</sup></td><td>44.68 ± 1.58 <sup>ab</sup></td><td>46.47 ± 2.48 <sup>b</sup></td><td>48.66 ± 1.06 <sup>b</sup></td></tr><tr><td>Hepatosomatic index ‡</td><td>0.77 ± 0.12</td><td>0.87 ± 0.11</td><td>0.84 ± 0.19</td><td>0.84 ± 0.9</td></tr><tr><td>Viscerosomatic index ‡</td><td>8.96 ± 1.37</td><td>9.77 ± 1.34</td><td>9.09 ± 0.58</td><td>8.74 ± 0.85</td></tr></tbody></table>

† data generated from pooled fish (n = 3 tanks of 30 fish each); ‡ data generated at end of trial from fish used for the metabolomics approach (n = 9 individual fish); * For final body weight, daily feed intake, feed conversion ratio (FCR) and nitrogen efficiency ratio, means accompanied by the same letter are not significantly different (Tukey's test, *P* > 0.05).

### 2.2. Annotation and Integration of Feed Soluble Compounds and Fish Metabolome NMR Spectra

In ¹H-NMR spectra of feed extracts, 28 compounds were putatively annotated (Table S2). In spectra of plasma and tissue extracts, a total of 84 metabolites were annotated: 23 in plasma, 36 in liver and 25 in muscle (Table S2). Twenty common metabolites were found in plasma and liver, 16 in plasma and muscle and 24 in liver and muscle. A Venn diagram details the common and specific metabolites within plasma and tissues (Figure S1). The 16 metabolites common to plasma, liver and muscle comprised 13 amino acids or their derivatives, two organic acids and one sugar.

### 2.3. Patterns of Compounds in Feed

To assess how compound relative intensities changed in feed with the inclusion of the insect ingredient, a hierarchical clustering analysis (HCA) was performed on compound intensity means
per diet. On the double HCA heatmap (Figure 1), the diets were ordinated from PB to INS15 and the feed compounds were grouped into two major clusters. The intensity changes of the compounds of each cluster are plotted in Figure S2. In Cluster 1 comprising eight compounds, the mean compound intensities decreased linearly from PB to INS15 feed. In Cluster 2 comprising 20 compounds, the mean compound intensities increased from PB to INS15 feed but less linearly than in Cluster 1.

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Hierarchical clustering analysis (HCA) heatmap of mean relative intensities of 28 compounds quantified using ¹H-NMR profiling in feed for control diet (PB) and three diets with increased amounts of insect ingredient (INS05, INS10, INS15).</figcaption></figure>

## 2.4. Fish Separation Based on Insect Amount in Diet

Based on the Kruskal–Wallis analysis of the intensities of 84 metabolites from plasma, liver or muscle NMR spectra, 32 metabolites were significantly affected by diet, 7 in plasma, 16 in liver and 9 in muscle (Figure S3). Thus, we analyzed the metabolome for specific changes in metabolites that could be related to linear trends observed for dietary compounds. We used a linear discriminant analysis to identify which metabolite intensities tended to increase or decrease linearly in the fish metabolome. A Partial Least Squares Discriminant Analysis with an Orthogonal Signal Correction (OSC-PLS-DA, Figure 2) was used to discriminate fish based on their metabolomes of plasma and tissues according to the four diets and to identify which metabolite variables were most representative of the fish separation between the four diets. The predictive accuracy and percentage of the variance explained by the diet factor for the OSC-PLS-DA were high: Q² = 0.90 and R²Y = 90.1%. The cross-validation *P*-value was significant (*P* = 0.036). As expected, the scores plot (Figure 2a) showed a clear group separation along the first component that concorded with the amount of insect protein extract. Despite the orthogonal signal correction, the diet groups were slightly tilted from a vertical axis. When integrating the scores (Figure 2a) and the loading plots (Figure 2b), overall plasma and muscle metabolites tended to have higher intensities in PB-fed fish, and liver metabolites tended to have higher intensities in INS15-fed fish.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** OSC-PLS-DA of 84 metabolites measured in fish plasma, liver or muscle, for fish fed PB, INS05, INS10 or INS15 diet for 84 days. Q² = 0.90, R²Y = 90.13%, P = 0.036. (a) Scores plot of individuals on first two components. Green circles, PB; light-blue triangles, INS05; blue diamond, INS10; dark-blue squares, INS15 diet. (b) Loading plots on first two components. Squares, plasma metabolites (P_); triangles, muscle metabolites (M_); diamonds, liver metabolites (L_). All metabolites with a loading value below -0.06 or above 0.06 on one of the first two components are annotated and colored according to their compound family: red, amino acids or dipeptides; green, organic acids; blue, sugars; black, alcohols or cyclitols; pink, amines or N-containing compounds; yellow, nucleobases or nucleosides; purple, vitamins.</figcaption></figure>
A threshold of 1.5 was chosen for the variable importance in the projection (VIP) scores to select the metabolite intensities that contributed most to fish separation according to the diet. Thirteen metabolite intensities had a VIP score higher than 1.5 (Figure 3). Among them, two were from plasma (histidine, alanine), three from muscle (taurine, histidine, acetate), and eight from liver (alanine, vitamin B5, beta-alanine, taurine, 1,2-propanediol, glutamine, choline, glutamate). All these metabolite intensities were significantly affected by the diet (Kruskal-Wallis analysis, *P* < 0.05), with the lowest *P*-value of 0.003 for plasma alanine, except liver choline intensity (*P* = 0.11). Taurine contributed to discriminating diets both in muscle and liver tissues, histidine both in plasma and muscle tissue, and alanine both in plasma and liver tissue. Glutamine and glutamate in liver tissue tended to discriminating diets.

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Variable importance in the projection (VIP) scores from the OSC-PLS-DA of 84 metabolites measured in fish plasma (P_), liver (L_) or muscle (M_), for fish fed PB, INS05, INS10 or INS15 diet for 84 days presented Figure 2. The scores values of the 13 variables with a score higher than 1.5 are presented.</figcaption></figure>

Noticeably, it was possible to plot and classify these metabolites for the four diets (Figure 4) according to two patterns paralleling the patterns highlighted for the feed compounds: gradual decrease (Figure 4a) or increase (Figure 4b) from PB to INS15. Intensities of liver taurine, muscle histidine, muscle acetate and plasma histidine were significantly reduced in INS-fed groups compared to the PB-fed group. Muscle taurine and liver alanine intensities were significantly lower in INS10 or INS15 groups than in PB and INS05 groups. Liver glutamate and liver 1,2-propanediol intensities were significantly higher in INS groups than in the PB group. Liver vitamin B5 and liver beta-alanine intensities were significantly higher in INS10 and INS15 groups than in PB and INS05 ones.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Whisker plots (n = 9 or 8) for 12 fish variables selected using OSC-PLS-DA presented in Figure 2 based on VIP score higher than 1.5 (Figure 3) for diet separation, and significantly affected by diet based on Kruskal-Wallis analysis (*P* < 0.05). (a) Fish variables showing a decreasing trend with level of insect ingredient amount in diet. (b) Fish variables showing an increasing trend with level of insect ingredient amount in diet. For each plot, medians of diets were compared according to Dunn's test (*P* < 0.05): diets accompanied by the same letters (a, b or c) above their corresponding whisker plots were not significantly different according to Dunn's test.</figcaption></figure>

Diet had a significant effect on the mean of the sum of six of the free essential amino acids in fish muscle (ANOVA, *P* = 5.4.10⁻⁵) measured using the intensities of the identified corresponding amino acids on the zgpr NMR spectra (Figure S4). The intensity sum of all these free essential amino acids identified and quantified were significantly lower in INS15 and INS10 fed groups compared to PB, while that of INS05 did not differ from PB (Tukey's test, *P* < 0.05).

**3. Discussion**

This study was conducted to assess the effect of supplementing an insect protein extract to a plant-based diet on the metabolome of rainbow trout. The hypothesis was that the numerous soluble compounds of the insect extract ingredient could act either as attractant to stimulate feed intake or metabolically to improve fish growth performances fed plant-based diets.
Based on the growth performances (Table 1), the insect protein-extract ingredient significantly increased final body weight and improved feed efficiency and nitrogen retention efficiency. However, the daily feed intake did not increase but significantly decreased from PB to INS15, demonstrating that the soluble compounds of the insect ingredient in the feed had no effect on the attractiveness of the diet. Other factors such as energy content or nutrient availability in feed usually drive feed intake [27]. The experimental diets were iso-energetic while free amino acid contents were higher in INS diets. Thus, it is suggested that the supply of nutrients as free form from INS diets covered the nutrient requirements faster than with PB diets. This contributed to the increased protein retention efficiency and probably explained the decreased feed intake also. Therefore, this offered a unique opportunity to analyze the changes induced by insect diets on metabolite patterns by assessing the potential metabolic effect of soluble compounds found in the diet.

### 3.1. Protein Synthesis-Oriented Metabolism Induced by Dietary Compounds

An interesting finding is the clear amino acid and protein metabolic status in trout, which could be deduced from the metabolite pattern of glutamate, glutamine in liver and alanine in liver and plasma combined with the decrease in free essential amino acids content in muscle. In carnivorous fish, amino acids not only constitute building blocks for protein synthesis in tissues but are also used as preferential energetic substrates together with fatty acids [28]. As energetic substrate, amino acids are catabolized by transamination in liver, except for the branched-chain amino acids which are mostly catabolized in muscle. Transamination transfers the amine function from an amino acid to an α-ketonic acid that generates glutamate and a novel α-ketonic acid. In the present study, the glutamate intensities in liver increased consistently with the ratio of insect extract in feed, which suggests a predominant amino acid catabolic activity for higher energy supply. Glutamate is then converted back to α-ketoglutarate by glutamate dehydrogenase, a step that produces free ammonia. However, even though teleosts can process and excrete high levels of ammonia directly, too high a concentration results in harmful effects [29]. Thus, ammonium can be fixed on glutamate to form glutamine through glutamine synthetase and alleviate such effects [30]. The concomitant rise in glutamine intensities observed in our study in the liver of fish fed the diet containing a high level of insect ingredient is consistent with this effect. In this case, increased glutamate and glutamine intensities in liver indicate an increase in amino acid catabolism for energy supply (Figure 5). This process might be fueled by the continuous supply of both free amino acids from the insect extract and the delayed release of amino acids from protein digestion. As we did not identify any changes in aspartate or asparagine intensity, the glutamate may be a privileged pathway to fuel the TCA cycle.

While there was a relative hepatic increase in glutamate and glutamine intensities, a decrease in alanine intensity was observed in liver and plasma. Alanine may be considered as one of the end-products of amino acid catabolism in muscle. The amino acids in excess in muscle during fasting or active swimming activity are further transaminated on pyruvate to produce alanine. Alanine is then transported through the blood to the liver, thus contributing to the Cahill cycle between muscle and liver (Figure 5). In fish fed insect extracts from 0% to 15% (PB to INS15), a decrease in alanine intensity was observed both in plasma and liver, which suggests that amino acid catabolism is reduced in the muscle of fish fed insect extract-supplemented diets.

The muscle metabolism seems to be preferentially oriented toward protein synthesis in fish fed diets supplemented with insect extract, as suggested by the reduced accumulation of free essential amino acids in muscle (Figure S4) and the increased whole-body nitrogen retention efficiency (Table 1). This protein metabolism state might explain the reduced levels of alanine observed in liver and plasma, which serve as a marker of reduced amino acid catabolism. This metabolism oriented towards protein synthesis in muscle might be favored by the significant supply of free amino acids and peptides from insect protein extracts, which are rapidly absorbed and transported to muscle compared to plant-based diets, in which the absorption of amino acids is generally delayed [31]. Furthermore, the increased glutamate and glutamine levels in liver indicate the increased use of amino acids to supply the energy required for protein synthesis in muscle.
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Simplified diagram of the relationships between liver, blood and muscle metabolites involved in Cahill cycle and protein synthesis, integrating hypotheses based on the present experiment. The increase of free amino acids from INS diet would rise the glutamate and glutamine contents in liver to supply energy to muscle. In muscle, the free amino acids are rather used for protein synthesis than for energetic purpose. The alanine content from Cahill cycle would consequently decrease.</figcaption></figure>

## 3.2. Metabolites as Markers of Alteration of Metabolism Induced by PB Diet

The decrease in taurine intensity in muscle and liver of fish fed the INS diet could partly be due to the profile of the soluble compounds of insect feed. Taurine is a precursor of biliary acids for lipid emulsification which facilitate their digestion. It is also involved in detoxification processes in liver and in osmoregulation, especially in muscle and heart [32]. Although rainbow trout can synthesize taurine *de novo*, it is considered as a semi-essential nutrient in rainbow trout, because a dietary taurine deficiency may degrade growth performance of teleosts [32]. It has been shown that trout hepatic and muscle taurine are a good marker of the level of taurine in the diet [33]. However, plant ingredients are not known to contain taurine, unlike animal feed ingredients [34]. Although insects may synthesize taurine, black soldier fly larvae meal contains very low amounts of taurine compared to fish meal [35]. Furthermore, taurine was not supplemented in our diets (Table 2), and no differences in taurine were observed in the feed profiling of soluble compounds. Thus, taurine in fish tissue must come exclusively from *de novo* synthesis. In rainbow trout, the main pathway for taurine synthesis is not yet clearly established, but most evidence indicates that synthesis goes through cysteine to form cysteine sulfinic acid, then hypotaurine, and finally taurine [32]. Cysteine is a non-essential sulfur-containing amino acid supplied either directly by ingredients in the diet or synthesized from methionine. However, diet formulation was done to ensure that total methionine and also cysteine contents were not different in the experimental diets. Thus, the relative predominance of taurine in the tissues of fish fed PB feeds could be due to the stimulation of taurine synthesis to ensure its specific functional properties [32]. The anti-nutrient factors of plant ingredients such as sucrose and stachyose oligosaccharides detected in PB diets are known to be osmotically active in the fish gut [36] where they alter both gut structure and function. Consequently, ion absorption might be altered, thus triggering the accumulation of taurine in liver and muscle. However, further information on the alteration of gut ion absorption in fish fed plant-based diet is required to strengthen this hypothesis.

Histidine also decreased in plasma and muscle in response to insect protein extract supplementation. Gomez-Requeni et al. (2004) [37] also found that free histidine in muscle gradually decreased with the reduction of plant ingredients in feed for a carnivorous fish. Thus, this seems to confirm that histidine changes in fish metabolome are not related to a specific alternative ingredient but rather to the decrease in plant protein ingredients. As histidine is a compound that shares the same pH buffering and metal-chelation properties as its related peptides anserine and carnosine, we suggest that an imbalanced pH homeostasis induced by the PB diet was progressively changed by insect protein extract supplementation. The increasing levels of beta-alanine in liver as the amount of
insect ingredient increased suggests that the status of carnosine, another related peptide involved in buffering capacity, is also modified.

Amongst the metabolites that demonstrate significant changes with insect protein extract complementation, some compounds such as vitamin B5 and 1,2-propanediol are interesting markers of insect ingredient composition. The annotation of the NMR signal corresponding to vitamin B5 requires further investigation, however the progressive improvement in vitamin status may contribute to the improvement of nutritional status in fish fed the INS diet. Although 1,2-propanediol could be considered as a potential energy substrate in fish, it is difficult to establish from its accumulation in liver of fish fed insect extract whether it plays any functional or metabolic role in fish.

# 4. Materials and Methods

## 4.1. Growth Trial and Sampling

Three hundred sixty female juvenile rainbow trout (49.10 ± 1.05 g average body weight) from INRAE experimental fish farm facilities (INRAE, UMR Nutrition, Metabolism and Aquaculture, F-64490 Lees-Athas, F-40360 Donzacq, France) were reared in 100 L tanks of a flow-through system supplied with water at constant temperature 17 ± 1 °C and constant oxygen level. Fish were split randomly into 12 groups of 30 fish. The fish were fed for 84 days either a full plant-based feed (PB) or plant-based feeds supplemented with 5%, 10% and 15% of insect ingredient (INS05, INS10 and INS15, respectively), i.e., three replicate tanks per diet. Fish were hand-fed twice a day until visual satiation. For each group, the biomass and the amount of distributed feed were measured every three weeks and the mean final body weight was calculated by dividing the total biomass weight by the number of fish alive in each tank. Nine fish per diet (3 per tank) were sampled 48 h after the last meal. Fish were sedated by immersion in a tank containing 10 mg/L benzocaine solution, then anaesthetized by immersion in a benzocaine solution of 30 mg/L. Then blood samples were collected with heparinized syringes and poured into 2 mL Eppendorf tubes. Each tube was centrifuged (3000× g, 5 min) and 800 µL of plasma were collected in a 2 mL Eppendorf tube, immediately stored at -20°C until the end of sampling, then further stored at -80 °C until sample preparation for NMR analysis. Just after the blood sampling, fish were euthanized by section of the spinal cord. The liver was collected and the gallbladder was delicately removed. The liver was rinsed in a saline solution (NaCl 9 g/L), dried on a wrapped paper then immediately weighed and frozen in liquid nitrogen. The whole digestive tract was removed manually and weighed. A portion of 2-3 g from the dorsal epaxial white muscle was sampled just below the dorsal fin after removal of the skin and the red muscle and frozen in liquid nitrogen. Liver and muscle samples were stored at -80 °C until sample preparation for NMR analysis. The ratio between liver mass to fish mass or gut mass to fish mass was calculated.

Both growth trial and sampling were in strict accordance with (i) EU legal frameworks concerning the protection of animals used for scientific research (Directive 2010/63/EU), (ii) the National Guidelines for Animal Care of the French Ministry of Research (decree n°2013-118, February 1st, 2013), and (iii) the local ethical committee "Comité d'Ethique en Expérimentation Animale Aquitaine Poissons Oiseaux" (C2EA-73). The trial did not need specific ethical approval because it involved standard rearing practices and diets formulated to cover the nutritional requirements of rainbow trout. All the staff of the experimental facility received training and personal authorization (N°B64 10 005).

## 4.2. *Diet and Fish Whole Body Composition Analysis.*

The feed preparation was described in Roques et al. (2018) [25]. All diets were isoproteic, isolipidic and isoenergetic. The control diet was exclusively based on plant ingredients (PB) so it was devoid of FM and fish oil, and was supplemented with graded amounts (0%, 5%, 10% and 15%) of an insect protein hydrolysate ingredient from defatted larvae of black soldier fly (*H. illucens*) (PB, INS05, INS10, INS15, respectively) (Table 2) described previously [25].
<table><caption>**Table 2.** Diet formulation (g.100 g⁻¹) of control diet (PB) and experimental diets (INS05, INS10, INS15). Ingredients and gross energy.</caption><thead><tr><th>Ingredients</th><th>PB</th><th>INS05</th><th>INS10</th><th>INS15</th></tr></thead><tbody><tr><td>DHA*-rich algae meal</td><td>6.84</td><td>6.84</td><td>6.84</td><td>6.84</td></tr><tr><td>Insect protein hydrolysate</td><td></td><td>5.00</td><td>10.00</td><td>15.00</td></tr><tr><td>Vegetable oils <sup>(a)</sup></td><td>18.10</td><td>16.95</td><td>16.00</td><td>14.95</td></tr><tr><td>Plant proteins <sup>(b)</sup></td><td>70.40</td><td>66.74</td><td>62.78</td><td>58.79</td></tr><tr><td>Rapeseed lecithin</td><td>1.00</td><td>1.00</td><td>1.00</td><td>1.00</td></tr><tr><td>Monocalcium phosphate</td><td>1.20</td><td>1.10</td><td>1.00</td><td>1.00</td></tr><tr><td>Phytase</td><td>0.02</td><td>0.02</td><td>0.02</td><td>0.02</td></tr><tr><td>Lysine 78%</td><td>0.50</td><td>0.50</td><td>0.50</td><td>0.50</td></tr><tr><td>DL-methionine 98%</td><td>0.65</td><td>0.56</td><td>0.57</td><td>0.61</td></tr><tr><td>Threonine 98%</td><td>0.20</td><td>0.20</td><td>0.20</td><td>0.20</td></tr><tr><td>Vitamin premix</td><td>0.30</td><td>0.30</td><td>0.30</td><td>0.30</td></tr><tr><td>Vitamin C monophosphate 35</td><td>0.04</td><td>0.04</td><td>0.04</td><td>0.04</td></tr><tr><td>Mineral premix</td><td>0.30</td><td>0.30</td><td>0.30</td><td>0.30</td></tr><tr><td>Liquid choline</td><td>0.15</td><td>0.15</td><td>0.15</td><td>0.15</td></tr><tr><td>Antioxidant</td><td>0.15</td><td>0.15</td><td>0.15</td><td>0.15</td></tr><tr><td>Antifungal</td><td>0.15</td><td>0.15</td><td>0.15</td><td>0.15</td></tr><tr><td><strong>Energy (kJ.g<sup>-1</sup> dry matter)</strong></td><td>24.03</td><td>24.50</td><td>24.59</td><td>24.71</td></tr></tbody></table>

*DHA, docosahexaenoic acid; (a) Vegetal oils including rapeseed oil and linseed oil; (b) Plant proteins including wheat gluten, hydrolyzed wheat gluten, pea protein, faba bean protein concentrate, soy concentrate, soybean meal, rapeseed meal, peeled faba bean and wheat.

Proximate analysis of the experimental diets and measurement of whole-body composition were performed as follows. Dry matter was analyzed by drying the samples to constant weight at 105 °C for 24 h. Crude protein was determined using the Kjeldahl method after acid digestion and estimated by multiplying total nitrogen by 6.25. Crude lipid content was quantified by petroleum diethyl ether extraction using the Soxhlet method. Gross energy content was determined in an adiabatic bomb calorimeter (IKA). Ash was examined by combustion in a muffle furnace at 550 °C for 16 h.

Daily feed intake as percentage of body weight was calculated as follows:

$$ \frac{F_d \times 100}{\left( \frac{BW_i + BW_f}{2} \right) \times n_j} $$

Where $F_d$ is the feed consumed during the growth period, $BW_i$ and $BW_f$, the initial and final average body weight, respectively and $n_j$ the number of days of trial. The feed conversion ratio, FCR, represents the efficiency of the transformation of feed to biomass and was calculated as follows:

$$ \frac{F_d}{BW_f + BW_d - BW_i} $$

Where $BW_d$ is the total biomass of the dead fish.

Nitrogen intake was calculated as follows: feed intake x nitrogen content in diet. Nitrogen gain corresponds to the difference of nitrogen in fish (fish body weight x nitrogen content in fish) between the end of the experiment and the beginning divided by the duration of experiment in days. The nitrogen efficiency ratio is the ratio between nitrogen gain in fish and nitrogen intake from feed during the experiment.

### 4.3. Sample preparation

For feed samples, freeze-dried feed powders were extracted according to Roques et al. (2018) [25] in five replicates per feed. Briefly, a series of hot ethanol/water extractions was performed, followed by steps of lyophilization and solubilization with deuterated phosphate buffer solution and pH adjustment to apparent pH 6.00 ± 0.02 by means of a titration robot (BTpH, Bruker, Kalrsruhe,
Germany). Lyophilized pH-adjusted feed extracts were dissolved again in deuterated water containing deuterated ethylene tetra-acetic acid disodium salt (EDTA-d₁₂) to chelate paramagnetic cations and (trimethylsilyl)propionic-2,2,3,3-d₄ acid sodium salt (TSP) was used for chemical shift calibration just before NMR spectra acquisition.

All fish were still juvenile female, except one male that was discarded from sample preparation. Plasma samples were individually thawed on an ice bed for a minimum of 2 h and no more than 2 h and 20 min, then 300 µL were transferred to an Eppendorf tube and diluted with 300 µL D₂O (99.8%, Eurisotop, Gif sur Yvette, France). A simple dilution with a lock substance is suitable in liquid ¹H-NMR to efficiently assess metabolites as a specific acquisition sequence (see Section 4.5) can be used to reduce the signals from macromolecules [22]. The tube was vortexed for 5 s maximum and 550 µL were transferred to a 5-mm NMR tube (Wilmad, Vineland, USA). No TSP was added to the plasma sample because proteins of the native plasma could bind to TSP and shift the -CH₃ resonance from 0 ppm. NMR acquisition was performed immediately after each sample preparation.

For liver and muscle samples, polar metabolites were extracted using phase separation based on methanol/chloroform/water recommended by Lin et al. (2007) [38]. In our case however for safety reasons the chloroform was substituted by dichloromethane which share similar extraction efficiency. The extraction solvent was a mix of dichloromethane (Honeywell Riedel-de Haën, Puriss ACS, Charlotte, North Carolina, USA) and methanol (Carlo Erba, RPE ACS, Chaussée du Vexin, France) (2:1 v:v) with 0.01% of butylhydroxytoluene as antioxidant (Sigma-Aldrich, Saint-Louis, Missouri, USA). A sample of 1.4 ± 0.2 g frozen tissue was weighed and immediately transferred to a 50-mL Falcon tube on an ice bed. Fifteen mL of the extraction solvent were poured into the 50-mL Falcon tube containing the frozen tissue. The tissue was ground with an Ultra-Turrax T25 disperser (IKA, Staufen, Germany) for 1 min at 9500 rpm, then the homogenate was filtered with two overlaid Whatman filter papers on the Büchner through vacuum filtration. The residual powder was resuspended in 15 mL of extraction solvent and the process repeated once. The combined filtered solutions were mixed with 7.5 mL of aqueous NaCl (0.73% w/v) to enhance the phase separation and were allowed to decant for 10 min minimum. They were then centrifuged at 2200× g for 10 min (AVANTI J-20 XP centrifuge, Beckman Coulter, Brea, California, USA). Fifteen mL of the polar phase were then collected and stored at -20°C until preparation for NMR acquisition.

To prevent a biased estimation of metabolite content, the volume of the polar phase collected was carefully adjusted to the weight of frozen tissue extracted in order to obtain the same final dried matter in each polar extract. The polar phase was transferred to a glass tube (Pyrex disposable culture tube, Mexico), then dried under a moderate nitrogen flow (ALPHAGAZ 1 Azote SMARTOP, Air Liquide, Paris, France) until visual dryness. Dry extracts were solubilized in 140 µL of a deuterated 200 mM potassium phosphate buffer solution containing 2 mM EDTA-d₁₂ (Sigma-Aldrich, Saint-Louis, Missouri, USA), vortexed for a few seconds, then transferred to a 2-mL plastic Eppendorf tube with a Pasteur pipette. Each glass tube was rinsed with 230 µL, then 330 µL D₂O. The buffered extract was vortexed for a few seconds, then pH was adjusted to apparent pH 6.00 with NaOD or DCl using a titration robot (BTpH, Bruker, Kalrsuhe, Germany). Finally, 550 µL of pH-adjusted extract were transferred to a 5 mm NMR tube (Wilmad, Vineland, USA) with 5 µL of TSP-d₄ (98% D) from Sigma-Aldrich (Paris, France).

### 4.4. Total Amino Acids in Feeds

Total amino acid profiles of feeds were determined by targeted ion-exchange HPLC-UV analyses of feed hydrolysates performed by Upscience (Vannes, France) according to the NF standards EN ISO 16634-1 and are presented as % of feed.

### 4.5. NMR Acquisition

NMR acquisition was performed in 5-mm NMR tubes using a Bruker Avance III 500 MHz spectrometer (Bruker, Wissembourg, France) equipped with a sample holder and an ATMA broadband inverse probe. 1D NMR acquisition parameters for feed extracts, have been detailed previously [25].
For plasma, 1D Carr-Purcell-Meiboom-Gill (CPMG) ¹H-NMR spectra with presaturation (cpmgpr1d) were acquired with 8 dummy scans before 64 scans of 32K data points, a 90° pulse angle, a 6000 Hz spectral width with 2.73 s acquisition time and 2 s recycle delay. A 200 µs echo time and 150 loops were used as parameters for the CPMG experiment. A 5.10⁻⁵ W irradiation power was used for presaturation.

For liver and muscle, 1D ¹H-NMR spectra of tissue extracts with presaturation (zgpr) were acquired with 32 scans of 64K data points, a 90° pulse angle, a 6000 Hz spectral width with 5.46 s acquisition time and 10 and 25 s recycle delay for liver and muscle, respectively. A 7.10⁻⁵ and 8.10⁻⁵ W irradiation power was used for presaturation for liver and muscle, respectively. Two-dimensional experiments, J-resolved (Jres), ¹H-¹H homonuclear correlated spectroscopy (COSY) and ¹H-¹³C heteronuclear single-quantum correlation (HSQC), were performed on selected samples to help for the annotation of ambiguous signals.

## 4.6. NMR Spectra Processing

Spectra were processed in groups of sample type with the NMRProcFlow tool [39]. Each spectrum was first zero-filled to four times their data points or twice for muscle. Free induction decay (FID) was then multiplied by an exponential window function with a line broadening of 0.3 Hz and automatically zero-order phased. For feed and tissue extracts, the spectra were calibrated on -CH₃ of TSP-d₄ (δ 0 ppm), and for plasma on -αCH of glucose (δ 5.22 ppm). A representative 1D ¹H-NMR spectrum of feed extracts, plasma, and liver or muscle extracts was annotated based on Roques et al. (2018) [25], Kullgren et al. (2010) [40], Shen et al. (2018) [41] and chemical shift, signal multiplicity and intensity ratio were compared to an in-house database, a public database (HMDB, http://www.hmdb.ca) and the ChenomX NMR Suite library 8.3 (ChenomX Inc., Edmonton, Canada). Metabolite annotation followed the recommended guidelines of Sumner et al. (2007) [42]. All full spectra of feed, plasma or tissues were manually corrected for baseline and uncontrolled peak shifts. A variable size spectral region, or bin, was defined on one well-resolved resonance for the relative quantification of each metabolite, on the full pattern whenever possible or on resolved peaks of the pattern, as detailed in Table S1. These spectral regions were then normalized to constant sum for each dataset. Each integrated spectral region was named after the compound/metabolite it corresponded to, preceded by the letter of the sample type (F: feed, P: plasma, L: liver, M: muscle).

## 4.7. Statistical Analysis

For the physiological parameters of fish, a one-way ANOVA was carried out with R statistical software (v.3.6.1, R Development Core team, 2008) to assess the effect of diet on fish final body weight BW, feed conversion, feed intake and nitrogen efficiency. When the *P*-value was lower than 0.05, Tukey's test was carried out for diet mean comparisons. Similarly, the Kruskal-Wallis test was used to assess the effects of hepatosomatic and viscerosomatic index on the fish sampled for metabolomics analysis (n = 9 per diet).

For feed characterization, before performing the hierarchical clustering analysis (HCA) of feed soluble compound ¹H-NMR data, the normalized values of bins were mean-centered and scaled to unit variance. The normalized and scaled values of all bins were used to calculate the relative mean intensity of compound per diet (n = 5 for experimental diet except for PB feed n = 4, where one sample was removed due to contamination detected during extraction). The HCA was then computed using R statistical software with Pearson correlation distance and average linkage.

For fish metabolomic data, prior to the Orthogonal Signal Correction Partial Least Squares Discriminant Analysis (OSC-PLS-DA), the data from plasma, liver and muscle were combined into one dataset. One sample of plasma was missing so the Non-linear Iterative Partial Least Squares (NIPALS) algorithm was used to impute the corresponding values. The bins were then mean-centered and scaled to unit variance. The imputation and OSC-PLS-DA were performed with BioStatFlow web software (v.2.9 http://www.biostatflow.org) on the combined dataset. A 200-fold cross validation (leave one out) was used to assess the significance (*P*-value) of the prediction and separation parameters (Q² and R²Y). A threshold of 1.5 was chosen for the selection of variable
importance in the projection (VIP). The diet effect was then tested on these variables with the Kruskal–Wallis test (*P* < 0.05) and Dunn's test (*P* < 0.05) for mean comparisons performed with R statistical software.

Based on the NMR annotation, we estimated the sum of free essential amino acids by calculating the sum of relative intensities of six of the free essential amino acids identified in the muscle for each fish to assess if fish tended to accumulate the free essential amino acids. The mean of the sum of free essential amino acids was calculated for each diet group (n = 9, except n = 8 for group INS15 where a male was discarded). A one-way ANOVA was used to determine whether diet had a significant effect on the mean of the sum of free essential amino acids in muscle. When the *P*-value was lower than 0.05, Tukey's test was performed for diet mean comparisons.

**5. Conclusions**

Integrative ¹H-NMR-metabolomics on both diet and fish have helped to identify the main underlying mechanisms contributing to fish performance. It thus appears that the insect (H. illucens) protein-hydrolysate tested in this work is a valuable alternative protein source to substitute plant- ingredients in fish diet for the improvement of fish performance in fish. Indeed, this ingredient increased the protein efficiency ratio without any increase in feed intake. The supply of higher level of free dietary amino acids and peptides by insect protein extract provides substrates for both protein synthesis-oriented metabolism and the associated energy requirement. This was strengthened by the glutamate and glutamine increase in liver, consistently with the reduced accumulation of essential free amino acids in muscle and reduced hepatic and circulating alanine. Likewise, dietary modulations could be also viewed through significant improvement of important metabolite status, as with taurine and histidine. Such an integrative approach improves the understanding of the metabolic pathways and processes involved in fish metabolism. It is a promising tool to improve new fish feeds formulated with alternative ingredients.

**Supplementary Materials:** The following are available online at www.mdpi.com/xxx/s1, Figure S1: Venn diagram of metabolites identified in plasma and tissues; Figure S2: Mean of all compounds of each major cluster of Figure 1 HCA for four diets.; Figure S3: Metabolites in plasma, liver or muscle for which diet had statistically significant effect; Figure S4: Mean of sum of free amino acids calculated from ¹H-NMR integration in muscle; Figure S5: Representative 500 MHz ¹H-NMR spectra of plasma of trout fed INS15 diets; Figure S6: Representative 500 MHz ¹H-NMR spectra of liver polar extract of trout fed INS15 diets; Figure S7: Representative 500 MHz ¹H-NMR spectra of muscle polar extract of trout fed INS15 diets; Table S1: Amount of feed consumed during fish growth period, nitrogen intake and nitrogen gain; Table S2: Annotation table of feed compounds and metabolites in plasma, liver and muscle; Table S3: MSI levels for annotated compounds and metabolites; Table S4: Total amino acid (free amino acids and protein amino acids) concentrations in the four feeds.

**Author Contributions:** S.R., B.F., and A.M. conceptualized and designed the experiments; Y.M. and L.G. supplied the resources for the experiment; S.R. and C.D. performed the experiment and validated the data; S.R., A.M. and D.J. performed the statistical analyses; S.R., B.F., and S.S.-C. interpreted the data; S.R. prepared and wrote the original draft; S.R., B.F., A.M., C.D., S.S.-C., and N.R. wrote, reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by FUI 2014 NINAqua, grant number 15008485-056 Nouvelle-Aquitaine Regional Council funding and MetaboHUB, grant number ANR-11-INBS-0010. Simon Roques was funded by "Phileo by Lesaffre" and the ANRT (Association Nationale Recherche Technologie).

**Acknowledgments:** We warmly thank Frédéric Terrier and his team from the INRAE experimental fish farm facilities at Donzacq (France) for managing the trial so well and helping us during sampling, Vincent Véron and Marianne Cluzeaud for their precious help during sampling, Laurence Larroquet for her help in sample preparation, and Ray Cooke for language editing.

**Conflicts of Interest:** S.R. was financed by "Phileo by Lesaffre". The private members of the FUI project had no role in either the analyses and interpretation of data or in the discussion of the manuscript.

**References**
1. FAO. La Situation Mondiale des Pêches et de l'aquaculture 2018. In *Atteindre Les Objectifs de Développement Durable*; FAO: Rome, Italy, 2018; ISBN 978-92-5-130692-5.

2. Francis, G.; Makkar, H.P.S.; Becker, K. Antinutritional factors present in plant-derived alternate fish feed ingredients and their effects in fish. *Aquaculture* **2001**, *199*, 197-227, doi:10.1016/S0044-8486(01)00526-9.

3. Tacon, A.G.J.; Hasan, M.R.; Métian, M. Demand and Supply of Feed Ingredients for Farmed Fish and Crustaceans. In *Trends and Prospects; Food and Agriculture Organization of the United Nations*: Rome, Italy, 2011; ISBN 978-92-5-106933-2.

4. Magalhães, R.; Sánchez-López, A.; Leal, R.S.; Martínez-Llorens, S.; Oliva-Teles, A.; Peres, H. Black soldier fly (Hermetia illucens) pre-pupae meal as a fish meal replacement in diets for European seabass (Dicentrarchus labrax). *Aquaculture* **2017**, *476*, 79-85, doi:10.1016/j.aquaculture.2017.04.021.

5. Stadtlander, T.; Stamer, A.; Buser, A.; Wohlfahrt, J.; Leiber, F.; Sandrock, C. *Hermetia illucens* meal as fish meal replacement for rainbow trout on farm. *J. Insects Food Feed* **2017**, *3*, 165-175, doi:10.3920/JIFF2016.0056.

6. Lock, E.R.; Arsiwalla, T.; Waagbø, R. Insect larvae meal as an alternative source of nutrients in the diet of Atlantic salmon (*Salmo salar*) postsmolt. *Aquac. Nutr.* **2016**, *22*, 1202-1213, doi:10.1111/anu.12343.

7. St-Hilaire, S.; Sheppard, C.; Tomberlin, J.K.; Irving, S.; Newton, L.; McGuire, M.A.; Mosley, E.E.; Hardy, R.W.; Sealey, W. Fly prepupae as a feedstuff for rainbow trout, *Oncorhynchus mykiss*. *J. World. Aquac. Soc.* **2007**, *38*, 59-67, doi:10.1111/j.1749-7345.2006.00073.x.

8. Barragan-Fonseca, K.B.; Dicke, M.; van Loon, J.J.A. Nutritional value of the black soldier fly (Hermetia illucens L.) and its suitability as animal feed—A review. *J. Insects Food Feed* **2017**, *3*, 105-120, doi:10.3920/JIFF2016.0055.

9. Henry, M.; Gasco, L.; Piccolo, G.; Fountoulaki, E. Review on the use of insects in the diet of farmed fish: Past and future. *Anim. Feed Sci. Technol.* **2015**, *203*, 1-22, doi:10.1016/j.anifeedsci.2015.03.001.

10. Aksnes, A.; Hope, B.; Jönsson, E.; Björnsson, B.T.; Albrektsen, S. Size-fractionated fish hydrolysate as feed ingredient for rainbow trout (*Oncorhynchus mykiss*) fed high plant protein diets. I: Growth, growth regulation and feed utilization. *Aquaculture* **2006**, *261*, 305-317, doi:10.1016/j.aquaculture.2006.07.025.

11. Wei, Y.; Liang, M.; Mai, K.; Zheng, K.; Xu, H. 1H NMR-based metabolomics studies on the effect of size-fractionated fish protein hydrolysate, fish meal and plant protein in diet for juvenile turbot (*Scophthalmus maximus* L.). *Aquac. Nutr.* **2017**, *23*, 523-536, doi:10.1111/anu.12420.

12. Refstie, S.; Olli, J.J.; Standal, H. Feed intake, growth, and protein utilisation by post-smolt Atlantic salmon (*Salmo salar*) in response to graded levels of fish protein hydrolysate in the diet. *Aquaculture* **2004**, *239*, 331-349, doi:10.1016/j.aquaculture.2004.06.015.

13. Hou, Y.; Wu, Z.; Dai, Z.; Wang, G.; Wu, G. Protein hydrolysates in animal nutrition: Industrial production, bioactive peptides, and functional significance. *J. Anim. Sci. Biotechnol.* **2017**, *8*, 24, doi:10.1186/s40104-017-0153-9.

14. Wang, J.; Yan, X.; Lu, R.; Meng, X.; Nie, G. Peptide transporter 1 (PepT1) in fish: A review. *Aquac. Fish.* **2017**, *2*, 193-206, doi:10.1016/j.aaf.2017.06.007.

15. Verri, T.; Terova, G.; Dabrowski, K.; Saroglia, M. Peptide transport and animal growth: The fish paradigm. *Biol. Lett.* **2011**, *7*, 597-600, doi:10.1098/rsbl.2010.1164.

16. Glencross, B.D.; Booth, M.; Allan, G.L. A feed is only as good as its ingredients: A review of ingredient evaluation strategies for aquaculture feeds. *Aquac. Nutr.* **2007**, *13*, 17-34, doi:10.1111/j.1365-2095.2007.00450.x.

17. Geay, F.; Ferraresso, S.; Zambonino-Infante, J.L.; Bargelloni, L.; Quentel, C.; Vandeputte, M.; Kaushik, S.; Cahu, C.L.; Mazurais, D. Effects of the total replacement of fish-based diet with plant-based diet on the hepatic transcriptome of two European sea bass (*Dicentrarchus labrax*) half-sibfamilies showing different growth rates with the plant-based diet. *BMC Genomics* **2011**, *12*, 522, doi:10.1186/1471-2164-12-522.

18. Panserat, S.; Hortopan, G.A.; Plagnes-Juan, E.; Kolditz, C.; Lansard, M.; Skiba-Cassy, S.; Esquerré, D.; Geurden, I.; Médale, F.; Kaushik, S.; et al. Differential gene expression after total replacement of dietary fish meal and fish oil by plant products in rainbow trout (*Oncorhynchus mykiss*) liver. *Aquaculture* **2009**, *294*, 123-131, doi:10.1016/j.aquaculture.2009.05.013.

19. Casu, F.; Watson, A.M.; Yost, J.; Leffler, J.W.; Gaylord, T.G.; Barrows, F.T.; Sandifer, P.A.; Denson, M.R.; Bearden, D.W. Investigation of graded-level soybean meal diets in red drum (*Sciaenops ocellatus*) using NMR-based metabolomics analysis. *Comp. Biochem. Physiol. Part D Genom. Proteom.* **2018**, *29*, 173-184, doi:10.1016/j.cbd.2018.11.009.
20. Schock, T.B.; Newton, S.; Brenkert, K.; Leffler, J.; Bearden, D.W. An NMR-based metabolomic assessment of cultured cobia health in response to dietary manipulation. *Food Chem.* **2012**, *133*, 90–101, doi:10.1016/j.foodchem.2011.12.077.

21. Abro, R.; Moazzami, A.A.; Lindberg, J.E.; Lundh, T. Metabolic insights in Arctic charr (*Salvelinus alpinus*) fed with zygomycetes and fish meal diets as assessed in liver using nuclear magnetic resonance (NMR) spectroscopy. *Int. Aquat. Res.* **2014**, *6*, 6–63, doi:10.1007/s40071-014-0063-9.

22. Gatesoupe, F.-J.; Fauconneau, B.; Deborde, C.; Madji Hounoum, B.; Jacob, D.; Moing, A.; Corraze, G.; Médale, F. Intestinal microbiota in rainbow trout, *Oncorhynchus mykiss*, fed diets with different levels of fish-based and plant ingredients: A correlative approach with some plasma metabolites. *Aquac. Nutr.* **2018**, *24*, 1563–1576, doi:10.1111/anu.12793.

23. Asakura, T.; Sakata, K.; Yoshida, S.; Date, Y.; Kikuchi, J. Noninvasive analysis of metabolic changes following nutrient input into diverse fish species, as investigated by metabolic and microbial profiling approaches. *PeerJ* **2014**, *2*, e550, doi:10.7717/peerj.550.

24. Mekuchi, M.; Sakata, K.; Yamaguchi, T.; Koiso, M.; Kikuchi, J. Trans-omics approaches used to characterise fish nutritional biorhythms in leopard coral grouper (*Plectropomus leopardus*). *Sci. Rep.* **2017**, *7*, 9372, doi:10.1038/s41598-017-09531-4.

25. Roques, S.; Deborde, C.; Richard, N.; Sergent, L.; Kurz, F.; Skiba-Cassy, S.; Fauconneau, B.; Moing, A. Characterizing alternative feeds for rainbow trout (*O. mykiss*) by ¹H NMR metabolomics. *Metabolomics* **2018**, *14*, 305, doi:10.1007/s11306-018-1454-5.

26. Jasour, M.S.; Wagner, L.; Sundekilde, U.K.; Larsen, B.K.; Rasmussen, H.T.; Hjermitslev, N.H.; Hammershøj, M.; Dalsgaard, A.J.T.; Dalsgaard, T.K. Fishmeal with different levels of biogenic amines in aquafeed: Comparison of feed protein quality, fish growth performance, and metabolism. *Aquaculture* **2018**, *488*, 80–89, doi:10.1016/j.aquaculture.2018.01.030.

27. De la Higuera, M. Effects of Nutritional Factors and Feed Characteristics on Feed Intake. In *Food Intake in Fish*; Houlihan, D., Boujard, T., Jobling, M., Eds.; Blackwell science: Oxford, UK, 2001; pp 250–268.

28. Kaushik, S.J.; Seiliez, I. Protein and amino acid nutrition and metabolism in fish: Current knowledge and future needs. *Aquac. Res.* **2010**, *41*, 322–332, doi:10.1111/j.1365-2109.2009.02174.x.

29. Wright, P.A. Nitrogen excretion: Three end products, many physiological roles. *J. Exp. Biol.* **1995**, *198*, 273–281.

30. Ip, Y.K.; Chew, S.F. Ammonia Production, Excretion, Toxicity, and Defense in Fish: A Review. *Front. Physiol.* **2010**, *1*, 134, doi:10.3389/fphys.2010.00134.

31. Yamamoto, T.; Unuma, T.; Akiyama, T. Postprandial changes in plasma free amino acid concentrations of rainbow trout fed diets containing different protein sources. *Fish. Sci.* **1998**, *64*, 474–481, doi:10.2331/fishsci.64.474.

32. Salze, G.P.; Davis, D.A. Taurine: A critical nutrient for future fish feeds. *Aquaculture* **2015**, *437*, 215–229, doi:10.1016/j.aquaculture.2014.12.006.

33. Yokoyama, M.; Nakazoe, J.-I. Accumulation and excretion of taurine in rainbow trout (*Oncorhynchus mykiss*) fed diets supplemented with methionine, cystine and taurine. *Comp. Biochem. Physiol. A Physiol.* **1992**, *102*, 565–568, doi:10.1016/0300-9629(92)90210-H.

34. Kataoka, H.; Ohnishi, N. Occurrence of Taurine in Plants. *Agric. Biol. Chem.* **2014**, *50*, 1887–1888, doi:10.1080/00021369.1986.10867664.

35. Belghit, I.; Liland, N.S.; Gjesdal, P.; Biancarosa, I.; Menchetti, E.; Li, Y.; Waagbø, R.; Krogdahl, Å.; Lock, E.-J. Black soldier fly larvae meal can replace fish meal in diets of sea-water phase Atlantic salmon (*Salmo salar*). *Aquaculture* **2019**, *503*, 609–619, doi:10.1016/j.aquaculture.2018.12.032.

36. Krogdahl, Å.; Penn, M.; Thorsen, J.; Refstie, S.; Bakke, A.M. Important antinutrients in plant feedstuffs for aquaculture: An update on recent findings regarding responses in salmonids. *Aquac. Res.* **2010**, *41*, 333–344, doi:10.1111/j.1365-2109.2009.02426.x.

37. Gómez-Requeni, P.; Mingarro, M.; Calduch-Giner, J.A.; Médale, F.; Martin, S.A.M.; Houlihan, D.F.; Kaushik, S.; Pérez-Sánchez, J. Protein growth performance, amino acid utilisation and somatotropic axis responsiveness to fish meal replacement by plant protein sources in gilthead sea bream (*Sparus aurata*). *Aquaculture* **2004**, *232*, 493–510, doi:10.1016/S0044-8486(03)00532-5.

38. Lin, C.Y.; Wu, H.; Tjeerdema, R.S.; Viant, M.R. Evaluation of metabolite extraction strategies from tissue samples using NMR metabolomics. *Metabolomics* **2007**, *3*, 55–67, doi:10.1007/s11306-006-0043-1.
39. Jacob, D.; Deborde, C.; Lefebvre, M.; Maucourt, M.; Moing, A. NMRProcFlow: A graphical and interactive tool dedicated to 1D spectra processing for NMR-based metabolomics. *Metabolomics* **2017**, *13*, 36, doi:10.1007/s11306-017-1178-y.

40. Kullgren, A.; Samuelsson, L.M.; Larsson, D.G.J.; Bjornsson, B.T.; Bergman, E.J. A metabolomics approach to elucidate effects of food deprivation in juvenile rainbow trout (*Oncorhynchus mykiss*). *Am. J. Physiol. Regul. Integr. Comp. Physiol.* **2010**, *299*, R1440–R1448, doi:10.1152/ajpregu.00281.2010.

41. Shen, G.; Huang, Y.; Dong, J.; Wang, X.; Cheng, K.-K.; Feng, J.; Xu, J.; Ye, J. Metabolic effect of dietary taurine supplementation on Nile tilapia (*Oreochromis nilotictus*) evaluated by NMR-based Metabolomics. *J. Agric. Food Chem.* **2018**, *66*, 368–377, doi:10.1021/acs.jafc.7b03182.

42. Sumner, L.W.; Amberg, A.; Barrett, D.; Beale, M.H.; Beger, R.; Daykin, C.A.; Fan, T.W.-M.; Fiehn, O.; Goodacre, R.; Griffin, J.L.; et al. Proposed minimum reporting standards for chemical analysis Chemical Analysis Working Group (CAWG) Metabolomics Standards Initiative (MSI). *Metabolomics* **2007**, *3*, 211–221, doi:10.1007/s11306-007-0082-2.

<figure><img src="image_8.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).