Black Soldier Fly Larvae Meal as Alternative to Fish Meal for Aquaculture Feed

: Hermetia illucens meal (HIM) as ingredient in feed represents a way to achieve more sustainable food production. The aim was to characterize the chemical, microbiological and organoleptic characteristics of four diets for Sparus aurata , isoenergetic and isoproteic, containing 0%, 25%, 35% and 50% of HIM in substitution of ﬁsh meal (FM). Analyses were carried out using gas chromatography for fatty acids and amino acids, ICP-OES for minerals and liquid chromatography for aﬂatoxins and following International Organization for Standardization methods for microbial ﬂora. E-sensing analysis of the diets was evaluated using an artiﬁcial sensory platform (E-eye, E-nose and E-tongue). The chemical results were submitted to a one-way ANOVA while Principal Component Analysis (PCA) of the e-sensing data was performed. No signiﬁcant differences were observed for polyunsaturated fatty acids, thrombogenic and peroxidation indices among the diets. The replacement of FM with HIM increased the content of lysine, methionine, isoleucine, leucine, threonine and valine, while phosphorus, calcium and sodium content decreased ( p < 0.01) as the percentage of HIM increased. Lead was signiﬁcantly below the maximum level set by the EU regulation. The diets showed good hygienic and sanitary quality. The artiﬁcial senses permitted distinguishing color, odor and taste among the diets. Data allow considering Hermetia illucens as alternative protein source in ﬁsh nutrition.


Introduction
Among the costs of the aquaculture farms, feeding represents the largest portion (about 60%) of the economic balance; therefore, the economic success of fish production sector is mainly linked to the use of low-cost nutritionally balanced diets [1]. Therefore, the goal of aquaculture nutritionists, as well as fish farmers, is to obtain a very good conversion ratio to cover the cost of feeding. Proteins with a well-balanced presence of essential and non-essential amino acids are the most important nutrients for the maintenance, growth and feed efficiency of fish. They are the nutrients with the highest cost, and, therefore, their inclusion in the fish feed plays an important role in the overall feed costs. Fish meal, characterized by a high protein content, an excellent amino acid profile, a low carbohydrate level and a high digestibility, is considered the most important feeding source in fish nutrition [2]. However, the significant increase in the price of fish meal, the high protein nutritional needs of fish and the alternating availability in fish meal supply, in recent decades, have led to the study of alternative protein sources in aquafeeds [3].
In relation to the increasing demand for protein sources in animal feeding, there is a great interest towards unconventional sources. In recent years, industrial by-products, co-products, insects, seaweed ingredients [4] and ex-food or former food products have been investigated as alternative ingredients for livestock and aquaculture feeding.
Interest in insects as feed ingredient for terrestrial and aquatic animals continues to grow [5][6][7][8]. At present, the exploiting of insects as feed ingredients is not in direct competition with food production. The use of insect meal as a substitute of fish meal seems to represent interesting perspectives for limiting the environmental impact of some production system, such as aquaculture, and for contributing to a circular economy and a "zero waste" society [9].
Among the insect species used as unconventional protein source for fish feeding, Hermetia illucens L. is of the most interesting source for its sustainability related to its capacity to convert organic waste material into biomass containing proteins (40-45%) with high biological value [10], fat (30-35%) with fatty acids of nutritional interest and ash (11-15%) with high mineral concentrations and a high Ca/P ratio [11]. In terms of protein quality, Hermetia illucens larvae contain a favorable essential amino acid profile closer to fish meal than that of soybean meal [5]. The mineral profile and fatty acid composition of the Hermetia illucens larvae were found to be influenced by the diet [12][13][14]. Liland et al. [14] reported that Hermetia illucens larvae are unable to synthesize polyunsaturated fatty acids, and, therefore, the presence of linoleic acid and alpha-linolenic acids, as well as of eicosapentaenoic and docosahexaenoic in the larvae most likely originates from the substrate. Furthermore, these authors observed that the presence of seaweed in the feeding media can enrich the larvae with macro-and microelements, making the insects a good source of minerals [14].
Nevertheless, Weththasinghe et al. [15] observed a linear decrease in protein and lipid digestibility, protein efficiency ratio and lipid retention in extruded diets for Atlantic salmon as the level of dietary Hermetia illucens meal increases.
The use of Hermetia illucens larvae meal (HIM) as a component of feed is a way to achieve more sustainable food production. According to current regulation in Europe [16], the use of HIM is mainly permitted in aquaculture. The optimal level of dietary substitution of fish meal for HIM varies considerably across studies, ranging from 25% to 100%, probably in relation to the different quality of larvae meal, fish species and diet formulation. Furthermore, to obtain a higher-quality fish, it is important to provide good quality and pathogen-free feeds to the fish [17]. In this view, the extrusion method destroys undesirable microbial flora, enzymes and anti-nutritional factors [18][19][20] and improves the nutritional value of final products and the apparent absorption of minerals [21,22].
Fish have strongly developed chemosensory and chemical signaling systems due to their living in an aquatic environment. The olfactory and gustatory systems comprise the main chemosensory pathways [23]. The combination of artificial senses (Electronic Nose, an Electronic Tongue and an Electronic Eye) proves to be a powerful tool to distinguish different organoleptic profiles related to the different chemical compositions of aquaculture feeds [24].
To further characterize the feasibility of HIM as a unconventional protein source in aquafeeds, the aim of this study was to characterize the chemical, microbiological and mineral composition of diets for Spaurus aurata L. containing HIM as a partial replacement for fish meal (FM). The effects of inclusion of HIM on the organoleptic characteristics of the diets were also investigated with the aim of providing detailed sensory information on the HIM-containing fish meal useful for feed industry.

Materials and Methods
Four experimental diets were formulated to satisfy the nutritional needs of Sparus aurata. The diets were isoenergetic (about 22 MJ/kg gross energy), isonitrogenous (about 43 g/100 g, as fed) and isolipidic (about 19 g/100 g, as fed). A basal diet (HIM0) with fish meal (FM), as exclusive protein source of animal origin, was prepared. FM was partially replaced with defatted Hermetia illucens meal at 25%, 35% and 50% (as fed basis) in the other three diets (HIM25%, HIM35%, HIM50%) corresponding to the inclusion levels of Sustainability 2021, 13, 5447 3 of 17 0%, 7.9%, 11% and 15.7% of HIM, respectively. The other components of the formulas were adapted to obtain diets with the same energetic content.
Diets were prepared by SPAROS Lda (Olhao, Portugal); all dietary ingredients were ground, mixed and extruded using die with 4 mm diameter; the oils were added using a vacuum coating technology. The ingredients of the diets and the proximate composition of the diets (HIM0, HIM25%, HIM35% and HIM50%) are reported in Table 1.

Fatty Acid Analysis
Triplicate feed samples were analyzed for the fatty acid composition. Each sample (ca. 2.5 g) added with sodium sulfate (1 g) was pounded manually. The lipids were extracted for 6 h with petroleum ether by a Soxtec™ 8000 Extraction system (FOSS, Padua, Italy). The fatty acid methyl esters (FAMEs) were produced from aliquots of lipids. In detail, sulfuric acid-methanol (1:9, v/v) reagent (2 mL) was added to the extracted lipid samples, and they were then heated at 100 • C for 1 h [25]. FAMEs were analyzed by a Trace 1310 chromatograph (Thermo Fisher Scientific, Milan, Italy) equipped with a flame ionization detector (FID) and a fused silica capillary column (30 m × 0.25 mm I.D., 0.25 µm film thickness) (Omegawax 250; Supelco, Bellefonte, PA, USA) maintained at 100 • C for 5 min, from 100 to 240 • C at 4 • C/min and final isotherm of 240 • C (20 min). Injector and detector temperatures were 250 • C. Injection volume and split ratio were 0.5 µL and 1:50, respectively. The carrier gas was helium at a flow rate of 1 mL/min. Data acquisition was carried out by a Chromeleon Software (Thermo Fisher Scientific, Milan, Italy). The identification of individual compounds was carried out by comparing their retention times with those of standards (mix 37 FAMEs, Supelco, Inc., Bellefonte, PA, USA). The results were expressed as g/100 g of the total fatty acids identified.

Amino Acid Analysis
Triplicate feed samples were analyzed for the amino acid composition. For the amino acid analysis, protein hydrolysis and derivatization were performed prior to the separation by GC-FID. Each sample (about 0.25 g) was hydrolyzed in 10 mL of a HCl solution (6 M) at 110 • C for 24 h. During the acid hydrolysis, the asparagine and glutamine were converted to aspartic and glutamic acids [29]; therefore, they were calculated as the sum of the aspartic acid plus asparagine and of the glutamic acid plus glutamine. For cysteine analysis, prior to the acid hydrolysis, a preliminary oxidation was performed for the deamination. Each sample was treated with formic acid:hydrogen peroxide (1:20, v/v) reagent (2 mL) for 30 min at room temperature. Then, a hydrolysis step with a HCl solution (6 M) was performed [30,31]. For the tryptophan analysis, each sample was hydrolyzed in 10 mL of a NaOH solution (4 M) at 112 • C for 16 h, and, after the hydrolysis, each sample was cooled and neutralized with acetic acid [32]. For the chromatographic analysis, procedures for purification, pre-column derivatization and qualitative and quantitative analyses of each amino acid were performed using the EZ:Faast Kit (Phenomenex, Torrance, CA, USA). A Trace 1310 chromatograph (Thermofisher, Waltham, MA, USA) was used, with a flame ionization detector (FID) and a ZB-AAA Amino Acid column (10 m × 0.25 mm ID); the oven temperature was programmed from 110 to 320 • C at 32 • C/min, with a final isotherm of 320 • C (1 min). Injector and detector temperatures were 250 and 320 • C, respectively. Injection volume and split ratio were 2.5 µL and 1:15, respectively.

Aflatoxin Analysis
The analysis of aflatoxins was carried out following the EN ISO method [33] and the criteria suggested by European regulation [34]. For the extraction of the aflatoxins, each sample (about 25 g) of the experimental diet was treated with 125 mL of a methanol:water mixture (7:3, v/v) and 5 g of NaCl. The extract was filtered, diluted with water and passed through an immunoaffinity column (Vicam) containing specific antibodies for aflatoxins B1, B2, G1 and G2. The aflatoxins were isolated, purified and concentrated on column and then recovered with methanol. The above procedure was performed in triplicate on each experimental diet. Aflatoxins were measured by a RP-HPLC coupled to a fluorescence detector (RF), using a post-column derivatization (Kobra cell system). Chromatographic separation was performed using a Nexera LC System (Shimadzu, Milan, Italy), equipped with a Luna column C18 250 × 4.6 mm (lenght × i.d.) and 5 µm of particle size (Phenomenex, Torrance, CA, USA). The mobile phase for the isocratic separation was a mixture of water:acetonitrile:methanol (3:1:1, v/v/v) with 0.35 mL of nitric acid 4 M and 120 mg/L of potassium bromide. The RF wavelength was set at 365 nm for excitation and 435 nm for emission. The injection volume was 50 µL and the amount of aflatoxins was calculated using an external standard method. The LC method for aflatoxin analyses has been validated for the simultaneous chromatographic determination of total (B1, B2, G1 and G2) and B1 aflatoxins.
A calibration curve was constructed for each aflatoxin; the linearity was also tested in the range of 0.05-22 µg/kg, providing a correlation coefficient (R2) of ≥0.9996. The Limit Of Detection (LOD, between 0 and 0.05 µg/kg) and the Limit Of Quantification (LOQ, between 0.05 and 0.13 µg/kg) were calculated by the signal-to-noise (S/N) ratio, which should be greater than 3 and 10, respectively, according to the IUPAC criteria.

Mineral Element Analysis
Triplicate feed samples were analyzed for the mineral composition. Each sample was weighed (about 0.5 g) into an acid-prewashed PTFE vessels and 7 mL of HNO3 at 65% and a Rhenium internal standard (1 mL) were added. The mixture was digested with 1 mL of H 2 O 2 at 30% using a closed-vessel microwave digestion system (Ethos 1, Milestone, Bergamo, Italy). To validate the analytical method, a Standard Reference Material of spinach leaves (SRM, NIST-1570a) obtained from the National Institute of Standards and Technology (Gaithersburg, MD, USA), was digested using the same analytical method described for the feed samples. An Avio200 ICP-OES instrument (Perkin Elmer, Waltham, MA, USA) equipped with a vertical DualView optical system and a S10 autosampler (Perkin Elmer, Waltham, MA, USA) was used to analyze the mineral content. Table 2 shows the recommended analytical lines length used to perform element analyses; the Argon line at 420.069 nm was used as an internal standard. The applied operational conditions are listed in Table 3. Data were processed using a PerkinElmer Syngistix™ for ICP software (Perkin Elmer, Waltham, MA, USA).  Table 3. Operational conditions of the ICP-OES.

Parameter Conditions
Radiofrequency power (W) 1500 Plasma gas flow (L/min) 9 Auxiliary gas (L/min) 0.2 Nebulizer gas (L/min) 0.7 Sample uptake (mL/min) 1 The position of the torch was optimized prior to the analytical phase using the optical optimization procedure of Syngistix™ ICP software with Mn analytical line. All the quantitative measurements were made against external calibration curves constructed from a standard solution of 0.05, 0.25 and 1 ppm of Perkin Elmer (Waltham, MA, USA) for ICP analysis. A Milli-Q ultrapure (Merck Millipore, Merck KGaA, Darmstadt, Germany) water system was used to produce water at 1.8 MΩ/cm for the preparation of solutions and to dilute samples as needed. The calibration curves for all elements were established using the calibration blank and the reagent blank, and all of them resulted with correlation coefficients (r 2 ) better than 0.999; the Detection Limits (DLs) of this procedure were determined by Sustainability 2021, 13, 5447 6 of 17 analyzing a matrix blank, which consisted of the same reagents and quantities as those used for sample preparation.

E-Sensing Analysis
The feed samples were analyzed using an artificial sensory platform consisting of an E-eye, E-nose and E-tongue.
E-eye: The image was acquired with an artificial vision system (Iris visual analyzer 400, Alpha MOS, Toulouse, France) equipped with a high-resolution charge-coupled device (CCD) camera with 16 million colors. Image acquisition was performed using a top illumination and a white tray at the bottom to easily remove the background contribution (threshold selection applied: R 0-145, G 0-121, B 0-109). Each sample was ground, and the powder was placed and flattened on a plastic Petri dish (diameter 92 mm, height 7 mm). For each sample, 15 images on 15 freshly prepared dishes were acquired. Color spectra were calculated by selecting only the contributions greater than 1%. E-nose: Odor analysis was performed by an electronic nose device (FOX 4000, Alpha MOS, Toulouse, France) equipped with 18 MOS (metal-oxide semiconductor) gas sensors and an automatic headspace sampler (HS100). For each sample, 15 replicates were prepared by placing 2 g of freshly ground fish feed into 10 mL headspace sealed vials. All parameters of the instrument are reported in Table 4. Volume 500 µL Speed 500 µL/s E-tongue: Artificial taste analysis was performed using a commercially available electronic tongue (Astree, Alpha MOS, Toulouse, France) equipped with a set of seven potentiometric sensors (ANS, PKS, CTS, NMS, CPS, ANS and SCS), an Ag/AgCl reference electrode (Metrohm, Pte Ltd., Singapore), a mechanical stirrer and a 48-position autosampler. Five grams of each sample were ground and placed in 50 mL of deionized water for 15 min and centrifuged at 3000 rpm for 30 min. The solution was filtered and placed in a 25 mL beaker for the analysis. Single sample analysis was repeated 30 times to obtain the most stable sensor response and the last 15 measurements were considered to perform data processing. The signal was acquired every second for 120 s and the average intensity of the last 20 s was measured. Prior to measurement, the sensors were conditioned using one of the samples as a standard.

Statistical Analysis
The chemical data were analyzed by a one-way ANOVA, using the XLSTAT statistical package [43]. The percentage integration of insect meal (HIM0, HIM25%, HIM35% and HIM50%) was used as a fixed effect. Separation of means was assessed by Tukey's test, and differences were significant if p < 0.05.
A Principal Component Analysis (PCA) of the sensory profile data was performed by Alpha Soft V12.4 (Alpha-MOS, Toulouse, France) to evaluate the discrimination ability between the four experimental diets. The effectiveness of discrimination was assessed by evaluating the discrimination index (DI), which gives the quality of discrimination through an indication of the surface between groups. The DI is calculated automatically by the instrument's software according to the following formula: The DI reaches a maximum value of 100 when the groups are completely resolved.

Fatty Acid Profile
The fatty acid profile of the four experimental diets is shown in Table 5. Saturated fatty acids did not show any significant (p > 0.05) differences, with the exception of the lauric acid (C12:0), which showed significantly (p < 0.01) higher values in the HIM35% and HIM50% diets than those observed in the HIM0 and HIM25% diets and for the palmitic acid (C16:0) which showed a significant (p < 0.05) higher value in the HIM0 diet than those observed in the diets containing Hermetia illucens meal. Monounsaturated fatty acids show similar content among the experimental diets as well as the polyunsaturated fatty acids, of both the omega 3 and omega 6 series. Table 6 shows the fatty acid classes; the sum of saturated, monounsaturated and polyunsaturated fatty acids; and some indices of nutritional interest, namely the atherogenic (AI) and thrombogenic (TI) indices, the peroxidation index (PI) and the hypocholesterolaemic/hypercholesterolaemic ratio. No significant (p > 0.05) difference was observed among the fatty acid classes, with the exception of the saturated fatty acids which show the highest (p < 0.05) level in the control diet (HIM0). The sum of the eicosapentaenoic (EPA) and docosahexaenoic (DHA) acids was similar among the diets. Similar values were recorded for AI, TI and PI. The H/H ratio showed significantly (p < 0.05) higher values in the HIM35% and HI50% diets than that recorded in the HIM0 diet while the HIM25% diet showed a value of the H/H ratio similar to those of the other diets. HIM0, fish meal; HIM25%, HIM35% and HIM50%, Hermetia illucens meal at 25%, 35% and 50% substitution rate of fish meal, respectively. # The concentration of fatty acid is expressed as g/100 g, considering 100 g the sum of the areas of all FAME identified. Mean values with different letters within the same row are significantly different, A-C at p < 0.01 and a and b at p < 0.05.  Table 7 reports the amino acid composition of the experimental diets. Twenty amino acids were identified and quantified; ten of these belong to the indispensable amino acids and ten to dispensable ones. Among the indispensable amino acids, six amino acids (isoleucine, leucine, lysine, methionine, threonine and valine) showed significantly (p < 0.01) higher values in the diets containing insect meal, while histidine showed significantly (p < 0.01) higher values in the HIM35% and HIM50% diets than those observed in the HIM0 and HIM25% diets. Arginine, phenylalanine and tryptophan showed the highest (p < 0.01) values in the control diet (HIM0). Among the dispensable amino acids, hydroxylisine, hydroxyproline and tyrosine showed similar values among the experimental diets. Glutamic acid plus glutamine showed significantly (p < 0.01) higher values in all diets containing the Hermetia illucens meal, while proline and serine showed significantly (p < 0.01) higher values in the diets in which FM has been replaced with HIM at 35% and 50%. The levels of aspartic acid plus asparagine and that of glycine were significantly (p < 0.01) higher in the HIM25% and HIM35% diets than those observed in the HIM0 and HIM50% diets. Cysteine showed the highest (p < 0.05) value in the HIM25% diet and alanine the highest (p < 0.01) level in the HIM0 diet.

Mineral Element Profile
In Table 8, the average values of minerals in the experimental diets are reported as: macroelements (phosphorus, calcium, potassium, magnesium and sodium), whose needs by the body are in large amounts, microelements (copper, zinc, manganese, iron, etc.), whose needs by the body are in small amounts [44] and toxic metal (lead). The amount of phosphorus, calcium and sodium decreased significantly (p < 0.01) by increasing the amount of HIM to replace fish meal. The presence of insect meal over 25% as a substitute for fish meal (HIM35% and HIM50%) resulted in significantly (p < 0.01) lower levels of potassium, although the HIM35% and HIM50% diets showed similar (p > 0.05) values between them. The Ca/P ratio in the four experimental diets was calculated due to the antagonist effect of these two macroelements [45]. Significantly (p < 0.01) higher levels for the Ca/P ratio were observed in the HIM0 and HIM25% diets than those of the HIM25% and HIM50% diets. Magnesium was significantly (p < 0.01) higher in the HIM25% diet than in HIM0, HIM35% and HIM50% diets, while these diets did not show (p > 0.05) differences among them. Iron, copper and zinc showed significantly (p < 0.01) higher levels in the HIM35% diet than those of the HIM0, HIM25% and HIM50% diets, except the HIM25% diet, which showed a similar content to the HIM50% diet. A significantly (p < 0.01) higher level of manganese was observed in the HIM25% diet than those of the HIM0, HIM35% and HIM50% diets. The increase in the level of HIM to replace fish meal resulted in a significant (p < 0.05) increase in the amount of chromium. Boron was significantly (p < 0.01) higher in the HIM25% and HIM35% diets compared to the other diets; lead showed the highest (p < 0.01) value in the HIM50% diet (Table 8).

Microbiological Profile
No significant difference on the load of each microbiological parameter was observed in the experimental diets. For all diets, only a few aerobic colonies with a load below 70 CFU/g (HIM0 = 65 CFU/g, HIM25% = 70 CFU/g; HIM35% = 50 CFU/g; HIM50% = 60 CFU/g) was observed. The counts of Enterobacteriaceae, coliforms, Clostridium spp., yeasts and molds were always under the LOD. Enterobacteriaceae, Salmonella spp., L. monocytogenes and Clostridium spp. were not detected.

E-Sensing Profile
Concerning data provided by the artificial sensory platform, the first step was to separately perform PCAs on the data from E-nose and E-tongue sensors and on the colors of the E-eye code. The next step was to look for the most effective way to combine the data provided by the E-eye, E-nose and E-togue in order to improve the discrimination capability. An intermediate fusion level was adopted in this study. The sensor data with the highest discrimination power were chosen; in particular, the data of four E-nose sensors (LY2/G, LY2/AA, P30/1 and T40/1), three E-tongue sensors (AHS, CTS and NMS) and four colors extracted from the E-eye (codes 1620, 1890, 1891 and 2147) were chosen. These datasets were reduced, due to the different data size, and a new PCA was performed and the Discrimination Index (DI) calculated. The result is shown in Figure 1.

E-Sensing Profile
Concerning data provided by the artificial sensory platform, the first step was to separately perform PCAs on the data from E-nose and E-tongue sensors and on the colors of the E-eye code. The next step was to look for the most effective way to combine the data provided by the E-eye, E-nose and E-togue in order to improve the discrimination capability. An intermediate fusion level was adopted in this study. The sensor data with the highest discrimination power were chosen; in particular, the data of four E-nose sensors (LY2/G, LY2/AA, P30/1 and T40/1), three E-tongue sensors (AHS, CTS and NMS) and four colors extracted from the E-eye (codes 1620, 1890, 1891 and 2147) were chosen. These datasets were reduced, due to the different data size, and a new PCA was performed and the Discrimination Index (DI) calculated. The result is shown in Figure 1.  The first two principal components calculated (PC1 and PC2) explain 99.4% of the total variance and show a DI of 93, highlighting a clear difference among the four diets, in relation to the substitution of fish meal with HIM. The combination of sensor responses with the highest discrimination power of the three artificial senses completely separates the four feed groups. Moreover, the loading plot helps identify the variables accountable for clustering of our dataset, by providing a numerical value which represents the contribution of each original variable to the score plot. Along the PC1 axis, from left to right, the groups are separated mainly by color codes 1890, 1891 and 2147, which are the most distinctive for the HIM0 group. On PC1, the volatile component of the different groups is also distinguished. Figure 1 shows the contribution of LY2/G, LY2/AA, P301 and T401. These sensors show a clear separation between HIM35% and HIM50% diets and the HIM0 and HIM25% diets. Regarding the taste profile, the three selected sensors (AHS, CTS and NMS) help to improve discrimination between groups. In particular, the AHS and CTS sensors have an important contribution on the HIM0 and HIM25% diets, while the NMS sensor mainly distinguishes the groups with higher insect meal content (HIM35% and HIM50%).

Discussion
Regarding the origin of fatty acids in the larvae, some studies report the possibility of modifying the fatty acid profile through the diet, while others the possibility of an endogenous synthesis of fatty acids. Knowledge of the fatty acid pathway is of great importance if HIM is used as an ingredient in animal feed. Our results are similar to the fatty acid profile determined by Belghit et al. [9] in four experimental diets formulated for Atlantic salmon with an increasing substitution of fish meal with HIM. St-Hilaire et al. [13] observed a reduction of alpha-linolenic acid (C18:3n3), eicosapentaenoic acid (C20:5n3) and docosahexaenoic acid (C22:6n3) in fish fed a diet containing HIM. This could represent a problem for producers and consumers of the seafood supply chain [46]. However, our results show similar content of the alpha-linolenic acid, eicosapentaenoic acid and docosahexaenoic acid among the diets, despite the increasing addition of HIM in the feed. The content of these fatty acids, essential for the growth and development in fish [47] and associated with disease prevention and health promotion for humans, thanks to the production of anti-inflammatory eicosanoids [48], could be associated with the use, in feeding of the larvae, of a substrate containing fish offal and algae [12,14], since Hermetia illucens is not able to synthesize PUFA [14]. In fact, only plants and marine algae possess the enzymes necessary for the synthesis of linoleic acid (C18: 2n6) and alpha-linolenic acid (the precursor of EPA and DHA) [49]; therefore, similar concentrations of n6-PUFA (mainly linoleic acid) and n3-PUFA among the four diets confirm our hypothesis on the influence of the feeding substrate of Hermetia illucens larvae [46]. Another interesting result concerned the similar values recorded for the health lipid indices (atherogenic and thrombogenic indices), the peroxidation index and the H/H ratio. The atherogenic and thrombogenic indices (AI and TI) take into account, in the formulas used for their calculation, the contribution that each fatty acid has on human health and, in particular, on the probability of influencing the incidence of cardiovascular diseases [26]. The decrease in PI, e.g., during feed storage, indicates the oxidative degradation of PUFA in primary and secondary oxidation products. This process results in a loss of shelf-life, nutritional value and feed safety, as well as reduced consumer acceptability [50]. Finally, the low H/H ratio observed in diets containing HIM suggests some positive nutritional effects on animals.
As regards the indispensable amino acids, the substitution of fish meal with Hermetia illucens increased the lysine (EAA) and methionine content in all the levels of inclusion (HIM25%, HIM35% and HIM50%). In fact, Hermetia illucens meal has high levels of amino acids such as methionine and lysine [30]. Furthermore, isoleucine, leucine, threonine and valine showed higher levels in all experimental diets containing Hermetia illucens meal; this fact is nutritionally interesting, as fish cannot synthesize these amino acids de novo [51]. Among the dispensable amino acids (NEAA), the inclusion of HIM at 35% and 50% levels resulted in an increase in proline and tyrosine (NEAA), as observed by De Marco et al. [30] and by Iaconisi et al. [51].
Minerals are present in low quantities in the feed despite having important metabolic roles [22,44]. Mineral elements in fish are involved in many biochemical processes as enzyme cofactors and activators, in the formation of skeletal muscles, in the transmission of nerve impulses and in the acid-base chemical balance [44,52,53]. The low levels of calcium, potassium and sodium observed in the formulation containing Hermetia illucens meal do not represent a problem in the feeding of fish that can readily obtain these minerals from the surrounding water environment [44,53]. On the contrary, the phosphorus requirement must be satisfied through the diet as the fish cannot readily obtain it from the surrounding aquatic environment, which contains a low amount of it [54,55]. Phosphorus and calcium are generally combined together in the fish body, so an adequate Ca/P ratio in the diet is crucial to ensure healthy bone formation and growth performance of fish [45,52]. In relation to their antagonistic effect [56], the dietary Ca/P ratio has to be taken into consideration, as an excess of Ca or P can cause mineral imbalances that can affect the absorption of other minerals (e.g., an excess of Ca into the diet can reduce the absorption of Zn, Fe and Mn), causing some important consequences for bone development, which can be adversely affected when this ratio increases, causing anomalies in mineral homeostasis and bone mass [57] or environmental contamination (excess P is excreted). In our study, the Ca/P ratio of all diets was close to 1, as recommended for fish [58,59].
As for the microelements, such as zinc, essential for the growth and development in fish, cofactor of many enzymes and necessary for the activity of the antioxidant enzyme superoxide dismutase [44,60]; iron, necessary for the blood and muscle pigments [44] and actively involved in oxidation-reduction reactions [60]; and copper, involved in enzymatic activities and oxygen transportation [60], the highest levels were observed in the diet containing a quote of 35% of defatted Hermetia illucens meal. This result would seem to highlight that the HIM35% diet is the best formula from a nutritional point of view, considering also that fish cannot obtain adequate quantities of these microelements from water [59,61] and that they must receive them through the diet [55].
The marine environment can convey toxic metals, naturally contained or introduced by various human activities through seafood [62]. Fish, both freshwater and marine water, that are at the top of the food chain, are extremely sensitive to exposure to lead (Pb), which is a highly toxic metal [63]. Its toxicity depends on various factors such diet and environment [63]. It seems that Sparus aurata is much more sensitive to waterborne lead exposure than in the diet [64]. However, considering the toxicity of lead, the European regulation [65] has set a maximum level of lead for complete animal feed of 5 ppm, which is considerably higher than the values found in our experimental diets.
Fish feeds become rapidly colonized by environmental microbes. Furthermore, all insects are colonized by microorganism and the insect microbiota is generally different from microorganism in the external environment, including ingested food [66,67]. Usually, the presence in commercial fish feeds of spoilage and pathogen bacteria such as Salmonella spp. and Escherichia coli can be related to poor hygienic storage condition, environmental contamination and problems in the extrusion treatment [68,69]. Our results on the microbiological profile of all diets confirm the effectiveness of the extrusion treatment in reducing the microbial load. The low charge observed for aerobic bacteria could be related to a secondary contamination after extrusion treatment. As recent studies have shown that aflatoxin contamination of animal feeds is a frequent issue, the European regulation has set maximum residue limits for total aflatoxins (B1, B2, G1 and G2) and aflatoxin B1 in animal feed [70]. To confirm the effectiveness of extrusion treatment and the good environmental condition during the storage, the total aflatoxins (B1, B2, G1 and G2) and aflatoxin B1 were found to be below the limit of detection in all experimental diets [71]. This is of particular interest because, in aquaculture, their presence in the diet can destroy the availability of certain nutrients, such as vitamin C and thiamine, reducing the immune defense of fish [72]. In addition, the contamination by mycotoxins of fish feed has been reported to cause intoxications [69], tissue abnormalities or liver injury, liver tumor, decreased growth rate and appetite [73] with a decrease of production efficiency and weight of the caught product and an increase in medical costs [74].
The E-nose and E-tongue sensors are non-specific and partially cross-sensitive; therefore, the value of each sensor is not directly related to each other. To better understand the data obtained by the artificial senses and improve the discrimination capability, it is necessary to optimize the analysis through a fusion process often categorized in a three-level model which includes the distinction of low, intermediate and high fusion levels [75,76]. In our study, the use of an intermediate fusion level helped improve the discrimination power. First, some relevant features were extracted from each data source separately, and then they were concatenated into a single array, which was used for multivariate classification and regression [77]. Selected color codes and E-nose and E-tongue sensors identified differences between groups in relation to the percentage of HIM integration. The volatile component was mainly distinguished by LY2/G, LY2/AA, P301 and T401, which are sensors primarily sensitive to volatile organic compounds derived from proteins and lipids such as nitrogencontaining compounds, hydrocarbons and aldehydes [78,79]. These sensors showed a clear separation of the HIM35% and HIM50% groups from HIM0 and HIM25%, in relation to the amino acidic content of the diets. Our observations appear interesting as the use of alternative ingredients in fish diets, in relation to their chemical composition and their level of inclusion, can reduce the acceptability of feeds, even if nutritionally balanced. In fact, different chemical substances can influence fish feeding behavior by acting as attractants through smell or taste [80]. Fish, in general, have a well-defined olfactory sensitivity to amino acids [81][82][83][84][85], which help them locate and identify food. Glycine and alanine are potent odorants that can stimulate feeding behavior by increasing food intake [86,87]. Our results show that HIM-integrated feeds had greater amounts of glycine and alanine com-pared to HIM0 diet; this probably contributed to the olfactory discrimination performed by the E-nose sensors. Regarding the taste profile, the three sensors selected mainly represent umami taste (NMS), sourness (AHS) and saltiness (CTS). The results of PCA show that the HIM0 and HIM25% diets have a more salty and sour taste, probably related to their high sodium content, while the diets containing insect meal showed higher percentages of glutamic and aspartic acids, representative of the umami taste [88], than those of the HIM0 diet.

Conclusions
Data suggest that inclusion of Hermetia illucens meal positively influenced the hypocholesterolaemic/hypercholesterolaemic ratio and the content of indispensable amino acids and microelements. The microbiological quality of all diets testifies to the good practices of hygiene and sanitation applied during the production processes of fish feeds. The Esensing analysis permitted distinguishing color, odor and taste in the four feed groups. The combination of sensor responses (E-eye, E-nose and E-tongue) proves to be a powerful tool for discriminating different organoleptic profiles linked to different chemical compositions of experimental diets.
This study represents a part of a larger investigation aimed at evaluating the suitability of HIM addition in the Sparus aurata diet through the study of the productive performance, the chemical and organoleptic characteristics of the fillets and the possible development of intestinal inflammation.