Influence of Dietary Long-Chain Polyunsaturated Fatty Acids and ω6 to ω3 Ratios on Head Kidney Lipid Composition and Expression of Fatty Acid and Eicosanoid Metabolism Genes in Atlantic Salmon (Salmo salar)

The interaction of dietary eicosapentaenoic acid and docosahexaenoic acid (EPA+DHA) levels with omega-6 to omega-3 ratios (ω6:ω3), and their impact on head kidney lipid metabolism in farmed fish, are not fully elucidated. We investigated the influence of five plant-based diets (12-week exposure) with varying EPA+DHA levels (0.3, 1.0, or 1.4%) and ω6:ω3 (high ω6, high ω3, or balanced) on tissue lipid composition, and transcript expression of genes involved in fatty acid and eicosanoid metabolism in Atlantic salmon head kidney. Tissue fatty acid composition was reflective of the diet with respect to C₁₈ PUFA and MUFA levels (% of total FA), and ω6:ω3 (0.5–1.5). Fish fed 0.3% EPA+DHA with high ω6 (0.3% EPA+DHA↑ω6) had the highest increase in proportions (1.7–2.3-fold) and in concentrations (1.4-1.8-fold) of arachidonic acid (ARA). EPA showed the greatest decrease in proportion and in concentration (by ~½) in the 0.3% EPA+DHA↑ω6 fed fish compared to the other treatments. However, no differences were observed in EPA proportions among salmon fed the high ω3 (0.3 and 1.0% EPA+DHA) and balanced (1.4% EPA+DHA) diets, and DHA proportions were similar among all treatments. Further, the transcript expression of elovl5a was lowest in the 0.3% EPA+DHA↑ω6 fed fish, and correlated positively with 20:3ω3, 20:4ω3 and EPA:ARA in the head kidney. This indicates that high dietary 18:3ω3 promoted the synthesis of ω3 LC-PUFA. Dietary EPA+DHA levels had a positive impact on elovl5a, fadsd5 and srebp1 expression, and these transcripts positively correlated with tissue ΣMUFA. This supported the hypothesis that LC-PUFA synthesis is positively influenced by tissue MUFA levels in Atlantic salmon. The expression of pparaa was higher in the 0.3% EPA+DHA↑ω6 compared to the 0.3% EPA+DHA↑ω3 fed fish. Finally, significant correlations between head kidney fatty acid composition and the expression of eicosanoid synthesis-related transcripts (i.e., 5loxa, 5loxb, cox1, cox2, ptges2, ptges3, and pgds) illustrated the constitutive relationships among fatty acids and eicosanoid metabolism in salmon.

Introduction

Long-chain polyunsaturated fatty acids (LC-PUFA) such as eicosapentaenoic (EPA, 20:5ω3), docosahexaenoic (DHA, 22:6ω3) and arachidonic (ARA, 20:4ω6) acids play important roles in fish growth, metabolism (Sargent et al., 2002), neural development, health and reproduction (Sargent et al., 1999; Tocher, 2010, 2015). They are also important structural components of cell membranes (Calder, 2013), and act as precursors to eicosanoid metabolites such as prostaglandins, leukotrienes, thromboxanes, docosanoids, maresins and resolvins which regulate inflammatory and immune response processes (Calder, 2007, 2013; Arts and Kohler, 2009; Martinez-Rubio et al., 2013; Holen et al., 2015; Montero et al., 2019). The levels of EPA and DHA in fish feeds showed a large reduction (from ~ 20 to 5% of dietary lipid) over the past ~3 decades (Einen and Thomassen, 1998; Ytrestøyl et al., 2015; Sprague et al., 2016), due to decreasing global availability and rising market price of their fish oil (FO) sources (Turchini et al., 2009; Tacon and Metian, 2015; Froehlich et al., 2018). The reduction in FO and fish meal (FM), and the concomitant rise in incorporation of plant-based ingredients have influenced fish health (e.g., liver steatosis, altered phagocytic activity, and modulation of immune-related transcript expression) (Montero et al., 2003; Ruyter et al., 2006; Caballero-Solares et al., 2017), and their nutritional quality for human consumers (Sprague et al., 2016; De Roos et al., 2017).

Salmonids have a capacity to desaturate and elongate the precursors 18:3ω3 and 18:2ω6 to LC-PUFA metabolites (Tocher et al., 2000; Hixson et al., 2014; Katan et al., 2019). The enzymes involved in the bioconversion of both ω3 and ω6 pathways are the fatty acyl desaturases (FADS2D5, FADS2D6a, FADS2D6b and FADS2D6c: Zheng et al., 2005; Monroig et al., 2010) and fatty acyl elongases (ELOVL5a, ELOVL5b, ELOVL2, ELOVL4a and ELOVL4b: Hastings et al., 2004; Morais et al., 2009; Carmona-Antoñanzas et al., 2011; Zhao et al., 2019). Hepatic LC-PUFA synthesis in Atlantic salmon (Salmo salar) is influenced by precursor (18:2ω6 and/or 18:3ω3) availability, as well as the concentration of LC-PUFA (i.e., EPA, ARA, DHA) (Jordal et al., 2005; Glencross et al., 2015; Katan et al., 2019). This pathway is regulated at the transcriptional level and is mediated by transcription factors [e.g., liver X receptor (LXR), sterol regulatory element binding protein (SREBP) 1 and 2, and peroxisome proliferator-activated receptor (PPAR) α, β and γ (Carmona-Antoñanzas et al., 2014; Hixson et al., 2017; Katan et al., 2019; Emam et al., 2020)]. However, the constitutive regulation of this pathway in salmon head kidney is less well-understood (Betancor et al., 2014), and requires further investigation given the central role of this organ in haematopoetic and immune processes (Tort et al., 2003; Zapata et al., 2006; Gjøen et al., 2007).

Plant oils used in aquafeeds generally contain high levels of ω6 FA, and as such, these oils and farmed fish consuming them may provide inadequate ratios of ω6 to ω3 (ω6:ω3) fatty acids for the consumer (Pickova and Mørkøre, 2007; Weaver et al., 2008; Young, 2009). Further, ω6:ω3 plays an important role in fish immune response [through their conversion into pro- and anti-inflammatory eicosanoids (Furne et al., 2013; Holen et al., 2018)] and FA metabolism [e.g., impacting LC-PUFA synthesis and transcript expression of genes involved in this pathway (Torstensen and Tocher, 2010; Vagner and Santigosa, 2011; Katan et al., 2019)]. Previous studies revealed that diets with varying ω6:ω3 ratio (i.e., ~0.3–4.0) influenced the transcription of inflammation- and FA metabolism-related genes in Atlantic salmon head kidney (Martinez-Rubio et al., 2013; Holen et al., 2018). However, the interaction of dietary ω6:ω3 with EPA+DHA levels, and their nutritional impacts on eicosanoid and fatty acid metabolism in salmon, have not been fully elucidated.

Emam et al. (2020) recently showed that diets with different combinations of ω6:ω3 (high ω6, high ω3, or balanced) and EPA+DHA levels (0.3, 1.0, or 1.4%) affected muscle lipid composition, and the hepatic transcript expression of fatty acid synthesis-related genes in Atlantic salmon. The current study used the same diets and fish population as in Emam et al. (2020), with the aim to investigate how diets with low FO inclusion levels impact head kidney lipid composition, and the transcript expression of genes involved in fatty acid and eicosanoid metabolism. Furthermore, correlation analyses were used to relate head kidney lipid composition with transcript expression.

Materials and Methods

Experimental Diets and Animals

Five experimental diets were manufactured by Cargill Canada (formerly EWOS Canada; Surrey, BC, Canada), top-coated with different oil mixes at the Chute Animal Nutrition Centre (Dalhousie University, Truro, NS, Canada), and stored at −20°C until needed. Experimental diets contained different levels of FO, poultry fat, and vegetable oils (VO; i.e., soy oil, linseed oil, and rapeseed oil) to generate three levels of EPA+DHA (0.3, 1.0, and 1.4%, as formulated) and contrasting ω6:ω3 ratios (high ω6, high ω3, and balanced). The diets were as follows: 0.3% EPA+DHA with higher ω6 (0.3% EPA+DHA↑ω6), 0.3% EPA+DHA with higher ω3 (0.3% EPA+DHA↑ω3), 1.0% EPA+DHA with higher ω6 (1%EPA+DHA↑ω6), 1.0%EPA+DHA with higher ω3 (1%EPA+DHA↑ω3), and 1.4%EPA+DHA with a more balanced ω6:ω3 ratio (1.4% EPA+DHA/balanced). All diets were formulated to be isonitrogenous and isoenergetic and to meet the nutritional requirements of salmonids (National Research Council (NRC), 2011). Dietary formulations and their lipid compositions were published previously (Emam et al., 2020). However, as they pertain to the current study, they are included as Supplementary Material here (Supplementary Tables 1, 2).

Atlantic salmon smolts were transported from a regional salmon farm to the Dr. Joe Brown Aquatic Research Building (Ocean Sciences Centre, Memorial University of Newfoundland, Canada; October 2016). After their arrival, fish were graded by size in order to select the most uniform population, PIT (Passive Integrated Transponder; Easy AV, Avid Identification Systems, Norco, CA, USA)-tagged for individual identification, and kept in 3,800-litre tanks until the beginning of the feeding trial. Then, salmon post-smolts [210 ± 44 g mean initial weight ± standard deviation (SD)] were randomly distributed into twenty 620-litre tanks (40-41 fish tank⁻¹), and subjected to an 8-week acclimation period during which they were given a commercial diet (EWOS Dynamic S, Cargill Inc., MN, USA). After the completion of the acclimation period, fish were switched from the commercial feed and fed with the experimental diets (4 tanks diet⁻¹) for 12 weeks. At all stages, fish were fed to satiation overnight using automatic feeders (AVF6 Vibrating Feeders, Pentair Aquatic Eco-Systems, Inc., Florida, USA). Uneaten pellets were collected every morning, and the daily amount of feed was adjusted based on the number of uneaten pellets. All tanks were supplied with 12°C flow-through filtered seawater at 12 l min⁻¹; the dissolved oxygen level was ~10 mg l⁻¹, and the photoperiod was maintained at 24-h light. Mortalities were weighed and recorded throughout the trial. Ethical treatment of fish in this experiment was carried out in accordance with the guidelines of the Canadian Council on Animal Care, and approved by the Institutional Animal Care Committee of Memorial University of Newfoundland (protocol # 16-75-MR).

Sample Collection

Salmon samples were collected from 10 opportunistically netted fish (each from a different tank) at week 0 (the day before fish were exposed to the experimental diets), and from 5 fish per tank at week 12. Sampling occurred after the fish were starved for 24 h and euthanized with an overdose of MS-222 (400 mg l⁻¹; Syndel Laboratories, Vancouver, BC, Canada). For lipid analyses, head kidney samples were collected in 15 ml glass test tubes that had been rinsed three times with methanol, followed by three rinses with chloroform. Samples were stored on ice during sampling, their wet weights recorded, and they were covered with 2 ml of chloroform (HPLC-grade). Finally, the glass test tubes were filled with nitrogen, the Teflon-lined caps were sealed with Teflon tape, and samples were stored at −20°C. For gene expression analyses, head kidney (50–100 mg sample⁻¹) tissues were collected in 1.5 ml nuclease-free tubes, flash-frozen in liquid nitrogen, and stored at −80°C until RNA extractions were performed.

Lipid Extractions

Lipid samples were extracted according to Parrish (1999). Samples were homogenized in a 2:1 mixture of ice-cold chloroform:methanol using a Polytron PCU-2-110 homogenizer (Brinkmann Instruments, Rexdale, ON, Canada). Chloroform-extracted water was added to bring the ratio of chloroform:methanol:water to 8:4:3. Each sample was sonicated for 4 min in an ice bath and centrifuged at 3,000 rpm for 2 min at room temperature. The bottom, organic layer was removed using a double pipetting technique, placing a 2 ml lipid-clean glass Pasteur pipette inside a 1 ml glass pipette, in order to remove the organic layer without disturbing the top, aqueous layer. Chloroform was then added back to the extraction test tube, and the entire procedure was repeated three times. All organic layers were pooled into a separate lipid-clean vial. Finally, samples were concentrated under a flow of nitrogen gas.

Lipid Class Separation

Iatroscan Mark 6 TLC-FID (thin-layer chromatography-flame ionization detector) (Mitsubishi Kagaku Iatron, Inc., Tokyo, Japan), and a three-step development method were used to determine the lipid class composition (Parrish, 1987). Lipid extracts and standards were applied to the Chromarods and focused to a narrow band using 100% acetone. The first development system was hexane:diethyl ether:formic acid (98.95:1.0:0.05). The rods were developed for 25 min, removed from the system for 5 min to dry, and replaced for 20 min. The second development was for 40 min in hexane:diethyl ether:formic acid (79:20:1). The final development system had two steps–the first was in 100% acetone for two 15 min time periods, followed by two 10 min periods in chloroform:methanol:chloroform extracted water (5:4:1). Before using each solvent system, the rods were dried in a constant humidity chamber. After each development system, the rods were partially scanned in the Iatroscan, and the data were collected using Peak Simple software (SRI Instruments, version 3.67, Torrance, CA, USA). The Chromarods were calibrated using standards from Sigma Chemicals (Sigma Chemicals, St. Louis, MO, USA).

Fatty Acid Methyl Ester (FAME) Derivatization

Lipid extracts (50 μl) were transferred into lipid-clean 15 ml glass vials and concentrated under a flow of nitrogen to complete dryness. Then, 1.5 ml of methylene chloride and 3 ml of Hilditch reagent (1.5 sulfuric acid: 98.5 anhydrous methanol) were added to each sample vial, and this was followed by a brief vortexing and 4 min sonication (Fisher Scientific FS30, Pittsburgh, PA, USA). Vials were then filled with nitrogen, capped and heated at 100°C for 1 h. Subsequently, 0.5 ml of saturated sodium bicarbonate solution and 1.5 ml of hexane were added to each vial, followed by a brief vortexing, and then the upper, organic layer was removed into a separate lipid-clean glass vial. Each sample was then dried and refilled with ~0.5 ml of hexane. Finally, vials were filled with nitrogen, capped, sealed with Teflon tape, and sonicated for 4 min to resuspend the fatty acids. All solvents were OmniSolv grade (VWR International, Mississauga, Ontario, Canada).

All FAMEs were analysed on a HP 6890 GC-FID (gas chromatograph-flame ionization detector) system equipped with a 7683 autosampler (Agilent Technologies Canada Inc., Mississauga, Ontario, Canada). Additional details regarding the GC conditions (e.g., column length, diameter, temperature, carrier gas) are described in Hixson et al. (2017). Fatty acid peaks were identified against known standards (PUFA 1, PUFA 3, BAME and a Supelco 37 component FAME mixture: Sigma-Aldrich Canada Ltd., Oakville, Ontario, Canada). Finally, chromatograms were integrated using the Varian Galaxie Chromatography Data System, version 1.9.3.2 (Walnut Creek, CA, USA). FA data were expressed as area percentage of FAME.

RNA Extraction, DNase Treatment, Column Purification and cDNA Synthesis

Head kidney samples were homogenized in TRIzol reagent (Invitrogen, Carlsbad, CA, USA), using the TissueLyser II system at 25 Hz for 2.5 min, using 5 mm stainless steel beads (QIAGEN, Mississauga, ON, Canada), and subjected to RNA extraction. This was followed by DNaseI treatment and column purification using the RNase-free DNase Set (QIAGEN) and the RNeasy Mini Clean-up Kit (QIAGEN), respectively. All procedures were conducted according to manufacturer instructions and as described in Xue et al. (2015). RNA integrity was verified by 1% agarose gel electrophoresis, and RNA purity and quantity were assessed by Nanodrop UV spectrophotometry (Thermo Fisher Scientific, Waltham, MA, USA). DNaseI-treated and column-purified RNA samples had A260/280 and A260/230 ratios of 1.8–2.3 (data not shown). All cDNAs were prepared by reverse transcription of 1 μg of DNaseI-treated, column-purified total RNA for each sample, using 1 μl of random primers (250 ng; Invitrogen), 1 μl of dNTPs (0.5 mM final concentration; Invitrogen), 4 μl of 5X first-strand buffer (1×final concentration; Invitrogen), 2 μl of DTT (10 mM final concentration; Invitrogen), and 1 μl of Moloney murine leukemia virus (M-MLV) reverse transcriptase (RT) (200 U; Invitrogen) at 37°C for 50 min. The total reaction volume was 20 μl, following manufacturer instructions, and as described in Xue et al. (2015). Finally, all cDNAs were diluted 40 times with nuclease-free water (Invitrogen) prior to real-time quantitative polymerase chain reaction (qPCR).

QPCR Analysis

Expression of 25 genes of interest (GOI) (Table 1) related to FA and eicosanoid metabolism was quantified by qPCR using head kidney cDNA templates from fish fed the five dietary treatments at week 12. For each dietary treatment, eight individuals (two from each tank) were used in the qPCR study. Only head kidney samples from fish that had specific growth rates within 1.5 standard deviations below and above the mean value of each tank were selected for the qPCR study, in order to reduce biological variability. The sequences of all primer pairs used in qPCR analyses, GenBank accession numbers of sequences used for primer design, and other details are presented in Table 1. Notably, primers for the transcripts fadsd6a, elovl2, and elovl5b (GenBank accession numbers AY458652, FJ237532 and FJ237532, respectively) failed quality testing due to low transcript expression levels, and thus were not included in the qPCR study. Each primer pair was quality-tested to verify that a single product was amplified with no primer dimers, and included standard curves and dissociation curves, as described in Rise et al. (2010) and Booman et al. (2011). In brief, the amplification efficiency (Pfaffl, 2001) of each primer pair was determined using a 5-point 1:3 dilution series starting with cDNA representing 10 ng of input total RNA. Two pools were generated (i.e., 0.3% EPA+DHA↑ω6 and 0.3% EPA+DHA↑ω3), and each pool consisted of 8 fish (with each fish contributing an equal quantity to the pool). The reported primer pair amplification efficiencies were an average of the two pools, except if one pool showed poor efficiency or spacing due to low expression levels (i.e., the 0.3% EPA+DHA↑ω3 pool was used for fadsd6b and fadsd6c, while the 0.3% EPA+DHA↑ω6 pool was used for lxrb).

TABLE 1

Table 1. qPCR primers.

To select the most suitable normalizer genes, six candidate normalizers were tested based on our previous qPCR studies (rpl32, eef1α-1, eef1α-2, actb, eif3d, pabpc1) (Xue et al., 2015; Caballero-Solares et al., 2017), and salmon literature on reference genes (actb, eef1α-1, eef1α-2) (Olsvik et al., 2005). Their qPCR primers were quality-tested as mentioned above. Half of the fish population involved in the qPCR study was utilized for normalizer testing (i.e., 20 total which consisted of 4 fish per treatment). Cycle threshold (C_T) values were measured using cDNA corresponding to 5 ng of input total RNA. Expression stability was then analysed using the geNorm algorithm (Vandesompele et al., 2002). Eif3d and pabpc1 were shown to be the most stable (i.e., geNorm M-values of 0.23 and 0.22, respectively) among the six candidate genes, and therefore were selected as normalizers.

All qPCR amplifications were performed in a total reaction volume of 13 μl and consisted of 4 μl of cDNA (5 ng input total RNA), 50 nM each of forward and reverse primer, 1×Power SYBR Green PCR Master Mix (Applied Biosystems), and nuclease-free water (Invitrogen). QPCR reactions, including no-template controls, were performed in technical triplicates using the ViiA 7 Real-Time PCR System (384-well format) (Applied Biosystems, Foster City, CA, USA) and the Power SYBR Green I dye chemistry. The Real-Time analysis program consisted of 1 cycle of 50°C for 2 min, 1 cycle of 95°C for 10 min, followed by 40 cycles (of 95°C for 15 s and 60°C for 1 min), with the fluorescence signal data collection after each 60°C step. When a C_T value within a triplicate was >0.5 cycle from the other two values, it was considered to be an outlier, discarded and the average C_T of the remaining two values was calculated. The relative quantity (RQ) of each GOI was calculated using a qBase relative quantification framework (Hellemans et al., 2007; Booman et al., 2014) with primer amplification efficiencies incorporated (Table 1). The expression levels of each GOI were normalized to both normalizer genes, and the sample with the lowest normalized expression was used as the calibrator sample (i.e., RQ = 1.0) for each GOI (Rise et al., 2015). Transcript expression data are presented as RQ values relative to the calibrator. QPCR fold-change values were calculated by dividing the mean RQ value of each dietary treatment by that of the 0.3% EPA+DHA↑ω6 fish.

Statistical Analyses

Tissue Lipid Composition and qPCR Data

Nested general linear models were used with tank nested in diet (Minitab 17 Statistical Software, State College, PA, USA). This was followed by Tukey post-hoc tests (p < 0.05) to identify significant differences among treatments at week 12 (Minitab 17 Statistical Software). However, when a significant tank effect was identified (p < 0.05), a one-way ANOVA followed by Tukey post-hoc tests was performed (Minitab 17 Statistical Software, State College, PA, USA). To identify significant differences in head kidney lipid composition between week 0 (i.e., initial) and week 12, a one-way ANOVA followed by Tukey post-hoc tests was used. To show the effects of dietary EPA+DHA (i.e., 0.3 and 1%) and ω6:ω3 (i.e., high ω6 and high ω3) factors on the transcript expression of each GOI, a two-way ANOVA was performed (Minitab 17 Statistical Software). The 1.4% EPA+DHA/balanced treatment was not included in this analysis since it was not formulated following the 2 × 2 factorial design of the other four treatments (i.e., 0.3 or 1.0% EPA+DHA vs. high ω6 or high ω3).

For the qPCR analysis, each dietary treatment was tested for outliers using Grubb's test (p < 0.05). In total, 7 RQ values were identified as statistical outliers in the entire dataset (i.e., out of 1,000 values comprising all samples and all GOIs), and excluded from the study. All GOIs had a sample size of 7–8 per dietary treatment. Finally, residuals were tested to verify normality, independence and homogeneity of variance. Normality was examined using the Anderson-Darling test. If the test failed (p < 0.05), a one-way ANOVA on ranks was performed, and followed by a Kruskal-Wallis test (SigmaPlot, Systat Software, Inc., Version 13, San Jose, CA, USA). In all cases, differences were considered statistically significant when p < 0.05.

Multivariate and Pearson Correlation Analyses

Principal coordinate analysis (PCoA) was used to describe head kidney lipid and FA composition in the five dietary treatments (PRIMER, Plymouth Routines in Multivariate Ecological Research; PRIMER-E Ltd., version 6.1.15, Ivybridge, UK). A similarity matrix was used, and the first two PCO axes (i.e., PCO1, PCO2) were plotted. SIMPER (Similarity of Percentages Analysis) was performed to quantify differences among treatments in lipid and FA data, while PERMANOVA (Permutational Multivariate Analysis of Variance) was used to perform pairwise tests between treatments (9,999 permutations). In all cases, the non-parametric Bray-Curtis similarity was used.

Pearson correlation analysis was performed to identify relationships between transcript expression (i.e., qPCR RQ data) and tissue lipid composition (i.e., % FA and lipid classes) in the head kidney, using individual fish. All GOIs and lipid classes [i.e., triacylglycerols (TAG), sterols (ST) and phospholipids (PL)] were used in the correlation analysis. However, transcripts and lipid classes with no significant correlations were removed. Individual saturated fatty acids (SFA) and monounsaturated fatty acids (MUFA) that accounted for >5% of the total FA, and all ω6 and ω3 FA were included in the analysis. IBM SPSS (IBM SPSS Statistics, Version 25, Armonk, New York, USA) was used for the correlation analyses. In order to separately group transcripts and lipid composition (% fatty acids and lipid classes), hierarchical clustering was used with group average in PRIMER (Version 6.1.15, Ivybridge, UK).

Results

Head Kidney Lipid Composition

The predominant lipid class in the head kidney was TAG (46.6–52.0%), followed by PL (22.8–28.6%), and ST (18.4–23.7%) (Table 2). No significant differences were observed among dietary treatments in any lipid class or in total lipids (mg g⁻¹) at week 12 (p = 0.09–0.72). ST showed a decrease from week 0 (initial) to week 12 in all treatments (p < 0.0001, Table 2).

TABLE 2

Table 2. Lipid and fatty acid composition (%) of Atlantic salmon head kidney^a before (Initial) and after 12 weeks of feeding diets with different ratios of ω6:ω3 and levels of EPA+DHA.

Individual SFA showed differences among treatments (i.e., 14:0 and 18:0; p < 0.0001 and 0.001, respectively) at week 12. The 1.4% EPA+DHA/balanced treatment had higher 14:0 than the 0.3% EPA+DHA↑ω6, 0.3% EPA+DHA↑ω3 and 1.0% EPA+DHA↑ω6 fed fish. Lower 18:0 levels were observed in the 1.4% EPA+DHA/balanced compared to the high ω6 (i.e., with 0.3% and 1.0% EPA+DHA) treatments. However, 16:0 and total SFA (ΣSFA) were similar among dietary treatments (p = 0.60 and 0.45, respectively) at week 12 (Table 2).

Salmon fed the 1.4% EPA+DHA/balanced diet had the highest level of total MUFA (ΣMUFA), while the 1% EPA+DHA↑ω3 fed fish showed higher ΣMUFA than the 0.3% EPA+DHA↑ω6 fish at week 12 (p < 0.0001). These differences were also observed in some individual MUFA (i.e., 18:1ω7, 20:1ω11, 20:1ω9, 22:1ω11; p < 0.0001). The 0.3% EPA+DHA↑ω3 fish had higher 18:1ω9 compared to fish fed the high ω6 diets (i.e., 0.3 and 1% EPA+DHA), and this FA was higher in the 1.4% EPA+DHA/balanced than the other treatments (not including 0.3% EPA+DHA↑ω3; p < 0.0001). All treatments showed a significant decrease in ΣMUFA, 18:1ω9 and other MUFA (i.e., 16:1ω7, 18:1ω7, 20:1ω9, 22:1ω11) when comparing week 12 to week 0 (p < 0.0001).

The levels of LC-PUFA precursors 18:2ω6 and 18:3ω3 reflected experimental diet levels (i.e., highest 18:2ω6 in the high ω6 fed salmon; highest 18:3ω3 in the high ω3 fed salmon; p < 0.0001; Table 2). Further, compared with week 0, 18:2ω6 and 18:3ω3 increased after 12 weeks of feeding high ω6 (i.e., 0.3% EPA+DHA↑ω6 and 1% EPA+DHA↑ω6) and high ω3 (0.3% EPA+DHA↑ω3 and 1% EPA+DHA↑ω3) diets, respectively (p < 0.0001; Table 2). The intermediate ω6 (18:3ω6, 20:2ω6 and 20:3ω6) and ω3 PUFA (18:4ω3, 20:3ω3 and 20:4ω3), and the ω6:ω3 ratios were different among treatments at week 12 (p < 0.0001) and responded to dietary changes in ω6 and ω3 PUFA. The proportion of 20:4ω6 was highest in fish fed 0.3% EPA+DHA↑ω6 diet (p < 0.0001), while no significant differences were observed among the remaining four treatments at week 12. In contrast, the 0.3% EPA+DHA↑ω6 fish showed the lowest 20:5ω3 (p < 0.0001), and similar levels of this ω3 LC-PUFA were observed among the remaining four treatments (p = 0.18) at week 12. Further, after 12 weeks of feeding, no differences were shown among dietary treatments in 22:6ω3 (p = 0.53). Compared with week 0, 20:4ω6 showed an increase after 12 weeks of feeding high ω6 (i.e., 0.3 and 1% EPA+DHA) and the 0.3% EPA+DHA↑ω3 diet, while 20:5ω3 increased in the 1% EPA+DHA↑ω3 and 1.4% EPA+DHA/balanced fed fish. All 5 treatments increased in 22:6ω3 at week 12 compared to week 0 (p < 0.0001). Fish fed the 0.3% EPA+DHA↑ω6 diet had the highest ω6:ω3, and the high ω6 fish (i.e., 0.3 and 1% EPA+DHA) had higher ω6:ω3 compared to high ω3 (i.e., 0.3 and 1% EPA+DHA) and 1.4% EPA+DHA/balanced fed fish at week 12 (p < 0.0001). The DHA:EPA ratio was higher in the 0.3% EPA+DHA↑ω6 fed fish compared with the other four treatments at week 12 (p < 0.0001; Table 2). Finally, EPA:ARA was highest in the 1% EPA+DHA↑ω3 and the 1.4% EPA+DHA/balanced treatments, and lowest in the 0.3% EPA+DHA↑ω6 fed fish (p < 0.0001).

PCoA showed a clear separation between the high ω6 and high ω3 treatments (i.e., 0.3% EPA+DHA↑ω6 and 1% EPA+DHA↑ω6 on top, and 0.3% EPA+DHA↑ω3 and 1% EPA+DHA↑ω3 on the bottom; Figure 1). The 1.4% EPA+DHA/balanced fed fish clustered with the high ω3 treatments. PCO1 and PCO2 accounted for 49.8 and 30.3% of the variability, respectively. PERMANOVA pairwise tests showed that the 0.3% EPA+DHA↑ω3 were not different from 1% EPA+DHA↑ω3 fed fish [p(perm) = 0.12]. However, all other treatments were significantly different from each other [p(perm) = 0.0001–0.003] (Supplementary Table 3). PCoA vectors showed that fish fed high ω6 diets (particularly 0.3% EPA+DHA↑ω6) were associated with 18:2ω6, 20:2ω6, 20:3ω6, 20:4ω6, and Σω6, while the high ω3 and 1.4% EPA+DHA/balanced diet fed fish were associated with 18:3ω3, 20:3ω3, 20:4ω3, 20:5ω3, Σω3, EPA:ARA, the MUFA 18:1ω7, 16:1ω7, 18:1ω9 and 20:1ω9, as well as ΣMUFA (Figure 1). The lipid classes ST and PL clustered with the high ω6 and high ω3 (including 1.4% EPA+DHA/balanced) treatments, respectively.

FIGURE 1

Figure 1. Principal coordinates analysis (PCoA) of lipid and FA composition (%) (r² > 0.45) in the head kidney of Atlantic salmon after 12 weeks of feeding diets with different ratios of ω6:ω3 and levels of EPA+DHA. PL, ST, and TAG represent the lipid classes phospholipid, sterol and triacylglycerol, respectively. ΣMUFA and EPA:ARA represent total monounsaturated fatty acids and the ratio of eicosapentaenoic (20:5ω3) to arachidonic (20:4ω6) acid, respectively.

SIMPER analysis indicated that the dissimilarity between 0.3% EPA+DHA↑ω6 and 0.3% EPA+DHA↑ω3, and between 0.3% EPA+DHA↑ω6 and 1% EPA+DHA↑ω6 were mainly caused by Σω6, Σω3, and the lipid classes TAG, PL and ST (Supplementary Table 4). The main FA contributors to the dissimilarities between the remaining high ω6 and high ω3 comparisons were Σω6 and 18:2ω6. Finally, the FA dissimilarity between 1.4% EPA+DHA/balanced and the high ω6 treatments was mainly driven by Σω6, ΣPUFA and ΣMUFA, while ΣPUFA contributed to the FA dissimilarity between 1.4% EPA+DHA/balanced and the high ω3 treatments (Supplementary Table 4). This analysis also showed that the highest dissimilarity in lipid composition profile was between 0.3% EPA+DHA↑ω6 and 1.4% EPA+DHA/balanced fed fish (average dissimilarity = 18.7%), while the lowest dissimilarity was between the two high ω3 (i.e., 0.3% EPA+DHA↑ω3 and 1% EPA+DHA↑ω3; average dissimilarity = 10.0%) and the two high ω6 (i.e., 0.3% EPA+DHA↑ω6 and 1% EPA+DHA↑ω6; average dissimilarity = 11.0%) treatments.

QPCR Analysis of Transcripts Involved in Fatty Acid Metabolism and Eicosanoid Synthesis

Dietary Impacts on Transcript Expression

One-way ANOVA indicated that salmon fed 0.3% EPA+DHA↑ω3, 1.4% EPA+DHA/balanced, and the 1% EPA+DHA diets had higher expression levels (1.31–1.40-fold; p = 0.01) of the transcript elovl5a compared with the 0.3% EPA+DHA↑ω6 fed fish (Figure 2B and Supplementary Table 5). The head kidney transcript expression of fadsd5 was higher in the 1% EPA+DHA↑ω3 compared with the 0.3% EPA+DHA↑ω6 fed fish (1.62-fold; p = 0.04; Figure 2C and Supplementary Table 5), while srebp1 mRNA levels were higher in the 1% EPA+DHA↑ω3 and 1.4% EPA+DHA/balanced treatments compared with the 0.3% EPA+DHA↑ω6 fed fish (1.48- and 1.47-fold, respectively; p = 0.02; Figure 2D and Supplementary Table 5). The transcript pparaa showed lower expression in the 0.3% EPA+DHA↑ω3 compared with the 0.3% EPA+DHA↑ω6 fed fish (0.37-fold; p = 0.02; Figure 2G and Supplementary Table 5). The prostaglandin synthesis-related transcript ptges3 showed an increasing trend (although not significant) in the 0.3% EPA+DHA↑ω3 compared with the 0.3% EPA+DHA↑ω6 fed fish (1.30-fold; p = 0.06; Supplementary Table 5). Finally, the transcript expression of elovl4b, srebp2 and lxrb (Figures 2A,E,F, respectively, and Supplementary Table 5) was not significantly impacted by diet.

FIGURE 2

Figure 2. Head kidney qPCR transcripts related to fatty acid metabolism (A–G) and eicosanoid synthesis (H) in salmon fed diets with different ratios of ω6:ω3 and levels of EPA+DHA for 12 weeks. Transcripts identified as statistically significant in the two-way ANOVA analysis are presented (A–I). Transcript expression data of these and other genes of interest (GOIs) are presented in Supplementary Table 5. Transcript expression values presented as mean relative quantity (RQ) ± SE (n = 7–8). Different letters above error bars indicate significant differences among treatments. When differences were not statistically significant (one-way ANOVA and Tukey post-hoc tests; p > 0.05), p-values are shown above error bars. Values below each GOI represent fold-changes relative to the 0.3% EPA+DHA↑ω6 treatment (see Materials and Methods). (I) Two-way ANOVA results. Significant differences (p < 0.05) are shown in bold font. Two-way ANOVA analysis of the entire set of GOIs is presented in Supplementary Table 6.

Effects of EPA+DHA and ω6:ω3 Factors on Transcript Expression

Two-way ANOVA analysis showed a significant interaction between the dietary factor EPA+DHA and ω6:ω3, impacting the transcript expression of elovl4b (Figure 2I and Supplementary Table 6). The transcript expression of elovl5a was significantly affected by both EPA+DHA and ω6:ω3, while fadsd5, srebp1, srebp2 and lxrb were significantly affected by dietary EPA+DHA levels. Finally, the transcript expression of pparaa and cox1 was significantly affected by the factor ω6:ω3 (Figure 2I and Supplementary Table 6).

Hierarchical Clustering and Pearson Correlation Analyses

Hierarchical Clustering of Head Kidney Transcript Expression and Lipid Composition

Hierarchical clustering of the qPCR analysed transcripts identified three separate clusters (Figure 3). Cluster I consisted of the transcript pparb2a, and the eicosanoid metabolism-related transcripts ptges3, ptges2, 5loxa, 5loxb, and cox1. Cluster II included pparaa and lxra, while cluster III included the lipid metabolism-related transcripts srebp2, elovl5a, srebp1, acox1, fadsd5, and the eicosanoid synthesis-related cox2 and pgds. Cluster analysis of head kidney lipid composition also showed three clusters. Cluster I was composed of the ω3 FA (i.e., 20:5ω3, 22:6ω3, 20:4ω3, 18:4ω3, 18:3ω3, 20:3ω3, and 22:5ω3), Σω3, 18:1ω9, ΣMUFA, EPA:ARA, and the lipid class TAG. Cluster II included 18:0, 16:0, ΣSFA and 22:4ω6. Cluster III showed the ω6 FA (i.e., 18:3ω6, 20:2ω6, 18:2ω6, 22:5ω6, 20:3ω6, 20:4ω6), Σω6, ω6:ω3, as well as ΣPUFA and the lipid class ST (Figure 3).

FIGURE 3

Figure 3. Hierarchical clustering and Pearson correlation matrix of head kidney transcript expression [qPCR relative quantity (RQ) values] and lipid composition (% fatty acids and lipid classes) in Atlantic salmon fed diets with different ratios of ω6:ω3 and levels of EPA+DHA for 12 weeks. Transcripts, fatty acids, and lipid classes that were not significantly correlated were excluded (see Materials and Methods). Statistically significant (p < 0.05) correlation coefficients are shown. Red cells signify negative relationships, and green cells signify positive relationships. ΣSFA, ΣMUFA, and ΣPUFA represent total saturated, monounsaturated, and polyunsaturated fatty acids, respectively. 20:5ω3, 22:5ω3, 22:6ω3, and 20:4ω6 represent EPA, DPA, DHA, and ARA, respectively. TAG and ST represent triacylglycerol and sterols, respectively.

Correlations Between Lipid Composition and Fatty Acid Metabolism Transcripts

The head kidney transcript expression of pparb2a was correlated positively with 20:4ω3 (p = 0.031; Figure 3), while pparaa showed negative correlations with several ω3 FA (i.e., 18:3ω3, 18:4ω3, 20:3ω3, 20:4ω3) and Σω3 (p = 0.002–0.044), and positive correlations with several ω6 FA (i.e., 18:2ω6, 20:2ω6, 20:3ω6, 20:4ω6, 22:5ω6), Σω6, and the ratio of ω6:ω3 (p = 0.009–0.038). The five transcripts lxra, srebp2, elovl5a, srebp1, and fadsd5 showed positive correlations with EPA:ARA (p = 0.001–0.032; Figure 3). The mRNA levels of elovl5a also showed positive correlations with 20:3ω3, 20:4ω3 and ΣMUFA (p = 0.013–0.044), while srebp1 correlated positively with 22:5ω3 and ΣMUFA and ΣSFA (p = 0.006–0.047). The transcript fadsd5 correlated positively with ΣMUFA, ΣSFA, and the FA 16:0 and 18:0 (p = 0.001–0.032). The same five transcripts (i.e., lxra, srebp2, elovl5a, srebp1, and fadsd5) and acox1 correlated negatively with ω6 FA; lxra with 20:2ω6, 20:3ω6 and Σω6 (p = 0.021–0.042), and srebp2 with 18:3ω6, 20:3ω6, 20:4ω6, 22:5ω6, as well as Σω6 and ΣPUFA (p = 0.0001–0.041). Elovl5a, srebp1, and fadsd5 correlated negatively with 18:2ω6, 18:3ω6, 20:2ω6, 20:3ω6, 20:4ω6, 22:5ω6, Σω6, ω6:ω3, and ΣPUFA (p = 0.0001–0.028). The transcript acox1 showed negative correlations with 18:3ω6 and 22:5ω6 (p = 0.011 and 0.016, respectively; Figure 3).

Correlations Between Lipid Composition and Eicosanoid Metabolism Transcripts

The transcript expression of ptges3 was correlated positively with 18:4ω3, 20:4ω3 and 18:0 (p = 0.012–0.049; Figure 3), while ptges2 showed a negative correlation with 20:5ω3 (p = 0.041). The transcript expressions of 5loxa and 5loxb showed positive correlations with 20:5ω3 and 22:6ω3 (p = 0.002-0.045). 5loxb also correlated positively with 20:4ω3 and Σω3 (p = 0.042 and 0.017, respectively), and negatively with ω6:ω3 (p = 0.044), while 5loxa correlated positively with ST and 20:4ω6 (p = 0.012 and 0.006, respectively), and negatively with 18:4ω3, TAG, 18:1ω9 and ΣMUFA (p = 0.006-0.038; Figure 3). The mRNA levels of cox1 showed positive correlations with the ω3 FA: 18:3ω3, 18:4ω3, 20:3ω3, 20:4ω3 and Σω3 (p = 0.009–0.043). The transcript cox2 correlated positively with 16:0, ΣSFA and 22:4ω6 (p = 0.007–0.029), and negatively with 18:3ω6, 22:5ω6 and ΣPUFA (p = 0.013–0.037). Finally, pgds correlated negatively with ΣPUFA (p = 0.045; Figure 3).

Discussion

Head Kidney Lipid Composition

The head kidney was mainly composed of TAG (46.6–52.0%), and this was followed by the lipid classes PL (22.8–28.6%), and ST (18.4–23.7%) (Table 2). Our lipid class (i.e., TAG and ST) proportions are in the range of what was found in earlier studies on Atlantic salmon head kidney (Martinez-Rubio et al., 2013; Foroutani et al., 2020). Similar to our findings, Foroutani et al. (2020) reported that TAG was the predominant lipid class in Atlantic salmon head kidney. Head kidney is a principal component of the immune system and the major site of haematopoiesis in fish (Tort et al., 2003; Zapata et al., 2006). Therefore, high levels of TAG may allow fish to fulfill the energy required to fuel cell differentiation processes and other immune functions. Interestingly, Thompson et al. (1995) reported that TAG proportions in Atlantic salmon head kidney leucocytes, macrophages, T- and B-cells varied between ~14 and 20%, and indicated that TAG may serve as an energy store and/or a source of free fatty acids that are used in cellular functions in salmon. Our study also showed that lipid class composition [when expressed as % or as concentration (mg g⁻¹ wet weight)] was not significantly affected by diet (Table 2, Supplementary Table 7). This result is in agreement with Betancor at al. (2014) who reported that dietary variation in LC-PUFA did not impact head kidney lipid classes in Atlantic salmon. However, the current study revealed that PL proportions were correlated positively with 20:5ω3 and negatively with 18:2ω6, 20:2ω6, and ω6:ω3 in the head kidney (Supplementary Table 8). Further, PCoA indicated associations between PL and the ω3 PUFA 22:6ω3, 20:5ω3, and Σω3 (Figure 1). These data suggest a preferential incorporation into the membrane of ω3 over ω6 PUFA. This can potentially increase the availability of anti-inflammatory eicosanoids in the head kidney (Martinez-Rubio et al., 2013). However, this can only be postulated, as fatty acid composition of the PL portion of the head kidney was not determined in our study.

The fatty acid composition of the head kidney reflected that of the experimental diets. For example, LC-PUFA precursor composition (i.e., 18:2ω6 and 18:3ω3) and ω6:ω3 ratio of the head kidney reflected the dietary composition (i.e., highest 18:2ω6 and ω6:ω3 in the high ω6 fed fish; highest 18:3ω3 in the high ω3 fed fish). Previous studies also indicated that head kidney lipid composition of Atlantic salmon (Gjøen et al., 2004; Foroutani et al., 2020), gilthead seabream (Sparus aurata) (Montero et al., 2003) and Nile tilapia (Oreochromis niloticus) (Chen et al., 2016) reflected dietary variation in 18:2ω6 and 18:3ω3. Gjøen et al. (2004) demonstrated that dietary replacement of FO with soy oil (i.e., 50%) modified head kidney ω6:ω3 (from 0.3 to 1.1), which reflected the dietary ω6:ω3. These ratios are comparable to the ω6:ω3 detected in our study (i.e., 0.5–1.5). At week 12, the 0.3% EPA+DHA↑ω6 fed fish had 1.7–2.3-fold more 20:4ω6 in terms of proportions and 1.4–1.8-fold in terms of concentrations, when compared with the other treatments. However, 20:5ω3 proportions (0.6–0.7-fold) and concentrations (0.4–0.6-fold) were lowest in this group, compared to the other treatments. Additionally, both proportions and concentrations of 20:5ω3 were similar between the 0.3% EPA+DHA↑ω3 and the 1% EPA+DHA↑ω3 fed fish. A similar effect on head kidney 20:4ω6 proportions was shown in a previous study where soy oil fed Atlantic salmon had 3-fold higher 20:4ω6 than FO fed fish (Gjøen et al., 2004). Furthermore, the differences in 20:4ω6 and 20:5ω3 proportions are in line with our previous study with Atlantic salmon muscle (Emam et al., 2020). Taken together, these data highlight the importance of dietary LC-PUFA precursors (i.e., 18:2ω6 and 18:3ω3) when salmon are fed with lower inclusion of EPA and DHA (0.3%), and their impact on head kidney LC-PUFA. Although LC-PUFA synthesis was not quantified in our study, the fact that intermediate ω6 and ω3 PUFAs were also significantly different among dietary treatments, and their head kidney proportions were higher than those of the diet, suggest that the LC-PUFA synthesis pathway was activated. Interestingly, 22:6ω3 proportions were not affected by diet at week 12, and 22:6ω3 concentrations (within the high ω3 and high ω6 treatments) did not vary with increasing dietary EPA+DHA levels. This suggests that in addition to elongation and desaturation, selective retention had also occurred. Indeed, earlier studies indicated higher 22:6ω3 retention in Atlantic salmon with decreasing dietary EPA+DHA (Glencross et al., 2014; Bou et al., 2017). However, the fact that 1% EPA+DHA/high↑ω3 and 1.4% EPA+DHA/balanced treatments showed higher concentration of 22:6ω3 than the 0.3% EPA+DHA/high↑ω6 fed fish suggests that both dietary EPA+DHA and ω6:ω3 influenced head kidney DHA content.

Salmon fed the 1.4% EPA+DHA/balanced diet had the highest level of ΣMUFA, and the 1% EPA+DHA↑ω3 fed fish showed higher ΣMUFA than the 0.3% EPA+DHA↑ω6 fed fish. These differences in ΣMUFA, and the observed changes in individual MUFA (i.e., 16:1ω7, 18:1ω7, 18:1ω9 and 20:1ω9) were also a reflection of the diets. Notably, the 1.4% EPA+DHA/balanced diet contained rapeseed oil which is rich in 18:1ω9 (Orsavova et al., 2015), and this may have influenced the observed changes in tissue MUFA. Emam et al. (2020) reported a similar response in salmon muscle, and showed that hepatic MUFA synthesis-related transcript scdb was positively correlated with muscle 20:1ω7 and 20:1ω9, among other MUFA. This suggests that head kidney MUFA differences between treatments in our study were also related to changes in MUFA synthesis. However, further analyses are required in order to test this hypothesis (e.g., transcript expression of MUFA synthesis-related genes in the head kidney).

Transcript Expression of Lipid Metabolism-Related Genes

Transcripts Involved in LC-PUFA Synthesis

Transcript expression levels of key genes encoding enzymes involved in LC-PUFA synthesis were modulated by diet in our study. Salmon fed with 0.3% EPA+DHA↑ω6 had the lowest expression of elovl5a, while the remaining four treatments had similar levels of this transcript (Figure 2B and Supplementary Table 5). The transcript expression of other genes encoding elongases (elovl4b) and desaturases (fadsd5, fadsd6b and fadsd6c) was also similar among 0.3% EPA+DHA↑ω3, 1% EPA+DHA (both high ω3 and high ω6) and 1.4% EPA+DHA/balanced fed fish (Supplementary Table 5). Further, elovl5a was correlated positively with ω3 (20:3ω3 and 20:4ω3), while both elovl5a and fadsd5 were correlated negatively with ω6 PUFA (18:2ω6, 18:3ω6, 20:2ω6, 20:3ω6, 20:4ω6, 22:5ω6) in the head kidney. These transcript expression changes along with the fatty acid data suggest that 20:5ω3 synthesis was lower in fish fed lower dietary levels of 18:3ω3 and EPA+DHA (i.e., 0.3%), and further support the notion that higher dietary 18:3ω3 promoted the synthesis of ω3 LC-PUFA in the head kidney. Indeed, LC-PUFA synthesis in salmonids is influenced by the availability of tissue LC-PUFA (i.e., EPA, DHA, ARA) and the precursors 18:3ω3 and 18:2ω6 (Jordal et al., 2005; Jump et al., 2005; Glencross et al., 2015; Katan et al., 2019). Fadsd5 expression was higher with 1% EPA+DHA↑ω3 compared with 0.3% EPA+DHA↑ω6 fed fish, while srebp1 showed higher expression levels in the 1% EPA+DHA↑ω3 and 1.4% EPA+DHA/balanced treatments compared with the 0.3% EPA+DHA↑ω6 fed fish. These transcripts, as well as srebp2 and lxrb, were significantly positively impacted by dietary EPA+DHA levels (Figures 2C,D,I). Our results are in contrast with previous studies reporting up-regulation of hepatic fadsd5 in Atlantic salmon fed plant-based (containing lower EPA and DHA) compared to FO diets (Zheng et al., 2004; Morais et al., 2011; Xue et al., 2015). However, our data is in agreement with Betancor et al. (2014) who reported up-regulation of head kidney srebp1 and a ~1.5-fold increase in fadsd5 in Atlantic salmon fed higher dietary levels of EPA+DHA (2% compared to 1% of diet). FO diets used in the former studies contained higher proportions of EPA+DHA (~11–24% of total FA) and higher inclusions of FO (~14–30% of diet) compared to our study (EPA+DHA proportions and FO inclusions were 2.8–8.8% and 0.1–6.8%, respectively). These as well as organ (i.e., liver vs. head kidney) differences may explain the discrepancies between studies. Further, dietary ratios of ALA:LNA were shown to influence LC-PUFA accumulation in salmonids and African catfish muscle tissues (Colombo et al., 2018; Sourabié et al., 2019), and this may have played a role in the regulation of LC-PUFA synthesis in our study. Finally, the positive impact of dietary EPA+DHA levels in our study could be related to MUFA levels, as the transcript expression of several lipid metabolism genes (elovl5a, fadsd5, and srebp1) showed a positive correlation with head kidney ΣMUFA (Figure 3). Interestingly, Emam et al. (2020) showed that hepatic elovl5a and srebp2 were positively correlated with the levels of individual MUFA in salmon muscle. Indeed, an earlier study with Atlantic salmon revealed that LC-PUFA synthesis is stimulated by dietary MUFA levels, and suggested that the sparing effect of SFA and MUFA on ω3 LC-PUFA could be the underlying mechanism (Emery et al., 2016).

The transcript expression of head kidney pparaa was higher in the 0.3% EPA+DHA↑ω6 compared with the 0.3% EPA+DHA↑ω3 fed fish (Figure 2G), and correlation analysis indicated that this transcript was positively correlated with ω6 PUFA (18:2ω6, 20:2ω6, 20:3ω6, 20:4ω6, 22:5ω6), Σω6 and ω6:ω3, and negatively with ω3 PUFA (18:3ω3, 18:4ω3, 20:3ω3, 20:4ω3) and Σω3 (Figure 3). Ppara is a transcription factor that regulates numerous genes involved in lipid-related processes such as FA oxidation, as well as bile acid and triacylglycerol metabolism (Kersten, 2014). Further, several studies in fish demonstrated involvement of Ppara in LC-PUFA synthesis (Dong et al., 2017; Zhu et al., 2019). Thus, our data suggest that pparaa may have played a regulatory role in ω6 PUFA synthesis, particularly in the 0.3% EPA+DHA↑ω6 fed fish. Emam et al. (2020) showed that hepatic pparaa expression was higher in the 0.3% EPA+DHA↑ω3 than the 0.3% EPA+DHA↑ω6 and 1% EPA+DHA treatments (i.e., high ω3 and high ω6), and also a positive correlation between liver pparaa and muscle ω3 PUFA. These contrasting results suggest that pparaa response in Atlantic salmon is tissue dependent. Finally, the fact that the transcript expression of pparaa was significantly impacted by dietary ω6:ω3 (Figure 2I), and the correlation with head kidney ω6:ω3 (Figure 3), suggest that this transcript may have been involved with inflammatory processes in our study. This is in line with previous studies with mammals (Jones et al., 2002; Becker et al., 2008; Varga et al., 2011) and Atlantic salmon (Martinez-Rubio et al., 2013) which showed that Ppara and ppara play a role in the resolution of inflammation. More studies are required in order to investigate the dietary-induced impacts of ppara on inflammatory processes in Atlantic salmon head kidney.

Transcripts Involved in Eicosanoid Metabolism

Our qPCR results (Supplementary Table 5) are in line with previous studies which indicated no differences in head kidney constitutive transcript expression of several genes related to eicosanoid synthesis (i.e., 5lox, cox2) in fish exposed to varying EPA levels [i.e., 0–20 μM (Furne et al., 2013), 1.1 and 19 mg g⁻¹ lipid (Salini et al., 2016)] or ω6:ω3 ratios (0.7–4.1) (Holen et al., 2018). However, Holen et al. (2018) reported that Atlantic salmon fed with soybean oil as a FO replacement had higher expression of pgds and pges in head kidney leucocytes compared to fish fed palm and rapeseed oils. Differences in dietary and tissue EPA:ARA, EPA+DHA, the tissues/cells examined (head kidney vs. isolated leucocytes), and variation between in vivo and in vitro models could have contributed to these discrepancies between studies. Further, our study revealed that head kidney transcript expression of cox1 was negatively impacted by dietary ω6:ω3 levels (Figures 2H,I), and positively correlated with head kidney ω3 fatty acids (i.e., 18:3ω3, 18:4ω3, 20:3ω3, 20:4ω3, and Σω3; Figure 3). The enzyme COX1 is constitutively expressed in several tissues (e.g., liver, kidney, spleen, gill, muscle, gut) and has maintenance and homeostatic functions (Tapiero et al., 2002; Olsen et al., 2012). Our data suggest that cox1 expression could result in production of anti-inflammatory prostaglandins in Atlantic salmon head kidney at the constitutive level. Further, differences among treatments in head kidney EPA:ARA (Table 2) may have altered the production of eicosanoids via cox1 transcription. Indeed, changes in cell membrane EPA:ARA could alter eicosanoid production and the inflammatory response in vertebrates (Calder, 2003; Arts and Kohler, 2009; Salini et al., 2016). However, since eicosanoid levels were not measured in our study, these hypotheses could not be tested. Further, the implications of FA-transcript correlations on eicosanoid production could not be assessed.

In addition to cox1, transcript expression of 5loxa, 5loxb, and ptges3 was also positively correlated with head kidney ω3 LC-PUFA. However, 5loxa showed positive correlations with both ω3 (i.e., 20:5ω3, 22:6ω3) and ω6 (20:4ω6) LC-PUFA. The latter result is in line with Katan et al. (2019) who reported that hepatic 5loxa showed a positive relationship (p = 0.059) with liver EPA+ARA in Atlantic salmon fed different mixes of plant oils as a FO replacement. However, this is in contrast with Caballero-Solares et al. (2020) who showed that hepatic 5loxa correlated negatively with liver EPA and DHA in salmon fed terrestrial plant ingredients as FO and FM alternatives. Discrepancies between studies suggest that the interaction of 5loxa transcription with tissue LC-PUFA in salmon depends on tissue (i.e., head kidney vs. liver) and/or dietary inputs (e.g., protein and lipid sources). Our study also showed positive correlations between ω3 PUFA (i.e., 18:4ω3, 20:4ω3) and ptges3 expression, and negative correlations between EPA and ptges2. In mammals, both PTGES2 and PTGES3 convert PGH₂ into PGE₂ (Neuman et al., 2017; Xu et al., 2019). Wei et al. (2016) showed that muscle ptges3 expression was induced in pigs fed with linseed-enriched diets for 60 days. However, previous studies reported no impacts of dietary fatty acids on the transcript expression of ptges2 and ptges3 in bovine oocytes (Ponter et al., 2012) and Eurasian perch (Perca fluviatilis) liver, brain and intestine (Geay et al., 2015). Clearly, more research is required to elucidate the influence of dietary and tissue PUFA on the transcript expression of ptges2 and ptges3 in fish.

Cox2 showed a positive correlation with 16:0, ΣSFA and 22:4ω6, and negative correlations with ΣPUFA, 18:3ω6 and 22:5ω6 in the current study (Figure 3). This finding supports previous human macrophage studies which indicated that SFA induced cox2 expression (Lee et al., 2001; Rocha et al., 2016). Astiz et al. (2012) showed that addition of 22:4ω6 and 22:5ω6 did not impact the production of COX2 in rat testis cells. However, a previous study on mouse brain indicated that reduction in the levels of ARA metabolites 22:4ω6 and 22:5ω6 reflected a shift of ARA from LC-PUFA biosynthesis prostaglandin or leukotriene synthesis (McNamara et al., 2009). Thus, the negative correlation between 22:5ω6 and cox2 could be a result of higher availability of ARA for COX-2 mediated prostaglandin synthesis. Finally, similarly to cox2, pgds correlated negatively with ΣPUFA. This latter finding is in line with Katan et al. (2019) who showed that hepatic pgds was negatively correlated with liver EPA+ARA+DHA in Atlantic salmon. However, Holen et al. (2015) reported that the combination of EPA+ARA+DHA did not modulate cox2 transcript expression in Atlantic salmon head kidney cells. In summary, further examination is needed in order to elucidate the influence of ARA metabolites (i.e., 22:4ω6 and 22:5ω6) on cox2 expression, and the interactive impacts of PUFA on the transcript expression of pgds and cox2 in Atlantic salmon head kidney.

Conclusion

This study investigated head kidney lipid composition and transcript expression of genes involved in fatty acid and eicosanoid metabolism in Atlantic salmon fed varying dietary levels of EPA+DHA (0.3, 1.0, and 1.4%) and ω6:ω3 ratios (high ω6, high ω3, and balanced) for 12 weeks. TAG was the predominant lipid class in all treatments, regardless of diet. This suggests that this lipid class could play an important role in immune and/differentiation processes in Atlantic salmon head kidney. Head kidney fatty acid composition was reflective of the diet with respect to C₁₈ PUFA and MUFA levels (% of total fatty acids), and responded to ω6:ω3 variation. Proportions of 20:5ω3 were similar among 0.3% EPA+DHA↑ω3, 1% EPA+DHA (both high ω3 and high ω6 treatments) and 1.4% EPA+DHA/balanced fed fish, although dietary 20:5ω3 varied by 2.5- to 3-fold. Tissue fatty acid composition changes agreed with positive correlations between head kidney 20:3ω3 and 20:4ω3 and elovl5a transcript levels. This suggested that high dietary 18:3ω3 promoted the synthesis of ω3 LC-PUFA in salmon fed lower dietary EPA+DHA levels (0.3%). Further, head kidney transcript expression of several genes involved in FA metabolism (e.g., elovl5a, fadsd5 and srebp1) was positively impacted by dietary EPA+DHA levels, and showed a positive correlation with head kidney ΣMUFA. This supported the hypothesis that LC-PUFA synthesis is stimulated by increasing MUFA levels in Atlantic salmon.

Several eicosanoid synthesis-related transcripts were significantly correlated with head kidney fatty acid composition. The head kidney transcript expression of 5loxb correlated positively with ω3 LC-PUFA and Σω3, and correlated negatively with ω6:ω3. Cox1 and ptges3 transcript expression correlated positively with ω3 PUFA, while cox1 was negatively impacted by dietary ω6:ω3 levels. We hypothesize that these transcript expression changes may have increased the production of anti-inflammatory eicosanoids in salmon head kidney. Finally, our study revealed significant positive and negative correlations of cox2 expression with SFA (i.e., 16:0, ΣSFA) and 22:5ω6, respectively. Given the important role of the head kidney in salmon immune responses and the interaction with dietary LC-PUFA, future studies should elucidate how ω6:ω3 and EPA+DHA levels affect LC-PUFA synthesis (e.g., stable isotope and FA mass balance methods) and eicosanoid production (e.g., LC-MS/MS) in this organ.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The animal study was reviewed and approved by Animal Care Committee of Memorial University of Newfoundland.

Author Contributions

MR, CP, and RT contributed to conception and design of the study. TK was involved in all aspects of the study, including experimental design, data and statistical analyses, and writing of the first manuscript draft. TK, XX, and AC-S took part in sampling. XX assisted with analyses related to the qPCR study. AC-S assisted with the Pearson correlation matrix. All authors contributed to original manuscript editing, revisions, and they read and approved the submitted version.

Funding

This study was part of the Integrated Pathogen Management of Co-infection in Atlantic Salmon (IPMC) project (Genomic Applications Partnership Program, GAPP #6607) which was funded by the Government of Canada through Genome Canada and Genome Atlantic. The IPMC project was also funded by InnovateNL (Government of Newfoundland and Labrador Department of Tourism, Culture, Industry and Innovation; Leverage R&D award #5401-1019-108), and Cargill Innovation (formerly EWOS Innovation). Furthermore, MR was supported by a Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant (341304-2012), and TK was supported by a Post-graduate Scholarship-Doctoral (PGS D) from NSERC, Ocean Industries Student Research Award (OISRA) from Innovate NL, and SGS fellowship from Memorial University Newfoundland.

Conflict of Interest

RT was employed by the company Cargill, but did not have any additional role in the study design, data collection, analysis and interpretation, writing of the manuscript, or in the decision to submit the manuscript for publication.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We would like to acknowledge Cara Kirkpatrick (Genome Atlantic) for managing the project, Dr. Dominic Nanton (Cargill) for feed formulations, Chute Animal Nutrition Center team (Dalhousie University) for feed processing, and Dr. Joe Brown Aquatic Research Building (JBARB) staff for their assistance with fish husbandry and sampling. We would like to thank Jeanette Wells for training, as well as Maryam Beheshti Foroutani, Kathleen S. Parrish and Hoda Rajabi for their assistance with the laboratory work. We are also grateful to the sampling teams in the MR and CP laboratories.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2020.602587/full#supplementary-material

References

Arts, M. T., and Kohler, C. C. (2009). “Health and condition in fish: the influence of lipids on membrane competency and immune response,” in Lipids in Aquatic Ecosystems, eds M. Kainz, M. Brett, M. Arts (New York, NY: Springer), 237–256. doi: 10.1007/978-0-387-89366-2_10

CrossRef Full Text | Google Scholar

Astiz, M., De Catalfo, G. H., De Alaniz, M. J. T., and Marra, C. A. (2012). Exogenous arachidonate restores the dimethoate-induced inhibition of steroidogenesis in rat interstitial cells. Lipids 47, 557–569. doi: 10.1007/s11745-012-3669-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Becker, J., Delayre-Orthez, C., Frossard, N., and Pons, F. (2008). Regulation of peroxisome proliferator-activated receptor-α expression during lung inflammation. Pulm. Pharmacol. Ther. 21, 324–330. doi: 10.1016/j.pupt.2007.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Betancor, M. B., Howarth, F. J. E., Glencross, B. D., and Tocher, D. R. (2014). Influence of dietary docosahexaenoic acid in combination with other long-chain polyunsaturated fatty acids on expression of biosynthesis genes and phospholipid fatty acid compositions in tissues of post-smolt Atlantic salmon (Salmo salar). Comp. Biochem. Physiol. Part B Biochem. Mol. Biol. 172, 74–89. doi: 10.1016/j.cbpb.2014.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Booman, M., Borza, T., Feng, C. Y., Hori, T. S., Higgins, B., Culf, A., et al. (2011). Development and experimental validation of a 20K Atlantic cod (Gadus morhua) oligonucleotide microarray based on a collection of over 150,000 ESTs. Mar. Biotechnol. 13, 733–750. doi: 10.1007/s10126-010-9335-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Booman, M., Xu, Q., and Rise, M. L. (2014). Evaluation of the impact of camelina oil-containing diets on the expression of genes involved in the innate anti-viral immune response in Atlantic cod (Gadus morhua). Fish Shellfish Immunol. 41, 52–63. doi: 10.1016/j.fsi.2014.05.017

CrossRef Full Text | Google Scholar

Bou, M., Berge, G. M., Baeverfjord, G., Sigholt, T., Østbye, T.-K., Romarheim, O. H., et al. (2017). Requirements of n-3 very long-chain PUFA in Atlantic salmon (Salmo salar L): effects of different dietary levels of EPA and DHA on fish performance and tissue composition and integrity. Br. J. Nutr. 117, 30–47. doi: 10.1017/S0007114516004396

CrossRef Full Text | Google Scholar

Caballero-Solares, A., Hall, J. R., Xue, X., Eslamloo, K., Taylor, R. G., Parrish, C. C., et al. (2017). The dietary replacement of marine ingredients by terrestrial animal and plant alternatives modulates the antiviral immune response of Atlantic salmon (Salmo salar). Fish Shellfish Immunol. 64, 24–38. doi: 10.1016/j.fsi.2017.02.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Caballero-Solares, A., Xue, X., Cleveland, B. M., Foroutani, M. B., Parrish, C. C., Taylor, R. G., et al. (2020). Diet-induced physiological responses in the liver of Atlantic salmon (Salmo salar) inferred using multiplex PCR platforms. Mar. Biotechnol. 22, 511–525. doi: 10.1007/s10126-020-09972-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Calder, P. C. (2003). Long-chain n-3 fatty acids and inflammation: potential application in surgical and trauma patients. Brazilian J. Med. Biol. Res. 36, 433–446. doi: 10.1590/S0100-879X2003000400004

PubMed Abstract | CrossRef Full Text | Google Scholar

Calder, P. C. (2007). Immunomodulation by omega-3 fatty acids. Prostaglandins Leukot. Essent. Fat. Acids 77, 327–335. doi: 10.1016/j.plefa.2007.10.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Calder, P. C. (2013). n-3 fatty acids, inflammation and immunity: new mechanisms to explain old actions. Proc. Nutr. Soc. 72, 326–336. doi: 10.1017/S0029665113001031

PubMed Abstract | CrossRef Full Text | Google Scholar

Carmona-Antoñanzas, G., Monroig, Ó., Dick, J. R., Davie, A., and Tocher, D. R. (2011). Biosynthesis of very long-chain fatty acids (C>24) in Atlantic salmon: Cloning, functional characterisation, and tissue distribution of an Elovl4 elongase. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 159, 122–129. doi: 10.1016/j.cbpb.2011.02.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Carmona-Antoñanzas, G., Tocher, D. R., Martinez-Rubio, L., and Leaver, M. J. (2014). Conservation of lipid metabolic gene transcriptional regulatory networks in fish and mammals. Gene 534, 1–9. doi: 10.1016/j.gene.2013.10.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C., Sun, B., Guan, W., Bi, Y., Li, P., Ma, J., et al. (2016). N-3 essential fatty acids in Nile tilapia, Oreochromis niloticus: effects of linolenic acid on non-specific immunity and anti-inflammatory responses in juvenile fish. Aquaculture 450, 250–257. doi: 10.1016/j.aquaculture.2015.08.005

CrossRef Full Text | Google Scholar

Colombo, S. M., Parrish, C. C., and Wijekoon, M. P. A. (2018). Optimizing long chain-polyunsaturated fatty acid synthesis in salmonids by balancing dietary inputs. PLoS ONE 13:e0205347. doi: 10.1371/journal.pone.0205347

PubMed Abstract | CrossRef Full Text | Google Scholar

De Roos, B., Sneddon, A. A., Sprague, M., Horgan, G. W., and Brouwer, I. A. (2017). The potential impact of compositional changes in farmed fish on its health-giving properties: is it time to reconsider current dietary recommendations? Public Health Nutr. 20, 2042–2049. doi: 10.1017/S1368980017000696

CrossRef Full Text | Google Scholar

Dong, X., Tan, P., Cai, Z., Xu, H., Li, J., Ren, W., et al. (2017). Regulation of FADS2 transcription by SREBP-1 and PPAR-α influences LC-PUFA biosynthesis in fish. Sci. Rep. 7, 1–11. doi: 10.1038/srep40024

PubMed Abstract | CrossRef Full Text | Google Scholar

Einen, O., and Thomassen, M. S. (1998). Starvation prior to slaughter in Atlantic salmon. II. Effects during ice storage on white muscle composition, texture and color of raw and sensory assessment of cooked fillet. Aquaculture 169, 37–53. doi: 10.1016/S0044-8486(98)00332-9

CrossRef Full Text | Google Scholar

Emam, M., Katan, T., Caballero-Solares, A., Taylor, R. G., Parrish, K. S., Rise, M. L., et al. (2020). Interaction between ω6 and ω3 fatty acids of different chain lengths regulates Atlantic salmon hepatic gene expression and muscle fatty acid profiles. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 375. doi: 10.1098/rstb.2019.0648

PubMed Abstract | CrossRef Full Text | Google Scholar

Emery, J. A., Norambuena, F., Trushenski, J., and Turchini, G. M. (2016). Uncoupling EPA and DHA in fish nutrition: dietary demand is limited in Atlantic salmon and effectively met by DHA alone. Lipids 51, 399–412. doi: 10.1007/s11745-016-4136-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Foroutani, M. B., Parrish, C. C., Wells, J., Taylor, R. G., and Rise, M. L. (2020). Minimizing marine ingredients in diets of farmed Atlantic salmon (Salmo salar): effects on liver and head kidney lipid class and fatty acid composition. Fish Physiol. Biochem. 46, 2331–2353. doi: 10.1007/s10695-020-00862-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Froehlich, H. E., Jacobsen, N. S., Essington, T. E., Clavelle, T., and Halpern, B. S. (2018). Avoiding the ecological limits of forage fish for fed aquaculture. Nat. Sustain. 1:298. doi: 10.1038/s41893-018-0077-1

CrossRef Full Text | Google Scholar

Furne, M., Holen, E., Araujo, P., Lie, K. K., and Moren, M. (2013). Cytokine gene expression and prostaglandin production in head kidney leukocytes isolated from Atlantic cod (Gadus morhua) added different levels of arachidonic acid and eicosapentaenoic acid. Fish Shellfish Immunol. 34, 770–777. doi: 10.1016/j.fsi.2012.11.044

PubMed Abstract | CrossRef Full Text | Google Scholar

Geay, F., Wenon, F. G. D., Tinti, J. M. E., and Tocher, S. N. M. M. D. R. (2015). Dietary linseed oil reduces growth while differentially impacting LC-PUFA synthesis and accretion into tissues in Eurasian Perch (Perca fluviatilis). Lipids 50, 1219–1232. doi: 10.1007/s11745-015-4079-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Gjøen, T., Kleveland, E. J., Moya-Falcón, C., Frøystad, M. K., Vegusdal, A., Hvattum, E., et al. (2007). Effects of dietary thia fatty acids on lipid composition, morphology and macrophage function of Atlantic salmon (Salmo salar L.) kidney. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 148, 103–111. doi: 10.1016/j.cbpb.2007.04.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Gjøen, T., Obach, A., Røsjø, C., Helland, B. G., Rosenlund, G., Hvattum, E., et al. (2004). Effect of dietary lipids on macrophage function, stress susceptibility and disease resistance in Atlantic salmon (Salmo salar). Fish Physiol. Biochem. 30, 149–161. doi: 10.1007/s10695-005-4318-7

CrossRef Full Text | Google Scholar

Glencross, B. D., De Santis, C., Bicskei, B., Taggart, J. B., Bron, J. E., Betancor, M. B., et al. (2015). A comparative analysis of the response of the hepatic transcriptome to dietary docosahexaenoic acid in Atlantic salmon (Salmo salar) post-smolts. BMC Genomics 16:684. doi: 10.1186/s12864-015-1810-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Glencross, B. D., Tocher, D. R., Matthew, C., and Bell, J. G. (2014). Interactions between dietary docosahexaenoic acid and other long-chain polyunsaturated fatty acids on performance and fatty acid retention in post-smolt Atlantic salmon (Salmo salar). Fish Physiol. Biochem. 40, 1213–1227. doi: 10.1007/s10695-014-9917-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Hastings, N., Agaba, M. K., Tocher, D. R., Zheng, X., Dickson, C. A., Dick, J. R., et al. (2004). Molecular cloning and functional characterization of fatty acyl desaturase and elongase cDNAs involved in the production of eicosapentaenoic and docosahexaenoic acids from α-linolenic acid in Atlantic salmon (Salmo salar). Mar. Biotechnol. 6, 463–474. doi: 10.1007/s10126-004-3002-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Hellemans, J., Mortier, G., De Paepe, A., Speleman, F., and Vandesompele, J. (2007). qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 8:R19. doi: 10.1186/gb-2007-8-2-r19

PubMed Abstract | CrossRef Full Text | Google Scholar

Hixson, S. M., Parrish, C. C., and Anderson, D. M. (2014). Changes in tissue lipid and fatty acid composition of farmed rainbow trout in response to dietary camelina oil as a replacement of fish oil. Lipids 49, 97–111. doi: 10.1007/s11745-013-3862-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Hixson, S. M., Parrish, C. C., Xue, X., Wells, J. S., Collins, S. A., Anderson, D. M., et al. (2017). Growth performance, tissue composition, and gene expression responses in Atlantic salmon (Salmo salar) fed varying levels of different lipid sources. Aquaculture 467, 76–88. doi: 10.1016/j.aquaculture.2016.04.011

CrossRef Full Text | Google Scholar

Holen, E., Araujo, P., Sissener, N. H., Rosenlund, G., and Waagbø, R. (2018). A comparative study: difference in omega-6/omega-3 balance and saturated fat in diets for Atlantic salmon (Salmo salar) affect immune-, fat metabolism-, oxidative and apoptotic-gene expression, and eicosanoid secretion in head kidney leukocytes. Fish Shellfish Immunol. 72, 57–68. doi: 10.1016/j.fsi.2017.10.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Holen, E., He, J., Espe, M., Chen, L., and Araujo, P. (2015). Combining eicosapentaenoic acid, decosahexaenoic acid and arachidonic acid, using a fully crossed design, affect gene expression and eicosanoid secretion in salmon head kidney cells invitro. Fish Shellfish Immunol. 45:695.e703. doi: 10.1016/j.fsi.2015.05.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, D. C., Ding, X., and Daynes, R. A. (2002). Nuclear receptor peroxisome proliferator-activated receptor α (PPARα) is expressed in resting murine lymphocytes. J. Biol. Chem. 277, 6838–6845. doi: 10.1074/jbc.M106908200

PubMed Abstract | CrossRef Full Text | Google Scholar

Jordal, A.-E. O., Torstensen, B. E., Tsoi, S., Tocher, D. R., Lall, S. P., and Douglas, S. E. (2005). Dietary rapeseed oil affects the expression of genes involved in hepatic lipid metabolism in Atlantic salmon (Salmo salar L.). J. Nutr. 135, 2355–2361. doi: 10.1093/jn/135.10.2355

PubMed Abstract | CrossRef Full Text | Google Scholar

Jump, D. B., Botolin, D., Wang, Y., Xu, J., Christian, B., and Demeure, O. (2005). Fatty acid regulation of hepatic gene transcription. J. Nutr. 135, 2503–2506. doi: 10.1093/jn/135.11.2503

PubMed Abstract | CrossRef Full Text | Google Scholar

Katan, T., Caballero-Solares, A., Taylor, R. G., Rise, M. L., and Parrish, C. C. (2019). Effect of plant-based diets with varying ratios of ω6 to ω3 fatty acids on growth performance, tissue composition, fatty acid biosynthesis and lipid-related gene expression in Atlantic salmon (Salmo salar). Comp. Biochem. Physiol. Part D Genomics Proteomics 30, 290–304. doi: 10.1016/j.cbd.2019.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Kersten, S. (2014). Integrated physiology and systems biology of PPARα. Mol. Metab. 3, 354–371. doi: 10.1016/j.molmet.2014.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, J. Y., Sohn, K. H., Rhee, S. H., and Hwang, D. (2001). Saturated fatty acids, but not unsaturated fatty acids, induce the expression of cyclooxygenase-2 mediated through Toll-like receptor 4. J. Biol. Chem. 276, 16683–16689. doi: 10.1074/jbc.M011695200

CrossRef Full Text | Google Scholar

Martinez-Rubio, L., Morais, S., Evensen, Ø., Wadsworth, S., Vecino, J. G., Ruohonen, K., et al. (2013). Effect of functional feeds on fatty acid and eicosanoid metabolism in liver and head kidney of Atlantic salmon (Salmo salar L.) with experimentally induced heart and skeletal muscle inflammation. Fish Shellfish Immunol. 34, 1533–1545. doi: 10.1016/j.fsi.2013.03.363

PubMed Abstract | CrossRef Full Text | Google Scholar

McNamara, R. K., Able, J., Jandacek, R., Rider, T., and Tso, P. (2009). Inbred C57BL/6J and DBA/2J mouse strains exhibit constitutive differences in regional brain fatty acid composition. Lipids 44, 1–8. doi: 10.1007/s11745-008-3244-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Monroig, Ó., Zheng, X., Morais, S., Leaver, M. J., Taggart, J. B., and Tocher, D. R. (2010). Multiple genes for functionalΔ 6 fatty acyl desaturases (Fad) in Atlantic salmon (Salmo salar L.): gene and cDNA characterization, functional expression, tissue distribution and nutritional regulation. Biochim. Biophys. Acta Mol. Cell Biol. Lipids 1801, 1072–1081. doi: 10.1016/j.bbalip.2010.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Montero, D., Kalinowski, T., Obach, A., Robaina, L., Tort, L., Caballero, M. J., et al. (2003). Vegetable lipid sources for gilthead seabream (Sparus aurata): effects on fish health. Aquaculture 225, 353–370. doi: 10.1016/S0044-8486(03)00301-6

CrossRef Full Text | Google Scholar

Montero, D., Torrecillas, S., Benítez-Dorta, V., Caballero, M. J., Izquierdo, M. S., and Zamorano, M. J. (2019). Effects of dietary vegetable oils on the expression of eicosanoid receptors genes in Senegalese sole (Solea senegalensis) intestine. Aquac. Rep. 15:100201. doi: 10.1016/j.aqrep.2019.100201

PubMed Abstract | CrossRef Full Text | Google Scholar

Morais, S., Monroig, O., Zheng, X., Leaver, M. J., and Tocher, D. R. (2009). Highly unsaturated fatty acid synthesis in Atlantic salmon: characterization of ELOVL5-and ELOVL2-like elongases. Mar. Biotechnol. 11, 627–639. doi: 10.1007/s10126-009-9179-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Morais, S., Pratoomyot, J., Taggart, J. B., Bron, J. E., Guy, D. R., Bell, J. G., et al. (2011). Genotype-specific responses in Atlantic salmon (Salmo salar) subject to dietary fish oil replacement by vegetable oil: a liver transcriptomic analysis. BMC Genomics 12:255. doi: 10.1186/1471-2164-12-255

PubMed Abstract | CrossRef Full Text | Google Scholar

National Research Council (NRC) (2011). Nutrient Requirements of Fish and Shrimp. Washington, DC: The National academies press.

Google Scholar

Neuman, J. C., Schaid, M. D., Brill, A. L., Fenske, R. J., Kibbe, C. R., Fontaine, D. A., et al. (2017). Enriching islet phospholipids with eicosapentaenoic acid reduces prostaglandin E2 signaling and enhances diabetic β-cell function. Diabetes 66, 1572–1585. doi: 10.2337/db16-1362

PubMed Abstract | CrossRef Full Text | Google Scholar

Olsen, R. E., Svardal, A., Eide, T., and Wargelius, A. (2012). Stress and expression of cyclooxygenases (cox1, cox2a, cox2b) and intestinal eicosanoids, in Atlantic salmon, Salmo salar L. Fish Physiol. Biochem. 38, 951–962. doi: 10.1007/s10695-011-9581-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Olsvik, P. A., Lie, K. K., Jordal, A.-E. O., Nilsen, T. O., and Hordvik, I. (2005). Evaluation of potential reference genes in real-time RT-PCR studies of Atlantic salmon. BMC Mol. Biol. 6:21. doi: 10.1186/1471-2199-6-21

PubMed Abstract | CrossRef Full Text | Google Scholar

Orsavova, J., Misurcova, L., Ambrozova, J. V., Vicha, R., and Mlcek, J. (2015). Fatty acids composition of vegetable oils and its contribution to dietary energy intake and dependence of cardiovascular mortality on dietary intake of fatty acids. Int. J. Mol. Sci. 16, 12871–12890. doi: 10.3390/ijms160612871

PubMed Abstract | CrossRef Full Text | Google Scholar

Parrish, C. C. (1987). Separation of aquatic lipid classes by chromarod thin-layer chromatography with measurement by latroscan flame ionization detection. Can. J. Fish. Aquat. Sci. 44, 722–731. doi: 10.1139/f87-087

CrossRef Full Text | Google Scholar

Parrish, C. C. (1999). Determination of total lipid, lipid classes, and fatty acids in aquatic samples. Lipids Freshw. Ecosyst. 635, 4–20. doi: 10.1007/978-1-4612-0547-0_2

CrossRef Full Text | Google Scholar

Pfaffl, M. W. (2001). A new mathematical model for relative quantification in real-time RT–PCR. Nucleic Acids Res. 29:e45. doi: 10.1093/nar/29.9.e45

PubMed Abstract | CrossRef Full Text | Google Scholar

Pickova, J., and Mørkøre, T. (2007). Alternate oils in fish feeds. Eur. J. Lipid Sci. Technol. 109, 256–263. doi: 10.1002/ejlt.200600222

CrossRef Full Text | Google Scholar

Ponter, A. A., Guyader-Joly, C., Nuttinck, F., Grimard, B., and Humblot, P. (2012). Oocyte and embryo production and quality after OPU-IVF in dairy heifers given diets varying in their n-6/n-3 fatty acid ratio. Theriogenology 78, 632–645. doi: 10.1016/j.theriogenology.2012.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Rise, M. L., Hall, J. R., Nash, G. W., Xue, X., Booman, M., Katan, T., et al. (2015). Transcriptome profiling reveals that feeding wild zooplankton to larval Atlantic cod (Gadus morhua) influences suites of genes involved in oxidation-reduction, mitosis, and selenium homeostasis. BMC Genomics 16:1016. doi: 10.1186/s12864-015-2120-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Rise, M. L., Hall, J. R., Rise, M., Hori, T. S., Browne, M. J., Gamperl, A. K., et al. (2010). Impact of asymptomatic nodavirus carrier state and intraperitoneal viral mimic injection on brain transcript expression in Atlantic cod (Gadus morhua). Physiol. Genomics 42, 266–280. doi: 10.1152/physiolgenomics.00168.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Rocha, D. M., Caldas, A. P., Oliveira, L. L., Bressan, J., and Hermsdorff, H. H. (2016). Saturated fatty acids trigger TLR4-mediated inflammatory response. Atherosclerosis 244, 211–215. doi: 10.1016/j.atherosclerosis.2015.11.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruyter, B., Moya-Falcón, C., Rosenlund, G., and Vegusdal, A. (2006). Fat content and morphology of liver and intestine of Atlantic salmon (Salmo salar): effects of temperature and dietary soybean oil. Aquaculture 252, 441–452. doi: 10.1016/j.aquaculture.2005.07.014

CrossRef Full Text | Google Scholar

Salini, M. J., Wade, N. M., Araújo, B. C., Turchini, G. M., and Glencross, B. D. (2016). Eicosapentaenoic acid, arachidonic acid and eicosanoid metabolism in juvenile barramundi Lates calcarifer. Lipids 51, 973–988. doi: 10.1007/s11745-016-4167-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Sargent, J., Bell, G., McEvoy, L., Tocher, D., and Estevez, A. (1999). Recent developments in the essential fatty acid nutrition of fish. Aquaculture 177, 191–199. doi: 10.1016/S0044-8486(99)00083-6

CrossRef Full Text | Google Scholar

Sargent, J. R., Tocher, D. R., and Bell, J. G. (2002). The lipids. Fish Nutr. 3, 181–257. doi: 10.1016/B978-012319652-1/50005-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Sourabié, A., Mandiki, S. N. M., Geay, F., Ahoulé, A. G., Naert, N., Toguyeni, A., et al. (2019). Tropical vegetable oils do not alter growth performance in African catfish through a high n-6 polyunsaturated fatty acids biosynthesis capacity. Lipids 54, 329–345. doi: 10.1002/lipd.12145

CrossRef Full Text | Google Scholar

Sprague, M., Dick, J. R., and Tocher, D. R. (2016). Impact of sustainable feeds on omega-3 long-chain fatty acid levels in farmed Atlantic salmon, 2006–2015. Sci. Rep. 6, 1–9. doi: 10.1038/srep21892

PubMed Abstract | CrossRef Full Text | Google Scholar

Tacon, A. G. J., and Metian, M. (2015). Feed matters: satisfying the feed demand of aquaculture. Rev. Fish. Sci. Aquac. 23, 1–10. doi: 10.1080/23308249.2014.987209

CrossRef Full Text | Google Scholar

Tapiero, H., Nguyen Ba, G., Couvreur, P., and Tew, K. D. (2002). Polyunsaturated fatty acids (PUFA) and eicosanoids in human health and pathologies. Biomed. Pharmacother. 56, 215–222. doi: 10.1016/S0753-3322(02)00193-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, K. D., Henderson, R. J., and Tatner, M. F. (1995). A comparison of the lipid composition of peripheral blood cells and head kidney leucocytes of Atlantic salmon (Salmo salar L.). Comp. Biochem. Physiol. Part B Biochem. 112, 83–92. doi: 10.1016/0305-0491(95)00051-9

CrossRef Full Text | Google Scholar

Tocher, D. R. (2010). Fatty acid requirements in ontogeny of marine and freshwater fish. Aquac. Res. 41, 717–732. doi: 10.1111/j.1365-2109.2008.02150.x

CrossRef Full Text | Google Scholar

Tocher, D. R. (2015). Omega-3 long-chain polyunsaturated fatty acids and aquaculture in perspective. Aquaculture 449, 94–107. doi: 10.1016/j.aquaculture.2015.01.010

CrossRef Full Text | Google Scholar

Tocher, D. R., Bell, J. G., Dick, J. R., Henderson, R. J., McGhee, F., Michell, D., et al. (2000). Polyunsaturated fatty acid metabolism in Atlantic salmon (Salmo salar) undergoing parr-smolt transformation and the effects of dietary linseed and rapeseed oils. Fish Physiol. Biochem. 23, 59–73. doi: 10.1023/A:1007807201093

PubMed Abstract | CrossRef Full Text | Google Scholar

Torstensen, B. E., and Tocher, D. R. (2010). “The Effects of fish oil replacement on lipid metabolism of fish,” in Fish Oil Replacement and Alternative Lipid Sources in Aquaculture Feeds, eds G. M. Turchini, W. K. Ng and D. R. Tocher (Boca Raton, FL: CRC Press), 405–437. doi: 10.1201/9781439808634-c13

CrossRef Full Text | Google Scholar

Tort, L., Balasch, J. C., and Mackenzie, S. (2003). Fish immune system. A crossroads between innate and adaptive responses. Inmunología 22, 277–286.

Google Scholar

Turchini, G. M., Torstensen, B. E., and Ng, W. (2009). Fish oil replacement in finfish nutrition. Rev. Aquac. 1, 10–57. doi: 10.1111/j.1753-5131.2008.01001.x

CrossRef Full Text | Google Scholar

Vagner, M., and Santigosa, E. (2011). Characterization and modulation of gene expression and enzymatic activity of delta-6 desaturase in teleosts: a review. Aquaculture 315, 131–143. doi: 10.1016/j.aquaculture.2010.11.031

CrossRef Full Text | Google Scholar

Vandesompele, J., De Preter, K., Pattyn, F., Poppe, B., Van Roy, N., De Paepe, A., et al. (2002). Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 3:RESEARCH0034. doi: 10.1186/gb-2002-3-7-research0034

PubMed Abstract | CrossRef Full Text | Google Scholar

Varga, T., Czimmerer, Z., and Nagy, L. (2011). PPARs are a unique set of fatty acid regulated transcription factors controlling both lipid metabolism and inflammation. Biochim. Biophys. Acta Mol. Basis Dis. 1812, 1007–1022. doi: 10.1016/j.bbadis.2011.02.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Weaver, K. L., Ivester, P., Chilton, J. A., Wilson, M. D., Pandey, P., and Chilton, F. H. (2008). The content of favorable and unfavorable polyunsaturated fatty acids found in commonly eaten fish. J. Am. Diet. Assoc. 108, 1178–1185. doi: 10.1016/j.jada.2008.04.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, H., Zhou, Y., Jiang, S., Huang, F., Peng, J., and Jiang, S. (2016). Transcriptional response of porcine skeletal muscle to feeding a linseed-enriched diet to growing pigs. J. Anim. Sci. Biotechnol. 7:6. doi: 10.1186/s40104-016-0064-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, G., Du, F., Wang, Y., Li, Y., Nie, Z., and Xu, P. (2019). Molecular characterization and differential expression of two prostaglandin E synthase 2 orthologs in coilia nasus. Pak. J. Zool. 51, 483–494. doi: 10.17582/journal.pjz/2019.51.2.483.494

CrossRef Full Text | Google Scholar

Xue, X., Hixson, S. M., Hori, T. S., Booman, M., Parrish, C. C., Anderson, D. M., et al. (2015). Atlantic salmon (Salmo salar) liver transcriptome response to diets containing Camelina sativa products. Comp. Biochem. Physiol. Part D Genomics Proteomics 14, 1–15. doi: 10.1016/j.cbd.2015.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Young, K. (2009). Omega-6 (n-6) and omega-3 (n-3) fatty acids in tilapia and human health: a review. Int. J. Food Sci. Nutr. 60, 203–211. doi: 10.1080/09637480903140503

PubMed Abstract | CrossRef Full Text | Google Scholar

Ytrestøyl, T., Aas, T. S., and Åsgård, T. (2015). Utilisation of feed resources in production of Atlantic salmon (Salmo salar) in Norway. Aquaculture 448, 365–374. doi: 10.1016/j.aquaculture.2015.06.023

CrossRef Full Text | Google Scholar

Zapata, A., Diez, B., Cejalvo, T., Gutierrez-de Frias, C., and Cortes, A. (2006). Ontogeny of the immune system of fish. Fish Shellfish Immunol. 20, 126–136. doi: 10.1016/j.fsi.2004.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, N., Monroig, Ó., Navarro, J. C., Xiang, X., Li, Y., Du, J., et al. (2019). Molecular cloning, functional characterization and nutritional regulation of two elovl4b elongases from rainbow trout (Oncorhynchus mykiss). Aquaculture 511:734221. doi: 10.1016/j.aquaculture.2019.734221

CrossRef Full Text | Google Scholar

Zheng, X., Tocher, D. R., Dickson, C. A., Bell, J. G., and Teale, A. J. (2004). Effects of diets containing vegetable oil on expression of genes involved in highly unsaturated fatty acid biosynthesis in liver of Atlantic salmon (Salmo salar). Aquaculture 236, 467–483. doi: 10.1016/j.aquaculture.2004.02.003

CrossRef Full Text | Google Scholar

Zheng, X., Torstensen, B. E., Tocher, D. R., Dick, J. R., Henderson, R. J., and Bell, J. G. (2005). Environmental and dietary influences on highly unsaturated fatty acid biosynthesis and expression of fatty acyl desaturase and elongase genes in liver of Atlantic salmon (Salmo salar). Biochim. Biophys. Acta Mol. Cell Biol. Lipids 1734, 13–24. doi: 10.1007/s11745-005-1355-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, K. C., Song, L., Guo, H. Y., Guo, L., Zhang, N., Liu, B. S., et al. (2019). Elovl4a participates in LC-PUFA biosynthesis and is regulated by PPARαβ in golden pompano Trachinotus ovatus (Linnaeus 1758). Sci. Rep. 9, 1–11. doi: 10.1038/s41598-019-41288-w

PubMed Abstract | CrossRef Full Text | Google Scholar