Optimizing the Timeframe to Produce Polychaetes (Hediste diversicolor) Enriched With Essential Fatty Acids Under Different Combinations of Temperature and Salinity

Polychaetes can be successfully employed to recover essential fatty acids (EFA) from wasted uneaten aquafeeds present in aquaculture effluents. The optimization of the timeframe required to produce premium ragworms (Hediste diversicolor) biomass rich in EFA is paramount to make available to the aquafeeds industry another alternative ingredient to fish meal and fish oil. The present study aimed to evaluate the potential enrichment of ragworms fatty acid (FA) profile when fed a commercial aquafeed during 10, 20, and 40 days (D10, D20, and D40) under different combinations of water temperature (20 and 25°C) and salinity (15, 20, and 25). Total FA incremented progressively overtime, with D40 polychaetes exhibiting average values ranging between 70 and 90 μg mg–1 DW. The average values of n-6 FA ranged between 13 and 17 μg mg–1 DW, while that of n-3 FA varied between 17 and 19 μg mg–1 DW at D40. No significant differences were found in the FA profile of H. diversicolor cultured under different combinations of temperature and salinity. The FA profile of cultured polychaetes exhibited between 28 and 31% dissimilarity from that of wild conspecifics and displayed a higher content of two essential n-3 FA: eicosapentaenoic (20:5 n-3, EPA) and docosahexaenoic acids (22:6 n-3, DHA) (values ranging between 9.6–11.2% and 4.3–5.0% of total FA, respectively). A higher similarity in FA profile was recorded between D40 polychaetes and aquafeed than with initially stocked or wild specimens. Palmitic (16:0), oleic (18:1 n-9), linoleic (18:2 n-6), eicosadienoic (20:2 n-6), EPA (20:5 n-3), and DHA (22:6 n-3) were the FA whose concentration exhibited the highest increment. Evidence of de novo FA biosynthesis was observed through the formation of some FA that were neither present in the initially stocked biomass, nor in the aquafeed supplied, such as 5,11-eicosadienoate (Δ 5,1120:2), 7,13,16-docosatrienoate (Δ 7,13,1622:3), dihomo-gamma-linolenic (20:3 n-6), eicosatrienoic (20:3 n-3) and eicosatetraenoic (20:4 n-3) acids. A plateau of total FA, n-6, and n-3 FA was not reached over the study period. Overall, the present study highlights the potential of H. diversicolor as an extractive species for integrated multi-trophic aquaculture (IMTA) applications.

Polychaetes can be successfully employed to recover essential fatty acids (EFA) from wasted uneaten aquafeeds present in aquaculture effluents. The optimization of the timeframe required to produce premium ragworms (Hediste diversicolor) biomass rich in EFA is paramount to make available to the aquafeeds industry another alternative ingredient to fish meal and fish oil. The present study aimed to evaluate the potential enrichment of ragworms fatty acid (FA) profile when fed a commercial aquafeed during 10, 20, and 40 days (D10, D20, and D40) under different combinations of water temperature (20 and 25 • C) and salinity (15, 20, and 25). Total FA incremented progressively overtime, with D40 polychaetes exhibiting average values ranging between 70 and 90 µg mg −1 DW. The average values of n-6 FA ranged between 13 and 17 µg mg −1 DW, while that of n-3 FA varied between 17 and 19 µg mg −1 DW at D40. No significant differences were found in the FA profile of H. diversicolor cultured under different combinations of temperature and salinity. The FA profile of cultured polychaetes exhibited between 28 and 31% dissimilarity from that of wild conspecifics and displayed a higher content of two essential n-3 FA: eicosapentaenoic (20:5 n-3, EPA) and docosahexaenoic acids (22:6 n-3, DHA) (values ranging between 9.6-11.2% and 4.3-5.0% of total FA, respectively). A higher similarity in FA profile was recorded between D40 polychaetes and aquafeed than with initially stocked or wild specimens. Palmitic (16:0), oleic (18:1 n-9), linoleic (18:2 n-6), eicosadienoic (20:2 n-6), EPA (20:5 n-3), and DHA (22:6 n-3) were the FA whose concentration exhibited the highest increment. Evidence of de novo FA biosynthesis was observed through the formation of some FA

INTRODUCTION
World aquaculture production reached an all-time record of 114.5 million tons in live weight in 2018 (USD 263.6 billion) (FAO, 2020). This figure includes the production of 54.2 million tons of fish (USD 139.7 billion) and 9.4 million tons of crustaceans (USD 69.3 billion), with 13 and 61% of these productions resulting from saltwater aquaculture (which includes marine and brackish water species), respectively (FAO, 2020). The majority of farmed species of both groups are produced using commercial aquafeeds (∼ 57 million tons in 2018) (FAO, 2020). Most marine teleosts (and other vertebrate species including humans) display a limited ability to perform de novo synthesis of n-3 highly unsaturated fatty acids (HUFA) due to the lack of desaturases 12 and 15, the enzymes responsible to produce polyunsaturated fatty acids (PUFA) from oleic acid (18:1 n-9) (Hastings et al., 2001(Hastings et al., , 2004Bell and Tocher, 2009;Monroig and Kabeya, 2018). As such, it is paramount to include in aquafeeds a balanced profile of HUFA and not their FA precursors (Sargent et al., 2003;Tocher, 2015Tocher, , 2003. Aquaculture of marine fishes and crustaceans thus depends on a rich supply of n-3 HUFA, including eicosapentaenoic (20:5 n-3, EPA) and docosahexaenoic (22:6 n-3, DHA) acids, through the formulation and supply of well-balanced aquafeeds commonly employing fishmeal and fish oil, two increasingly scarcer and costly resources for aquafeeds Metian, 2015, 2008;Cashion et al., 2017;FAO, 2020). The progressive reduction in the supply of these marine origin raw materials has been coupled with a surging demand driven by a fastgrowing aquaculture industry (farmed aquatic animals grew on average 5.3% per year between 2001 and 2018) (FAO, 2020). There is a global requirement of approximately 0.4 million metric tons of n-3 HUFA per year (Tocher, 2015). Formulated aquafeeds have to both satisfy the nutritional needs of cultured species and safeguard that at the end of a productive cycle farmed species display an optimal biochemical profile for human nutrition.
Polychaetes display a significant commercial value for the culture of vertebrates and invertebrates (García-Alonso et al., 2008;Olive, 1999). According to reported values, approximately 121.000 tons of polychaetes (with an estimated value of USD 8.39 billion) were harvested globally in 2015, with these figures being comparable to some of the most important world fisheries (Watson et al., 2016). It has already been acknowledged that the collection of polychaetes from the wild is likely insufficient to satisfy global market demands and that this practice drives a multitude of negative environmental impacts (Fidalgo e Costa et al., 2006;Pombo et al., 2018;Xenarios et al., 2018). The polychaete Hediste diversicolor O.F. Müller, 1776, popularly known as ragworm, is a highly valued bait species for sports fishing (Fidalgo e Costa et al., 2006;Carvalho et al., 2013;Sá et al., 2017;Pombo et al., 2018). It also plays a key dietary role on the nutrition and production of some fish and crustacean species (e.g., soles, shrimps and crabs), being often used to trigger gonad maturation and spawning (Lytle et al., 1990;Dinis, 1992;Luis and Ponte, 1993;Dinis et al., 1999;Meunpol et al., 2005). This omnivorous species is an active predator (Fauchald and Jumars, 1979) that is also able to perform deposit-feeding (Reise, 1979;Scaps, 2002) and filter-feeding (Riisgård, 1994;Vedel et al., 1994). Their deposit-feeding behavior allows them to consume considerable amounts of organic matter present in the substrate where it burrows (Reise, 1979;Scaps, 2002). The biomass of ragworms contains essential ingredients that can act as important supplements for aquafeeds formulation, such as amino acids and other odorants that elicit a feeding response for several fish species (e.g., Senegalense sole) (Velez et al., 2007;Wang et al., 2019). Ragworms also display high proteic and lipidic contents (54-60% and 11-22% of DW biomass, respectively) (Wang et al., 2019;Yousefi-Garakouei et al., 2019). This species has already been shown to perform de novo biosynthesis of some essential fatty acids (EFA) from acetyl coenzyme A, by using several fatty acids (FA) desaturase and enlogase enzymes; as such, it is usual to detect higher concentrations of PUFA and HUFA in ragworms biomass than on their diet (Vásquez et al., 2014;Santos et al., 2016). The abundance of EFA and ragworms ability to recycle n-3 HUFA (such as EPA and DHA) from uneaten aquafeeds and sludge from aquaculture effluents, that would otherwise be lost to the environment, make H. diversicolor a key extractive species for integrated multi-trophic aquaculture (IMTA) applications (Bischoff et al., 2009;Santos et al., 2016;Marques et al., 2018;Wang et al., 2019;Yousefi-Garakouei et al., 2019). A large variability in the total pool of FA determined for H. diversicolor biomass has been reported to date, with values ranging between 50 and 280 µg mg −1 DW for specimens fed commercial aquafeeds (Fidalgo e Costa et al., 2000;García-Alonso et al., 2008;Santos et al., 2016;Pajand et al., 2017;Wang et al., 2019) or between 24 and 110 µg mg −1 DW when cultured using aquaculture effluents (García-Alonso et al., 2008;Bischoff et al., 2009;Pajand et al., 2017;Marques et al., 2018;Wang et al., 2019;Yousefi-Garakouei et al., 2019). Also, the proportion of n-3 and n-6 FA reported to date for this species is also highly variable, with the above-mentioned works reporting values ranging between 5 and 33% for n-3 FA and 9-27% to n-6 FA. This variability is a consequence of several factors, such as the duration of experimental trials, the maturation stage of polychaetes at the beginning and during experiments, the composition of the supplied diet and the abiotic conditions experienced during culture (e.g., temperature, salinity, photoperiod). The main source of EFA in IMTA designs will always be the aquafeed supplied to upstream species. However, the majority of studies performed to date including polychaetes as extractive species under IMTA frameworks fails to characterize the proportion of feces and uneaten aquafeed present in the particulate fraction of nutrients being supplied. Additionally, only a few studies have evaluated the initial FA profile of polychaetes in order to assess how these biomolecules evolve when ragworms are supplied a certain pool of nutrients (Bischoff et al., 2009;Marques et al., 2018;Wang et al., 2019). It also remains to be clarified how the enrichment in EFA occurs in ragworms biomass over time under different combinations of water temperature and salinity, and whether it is possible to reach a plateau for enrichment in total FA, n-6 and n-3 FA. All this information is paramount to optimize the timeframe required to produce a premium polychaete biomass rich in EFA.
In order to shed some light over the above-mentioned questions, the present study aimed to evaluate the evolution of the FA profile of H. diversicolor fed a commercial aquafeed (with a well-known FA composition) during 10, 20, and 40 days (D10, D20 and D40) under different combinations of water temperature (20 and 25−C) and salinity (15, 20, and 25). These are optimal conditions for the culture of ragworms, which are within the range of temperature and salinity commonly employed by warmtemperate aquaculture systems operating with brackish water. A comparison between cultured and wild polychaetes over time was also performed to confirm that the evolution of FA profile of polychaetes fed on aquafeeds was not influenced by natural cycles. The feeding and growth performances of H. diversicolor cultured under the different combinations of water temperature and salinity were evaluated over time as well.

Experimental Set-Up
To evaluate the evolution of FA profile of H. diversicolor, five independent replicates were considered for each of the sixtemperature x salinity combinations, at each of the three-time sampling points considered (D10, D20, D40). Each replicate consisted of 10 polychaetes stocked on a glass flask (100 mm × 100 mm × 180 mm, ∼1-L volume) filled with 80 mm of sand (0.5-0.7 mm grain size) (approximately 1,000 ind. m −2 ) and aerated artificial sea water (see below for details). Five water baths were used for each of the temperatures being tested, with 100-W Eheim R thermocontrol 3612 aquarium heaters being used to control this parameter. Replicates of the different salinities were randomly distributed over these water baths. The different salinities were previously prepared using artificial sea water (prepared by mixing Red Sea R salt with tap water purified by a reverse osmosis unit), with the whole water volume of each glass flask being changed every 2 days. Constant water aeration was secured by air stones connected to an aeration pump (Hailea R vortex blower) and a 16:8 h light:dark photoperiod was employed using fluorescent white tubes connected to automatic timers. The schematic representation of the experimental set-up is represented in Figure 1.
Polychaetes were fed, once a day, a commercial grow-out diet for flatfish that present 62% of crude protein, 18% of crude fat, and 0.3% of crude fiber (WINFlat-SPAROS). Aquafeed was supplied ad libitum until polychaetes lost interest in the pellets and these remained untouched at the top of the sand in the glass flasks. Uneaten feed was never removed from the glass flasks, with its presence or absence being used to adjust the feed being supplied the next day.
Water parameters, namely dissolved oxygen (DO), pH, temperature and salinity were monitored daily using a manual probe (pH/Cond 3320, WTW, Weilheim, Germany), with ammonia and nitrites being determined once a week using colorimetric tests (Salifert Profi Test).

Polychaetes Stocking and Sampling
Wild specimens of H. diversicolor were collected from Ria de Aveiro coastal lagoon (Portugal) at Espinheiro Channel (40 • 38 8.4 N, 8 • 39 41.8 W) at the end of February 2019 and transported to the laboratory. An acclimation period of 24 h to experimental conditions was carried out and no feed was supplied during this period, so as not to modify the FA profile of wild polychaetes. The polychaetes were randomly selected and weighed individually and distributed per each of the five replicates of the different experimental conditions being tested. At the end of each sampling point (D10, D20, D40), polychaetes were depurated for 24 h in aerated glass flaks holding precombusted sterilized sand and artificial seawater (prepared as detailed above) to safeguard empty guts and avoid potential bias when performing FA analysis. The same temperature x salinity combination was maintained for each condition during the depuration process. Following depuration, all samples were freeze-dried and stored at −80 • C until further analysis.
Five composite samples of polychaetes (10 specimens per composite sample) were collected from the wild at the beginning of the experiment (February 2019) to characterize the FA profile of the initially stocked biomass (Initial). To clarify if the evolution of the FA profile of polychaetes being supplied with aquafeeds in the laboratory was not influenced by natural cycles, another five composite samples of wild specimens (10 specimens per composite sample) were collected in the same site at each sampling time point of the experiment (Wild10, Wild20, and Wild40) (February-March 2019). Depuration and storage processes for these samples were the same as detailed above.
Five composite samples of the aquafeed being used were also freeze dried and stored at −80 • C for FA analysis.

Lipid Extraction and Fatty Acid Analysis
The FA content of both polychaetes and aquafeeds was determined by screening FA methyl esters (FAME) analyzed by FIGURE 1 | Schematic representation of one of the five replicate water baths used to control temperatures (20 and 25 • C) housing randomly distributed glass flasks stocked with 10 Hediste diversicolor at different salinities (S-15, 20, and 25) to evaluate the evolution of their fatty acid profiles after being fed a commercial aquafeed during 10, 20, and 40 days (D10, D20, and D40).
gas chromatography-mass spectrometry (GC-MS). Before analysis, all freeze-dried samples were powdered and homogenized. Lipid extraction was performed following the Bligh and Dyer (1959) method using an initial sample biomass of 20 mg. Each sample was transferred to a glass tube, homogenized with 2,500 µL methanol (MeOH) and 1,250 µL dichloromethane (CH 2 Cl 2 ), sonicated for 1 min and then incubated on ice in an orbital shaker for 30 min. Afterward, each sample was centrifuged at 2,000 rpm for 10 min and the organic phase was collected. A volume of 1,250 µL CH 2 Cl 2 and 1,250 µL of Mili Q water was added to the total organic phase to promote phase separation, followed by centrifugation for 10 min at 2,000 rpm. Organic phase aliquots of lipid extract (74 µL for polychaetes and 57 µL for aquafeeds; total lipid extract: 3,000 µL) were collected for a new tube, previously washed in n-hexane. In order to define the final aliquot volume, a previous study was performed to calculate the total amount of lipid extract present in the different aliquots and analyzing them in GC-MS. Lipid extract aliquots were dried under a nitrogen gas stream, and posteriorly used for methylation. FAME were prepared by adding 1 mL of the internal standard 19:0 FA (1.5 µg mL −1 ) in n-hexane (99%) to the tube containing lipids. Subsequently, 200 µL of methalonic (MeOH) KOH solution (2M) were added and vigorously vortexed for 2 min. Following this procedure, 2 mL of NaCl saturated solution was added to the tube, and then centrifuged for 5 min at 2,000 rpm. Then 600 µL of the organic phase contained FAME were collected. To remove cholesterol from FAME, the solution was cleaned by solid-phase extraction. A column containing 0.1 g of silica was activated with 3 mL of n-hexane (99%) and the aliquot containing FAME (600 µL) was applied in the column. Subsequently, 3 mL of n-hexane 99%:diethyl ether (95:5 by volume) were eluted to recover the FAME. This fraction was dried under a nitrogen stream and stored at −20 • C until GC-MS analysis. FAME were then dissolved in 100 µL of n-hexane 99% and 2 µL of this solution were used for analysis on a GC-MS system (Agilent Technologies 5977 B GC/MSD) equipped with a DB-FFAP column (123-3232, J and W Scientific, Folsom, CA, United States) with the following specifications: 30 m length, an internal diameter of 0.32 mm and a film thickness of 0.25 µm. The system employed includes a Mass Selective Detector operating with an electron impact mode at 70 eV and scanning the mass range m/z 50-550 in a 1 s cycle in a full scan mode acquisition. The oven temperature program was as follows: (1) initial temperature of 80 • C for 2 min; (2) linear increase to 160 • C (25 • C min −1 ); (3) linear increase to 210 • C (2 • C min −1 ); (4) linear increase to 225 • C (20 • C min −1 ); and (5) standing at 225 • C for 20 min. Helium was used as the carrier gas (1 mL min −1 ). FA were integrated through Agilent's Masshunter Solutions Quantitative Analysis v10.0 automatic integration. Identification of FA was performed considering retention time and analysis of MS spectra in comparison with MS spectra of FA standards (Supelco 37 Component FAME Mix, ref. 47885-U, Sigma-Aldrich) and comparison with chemical databases (Wiley 275 library, AOCS lipid library, and NIST 2014 Mass spectral library). The FA content (expressed as µg of FA mg −1 dry weight, DW) was determined using calibration curves obtained from FAME certified standard mixture (Supelco R 37 Component FAME Mix, Sigma-Aldrich) and 19:0 FA as internal standard.
In the present study, PUFA were defined as all FA with two or three double bonds, while HUFA refers to all FA with four or more double bonds.
FA increments between day 1-10, day 11-20, and day 21-40 (D1-D10, D11-D20, and D21-D40) were determined by calculating the differences in concentration between samples collected at the beginning of the experiment and D10, between samples from D10 and D20, and between samples from D20 and D40, respectively.

Statistical Analysis
A two-way nested analysis of variance (ANOVA) was performed to determine the existence of significant differences in SGR, DGR and FR of H. diversicolor cultured under different combinations of water temperature and salinity. This analysis was performed for each independent group fed a commercial aquafeed during 10, 20, and 40 days (D10, D20, and D40) (n = 5). The factors being tested were as follows: "Temperature" (two fixed levels, 20 and 25 • C) and "Salinity" (three fixed levels, 15, 20, and 25) nested within temperature. Data were previously checked for normality (Anderson-Darling test) and homogeneity of variances (Bartlett's and Levene's tests for normal and non-normal distribution, respectively). When normality was not verified, the hypotheses were tested using non-parametric Kruskal-Wallis test. Significant differences were always considered at p < 0.05. These statistical analyses were performed using MINITAB 18 Software (State College, PA).
To determine the existence of significant differences in the FA profile of cultured H. diversicolor under the different combinations of water temperature and salinity tested in the present work, a two-way nested analysis of similarities (ANOSIM) was performed on a resemblance matrix produced using Bray Curtis similarity coefficient of data previously transformed using the formula log (x + 1). This analysis was performed for each independent group at D10, D20, and D40 days (n = 5). The factors being tested were as follows: "Temperature" (two levels, 20 • C and 25 • C) and "Salinity" (three levels, 15, 20, and 25) nested within temperature. The differences between FA profile of cultured polychaetes under different 16.06 ± 1.68 18:1 n-7 5.69 ± 0.55 20:1 n-13 1.10 ± 0.04 20:1 n-11 2.73 ± 0.25 20:1 n-9 1.20 ± 0.03 22:1 n-11 2.94 ± 0.24 MUFA 35.56 ± 3.46 18:2 n-6 (LA) 14.32 ± 1.49 18:3 n-3 (ALA) 3.13 ± 0.28 combinations of water temperature and salinity and the FA profile of wild conspecifics collected in the same sampling period (Wild10, Wild20, and Wild40) (n = 5) were performed using a one-way ANOSIM. When significant differences were recorded, a similarity percentages (SIMPER) analysis was performed to evaluate which FA contributed the most to the dissimilarities being recorded between samples, until a total of 50% cumulative dissimilarity was achieved. A principal coordinates analysis (PCO) plot was also used to compare the results of total FA concentration, FA classes concentration (saturated FA [SFA], monounsaturated FA [MUFA], PUFA, and HUFA) and n-3 and n-6 FA concentration from D40 polychaetes cultured under different combinations of temperature and salinity, Wild40 polychaetes, initially stocked polychaetes (Initial) and the aquafeed supplied (resemblance matrix produced using Bray Curtis similarity used as input). These statistical analyses were performed using PRIMER v6 with the PERMANOVA + add-on. For a detailed description of all the statistical methods employed please see Anderson et al. (2008).

Experimental Conditions
The average values of DO and pH inside the glass flask varied between 6.7-7.9 mg L −1 and 8.2-8.4 along the duration of the experiment, respectively. The average value of ammonia (NH 4 ) and nitrites (NO 2 ) varied between 0.4-1.2 and 0.5-1.4 mg L −1 , respectively. Temperature and salinity were kept stable in the predefined conditions (see Supplementary Table 1).

Feeding and Growth Performance of Polychaetes
The FR (%), SGR (% d −1 ), and DGR (g d −1 ) determined for polychaetes cultured under each combination of water temperature and salinity are summarized in
A progressive increment of total FA content over time was recorded in all combinations of water temperature and salinity evaluated in the present work, with values varying from 30 µg mg −1 DW (Initial) to 70-90 µg mg −1 DW (D40) (Figure 2). The proportion of PUFA incremented approximately 10% in the whole pool of FA over the 40 days of the experimental trial, with all FA classes (SFA, MUFA, PUFA, and HUFA) being present with a similar proportion at the end of the study (˜25% in all FA classes) (Supplementary Figure 1).

Fatty acid
PCO plots revealed two well-separated groups, one formed by D40 polychaetes and aquafeed, and the other one formed by Initial and Wild40 polychaetes (Figures 5A-C). A similarity above 90% was recorded between D40 polychaetes and aquafeed for total FA concentration, FA classes concentration (SFA, MUFA, PUFA, and HUFA) and n-3 and n-6 FA concentration. The two PCO axis explained more than 99% of the variation recorded between samples from different group.
The amount of commercial aquafeed supplied to polychaetes stocked under the different treatments performed allowed to estimate an increment of 25.0-37.0 µg mg −1 DW in the total FA present in polychaetes biomass. Concerning n-3 and n-6 FA, the increments recorded were 5.6-8.0 and 4.8-7.7 µg mg −1 DW in polychaetes biomass, respectively (Table 4). DHA and EPA increased between 2.0-2.6 and 2.3-3.2 µg mg −1 DW in polychaetes biomass, respectively.

DISCUSSION
In the present study it was possible to verify that H. diversicolor fed a commercial aquafeed during 10, 20, and 40 days (D10, D20, and D40) under different combinations of water temperature (20 and 25 • C) and salinity (15, 20, and 25) displayed similar SGR and DGR. This finding is in line with the ecophysiological traits of this species, which is very abundant in extreme environments, such as intertidal mudflats, exposed to major shifts in water temperature and salinity (Scaps, 2002;García-Alonso et al., 2008). The values of SGR reported in the present work were lower than those reported by previous studies testing commercial aquafeeds to culture this species (values ranging between 6.0 and 6.8% d −1 ) (Fidalgo e Costa et al., 2000;Santos et al., 2016). However, our values were higher than most of those reported when this species was fed on nutrients derived from  aquaculture effluents (values ranging between 1.2 and 3.0% d −1 ) (Bischoff, 2007;Pajand et al., 2017;Wang et al., 2019). The values of DGR determined for D10 polychaetes were similar to those reported by Santos et al. (2016) for the same species fed a commercial aquafeed for 60 days. The progressive decrease of DGR over time and the lower SGR recorded for D40 polychaetes, may be related with their sexual maturation, as detailed below. It is worth highlighting that D40 polychaetes incremented their average weight up to 4-times, confirming the potential of this species to be produced under an IMTA framework and to generate a considerable biomass in a short period of time. The progressive mortality verified in D10, D20, and D40 samples, which reached values ∼50%, may have been related to the high stocking density employed at the beginning of the experiment (1,000 ind. m −2 ). This assumption is supported by the results reported by Nesto et al. (2012), who verified that H. diversicolor exhibited a higher growth and survival when cultured at a density of 300 ind. m −2 than at 1,000 and 3,000 ind. m −2 . Total FA reported for D40 polychaetes in the present work are among the highest values recorded to date for H. diversicolor (see Table 5). Geographic location, abiotic conditions, diet FIGURE 7 | Schematic representation of putative biosynthetic pathways of some polyunsaturated fatty acids (PUFA) identified in Hediste diversicolor from the wild and after being fed on a commercial aquafeed for 40 days. Methyl-end desaturases (ω, green arrows); Front-end desaturases (red arrows); Elongation reactions (Elo, blue arrows). The " " refers to the carbon position at which the incipient double bound (unsaturation) locates within the methyl and front ends of fatty acyl chains. AdA, adrenic acid; ALA, α-linolenic acid; ARA, arachidonic acid; DGLA, dihomo-gamma-linolenic acid; DHA, docosahexaenoic acid; DPA, Docosapentaenoic acid; EPA, eicosapentaenoic acid; ETA, eicosatetraenoic acid; ETE, eicosatrienoic acid; LA, linoleic acid. Scheme adapted from Monroig and Kabeya (2018). availability and composition, and maturation state contribute to the shaping of FA profiles in this species (Luis and Passos, 1995;García-Alonso et al., 2008;Santos et al., 2016;Marques et al., 2018). Therefore, the interpretation of contrasting FA profiles reported in different studies is a challenging task. Initially stocked polychaetes, as well as wild conspecifics Wild10, Wild20, and Wild40, were collected during February-March, a transition period between winter and spring in the sampling area. Wild specimens exhibited a reddish-brown coloration (Figure 6A), a characteristic feature of immature worms (Scaps, 2002). Two spawning periods, one in spring and the other in early autumn, were identified in wild populations of H. diversicolor monitored at the same sampling site of the specimens used in the present study (Abrantes et al., 1999). According to Luis and Passos (1995), during winter and early spring months, in the presence of a suitable diet, sexually immature worms accumulate lipids in their body. Immediately before reproduction, a bulk of FA is made available through triacylglycerol depletion, being these essential during the rapid growth of oocytes (Luis and Passos, 1995). This polychaete species is therefore well-known to reflect in its body composition the FA profile of its diet (Fidalgo e Costa et al., 2000;Santos et al., 2016;Pajand et al., 2017;Marques et al., 2018;Wang et al., 2019). This feature, along with maturation, are likely among the factors which most contributed for the differences in FA composition reported in previous studies testing the production of H. diversicolor using commercial aquafeeds or wasted nutrients derived from aquaculture effluents. In the present work it was possible to monitor the effect of diet over time (by comparing Initial, D10, D20, and D40 samples). The FA which exhibited higher increments in polychaete biomass were exactly the ones presented in higher levels in the commercial aquafeed employed. With only a few exceptions, these FA are also among the most wellrepresented ones reported in previous studies addressing the culture of this species (Fidalgo e Costa et al., 2000;García-Alonso et al., 2008;Bischoff et al., 2009;Santos et al., 2016;Pajand et al., 2017;Marques et al., 2018;Wang et al., 2019;Yousefi-Garakouei et al., 2019).
An advantage of the present study in relation to previous works is that it monitored the evolution of FA profile of ragworms over time. This issue is of paramount importance if producers aim to enhance the amount of EFA of cultured ragworms by using enriched diets for a short period of time at the end of the production cycle to obtain premium polychaetes. In addition, producers may also select to enrich wild polychaetes biomass, rather than culturing these organisms over their whole life cycle. A progressive increment of total FA was recorded over time. In the first 20 days (comparing Initial, D10 and D20 polychaetes), a predominant increment of SFA and MUFA was evidenced, while in the last 20 days of the trial (comparing D20 and D40 polychaetes) an increment of PUFA and HUFA prevailed. These results may well reflect an adaptation of polychaetes to the new diet being provided and/or the experimental conditions employed, latter followed by metabolic shifts triggered by maturation. According to Nesto et al. (2012), diets similar to the one used in the present work, with high-protein content, despite originating the best results of daily biomass production, also triggered an earlier onset of gametogenesis. Despite the increments recorded, it is worth highlighting that a plateau of total FA, n-6 FA and n-3 FA was not reached during the 40 days of the experimental trials. This finding allows us to infer that it may still be possible to further enhance H. diversicolor biomass in EFA if the study had been extended over a longer period. Nonetheless, one must acknowledge that when this species reaches adult size (similar to that recorded at D40 polychaetes) it starts to reproduce and die.

CONCLUSION
The present study demonstrated how H. diversicolor biomass can be successfully enriched with n-6 and n-3 FA (including EPA and DHA) when provided a commercial aquafeed, even when exposed to different combinations of water temperature and salinity. This finding highlights the potential of using H. diversicolor as an extractive species in IMTA designs, thus allowing to recover valuable nutrients present in aquaculture effluents that would otherwise be wasted. The production of a ragworm biomass rich in EFA allows to provide the aquafeed industry with another alternative ingredient for fish meal and fish oil, at least for the formulation of premium maturation and finishing diets. The development of innovative aquaculture production models based on the integration of extractive species, such as polychaetes, to recover wasted nutrients derived from the aquaculture of fed species complies with sustainable aquaculture guidelines and is aligned with UN 2030 Sustainable Development Goals Agenda.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

AUTHOR CONTRIBUTIONS
DJ and RC contributed to conceptualization and design of the study. DJ and HI performed the experiment culturing polychaetes under different conditions and collected the biomass samples. DJ, FR, and MD designed and performed the analytical work to evaluate fatty acid profile. DJ performed the statistical analysis and wrote the first draft of the manuscript. RC, AL, and MD contributed with supervision, funding and resources acquisition. All authors contributed to manuscript revision, read, and approved the submitted version.

FUNDING
This work was supported by the project "VALORMAR-Full valorization of marine resources: potential, technological innovation and new applications" (24517, 10/SI/2016-Mobilizing R&TD Programs) co-funded by Compete 2020, Lisboa 2020, Algarve 2020, Portugal 2020 and European Union, through the European Regional Development Fund (ERDF), by project "AquaMMIn-Development and validation of a modular integrated multitrophic aquaculture system for marine and brackish water species" (MAR-02.01.01-FEAMP-0038) co-funded by Portugal 2020 and the European Union through Mar2020, the Operational Programme (OP) for the European Maritime and Fisheries Fund (EMFF) in Portugal and by the Integrated Programme of SR&TD "SmartBioR-Smart Valorization of Endogenous Marine Biological Resources Under a Changing Climate" (Centro-01-0145-FEDER-000018), and co-funded by Centro 2020 program, Portugal 2020, European Union, through the European Regional Development Fund. Daniel Jerónimo acknowledges FCT (Portuguese Foundation for Science and Technology) by financially supporting his Ph.D. grant (PD/BD/127989/2016) and Felisa Rey by her Junior Researcher contract (reference CEECIND/00580/2017). Thanks for financial support are also due to CESAM (UIDB/50017/2020 + UIDP/50017/2020), to FCT/MEC through national funds, and the cofounding by the FEDER, within the PT2020 Partnership Agreement and Compete 2020.