- Research article
- Open Access
Functional genomics reveals increases in cholesterol biosynthetic genes and highly unsaturated fatty acid biosynthesis after dietary substitution of fish oil with vegetable oils in Atlantic salmon (Salmo salar)
BMC Genomics volume 9, Article number: 299 (2008)
There is an increasing drive to replace fish oil (FO) in finfish aquaculture diets with vegetable oils (VO), driven by the short supply of FO derived from wild fish stocks. However, little is known of the consequences for fish health after such substitution. The effect of dietary VO on hepatic gene expression, lipid composition and growth was determined in Atlantic salmon (Salmo salar), using a combination of cDNA microarray, lipid, and biochemical analysis. FO was replaced with VO, added to diets as rapeseed (RO), soybean (SO) or linseed (LO) oils.
Dietary VO had no major effect on growth of the fish, but increased the whole fish protein contents and tended to decrease whole fish lipid content, thus increasing the protein:lipid ratio. Expression levels of genes of the highly unsaturated fatty acid (HUFA) and cholesterol biosynthetic pathways were increased in all vegetable oil diets as was SREBP2, a master transcriptional regulator of these pathways. Other genes whose expression was increased by feeding VO included those of NADPH generation, lipid transport, peroxisomal fatty acid oxidation, a marker of intracellular lipid accumulation, and protein and RNA processing. Consistent with these results, HUFA biosynthesis, hepatic β-oxidation activity and enzymic NADPH production were changed by VO, and there was a trend for increased hepatic lipid in LO and SO diets. Tissue cholesterol levels in VO fed fish were the same as animals fed FO, whereas fatty acid composition of the tissues largely reflected those of the diets and was marked by enrichment of 18 carbon fatty acids and reductions in 20 and 22 carbon HUFA.
This combined gene expression, compositional and metabolic study demonstrates that major lipid metabolic effects occur after replacing FO with VO in salmon diets. These effects are most likely mediated by SREBP2, which responds to reductions in dietary cholesterol. These changes are sufficient to maintain whole body cholesterol levels but not HUFA levels.
FO has been, and still remains, the principal lipid source for fish diets in intensive aquaculture. Because the wild fisheries, from which FO is obtained, are currently fished to the maximal sustainable levels there is pressure to utilise alternative lipid sources in aquaculture diets. Currently the most practical alternatives are VOs and recent studies have demonstrated that VO can be used to replace up to 75% of FO without significant effects on growth in Atlantic salmon . However the effects of VO-based diets on normal metabolism and physiology and ultimately fish health and welfare are not yet understood. A variety of potentially deleterious effects have been reported in VO-fed fish. These include cardiac lesions , liver histopathology , compromised immune function  and disruption of intestinal function [6–9]. There is good reason to suspect that feeding of VO may have a major impact on fish physiology since VO differ considerably in composition compared to FO. For example VOs are rich in shorter chain, C18 polyunsaturated fatty acids (PUFA) and devoid of n-3HUFA. Indeed, fatty acid compositions in fish, including salmon, fed VO are characterised by increased levels of C18 PUFA and decreased levels of n-3HUFA, which could compromise their nutritional value to the human consumer [2, 10].
Massively parallel gene expression profiling technologies such as cDNA or oligonucleotide microarrays are powerful tools for discovering genes which change their tissue or cellular expression levels in response to changed conditions and thence enable the physiological mechanisms underlying such changes to be elucidated. Until recently there have been few such resources for commercially important fish species. The recent development of high density cDNA microarrays for Atlantic salmon (GRASP; TRAITS-SGP; [11, 12], the most commercially valuable farmed fish species in Europe and the Americas, has opened the way for fundamental studies on diet-gene interactions and promises to greatly advance understanding of fish nutrition.
The aim of this study was to discover mechanisms for physiological adaptation to VO-based diets in fish. This was achieved by measuring the effects on hepatic gene expression using a high density cDNA microarray, and by complementary biochemical and compositional assays of Atlantic salmon smolts fed diets in which 100% of the FO was replaced with three VOs, rapeseed oil (RO), linseed oil (LO) or soybean oil (SO).
Growth and biometry
Slightly, but significantly, lower final weights were obtained with fish fed the LO diets compared to fish fed FO (Table 1). However, no effect of dietary oil was observed in the specific growth rate (SGR) obtained for individually pit-tagged fish. Dietary VOs had no effect on feed efficiency as measured by feed conversion ratio (FCR), hepato-somatic index (HSI) or viscero-somatic index (VSI) or condition factor (Table 1).
Proximate, lipid and fatty acid compositions
All the VO diets significantly increased the relative percentage of protein in the whole fish compared to fish fed the FO diet (Table 1). There was also a tendency for the VO diets to decrease whole body lipid, and increase the protein:lipid ratio although these parameters were only significant in fish fed LO (Table 1). Liver total lipid content was not significantly affected by dietary VOs, but there was a trend towards increased free fatty acids in fish fed LO and increased total neutral lipid and triacylglycerol in fish fed SO and LO (Table 2). Dietary VO had only a few significant effects on the lipid class composition of liver with the most obvious being fish fed LO and SO showing reduced proportions of phosphatidylethanolamine compared to fish fed FO (Table 2). Differences in phosphatidylcholine levels were also observed between RO and LO, although none of the VO diets showed significant differences compared to FO. Increases in total unknown neutral lipid were also apparent in RO and SO fed fish compared to FO. In general the fatty acid compositions of liver and flesh total lipid (Table 3) reflected the fatty acid compositions of the diets. Thus fish fed the FO diet were characterised by high levels of 16:0, 20:5n-3 (eicosapentaenoic acid; EPA) and 22:6n-3 (docosahexaenoic acid; DHA), with 20:4n-6 (arachidonic acid; ARA) as the major n-6 fatty acid. Feeding the VOs resulted in increased proportions of 18:1n-9 and total monoenes, 18:2n-6 and total n-6PUFA, and 18:3n-3, and decreased the proportions of 16:0 and total saturated fatty acids, ARA, EPA and DHA. Specifically, the major fatty acids in fish fed RO, SO and LO were 18:1n-9, 18:2n-6 and 18:3n-3, respectively. Fish fed the vegetable oils also showed increased proportions of the metabolites of these fatty acids, namely 20:3n-3 and 20:4n-3 in fish fed LO, 20:2n-6 and 20:3n-6 in fish fed SO, and 20:1n-9 in fish fed RO.
Liver enzyme activities
No effect of diet was detected on the activity of liver carnitine palmitoyl transferase 1 (CPT1) activity, although fatty acid β-oxidation was increased in LO and RO fed fish. HUFA synthesis was increased in liver microsomes from fish fed all the VOs compared to fish fed FO. The activity of malic enzyme (ME) was increased in fish fed the RO diet in comparison to fish fed FO (Fig. 1).
Microarray gene expression profiling
Comparison of each experimental diet to FO identified 428, 350 and 434 features from RO, LO and SO respectively which were significantly different (p ≤ 0.05). These lists from each VO were then compared and features which were significantly changed relative to FO in two or more VO diets (T-test, P ≤ 0.05) were retained. This resulted in a list of 132 features. Features were then eliminated which were changed in different directions in different VO diets, resulting in a list of 128 features, 88 of which could be tentatively identified according to BLAST results. These genes whose expression was changed, and whose identity could be indicated from BLAST similarity, in two or more VO diets compared to FO, are listed in Table 4.
Gene ontology (GO) terms for biological processes were available for 39 of the 88 features with significant BLAST similarities that were changed in VO diets. The biological process GO terms associated with these 39 genes were compared with the total GO list for the TRAITS/SGP, and tested for significant (p ≤ 0.05) overrepresentation as implemented in GeneSpring GX 7.3.1 (Agilent). Only GO categories represented by more than one gene from the final gene list were considered significant and are reported in Table 5. Overrepresented GO categories included lipid metabolism, sterol biosynthesis, isoprenoid biosynthesis, fatty acid desaturation, mRNA processing, and translational elongation. It should be noted that the current gene ontology annotations (biological process) for the TRAITS/SGP array are only available for 3921 of the 16950 features. Notably amongst those genes in the final list which were overrepresented after GO analysis, were genes showing the greatest increase in expression. These included Δ5/6 fatty acyl desaturases which function in HUFA biosynthesis (enriched GO:6636, fatty acid desaturation) and isopentenyl-diphosphate Δisomerase (IPP), lathosterol oxidase, and squalene epoxidase, all genes of the cholesterol biosynthesis pathway (enriched GO:8299, isoprenoid biosynthesis; GO6695, cholesterol biosynthesis). Levels of other genes associated with lipid metabolism (enriched GO: 6629) or transport were also changed. These included increases in apoliporotein B (apoB) and decreases in an apoB-like gene in LO and SO. Effects were seen in levels of expression of genes involved in transcription, RNA processing and translation (enriched GO:6414).
Many of the remaining genes from the final list could also be grouped into distinct biological process categories, although they were not identified as such by GO analysis. Changes in carbohydrate metabolism genes included an increase in the expression of glycogen phosphorylase in all VO, and glucose 6-phosphate dehydrogenase in SO and RO. Several changes in expression levels of genes of protein degradation associated with the ubiquitin-proteosome pathway were observed, and also in the protease calpain1 which was upregulated more than five fold in all VO. Signal transduction genes were also changed and MARCKS-like 1 was increased in all VO treatments and protein phosphatase 1β decreased in all VO treatments. Genes for proteins with extracellular locations which changed in VO included the protease inhibitors cystatin, spint2 and tissue factor pathway inhibitor 2, as well as decreased expression of collagen1A3. In addition there were 44 features changed by VO for which no identity could be confirmed by BLAST (not shown).
Quantitative PCR (QPCR)
Results for QPCR are presented in Table 6. Three genes of the cholesterol biosynthesis pathway were upregulated in VO fed fish. One of these, IPP, represented by two microarray features, was selected for QPCR, another 7-dehydrocholesterol reductase was also on the array but was not significantly changed (p = 0.07). A further two genes, mevalonate kinase (MK) and SREBP2 were selected for QPCR because they were not present on the microarray, but were identified in BLAST searches of Atlantic salmon EST collections. MK was upregulated in all three VO, whilst SREBP2 mRNA was significantly increased in LO and RO, although not in SO (p = 0.057). HMG-CoA reductase (HMGR), another gene in the cholesterol biosynthesis pathway was not changed, either on the basis of microarray data or by QPCR. Δ6 fatty-acyl desaturase mRNA was also confirmed by QPCR to increase in all VO diets. A number of other lipid metabolic genes, selected on the basis of the microarray data were also assayed by QPCR. QPCR confirmed the microarray result that carnitine palmitoyl transferase A (CPT1A) mRNA was not changed, whilst QPCR of apolipoprotein B mRNA (APOB) demonstrated upregulation in LO and RO, whereas the array results indicated upregulation in SO and LO. Three other genes (catalase, PPARβ and short chain alcohol dehydrogenase) whose expression was not different between FO and VO by array analysis were also confirmed not to be changed after QPCR. However PPARγ which showed altered expression levels on micorarray analysis was not changed as measured by QPCR.
Although there were some relatively small differences, no great effects of diet on growth were observed in the present trial. Nevertheless, the effects of diet on tissue fatty acid composition are very clear with VO having highly significant effects. These effects are as expected and as reported in several previous trials on salmon [10, 13, 14].
As a means of more fully investigating the possible effects of VO substitution, liver was chosen to undertake a gene expression profiling experiment using the recently developed TRAITS/SGP microarray for Atlantic salmon. Liver integrates a large part of the nutritional uptake from the diet and has particular roles in distributing dietary lipid, and synthesising lipid de novo as well as participating in detoxification and excretion pathways. The use of high density cDNA microarrays is now a well established method, and we sought to use the TRAITS/SGP array to discover genes, and indicate biochemical and physiological pathways whose expression patterns might change during adaptation to VO-based diets. In this DNA microarray hybridisation experiment, application of ANOVA indicated that 784 features were changed at p ≤ 0.05 across all diet groups. However, application of either multiple testing correction (Benjamini and Hochberg) or a post hoc test (Tukeys multiple comparison) reduced this list of features 50 fold. Nevertheless, it was apparent that amongst the features from the list generated by ANOVA which passed filters for multiple testing correction or post hoc testing, were genes for Δ5/Δ6 fatty acyl desaturases. Given that these desaturase genes were known from prior experiments to be up-regulated in VO fed salmon  we filtered genes according to the criteria of having expression changed in the same direction in two or more of the VO dietary conditions after individual t-tests. In support of this approach it is notable that of the 132 genes which were scored as being changed in two or more VO diets, only 4 showed expression changes in opposite directions, far less than would be expected by chance alone. We tested the validity of the list of commonly changed candidates by a combination of compositional, biochemical and independent gene expression measurement by QPCR. The overall consistency of gene expression and physiological data demonstrates the validity of our approach to interpreting the microarray results presented here.
The most consistently and highly up-regulated genes that could be identified by homology or had previously been characterised and were clearly overrepresented after GO analysis were those of the Δ5 and Δ6 fatty acyl desaturases. Desaturases are encoded by several genes in salmon and each of them were represented on the array. Every one of these showed significant up-regulation in two or more of the VO diets. QPCR confirmed that desaturase mRNAs were increased and in addition the biosynthesis of HUFA was increased in liver microsomes from VO-fed fish. However, despite this clear increase in desaturase mRNA and HUFA biosynthesis, the amount of HUFA in salmon tissue was still considerably reduced when FO was replaced by VO in the diet. These results confirm previous reports describing the expression of salmon Δ5/Δ6 desaturase, where feeding of VO has been shown to increase desaturase expression and HUFA biosynthesis [15, 16].
In addition to the increase in fatty acyl desaturase expression, increases in mRNA for three enzymes (IPP, squalene epoxidase and lathosterol oxidase) of the cholesterol biosynthetic pathway were also evident from microarray data and from QPCR for IPP. To test the hypothesis that the cholesterol biosynthetic pathway was up-regulated following VO feeding the Atlantic salmon EST collections were searched for genes which are known to participate in cholesterol biosynthesis, but were not present on the TRAITS/SGP array. Two genes were identified and assayed by QPCR, MK and SREBP2. MK is known to be highly regulated by cholesterol level in mammalian liver and SREBP2 is a master transcriptional regulator of both cholesterol and HUFA biosynthesis . Messenger RNA for both of these genes was increased by a similar amount to that seen for IPP. The activity of SREBP2 is controlled by cellular cholesterol status, such that reductions in membrane cholesterol trigger proteolytic cleavage of membrane bound SREBP2 and lead to its translocation to the nucleus where it drives the transcription of the genes for cholesterol and HUFA biosynthesis . Of the 12 genes involved in lipid metabolism which were changed on the microarray after feeding VO, 5 are known to be regulated by SREBPs (Δ5 desaturase, Δ6 deasturase, IPP, lathosterol oxidase and squalene epoxidase) . Another gene in the cholesterol biosynthesis pathway, 7-dehydrocholesterol reductase, which by QPCR was upregulated in all VO, was on present the TRAITS-SGP array, but was not present in the final output gene list. One other gene critical for cholesterol biosynthesis, HMG-CoA reductase was present on the array but remained unchanged and QPCR analysis confirmed that VO had no effect on levels of HMG CoA reductase mRNA. This is surprising since in mammals HMG-CoA reductase is regulated directly by cholesterol levels, by hormonal stimuli, and transcriptionally by SREBPs, and therefore our results suggest some differences in the regulation of this gene between fish and mammals. Nevertheless, taken together the effects on genes of the cholesterol biosynthesis pathway strongly suggests that salmon fed VO respond to reduced cholesterol levels in the diet. Although we did not measure cholesterol contents of the diets, a recent report, which quantifies cholesterol and total sterol levels in the raw ingredients used to formulate the diets, confirms the deficiency of cholesterol in VO diets . Moreover, it is known that in vegetable oils the bulk of sterol is present as sitosterol, stigmasterol, campesterol or brassicasterol depending on source . These plant sterols are not taken up by mammals, and indeed inhibit cholesterol uptake  and, although not proven, are unlikely to be taken up by fish. Therefore one of the major responses of salmon liver during adaptation to VO is up-regulation of cholesterol biosynthesis to compensate for deficiencies in the diet. Importantly, and in contrast to the situation with HUFA, lower levels of dietary cholesterol were fully compensated for, since no differences were observed in cholesterol content, in liver or in flesh between fish fed FO or VOs. From the data presented here the compensation for reduced dietary cholesterol is due to up-regulation of hepatic cholesterol biosynthetic pathways, although other pathways such as decreased bile acid synthesis might also contribute. It should be noted, however, that this experiment consisted of a 12 week trial but in a longer term trial, after feeding VO for up to 22 months, both liver cholesterol and plasma LDL cholesterol were lower compared to FO .
Further to the up-regulation of cholesterol biosynthesis genes, apoB-100 which is a component of low density lipoprotein (LDL) required for the transport of cholesterol from liver to other tissues, also showed increases in expression level, both from microarray and QPCR measurement. Furthermore, from the microarray results angiopoietin-related protein 3 (ANGPTL3) was down-regulated in VO fed fish. ANGPTL3 is a liver derived factor which suppresses lipid uptake from LDL by tissues and promotes fatty acid release from adipocytes . Thus the reduction in ANGPTL3 observed in salmon fed VO is consistent with the increase in apoB-100 and taken together indicates increases in lipid transport to tissues from liver via LDL. Interestingly an apoB-like protein, which is related but clearly distinct from apoB-100 was down-regulated in VO fed fish, and a protein which is suggested to act as an "anti-freeze" in fish blood , but which resembles apoA, another lipid transport molecule, is up-regulated in all VO. Several other genes which showed changes in expression levels in VO fed fish have also been linked to lipid transport. For example cystatin, which was increased in all VO fed fish, is involved in tissue lipid uptake, suggested to be related to its role in controlling elastin and collagen breakdown in the extracellular matrix . Notably in this regard there were also decreases, in two microarray features identified as a collagen subtype, COL1A3, in all VO fish.
Overall these changes indicate a change in the function of the liver with respect to cholesterol synthesis and distribution, with an increase in cholesterol biosynthesis and an increase the transport processes which distribute cholesterol from the liver to other tissues. In animals fed normal cholesterol in the diet, the majority of cholesterol in the liver, from endogenous or exogenous pathways, is transported in the bile, as cholesterol or as bile acids, and thence to the gut where, together with phosphatidylcholine, it assists in the uptake of dietary lipid. In this regard, levels of hepatic phosphatidylcholine and its precursor phosphatidylserine were one of the few lipid compositional measures to vary between dietary treatments. These compositional changes might be related to changes in SO and RO diets of phospholipase D mRNA, which encodes an enzyme which releases choline from PC. Changes in bile transport, if occurring, may have consequences for lipid digestion, uptake and transport in the intestine which might explain the accumulation of lipid droplets in intestinal cells previously been observed in salmonids fed VO [6–8]. Interestingly, similar intestinal lipid accumulation is observed in rodents with deficiencies in bile transport .
Also notable from the microarray data was an increase in mRNA for glycogen phosphorylase (the first step of glycogen utilisation) and also in expression of glucose-6-phosphate dehydrogenase. These systems are involved in NADPH production  and suggests that the requirement for reducing power for cholesterol and HUFA biosynthesis is provided by breakdown of glycogen and up-regulation of NADPH generating systems. Another enzyme of importance for biosynthetic NADPH production, ME, was also shown to have increased activity in RO fed fish.
Some reports have indicated an increase in β-oxidation in VO compared to FO fed fish [27, 28]. From our microarray data we did not see increases in the expression of any genes linked to mitochondrial β-oxidation, indeed mRNAs for components of mitochondrial oxidative phosphorylation (pyrophosphates and phosphate carrier proteins) were decreased. Furthermore, QPCR and enzymic activity measurement confirmed that CPT1, an important rate limiting step in mitochondrial β-oxidation, was not significantly altered by VO feeding. However, it is known that peroxisomal β-oxidation predominates in salmon liver  and so the mRNA expression level of a critical enzyme of peroxisomal activity, acyl-CoA oxidase  was also measured. Acyl-CoA oxidase mRNA was increased in salmon fed the LO and RO diets, and this was consistent with the increased peroxisomal β-oxidation observed in fish fed LO and RO.
Of course this experiment represents a "snapshot" of gene expression 20 hours after final feeding and in liver only. At different prandial states and in different tissues the expression of genes for metabolic enzymes may be different. For example, shortly after feeding, when energy is in abundance, it is possible that the pentose phosphate pathway might be fuelled by gluconeogensis rather than glycogen breakdown, and these differences between FO and VO diets may be masked. With regard to fatty acid β-oxidation, sea bream exhibit differences in mitochondrial and peroxisomal activity depending on whether they are measured at 6 or 24 h after feeding .
The microarray data also indicated down-regulation of PPARγ, although, after QPCR no effect of VO feeding was apparent. PPARs are critical transcriptional regulators of lipid metabolism and energy homeostasis, and are encoded by three genes in mammals known as PPARα, PPARβ and PPARγ. However, in salmon at six PPAR genes have so far been discovered and they share considerable homology [33, 34]. Therefore it is possible that the microarray results reflect cross-hybridisation of different, but related transcripts, although a previous microarray study has also reported a decrease in PPARγ mRNA following VO feeding in salmon . In general cDNA microarrays are unlikely to be specific tools for measuring the expression of multiple closely related genes and this problem is exacerbated in Atlantic salmon which have polyploid genomes, often possessing more members within each gene family than the majority of diploid vertebrates. Similar issues been pointed out for the Δ5 and Δ6 fatty acyl desaturase genes which are 95% identical at the nucleotide level in salmon .
Increased lipid content, which might be an indication of health problems, in the livers of salmon fed VO has been observed previously [21, 36]. In our study, although there were no significant increases in hepatic lipid, a trend for increased lipid was supported by higher TAG and free fatty acid in VO fed fish and levels of mRNA for a marker for intracellular lipid accumulation, adipophilin , were increased in SO and LO, supporting the idea that at least some VO diets cause increases in hepatic lipid accumulation.
Finally from the microarray results, although not directly addressed by QPCR or biochemistry, were changes in the expression of genes involved in RNA and protein processing and turnover. These changes may be a reflection of changes in cholesterol biosynthetic pathways and in lipid transport, reflecting a requirement to export cholesterol to tissues.
In conclusion the present study demonstrates the utility of transcriptomic tools for studies on non-model commercially important fish species. Although, there were limitations in the statistical power of the design due largely to the inherent variability amongst recently domesticated farmed Atlantic salmon, by using appropriate filtering of microarray data, independent gene expression tests and biochemical and compositional data we were able to physiologically anchor a significant proportion of the gene expression results. This combined approach leads to the hypothesis that the major lipid metabolic effects of replacing FO with VO in salmon diets are mediated by SREBP2, a transcription factor which responds to reductions in membrane cholesterol level caused by reduced cholesterol in VO based diets. Furthermore, the application of the TRAITS/SGP microarray indicated several genes significantly altered by vegetable oils which suggest mechanisms to explain previously observed health impairment in VO fed salmon.
Diets and animals
Four diets (4 mm pellets) with the same basal composition but coated with four different oils were prepared at the Skretting Technology Centre, Stavanger, Norway. The diets were formulated to satisfy the nutritional requirements of salmonid fish (U.S. National Research Council 1993), and to contain 33% fat and 47% protein, consisting of fish meal (55%), corn gluten (10%), wheat (8.3%), oil (26.3%), mineral and vitamin mixes (0.2%) and Carophyl Pink® (0.06%). The oils used were FO (Anchovy oil, Skretting, Stavanger, Norway), RO (Denofa AS, Fredrikstad, Norway), LO (NOBA Vetveredeling B.V., Hamburg, Germany) and SO (Denofa AS, Fredrikstad, Norway). Fatty acid and proximate compositions of the experimental diets are given in Table 7. One thousand Atlantic salmon post-smolts, average weight 132 g, were distributed randomly into four 2 × 2 m tanks (250/tank) at the Nutreco Aquaculture Research Centre, Lerang Research Station, Stavanger, Norway. In order to measure growth in individual salmon, twenty fish per tank were PIT-tagged by implanting a micro-transponder into the peritoneal cavity. The fish were conditioned to the new environment for 3 weeks before feeding the experimental diets. During the conditioning period, the fish received commercial diet containing FO (Atlantic 4 mm; Skretting). Feed intake was monitored during this period, and the trial did not start until appetite was at least 0.8% body weight. After the conditioning period, the fish were fed the experimental diets to satiation according to usual procedures at Lerang Research Station, for a period of 16 weeks. Feed intake was monitored throughout the trial and waste feed was collected from the effluent water from each tank by a wire mesh collector and dried. Feed given, waste feed and the resulting net feed intake were registered daily. The tanks were supplied with sea water at constant temperature (7.9°C ± 0.1°C). The oxygen level varied between 8 and 12 ppm, with an average of 9.9 ppm. A photoperiod of 18 h light and 6 h darkness was applied.
At the start and end of the trial, all the fish in each tank were anaesthetised with metacain (50 mg/L) and individually weighed and measured. At the end of the trial, the 20 fish per dietary treatment were sampled for compositional, enzymatic and gene expression analyses, with three whole fish frozen immediately on dry ice and subsequently stored at -20°C for whole body compositional (proximate) analyses. The other sampled fish were eviscerated and twelve were used for biometric determinations (hepato- and viscero-somatic indices) and for lipid analyses. Livers were taken from each fish, pooled in four pools of three fish each, and frozen immediately in liquid nitrogen. Samples of liver were dissected from the remaining five fish from each dietary treatment and these samples were individually used for enzymatic, microarray and quantitative PCR (QPCR) analysis. A portion of each liver was taken for enzymatic analyses and frozen in liquid nitrogen. For gene expression analysis, samples of 0.5 g of fresh tissue were rapidly disrupted in 5 ml of TriReagent using an Ultra-Turrax homogeniser (Fisher Scientific, Loughborough, U.K.), and immediately frozen in liquid nitrogen for storage prior to RNA extraction.
Proximate analysis, consisting of moisture content and lipid and protein content, was undertaken as previously described . Lipid analyses were conducted on liver from twelve (four pools of three) fish from each diet and homogenised into pooled "pates". Total lipid was extracted from diets or 1 g portions of tissue pates and separation of lipid classes was performed by high-performance thin-layer chromatography (HPTLC) and fatty acid analysis performed as described previously .
Enzymatic and metabolic assays
All enzymatic and metabolic assays were carried out on the same five fish randomly selected from each dietary treatment. HUFA biosynthesis was determined in liver microsomes as described by Leaver et al., 2006. Assays of malic enzyme (ME), carnitine acyltransferase I (CPT1) and peroxisomal fatty acid β-oxidation in tissue homogenates were carried out on liver homogenates as previously described [14, 38].
Microarray gene expression profiling
Liver was selected as the tissue on which to focus gene expression studies because of its central role in the integration and control of nutritional physiology and energy homeostasis. Total RNA was extracted from the same five fish from each dietary treatment as used for enzymatic determinations. RNA was prepared from TriReagent homogenates as described in the manufacturers protocol (Sigma). Two separate RNA extracts were prepared from each sample, one for microarray probe preparation and the second for QPCR. Each biological replicate was co-hybridised in a two dye experiment with a single pooled reference sample. The pooled reference sample comprised equal amounts of RNA from each of the 20 biological samples. RNA was reverse transcribed and labelled with either Cy3 (pooled reference) or Cy5 (FO, RO, SO and LO fed fish) fluors using the FAIRPLAY II cDNA indirect labelling kit (Stratagene) according to the manufacturer's instructions. Briefly 20 μg total RNA was reverse transcribed after being primed with oligo dT, which incorporated aminoallyl-dUTP into the newly synthesised cDNA strand. The RNA template was then hydrolysed using 1 M NaOH for 15 min and neutralised with 1 M HCl. The cDNA was precipitated overnight. cDNA pellets were washed in 80% ethanol and air dried before being resuspended in 5 μL 2× coupling buffer (Stratagene FAIRPLAYII reagent). Once the cDNA had fully dissolved (after at least 30 min) 5 μL of either Cy3 or Cy5 dye was added and the samples incubated in the dark for 30 min. The Cy3 and Cy5 dyes (GE Healthcare) were dissolved in DMSO prior to being added to the coupling buffer. To remove unincorporated dye, the labelled cDNA (total volume 10 μL) was passed through a SpinEX column (Qiagen). Dye incorporation was assessed by separating 0.5 μL of the sample on a mini-agarose gel and visualising fluorescent products on a fluorescence scanner (Typhoon Trio, GE Healthcare). No microarray pre-hybridisation step was required. For hybridisation the remainder of each labelled cDNA (7–9 μL; 16–30 pmol each dye) was added to 85 μL hybridisation buffer (UltraHyb, Ambion), 10 μL poly(A)80 (10 mg mL-1; Sigma) and 5 μL ultrapure BSA (10 mg mL-1; Ambion). The hybridisation mixture was heated to 95°C for 3 min, then cooled to 60°C before being applied to the microarray (TRAITS/SGP, see section 2.7). Hybridisations were performed on a Gene TAC Hyb Station (Genomic Solutions) for 16 h at 45°C. Slides were then automatically washed with 2×SSC, 0.5% SDS for 10 min at 60°C; 0.2×SSC, 0.5% SDS for 10 min at 42°C; and finally 0.2 × SSC 10 min at 42°C. Following manual rinsing in isopropanol and drying by centrifugation, hybridised slides were scanned at 10 μm resolution using a Perkin Elmer ScanArray Express HT scanner. BlueFuse software (BlueGnome) was then used to identify and quantify features. Following manual spot editing to remove obvious artefactual features and fusion of duplicate spot data (BlueFuse proprietary algorithm), the resulting intensity values and quality annotations were exported into the GeneSpring GX version 7.3.1 (Agilent Technologies) analysis platform. Data transformation, normalisation and quality filtering were as follows: 1) all raw spot intensity values less than 0.01 were set to 0.01 to facilitate comparisons of transformed data; 2) a 'per spot per chip' intensity dependent (Lowess) normalisation was undertaken using software defaults (20% smoothing/cutoff 10); 3) data were filtered using a BlueFuse spot confidence value ≥ 0.2 in ≥ 8 slides and BlueFuse spot quality of ≥ 0.5 in ≥ 8 slides to filter out features which were below reasonable quality thresholds Bluefuse uses Bayesian statistical models of the signal and noise distributions on the slide to estimate the confidence that should be placed in each result. Similarly Bluefuse integrates various spot quality parameters (eg circularity, uniformity etc) to provide a spot quality assessment. This gave a final list of 11024 genes which were eligible for statistical analysis and application of ANOVA indicated that 9964 features had sufficient data for testing across all treatments. Experimental annotations complied fully with MIAME guidelines  and data has been submitted to the ArrayExpress databank (accession E-TABM-478; ).
Atlantic salmon ESTs were tentatively identified using BLASTX to search the Genbank non-redundant protein sequence database (all non-redundant GenBank CDS translations + RefSeq Proteins + PDB + SwissProt + PIR + PRF databases) applying a cut-off value of ≤ e-10. ESTs which were not assigned using this procedure were then compared to the Atlantic salmon DFCI Atlantic salmon Gene Index http://compbio.dfci.harvard.edu/tgi/tgipage.html comprising ESTs and ESTs assembled into Tentative Consensus (TC) sequences using BLASTN. TCs or ESTs showing 99.9% identity over a continuous overlap of >100 nucleotides were then used to search the Genbank non-redundant protein sequence database applying a cut-off value of ≤ e-10.
Quantitative PCR (QPCR)
QPCR was performed using an RNA extract from the same five animals that had been taken for enzymic and microarray analysis. Complementary DNA (cDNA) was synthesized from total salmon hepatic RNA using a commercially available kit according the manufacturers instructions (Reverse-iT Max 1st strand synthesis kit, ABgene). Briefly each reaction of 20 μL contained 1.5 μg of total RNA, 300 ng of random hexamers and 125 ng of anchored oligo-dT. Following cDNA synthesis at 42°C for 1 hour reactions were stopped by heating at 75°C for 10 min and cDNA diluted to 500 μL total volume with water. Real-time PCR was performed using a Quantica machine (Techne). Quantitative PCR analysis for each gene was performed in triplicate in a total volume of 20 μL containing 5 μL cDNA (equivalent to 15 ng of input RNA), 100 nM of each primer and 10 μL of Absolute QPCR SYBR Green Mix 2× (ABgene). For each target gene, forward and reverse primers were chosen from the available EST sequences by using Primer3 software. Target genes were selected by consideration of microarray data and by BLAST comparison of Atlantic salmon ESTs deposited in Genbank with sequences of mammalian genes of interest. Primer pairs and EST information are provided in Table 8. Thermal cycling was initiated with incubation at 95°C for 15 min in order to activate the Thermo-Start® DNA Polymerase present in the mix. After this initial step, forty-five cycles of PCR were performed. Each PCR cycle consisted of heating for 15 s at 95°C for denaturing, and then for 15 s at 60°C and 30 s at 72°C for annealing and extension. Cycle threshold (CT) values corresponded to the number of cycles at which the fluorescence emission monitored in real time exceeded the threshold limit.
Melting curve analysis was performed to indicate the production of a single product in these reactions. Standard curves were established for each gene by plotting CT values against the log10 of five different dilutions (in triplicate) of cDNA sample solutions. In addition a subset of randomly selected samples generated by each primer pair was analysed by agarose gel electrophoresis and sequencing to confirm the identity of amplicons. Real-time efficiency was determined for each gene from the slopes given by Quantsoft software, applying the equation E = 10(-1/slope). The calculated relative expression ratio of each gene was based on the PCR efficiency (E) and CT of sample compared with control, and expressed in comparison to the reference genes, glyceraldehyde phosphate dehydrogenase (GAPDH), β-actin, and unidentified EST Acc No. AJ425111 (REST© software, ). Statistically significant differences in gene expression between the control (FO) and samples (LO, SO and RO) were evaluated in group means by randomization tests  using REST© software. Five thousand random allocations were applied and differences were considered to be significant at P ≤ 0.05.
Materials for enzyme assays: [1-14C] Palmitoyl CoA (50–55 mCi/mmol) and [methyl-3H} L-carnitine (60–86 Ci/mmol) were obtained from GE Healthcare. BHT, Carnitine, coenzyme A, DTT, FAF-BSA, glutathione, HEPES, KCN, malate, N-acetylcysteine, NADH, NADP, nicotinamide, palmitoyl-CoA, perchloric acid, silver nitrate and TriReagent were obtained from Sigma.
Materials for lipid analyses: HPTLC (10 cm × 10 cm × 0.15 mm) and TLC (20 cm × 20 cm × 0.25 mm) plates, precoated with silica gel 60 (without fluorescent indicator) were obtained from Merck. All solvents were HPLC grade and were obtained from Fisher Scientific.
Unless otherwise stated, all data are presented as means ± SD (n value as stated). The effects of dietary treatment on biometry, composition and enzyme activity were analysed by one-way analysis of variance (ANOVA) followed, where appropriate, by Tukey's comparison test. Percentage data and data which were identified as non-homogeneous (Bartlett's test) were subjected to arcsine transformation before analysis. Differences were regarded as significant when P ≤ 0.05. Statistical treatment of microarray and QPCR data is described within the specific Methods and Results sections.
The probability that a particular biological process GO term was enriched in the input GO list (All GO annotations in the TRAITS-GS array) compared to the output GO list (all GO annotations in the experimentally altered gene list) was calculated using a hypergeometric distribution model as implemented by GeneSpring GX 7.3.1. The GO annotation of the TRAITS-GS microarray has been described by Taggart et al., 2008 
Sargent JR, Tacon AGJ: Development of farmed fish: a nutritionally necessary alternative to meat. Proceedings of the Nutrition Society. 1999, 58: 377-383.
Torstensen BE, Bell JG, Rosenlund G, Henderson RJ, Graff IE, Tocher DR, Lie O, Sargent JR: Tailoring of Atlantic salmon (Salmo salar L.) flesh lipid composition and sensory quality by replacing fish oil with a vegetable oil blend. Journal of Agricultural and Food Chemistry. 2005, 53: 10166-10178. 10.1021/jf051308i.
Bell JG, McVicar AH, Park MT, Sargent JR: High Dietary Linoleic-Acid Affects the Fatty-Acid Compositions of Individual Phospholipids from Tissues of Atlantic Salmon (Salmo Salar) - Association with Stress Susceptibility and Cardiac Lesion. Journal of Nutrition. 1991, 121: 1163-1172.
Bell JG, Tocher DR, Farndale BM, McVicar AH, Sargent JR: Effects of essential fatty acid-deficient diets on growth, mortality, tissue histopathology and fatty acid compositions in juvenile turbot (Scophthalmus maximus). Fish Physiology and Biochemistry. 1999, 20: 263-277. 10.1023/A:1007743532618.
Mourente G, Good JE, Bell JG: Partial substitution of fish oil with rapeseed, linseed and olive oils in diets for European sea bass (Dicentrarchus labrax L.): effects on flesh fatty acid composition, plasma prostaglandins E-2 and F-2 alpha, immune function and effectiveness of a fish oil finishing diet. Aquaculture Nutrition. 2005, 11: 25-40. 10.1111/j.1365-2095.2004.00320.x.
Olsen RE, Myklebust R, Ringo E, Mayhew TM: The influences of dietary linseed oil and saturated fatty acids on caecal enterocytes in Arctic char (Salvelinus alpinus L.): a quantitative ultrastructural study. Fish Physiology and Biochemistry. 2000, 22: 207-216. 10.1023/A:1007879127182.
Olsen RE, Dragnes BT, Myklebust R, Ringo E: Effect of soybean oil and soybean lecithin on intestinal lipid composition and lipid droplet accumulation of rainbow trout, Oncorhynchus mykiss Walbaum. Fish Physiology and Biochemistry. 2003, 29: 181-192. 10.1023/B:FISH.0000045708.67760.43.
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: 253-271. 10.1016/S0044-8486(01)00852-3.
Caballero MJ, Izquierdo MS, Kjorsvik E, Montero D, Socorro J, Fernandez AJ, Rosenlund G: Morphological aspects of intestinal cells from gilthead seabream (Sparus aurata) fed diets containing different lipid sources. Aquaculture. 2003, 225: 325-340. 10.1016/S0044-8486(03)00299-0.
Bell JG, Tocher DR, Henderson RJ, Dick JR, Crampton VO: Altered fatty acid compositions in Atlantic salmon (Salmo salar) fed diets containing linseed and rapeseed oils can be partially restored by a subsequent fish oil finishing diet. Journal of Nutrition. 2003, 133: 2793-2801.
von Schalburg KR, Rise ML, Cooper GA, Brown GD, Gibbs AR, Nelson CC, Davidson WS, Koop BF: Fish and chips: Various methodologies demonstrate utility of a 16,006-gene salmonid microarray. BMC Genomics. 2005, 6: 126-10.1186/1471-2164-6-126.
Taggart JB, Bron JE, Martin SAM, Seear P, Hoyheim B, Talbot R, Villeneuve L, Sweeney GE, Houlihan DF, Secombes CJ, Tocher DR, Teale AJ: A description of the origin, designs and performance of the TRAITS/SGP Atlantic salmon (Salmo salar) cDNA microarray. Journal of Fish Biology. 2008,
Rosenlund G, Obach A, Sandberg MG, Standal H, Tveit K: Effect of alternative lipid sources on long-term growth performance and quality of Atlantic salmon (Salmo salar L.). Aquaculture Research. 2001, 32: 323-328. 10.1046/j.1355-557x.2001.00025.x.
Torstensen BE, Lie O, Froyland L: Lipid metabolism and tissue composition in Atlantic salmon (Salmo salar L.) - Effects of capelin oil, palm oil, and oleic acid-enriched sunflower oil as dietary lipid sources. Lipids. 2000, 35: 653-664. 10.1007/s11745-000-0570-6.
Zheng XZ, Torstensen BE, Tocher DR, Dick JR, Henderson RJ, Bell JG: Environmental and dietary influences on highly unsaturated fatty acid biosynthesis and expression of fatty acyl desaturase and elongase genes in liver of Atlantic salmon (Salmo salar). Biochimica et Biophysica Acta-Molecular and Cell Biology of Lipids. 2005, 1734: 13-24. 10.1016/j.bbalip.2005.01.006.
Tocher DR, Bell JG, McGhee F, Dick JR, Fonseca-Madrigal J: Effects of dietary lipid level and vegetable oil on fatty acid metabolism in Atlantic salmon (Salmo salar L.) over the whole production cycle. Fish Physiology and Biochemistry. 2003, 29: 193-209. 10.1023/B:FISH.0000045722.44186.ee.
Horton JD, Shah NA, Warrington JA, Anderson NN, Park SW, Brown MS, Goldstein JL: Combined analysis of oligonucleotide microarray data from transgenic and knockout mice identifies direct SREBP target genes. Proceedings of the National Academy of Sciences of the United States of America. 2003, 100: 12027-12032. 10.1073/pnas.1534923100.
Tocher DR, Bendiksen EA, Campbell PJ, Bell JG: The role of phospholipids in nutrition and metabolism of teleost fish. Aquaculture. 2008,
Phillips KM, Ruggio DM, Toivo JI, Swank MA, Simpkins AH: Free and esterified sterol composition of edible oils and fats. Journal of Food Composition and Analysis. 2002, 15: 123-142. 10.1006/jfca.2001.1044.
Ostlund RE: Phytosterols and cholesterol metabolism. Current Opinion in Lipidology. 2004, 15: 37-41. 10.1097/00041433-200402000-00008.
Jordal AEO, Lie O, Torstensen BE: Complete replacement of dietary fish oil with a vegetable oil blend affect liver lipid and plasma lipoprotein levels in Atlantic salmon (Salmo salar L.). Aquaculture Nutrition. 2007, 13: 114-130. 10.1111/j.1365-2095.2007.00455.x.
Kersten S: Regulation of lipid metabolism via angiopoietin-like proteins. Biochemical Society Transactions. 2005, 33: 1059-1062. 10.1042/BST20051059.
Deng GJ, Laursen RA: Isolation and characterization of an antifreeze protein from the longhorn sculpin, Myoxocephalus octodecimspinosis. Biochimica et Biophysica Acta-Protein Structure and Molecular Enzymology. 1998, 1388: 305-314. 10.1016/S0167-4838(98)00180-0.
Bengtsson E, To F, Hakansson K, Grubb A, Branen L, Nilsson J, Jovinge S: Lack of the cysteine protease inhibitor cystatin C promotes atherosclerosis in apolipoprotein E-deficient mice. Arteriosclerosis Thrombosis and Vascular Biology. 2005, 25: 2151-2156. 10.1161/01.ATV.0000179600.34086.7d.
Voshol PJ, Minich DM, Havinga R, Elferink RPJO, Verkade HJJ, Groen AK, Kuipers F: Postprandial chylomicron formation and fat absorption in multidrug resistance gene 2 P-glycoprotein-deficient mice. Gastroenterology. 2000, 118: 173-182. 10.1016/S0016-5085(00)70426-4.
Pollak N, Dolle C, Ziegler M: The power to reduce: pyridine nucleotides - small molecules with a multitude of functions. Biochemical Journal. 2007, 402: 205-218. 10.1042/BJ20061638.
Stubhaug I, Lie O, Torstensen BE: Fatty acid productive value and beta-oxidation capacity in Atlantic salmon (Salmo salar L.) fed on different lipid sources along the whole growth period. Aquaculture Nutrition. 2007, 13: 145-155. 10.1111/j.1365-2095.2007.00462.x.
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. Biochimica et Biophysica Acta-Molecular and Cell Biology of Lipids. 2005, 1734: 277-288. 10.1016/j.bbalip.2005.04.003.
Froyland L, Lie O, Berge RK: Mitochondrial and peroxisomal beta-oxidation capacities in various tissues from Atlantic salmon Salmo salar. Aquaculture Nutrition. 2000, 6: 85-89. 10.1046/j.1365-2095.2000.00130.x.
Morais S, Knoll-Gellida A, Andre M, Barthe C, Babin PJ: Conserved expression of alternative splicing variants of peroxisomal acyl- CoA oxidase 1 in vertebrates and developmental and nutritional regulation in fish. Physiological Genomics. 2007, 28: 239-252. 10.1152/physiolgenomics.00136.2006.
Diez A, Menoyo D, Perez-Benavente S, Calduch-Giner JA, de Celis SVR, Obach A, Favre-Krey L, Boukouvala E, Leaver MJ, Tocher DR, Perez-Sanchez J, Krey G, Bautista JM: Conjugated linoleic acid affects lipid composition, metabolism, and gene expression in Gilthead sea bream (Sparus aurata L). Journal of Nutrition. 2007, 137: 1363-1369.
Desvergne BA, Michalik L, Wahli W: Transcriptional regulation of metabolism. Physiological Reviews. 2006, 86: 465-514. 10.1152/physrev.00025.2005.
Leaver MJ, Boukouvala E, Antonopoulou E, Diez A, Favre-Krey L, Ezaz MT, Bautista JM, Tocher DR, Krey G: Three peroxisome proliferator-activated receptor isotypes from each of two species of marine fish. Endocrinology. 2005, 146: 3150-3162. 10.1210/en.2004-1638.
Leaver MJ, Ezaz MT, Fontagne S, Tocher DR, Boukouvala E, Krey G: Multiple peroxisome proliferator-activated receptor beta subtypes from Atlantic salmon (Salmo salar). Journal of Molecular Endocrinology. 2007, 38: 391-400. 10.1677/JME-06-0043.
Jordal AE, Torstensen BE, Tsoi S, Tocher DR, Lall SP, Douglas SE: Dietary rapeseed oil affects the expression of genes involved in hepatic lipid metabolism in Atlantic salmon (Salmo salar L.). J Nutr. 2005, 135: 2355-2361.
Ruyter B, Moya-Falcon C, Rosenlund G, Vegusdal A: Fat content and morphology of liver and intestine of Atlantic salmon (Salmo salar): Effects of temperature and dietary soybean oil. Aquaculture. 2006, 252: 441-452. 10.1016/j.aquaculture.2005.07.014.
Motomura W, Inoue M, Ohtake T, Takahashi N, Nagamine M, Tanno S, Kohgo Y, Okumura T: Up-regulation of ADRP in fatty liver in human and liver steatosis in mice fed with high fat diet. Biochemical and Biophysical Research Communications. 2006, 340: 1111-1118. 10.1016/j.bbrc.2005.12.121.
Leaver MJ, Tocher DR, Obach A, Jensen L, Henderson RJ, Porter AR, Krey G: Effect of dietary conjugated linoleic acid (CLA) on lipid composition, metabolism and gene expression in Atlantic salmon (Salmo salar) tissues. Comparative Biochemistry and Physiology A-Molecular & Integrative Physiology. 2006, 145: 258-267. 10.1016/j.cbpa.2006.06.034.
Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC, Gaasterland T, Glenisson P, Holstege FCP, Kim IF, Markowitz V, Matese JC, Parkinson H, Robinson A, Sarkans U, Schulze-Kremer S, Stewart J, Taylor R, Vilo J, Vingron M: Minimum information about a microarray experiment (MIAME) - toward standards for microarray data. Nature Genetics. 2001, 29: 365-371. 10.1038/ng1201-365.
ArrayExpress. 2008, [http://www.ebi.ac.uk/arrayexpress]
Primer3. 2008, [http://frodo.wi.mit.edu/cgi-bin/primer3/primer3_www.cgi]
wwwGene-Quantificationinfo. 2008, [http://www.gene-quantification.info/]
Pfaffl MW, Horgan GW, Dempfle L: Relative expression software tool (REST (c)) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Research. 2002, 30:
Martin SAM, Taggart JB, Seear P, Bron JE, Talbot R, Teale AJ, Sweeney GE, Hoyheim B, Houlihan DF, Tocher DR, Zou J, Secombes CJ: Interferon type I and type II responses in an Atlantic salmon (Salmo salar) SHK-1 cell line by the salmon TRAITS/SGP microarray. Physiological Genomics. 2007, 32: 33-44. 10.1152/physiolgenomics.00064.2007.
This majority of this work was supported by the European Union (project number: Q5RS-2000-30360) and by a Marie Curie Fellowship to LV (project number: EIF-024589) and also partly by the University of Stirling through provision of a faculty research development fund award and by BBSRC UK (project number: 98/EGA17674).
MJL conceived and co-ordinated the study, integrated the data, performed biochemical studies and drafted the manuscript. LANV performed microarray sample preparation and hybridisations and QPCR. AO participated in experimental design, formulated diets and performed analysis of growth parameters. LJ performed dietary trial and co-ordinated sampling procedures. JEB participated in the analysis of microarray data and other statistical analyses. DRT performed lipid analyses and fatty-acyl desaturation assays. JBT participated in the co-ordination of microarray studies and experimental design. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.