Protein, Essential Amino Acid, and Fatty Acid Composition of Five Target Fishery Species of Central Mediterranean Sea

Simple Summary Fishery products are an important part of the Mediterranean diet and constitute a high source of nutrients for consumers, contributing to the prevention of several health diseases. This work analysed the protein, amino acids, and fatty acids of Mediterranean target species (four fish species and one crustacean species) which are mostly consumed in Sicily. The nutritional composition of the selected species together with the evaluation of nutrient indices showed the high quality of the sampled fish. The excellent profile of the fatty acid of “blue fish” species highlights the opportunity for local people to benefit from a highly nutritional product originating from local and sustainable fisheries and with an affordable commercial value. Abstract The protein, essential amino acid, and fatty acid composition of European pilchard (Sardina pilchardus), European hake (Merluccius merluccius), surmullet (Mullus surmuletus), red mullet (Mullus barbatus), and deep water rose shrimp (Parapenaeus longirostris) from the central Mediterranean Sea were investigated. All the species showed an essential amino acid content of about 50% of total amino acids, while the protein and total fatty acids content varied from 19.9 to 24.8% and from 1.4 to 5.1%, respectively. The fatty acid profile mainly followed the order SFA (39.1–52.6%) > PUFA (21.0–39.3%) > MUFA (15.6–24.3%). Palmitic and stearic acids were predominant among saturated fatty acids (38–52% and 21–25%, respectively), while palmitoleic and oleic acids were the most represented of the total monounsaturated acids (10–21% and 55–68%, respectively). All the species, as expected, showed a more significant proportion of n-3 PUFA (EPA + DHA) of about 81–93% of the total PUFA, with the highest values was found in European pilchard. Also, several fat quality index values, such as n-6/n-3 ratio, PUFA/SFA, the index of atherogenicity (IA), the index of thrombogenicity (IT), the hypocholesterolemic/hypercholesterolemic ratio (HH), and fish lipid quality/flesh lipid quality (FLQ) were calculated to assess the nutritional quality. All the obtained results, along with the fat quality indexes, indicated the excellent nutritional values of the selected species.


Introduction
Fishery products are largely considered healthy foods, representing the primary source of omega-3 polyunsaturated fatty acids (n-3 PUFAs) and contributing to the intake of numerous nutrients [1][2][3].This awareness/perception is deeply rooted in the central Mediterranean communities, where fishery products (bony fish, crustaceans, cephalopods, and elasmobranchs) are a regular part of the well-known "Mediterranean diet" [4].The protein and lipid contents of the edible part of fishery products make them nutritionally valuable and a valid substitute for other meals based on animal sources (such as red meat, pork, poultry, and seafood products from intensive aquaculture).
Several research studies have highlighted the importance of the effect of fish proteins in maintaining significant body functions and preventing health diseases related to metabolic syndromes, inflammation, or insulin resistance [5][6][7][8].Fish proteins are also rich in essential amino acids (EAAs) such as Valine, Leucine, Isoleucine, Threonine, Phenylalanine, Lysine, and Tryptophan [9,10].EAA intake is very important since they have a high biological value, regulating key metabolic pathways and being associated with the prevention and treatment of several organism's diseases and dysfunctions, such as obesity and diabetes, infertility, growth limitation, and intestinal, neurological, and renal dysfunctions [11][12][13].
PUFAs of the n-3 group, represented by eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA), are well known to exert several beneficial effects on human health, preventing cardiovascular, autoimmune, and neural diseases [14][15][16].Thus, seafood consumption has been demonstrated to exert health benefits in humans despite the potential risks due to contaminant bioaccumulation, especially in fish tissues [17][18][19], and/or the emerging threat of marine micro-and nanoplastic pollution [20].
Nevertheless, the higher consumption of processed food and vegetable oils in the modern Western diet results in a higher intake of n-6 PUFA [21,22].This higher intake of n-6 PUFA is known to be responsible for many health diseases [23,24], such as an increased incidence of neurodegenerative and neurological disorders [25].Therefore, increasing n-3 and reducing n-6 consumption to reduce the n-6/n-3 ratio between 5:1 and 10:1 is recommended for beneficial effects regarding human health [24,26,27].The fatty acid composition and content of marine fish can vary depending on various factors, such as diet, age, size, reproductive cycle, and environmental conditions [27].Thus, the knowledge of fatty acid content in different species, together with the evaluation of several nutritional indices, could be useful for consumers to identify highly healthy species, enabling them to make the best choices for their health.Some low-valued fishery products, as well as by-catches of massive fisheries [28,29], have the advantage of being, at the same time, highly nutritive food and easily accessible at low or no cost for people with low or zero income [11].This may be a strong tool to reduce worldwide malnutrition [30], providing at least 20% of protein intake for a third of the world's population that lives under poverty [31].Under this framework and based on landed production data provided by the DCF (Data Collection Framework; https: //datacollection.jrc.ec.europa.eu,accessed on 26 June 2020), we determine the nutritional composition (i.e., total protein and fatty acids) of five target species of the Italian commercial fishery from different markets in order to evaluate the advantages for human health and the differences between species coming from different sampling areas, consequently determining the most suitable for human consumption.

Samples Strategy and Collection
Five target species have been selected based on their total landing weight and/or related economic value.Specifically, one small pelagic species such as European pilchard (Sardina pilchardus), three demersal bony fish species such as European hake (Merluccius Merluccius), surmullet (Mullus surmuletus), and red mullet (Mullus barbatus), and one crustacean species such as deep water rose shrimp (Parapenaeus longirostris) were directly collected with the support of professional local fishermen from some major Sicilian fish landing ports (southern Italy) (Figure 1).Six ports were identified as follows: Porticello (PT) and S. Agata (SA) on the northern coast (the southern Tyrrhenian Sea), Catania (CT) on the east coast (Ionian Sea), and Mazara del Vallo (MV), Sciacca (SC), and Pozzallo (PZ) on the southern coast (Strait of Sicily Sea) (Figure 1), with the first aim of evaluate potential differences in the nutritional values of the selected species from different areas.After their collection, samples were packaged in a sealed polyethylene box and quickly frozen and stored at −20 • C. Once in the laboratory, biometric parameters, such as total length (cm) for fish and carapace length (cm) for crustaceans and weight (g), were measured for each specimen (Table 1).Specimens were dissected after thawing by means of stainless-steel scissors in order to avoid contamination.The edible part (muscle), cleaned by skin or carapace, was collected, and different pools (each composed between 1 and 30 organisms, depending on the species) were obtained.Finally, samples were homogenised for a few seconds with a knife mill (Retsch-grindomix GM 200) and set for biochemical analysis.
tacean species such as deep water rose shrimp (Parapenaeus longirostris) were directly collected with the support of professional local fishermen from some major Sicilian fish landing ports (southern Italy) (Figure 1).Six ports were identified as follows: Porticello (PT) and S. Agata (SA) on the northern coast (the southern Tyrrhenian Sea), Catania (CT) on the east coast (Ionian Sea), and Mazara del Vallo (MV), Sciacca (SC), and Pozzallo (PZ) on the southern coast (Strait of Sicily Sea) (Figure 1), with the first aim of evaluate potential differences in the nutritional values of the selected species from different areas.After their collection, samples were packaged in a sealed polyethylene box and quickly frozen and stored at −20 °C.Once in the laboratory, biometric parameters, such as total length (cm) for fish and carapace length (cm) for crustaceans and weight (g), were measured for each specimen (Table 1).Specimens were dissected after thawing by means of stainless-steel scissors in order to avoid contamination.The edible part (muscle), cleaned by skin or carapace, was collected, and different pools (each composed between 1 and 30 organisms, depending on the species) were obtained.Finally, samples were homogenised for a few seconds with a knife mill (Retsch-grindomix GM 200) and set for biochemical analysis.

Protein and Amino Acids Content
Considering that the protein content of fish flesh, unlike the fat content, is very constant, independent of seasonal variations and feeding and reproductive cycles, and that only small differences between species are present, protein measurements were performed on a singular pool of each species from different areas consisting of homogenised tissue [32].The protein determination was calculated by nitrogen values (%) determined by a Thermo Electron Flash EA 1112.The fresh samples (ca.0.5 mg) were packed into tin capsules and run against an internal standard (acetanilide with C = 71.09%,N = 10.36%).(*) total length from the tip of the rostrum to the posterior edge of the carapace.

Protein and Amino Acids Content
Considering that the protein content of fish flesh, unlike the fat content, is very constant, independent of seasonal variations and feeding and reproductive cycles, and that only small differences between species are present, protein measurements were performed on a singular pool of each species from different areas consisting of homogenised tissue [32].The protein determination was calculated by nitrogen values (%) determined by a Thermo Electron Flash EA 1112.The fresh samples (ca.0.5 mg) were packed into tin capsules and run against an internal standard (acetanilide with C = 71.09%,N = 10.36%).The crude protein was determined according to the Association of Official Analytical Chemists by converting the nitrogen content (N × 6.25) [33].Based on six analyses of the homogenised fish tissue, the repeatability test was carried out and the value was below 2%.The amino acid composition in fresh tissue was determined by ethyl chloroformate (ECF) derivatisation reagent [34] and using a GC/FID instrument (GC TRACE 1310 and Triplus 100LS Autosampler, purchased by Thermo Fisher Scientific, San Jose, CA, USA).Standard solutions were prepared to calibrate the instrument method and perform the calibration curves using the standard mixture 79248 (by Sigma Aldrich, Milan, Italy) composed of 20 amino acids.This standard was appropriately diluted with distilled water at several levels for each amino acid at 0.025 mM, 0.05 mM, 0.125 mM, 0.25 mM, and 0.5 mM, and the samples were treated as described below.About 10 mg of the freeze-dried samples were hydrolysed at 110 • C for 24 h with 200 µL of 6 M HCl in an autoclavable vial.After the vial was cooled to room temperature and the solution was dried under a gentle N 2 current flow, the residues were re-dissolved in 150 µL of distilled water and 150 µL of chloroform.After manual agitation using a glass capillary directly in the vial and waiting for the formation of two phases, the chloroform phase was discarded, and 100 µL of the water solution was transferred into a second micro-vial to derivatise the solubilised amino acids.This 100 µL of solution of each sample or 100 µL of each level of the amino acid standard solution was treated as follows: (i) mixed with 50 µL of a 4:1 solution of ethanol and pyridine; (ii) added 10 µL of ECF which was agitated slowly due to the released CO 2 gas, causing effervescence and consequently losing part of the sample; (iii) added 50 µL of a 1% ethyl chloroformate solution in chloroform; and (iv) 50 µL of the saturated NaHCO 3 solution was added to neutralise the solution.At the end of this procedure, each sample consisted of an aqueous top layer and a chloroform bottom layer.The lower organic layer was anhydrified by adding anhydrous sodium sulphate and was transferred to a fresh vial for the following GC/FID analysis.The amino acid was determined with a GC/FID Trace 1310 instrument (by Thermo Fisher Scientific) using a capillary column TR-1701 (Thermo Fisher Scientific) of 30 m × 0.25 mm × 0.25 µm and helium (He) carrier gas at a flow rate of 1.0 mL•min −1 ; the PTV injector was set in split mode (10:1) at 280 • C; and the FID detector temperature was 300 • C. The oven thermal gradient programme was set as follows: starting at 100 • C and held for 2 min, then ramped at 10 • C•min −1 to 270 • C, followed by a ramp, and held for 10 min.As for protein determination, the repeatability test for amino acid measurement was run, and it was between 8 and 12%.

Fatty Acids Profile Determination
The fatty acid (FA) profile was determined as fatty acid methyl esters (FAMEs) based on the lipid extraction reported by the Bligh and Dyer method [35].The fresh fish tissue of each pool per site was treated following these steps: about 1 g of the thawed sample was weighed directly in a glass vial; 10 mL of a chloroform/methanol solution (1:1 v/v) was pooled, capped, and placed in an ultrasound bath for 30 min; and in the same vial, another aliquot of 5 mL of chloroform was added and placed another time in an ultrasound bath.The solution containing the fish tissue fragments was filtered using a fast filter paper and pooled in a second vial.The residual was washed with another 5 mL of chloroform, and the paper filter was squeezed to quantitatively separate the lipid extracted from the tissue fragment and pool it in the previous second vial.After that, a chloroform clean-up procedure was followed, with about 5 mL of 1% saline water to remove methanol, polar substances, and co-extracted water-soluble impurities.After two layers had been formed, the upper water solution was discarded, and the chloroform solution containing lipid extract was anhydrified by adding a spoonful (about 0.5 g) of anhydrous sodium sulphate (Na 2 SO 4 ) and was then dried using a vacuum multivapor.Finally, 1 mL of hexane was added to re-dissolve the lipid extracted using a vortex mixer for 30 s.The FAMEs were obtained by a trans-esterification reaction: 0.2 mL of 2 M alcoholic potassium hydroxide (KOH) was added to the hexane solution directly into the previous vial, agitated by vortex for 30 s, and allowed at least 30 min to complete the reaction.The same vial was placed in the TriPlus autosampler of the GC/FID Trace 1310 instrument (Thermo Fisher Scientific).The fatty acid was determined with a GC/FID using a capillary column TR-FAME (Thermo Fisher Scientific) of 100 m × 0.25 mm × 0.1 µm and helium (He) carrier gas at a flow rate of 1.5 mL min −1 ; the PTV injector was set in split mode (20:1) at 280 • C; and the FID detector temperature was 300 • C. The oven thermal gradient programme was set as follows: starting at 100 • C (held for 2 min), then ramped at 2 • C min −1 to 235 • C (held for 5 min), followed by a ramp to 260 • C ramped at 40 • C (held for 5).The determination of FAMEs in the samples was identified by comparing the retention time for each fatty acid with the retention time measured using a Supelco MIX37 Component FAMEs mix standard (Supelco, Bellefonte, PA, USA).The instrument's calibration to quantify each fatty acid was developed using four levels of Mix37, as mentioned above, appropriately diluted with hexane.
The sum of each fatty acid determined the total FAs, and the profile was expressed as a percentage with respect to the total FAs.

Nutritional Quality Indexes
Different nutritional indexes based on fatty acid compounds were used to evaluate the lipid quality [36] of the selected species.The n-6/n-3 ratio was calculated as Σ n-6 PUFA/Σ n-3 PUFA; PUFA/SFA refers to the fraction of PUFAs over SFAs; the index of atherogenicity (IA) and thrombogenicity (IT), the Hypocholesterolemic/hypercholesterolemic ratio (HH), and the fish lipid quality (FLQ) were calculated as follows [37][38][39]: where MUFA stands for monounsaturated fatty acids; PUFA stands for polyunsaturated fatty acids; UFA stands for unsaturated fatty acids; and FA stands for fatty acids

Statistical Analysis
Significant differences in the total protein and essential amino acid percentages were tested among species using Kruskal-Wallis tests.A post hoc Dunn test was performed in the case of statistically significant results to determine which species differed from each other.
Nonmetric multidimensional scaling analysis (NMDS), based on the Euclidean distance matrix, was applied to categories of fatty acids (SFA, MUFA, PUFA n-3, PUFA n-6) to assess if the composition of the nutritional parameters were linked to particular species or fish markets.The final result was the transformation of the original multivariate dataset into a space with a reduced number of dimensions, visualised in an ordination diagram with two axes.The performance of the final model was measured through the "stress value" [40], which indicates the difference between the dissimilarities in the reduced dimension when compared to the complete multidimensional space.The correlation between the variables and axes was represented in the ordination plot by means of arrows with a length proportional to its strength.The closeness between points indicated how similar the fish individuals were among them based on fatty acid composition.
The calculated indices were plotted through boxplots by species and fish markets to assess possible different nutritional qualities.The significance of these differences was performed using Kruskal-Wallis tests.All statistical analyses were implemented using the statistical software R 4.1.0[41].

Protein Content and Essential Amino Acids
The average percentage of total proteins ranged from 19.9 to 24.8% of wet weight (w.w.) registered in European hake and European pilchard, respectively (Table 2).No statistically significant difference in amino acid percentage was found among species (KW = 4.81; p-value > 0.05), while the percentage of total proteins differed significantly between European hake and European pilchard (Dunn test = −3.36;p-value < 0.05) and between deep water rose shrimp and European pilchard (Dunn-test = −3.26;p-value < 0.05).The Daily Recommended Dietary Allowance (RDA) of proteins indicated for an adult (60 kg) by the World Health Organisation (WHO) is 0.8 g protein/kg body weight [42].Considering a 150 g portion of the selected species, the contribution to the recommended daily protein amount is 63% for European hake and deep water rose shrimp, 72% for red mullet and surmullet, and 78% for European pilchard.High-quality proteins are readily digestible and contain the dietary EAAs in amounts that correspond to human needs [42].We reported the content of EAAs measured in the selected species because it is the main determinant of the protein nutritional quality [43].The average fraction of EAAs in analysed seafood samples constitutes from 48 to 52% of the total amino acidic content.Among the analysed EAAs, leucine was the predominant fraction (2.2-3.2%w.w.).Our results indicate that all the selected species constitute a rich source of healthy protein and essential amino acids which is able to reach the recommended daily dose.
Among the monounsaturated fatty acids, the palmitoleic (C16:1 n-7) and oleic (18:1 n-9) acids were the most represented, accounting for 10-21% and 55-68% of total MUFA, respectively, as reported from various studies on marine species [49,50,52] (Tables S1, S2 and 3).The total PUFA varied from 21.0% to 39.3%, as registered in red mullet and European pilchard, respectively.Eicosapentaenoic (20:5 n-3, EPA) and docosahexaenoic (22:6 n-3, DHA) were the most abundant, accounting for 63-85% of the total PUFA.Generally, DHA assumed higher values than EPA in all the species analysed (Table 3).The highest value of EPA was detected in deep water rose shrimp (11.59% of FA) and the lowest in European hake (4.91%), whereas the highest content of DHA was found in European hake (26.28%) and the lowest in red mullet (7.57%).These results are in agreement with other works focused on the same species [49,53,54].The European Food Safety Authority (EFSA) suggests a daily intake of EPA and DHA consumption between 0.25 and 0.5 g in order to protect against the risk of cardiovascular diseases [55].The highest EPA and DHA content (mg g w.w.) was found in European pilchard, followed by deep water rose shrimp, European hake, surmullet, and red mullet (Figure 2, bottom).Based on data in Table 3, and considering a fish portion of 150 g, the EPA + DHA values were 0.48, 0.64, 0.67, 0.73, and 2.56 g per portion of red mullet, surmullet, European hake, deep water rose shrimp, and European pilchard, respectively.Figure 2, top, shows the proportion of the total fatty acid composition in the analysed species following the order SFA (39.1-52.6%)> PUFA (21-39.31%)> MUFA (15.63-24.34%).Only the red mullet showed a content of MUFA slightly higher than PUFA (24.34 and 21%, respectively), as found by Prato and Biandolino [46].Both n-6 and n-3 PUFA play crucial functions in regulating several biochemical and physiological processes with opposite effects [56].Since they are not synthesised by the organism, they must be assumed from the diet.Thus, an unbalanced fatty acid intake can result in the occurrence of several diseases [24,57].
The n-3 PUFA is a vital structural constituent of the phospholipid cell membranes participating in cell regulation and playing a crucial role in preventing many disorders such as cardiovascular, atherosclerosis, heart attack, depression, cancer, neural, and autoimmune diseases [13,15,24,[58][59][60][61][62][63][64][65].On the contrary, previous studies found that a higher intake of n-6 PUFA is often associated with adverse effects on human health, such as inflammatory processes, cancer, autoimmune, and cardiovascular disorders [66][67][68][69].All the analysed species contain a higher proportion of n-3 PUFA with respect to n-6 PUFA, representing 81-93% of the total PUFA.The highest content of n-3 PUFA was registered in European pilchard (36.27%) and the lowest in red mullet (17.37%).
The statistical approach using the NMDS model allowed for the identification of differences in fatty acid composition among the species and sampling areas.The stress value was smaller than 10%, indicating a good representation of the original data.All variables were significantly correlated with the axes.Scores from NMDS along the two first dimensions are graphically represented in Figure 3 Some species such as deep water rose shrimp, red mullet, and European hake are clustered together within the plot, indicating similar characteristics regardless of their origin.Differently, the surmullet and European pilchard are spread across multiple quadrants.Specifically, the group formed by surmullet from MV is located in a different quadrant (I quadrant) when compared to specimens from the other areas, therefore indicating distinct characteristics.Similarly, the group consisting of the European pilchard from CT diverges from the specimens of the same species from different areas.These results are probably due to both the intrinsic characteristics of the species and the availability of food in the area.Indeed, although marine fish could synthesise PUFAs (EPA and DHA) through characteristic enzymatic pathways, this enzymatic activity is low, and their content of PUFA appears to be closely related to diet [70,71].The primary production of PUFA can occur both in pelagic and benthic algae [72,73], and microalgae constitute an important source of these marine bioactive materials [74].Furthermore, oceanographic features play a key role in the distribution of primary producers, as well as in habitat characterisation and the availability of infauna, affecting the diet of marine species.The variability of the nutritional composition of fish depends on several biological parameters, such as different age, size, and the maturation status of the fish, and by different habitat characteristics, such as water temperature, food availability, and fishing impact [47,51,[75][76][77][78][79][80][81].The combination of these features easily explains the different fatty acid distributions in specimens of the same species coming from different fishing areas of the Mediterranean Sea.

Nutritional Quality Indices
Several fat quality indexes (Table 4), such as the n-6/n-3 ratio, PUFA/SFA, the index of atherogenicity (IA), the index of thrombogenicity (IT), the hypocholesterolemic/hypercholesterolemic ratio (HH), and fish lipid quality/flesh lipid quality (FLQ), can be useful to assess the contribution of FAME towards human health [37,38,[82][83][84].In this study, we found Some species such as deep water rose shrimp, red mullet, and European hake are clustered together within the plot, indicating similar characteristics regardless of their origin.Differently, the surmullet and European pilchard are spread across multiple quadrants.Specifically, the group formed by surmullet from MV is located in a different quadrant (I quadrant) when compared to specimens from the other areas, therefore indicating distinct characteristics.Similarly, the group consisting of the European pilchard from CT diverges from the specimens of the same species from different areas.These results are probably due to both the intrinsic characteristics of the species and the availability of food in the area.Indeed, although marine fish could synthesise PUFAs (EPA and DHA) through characteristic enzymatic pathways, this enzymatic activity is low, and their content of PUFA appears to be closely related to diet [70,71].The primary production of PUFA can occur both in pelagic and benthic algae [72,73], and microalgae constitute an important source of these marine bioactive materials [74].Furthermore, oceanographic features play a key role in the distribution of primary producers, as well as in habitat characterisation and the availability of infauna, affecting the diet of marine species.The variability of the nutritional composition of fish depends on several biological parameters, such as different age, size, and the maturation status of the fish, and by different habitat characteristics, such as water temperature, food availability, and fishing impact [47,51,[75][76][77][78][79][80][81].The combination of these features easily explains the different fatty acid distributions in specimens of the same species coming from different fishing areas of the Mediterranean Sea.

Nutritional Quality Indices
Several fat quality indexes (Table 4), such as the n-6/n-3 ratio, PUFA/SFA, the index of atherogenicity (IA), the index of thrombogenicity (IT), the hypocholesterolemic/ hypercholesterolemic ratio (HH), and fish lipid quality/flesh lipid quality (FLQ), can be useful to assess the contribution of FAME towards human health [37,38,[82][83][84].In this study, we found a n6/n3 ratio between 0.06 in European hake and 0.19 in red mullet, thus being lower than the value of 4.0 recommended by the UK Department of Health [84].This result indicates that the selected species contribute to a proper intake of these fatty acids (Table 4).The PUFA/SFA ratio reflects the effects of PUFA and SFA on human health [82].A minimum value ratio of 0.45 is recommended for one's diet since a higher value may result in increased blood levels of cholesterol [84].The PUFA/SFA ratio measured in the selected species was higher than the recommended value in all the species analysed except for red mullet (0.41), probably because of his lower PUFA content (21% of FAME) (Tables 2 and 3).However, it should be kept in mind that these indices may underestimate the quality of fatty acids, providing simplistic guidance.Generally, the ratio of PUFA/SFA or n-6/n-3 fails to reflect the effects of MUFAs [82].Therefore, additional indices capable of accounting for the functional effects of fatty acids are used.The indexes IA, IT, HH, and FLQ were calculated in the lipid fraction.IA, IT, and HH are associated with the potential platelet aggregation and subsequent adverse risks at the cardiovascular system [85].Lower IA and IT values are desirable to prevent cardiovascular disorder [37].In the selected species, the IA index varied from 0.58 for deep water rose shrimp to 0.92 for red mullet and the IT values ranged from 0.27 for European pilchard to 0.60 for red mullet (Table 3), thus remaining within the range of the expected values [27].On the contrary, a higher value of HH is more beneficial for human health [86].The measured HH index ranged from 1.27 (red mullet) to 3.37 (European pilchard), the latest exceeding values reported for marine fish (0.87-2.46) [27,87].The FLQ index is considered more suitable for fish or marine products due to their higher proportions of EPA and DHA [83].The obtained FLQ values fall in the range of variability measured for different marine species [83] (and references therein).The highest values were registered in European pilchard (33.8) and European hake (31.2) (Table 3), confirming the well-known value of these species relating to their higher n-3 PUFA content [49,51,88].

Conclusions
This study represents the first approach to characterise the nutrient parameters of most widely commercialised fish species in Sicily (Italy, Mediterranean Sea).The variability of fatty acid concentrations, both at the intra-and inter-specific levels, reflect the different biological characteristics of the analysed species and, presumably, the different conditions of their feeding areas.Although this study presents some limitations, such as the absence of minerals and vitamins analysis, our results demonstrate the high nutritional values of the selected seafood and especially highlight the healthy fatty acid nutrient profile of the "blue fish" species such as European pilchard, which shows the highest total protein amount, the highest n3-PUFA content, and the best nutritional quality indices.This study encourages the consumption of low-cost marine target species of local, artisanal, and sustainable fisheries which are characteristic of the study area, highlighting the opportunity for consumers to take advantage of high-nutrient resources with affordable commercial value.

Figure 1 .
Figure 1.Chart of the Mediterranean Sea indicating the Sicilian fish markets.

Figure 1 .
Figure 1.Chart of the Mediterranean Sea indicating the Sicilian fish markets.
. The samples are clearly clustered in the four quadrants of the plot.One cluster (I quadrant) is composed of European hake from all sites except PT, surmullet sampled in MV, and European pilchard from PZ.This cluster does not have a predominance of a particular fraction of fatty acids, but it is characterised by low MUFA and n-6 PUFA values.The second cluster (II quadrant), characterised by high SFA values, is composed of European hake from PT and surmullet and red mullet from SA.The third cluster (III quadrant), with high MUFA and n-6 PUFA values, is composed of red mullet from SC and deep water rose shrimp from CT and SC.The fourth cluster (IV quadrant) includes European pilchard individuals collected in CT, PT, and MV and deep water rose shrimp from PT and SA.This cluster is dominated by high levels of n-3 PUFA and low values of SFA, MUFA, and n-6 PUFA.Moreover, some samples are located along the negative side of the horizontal axis because they are dominated by SFA, MUFA, and n-6 PUFA, with low n-3 PUFA values.This group includes surmullet from CT and PT and red mullet from SC and PZ.Animals 2024, 14, x FOR PEER REVIEW 10 of 15

Figure 3 .
Figure 3. Ordination plot generated by NMDS model applied to fatty acid content in specimens from different fish markets.The points are represented with different colours and shapes, corresponding to a specific fish market and species, respectively.

Figure 3 .
Figure 3. Ordination plot generated by NMDS model applied to fatty acid content in specimens from different fish markets.The points are represented with different colours and shapes, corresponding to a specific fish market and species, respectively.

Table 1 .
Biometrics data, number of individuals per site, and habitat of the species.

Table 1 .
Biometrics data, number of individuals per site, and habitat of the species.

Table 2 .
Total protein content, total EAAs, and individual EAAs (% w.w.) of the selected species and relative sampling area.

Table 3 .
Mean ± standard deviation of total fatty acid content (% w.w.) and fatty acid composition (% of TFA) in the selected species from different sampling areas.

Table 4 .
Nutritional quality indexes of the lipid fraction of the selected species.