Abstract
Most commercial products available for Artemia sp. enrichment include high levels of docosahexaenoic acid (DHA) to boost its nutritional value, though with limited success. In this regard, the present study evaluated the alternative utilization of eicosapentaenoic acid (20 : 5n-3; EPA) oils to improve the n-3 long-chain highly unsaturated fatty acid content (n-3 LC-HUFA) in enriched Artemia sp. to feed greater amberjack (Seriola dumerili) larvae. Five experimental emulsions containing increasing levels of EPA from 0.8% to 60% of total fatty acids (TFA) and n-3 LC-HUFA (1.3%–70.6% TFA) were formulated. Each diet was fed to greater amberjack larvae (17–35 days posthatch (dph)) in three replicate 200-L tanks. The dietary EPA supplementation significantly improved larval growth during the feeding trial and survival at 35 dph (). In addition, larval fatty acid profiles showed a positive correlation with dietary EPA. Finally, despite the sum of total skeletal anomalies and column anomalies were not significantly affected by dietary EPA, increasing dietary EPA, and n-3 LC-HUFA tended to reduce the incidence of these types of anomalies in greater amberjack larvae at 35 dph. Based on these results, S. dumerili larvae could be successfully grown with low DHA but high EPA-rich oil enrichment products when n-3 LC-HUFA content in Artemia sp. is maintained in sufficient amounts.
1. Introduction
The greater amberjack Seriola dumerili (Risso, 1810) is an excellent candidate for marine fish species aquaculture diversification due to its high consumer acceptability worldwide, high market value, and excellent flesh quality [1–3]. Although the aquaculture of Seriola sp. species has increased in recent years from 160,000 metric tons in 2015 to 167,000 metric tons in 2019 [4], their production is still constrained by the reduced number of commercial hatcheries and juveniles availability. Seriola juveniles for on-growing operations are mostly collected from the wild in Asian countries, while are produced in hatcheries in Europe, Australia, and America. In this regard, the limited knowledge of zootechnics and nutritional requirements during larval stages is one of the major causes for the poor hatchery survival and juvenile quality obtained under aquaculture sustainable methodologies. Up to date, the most common feeding protocol for greater amberjack larvae produced in hatcheries is based on the use of enriched rotifers (Brachionus sp.) from the first feeding, followed by the use of brine shrimp (Artemia sp. nauplii and meta nauplii) appropriately enriched with commercial or experimental products to finally cofeed during the weaning stages until complete replacement of live feed at later stages [5–11]. Thus, hatchery-reared juveniles are usually characterized by the high number of weak fish with low resistance to manipulation at early stages [10], high size variability, and high incidence of skeleton anomalies in comparison to the wild-captured ones [11]. Although some attempts are still underway to introduce the use of marine copepods as food for the first feeding of greater amberjack, the complexity and costs of copepod production limit this approach to a reduced number of commercial and experimental stations with specific knowledge and capacity to produce copepods at a large scale. Nevertheless, with an adequate formulation of the enrichment products and the knowledge and quantification of the species-specific requirements, standard live preys (rotifers and Artemia sp.) can be used as effectively as copepods for greater amberjack culture. One of the most important nutritional differences and benefits of copepods in comparison to rotifers and Artemia sp. is related to their optimal lipid classes and fatty acid profiles, which usually cover the nutritional needs of most marine fish larvae [12, 13]. In this regard, it is well known that one of the most important factors affecting the performance of marine fish larvae is the n-3 long-chain highly unsaturated fatty acid content (n-3 LC-HUFA) composition of their feed [14]. Thus, depending on the enzymatic capacity of the different species for desaturation and elongation to convert the C:18 polyunsaturated fatty acids, α-linolenic acid (18 : 3n–3), into the biologically active long-chain (C:20–24) polyunsaturated fatty acids, they must be included or not in the larval diets to fulfill their nutritional requirements [14]. Specifically, docosahexaenoic acid (DHA, 22 : 6n-3) and eicosapentaenoic acid (EPA, 20 : 5n-3) are essential components of cellular membranes, modulate physiological processes, including membrane transport, receptors, and enzymatic activities [15], and play essential roles in fish larval development in many marine species [10, 16–18]. However, to date, a limited number of studies have addressed the fatty acids requirements for greater amberjack during live feed stages, mostly focusing on DHA enrichment [8, 10]. As stated above, the feeding sequence of S. dumerili larvae usually includes Artemia sp. nauplii, enriched or not, from as early stages as 10 days posthatch (dph) [10, 11]. However, the fatty acid composition of Artemia sp., even after enrichment, is usually deficient in n-3 LC-HUFA, particularly in DHA [19, 20]. Major causes of this effect are associated with enrichment protocols themselves, that is, accelerating autoxidation of n-3 LC-HUFA due to high water temperatures (28°C) and vigorous aeration [21], while others are intrinsic to Artemia sp. metabolism, such as undesired metabolic retroconversion of DHA to EPA, its mobilization to neutral lipid fraction and major accumulation as triglycerides [20, 22].
In this regard, most of the commercial products available for Artemia sp. enrichment include high levels of DHA, trying to boost n-3 LC-HUFA and particularly DHA in Artemia sp. while maintaining a high (>2) DHA/EPA ratio. With this perspective, the objective of the present study was to determine the potential utilization of purified EPA oils as a substitute for DHA during Artemia sp. enrichment to cover the nutritional requirements of n-3 LC-HUFA in S. dumerili larvae, evaluating its effects on survival, growth, stress resistance/tolerance, skeletal anomalies development, and tissue composition.
2. Materials and Methods
2.1. Ethic Statement
All experimental protocols described in this manuscript fully comply with the recommendations in the Guide for Care and Use of Laboratory Animals of the European Union Council (2010/63/EU) and Spanish National Legislation (Spanish Royal Decree 53/2013) which regulate animal usage in experimentation and/or other scientific purposes, and with the general guidelines approved in the EU-project DIVERSIFY (KBBE.2013.1.2–09 KBBE.2013.1.2–09, grant agreement no. 603121).
2.2. Larval Rearing
Greater amberjack larvae of 17 dph, (n = 1,000 per tank; total length (TL) 6.39 ± 0.44 mm; fresh weight 2.94 ± 0.57 mg; mean ± standard deviation) were randomly distributed into 15 light gray, cylindrical fiberglass tanks (five triplicate treatments) of 200 L. All tanks were equipped with continuous aeration and supplied with filtered UV-sterilized seawater (salinity 37 psu) in a flow-through system. Water exchange was gradually increased from 0.80 L min−1 at 17 dph to 1.60 L min−1 at 25 dph and finally to 3.30 L min−1 from 30 dph to the end of the feeding trial at 35 dph, ensuring good water quality and removal of nonfeeding preys from 1 day to the next. The average seawater temperature and dissolved oxygen during the experimental period were 24.15 ± 0.35°C and 6.55 ± 0.41 ppm, respectively. An artificial fluorescent light above each tank provided a surface light intensity ranging between 1,000 and 1,500 lux (Digital Lux Tester YF-1065; Powertech Rentals, Osborne, Australia) at the center of each rearing tank, with a photoperiod of 12 hr light/12 hr darkness cycle with lights on from 08:30 a.m. to 8:30 p.m. From 17 to 22 dph, there was an overlap between rotifers (nonenrich) and enriched Artemia sp. with a gradual reduction in the amounts of nonenriched rotifers from 5,000 to 0 individuals L−1 and a progressive increasing density of enriched Artemia sp. (125–500 individuals L−1) provided two times per day (9:00–15:00 hr) during the same period. From 23 to 35 dph, larvae were exclusively fed with 500 individuals L−1 of enriched Artemia sp. two times per day (9:00–15:00 hr), which were offered following complete clearing of the tank between feedings.
2.3. Growth and Survival
All larval samplings were carried out randomly from each experimental tank 1 hr before the first feed at 17 dph and after night starvation at 35 dph. Larvae were euthanized with an overdose of clove oil containing 87% eugenol (Guinama®, Valencia, Spain) and TL was measured with a profile projector (Mitutoyo PJ-A3000, Kanagawa, Japan). Additionally, the fresh body weight (FW) of 30 larvae/tank was determined initially (17 dph) and finally at 35 dph. Larval survival was determined by hand counting of all the remaining alive larvae in each tank at the end of the experiment at 35 dph.
2.4. Stress Test
At 35 dph, 30 larvae per tank were used to conduct an acute air stress test, by handling them out of the water in a scoop net for 30 s. This test was repeated with another 30 larvae, but increasing the exposure time to 60 s. After the test, each set of larvae was returned to separate buckets with aerated seawater to determine survival after 24 hr [10, 23].
2.5. Experimental Enrichments
2.5.1. Preparation
To determine the capacity of EPA inclusion in enrichments to raise the n-3 LC-HUFA content in Artemia sp., five experimental enrichments (EPA 0–4) with increasing EPA content (0%–60%) were formulated (Table 1). The enrichments were prepared by mixing increasing amounts of a commercial oil (Incromega EPA 500TG; Croda, Barcelona, Spain) containing 63% EPA and 8% DHA in total fatty acids (TFA) with a second commercial oil (Oleic Oil; Sigma–Aldrich, Madrid, Spain) containing 77% oleic acid (OA, 18 : 1n-9) in TFA. Additionally, every enrichment preparation was mixed with constant amounts of soya lecithin (Korot SL, Alcoy, Spain) containing 54% linoleic acid (18 : 2n-6) and trace amounts of EPA and DHA as emulsifier. Also, enrichment preparations were fortified with 3 g kg−1 vitamin E (DL-α-tocopherol acetate; Sigma–Aldrich, Madrid, Spain) and 2.5 g kg−1 vitamin C (L-ascorbic acid, Asc; Sigma–Aldrich, Madrid, Spain) to prevent the oxidation according to Betancor et al. [24], Hamre et al. [13], and Roo et al. [10]. Once prepared, experimental enrichments were stored in a fridge at 4°C until used. During the feeding trial, three subsamples of each experimental enrichment (ca. 100 g), newly hatched Artemia sp. (ca. 100 g), and enriched Artemia sp. (ca. 100 g) were collected and analyzed for proximate and fatty acid composition (Table 2).
2.5.2. Artemia sp. Enrichment
Brine shrimp Artemia sp. cysts (Sept Art type; INVE Aquaculture, Dendermonde, Belgium) were incubated for 24 hr in 150 L cylindrical–conical tanks with highly aerated seawater (37 g L−1) at 28°C. Newly hatched nauplii were harvested, washed, and transferred to a 10 L bucket (one bucket for each experimental emulsion) filled with clean seawater and provided with moderated aeration, oxygen supply automatically controlled to maintain 6.0 ± 0.2 ppm and 28°C temperature. To enrich the Artemia sp., 1.1 ml of each experimental enrichment was stirred for 1 min with 100 ml of fresh water and added to each biker. The same enrichment conditions were utilized for all the experimental enrichments assayed, including enrichment time (18 hr), Artemia sp. density (500,000 individuals L−1) and amount of experimental enrichment.
2.6. Lipid Content and Fatty Acid Analysis
A sample of larvae (ca. 50 individuals) from each experimental tank at 35 dph was collected for biochemical analysis. All samples were flushed with N2 and kept frozen at −80°C until analysis was carried out. Total lipids were extracted in chloroform : methanol (2 : 1, v : v) using the method of Folch et al. [25]. Fatty acids were prepared by transesterification using the method of Christie and Han [26]. Separation and identification of the fatty acids was realized with gas chromatography (GC Thermo Finnigan Fucus GC, Milan, Italy) under the conditions reported by Izquierdo et al. [27]. Dry matter, lipids, and ash content were determined using the methods of analysis of the Association of Official Analytical Chemists [28]. Results of total lipids, TFA content in the enrichment emulsions, enriched Artemia sp. and larvae for each treatment are shown in Tables 1–3, respectively. Data are expressed as percentage of total fatty acids (% TFA).
2.7. Skeletal Anomalies
At 35 dph, 100 larvae per tank were euthanized with clove oil and fixed in 10% buffered formaldehyde, cleared, stained with alizarin red, and immediately photographed to evaluate skeletal anomalies occurrence [29]. Anomalies were studied in the different axial column regions according to the methods described by Boglione et al. [30].
2.8. Statistics
All the data were statistically analyzed using SPSS Statistical Software System 15.0 (SPSS, http://www.spss.com). Data were presented as mean + standard deviation. Effects of the dietary levels of EPA were analyzed by regression and analysis of variance at a significant level of 5%. All variables were checked for normality and homogeneity of variance using the Kolmogorov–Smirnoff and the Levene tests, respectively [31].
The nonparametric Kruskal–Wallis test was accomplished in the analysis of survival after stress test data as a result of heterogeneous variance. To evaluate the differences in skeletal anomalies, logarithmic linear statistical analysis was performed [31].
3. Results
3.1. Experimental Emulsions and Artemia sp. Enrichment
Formulations and FA compositions of the experimental emulsions (EPA 0–4) are listed in Table 1. The use of OA oil to compensate for the reduction of EPA-rich oil (EPA 500TG) resulted in variations in saturated, monoenoic, total n-3, total n-3 LC-HUFA, total n-9, and different FA ratios such as DHA/EPA between the experimental emulsions. The lowest level of EPA in the experimental emulsions was present in EPA-0 with 0.8% in TFA, while the EPA content in EPA 1–4 increased from 14.2% to 60.1% in TFA. The DHA and arachidonic acid (ARA) in the experimental emulsions increased in accordance with EPA, from 0.3% to 6.7% and from 0% to 3.5% in TFA, respectively. Artemia sp. enriched with different experimental emulsions resulted in five enriched Artemia sp. treatments (Table 2). Thus, EPA and n-3 LC-HUFA content in Artemia sp. mirrored EPA and n-3 LC-HUFA content in the experimental emulsion, ranging from 1.08% to 22.9% and from 3.2% to 29.0%, respectively. Other fatty acids such as DHA (ranged from 0.14% to 3.01% TFA) and ARA (ranged from 0.39% to 1.72% TFA) were also correlated with their emulsion content. On the other hand, the DHA/EPA ratio was kept quite constant (average 0.1) for all the experimental treatments (Table 2).
3.2. Larval Growth and Survival
Larval growth was significantly improved with the increase of dietary EPA and n-3 LC-HUFA during the feeding trial (Figures 1(a) and 1(b)). Mean values for all the dietary treatments were 11.72 ± 0.47 mm in TL and 27.13 ± 3.21 mg in FW at 35 dph. Larval fed EPA-2 and EPA-3 showed significant () better growth in terms of TL than those larvae fed EPA-0, but similar to EPA-1 and EPA-4 (Figure 1(a)). Growth in terms of fresh weight followed a similar pattern to that in TL. Thus, larvae fed EPA-3 showed significantly () higher FW than EPA-0 and EPA-1 but like EPA-2, and EPA-4 () (Figure 1(b)). The relationship between final TL and Artemia sp. EPA content was described by the second-order polynomial regression: y = 10.893 + 0.214x − 0.080 x2 (y = Total length; x = Dietary EPA (% TFA); R2 = 0.965; ). Also, the final fresh weight was described by the equation: y = 21.503 + 1.460x – 0.0544 x2 (y = Fresh weight; x = Dietary EPA (% TFA); R2 = 0.989; ). Under the experimental conditions applied, both growth models suggest that maximum growth would be achieved in the range of dietary EPA and n-3 LC-HUFA concentrations tested, between 8%–16% TFA and 11%–21% TFA, respectively. Specifically, the highest growth values would be achieved at 13.3% EPA and 17.7% n-3 LC-HUFA content in Artemia sp. Larval survival was significantly () improved by dietary EPA at 35 dph, the lowest survival was recorded in those larvae receiving the lowest EPA and n-3 LC-HUFA in the Artemia sp. (EPA-0). Hence, larval resistance to stress test (30 s), determined as the survival rate 24 hr after handling, was significantly affected by dietary EPA and n-3 LC-HUFA at 35 dph, being the lowest survival recorded in those larvae fed the lowest EPA and n-3 LC-HUFA in the diet (EPA-0) (Figures 2(a) and 2(b)), but no survivors were recorded at 60 s stress test in any treatment assayed.

(a)

(b)

(a)

(b)
3.3. Proximate and Fatty Acid Composition in S. dumerili Larvae
The results of the total lipid and fatty acid composition of S. dumerili larvae at 35 dph are shown in Table 3. Total lipids, ash content, and moisture were similar for the larvae fed any of the Artemia sp. treatments. Regarding the larvae FA composition, n-3 LC-HUFA were significantly () affected by increased levels of dietary EPA in EPA-0 and EPA-1 treatment in comparison to EPA-4. Similarly, the EPA level in larvae fed EPA-0 was significantly () different than EPA-4 fed ones. DHA contents in larval tissues were also affected by dietary treatments, with the lowest values for those larvae fed EPA-0, which was significantly () lower than those larvae fed EPA-4. On the other hand, ARA content or DHA/EPA ratios in larval tissues were not affected () by the different treatments.
3.4. Larval Quality
The incidence of total acute skeleton anomalies was not significantly affected by dietary EPA and n-3 LC-HUFA with an average value for all dietary treatments of 22.63% ± 6.25% (Figure 3(a)). Regardless of the EPA and n-3 LC-HUFA content, most of the acute skeleton anomalies were major alterations in the hemal region, related to column anomalies, with scoliosis being the most important recorded affection (Figure 3(b)). The relationship between total acute skeleton anomalies and Artemia sp. EPA content was described by the second-order polynomial regression: y = 32.269 + (−1.654x) + 0.043x2 (y = Total anomalies; x = Dietary EPA (% TFA); R2 = 0.9854; ) (Figure 3(a)). Similarly, the relationship between the sum of column anomalies and Artemia sp. EPA content was described by the same model, represented as: y = 32.163 + (−1.658x) + 0.043x2 (y = Total column anomalies; x = Dietary EPA (% TFA); R2 = 0.9854; ) (Figure 3(b)), suggesting that increasing dietary EPA and n-3 LC-HUFA tends to reduce the incidence of these anomalies during the live stage of development evaluated. Under the experimental conditions and the regression models applied, these results suggest that the best larval quality in terms of skeleton anomalies occurrence was achieved at EPA levels between 16% and 21% TFA with a maximum around 17% EPA content in Artemia sp., which matched with the data obtained to reach the best growth.

(a)

(b)
4. Discussion
Present results demonstrated successful EPA and n-3 LC-HUFA enrichment of Artemia sp. nauplii with the experimental emulsions formulated. Gradual increase in EPA and n-3 LC-HUFA levels in Artemia sp. improved greater amberjack larval growth. Thus, the highest growth was obtained when larvae were fed Artemia sp. containing n-3 LC-HUFA in a range of 11%–21% TFA (2.3–5.3 g 100 g−1n-3 LC-HUFA dry weight (dw)), with EPA content between 8% and 16% TFA (0.8–2.2 g 100 g−1 EPA dw) despite a low DHA level (0.87%–2.22% TFA; 0.1–0.5 g 100 g−1 DHA dw) and a low DHA/EPA ratio (average 0.1). Furthermore, the optimal n-3 LC-HUFA levels based on EPA provision in this study is consistent with the optimum dietary n-3 LC-HUFA levels determined to obtain maximum growth and survival in a previous study conducted by our team but using rich DHA oils as the main n-3 LC-HUFA source [10] These data also agreed with n-3 LC-HUFA levels reported by Takeuchi et al. [32] and Ishizaki et al. [33] in other carangid species such as striped jack (Pseudocaranx dentex) and yellowtail (Seriola quinqueradiata). These authors did not find significant differences in either growth or survival when enriched Artemia sp. was used with both DHA-rich oil or EPA ones with n-3 LC-HUFA levels in the range of 11%–25% TFA and 13%–21% TFA, respectively. Hence, Matsunari et al. [7] reported better growth and survival in younger amberjack larvae (3–10 dph) when n-3 LC-HUFAs are in the range of 15%–22% TFA but using DHA-rich oils. In the same way, present data showed a reduction in greater amberjack larval growth and survival when fed on Artemia sp. enriched with n-3 LC-HUFA content lower than 11% TFA. In the aforementioned trials with greater amberjack, yellowtail, and striped jack, all of those species performed significantly better when fed live preys containing n-3 LC-HUFA levels over this value, which could be considered as a minimum amount to be raised during live prey enrichment for greater amberjack larval production [10]. Additionally, the negative effect on larval growth when the level of n-3 LC-HUFA is higher than 21% might be an indication of an excess of n-3 LC-HUFA beyond larval requirements. Conversely, increased levels of DHA and EPA and DHA/EPA ratios in live preys are reported as nutritional indicators to promote marine fish larval growth, survival, stress resistance, pigmentation, and larval quality [15, 34, 35]. Thus, target values of DHA/EPA over 1.0 and even close to 2.0 or higher are usually recommended during live prey enrichment, taking as reference values reported in natural live preys such as copepods [12]. However, the intrinsic capacity of Artemia sp. to metabolize DHA during the enrichment and postenrichment process, led to the retroconversion of DHA to EPA [19, 22], thus reducing the DHA/EPA ratio and the nutritional value of this live prey. In this study, the utilization of EPA-rich oil, as a main source of n-3 LC-HUFA, despite elevating the n-3 LC-HUFA content in Artemia sp., led to a reduction in DHA/EPA ratio to values close to 0.1, but with an efficient utilization of the fatty acids provided in the experimental emulsion, with no apparently harmful effects on larval growth, survival, resistance to stress test or quality when sufficient n-3 LC-HUFA were provided. Most published studies report that DHA is more limiting for growth, survival, and nervous system and schooling behavior development than EPA [15, 36]. On the contrary, other studies identify high levels of dietary DHA as a cause of muscular dystrophy [37] or spinal anomalies in larvae [38] due to the high peroxidation risk of DHA and the formation of toxic oxidized compounds, if they are not properly protected with antioxidant substances [18, 35].
In this regard, the relationship between dietary n-3 LC-HUFA and its role in bone growth and development has been previously addressed [17, 35]. Present data showed a 10% reduction in total acute and sum of column skeleton anomalies when n-3 LC-HUFA content was increased from 3% to 11% TFA, suggesting the beneficial effect of n-3 LC-HUFA supplementation and confirming the appropriated FA antioxidant protection during experimental enrichment formulation.
5. Conclusion
Based on these results, the utilization of EPA-rich oils to increase n-3 LC-HUFA content in Artemia sp. seems to be an effective option to improve Artemia sp. nutritional quality. Additionally, the utilization of this experimental enrichment with a n-3 LC-HUFA content between 11% and 21% TFA and current methodology avoids Artemia sp. nauplii retroconversion of DHA to EPA, providing a sufficient amount of n-3 LC-HUFA for greater amberjack larvae during Artemia sp. feeding stage, with no harmful effect over larval growth, survival, and stress resistance.
Data Availability
Data used to support the findings of this study are available upon request.
Conflicts of Interest
The authors declare that they have no known competing financial interests or personal relationships that could have influenced the work reported in this paper.
Authors’ Contributions
Javier Roo did the experimental conceptualization, data curation, formal analysis, original draft writing, funding acquisition, and project administration. Daniel Montero draft review and editing. Quirós-Pozo Raquel did the data curation, formal analysis, and original draft revision. Christian Monzón-Rivero did the data curation, formal analysis. Marisol Izquierdo López funding acquisition, project administration, supervision, draft review, and editing.
Acknowledgments
We thank Mr. Ruben Rodriguez for his support with live prey production and larval rearing. This study was undertaken within the project “Exploring the biological and socioeconomic potential of new/emerging candidate fish species for expansion of the European aquaculture industry (DIVERSIFY; GA603121)” funded by the European Commission. Lab analyses were co-founded by the THINKINGAZUL program supported by the Spanish Ministry of Science and Innovation with funds from the European Union NextGenerationEU (PRTR-C17.I1). Additionally, publication fees have been covered by the project ULPGC EXCELLENCE, funded by the Ministry of Universities, Science, Innovation and Culture of the Canary Islands Government.
Supplementary Materials
Supplementary 1. Graphic summary that synthesizes the contents of the main results of the study in an illustrated and concise way.
Supplementary 2. Larval rearing sequence and sampling points utilized in the study.