Skip to main content

Genetic architecture and genomic selection of fatty acid composition predicted by Raman spectroscopy in rainbow trout

Abstract

Background

In response to major challenges regarding the supply and sustainability of marine ingredients in aquafeeds, the aquaculture industry has made a large-scale shift toward plant-based substitutions for fish oil and fish meal. But, this also led to lower levels of healthful n−3 long-chain polyunsaturated fatty acids (PUFAs)—especially eicosapentaenoic (EPA) and docosahexaenoic (DHA) acids—in flesh. One potential solution is to select fish with better abilities to retain or synthesise PUFAs, to increase the efficiency of aquaculture and promote the production of healthier fish products. To this end, we aimed i) to estimate the genetic variability in fatty acid (FA) composition in visceral fat quantified by Raman spectroscopy, with respect to both individual FAs and groups under a feeding regime with limited n-3 PUFAs; ii) to study the genetic and phenotypic correlations between FAs and processing yields- and fat-related traits; iii) to detect QTLs associated with FA composition and identify candidate genes; and iv) to assess the efficiency of genomic selection compared to pedigree-based BLUP selection.

Results

Proportions of the various FAs in fish were indirectly estimated using Raman scattering spectroscopy. Fish were genotyped using the 57 K SNP Axiom™ Trout Genotyping Array. Following quality control, the final analysis contained 29,652 SNPs from 1382 fish. Heritability estimates for traits ranged from 0.03 ± 0.03 (n-3 PUFAs) to 0.24 ± 0.05 (n-6 PUFAs), confirming the potential for genomic selection. n-3 PUFAs are positively correlated to a decrease in fat deposition in the fillet and in the viscera but negatively correlated to body weight. This highlights the potential interest to combine selection on FA and against fat deposition to improve nutritional merit of aquaculture products. Several QTLs were identified for FA composition, containing multiple candidate genes with indirect links to FA metabolism. In particular, one region on Omy1 was associated with n-6 PUFAs, monounsaturated FAs, linoleic acid, and EPA, while a region on Omy7 had effects on n-6 PUFAs, EPA, and linoleic acid. When we compared the effectiveness of breeding programmes based on genomic selection (using a reference population of 1000 individuals related to selection candidates) or on pedigree-based selection, we found that the former yielded increases in selection accuracy of 12 to 120% depending on the FA trait.

Conclusion

This study reveals the polygenic genetic architecture for FA composition in rainbow trout and confirms that genomic selection has potential to improve EPA and DHA proportions in aquaculture species.

Peer Review reports

Background

Over the last several decades, aquatic food production has evolved away from capture fisheries toward the culture of increasing numbers of farmed fish species [1]. However, aquaculture operations face major challenges regarding the sustainability of the feed used, particularly with regard to marine ingredients. To address these concerns, fish oil and fish meal—the traditional sources of proteins and lipids in aquafeeds—have been largely substituted (60–80%) with plant-based ingredients. However, compared to fish oils, vegetable oils differ in their fatty acid composition: they are rich in oleic acid (OA; C18:1 n-9), linoleic acid (LA; C18:2 n-6), and alpha-linolenic acid (ALA; C18:3 n-3) and contain very low proportion or no n-3 long-chain polyunsaturated fatty acids (n-3 LC PUFAs). In farmed trout, one result of this dietary shift has been reductions in the levels of n-3 LC PUFAs, especially eicosapentaenoic (EPA; C20:5 n-3) and docosahexaenoic (DHA; C22:6 n-3) acids, in flesh [2]. n-3 LC PUFAs are known to have beneficial effects on human health, including the prevention of a range of cardiovascular and inflammatory diseases and neurological disorders [3, 4]. Freshwater fish are theoretically capable of biosynthesising DHA and EPA via desaturation and elongation of the ALA found in some vegetable oils [5]. In practice, however, such bioconversion is insufficient to compensate for the lack of dietary n-3 LC PUFAs, resulting in a significant reduction in levels of these healthful fatty acids (FAs) in fish tissues [5, 6]. Especially, in rainbow trout, two fads2 genes encoding proteins with delta 5 and delta 6 desaturase activities and two elongase enzymes Elovl5 and Elovl 2 have been isolated and functionally characterised [7, 8]. Rainbow trout is dependent on Elovl2 for 22:5 n-3 to 24:5 n-3 synthesis and ultimately DHA synthesis [8]. When diets are high in alpha-linolenic acid (ALA, 18: 3 n-3) with no added EPA or DHA, fads2, Elovl5 and Elovl2 are most highly expressed in rainbow trout liver [9]. PUFAs are also important in the fish life cycle, most notably for their roles in reproduction, egg quality, and offspring development [5], and have effects on nutritional quality of fish flesh for human consumption [6, 10]. The relative amounts of EPA and DHA formed are determined by the activities of desaturase and elongase enzymes, which are themselves influenced by several factors, such as nutrition, the environment, physiology, and genetics [11]. To date, numerous studies have examined the characteristics of FA metabolism and the effect of dietary oil sources on the FA composition of farmed fish in attempts to find solutions for this problem [12,13,14,15,16,17,18].

To meet consumer demand, producers are constantly making improvements in husbandry techniques, nutrition, and genetic management. One potential strategy for meeting the demand for farmed fish without compromising nutritional value could be to combine genetic selection with changes in commercial feed formulations. Recent investigations of the genetic variability underlying lipid deposition and metabolism in fish have identified a highly heritable genetic component that governs the capacity to synthesise and/or deposit LC-PUFAs [19,20,21,22]. In trout, studies have revealed the efficiency of divergent selection for total fat content [17, 23,24,25] and potential interactions between dietary lipid level and genetic selection for body fat [26, 27]. Thus, to counteract the diet-based decline in PUFA content in flesh, one solution may be breeding programmes that select fish with better abilities to retain and/or synthesise PUFAs. To date, studies have illustrated the potential of selective breeding to increase n-3 LC PUFA levels in salmon [20, 28,29,30], yellow croaker [22], tilapia [31], Asian sea bass [32] and common carp [33]. Similar work has also been conducted in trout [30], but as yet, we lack knowledge on the genetic architecture that shapes the relative proportions of individual FAs in this species. A common strategy for this purpose is the genome-wide association study (GWAS), which has been used to identify the genetic regions and loci significantly associated with FA composition in species such as cattle (in meat [34,35,36] and milk [37]) and pigs [38,39,40,41], common carp [42], Asian seabass [32, 43], tilapia [44] and Atlantic salmon [28]. In rainbow trout, though, the genetic parameters of FA composition remain unknown, and no GWAS has been performed to increase our knowledge of the genetics of n-3 LC PUFA composition. Furthermore, the relationships between FA composition and traits related to lipid deposition, weight, yield, or quality have not yet been characterised. Phenotype-based research has shown that total fat content (as measured with a Fatmeter) in the muscle increases with the growth and development of fish [45, 46], but the genetic correlations between production traits and FA composition—which would play a crucial role in optimising the efficiency of breeding programmes—remain unknown.

In the literature, the technology most commonly used for FA characterisation is chemical extraction and gas chromatography. However, such analyses are expensive, invasive, and time-consuming, and are thus not easily applicable to a breeding programme. In this context, what is needed are alternative methods that are affordable and potentially non-invasive, which could be used to estimate the proportions of different FAs (group or individual FA) in a population of fish that is large enough to enable effective analyses of the genetic architecture of traits. One potential approach could be the use of Raman scattering spectroscopy, a rapid, non-destructive method for molecular characterisation based on vibrational spectrometry. Specifically, this method enables the qualitative and quantitative characterisation of molecules through analyses of spectral bands with respect to the fundamental vibrational modes of their chemical bonds. This technology has already been used to determine total lipid concentration in minced salmon flesh [47] and FA composition in pork adipose tissue [48], as well as in fish oils [49, 50] and in Atlantic salmon flesh [51].

As a step toward the long-term objective of increasing the nutritional value of farmed rainbow trout and ensuring high levels of n-3 PUFAs in the flesh, the main goal of the present study was to investigate the genetics of n-3 LC PUFA composition. Specifically, we used a commercially selected population of rainbow trout to (1) estimate the genetic variability of the FA composition of visceral adipose tissue, as indirectly estimated by Raman spectroscopy; (2) analyse genetic and phenotypic correlations among different FA traits and between FAs and production traits linked to yields and fat deposition, (3) detect QTLs associated with FA proportions and identify candidate genes present within those regions, and (4) estimate the efficiency of genomic selection (GS) compared to pedigree-based BLUP selection using phenotypes predicted by Raman spectroscopy. In this study, we targeted quantification of FA proportions in adipocytes from visceral fat. A reason was that Raman spectroscopy performed on adipocytes is less tedious and expensive than the successive mincing and lyophilisation of flesh reported in previous studies. A second one, was that the procedure proposed could also be usable as highthroughput phenotyping technology by biopsy in a non-destructive manner on live candidates as Raman spectroscopy required a limited size of sample (< 1 g). The methodology of the study encompassing gaz chromatography, Raman spectroscopy, MRI, microwaves and genotyping for calibration step and commercially selected fish was illustrated in Fig. 1.

Fig. 1
figure 1

Description of the methodology of the study encompassing gaz chromatography, Raman spectroscopy, phenotyping and genotyping for calibration and commercial selected fish

Results

Basic characteristics of fatty acid composition

Descriptive statistics of the fish fatty acids (FAs) proportions predicted by Raman spectroscopy in this study are presented in Table 1 and Fig. 2.

Table 1 Summary statistics for fatty acid traits in rainbow trout
Fig. 2
figure 2

Bar plot showing fatty acid prediction in percentage for SFA, MUFA (OA and others), PUFA (n-3 PUFA (ALA; EPA; DHA); n-6 PUFA (LA; ARA) and others) in rainbow trout

In this study, the most abundant individual fatty acid was oleic acid (C18:1, 43.76%), which represented about 90% of all monounsaturated fatty acids (MUFAs) present in visceral fat. Saturated fatty acids (SFAs) made up 22.69% of total FA content, while polyunsaturated fatty acids (PUFAs) represented 27.05%. This latter group was composed of approximately 33% n-3 PUFAs (sum of omega-3 fatty acids) and 64% n-6 PUFAs (sum of omega-6 fatty acids). Of the individual n-6 PUFAs, the most abundant was linoleic acid (LA, 16.08%), which represented ca. 60% of all PUFAs and about 93% of n-6 PUFAs. Of the individual n-3 PUFAs, alpha-linolenic acid (ALA) had the highest concentration, accounting for 50% of n-3 PUFAs and about 17% of total PUFA content. DHA was the second most abundant n-3 PUFA, representing ca. 22% of n-3 PUFAs and 7% of total PUFAs. Taken together, the sum of EPA + DHA corresponded to 35% of n-3 PUFAs and 11% of PUFAs as a whole. Overall, the FA composition of fish was similar to that of the feed (Additional file 1).

Among individuals, we observed relatively little variability in the estimated proportions of different types of FAs (coefficient of variation (CV) ranging from 2 to 5% for SFAs, MUFAs, and PUFAs). Instead, much more variation was detected among individuals in the estimated proportions of individual FAs (CV ranging from 3 to 31%), with the largest CV found for DHA. In general, proportions of n-3 PUFAs were more variable among fish than those of n-6 PUFAs (CVs of 13% and ca. 3%, respectively).

Heritability estimates and correlations among FAs

Figure 3 depicts heritability (h2) estimates based on pedigree information for individual FAs and FA groups, as well as phenotypic correlations (rp) among the 12 traits.

Fig. 3
figure 3

Genetic parameters of fatty acid traits in rainbow trout. Heritabilities (± standard error) in bold red on the diagonal, genetic correlations (± standard error) in the upper triangle, phenotypic correlations in the lower triangle

In general, heritability estimates were low and varied from 0.02 ± 0.03 to 0.24 ± 0.05. With respect to FA groups, higher heritability was estimated for PUFAs (h2 = 0.16 ± 0.05) compared to SFAs (h2 = 0.08 ± 0.04) and MUFAs (h2 = 0.12 ± 0.04). Of the two types of PUFAs, n-6 PUFAs had the higher heritability (h2 = 0.24 ± 0.05). For individual fatty acids, the highest heritabilities were estimated for arachidonic acid (ARA), LA, and EPA (0.21 ± 0.05, 0.18 ± 0.05, and 0.16 ± 0.05, respectively).

Our analysis of phenotypic correlations among traits revealed that, as a group, values for PUFAs as a whole were more strongly correlated to those of n-3 PUFAs than of n-6 PUFAs. Concentrations of ARA, an individual n-6 PUFA, were highly correlated to those of the larger PUFA group, n-3 PUFAs, and the individual FAs EPA, DHA, and sum EPA + DHA. The two types of n-3 LC PUFAs, EPA and DHA, were highly correlated with each other. EPA and sum EPA + DHA were more strongly correlated to n-6 PUFAs, LA, and ARA (rp > 0.5) than to n-3 PUFAs and ALA (rp < 0.5).

When we examined the genetic correlations (rg) among FAs, the results were similar to the phenotypic correlations. From a group perspective, SFAs were highly positively genetically correlated to MUFAs (rg = 1 ± 0) and highly negatively genetically correlated to PUFAs (− 0.99 ± 0.22). EPA was highly genetically correlated to LA (0.89 ± 0.08). DHA and the sum of EPA + DHA were positively correlated to LA (rg = 1 ± 0.13, 0.80 ± 0.66, respectively), while the sum of EPA + DHA was also highly correlated to ARA (rg = 1 ± 0.53). However, for the majority of FAs the standard errors of the genetic correlations were high (up to ±2.86 for ALA) and some parameters did not converge, thus some results would need confirmation.

Correlations between fatty acid proportions and other traits

Genetic correlations between production traits and fatty acids are presented in Fig. 4 and phenotypic correlations can be found in Additional file 2.

Fig. 4
figure 4

Estimates of genetic correlations between proportions of fatty acids and production or quality traits

The growth parameters BW and K and fat content parameters Fat, MRI_F%, MRI_F_sc%, and MRI_F_F% were moderately and positively genetically correlated with SFAs and MUFAs (0.2 ± 0.32 to 0.64 ± 0.36). Instead, generally speaking those same production traits tended to be negatively correlated with all PUFA groups and individual PUFAs (− 0.11 ± 0.6 to − 0.93 ± 1.25). BW was highly positively correlated with SFAs and MUFAs (0.64 ± 0.36 and 0.49 ± 0.22, respectively) and highly negatively correlated with PUFAs, n-3 PUFAs, and n-6 PUFAs (− 0.6 ± 0.15, − 0.78 ± 0.73, and − 0.48 ± 0.14, respectively). Fat percentage correlated positively with SFAs and MUFAs (0.55 ± 0.25 and 0.25 ± 0.21), and negatively with PUFAs as a group (− 0.44 ± 0.16) as well as all individual PUFAs (− 0.11 ± 0.6 to − 0.74 ± 0.91).

PUFAs, n-3 PUFAs, and n-6 PUFAs were positively genetically correlated with yield traits (HGCarc% and Carc%), with the highest positive correlation found between Carc% and n-6 PUFAs (0.47 ± 0.13). These yield traits tended to be negatively correlated with SFAs and MUFAs, with the highest negative correlation found between Carc% and SFAs (− 0.44 ± 0.36).

DHA and the sum of EPA + DHA showed a moderate, positive correlation with HGCarc% and Carc% (0.17 ± 0.48 to 0.25 ± 0.57).

We must reiterate here that the standard errors of the genetic correlations were generally quite high; estimates of genetic correlations were thus relatively imprecise.

QTL detection and identification of candidate genes

GWAS analyses were performed for the 12 individual FAs or FA groups using approximately 30,000 SNPs. We identified 10 QTLs with evidence effects on the fatty acids proportions as well as 13 putative QTLs with only putative evidence (Fig. 5 and Table 2). Table 2 presents all regions detected for each FA trait and the percentage of total additive genetic variation explained by each. We will first describe the regions for which we detected strong evidence of an effect, then the regions associated with more than one FA, and finally the regions that explained more than 1% of the total additive variance of a trait. Regions associated with only one FA and that accounted for less than 1% of genetic variance are included in Table 2 but will not be described here. The candidate genes identified within the different QTL regions are detailed in Additional file 3.

Fig. 5
figure 5

Manhattan plot of QTLs detected using a Bayesian GWAS for fatty acid traits. SNPs are plotted on the x-axis according to their position on each chromosome. The y-axis represents twice the natural logarithm of the Bayes factor associated with each trait. Red and blue dashed horizontal lines correspond to the threshold value of logBF ≥6.0 for definitive evidence of a QTL and the threshold value of logBF ≥5.0 for a putative QTL, respectively

Table 2 Summary statistics for GWAS for fatty acid traits in rainbow trout using BayesCπ method

Within a single region on Omy7, we identified with evidence effect a QTL for n-6 PUFAs and putative QTLs for LA and EPA. This region explained the highest percentage of genetic variance for any of the traits investigated—1.67, 0.86, and 1.02%, respectively—and also had the largest credibility interval, between 2653 kb and 3705 kb wide. The peak SNPs for n-6 PUFAs, LA, and EPA were located, respectively, at 8.34 Mb (coiled coil domain-containing protein 122-like), 8.76 Mb (intergenic region between tnika (traf2 and NCK-interacting protein kinase-like) and slc7a14b (probable cationic amino acid transporter) genes), and 9.88 Mb (intergenic region between ino1a (inositol-3-phosphate synthase 1-A-like) and gnao1a (guanine nucleotide-binding protein subunit alpha-11-like) genes) respectively.

Between 63.06 Mb and 64.29 Mb on Omy 1, we detected QTLs for n-6 PUFAs and LA, and putative QTLs for MUFAs and EPA. Depending on the trait, this region explained between 0.59 and 1.08% of the genetic variance and had a credibility interval that was between 688 kb and 1232 kb wide. The peak SNP for n-6 PUFAs and LA, at 64.29 Mb, was located in the intergenic region between tlx1 (T-cell leukaemia homeobox protein 1 isoform X1) and lbx1 (transcription factor LBX1-like). The peak SNP for MUFAs and EPA, at 63.63 Mb, was located in the PH and SEC7 domain-containing protein 1-like isoform region.

A region on Omy3 (located between 75.82 Mb and 75.96 Mb) contained QTLs for PUFAs and ARA, but these explained less than 0.5% of the genetic variance in each trait. The peak SNP for these QTLs was located in the intergenic region between tmtops2a (vertebrate ancient opsin-like) and jam2b (junctional adhesion molecule B-like isoform X1).

A region on Omy14 (located between 2.57 and 2.73 Mb) contained a QTL for n-6 PUFAs and putative QTLs for LA and EPA. This region explained a low percentage of the genetic variance in these traits—0.42, 0.45, and 0.59%, respectively—and had credibility intervals that spanned from 50 to 165 kb. The peak SNP for this region was located at 2.73 Mb, in the intergenic region between rims2b (regulating synaptic membrane exocytosis protein 2-like) and rims2a (regulating synaptic membrane exocytosis protein 2-like isoform X1).

On Omy26, one region was identified as a putative QTL for both n-6 PUFAs and LA. This region explained 0.56 and 0.19% of genetic variance, respectively, and was located between 24.05 Mb and 25.30 Mb. The peak SNP for this QTL was located at 24.05 Mb, in the intergenic region between rbtn1 (rhombotin-1) and rergl (ras-related and estrogen-regulated growth inhibitor-like protein).

Two other putative QTLs were detected that each explained more than 1% of the genetic variance associated with PUFAs as a group. The first was located on Omy7, had a credibility interval of 1239 kb, and explained 1.29% of the genetic variance for this trait. The peak SNP for this QTL was located at 63.78 Mb, in the znf1007 (zinc finger protein 595-like) gene. The second putative QTL was located on Omy22 between 4.94 and 8.31 Mb, and explained 1.6% of genetic variance. The peak SNP for this QTL was located in the lmo7 (LIM domain only protein 7) gene at 6.74 Mb.

No significant QTL was found for SFAs, n-3 PUFAs, OA, ALA, DHA, DPA, and EPA + DHA, which suggests that the underlying genetic architecture for these traits is very polygenic and affected by multiple loci with small effects in this population of rainbow trout.

Genomic selection

The accuracies of estimated breeding values (EBVs) and genomic estimated breeding values (GEBVs) are shown in Table 3 for all FAs; data on the efficiency of genomic selection, along with the inflation coefficients, can be found in Table 3.

Table 3 Mean and standard deviation (in brackets) over 40 replicates of the selection accuracy (Acc) and inflation coefficient (b) of EBVs and GEBVs for fatty acid traits

Regardless of the trait in question, all inflation coefficients were statistically indistinguishable from 1 (see Table 3) and nearly identical between pedigree-based (BLUP) and genomic (GBLUP) approaches.

GEBVs were more accurate than EBVs for all of the predicted traits. The accuracy of genomic selection was approximately 45% higher than that of pedigree-based selection for all FAs. The mean accuracy of GEBVs varied from 0.34 for ALA to 0.70 for n-6 PUFAs, whereas the corresponding accuracies of EBVs were 0.30 and 0.51, respectively. The highest increase in accuracy for a GEBV with respect to the corresponding EBV was obtained for DHA (+ 119.8%), while the lowest was obtained for ARA (+ 11.8%).

Discussion

In this study, we estimated the genetic parameters of FA composition predicted by Raman spectroscopy and the correlations between FA traits of visceral adipocytes and those associated with production and quality traits (body weight, carcass and filet yields, and lipid deposition) in a commercial population of rainbow trout. We identified quantitative trait loci and evaluated the potential for genomic selection aimed at increasing the abundance of long-chain omega-3 polyunsaturated fatty acids in these fish.

Raman phenotyping and sampling of visceral adipose tissue

In the majority of previous studies on aquatic species, genetic variability in FA composition has been evaluated using chemical analyses and gas chromatography [19, 20, 31, 33, 52, 53]. Due to the high cost of these analyses, such studies have typically relied on relatively small datasets, which limits the usefulness of this approach for rigorous genetic evaluation and selection. The inclusion of FA traits in breeding programmes would thus depend greatly on the availability of alternative methods for FA analysis. The development of breeding programs could also benefit from alternative non-destructive technologies to predict FA on live candidates. For this reason, the present study investigated the use of Raman spectroscopy, which is less expensive than GC-based analyses and would therefore enable the construction of the larger datasets (thousands of individuals) necessary for study of genetic architecture. As reported previously, this technology requires also a small size of tissue (few μm3) that could be frozen. Thus it could be a non-lethal biopsy applied on live candidates, making it affordable in breeding programs.

The accuracy of this technique is determined by calibration equations that are constructed using GC-based FA quantification [Prado E, Eklouh-Molinier C, Enez F, Causeur D, Blay C., Dupont-Nivet M, Labbé L, Petit V, Moreac A, Taupier G, Haffray P, Bugeon J, Corraze G, Nazabal, V. Prediction of fatty acids composition in the rainbow trout Oncorhynchus mykiss by using Raman micro-spectroscopy. Submitted]. For the majority of FAs (MUFAs, PUFAs, n-6 PUFAs, OA, LA, ALA, EPA, DHA, and EPA + DHA), the quality of prediction was high (correlation coefficient (R2) ≥ 0.75, as estimated using ridge regression methods). Only three of the FA traits tested (SFAs, ARA, and n-3 PUFAs) demonstrated weaker predictive power, with R2 < 0.75 [50]. One potential disadvantage of this technique, however, could be its limit of detection, which is around 1%; some of the values for individual FAs in the current study were close to this percentage. Curiously, though, the predictive power for the SFA group was quite limited (R2 = 0.42), even though the proportions calculated for this group were well above the detection threshold (around 22%). Further investigations are currently underway to study these results more thoroughly.

Another potential point of contention could be the use of adipocytes from visceral fat instead of minced fillet or adipocyte cells from myosepta. Major reasons of this choice were detailed in the introduction. Adipocytes are located in various adipose tissues (mostly visceral, subcutaneous, liver, red and white muscles, brain, pancreas, mandible, cranium, and tail fin); they play an important role in the long-term storage of FAs as triglycerides (lipogenesis) and their subsequent release (lipolysis) in the blood for use in cell growth, endocrine regulation, reproduction, or as energy sources after β-oxidation [54, 55]. The processes of adipocyte multiplication (hyperplasia) or growth in size (hypertrophy), as well as lipolysis, are regulated by a variety of enzymes, transcription factors, and hormones that also interact with the enzymes fatty acyl elongase (Elovl) and desaturase (Fad) to yield the final edible FA composition of the muscle. Genetic variability in the FA composition of the fillet results then from all the previous physiological process acting in interaction with important environmental factors such as swimming activity, water temperature, feed composition, or feeding practices. At this time, there is no published information about potential differences in FA composition between adipocytes in different locations in the fillet. In targeting visceral adipocyte composition, we aimed to evaluate genetic differences in lipogenesis- and lipolysis-associated processes while minimising, as much as possible, local interaction with other cell types or tissues. This may limit the applicability of our results to this type of tissue. The proportion of polar lipids is higher in muscle (around 25%) than in adipose tissue (5–10%). The n-3 long-chain PUFAs are mostly associated with sn-2 positions of phospholipids, so we can reasonably hypothesise that the relative proportion of n-3 LC-PUFA (in % of total FA) would be higher in muscle than in adipose tissue and then can be better predicted by Raman spectroscopy due to their higher abundance (in proportion). Our study provides a starting point—and a new technical methodology—for the characterisation of potential genetic variation of FA composition and will need further investigation to compare muscle and visceral FA composition.

Heritabilities for fatty acid composition and correlations

To our knowledge, this study is the first to report pedigree heritability estimates for all FA groups and individual FAs in rainbow trout. An earlier study in the same species estimated heritability for EPA (0.61 ± 0.17) and DHA (0.77 ± 0.11), but with a limited study population (220 fish from 44 families) [52]. Another study on Atlantic salmon reported a similarly high estimate for the n-3 PUFA group (0.77 ± 0.14) [19]. These estimates are very different than those we obtained for EPA (0.17 ± 0.05) and DHA (0.03 ± 0.03), probably due to differences in population structure, phenotyping method, and/or in the age and diet of the studied fish. In our case, the contrast between the limited proportions of PUFAs (15% of total FAs), EPA (1–2.3% of total FAs), and DHA (around 1% of total FAs) and the much-larger proportions of MUFAs (47–56% of total FAs) and SFAs (14–27% of total fat) may have hindered our ability to accurately estimate the underlying genetic variation. However, in general, our results are consistent with those documented from many other aquatic species. For example, in a study of Atlantic salmon whose diet contained fish meal and fish oil, and thus higher PUFA concentrations than in the present study, Horn et al. reported heritabilities that were similar to ours, in the range of 0.09 ± 0.06 to 0.26 ± 0.08 for individual n-3 PUFAs [20]. Although there were large differences between the two studies in the FA composition of the feed (EPA + DHA made up 3% of total dietary fat in the trout feed used here, compared with 17% in their salmon feed; trout feed was 25–30% fat versus 36% for salmon), this does not seem to have influenced the heritability estimates. For the shrimp Litopenaeus vannemei, published heritability estimates were even lower, ranging from 0 to 0.19 ± 0.07 for all fatty acids, 0.07 ± 0.05 for EPA, and 0.12 ± 0.06 for DHA [53]. In tilapia, heritability estimates ranged from 0 to 0.39 ± 0.11, with 0.10 ± 0.08 for EPA and 0.004 ± 0.07 for DHA [31]; in common carp, estimated heritabilities ranged from 0.03 ± 0.10 to 0.37 ± 0.22, with only three estimates different from zero: PUFAs (0.29 ± 0.17), n-3 PUFAs (0.37 ± 0.22), and EPA (0.34 ± 0.20) [33]. Thus, our estimates for the heritability of FA groups or individual FAs are in the same range as previous results obtained from other aquaculture species. Although these estimates were globally low, they do point to a genetic component in rainbow trout for FA composition and abundance, particularly for n-6 PUFAs, which could be exploited through genetic selection. The higher heritabilities we observed for n-6 PUFAs, both as a group and for individual FAs, were most likely due to the improved prediction accuracy for these FAs compared to SFAs, MUFAs, or n-3 PUFAs (both as a group or as individuals). It could be that SFAs or MUFAs have low heritabilities because they are mainly modulated by the diet, which in this study contained very high proportions of these FAs. Additional information from multiple species may help us to better understand the genetic determinism in the FA transformation cascade and the respective influences of the elongase and desaturase enzymes in different contexts.

Strong genetic correlations, whether positive or negative, between different fatty acids is favourable to select the trait with the highest heritability. The genetic correlations estimated in the present study among FAs and traits associated with body weight, yield, and quality enable us to make predictions about changes that may occur as a result of selection on a given trait. Generally, our results indicated that improvement in one FA group is likely to cause unfavourable changes in other groups. In our study, the negative genetic correlations between PUFAs and SFAs/MUFAs demonstrate the trade-offs between groups. Generally, the genetic associations were synergic within FA groups but antagonistic between groups, which arises as a result of constraints inherent in the underlying biosynthesis pathways of these FAs. Additionally, even though our dataset of 1382 fish was larger than those used in previous studies, only some of the bivariate models for individual fatty acids converged (due to very high correlations between those pairs of traits) and all genetic correlations had large standard errors. Thus, we based our interpretation of the results on group traits rather than on individual estimates.

Another important contribution of this study is our characterisation of the genetic correlations between FAs and traits related to body weight, yield, and quality. In rainbow trout, and in fish in general, the main selection objectives are faster growth, increased yields, and disease resistance. Given the genetic correlations estimated here, genetic selection for increased BW or fat content would be expected to increase the proportions of SFAs and MUFAs and decrease the proportion of PUFAs; such a reduction in the abundance of healthful n-3 PUFAs could diminish the nutritional quality of flesh. Similar results have been observed in studies of common carp [33], Atlantic salmon [19], and tilapia [31], in which unfavourable genetic correlations were reported between performance traits and EPA or n-3 PUFAs. The fact that the same tendency we observed in visceral adipocytes has been reported in other studies based on fillets indicates that adipocytes probably exhibit similar relative variations in FA composition regardless of their origin.

One explanation of why faster growth or yield improvement could drive differential lipid composition may also be linked with nutrition. It is well known that the FA composition of diets strongly affects the composition of lipids in fish [16, 56,57,58]. In a 2018 study, trout that were selected for faster growth tended to have a higher feed intake; the fish then stored unutilised energy as MUFAs, which shifted the balance of FAs in flesh away from PUFAs [33]. Here, we identified a moderate positive genetic correlation between fillet yield (predicted by HGCarc%) or carcass yield (Carc%) and certain FA traits (DHA and sum EPA + DHA). This could be linked to the fact that, compared to oleic acid (18:1 n-9), EPA and DHA lower the accumulation of triglycerides in adipocytes, thus limiting the development of adipose tissue in fish [54, 59]. Carcass yield is inversely correlated to visceral yield (with increased visceral tissue, carcass yield declines), so it is likely that as the amount of visceral tissue increases, we would see a corresponding decline in EPA and DHA content. Our results give support to this hypothesis on the link between EPA and DHA and oleic acid storage and the interest to limit fat deposition in the muscle and/or the visceral fat tissue in selecting leaner fish with higher fillet yield. Thus selecting for increasing carcass yield should be a good strategy which should preserve or increase EPA and DHA content, therefore maintaining a favourable balance of FAs and the accompanying nutritional benefits for humans [55]. Such an approach might be preferable to one based only on growth, which is negatively correlated to n-3 LC PUFA content. However, selecting for both carcass yield and EPA/DHA content would reinforce the impacts on FA composition.

New QTLs and candidate genes for fatty acid composition

All previous studies on the genetic basis of FA composition in aquaculture species have described a polygenic architecture, with a few QTLs responsible for moderate levels of genetic variation in FAs [28, 32, 33, 44, 60].

Our study identified QTLs on Omy7 that were associated with n-6 PUFAs, LA, and EPA; contained within this region was the gene onmy-cd8a (T-cell surface glycoprotein CD8 alpha precursor), which is indirectly related to FA activity. In Atlantic salmon, a gene in the same family (LOC106581970, T-cell surface glycoprotein CD3 zeta chain-like, on chromosome 21) was also identified as a potential candidate gene with links to EPA and DHA [28]. It has been hypothesised that EPA and DHA are able to modulate T-cell activation to exert anti-inflammatory influences [61,62,63]. Within this same region on Omy7, we also identified mrps9 (28S ribosomal protein S9, mitochondrial isoform) near the peak SNP. Members of the MRPS family are involved in the synthesis of protein inside mitochondria [64], with one of their roles being the induction of apoptosis by SFAs in several cell types [65, 66]. Another member of this family, MRPS30 (mitochondrial ribosomal protein S30), was highlighted by a GWAS analysis on FA composition in sheep [64]; this gene was associated with myristic acid (C14:0) content and the ratio of n-6 to n-3 PUFAs.

On Omy1, we detected QTLs associated with n-6 PUFAs, MUFAs, LA, and EPA; close to the peak SNP in this region, we identified three interesting genes with indirect links to FAs: apmap, accs1, and abhd12. Apmap (adipocyte plasma membrane associated protein) has arylesterase and strictosidine synthase activity and may play a role in adipocyte differentiation [67], while acss1 (acetyle coenzyme A synthetase 2-like) catalyses the synthesis of acetyl-CoA from short-chain FAs [66]. Abhd12 (Lysophospholipase ABHD12) is known to be involved in immune and neurological processes and plays a role in the regulation of lysophospatidylserine pathways and related with very-long chain lipids [68,69,70].

In the ARA-associated QTL on Omy10, the gene acsl3a (long-chain fatty acid CoA ligase 3-like) was present close to the peak SNP [71]. This gene belongs to the same family as Acyl-coenzyme A (CoA) synthetase 1 (acsl1), a well-studied obesogenic gene involved in FA metabolism that is associated with high caloric food intake in mice and humans [72]. In general, members of this family appear to play similar roles: acsl5 has been implicated in lipid biosynthesis and FA degradation [73], and acsl1 and acsl5 were linked with ARA and ARA/ALA, respectively, in a study of common carp [42]. Another gene of interest present in this region was mogat3b (monoacylglycerol O-acyltransferase 3b), which is predicted to have diacylglycerol O-transferase activity; it also plays a role in triglyceride biosynthesis and may be involved in the absorption of dietary fat [74,75,76]. An important paralog of this gene is dgat2 (diacylglycerol o-transferase 2), which has been implicated in the catalysis of the final stage of triacylglycerol biosynthesis [77] and linked to adipogenesis [78]. This gene was also highlighted in a study on sheep, in a QTL associated with SFAs, C18:0, C16:1, and MUFAs [79].

In a QTL associated with PUFAs on Omy12, we found the gene LOC110538527, which encodes the butyrophilin subfamily 1 member A1 protein. Butyrophilin is the main protein associated with milk fat droplets and milk quality in cattle [80], and several members of this family were highlighted in a GWAS analysis of FA composition in beef cattle [34].

The candidate gene TBC1D4 was found close to the peak SNP in a QTL on Omy22 associated with n-6 PUFAs. Members of the TBC1 domain family have been associated with insulin and FA composition in previous studies of humans and mice, and warrant further investigation [81,82,83]. In pigs, this gene family, and more specifically, TBC1D1, was also reported to have putative effects on three FA ratios (C16:1 n-7/C16:0, C18:1 n-9/C16:1 n-7, PUFA/MUFA), on PUFAs as a whole, as well as on n-3 and n-6 PUFAs [84].

Some previous studies of the genetic basis of FA content have identified candidate genes linked to FA metabolism, such as the elovl or fad genes in pigs [85] or elovl2 (involved in the conversion of DPA to DHA) in Atlantic salmon [28]. In our study, we did not detect any strong candidate genes that were directly involved in the bioconversion of PUFAs, n-3 PUFAs, EPA, or DHA in rainbow trout. However, we did note several candidate genes with links to EPA-, DHA-, or PUFA-related traits. Other candidate genes highlighted by our analysis included sdhaf4 (succinate dehydrogenase assembly factor 4, mitochondrial-like), free fatty acid receptor 2-like, atp11a (phospholipid-transporting ATPase IH isoform1), plpp2a (phospholipid phosphatase 2-like), isocitrate dehydrogenase, pemt (phosphatidylethanolamine N-methyltransferase), inpp4ab (Inositol polyphosphate-4-phosphatase type I Ab), mtmr2 (Myotubularin-related protein 2), ebpl (EBP-like), dhdds (alkyl transferase), and sstr5 (Somatostatin receptor 5). Several factors could explain the lack of association with n-3 PUFA or n-6 PUFA bioconversion pathway genes. One explanation could be that the SNPs on the array are not located in gene regions that influence lipid metabolism, or that any such genes located near the SNPs remain uncharacterised. It could also be that the diversity of EPA- and DHA-pathway genes is low, which would reduce our ability to detect QTLs. Finally, GWAS efficiency is influenced by the accuracy of phenotypic recording. It is possible that our decision to predict FA composition using Raman spectroscopy, which is less accurate than gas chromatography, might have influenced the result.

Genomic selection for fatty acids

Unlike pedigree-based BLUP selection, in which candidates from the same family have the same EBV, genomic selection (GS) is able to differentiate among candidates within a family and may thus yield improved estimates of breeding values and better selection efficiency. In particular, GS is more efficient than pedigree-based selection for traits which cannot be measured directly on selection candidates and for traits with limited heritability, such as FA proportions. Furthermore, GS has the additional advantage of requiring fewer phenotypes, since only the reference population must be phenotyped. This is especially beneficial for traits that are complex and expensive to measure, such as FA composition. In our study, GS for FA traits using a reference population of 1100 individuals was estimated to improve accuracy by 12 to 120% compared to BLUP selection. In the literature to date, only one study has compared the accuracy of selection for FA content in fish (in this case, Atlantic salmon) between pedigree- and genomic-based methods [86]. Those authors reported low to moderate GS accuracies (0.27 to 0.61), which are similar to the estimates we found here (0.34 to 0.70). For DHA, the pedigree and genomic accuracies were 0.33 and 0.41, respectively (a 26% gain in accuracy with GS, compared with 120% in our study), while for EPA, these values were 0.37 and 0.32, respectively (a higher accuracy for pedigree prediction; − 14% change in accuracy for GS compared to the 23% gain estimated in our study). In studies of cattle, the accuracy of GS was less than 0.40 for the majority of FAs [34, 36] and the reliability of genomic prediction for milk FA composition from three Holstein populations was also less than 0.40 [87]; despite this, though, genomic prediction was always superior to pedigree prediction. In our study, the two lowest values of genomic accuracy were obtained for n-3 PUFAs and ALA, which were also the two FAs with the lowest heritabilities. Compared to traditional pedigree-based selection, GS seems to demonstrate particular potential for the healthy, desirable FAs EPA, DHA, and EPA + DHA (gains of + 23%, + 120%, and + 70%, respectively).

Genotyping and data collection are costly, and the relative advantage of using SNP data in selection ultimately depends on how these costs are offset by the value of the improvement in traits of interest. The higher genetic gains enabled by GS with respect to pedigree selection may partially cover the extra costs of genotyping [88]. Because of this, it may be feasible to incorporate GS for FA traits in breeding programmes in order to select fish with a superior genetic basis for FA composition. The major advantage of GS is the ability to differentiate among candidates within a family, which here resulted in a gain in accuracy between 12% and 120%; this advantage is particularly pronounced for traits that cannot be measured directly on a selection candidate. Furthermore, the cost of implementation of GS could be reduced by optimising the density of the SNP panel in use. Several studies have evaluated the accuracy of predictions made using more cost-effective lower-density SNP chips, and, for the traits examined, the use of 500-SNP panels still provided predictive accuracy that was higher than that of BLUP [89, 90].

Conclusions

In summary, this work provides new insights into the genetics of fatty acid traits in rainbow trout. Our results reveal that fatty acid proportions are very polygenic traits, but, under the conditions investigated here, most appear to be moderately heritable, with the highest heritability observed for n-6 PUFAs (0.24). We detected several genomic regions that explained up to 2% of the genetic variance in proportions of MUFAs, PUFAs, n-6 PUFAs, LA, ARA, and EPA. When we investigated these regions, we identified several genes (mrps9, mogat3b, TBC1D4, acsl3, onmy-cd8a, butyrophiulin family, apmap, acss1, and abhd12) that can be indirectly implicated in fatty acid metabolism. These genes represent good candidates for further functional validation to decipher the biological mechanisms underlying variation in fatty acid traits in rainbow trout. The work also provides new information allowing to propose that selection on carcass or fillet yields should improve n-3 LC PUFA composition of the fillet, but a combined approach based on both yield and FA composition in the fillet should further increase the efficiency of selection. Finally, our analyses indicate that, with a reference population of about 1100 individuals, the implementation of genomic selection in a breeding programme for fatty acid traits would enable a gain in accuracy of 12–120% compared to standard pedigree-based selection. These results suggest that genomic evaluation is a feasible strategy for selecting trout with superior genetic merit for traits related to production, quality, and fatty acid composition.

Methods

Fish production and trait recording

The fish used in this study were derived from a commercially selected line from the Sources de l’Avance breeding company, a subsidiary of Aqualande Group (Pissos, France). The line has been previously selected for 9 generations with a multi-trait selection combining mass selection on growth and carcass yield assisted by ultrasound and morphometry and sib selection on carcass yield and fillet yield [91].

Through 10 factorial crosses, 84 dams were crossed with 99 neomales (sex-reversed females used as sires) on the same day to create 831 families of rainbow trout. A piece of fin was sampled from each parent for DNA extraction and subsequent genomic analysis. Trout were reared under commercial conditions in the “Viviers de la Houtine” growing farm (Belin-Beliet, France); further details on rearing can be found in [46]. Fish were fed to satiation using extruded commercial feed: Neo start (17% lipids) and Neo CDC (23% lipids) (Le Gouessant, Lamballe, France) during the first stage, then Extra CDC AQL G25 (25% lipids) (Le Gouessant, Lamballe, France) and Viva Pro 7F NAT29 (30% lipids) (Aqualia, Arue, France) until the end of the experiment. The composition of the feed mimicked that used in commercial operations: the majority of fish meal and fish oil had been substituted with plant-based ingredients, so that PUFA levels were minimal (detailed in Additional file 1). Fish were reared following standard practices, and measurements were performed only on slaughtered fish; there was thus no need to consult an ethics committee. At 469 days post-fertilisation (dpf), fish were individually tagged with RFID transponders and fin samples were collected and preserved in 95% ethanol for DNA extraction and genomic analysis. At the end of the growing period (between 503 and 506 dpf), data were collected from 1410 fish randomly sampled. Fish were humanely killed by a blow to the head and bled by cutting the gills in an ice water bath, in accordance with good animal slaughtering practices. Post-mortem data collection and processing were accomplished as quickly as possible to ensure data accuracy. All measured traits were defined according to the ATOL (Animal Trait Ontology for Livestock) database, available online (https://www.atol-ontology.com/en/). From each fish, the following traits were measured: body weight (BW, ATOL_0000351), body length (BL, ATOL_0001658), head weight (HeadW, ATOL_0001545), headless gutted carcass weight (HGCarcW, ATOL_0002260), and viscera weight (ViscW, ATOL_0002258). These traits were combined to calculate three synthetic traits or processing yields: Fulton coefficient condition, calculated as K = BW(g)100/BL3 (cm) (ATOL_0001653); headless gutted carcass yield (HGCarc%, ATOL_0002261), an indirect predictor of fillet yield (rg = 0.97 [92];); and gutted carcass yield (Carc%, ATOL_0000548). Total fat content in muscle (Fat, ATOL_0001663) was estimated indirectly using a Fish Torry Fat-meter® positioned on the skin as described in [93]. For each fish, one steak was cut in front of the dorsal fin and photographed using a digital camera (Canon EOS 1000 D 10 M Pixels, adapted with a shooting tent (Literoom Photoflex©) to avoid specular reflection and a copy stand to fix the camera. The steaks were then packed in individual plastic bags and frozen at − 20 °C until magnetic resonance imaging (MRI) analysis. The digital pictures were analysed using a modification of the method described by Marty-Mahe et al. (2003) [94]. First, the image was L*a*b* transformed, then colour image segmentation was performed with Visilog 7.3 for Windows© to quantify areas of the steak (ATOL_0005553). Flesh colour (ATOL_0001017) (including myosepta) was expressed in the L*, a*, b* system—representing luminosity, redness and yellowness, respectively—as recommended by the CIE (CIE 1976) [95].

The process of magnetic resonance image formation is based on a combination of permanent and radio-frequency (RF) magnetic fields [96]. An automatic image analysis scheme can then be used to distinguish flesh from subcutaneous fat tissue. The method for mapping the distribution of fat is detailed in [97]. For each steak, MRI measurements were used to determine total fat content (MRI_F%), subcutaneous fat proportion (MRI_F_SC%), and fat content in flesh (IRM_F_F%). The genetic architecture of all of these traits has already been analysed, and further details on trait measurement can be found in additional Table 4 and Blay et al. (2021) [46]. Data from that study are used here to estimate the phenotypic and genetic correlations of production traits with fatty acids.

Prediction of fatty acid composition

The FA composition of adipocytes from visceral fat (ATOL_0000074) was predicted by Raman spectroscopy (LabRAM HR800, Horiba Scientific) following the calibration method developed by [Prado E, Eklouh-Molinier C, Enez F, Causeur D, Blay C., Dupont-Nivet M, Labbé L, Petit V, Moreac A, Taupier G, Haffray P, Bugeon J, Corraze G, Nazabal, V. Prediction of fatty acids composition in the rainbow trout Oncorhynchus mykiss by using Raman micro-spectroscopy. Analytical Chimica Acta, submitted], with gas chromatography (Eurofins Analytics, Nantes) as a reference method [98]. The detection limit for gas chromatography was estimated to be 0.05% of total FAs. This is detailed in [Prado E, Eklouh-Molinier C, Enez F, Causeur D, Blay C., Dupont-Nivet M, Labbé L, Petit V, Moreac A, Taupier G, Haffray P, Bugeon J, Corraze G, Nazabal, V. Prediction of fatty acids composition in the rainbow trout Oncorhynchus mykiss by using Raman micro-spectroscopy, submitted] but briefly (see also Fig. 1), the calibration was made from 1) 259 individuals divided into two groups fed with two different diets, one enriched in marine fish oil and fish meal and one substituted with plant ingredients 2) nine additional individuals from the 1410 individuals of this study. Visceral adipose tissues were analysed by gas chromatography and also by Raman spectroscopy for 268 samplesand only with Raman for the rest of the samples (1401).

Visceral fat (2 g), mostly composed of adipocytes, was collected from the “front” lobe of visceral adipose tissue from 1410 commercial trout and preserved in liquid nitrogen. Raman measurements were acquired on 1410 sibs (including the 9 used for the calibration) using a 10x objective and an excitation wavelength of 785 nm. Two spectral ranges were recorded: 550 to 1800 cm− 1 and 2610 to 3100 cm− 1. The calibration by means of a ridge regression model was used to predict FA composition in the 1410 rainbow trout. Specifically, the proportions of saturated fatty acids (SFAs), mono-unsaturated fatty acids (MUFAs), poly-unsaturated fatty acids (PUFAs), omega-3 (n-3 PUFAs) and omega-6 (n-6 PUFAs) fatty acids, and several individual fatty acids were calculated as a percentage of total FA content. The acquisition parameters, pre-processing treatment, and statistical analyses are described in detail by Prado et al. (2021) [Prado E, Eklouh-Molinier C, Enez F, Causeur D, Blay C., Dupont-Nivet M, Labbé L, Petit V, Moreac A, Taupier G, Haffray P, Bugeon J, Corraze G, Nazabal, V. Prediction of fatty acids composition in the rainbow trout Oncorhynchus mykiss by using Raman micro-spectroscopy, submitted.].

Genotyping

Fish were genotyped using the 57 K SNP Axiom® Trout Genotyping array from Thermofisher [99] at the INRAE genotyping platform Gentyane (Clermont-Ferrand, France).

SNPs were subjected to several steps of quality-control filtering, as described in D’Ambrosio et al. (2019) [100], in particular to remove SNPs with probe polymorphism and multiple locations on the genome. From the initial set of 57 K SNPs, we retained those with a call rate higher than 0.97, no significant deviation from Hardy-Weinberg equilibrium (p-value > 0.0001), and a minor allele frequency (MAF) higher than 0.05 for further analysis. In addition, samples in which less than 90% of SNPs were genotyped were removed from the analysis. Finally, all missing SNP genotypes of the remaining individuals were imputed using FImpute software 2.0 [101]. A total of 29,652 SNPs passed quality-control filtering. Following DNA quality control, our final dataset contained genotypes for 1382 of the 1410 fish that had been phenotyped for FA traits.

Estimation of genetic parameters

Heritability (h2) and phenotypic and genetic correlations (rp, and rg, respectively) were estimated using the restricted maximum likelihood method (AIREML) and BLUPF90 software [102]. Univariate analyses were performed to estimate the heritability of all traits. Bivariate analyses were performed to estimate genetic correlations between traits, using the following animal model:

$$ {Y}_i=\upmu +{\mathrm{a}}_i+{\mathrm{e}}_i $$

where Yi is the performance of the ith animal, μ is the overall mean of the population, ai is the additive effect of the ith animal, and ei is the residual random error term. The pedigree under consideration contained 17,235 animals over 9 generations. The maternal variance was not significantly different from zero and thus maternal effects were not included in the final models for all traits. All fish were female and reared in the same raceway, so there was no need for a fixed effect to be included in the model.

Heritability estimates were calculated as the ratio of additive genetic variance (Va) divided by the total phenotypic variance (Vp). For the genetic correlation among FAs, data were transformed using a logarithmic function to improve the convergence of estimates.

GWAS and QTL detection

GWAS was performed using BayesCπ [103], a Bayesian Stochastic Search Variable Selection approach, implemented in BESSiE software (version 1.0) [104].

In the BayesCπ model, only a certain proportion of SNPs (π) are assumed to have a non-zero effect on the phenotype. The marker effects are estimated through an MCMC algorithm that considers a mixture of markers, of which proportion π have effects that follow a normal distribution N (0, σ2a) and proportion 1 – π have zero effect. The general model used is:

$$ {Y}_i=\mu +\sum \limits_{j=1}^n{\delta}_{jk}{a}_j{g}_{ij}+{\varepsilon}_{ik} $$

with Y the phenotype observed for the ith individual, μ the overall mean in the population, n the total number of SNPs in the analysis, aj the additive effect of the reference allele for the jth SNP, gij the genotype for individual i (coded as 0, 1, or 2), and εik the residual effect for the ith individual in the kth iteration. The vector of residual effects is normally and independently distributed, \( \boldsymbol{\varepsilon} \sim N\left(0,\boldsymbol{I}{\sigma}_e^2\right), \) with σ2e the residual variance.

A total of 200,000 cycles were used, with a burn-in period of 5000 cycles. Results were saved every 20 cycles. In order to check convergence, the MCMC algorithm was initiated three times with three different chains for the random number generator. Convergence was assessed by visual inspection of plots of the posterior density of genetic and residual variances and by high correlations (r > 0.99) between the genomic estimated breeding values (GEBVs) estimated from the different chains of the MCMC algorithm.

At each cycle k, the decision to include the jth SNP in the model depended on the indicator variable δjk: if δjk was equal to 1, the effect of the jth SNP was estimated as aj, while if δjk was equal to 0, no effect was estimated. This indicator variable was sampled from a binomial distribution with a probability π that δjk was equal to 1 (i.e. the SNP has a non-zero effect) and a probability 1-π that δjk was equal to 0. The proportion 1-π was sampled from a beta distribution, B(α, β), with α = 300 and β = 29,652; the value of π was kept almost constant at 1%, corresponding to approximately 300 SNPs selected at each cycle from the 29,652 markers. This value of π was considered to be a good compromise in our variable selection algorithm between the high degree of polymorphism of the traits under study and the limited number of individuals (n = 1382) in our dataset.

The degree of association between each SNP and a given phenotype was assessed using the Bayes Factor (BF): \( \mathrm{BF}=\frac{\raisebox{1ex}{${\mathrm{P}}_{\mathrm{i}}$}\!\left/ \!\raisebox{-1ex}{$\left(1-{\mathrm{P}}_{\mathrm{i}}\right)$}\right.}{\raisebox{1ex}{$\uppi $}\!\left/ \!\raisebox{-1ex}{$\left(1-\uppi \right)$}\right.} \), where Pi is the probability that the ith SNP has a non-zero effect.

Following Kass and Raftery (1995) [105], QTLs were evaluated based on calculations of 2*ln (BF) (twice the natural logarithm of the BF); the threshold 2*ln (BF) ≥ 6 was considered evidence for a QTL. As proposed by Michenet et al. (2016) [106], a credibility interval was constructed around the peak SNP that encompassed, within a sliding window of 1 Mb on either side, all neighbouring SNPs for which 2*ln (BF) ≥ 3. If 5 ≤ 2*ln (BF) < 6 for a peak, the QTL was considered putative, and was included in further analyses only if it explained at least 0.5% of the genetic variance for a trait or if it was linked with more than one trait.

All candidate genes that were located within the confidence or credibility intervals established using the BayesCπ approach are listed in Additional file 3, with annotation from the NCBI Oncorhynchus mykiss genome assembly release 100 (GCF_002163495.1), and gene symbols from Lallias et al. (2020) [107].

Selection efficiency

To assess the relative efficiencies of a pedigree-based selection approach (best linear unbiased predictor; BLUP) and a genomic selection approach (GBLUP), estimated breeding values (EBV and GEBV, respectively) were derived via consideration of either the pedigree relationship matrix (A) or the genomic matrix (G) using the software package BLUPf90 [102].

For pedigree-based evaluation, the following BLUP model was applied for the estimation of breeding values:

$$ {Y}_i=\upmu +{Z}_i{\mathrm{a}}_i+{\mathrm{e}}_i $$

where Yi is the performance of the ith animal, μ is the overall mean of the population, and ai and ei are the vectors of additive genetic effects and residual effects that explain the performance of all phenotyped animals, respectively. Z is the incidence matrix for ai. In our model, vector ai corresponded to the breeding values of 17,235 individuals related through the pedigree relationship matrix A of the 1382 phenotyped fish.

For genomic selection, the genomic relationship matrix G was used in place of the pedigree matrix A [108]. The following GBLUP model was applied for the estimation of genomic breeding values (GEBV):

$$ {Y}_i=\upmu +{Z}_i{\mathrm{g}}_i+{\mathrm{e}}_i $$

with the vector gi corresponding to the breeding values of 1382 phenotyped and genotyped individuals related through the genomic relationship matrix G.

To assess the accuracy of (G) EBVs, 40 replicates of Monte Carlo ‘leave-one-group-out’ cross validation tests [109] were performed. For each replicate, 314 fish from the 1382 phenotyped and genotyped individuals were randomly chosen for the validation set and 1068 fish were chosen for the training set. The phenotypes recorded in the validation population were then masked and breeding values were estimated using (G) BLUP model.

For each replicate, the selection accuracy (Acc) was computed as:

$$ Acc=r\left((G){EBV}_v,y\right)/\surd h2 $$

where r represents the correlation between (G) EBVv and y, (G) EBVv represents the (genomic) breeding values of individuals belonging to the validation dataset, y is their phenotype, and h2 is the heritability estimated using pedigree information. In addition to evaluating selection accuracy, we also assessed the quality of BLUP and GBLUP evaluation by deriving the inflation coefficient of EBVs as a measure of selection bias. The inflation coefficient is the slope of the regression of the phenotypes on the (G)EBVs. In the absence of selection bias, this coefficient is expected to be equal to 1; in case of EBV over-dispersion (inflation), the coefficient is below 1 and in the case of EBV under-dispersion the value is above 1.

The values obtained for selection accuracy and the inflation coefficient were averaged over the 40 replicates.

Availability of data and materials

The data that support the findings of this study belong to the breeding company Aqualande; they were used under license for the current study, and so are not publicly available. The data can be made available for the purpose of reproduction of the results on request via a material transfer agreement and with the permission of Aqualande. Requests to access the datasets should be directed to Mathilde Dupont-Nivet or Vincent Petit.

References

  1. 1.

    FAO. La situation mondiale des pêches et de l’aquaculture 2020. La durabilité en action. Rome: FAO; 2020. https://doi.org/10.4060/ca9229fr.

    Book  Google Scholar 

  2. 2.

    Torstensen BE, Bell JG, Rosenlund G, Henderson RJ, Graff IE, Tocher DR, et al. Tailoring of Atlantic salmon (Salmo salar L.) flesh lipid composition and sensory quality by replacing fish oil with a vegetable oil blend. J Agric Food Chem. 2005;53(26):10166–78. https://doi.org/10.1021/jf051308i.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Calder PC, Yaqoob P. Omega-3 polyunsaturated fatty acids and human health outcomes. BioFactors. 2009;35(3):266–72. https://doi.org/10.1002/biof.42.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Ruxton CHS, Calder PC, Reed SC, Simpson MJA. The impact of long-chain n −3 polyunsaturated fatty acids on human health. Nutr Res Rev. 2005;18(1):113–29. https://doi.org/10.1079/nrr200497.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Tocher DR. Metabolism and functions of lipids and fatty acids in teleost fish. Rev Fish Sci. 2003;11(2):107–84. https://doi.org/10.1080/713610925.

    CAS  Article  Google Scholar 

  6. 6.

    Turchini GM, Torstensen BE, Ng WK. Fish oil replacement in finfish nutrition. Rev Aquac. 2009;1(1):10–57. https://doi.org/10.1111/j.1753-5131.2008.01001.x.

    Article  Google Scholar 

  7. 7.

    Abdul Hamid NK, Carmona-Antoñanzas G, Monroig Ó, Tocher DR, Turchini GM, Donald JA. Isolation and functional characterisation of a fads2 in rainbow trout (Oncorhynchus mykiss) with Δ5 desaturase activity. PLoS One. 2016;11(3):e0150770. https://doi.org/10.1371/journal.pone.0150770.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Gregory MK, James MJ. Rainbow trout (Oncorhynchus mykiss) Elovl5 and Elovl2 differ in selectivity for elongation of omega-3 docosapentaenoic acid. Biochim Biophys Acta. 2014;1841(12):1656–60. https://doi.org/10.1016/j.bbalip.2014.10.001.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Gregory MK, Collins RO, Tocher DR, James MJ, Turchini GM. Nutritional regulation of long-chain PUFA biosynthetic genes in rainbow trout (Oncorhynchus mykiss). Br J Nutr. 2016;115(10):1721–9. https://doi.org/10.1017/S0007114516000830.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Grigorakis K. Compositional and organoleptic quality of farmed and wild gilthead sea bream (Sparus aurata) and sea bass (Dicentrarchus labrax) and factors affecting it: a review. Aquaculture. 2007;272(1-4):55–75. https://doi.org/10.1016/j.aquaculture.2007.04.062.

    Article  Google Scholar 

  11. 11.

    Fonseca-Madrigal J, Bell JG, Tocher DR. Nutritional and environmental regulation of the synthesis of highly unsaturated fatty acids and of fatty-acid oxidation in Atlantic salmon (Salmo salar L.) enterocytes and hepatocytes. Fish Physiol Biochem. 2006;32(4):317–28. https://doi.org/10.1007/s10695-006-9109-2.

    CAS  Article  Google Scholar 

  12. 12.

    Zheng X, Tocher DR, Dickson CA, Bell JG, Teale AJ. 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. 2004;236(1-4):467–83. https://doi.org/10.1016/j.aquaculture.2004.02.003.

    CAS  Article  Google Scholar 

  13. 13.

    Stubhaug I, Tocher DR, Bell JG, Dick JR, Torstensen BE. Fatty acid metabolism in Atlantic salmon (Salmo salar L.) hepatocytes and influence of dietary vegetable oil. Biochim Biophys Acta - Mol Cell Biol Lipids. 2005;1734(3):277–88. https://doi.org/10.1016/j.bbalip.2005.04.003.

    CAS  Article  Google Scholar 

  14. 14.

    Johansson L, Kiessling A, Kiessling KH, Berglund L. Effects of altered ration levels on sensory characteristics, lipid content and fatty acid composition of rainbow trout (Oncorhynchus mykiss). Food Qual Prefer. 2000;11(3):247–54. https://doi.org/10.1016/S0950-3293(99)00073-7.

    Article  Google Scholar 

  15. 15.

    Turchini GM, Francis DS. Fatty acid metabolism (desaturation, elongation and β-oxidation) in rainbow trout fed fish oil- or linseed oil-based diets. Br J Nutr. 2009;102(1):69–81. https://doi.org/10.1017/S0007114508137874.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Lazzarotto V, Médale F, Larroquet L, Corraze G. Long-term dietary replacement of fishmeal and fish oil in diets for rainbow trout (Oncorhynchus mykiss): effects on growth, whole body fatty acids and intestinal and hepatic gene expression. PLoS One. 2018;13(1):1–25. https://doi.org/10.1371/journal.pone.0190730.

    CAS  Article  Google Scholar 

  17. 17.

    Kamalam BS, Médale F, Larroquet L, Corraze G, Panserat S. Metabolism and fatty acid profile in fat and lean rainbow trout lines fed with vegetable oil: effect of carbohydrates. PLoS One. 2013;8(10):e76570. https://doi.org/10.1371/journal.pone.0076570.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Yıldız M, Eroldoğan TO, Ofori-Mensah S, Engin K, Baltacı MA. The effects of fish oil replacement by vegetable oils on growth performance and fatty acid profile of rainbow trout: re-feeding with fish oil finishing diet improved the fatty acid composition. Aquaculture. 2018;488:123–33. https://doi.org/10.1016/j.aquaculture.2017.12.030.

    CAS  Article  Google Scholar 

  19. 19.

    Leaver MJ, Taggart JB, Villeneuve L, Bron JE, Guy DR, Bishop SC, et al. Heritability and mechanisms of n- 3 long chain polyunsaturated fatty acid deposition in the flesh of Atlantic salmon. Comp Biochem Physiol - Part D Genomics Proteomics. 2011;6(1):62–9. https://doi.org/10.1016/j.cbd.2010.04.002.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Horn SS, Ruyter B, Meuwissen THE, Hillestad B, Sonesson AK. Genetic effects of fatty acid composition in muscle of Atlantic salmon. Genet Sel Evol. 2018;50(1):1–12. https://doi.org/10.1186/s12711-018-0394-x.

    CAS  Article  Google Scholar 

  21. 21.

    Morais S, Pratoomyot J, Taggart JB, Bron JE, Guy DR, Bell JG, et al. Genotype-specific responses in Atlantic salmon (Salmo salar) subject to dietary fish oil replacement by vegetable oil: a liver transcriptomic analysis. BMC Genomics. 2011;12(1):255. https://doi.org/10.1186/1471-2164-12-255.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Dong L, Xiao S, Wang Q, Wang Z. Comparative analysis of the GBLUP, emBayesB, and GWAS algorithms to predict genetic values in large yellow croaker (Larimichthys crocea). BMC Genomics. 2016;17(1):1–11. https://doi.org/10.1186/s12864-016-2756-5.

    CAS  Article  Google Scholar 

  23. 23.

    Kamalam BS, Panserat S, Aguirre P, Geurden I, Fontagné-Dicharry S, Médale F. Selection for high muscle fat in rainbow trout induces potentially higher chylomicron synthesis and PUFA biosynthesis in the intestine. Comp Biochem Physiol - A Mol Integr Physiol. 2013;164(2):417–27. https://doi.org/10.1016/j.cbpa.2012.11.020.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Kamalam BS, Medale F, Kaushik S, Polakof S, Skiba-Cassy S, Panserat S. Regulation of metabolism by dietary carbohydrates in two lines of rainbow trout divergently selected for muscle fat content. J Exp Biol. 2012;215(15):2567–78. https://doi.org/10.1242/jeb.070581.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Quillet E, Le Guillou S, Aubin J, Fauconneau B. Two-way selection for muscle lipid content in pan-size rainbow trout (Oncorhynchus mykiss). Aquaculture. 2005;245(1-4):49–61. https://doi.org/10.1016/j.aquaculture.2004.12.014.

    CAS  Article  Google Scholar 

  26. 26.

    Kolditz C, Borthaire M, Richard N, Corraze G, Panserat S, Vachot C, et al. Liver and muscle metabolic changes induced by dietary energy content and genetic selection in rainbow trout (Oncorhynchus mykiss). Am J Physiol - Regul Integr Comp Physiol. 2008;294(4):R1154–64. https://doi.org/10.1152/ajpregu.00766.2007.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Kolditz CI, Paboeuf G, Borthaire M, Esquerré D, SanCristobal M, Lefèvre F, et al. Changes induced by dietary energy intake and divergent selection for muscle fat content in rainbow trout (Oncorhynchus mykiss), assessed by transcriptome and proteome analysis of the liver. BMC Genomics. 2008;9(1):506. https://doi.org/10.1186/1471-2164-9-506.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Horn SS, Ruyter B, Meuwissen THE, Moghadam H, Hillestad B, Sonesson AK. GWAS identifies genetic variants associated with omega-3 fatty acid composition of Atlantic salmon fillets. Aquaculture. 2020;514:734494. https://doi.org/10.1016/j.aquaculture.2019.734494.

    CAS  Article  Google Scholar 

  29. 29.

    Bell JG, Pratoomyot J, Strachan F, Henderson RJ, Fontanillas R, Hebard A, et al. Growth, flesh adiposity and fatty acid composition of Atlantic salmon (Salmo salar) families with contrasting flesh adiposity: effects of replacement of dietary fish oil with vegetable oils. Aquaculture. 2010;306(1-4):225–32. https://doi.org/10.1016/j.aquaculture.2010.05.021.

    CAS  Article  Google Scholar 

  30. 30.

    Schlechtriem C, Bron JE, Tocher DR. Inter-individual variation in total fatty acid compositions of flesh of Atlantic salmon smolts-fed diets containing fish oil or vegetable oil. Aquac Res. 2007;38(10):1045–55. https://doi.org/10.1111/j.1365-2109.2007.01759.x.

    CAS  Article  Google Scholar 

  31. 31.

    Nguyen NH, Ponzoni RW, Yee HY, Abu-Bakar KR, Hamzah A, Khaw HL. Quantitative genetic basis of fatty acid composition in the GIFT strain of Nile tilapia (Oreochromis niloticus) selected for high growth. Aquaculture. 2010;309(1-4):66–74. https://doi.org/10.1016/j.aquaculture.2010.08.034.

    CAS  Article  Google Scholar 

  32. 32.

    Xia JH, Lin G, He X, Yunping B, Liu P, Liu F, et al. Mapping quantitative trait loci for Omega-3 fatty acids in Asian seabass. Mar Biotechnol. 2014;16(1):1–9. https://doi.org/10.1007/s10126-013-9524-1.

    CAS  Article  Google Scholar 

  33. 33.

    Prchal M, Vandeputte M, Gela D, Doležal M, Buchtová H, Rodina M, et al. Estimation of genetic parameters of fatty acids composition in flesh of market size common carp (cyprinus Carpio L.) and their relation to performance traits revealed that selective breeding can indirectly affect flesh quality. Czech J Anim Sci. 2018;63(7):280–91. https://doi.org/10.17221/30/2018-CJAS.

    CAS  Article  Google Scholar 

  34. 34.

    Chen L, Ekine-Dzivenu C, Vinsky M, Basarab J, Aalhus J, Dugan MER, et al. Genome-wide association and genomic prediction of breeding values for fatty acid composition in subcutaneous adipose and longissimus lumborum muscle of beef cattle. BMC Genet. 2015;16(1):135. https://doi.org/10.1186/s12863-015-0290-0.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Lemos MVA, Chiaia HLJ, Berton MP, Feitosa FLB, Aboujaoud C, Camargo GMF, et al. Genome-wide association between single nucleotide polymorphisms with beef fatty acid profile in Nellore cattle using the single step procedure. BMC Genomics. 2016;17(1):213. https://doi.org/10.1186/s12864-016-2511-y.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Zhu B, Niu H, Zhang W, Wang Z, Liang Y, Guan L, et al. Genome wide association study and genomic prediction for fatty acid composition in Chinese Simmental beef cattle using high density SNP array. BMC Genomics. 2017;18(1):464. https://doi.org/10.1186/s12864-017-3847-7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Bouwman AC, Bovenhuis H, Visker MHPW, van Arendonk JAM. Genome-wide association of milk fatty acids in Dutch dairy cattle. BMC Genet. 2011;12(1):1–12. https://doi.org/10.1186/1471-2156-12-43.

    CAS  Article  Google Scholar 

  38. 38.

    Puig-Oliveras A, Revilla M, Castelló A, Fernández AI, Folch JM, Ballester M. Expression-based GWAS identifies variants, gene interactions and key regulators affecting intramuscular fatty acid content and composition in porcine meat. Sci Rep. 2016;6(1):31803. https://doi.org/10.1038/srep31803,.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    van Son M, Enger EG, Grove H, Ros-Freixedes R, Kent MP, Lien S, et al. Genome-wide association study confirm major QTL for backfat fatty acid composition on SSC14 in Duroc pigs. BMC Genomics. 2017;18(1):1–13. https://doi.org/10.1186/s12864-017-3752-0.

    CAS  Article  Google Scholar 

  40. 40.

    Yang B, Zhang W, Zhang Z, Fan Y, Xie X, Ai H, et al. Genome-wide association analyses for fatty acid composition in porcine muscle and abdominal fat tissues. PLoS One. 2013;8(6):e65554. https://doi.org/10.1371/journal.pone.0065554.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Zhang Y, Zhang J, Gong H, Cui L, Zhang W, Ma J, et al. Genetic correlation of fatty acid composition with growth, carcass, fat deposition and meat quality traits based on GWAS data in six pig populations. Meat Sci. 2019;150:47–55. https://doi.org/10.1016/j.meatsci.2018.12.008.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Zhang H, Xu P, Jiang Y, Zhao Z, Feng J, Tai R, et al. Genomic, transcriptomic, and epigenomic features differentiate genes that are relevant for muscular polyunsaturated fatty acids in the common carp. Front Genet. 2019;10:217. https://doi.org/10.3389/fgene.2019.00217.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Wang L, Chua E, Sun F, Wan ZY, Ye B, Pang H, et al. Mapping and validating QTL for fatty acid compositions and growth traits in Asian seabass. Mar Biotechnol. 2019;21(5):643–54. https://doi.org/10.1007/s10126-019-09909-7.

    CAS  Article  Google Scholar 

  44. 44.

    Lin G, Wang L, Te Ngoh S, Ji L, Orbán L, Yue GH. Mapping QTL for Omega-3 content in hybrid saline Tilapia. Mar Biotechnol. 2018;20(1):10–9. https://doi.org/10.1007/s10126-017-9783-3.

    CAS  Article  Google Scholar 

  45. 45.

    Rye M, Gjerde B. Phenotypic and genetic parameters of body composition traits and flesh colour in Atlantic salmon, Salmo salar L. Aquac Res. 1996;27(2):121–33. https://doi.org/10.1111/j.1365-2109.1996.tb00976.x.

    Article  Google Scholar 

  46. 46.

    Blay C, Haffray P, Bugeon J, d’Ambrosio J, Dechamp N, Collewet G, et al. Genetic parameters and genome-wide association studies of quality traits characterized using imaging technologies in rainbow trout, Oncorhynchus mykiss. Front Genet. 2021;12:219. https://doi.org/10.3389/fgene.2021.639223.

    CAS  Article  Google Scholar 

  47. 47.

    Difford GF, Horn SS, Dankel KR, Ruyter B, Dagnachew BS, Hillestad B, et al. The heritable landscape of near-infrared and Raman spectroscopic measurements to improve lipid content in Atlantic salmon fillets. Genet Sel Evol. 2021;53(1):1–11. https://doi.org/10.1186/s12711-021-00605-6.

    CAS  Article  Google Scholar 

  48. 48.

    Olsen EF, Rukke EO, Flåtten A, Isaksson T. Quantitative determination of saturated-, monounsaturated- and polyunsaturated fatty acids in pork adipose tissue with non-destructive Raman spectroscopy. Meat Sci. 2007;76(4):628–34. https://doi.org/10.1016/j.meatsci.2007.02.004.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Killeen DP, Marshall SN, Burgess EJ, Gordon KC, Perry NB. Raman spectroscopy of fish oil capsules: polyunsaturated fatty acid quantitation plus detection of ethyl esters and oxidation. J Agric Food Chem. 2017;65(17):3551–8. https://doi.org/10.1021/acs.jafc.7b00099.

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    Bekhit MY, Grung B, Mjøs SA. Determination of omega-3 fatty acids in fish oil supplements using vibrational spectroscopy and chemometric methods. Appl Spectrosc. 2014;68(10):1190–200. https://doi.org/10.1366/13-07210.

    CAS  Article  PubMed  Google Scholar 

  51. 51.

    Afseth NK, Wold JP, Segtnan VH. The potential of Raman spectroscopy for characterisation of the fatty acid unsaturation of salmon. Anal Chim Acta. 2006;572(1):85–92. https://doi.org/10.1016/j.aca.2006.05.013.

    CAS  Article  PubMed  Google Scholar 

  52. 52.

    Overturf K, Welker T, Barrows F, Towner R, Schneider R, Lapatra S. Variation in rainbow trout, Oncorhynchus mykiss, to biosynthesize Eicosapentaenoic acid and docosahexaenoic acid when reared on plant oil replacement feeds. J World Aquac Soc. 2013;44(3):326–37. https://doi.org/10.1111/jwas.12041.

    CAS  Article  Google Scholar 

  53. 53.

    Nolasco-Alzaga HR, Perez-Enriquez R, Enez F, Bestin A, Palacios-Mechetnov E, Haffray P. Quantitative genetic parameters of growth and fatty acid content in the hemolymph of the Whiteleg shrimp Litopenaeus vannamei. Aquaculture. 2018;482:17–23. https://doi.org/10.1016/j.aquaculture.2017.09.015.

    CAS  Article  Google Scholar 

  54. 54.

    Salmerón C. Adipogenesis in fish. J Exp Biol. 2018;121(Suppl_1). https://doi.org/10.1242/jeb.161588.

  55. 55.

    Tocher DR. Omega-3 long-chain polyunsaturated fatty acids and aquaculture in perspective. Aquaculture. 2015;449:94–107. https://doi.org/10.1016/j.aquaculture.2015.01.010.

    CAS  Article  Google Scholar 

  56. 56.

    Caballero MJ, Obach A, Rosenlund G, Montero D, Gisvold M, Izquierdo MS. Impact of different dietary lipid sources on growth, lipid digestibility, tissue fatty acid composition and histology of rainbow trout, Oncorhynchus mykiss. Aquaculture. 2002;214(1-4):253–71. https://doi.org/10.1016/S0044-8486(01)00852-3.

    CAS  Article  Google Scholar 

  57. 57.

    Steffens W. Effects of variation in essential fatty acids in fish feeds on nutritive value of freshwater fish for humans. Aquaculture. 1997;151(1-4):97–119. https://doi.org/10.1016/S0044-8486(96)01493-7.

    CAS  Article  Google Scholar 

  58. 58.

    Turchini GM, Hermon KM, Francis DS. Fatty acids and beyond: fillet nutritional characterisation of rainbow trout (Oncorhynchus mykiss) fed different dietary oil sources. Aquaculture. 2018;491:391–7. https://doi.org/10.1016/j.aquaculture.2017.11.056.

    CAS  Article  Google Scholar 

  59. 59.

    Todorčević M, Vegusdal A, Gjøen T, Sundvold H, Torstensen BE, Kjær MA, et al. Changes in fatty acids metabolism during differentiation of Atlantic salmon preadipocytes; effects of n-3 and n-9 fatty acids. Biochim Biophys Acta - Mol Cell Biol Lipids. 1781;2008(6-7):326–35. https://doi.org/10.1016/j.bbalip.2008.04.014.

    CAS  Article  Google Scholar 

  60. 60.

    Wang W, Tan S, Luo J, Shi H, Zhou T, Yang Y, et al. GWAS analysis indicated importance of NF-κB signaling pathway in host resistance against motile Aeromonas septicemia disease in catfish. Mar Biotechnol. 2009;2019(3):335–47. https://doi.org/10.1007/s10126-019-09883-0.

    CAS  Article  Google Scholar 

  61. 61.

    Denys A, Hichami A, Khan NA. n-3 PUFAs modulate T-cell activation via protein kinase C-α and -ε and the NF-κB signaling pathway. J Lipid Res. 2005;46:752–8. https://doi.org/10.1194/jlr.M400444-JLR200.

    CAS  Article  PubMed  Google Scholar 

  62. 62.

    Lee JY, Plakidas A, Lee WH, Heikkinen A, Chanmugam P, Bray G, et al. Differential modulation of Toll-like receptors by fatty acids: Preferential inhibition by n-3 polyunsaturated fatty acids. J Lipid Res. 2003;44:479–86. https://doi.org/10.1194/jlr.M200361-JLR200.

    CAS  Article  PubMed  Google Scholar 

  63. 63.

    Rogero MM, Calder PC. Obesity, inflammation, toll-like receptor 4 and fatty acids. Nutrients. 2018;10(4):432. https://doi.org/10.3390/nu10040432.

    CAS  Article  PubMed Central  Google Scholar 

  64. 64.

    Kenmochi N, Suzuki T, Uechi T, Magoori M, Kuniba M, Higa S, et al. The human mitochondrial ribosomal protein genes: mapping of 54 genes to the chromosomes and implications for human disorders. Genomics. 2001;77(1-2):65–70. https://doi.org/10.1006/geno.2001.6622.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  65. 65.

    Wei Y, Wang D, Topczewski F, Pagliassotti MJ. Saturated fatty acids induce endoplasmic reticulum stress and apoptosis independently of ceramide in liver cells. Am J Physiol - Endocrinol Metab. 2006;291(2):E275–81. https://doi.org/10.1152/ajpendo.00644.2005.

    CAS  Article  PubMed  Google Scholar 

  66. 66.

    Tang DG, La E, Kern J, Kehrer JP. Fatty acid oxidation and signaling in apoptosis. Biol Chem. 2002;383(3-4):425–42. https://doi.org/10.1515/BC.2002.046.

    CAS  Article  PubMed  Google Scholar 

  67. 67.

    Ilhan A, Gartner W, Nabokikh A, Daneva T, Majdic O, Cohen G, et al. Localization and characterization of the novel protein encoded by C20orf3. Biochem J. 2008;414(3):485–95. https://doi.org/10.1042/BJ20080503.

    CAS  Article  PubMed  Google Scholar 

  68. 68.

    Ogasawara D, Ichu TA, Vartabedian VF, Benthuysen J, Jing H, Reed A, et al. Selective blockade of the lyso-PS lipase ABHD12 stimulates immune responses in vivo. Nat Chem Biol. 2018;14(12):1099–108. https://doi.org/10.1038/s41589-018-0155-8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  69. 69.

    Joshi A, Shaikh M, Singh S, Rajendran A, Mhetre A, Kamat SS. Biochemical characterization of the PHARC-associated serine hydrolase ABHD12 reveals its preference for very-long-chain lipids. J Biol Chem. 2018;293(44):16953–63. https://doi.org/10.1074/jbc.RA118.005640.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  70. 70.

    Kelkar DS, Ravikumar G, Mehendale N, Singh S, Joshi A, Sharma AK, et al. A chemical–genetic screen identifies ABHD12 as an oxidized-phosphatidylserine lipase. Nat Chem Biol. 2019;15(2):169–78. https://doi.org/10.1038/s41589-018-0195-0.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  71. 71.

    Gaudet P, Livstone MS, Lewis SE, Thomas PD. Phylogenetic-based propagation of functional annotations within the gene ontology consortium. Brief Bioinform. 2011;12(5):449–62. https://doi.org/10.1093/bib/bbr042.

    Article  PubMed  PubMed Central  Google Scholar 

  72. 72.

    Joseph R, Poschmann J, Sukarieh R, Too PG, Julien SG, Xu F, et al. ACSL1 is associated with fetal programming of insulin sensitivity and cellular lipid content. Mol Endocrinol. 2015;29(6):909–20. https://doi.org/10.1210/me.2015-1020.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  73. 73.

    Bowman TA, O’Keeffe KR, D’Aquila T, Yan QW, Griffin JD, Killion EA, et al. Acyl CoA synthetase 5 (ACSL5) ablation in mice increases energy expenditure and insulin sensitivity and delays fat absorption. Mol Metab. 2016;5(3):210–20. https://doi.org/10.1016/j.molmet.2016.01.001.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  74. 74.

    Cheng D, Nelson TC, Chen J, Walker SG, Wardwell-Swanson J, Meegalla R, et al. Identification of acyl coenzyme a:monoacylglycerol acyltransferase 3, an intestinal specific enzyme implicated in dietary fat absorption. J Biol Chem. 2003;278(16):13611–4. https://doi.org/10.1074/jbc.C300042200.

    CAS  Article  PubMed  Google Scholar 

  75. 75.

    Ma Z, Onorato JM, Chen L, Nelson DW, Yen CLE, Cheng D. Synthesis of neutral ether lipid monoalkyl-diacylglycerol by lipid acyltransferases®. J Lipid Res. 2017;58(6):1091–9. https://doi.org/10.1194/jlr.M073445.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  76. 76.

    Brandt C, McFie PJ, Stone SJ. Biochemical characterization of human acyl coenzyme a: 2-monoacylglycerol acyltransferase-3 (MGAT3). Biochem Biophys Res Commun. 2016;475(3):264–70. https://doi.org/10.1016/j.bbrc.2016.05.071.

    CAS  Article  PubMed  Google Scholar 

  77. 77.

    Yen CLE, Stone SJ, Koliwad S, Harris C, Farese RV. DGAT enzymes and triacylglycerol biosynthesis. J Lipid Res. 2008;49:2283–301. https://doi.org/10.1194/jlr.R800018-JLR200.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  78. 78.

    Keats EC, Dominguez JM, Grant MB, Khan ZA. Switch from canonical to noncanonical Wnt signaling mediates high glucose-induced adipogenesis. Stem Cells. 2014;32(6):1649–60. https://doi.org/10.1002/stem.1659.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  79. 79.

    Rovadoscki GA, Pertile SFN, Alvarenga AB, Cesar ASM, Pértille F, Petrini J, et al. Estimates of genomic heritability and genome-wide association study for fatty acids profile in Santa Inês sheep. BMC Genomics. 2018;19(1):1–14. https://doi.org/10.1186/s12864-018-4777-8.

    CAS  Article  Google Scholar 

  80. 80.

    Bhattacharya TK, Misra SS, Sheikh FD, Sukla S, Kumar P, Sharma A. Effect of butyrophilin gene polymorphism on milk quality traits in crossbred cattle. Asian-Australasian J Anim Sci. 2006;19(7):922–6. https://doi.org/10.5713/ajas.2006.922.

    CAS  Article  Google Scholar 

  81. 81.

    Cartee GD. Roles of TBC1D1 and TBC1D4 in insulin- and exercise-stimulated glucose transport of skeletal muscle. Diabetologia. 2015;58(1):19–30. https://doi.org/10.1007/s00125-014-3395-5.

    CAS  Article  PubMed  Google Scholar 

  82. 82.

    Maher AC, McFarlan J, Lally J, Snook LA, Bonen A. TBC1D1 reduces palmitate oxidation by inhibiting β-HAD activity in skeletal muscle. Am J Physiol - Regul Integr Comp Physiol. 2014;307(9):R1115–23. https://doi.org/10.1152/ajpregu.00014.2014.

    CAS  Article  PubMed  Google Scholar 

  83. 83.

    Treebak JT, Pehmøller C, Kristensen JM, Kjøbsted R, Birk JB, Schjerling P, et al. Acute exercise and physiological insulin induce distinct phosphorylation signatures on TBC1D1 and TBC1D4 proteins in human skeletal muscle. J Physiol. 2014;592(2):351–75. https://doi.org/10.1113/jphysiol.2013.266338.

    CAS  Article  PubMed  Google Scholar 

  84. 84.

    Zhang W, Yang B, Zhang J, Cui L, Ma J, Chen C, et al. Genome-wide association studies for fatty acid metabolic traits in five divergent pig populations. Sci Rep. 2016;6(1):24718. https://doi.org/10.1038/srep24718.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  85. 85.

    Zhang W, Zhang J, Cui L, Ma J, Chen C, Ai H, et al. Genetic architecture of fatty acid composition in the longissimus dorsi muscle revealed by genome-wide association studies on diverse pig populations. Genet Sel Evol. 2016;48(1):5. https://doi.org/10.1186/s12711-016-0184-2.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  86. 86.

    Horn SS, Meuwissen THE, Moghadam H, Hillestad B, Sonesson AK. Accuracy of selection for omega-3 fatty acid content in Atlantic salmon fillets. Aquaculture. 2020;519:734767. https://doi.org/10.1016/j.aquaculture.2019.734767.

    CAS  Article  Google Scholar 

  87. 87.

    Gebreyesus G, Bovenhuis H, Lund MS, Poulsen NA, Sun D, Buitenhuis B. Reliability of genomic prediction for milk fatty acid composition by using a multi-population reference and incorporating GWAS results. Genet Sel Evol. 2019;51(1):1–14. https://doi.org/10.1186/s12711-019-0460-z.

    Article  Google Scholar 

  88. 88.

    Sonesson AK, Meuwissen TH. Testing strategies for genomic selection in aquaculture breeding programs. Genet Sel Evol. 2009;41(1):37. https://doi.org/10.1186/1297-9686-41-37.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  89. 89.

    Vallejo RL, Silva RMO, Evenhuis JP, Gao G, Liu S, Parsons JE, et al. Accurate genomic predictions for BCWD resistance in rainbow trout are achieved using low-density SNP panels: evidence that long-range LD is a major contributing factor. J Anim Breed Genet. 2018;135(4):263–74. https://doi.org/10.1111/jbg.12335.

    CAS  Article  Google Scholar 

  90. 90.

    Vallejo RL, Leeds TD, Gao G, Parsons JE, Martin KE, Evenhuis JP, et al. Genomic selection models double the accuracy of predicted breeding values for bacterial cold water disease resistance compared to a traditional pedigree-based model in rainbow trout aquaculture. Genet Sel Evol. 2017;49(1):1–13. https://doi.org/10.1186/s12711-017-0293-6.

    CAS  Article  Google Scholar 

  91. 91.

    Haffray P, Bugeon J, Rivard Q, Quittet B, Puyo S, Allamelou JM, et al. Genetic parameters of in-vivo prediction of carcass, head and fillet yields by internal ultrasound and 2D external imagery in large rainbow trout (Oncorhynchus mykiss). Aquaculture. 2013;410-411:236–44. https://doi.org/10.1016/j.aquaculture.2013.06.016.

    Article  Google Scholar 

  92. 92.

    Haffray P, Bugeon J, Pincent C, Chapuis H, Mazeiraud E, Rossignol MN, et al. Negative genetic correlations between production traits and head or bony tissues in large all-female rainbow trout (Oncorhynchus mykiss). Aquaculture. 2012;368–369:145–52. https://doi.org/10.1016/j.aquaculture.2012.09.023.

    Article  Google Scholar 

  93. 93.

    Douirin C, Haffray P, Vallet JL, Fauconneau B. Determination of the lipid content of rainbow trout (Oncorhynchus mykiss) filets with the torry fish fat meter®. Sci Aliments. 1998;18:527–35.

    CAS  Google Scholar 

  94. 94.

    Marty-Mahe P, Loisel P, Brossard D. Color image segmentation to detect defects on fresh ham. Sixth Int Conf Qual Control by Artif Vis. 2003;5132:45. https://doi.org/10.1117/12.514936.

    Article  Google Scholar 

  95. 95.

    Robertson AR. The CIE 1976 color-difference formulae. Color Res Appl. 1977;2(1):7–11. https://doi.org/10.1002/j.1520-6378.1977.tb00104.x.

    Article  Google Scholar 

  96. 96.

    Lin W. Principles of magnetic resonance imaging: a signal processing perspective, vol. 19: SPIE Optical Engineering Press; 2005. https://doi.org/10.1109/memb.2000.870245.

  97. 97.

    Collewet G, Bugeon J, Idier J, Quellec S, Quittet B, Cambert M, et al. Rapid quantification of muscle fat content and subcutaneous adipose tissue in fish using MRI. Food Chem. 2013;138(2-3):2008–15. https://doi.org/10.1016/j.foodchem.2012.09.131.

    CAS  Article  PubMed  Google Scholar 

  98. 98.

    Thomas F, Jamin E, Wietzerbin K, Guérin R, Lees M, Morvan E, et al. Determination of origin of Atlantic salmon (Salmo salar): the use of multiprobe and multielement isotopic analyses in combination with fatty acid composition to assess wild or farmed origin. J Agric Food Chem. 2008;56(3):989–97. https://doi.org/10.1021/jf072370d,.

    CAS  Article  PubMed  Google Scholar 

  99. 99.

    Palti Y, Gao G, Liu S, Kent MP, Lien S, Miller MR, et al. The development an-d characterization of a 57K single nucleotide polymorphism array for rainbow trout. Mol Ecol Resour. 2015;15(3):662–72. https://doi.org/10.1111/1755-0998.12337.

    CAS  Article  PubMed  Google Scholar 

  100. 100.

    D’Ambrosio J, Phocas F, Haffray P, Bestin A, Brard-Fudulea S, Poncet C, et al. Genome-wide estimates of genetic diversity, inbreeding and effective size of experimental and commercial rainbow trout lines undergoing selective breeding. Genet Sel Evol. 2019;51(1):1–15. https://doi.org/10.1186/s12711-019-0468-4.

    Article  Google Scholar 

  101. 101.

    Sargolzaei M, Chesnais J, Schenkel FS. FImpute Users’ Guide, vol. 17; 2014.

    Google Scholar 

  102. 102.

    Misztal I, Tsuruta S, Strabel T, Auvray B, Druet T, Lee DH. BLUPF90 and related programs (BGF90), vol. 33. Montpellier: Commun. No. 28–07. 7th World Congr. Genet. Appl. to Livestick Prod; 2002. p. 2001–2.

    Google Scholar 

  103. 103.

    Habier D, Fernando RL, Kizilkaya K, Garrick DJ. Extension of the bayesian alphabet for genomic selection. BMC Bioinformatics. 2011;12(1):186. https://doi.org/10.1186/1471-2105-12-186.

    Article  PubMed  PubMed Central  Google Scholar 

  104. 104.

    Boerner V, Tier B. BESSiE: a software for linear model BLUP and Bayesian MCMC analysis of large-scale genomic data. Genet Sel Evol. 2016;48(1):63. https://doi.org/10.1186/s12711-016-0241-x.

    Article  PubMed  PubMed Central  Google Scholar 

  105. 105.

    Kass RE, Raftery AE. Bayes factors. J Am Stat Assoc. 1995;90(430):773–95. https://doi.org/10.1080/01621459.1995.10476572.

    Article  Google Scholar 

  106. 106.

    Michenet A, Barbat M, Saintilan R, Venot E, Phocas F. Detection of quantitative trait loci for maternal traits using high-density genotypes of blonde d’Aquitaine beef cattle. BMC Genet. 2016;17(1):88. https://doi.org/10.1186/s12863-016-0397-y.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  107. 107.

    Lallias D, Bernard M, Ciobotaru C, Dechamp N, Labbé L, Goardon L, et al. Sources of variation of DNA methylation in rainbow trout: combined effects of temperature and genetic background. Epigenetics. 2020. https://doi.org/10.1080/15592294.2020.1834924.

  108. 108.

    VanRaden PM. Efficient methods to compute genomic predictions. J Dairy Sci. 2008;91(11):4414–23. https://doi.org/10.3168/jds.2007-0980.

    CAS  Article  PubMed  Google Scholar 

  109. 109.

    Kuhn M, Johnson K, Kuhn M, Johnson K. Over-Fitting and Model Tuning. Appl. Predict. Model. 2013:61–92. https://doi.org/10.1007/978-1-4614-6849-3_4.

Download references

Acknowledgements

We gratefully acknowledge Jérôme Doerflinger, Thierry Espirat, and Mathilde Horat from Aqualande; Laurence Larroquet from INRAE; and Sophie Brard, Romain Morvezen, Pierre Patrice, and Maxime Reverchon from SYSAAF for their invaluable help in data collection. We thank Alain Moréac from Scanmat for his help for Raman data collection and Guylaine Collewet from INRAE for MRI data collection.

Funding

This study was supported by the European Maritime and Fisheries Fund and co-funded by FranceAgrimer (OmegaTruite project, n° P FEA470017FA1000008).

Author information

Affiliations

Authors

Contributions

P.H, M.D-N, G. C, J.B, V. N, and V.P. contributed to the conception and design of the study. C.B., E.P., N.D., J.B., F.E., G.C., and M.D-N collected the data. E. P, D. C, C.E-M, and F. E carried out Raman data collection and calibration. C.B. and F.E. organised the database. C.B., J.D’A., E.P., N.D., and F.E. performed data analysis. C.B., N.D., F.P., and M.D-N. collaborated on the design of the statistical methodology and interpretation of the results. F.P. and M.D-N. supervised the data analysis. C.B. wrote the first draft of the manuscript. M.D.-N, P.H., G.C., J.B., V.N. and F.P. contributed to manuscript revision. All authors read and approved the submitted version.

Corresponding author

Correspondence to Mathilde Dupont-Nivet.

Ethics declarations

Ethics approval and consent to participate

This study was conducted in accordance with EU Directive 2010–63-EU on the protection of animals used for scientific purposes. The fish were reared according to normal husbandry practices in the breeding programme of Les Sources de l’Avance. They were not subjected to practices likely to cause pain, suffering, distress, or lasting harm equivalent to, or higher than, that caused by the introduction of a needle, in accordance with good veterinary practice. As such, the experiment did not require approval by an Ethics Committee, in accordance with Article 2.5 of the Directive.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Fatty acid composition of the two final diets (% total FAs).

Additional file 2.

Estimates of phenotypic correlations between proportions of fatty acids and production or quality traits.

Additional file 3.

Candidate genes from the NCBI Oncorhynchus mykiss Annotation Release 100 (GCF_002163495.1) that were located within the confidence or credibility intervals estimated using BCπ analysis.

Additional file 4.

Summary statistics for 11 production and quality traits in rainbow trout.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Blay, C., Haffray, P., D’Ambrosio, J. et al. Genetic architecture and genomic selection of fatty acid composition predicted by Raman spectroscopy in rainbow trout. BMC Genomics 22, 788 (2021). https://doi.org/10.1186/s12864-021-08062-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-021-08062-7

Keywords

  • Fish
  • Genomic selection
  • QTL
  • GWAS
  • Fatty acid
  • Raman
  • Genetic correlations
  • Accuracy