Effects of genotype and dietary fish oil replacement with vegetable oil on the intestinal transcriptome and proteome of Atlantic salmon (Salmo salar)
BMC Genomics volume 13, Article number: 448 (2012)
Expansion of aquaculture requires alternative feeds and breeding strategies to reduce dependency on fish oil (FO) and better utilization of dietary vegetable oil (VO). Despite the central role of intestine in maintaining body homeostasis and health, its molecular response to replacement of dietary FO by VO has been little investigated. This study employed transcriptomic and proteomic analyses to study effects of dietary VO in two family groups of Atlantic salmon selected for flesh lipid content, 'Lean' or 'Fat'.
Metabolism, particularly of lipid and energy, was the functional category most affected by diet. Important effects were also measured in ribosomal proteins and signalling. The long-chain polyunsaturated fatty acid (LC-PUFA) biosynthesis pathway, assessed by fatty acid composition and gene expression, was influenced by genotype. Intestinal tissue contents of docosahexaenoic acid were equivalent in Lean salmon fed either a FO or VO diet and expression of LC-PUFA biosynthesis genes was up-regulated in VO-fed fish in Fat salmon. Dietary VO increased lipogenesis in Lean fish, assessed by expression of FAS, while no effect was observed on β-oxidation although transcripts of the mitochondrial respiratory chain were down-regulated, suggesting less active energetic metabolism in fish fed VO. In contrast, dietary VO up-regulated genes and proteins involved in detoxification, antioxidant defence and apoptosis, which could be associated with higher levels of polycyclic aromatic hydrocarbons in this diet. Regarding genotype, the following pathways were identified as being differentially affected: proteasomal proteolysis, response to oxidative and cellular stress (xenobiotic and oxidant metabolism and heat shock proteins), apoptosis and structural proteins particularly associated with tissue contractile properties. Genotype effects were accentuated by dietary VO.
Intestinal metabolism was affected by diet and genotype. Lean fish may have higher responsiveness to low dietary n-3 LC-PUFA, up-regulating the biosynthetic pathway when fed dietary VO. As global aquaculture searches for alternative oils for feeds, this study alerts to the potential of VO introducing contaminants and demonstrates the detoxifying role of intestine. Finally, data indicate genotype-specific responses in the intestinal transcriptome and proteome to dietary VO, including possibly structural properties of the intestinal layer and defence against cellular stress, with Lean fish being more susceptible to diet-induced oxidative stress.
Fish are important components of the human diet, being highly nutritious and valued as the main source of n-3 long-chain polyunsaturated fatty acids (LC-PUFA). These essential fatty acids, mainly eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA), have well-known health-promoting properties, including protection against a range of cardiovascular and inflammatory diseases, and neurological disorders . With population growth and increasing awareness of the importance of fish consumption as part of a healthy diet, worldwide demand for seafood continues to grow. However, as traditional fisheries are largely in decline, aquaculture must meet this demand. Aquaculture is the fastest-growing food production sector with an average annual growth rate of 6.6%, accounting for 46% of total fish supply . In the European and American continents, aquaculture production is largely dominated by salmonid species, mainly Atlantic salmon, and feeds for such carnivorous species have traditionally relied on fishmeal (FM) and fish oil (FO) from wild stocks. Recent estimates indicated that 88.5% of global production of FO was used by the aquaculture sector, with salmonid culture taking the largest share (56%) . With ever increasing demands for aquafeeds and reduction in fisheries landings, the availability of FO and FM seriously limits the growth of aquaculture production and there is an urgent need to find more sustainable alternatives.
Vegetable oils (VO) can replace FO in salmon feeds without compromising fish growth or condition although, at high levels of replacement, tissue levels of n-3 LC-PUFA are significantly reduced [4, 5]. The effects of FO replacement by VO are becoming well characterized in the hepatic transcriptome of salmonids [6–10], and other species . However, studies on intestinal transcriptome are few and restricted to effects of replacement of FM by plant proteins, particularly soybean meal [12, 13], given its potential to cause enteritis. It is now clear that intestine in salmonids is not simply a site for reacylation and packaging of dietary lipids but it also has important roles in fatty acid metabolism, including LC-PUFA biosynthesis [14, 15]. Furthermore, dietary VO can induce major histological changes in fish enterocytes, originating mostly from supranuclear lipid droplet formation, possibly due to altered reacylation mechanisms and decreased phospholipid synthesis [16–18]. In some cases, these accumulations were large enough to be deemed pathological [19, 20]. A recent study investigating effects of dietary FO replacement by VO on intestinal transcriptome in Atlantic cod indicated potential effects on lipid absorption and transport and suggested morphological and structural changes to the intestinal muscle layer . Furthermore, both this and a previous study on Atlantic salmon showed significant effects on expression of genes involved in cell proliferation and apoptosis . Therefore, there is indication that intestine may be affected by changes in lipid components of feed formulations. Given its crucial roles in nutrient absorption, protection against the entry of pathogens, and immune function , further attention is warranted and impacts of FO replacement require investigation in intestine, particularly in salmon where important changes in diet formulation are already being applied. This study is a large-scale analysis of the effects of replacement of dietary FO by VO on the transcriptome and proteome of Atlantic salmon intestine. Furthermore, given recent interest in evaluating genetic selection as a feasible strategy, in conjunction with changes in commercial feed formulation, to meet worldwide demand for farmed fish without compromising animal welfare or nutritional value [9, 10, 24], two groups of Atlantic salmon families, Lean and Fat, were studied to examine the potential effects of genetic background. This experiment was performed in parallel with another microarray study looking at effects in the hepatic transcriptome, analysing samples from the same individuals [9, 10], enabling a global and comprehensive assessment of the physiological and molecular effects of FO replacement by VO in Atlantic salmon, including potential interactions with genotype.
Two way-ANOVA of the cDNA array dataset returned 1409, 1626 and 862 significant genes for the factors diet, genotype and diet × genotype interaction, respectively. Detailed analysis was restricted to the top 100 most significant features, which were categorised according to biological function, based on mammalian homolog genes. Metabolism, particularly of lipid and energy, was the functional category most affected by diet accounting for 39-41% of the top 100 annotated genes (Table 1), and showing highest diet × genotype interaction (Additional file 1). Diet also impacted translation (17%; mostly ribosomal proteins) and signalling (17%). In contrast, genotype affected less markedly metabolism (21%), whereas structural proteins (25%; mostly different types of collagen) and proteins involved in the regulation of transcription (21%) predominated (Table 2).
Gene Ontology (GO) enrichment analysis was performed on the complete significant lists, enabling identification of GO terms significantly enriched in the input entity list, in comparison to the whole array, providing clues as to which biological processes might be particularly altered in the experimental conditions being compared. It revealed no significant enrichment of GO terms in the genotype list, while 20 and 7 GO terms were significantly enriched in the diet and interaction lists, respectively. GO terms enriched in the diet list included structural constituents of ribosome, structural molecule activity, cytosolic ribosome, cytosol, ribosomal subunit, translation, cellular biosynthetic process, gene expression, macromolecule and biopolymer biosynthetic process and other related terms. This was explained by the large number of ribosomal proteins, components of both the 40S and 60S subunits, which were down-regulated by dietary VO (Table 1). In contrast, several Δ6-desaturase (Δ6fad) clones showing a diet × genotype interaction (see Additional file 1) caused a significant enrichment of the GO terms oxidoreductase activity, stearoyl-CoA 9-desaturase activity, unsaturated fatty acid biosynthetic activity/metabolic processes and very long chain fatty acid (VLCFA) biosynthetic activity/metabolic processes.
RT-qPCR analysis of gene expression
The expression of several genes significantly affected or related to processes affected by the two factors in the microarray analysis was determined by RT-qPCR (Table 3). For diet, a reasonably good match was found for Δ5 fatty acyl desaturase (Δ5fad), NADH dehydrogenase subunit 1 (NADH1), proliferation-associated 2G4b (PA2G4), 60S acidic ribosomal protein (RP60S), proliferating cell nuclear antigen (PCNA; also known as cyclin) and cytochrome P450 1A (CYP1A), particularly in the Fat group where fold-changes were generally more pronounced and significant. No change in expression of uncoupling protein 2 (UCP2) with diet was measured while, for myosin heavy chain (MYO) and methylenetetrahydrofolate dehydrogenase (NADP+ dependent) 1-like (MTHFD1), RT-qPCR indicated a change opposite to that suggested by microarray. Regarding genotype, a good match was obtained for CYP1A, proteasome subunit beta type-8 precursor (PSMB8) and alpha 2 type I collagen (COL1A2), while transgelin 2 (TAGLN) expression did not differ between family groups, and for ATP-binding cassette sub-family A member 1 (ABCA1) there was an inverse change in expression.
As well as validation above, RT-qPCR was used to further analyse genes/pathways identified by microarray and published data as potentially interesting, including lipid metabolism, xenobiotic and oxidative stress, and apoptosis. One was LC-PUFA biosynthesis, given that Δ5fad was significantly affected by diet in the microarray analysis, with a stronger response in Fat fish, whereas Δ6fad showed a significant diet × genotype interaction (supplementary file 1) confirmed by RT-qPCR. The Δ6fad transcript was only significantly up-regulated in Fat fish fed VO, compared to FO, and Lean fish showed higher levels of Δ6fad expression than Fat fish when fed FO, while the opposite trend was noted when fish were fed VO. Fatty acyl elongases (elovl’s) were also quantified and their expression (except elovl5a) broadly followed that of fads: significantly up-regulated when dietary VO replaced FO in the Fat group. Additionally, elovl5b and, particularly, elovl2 showed a trend for increased expression (high fold-change) in Lean fish, compared to Fat fish, when fed FO, while an opposite trend was observed when salmon were fed VO. Although genes involved in fatty acid synthesis and oxidation showed few significant differences, expression of fatty acid synthase (FAS) was up-regulated in fish fed VO, but only in Lean fish. The expression of peroxisome proliferator-activated receptors (PPARs), involved in the regulation of multiple lipid metabolism genes, was determined but only PPARγ showed any significant change, being up-regulated in the Fat group when dietary VO replaced FO. Of the xenobiotic and oxidative metabolism genes assayed, apart from CYP1A, only catalase (CAT) was affected by diet (up-regulated in fish fed VO) and only significantly in the Lean family. In contrast, metallothionein A (MT-A) showed higher expression in Fat fish, but only when fed VO, while a marginal down-regulation was observed when comparing VO and FO-fed fish in the Lean group. Of genes related to apoptosis, CASP3B was up-regulated by VO in Lean fish whereas a similar fold-change was marginally non-significant in the Fat fish.
Intestine fatty acid composition
The levels of most fatty acids in pyloric caeca were affected by diet, whereas genotype had no significant effect (Table 4). However, some fatty acids also showed a significant diet × genotype interaction, indicating that the effect of diet depended on the genetic background of the fish. For instance, interactions were observed for some LC-PUFA as a result of higher levels being found in the Lean group, compared to Fat, when fish were fed VO, while the reverse was observed when fed FO. Another unexpected result was that similar levels of DHA in FO- and VO-fed Lean fish meant that, in spite of substantial differences in Fat fish fed the two diets, the effect of diet on DHA was marginally non-significant. Similarly, levels of EPA and 22:5n-3 between FO- and VO-fed fish were noticeably closer in the Lean group.
Of the protein spots identified as being differentially expressed between diets or genotypes (Additional files 2 and 3), only 17 and 29 could be excised and, of these, only 9 and 20, respectively, returned reliable identifications by peptide fragment fingerprinting (Figure 1). Proteins significantly up-regulated by dietary VO are likely implicated in xenobiotic/drug metabolism (epoxide hydrolase 2, EPHX2), protection from oxidative stress (hemopexin-like protein, HPX and peroxiredoxin-1, PRDX1) and induction of apoptosis and inflammatory responses (galectin-2, LGALS2). Those proteins down-regulated by dietary VO included proteins responsible for protein folding and involved in signalling (2-peptidylprolyl isomerase A, PPIA), actin-based motility (myosin light chain, MYL) and DNA replication, repair or transcription (histone H2A, H2A) (Table 5). Proteins affected by genotype encompass a variety of pathways, of which only a few are related to metabolism, namely carbohydrate (alpha-enolase, ENO1; pancreatic alpha-amylase precursor, AMY2; triosephosphate isomerase 1b, TPI1), folate (dihydropteridine reductase, DHPR) or retinol (retinol-binding protein II, RBP2) metabolism (Table 6). Other proteins may have potentially multiple roles but can broadly be assigned roles in response to oxidative and cellular stress (HPX; PRDX; heat shock protein 70, HSP70), oxygen transport (alpha globin, HBA), signal transduction (calreticulin, CALR), transcription/RNA repair (heterogeneous nuclear ribonucleoprotein A0, HNRNPA0; histone H2A, H2A), apoptosis (CASP3), cellular transport, potentially also associated with apoptosis (voltage-dependent anion channel 2–2, VDAC2; annexin A4, ANX4), and proteolysis (proteasome beta 1 subunit, PSMB1). As with the microarray analysis, a few proteins with a more structural function and particularly associated with tissue contractile properties were affected by genotype, showing lower levels in Lean fish. These included calponin-1 (CNN1) and transgelin (TAGLN), the latter which was also found to be significantly affected by microarray, albeit up-regulated in Lean fish. Most proteins significantly affected by genotype showed lower levels of expression in the Lean group, with the exception of ENO1, HSP70, TPI1, H2A and HBA.
Dietary plant ingredients can induce chronic intestinal inflammatory conditions in salmonids that can ultimately result in carcinogenesis . This extreme reaction is rare and usually associated with soy protein at high levels . Dietary n-3 LC-PUFA have important anti-inflammatory and anti-carcinogenic effects in mammalian intestine . Therefore, use of feeds containing high percentages of plant proteins combined with replacement of FO by VO, as is now prevalent in the industry, requires studies on dietary effects on intestinal transcriptomes and proteomes. However, interpretation of the data was difficult as the effects on dietary treatments and/or family groups were subtle, as also observed in liver transcriptome , and is typical of this type of experiment [6, 8, 11, 27]. Partly as a consequence, validation of the microarray data gave variable results, from perfect match, to opposite changes in a few, although effects observed in the microarray, with fold changes as low as 1.2 (e.g., NADH1 and PA2G4) were validated by RT-qPCR. In view of the whole genome duplication event that occurred in salmonids , gene expression studies are often more challenging due to the presence of highly similar genes which may hybridize with cDNA probes presenting low specificity, further complicated if similar transcripts, corresponding to duplicated genes, are differentially regulated . Nonetheless, the presence of several features related to specific processes in both the transcriptomic and proteomic analysis gave supporting evidence to the pathways likely differentially affected by dietary oil and genetic background related to flesh adiposity.
Intestinal LC-PUFA biosynthesis capacity is differentially affected by diet and genotype
Considering whether genetic selection for fish families showing better adaptation to more sustainable feeds might be a viable approach to develop aquaculture, one outcome of this investigation was to establish if effects of diet on expression of LC-PUFA biosynthesis genes depended on genotype, as shown in the liver transcriptome of these fish . This was not seen in the hepatic transcriptome of European sea bass families showing different growth rates when fed a vegetable diet but, in this, case similar LC-PUFA profiles were also noted in both genotypes in response to the vegetable diet . In both salmon tissues, differences in n-3 LC-PUFA content between fish fed FO or VO were smaller in the Lean family group. This was due to higher levels of n-3 LC-PUFA in Lean salmon, compared to Fat, when fish were fed VO, but higher amounts in the Fat family group when fed FO. However in liver, up-regulation of LC-PUFA biosynthesis when fish were fed VO was much larger in the Lean family group, whereas in intestine the same individuals only showed significant up-regulation in Fat fish. This appears contradictory but can be explained by the differential tissue n-3 LC-PUFA contents. Although the difference was smaller in Lean fish compared to Fat fish in both tissues, in liver there was still a considerable difference in n-3 LC-PUFA levels between fish fed FO or VO, while in intestine levels were similar. PUFA have important activities on transcription factors, either as direct ligands or through effects on membrane composition , affecting transcription of many genes involved in lipid metabolism, including desaturases and elongases [30, 31]. In salmon, regulation of genes of LC-PUFA biosynthesis that are known to respond to dietary composition, i.e., Δ5fad, Δ6fad, elovl5b and elovl2 [32–34], appear to show high plasticity and are likely under feedback regulation by tissue n-3 LC-PUFA. Both studies suggest that the Lean family group may show an enhanced response to low dietary n-3 LC-PUFA, with greater up-regulation of biosynthesis when fed VO. In contrast to liver, this response was sufficient in intestine to maintain tissue n-3 LC-PUFA, particularly DHA, at similar levels to fish fed FO. Considering that differences in desaturase expression between the Fat and Lean fish were only significant when FO but not VO, was fed, suggests that the likely mechanism is through negative feedback by high levels of n-3 LC-PUFA rather than positive feedback from low levels of LC-PUFA and/or higher levels of shorter chain precursors.
Other dietary effects on lipid metabolism
Transcriptional regulation of desaturases and elongases by LC-PUFA may involve both PPARα and sterol regulatory element binding protein 1c (SREBP-1c) [30, 31]. In liver, expression of Δ5fad, Δ6fad, elovl2, PPARα, and possibly PPARβ, appeared co-ordinately regulated by diet depending on genotype, while PPARγ was not affected . In intestine, however, expression of PPARα and PPARβ was not affected by either diet or genotype, while PPARγ was up-regulated by dietary VO, significantly in Fat fish. This suggests that dietary regulation of lipid metabolism genes in fish intestine might differ to mammals, where PPARα showed differential expression in response to dietary EPA and DHA in murine intestine . Reasons for differential regulation of PPARs between salmon liver and intestine are unclear, but may be due to different patterns of tissue expression. In plaice and seabream, there was no nutritional regulation (fasting/feeding) of PPARs in the intestine, where PPARγ was the dominant isotype, in contrast to liver where PPARα was dominant . PPARγ in both mammals and fish is predominantly expressed in adipose tissue and promotes adipocyte differentiation and lipid storage [36–38]. In mammals, PPARγ activates the expression of genes characteristic of mature adipocytes and adipogenesis, including FAS [37, 38] and hence the expression of PPARγ, up-regulated in salmon fed VO, might be related to increased expression of FAS. However, increased PPARγ expression was only significant in Fat fish whereas FAS was significantly up-regulated only in Lean salmon. As fish PPARγ is functionally the most different of the three isotypes compared to mammalian PPARs, and is expressed more widely in fish tissues that in mammals, other mechanisms and functions may underlie the observed regulation .
In this study, the hypotriglyceridemic effect of LC-PUFA, well established in mammals [39, 40], was also observed in salmon intestine. Lipogenesis was down-regulated in FO-fed fish, as demonstrated by decreased FAS expression (RT-qPCR) and the presence of a transcript containing a beta-ketoacyl synthase domain, a component of FAS (microarray). The differences in FAS expression were not as marked as in liver  and were only significant in Lean fish but, together with the LC-PUFA biosynthesis data, demonstrate the active role of salmon intestine in lipid metabolism . However, despite up-regulation of lipogenesis by dietary VO, lipid accumulation in enterocytes was lower than in fish fed FO, contrary to previous reports of VO promoting lipid accumulation in enterocytes [16–18].
In contrast, the hypotriglyceridemic effect of LC-PUFA did not involve the typical increase in β-oxidation , reported in mice intestine [36, 40]. As in liver, no changes were observed in the expression of β-oxidation genes carnitine palmitoyltransferase I (CPT1) and acyl-CoA oxidase (ACO) . Nonetheless, effects of dietary lipid on energy metabolism were observed in intestine. In particular, UCP and transcripts involved in the mitochondrial electron transport chain, including components of cytochrome c oxidase, NADH1 and ubiquinol-cytochrome c reductase complexes, and the mitochondrial metabolite transporter SCaMC-2, were slightly down-regulated by dietary VO, possibly suggesting reduced energetic metabolism. EPA may act as a mitochondrial proliferator in both rat and salmon liver [41, 42], which might also explain this result.
Vegetable oils as potential sources of contaminants in alternative diets for aquaculture?
The use of FO in feeds has also been questioned in relation to levels of persistent organic pollutants, POPs , whereas VO are generally considered safer alternatives to dietary FO . However, results suggested that VO had greater impact on xenobiotic metabolism in salmon intestine than FO, with the transcriptome and proteome both showing up-regulation of transcripts and proteins involved in detoxification and protection from oxidative stress in fish fed VO. This was surprising, considering expression of these genes has been associated with FO supplementation and higher levels of organic contaminants and/or increased peroxidative susceptibility of LC-PUFA [22, 45, 46]. In particular, and linked to a detoxification role, up-regulation of CYP1A transcript and epoxide hydrolase 2 (EPHX2) protein was found in fish fed VO. Cytochrome P450 1A metabolizes many exogenous and endogenous molecules, including pollutants (e.g., polychlorinated biphenyls-PCBs, polycyclic aromatic hydrocarbons-PAHs, dioxins and dibenzofurans) as well as several metabolic products (e.g., steroids, bile acids). Its expression is sensitive to low levels of contaminants and therefore it is a commonly used marker [47, 48]. Lipid peroxidation and antioxidant enzymes are also biomarkers for environmental xenobiotic contamination in fish, as the catalytic actions of detoxifying enzymes like CYP1A produce high levels of reactive oxygen species (ROS) [49, 50]. We observed up-regulation of CAT and of a selenoprotein transcript, as well as of HPX and PRDX1 proteins in salmon fed VO. CAT and PRDX1 catalyze the decomposition of hydrogen peroxide , hemopexin prevents heme-mediated oxidative stress , and selenoproteins have diverse roles, including selenium transport and antioxidant defense . Expression of other genes involved in protection from oxidative stress like GST and SOD were not affected, but MT-A was down-regulated in the Lean family group fed VO. MT-A, besides having an antioxidant role, also protects against heavy metal toxicity and maintains physiological zinc homeostasis  and hence could be responding to contaminants more abundant in FO (e.g., heavy metals) .
Previous analysis of contaminants in comparable feeds using the same standard FO and a similar VO blend, showed that levels of several POPs, including organochlorine pesticides, and heavy metals were substantially lower, but polycyclic aromatic hydrocarbons (PAH) were 10-fold increased in the VO diet compared to the FO diet . PAH derive from incomplete combustion of organic compounds and can be found in high concentrations in fats, including VO, through multiple routes of contamination, before, during or after oil processing . They are metabolized by both CYP1A1 and EPHX2, among other enzymes, and CYP1A1 is induced by PAH in mammals . Therefore, higher PAH levels in the VO diet might explain, at least partly, the results obtained in the present study although, unlike POPs, PAHs are not persistent and are readily eliminated from fish tissues .
High doses of PAH result in significant intestinal hyperplasia in fish, with an increase in cell proliferation and faster epithelial turnover . Previous studies on intestinal gene expression in fish indicated a reduction in cell proliferation or differentiation associated with dietary FO replacement by VO [21, 22], possibly due to lower levels of membrane LC-PUFA and reduced oxidative stress . In the present study, no major impact on cell proliferation was apparent in the intestinal transcriptome or proteome data. Two transcripts related to cell proliferation, PA2G4 and cyclin G1, were slightly down-regulated in fish fed VO, but in mammals these have opposing effects [60, 61] and, furthermore, two mammalian PA2G4 isoforms have been shown to have opposite effects in cellular proliferation  and hence results are inconclusive.
Previously, expression of caspases, effectors of controlled cell death or apoptosis, was affected by replacement of dietary FO by VO in fish [21, 22]. Apoptosis is particularly important in organs with high rates of cellular turnover such as intestine but, in addition to maintaining normal gut function, apoptosis may be affected by pathological or toxic conditions, including those induced by environmental chemical contaminants [63, 64]. In the present study, expression of CASP3B was up-regulated in salmon fed VO, particularly in the Lean family group and a similar, non-significant trend was observed for CASP6A/B. As ROS are important signalling molecules in apoptotic processes, these results could be linked to a cytotoxic effect causing increased oxidative stress in VO . Relevant to the above was the up-regulation of galectin 2 (LGALS2) in the proteome of salmon fed VO. Galectins are pleiotropic regulators of immune functions and are up-regulated by injury and infectious conditions, have well-recognized modulatory roles in mammalian intestinal inflammatory diseases, and their mode of action involves induction of apoptosis [65, 66]. The lack of major effects on cell proliferation and only slight up-regulation of CASP3 and LGALS2 suggests that any contaminant doses experienced by the fish were unlikely to have caused any serious morphophysiological damage in the intestine. As similar trends were not seen in the hepatic transcriptome of these individuals , this may suggest intestine can potentially metabolize and detoxify xenobiotics present in the diet. Furthermore, there were no growth or general performance issues with these fish . Therefore, the data do not imply abnormal gastrointestinal functions or effects on final product quality.
Effect of genotype in intestinal transcriptome and proteome
Contrary to diet, genotype did not have a major impact on metabolism genes, apart from transcripts related to the proteasomal degradation pathway including a strong down-regulation of PSMB8 in Lean fish, particularly fed VO. This gene has been recently found to have a molecular evolution history that suggests a very strong selective pressure for its functional dimorphism to be maintained in vertebrates . Two different alleles, A-type and F-type, can be found in basal vertebrate species, including Atlantic salmon. The PSMB8F lineage was lost in common ancestors of higher teleosts and tetrapods but was then independently revived de novo through the appearance of F-type alleles within the PSMB8A lineage. In this study we did not find evidence of significant differences between families groups for the A-type allele in the transcriptomic analysis as the probe showing significant variation between families in the microarray corresponded to the PSMB8F allele. Hence, it was also the F-type transcript that was validated by RT-qPCR, using type-F specific primers. However, further studies would be required to confirm this and to assess the functional significance of this result. On the other hand, expression of PSMB1 was down-regulated in the intestine proteome of Lean fish. Proteolysis through this pathway is essential for many cellular processes, including the cell cycle, signalling, cellular defence and responses to oxidative stress . Therefore, this response might be related to defence against cellular stress, as another difference between the two family groups was related to xenobiotic and oxidant metabolism. Apart from lower expression of a CYP1A transcript in Lean fish, two proteins with antioxidant roles, HPX and PRDX, were down-regulated in the proteome of Lean compared to Fat fish. Alpha globin, or haemoglobin alpha (HBA), a major component of blood and potent mediator of oxidative stress, can have both protective and damaging effects depending on complex interactions in H2O2-rich environments . However, given its opposite regulation to HPX, whose main role is to scavenge heme and protect from its toxic effects , up-regulation of HBA in Lean fish may indicate heme-mediated oxidative stress.
The apoptotic pathway may be differentially affected by genotype, with down-regulation of CASP3, VDAC2 and ANXA4 in the Lean family group, the latter two transport proteins having well-recognized roles in apoptosis [70, 71]. In contrast, heat shock proteins that protect against environmental stresses were increased in the intestine transcriptome and proteome of Lean salmon. This response could be associated with the observed changes in the ubiquitin-proteasome degradation system, as the systems have been functionally coupled in mammals. Thus, moderate exposure to a heat shock can cause a transient increase in intracellular proteolysis by the ubiquitin-proteasome pathway, followed by a phase of slower or even inhibited protein degradation [72, 73]. Furthermore, Pirkkala et al.  demonstrated transcriptional induction of heat shock genes when proteasome activity was down-regulated. However, judging by the fold-changes, these effects are only relevant when fish were fed VO, and hence could be more related to dietary changes. Collectively, the data may indicate higher sensitivity of Lean fish to environmental or endogenous stresses due to replacement of dietary FO by VO.
The predominant influence of genotype was in the expression of intestinal transcripts of structural proteins, particularly collagen alpha chains, but also osteonectin, TAGLN, troponin I and keratocan, which were up-regulated in Lean fish, whereas troponin C was down-regulated. Furthermore, CNN1 and TAGLN were down-regulated in the intestinal proteome in Lean fish. Collagen, the main component of connective tissue, helps to maintain the structural integrity of tissues , while osteonectin is an extracellular matrix glycoprotein with high affinity towards collagen and whose expression has been associated with remodelling processes in tissues, including human intestine during development/morphogenesis and in diseased mucosa . Troponin, TAGLN and CNN1 are all involved in actin binding, actin-myosin interaction and muscle contraction. The inverse regulation of troponins is not conflicting as they have different roles in actomyosin cross-bridge formation and contraction; binding of troponin C to Ca2+ induces conformational changes that counteract the inhibitory action of troponin I . Expression of TAGLN transcript and protein showed opposite effects but a lack of correlation between transcriptomic and proteomic data is not unprecedented . As discussed above, this result might also be explained by the presence of similar duplicated genes in Atlantic salmon that are regulated differently. Transcriptomic results were validated by RT-qPCR for COL1A2, although only significantly when fish were fed the VO diet, for which fold changes were higher. In addition, in the microarray results differences in expression of structural proteins between family groups were consistently more accentuated in fish fed VO which could suggest a cumulative effect of diet. Furthermore, MYO was up-regulated in fish fed VO compared to FO but only in Lean fish (RT-qPCR), and significant diet × genotype interactions were found for alpha-actinin 1, tubulin beta-2 chain and procollagen-lysine 2-oxoglutarate 5-dioxygenase 2 (Additional file 1), which were up-regulated in Lean salmon, compared to Fat, but only when fed VO. In cod, replacement of FO by VO resulted in changes in intestinal expression of structural genes with the potential to alter the structural and mechanical properties of the intestinal muscle layer, including a range of actin-binding transcripts .
The present study is the first investigation of the influence of genetic background of families with different flesh adiposity phenotypes on intestinal gene expression of a fish species. Effects were subtle and consequently their potential impacts difficult to fully assess. However, if genetic selection for families better adapted to alternative formulations is an approach taken in the future, the potential for genotype-specific differences being exacerbated when VO replaces dietary FO should be further examined to assess the consequences of these changes in intestinal gene expression.
Metabolic activity, particularly lipid and energy, of intestinal tissue responded to dietary lipid composition but was also affected by genotype. The LC-PUFA biosynthesis pathway, typically up-regulated when salmon are fed VO, was especially influenced by genetic background. The Lean fish showed an enhanced response to low dietary n-3 LC-PUFA and the expression of Δ5fad, Δ6fad, elovl5b and elovl2 in the intestine showed high plasticity and was reflected in tissue biochemical composition indicating that their transcriptional regulation might be under feedback control by n-3 LC-PUFA, mainly DHA. Lower n-3 LC-PUFA in VO increased lipogenesis in Lean salmon, assessed by expression of FAS, while β-oxidation appeared unaffected, although transcripts involved in mitochondrial respiratory or electron transport chains were down-regulated, suggesting reduced activity in fish fed VO. Higher expression of genes and proteins involved in xenobiotic metabolism (CYP1A and EPHX2), antioxidant defence (CAT, HPX and PRDX1), and apoptosis (Casp3B) were observed in VO-fed fish, suggesting they might be responding to higher levels of contaminants, particularly PAH, in the diet. However, the intestine appeared able to metabolize and detoxify xenobiotic substances potentially present in the diet without major deleterious effects. Nonetheless, the data suggest that further attention should be given to contaminants in VO in the future. On the other hand, the data indicate potential genotype-specific differences in the response of the intestinal transcriptome and proteome to dietary VO. These include potential changes in structural properties of the intestinal layer and defence against cellular stress suggesting the Lean group was more susceptible to diet-induced oxidative stress. Considering the possibility of selecting families better adapted to alternative diet formulations, and the central role of intestine as a major barrier to nutrients, contaminants and pathogens, greater attention should be given to this organ when evaluating the effects of diet and genotype.
Feeding trial and sampling
A dietary trial was conducted using two genetically characterized groups of Atlantic salmon (Salmo salar) post-smolts comprising full-sib families selected from a breeding program (Landcatch Natural Selection Ltd., Argyll, Scotland). The choice of the family groups was based on estimated breeding values (EBVs) of the parents for high or low flesh adiposity, assessed by Torry Fatmeter (Distell Industries, West Lothian, UK), a trait that was found to have a heritability ranging from 0.17 to 0.39 in this dataset. The two groups were created from four unrelated full-sib families; two families from the extreme lower end of the EBV distribution for flesh lipid content (Lean) and two families from the extreme upper end of the distribution (Fat). The average EBV for the lipid content of the Fat families was 2.00 percentage units higher than that of the Lean families, representing a standardised selection differential of 2.33 standard deviations. Assessment of the flesh and visceral lipid contents at the end of the trial confirmed differences in adiposity between the groups .
Two thousand fish of each group were stocked into eight 12 x 5m3 net pens at the Ardnish Fish Trials Unit (Marine Harvest Scotland, Lochailort, Scotland; 500 fish pen-1). Duplicate pens of each group were fed one of two experimental diets (Skretting ARC, Stavanger, Norway) containing 25-32% fish meal, 40-45% plant meals and 27.5-30% oil (percentages varying according to pellet size) supplied either as standard northern FO or as a VO blend comprising rapeseed, palm and Camelina oils in a ratio of 5:3:2. Diets were formulated to fully satisfy the nutritional requirements of salmonid fish  and contained similar levels of PUFA but different n-3 and n-6 PUFA contents, 25.3% and 4.6% in the FO diet and 13.4% and 17.1% in the VO diet, respectively . After 55 weeks, 25 fish per pen were sampled 24 h after the last meal. Fish were killed by a blow to the head following anaesthesia, and intestinal tissue (pyloric caeca) collected, immediately frozen in liquid nitrogen and stored at −70°C prior to analyses. Further details can be found in Bell et al. .
Lipid extraction and fatty acid analyses
Total lipid from 1 g of intestine of four fish per treatment was extracted and determined gravimetrically , and fatty acid methyl esters (FAME) prepared by acid-catalysed transesterification of total lipid . FAME were separated and quantified by gas chromatography as described in detail previously . Significant differences in intestinal fatty acid composition were determined by two-way ANOVA (p<0.05) using the SPSS 16.0 statistical package (SPSS Inc., Chicago IL, USA).
RNA extraction and purification
Intestinal tissue (0.2 g) from six individuals per experimental group was homogenised in 2mL TRI Reagent and total RNA isolated following manufacturer’s instructions (Ambion, Applied Biosystems, Warrington, U.K.). RNA quantity and quality (integrity and purity) were assessed by gel electrophoresis and spectrophotometry (NanoDrop ND-1000, Thermo Scientific, Wilmington, U.S.A.), and 100 μg of total RNA from each sample further cleaned by mini spin-column purification (RNeasy Mini Kit, Qiagen, Crawley, UK).
Microarray hybridizations, image processing and statistical analysis
The TRAITS/SGP (v.2.1) salmon 17k cDNA microarray described by Taggart et al.  was used (ArrayExpress accession: A-MEXP-1930). A dual-labelled experimental design was employed, with each sample being competitively hybridised against a common pooled-reference. The experiment comprised 2 genotypes (Lean/Fat) × 2 diets (VO/FO) × 6 biological replicates. Indirect labelling was employed for preparing the microarray targets. Antisense amplified RNA (aRNA) was produced from 500 ng of purified total RNA per sample using the Amino Allyl MessageAmpTM II aRNA Amplification Kit as per manufacturer’s instructions (Ambion, Applied Biosystems), followed by Cy3 (for samples) or Cy5 (for pooled reference) fluor (GE HealthCare, Buckinghamshire, U.K.) incorporation through dye-coupling reaction . Microarray hybridizations were performed in a Lucidea semi-automated system (GE Healthcare) without pre-hybridization. For each array, every labelled biological replicate and corresponding pooled reference (40 pmol each dye, c. 150 ng aRNA) were combined and added to the hybridization solution. Two post-hybridization automatic washes followed by six manual washes to a final stringency of 0.1× SSC (EasyDipTM Slide staining system; Canemco Inc., Quebec, Canada) were performed before scanning.
Scanning was performed at 10 μm resolution using an Axon GenePix 4200AL Scanner (MDS Analytical Technologies, Wokingham, Berkshire, U.K.). Laser power was constant (80%) and “auto PMT” was enabled to adjust each channel at less than 0.1% feature saturation and Cy3/Cy5 mean intensity close to one. BlueFuse software (BlueGnome, Cambridge, U.K.) was used to identify features and extract fluorescence intensity values from TIF images. The resulting fluorescence intensity data (BlueFuse proprietary algorithm) and quality annotations for the 17,102 gene features were exported into the GeneSpring GX version 10.0.2 analysis platform (Agilent Technologies, Wokingham, Berkshire, U.K.) after undergoing block Lowess normalization. All control features were excluded from subsequent analyses. Data transformation and quality filtering were as in Morais et al. [9, 10]. This gave a final list of 15,498 genes that were eligible for statistical analysis. Experimental annotations complied fully with minimum information about a microarray experiment (MIAME) guidelines  and experimental hybridisations are archived on the EBI ArrayExpress database (http://www.ebi.ac.uk/arrayexpress/) under accession number E-TABM-1173. Hybridization data were analysed in GeneSpring by two-way ANOVA, which examined the explanatory power of the variable diet and genotype and the interaction between the two, followed by Gene Ontology (GO) enrichment analysis of the significant lists of features, at a significance level of 0.05. No multiple test correction was employed, as previous analyses, confirmed by RT-qPCR, indicated that such corrections are over-conservative for this type of data [9, 32].
RT-qPCR gene expression analysis
Expression of selected genes, for microarray validation and to further examine biological processes of interest, was studied by reverse transcription quantitative real time PCR (RT-qPCR), with target qPCR primer sequences given in Additional file 2. In addition, amplification of two reference genes, cofilin-2 and elongation factor-1α (elf-1α), was performed. One μg of column-purified total RNA per sample was reverse transcribed into cDNA using the VersoTM cDNA kit (ABgene, Surrey, U.K.) using a mixture of random hexamers (400ng/μL) and anchored oligo-dT (500ng/μL) at 3:1 (v/v). Negative controls (no enzyme) were performed to check for genomic DNA contamination. A similar amount of cDNA was pooled from all samples and the remaining cDNA diluted 20-fold with water. RT-qPCR analysis used relative quantification with the amplification efficiency of each primer pair assessed by serial dilutions of the cDNA pool. Amplifications were carried out in duplicate using a Quantica machine (Techne, Cambridge, U.K.) in a final volume of 20 μl containing 2–8 μl diluted cDNA, 0.5 μM of each primer and 10 μl AbsoluteTM QPCR SYBR® Green mix (ABgene), with a systematic negative control (NTC-non template control). The qPCR profiles contained an initial activation step at 95°C for 15 min, followed by 30–40 cycles: 15 s at 95°C, 15 s at the specific primer pair annealing temperature (Ta; Additional file 2) and 15 s at 72°C. After amplification, a melt curve was performed confirming a single product in each reaction, RT-qPCR product sizes checked by agarose gel electrophoresis, and identity of amplicons confirmed by sequencing. Gene expression was analysed using the relative expression software tool (REST 2008, http://www.gene-quantification.info/), employing a pair wise fixed reallocation randomisation test (10,000 randomisations) with efficiency correction .
Protein extraction and labelling
Six intestine samples (100 mg) per treatment were rapidly disrupted by homogenization and sonication on ice in 1 ml of DIGE lysis/labeling buffer (7M Urea, 2M Thiourea, 4% CHAPS and 30mM Tris, pH 8.5) in the presence of 10 μl of a protease inhibitor cocktail (ref. P8340, Sigma-Aldrich) and 4μl of 250 mM EDTA. After centrifugation at 12,000 × g for 20 min at 4°C, the supernatant was recovered and protein concentration determined (Quick Start Bradford Protein Assay Kit, ref. # 500–0202, Bio-Rad). Protein (500 μg) was purified by precipitation (ReadyPrepTM 2-D Cleanup kit, Bio-Rad) and the pellet re-suspended in DIGE lysis/labeling buffer at ~5μg/μl.
Samples were labelled using CyDye DIGE fluors (5nmol labelling kit; GE Healthcare), following manufacturer’s instructions. Three of the experimental replicates of each treatment (50 μg protein with pH corrected to 8.0-8.5) were labelled individually with 400 pmol Cy3 and the remaining three with 400 pmol Cy5. In addition, equal amounts of all experimental samples were pooled and 600 μg of protein (50 μg per 2-D gel) were batch-labelled with Cy2. The three labelled samples, corresponding to two experimental samples and one internal reference pool, were then combined to have in each 2-D gel samples corresponding to fish fed either FO or VO within the same family group.
Two-dimensional (2-D) polyacrylamide gel electrophoresis
Rehydration buffer (ReadyPrep 2-D starter kit, Bio-Rad) containing 0.2% DTT was added to the pooled protein samples to a final volume of 450 μl, which were loaded onto Immobiline DryStrip pH 3–11 NL, 24 cm IPG strips (GE Healthcare) by passive rehydration at room temperature overnight in the dark. Proteins were separated in the first dimension by isoelectric focusing at 20°C (Ettan IPGphor, GE Healthcare), applying increasing voltage until 200 V for 4 h, increasing to 500 V over a period of 3 h, then keeping the applied tension at a constant 1000 V for 1 h, followed by a further increase to 8000 V over 90 min, maintaining this voltage for almost 9 h (total of 79,200 Vh). After isoelectric focusing the strips were equilibrated in two 40 min steps using 50mM Tris–HCl pH 8.8, 6M urea, 30% glycerol, 2% SDS buffer, to which 2 % DTT (w/v) and 2.8% iodoacetamide (GE Healthcare) were added to produce reducing and alkylating buffers, respectively. The strips were loaded onto a 12.5% acrylamide gel cast between low fluorescence glass cassettes (EttanDALTsix gel caster system, GE Healthcare). The strips were overlaid with ReadyPrep Overlay Agarose (Bio-Rad) and the six gel cassettes run in the EttanDALT system in two steps: at 60 mA, 80 V, 6 W (upper limits) for 1 h, and then 240 mA, 500 V, 78 W until the bromophenol blue dye front had run to ~1 cm above the bottom of the gels. Laemmli buffers (250mM Tris, 1.92M glycine, 1% SDS) were used in the lower (1X) and upper (2X) chambers, respectively.
Gel imaging and analysis
Labelled gels were scanned using a Typhoon TRIO (GE Healthcare) and Cy2, Cy3 and Cy5 images acquired using 520BP40, 580BP30 and 670BP30 laser emission filters, respectively, at 500 PMT and 100 μm resolution. Images were cropped to remove extraneous areas prior to analysis (ImageQuantTL, GE Healthcare), and image analysis performed using DeCyder V7.0 (GE Healthcare). The estimated number of spots for each co-detection procedure was set at 10,000 and an exclusion filter was applied to remove spots with a volume lower than 30,000. Differential expression of protein spots was examined by two-way ANOVA at a significance level of 0.05. After verifying that significant spots were well matched across the gels, two pick lists were generated with a total of 22 and 45 spots for the diet and genotype factors, respectively (Figure 1; Additional files 3 and 4).
Spot picking and protein identification by peptide fragment fingerprinting
Four preparative gels were run under the conditions described above but with higher amounts of protein (600 and 1000 μg). They were stained with colloidal Coomassie and, whenever possible, spots were excised and sequenced in the Mass Spectrometry Laboratory ITQB-UNL (Oeiras, Portugal), where in-gel digestion and extraction of the proteins from the gel was performed, followed by micropurification, and peptides identified by mass spectrometry (using an Applied Biosystems (ABI) 4800 MALDI TOF/TOF Analyzer). The search engine MASCOT was then used to identify and confirm protein ID’s from the peptide mass fingerprinting (PMF; MS) and peptide fragment fingerprinting (PFF; MS/MS) data.
Calder PC, Yaqoob P: Omega-3 polyunsaturated fatty acids and human health outcomes. Biofactors. 2009, 35: 266-272. 10.1002/biof.42.
FAO: The State of World Fisheries and Aquaculture (SOFIA) 2010. 2010, Food and Agriculture Organization of the United Nations, Rome
Tacon AGJ, Metian M: Global overview on the use of fish meal and fish oil in industrially compounded aquafeeds: Trends and future prospects. Aquaculture. 2008, 285: 146-158. 10.1016/j.aquaculture.2008.08.015.
Torstensen BE, Espe M, Sanden M, Stubhaug I, Waagbø R, Hemre G-I, Fontanillas R, Nordgarden U, Hevrøy EM, Olsvik P, Berntssen MHG: Novel production of Atlantic salmon (Salmo salar) protein based on combined replacement of fish meal and fish oil with plant meal and vegetable oil blends. Aquaculture. 2008, 285: 193-200. 10.1016/j.aquaculture.2008.08.025.
Bell JG, Pratoomyot J, Strachan F, Henderson RJ, Fontanillas R, Hebard A, Guy DR, Hunter D, Tocher DR: 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: 225-232. 10.1016/j.aquaculture.2010.05.021.
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.
Panserat S, Kolditz C, Richard N, Plagnes-Juan E, Piumi F, Esquerré D, Médale F, Corraze G, Kaushik S: Hepatic gene expression profiles in juvenile rainbow trout (Oncorhynchus mykiss) fed fishmeal or fish oil-free diets. Br J Nutr. 2008, 100: 953-967. 10.1017/S0007114508981411.
Panserat S, Hortopan GA, Plagnes-Juan E, Kolditz C, Lansard M, Skiba-Cassy S, Esquerre D, Geurden I, Médale F, Kaushik S, Corraze G: Differential gene expression after total replacement of dietary fish meal and fish oil by plant products in rainbow trout (Oncorhynchus mykiss) liver. Aquaculture. 2009, 294: 123-131. 10.1016/j.aquaculture.2009.05.013.
Morais S, Pratoomyot J, Taggart JB, Bron JE, Guy DR, Bell JG, Tocher DR: 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: 255-10.1186/1471-2164-12-255.
Morais S, Pratoomyot J, Torstensen BE, Taggart JB, Guy DR, Bell JG, Tocher DR: Diet x genotype interactions in hepatic cholesterol and lipoprotein metabolism in Atlantic salmon (Salmo salar) in response to replacement of dietary fish oil with vegetable oil. Br J Nutr. 2011, 106: 1457-1469. 10.1017/S0007114511001954.
Geay F, Ferraresso S, Zambonino-Infante JL, Bargelloni L, Quentel C, Vandeputte M, Kaushik S, Cahu CL, Mazurais D: Effects of the total replacement of fish-based diet with plant-based diet on the hepatic transcriptome of two European sea bass (Dicentrarchus labrax) half-sibfamilies showing different growth rates with the plant-based diet. BMC Genomics. 2011, 12: 522-10.1186/1471-2164-12-522.
Murray HM, Lall SP, Rajaselvam R, Boutilier LA, Blanchard B, Flight RM, Colombo S, Mohindra V, Douglas SE: A nutrigenomic analysis of intestinal response to partial soybean meal replacement in diets for juvenile Atlantic halibut, Hippoglossus hipploglossus, L. Aquaculture. 2010, 298: 282-293. 10.1016/j.aquaculture.2009.11.001.
Skugor S, Grisdale-Helland B, Refstie S, Afanasyev S, Vielma J, Krasnov A: Gene expression responses to restricted feeding and extracted soybean meal in Atlantic salmon (Salmo salar L.). Aquacult Nutr. 2011
Bell MV, Dick JR, Porter AEA: Pyloric ceca are a major site of 22:6n-3 synthesis in rainbow trout (Oncorhynchus mykiss). Lipids. 2003, 39: 39-44.
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: 317-328. 10.1007/s10695-006-9109-2.
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.
Olsen RE, Dragsnes BT, Myklebust R, Ringø R: Effect of soybean oil and soybean lecithin on intestinal lipid composition and lipid droplet accumulation of rainbow trout, Oncorhynchus mykiss Walbaum. Fish Physiol Biochem. 2003, 29: 181-192.
Olsen RE, Myklebust R, Kaino T, Ringø E: Lipid digestibility and ultrastructural changes in the enterocytes of Arctic char (Salvelinus alpinus L.) fed linseed oil and soybean lecithin. Fish Physiol Biochem. 1999, 21: 35-44. 10.1023/A:1007726615889.
Olsen RE, Myklebust R, Ringø 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 Physiol Biochem. 2000, 22: 207-216. 10.1023/A:1007879127182.
Morais S, Edvardsen RB, Tocher DR, Bell JG: Transcriptomic analyses of intestinal gene expression of juvenile Atlantic cod (Gadus morhua) fed diets with Camelina oil as replacement for fish oil. Comp Biochem Physiol B. -
Olsvik PA, Torstensen BE, Berntssen MHG: Effects of complete replacement of fish oil with plant oil on gastrointestinal cell death, proliferation and transcription of eight genes’ encoding proteins responding to cellular stress in Atlantic salmon Salmo salar L. J Fish Biol. 2007, 71: 550-568. 10.1111/j.1095-8649.2007.01521.x.
Teitelbaum JE, Allan Walker W: Review: the role of omega 3 fatty acids in intestinal inflammation. J Nutr Biochem. 2001, 12: 21-32. 10.1016/S0955-2863(00)00141-8.
Leaver MJ, Taggart JB, Villeneuve L, Bron JE, Guy DR, Bishop SC, Houston RD, Matika O, Tocher DR: 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: 62-69. 10.1016/j.cbd.2010.04.002.
Dale OB, Tørud B, Kvellestad A, Koppang HS, Koppang EO: From chronic feed-induced intestinal inflammation to adenocarcinoma with metastases in salmonid fish. Cancer Res. 2009, 69: 4355-4362. 10.1158/0008-5472.CAN-08-4877.
Burrells C, Williams PD, Southgate PJ, Crampton VO: Immunological, physiological and pathological responses of rainbow trout (Oncorhynchus mykiss) to increasing dietary concentrations of soybean proteins. Vet Immunol Immunopathol. 1999, 72: 277-288. 10.1016/S0165-2427(99)00143-9.
Kolditz CI, Paboeuf G, Borthaire M, Esquerré D, SanCristobal M, Lefèvre F, Médale F: 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: 506-10.1186/1471-2164-9-506.
Allendorf FW, Thorgaard GH: Tetraploidy and the evolution of salmonid fishes. Evolutionary Genetics of Fishes. Edited by: Turner BJ. 1984, Plenum Press, New York, 55-93.
Jump DB, Clarke SD, Thelen A, Liimatta M, Ren B, Badin M: Dietary polyunsaturated fatty acid regulation of gene transcription. Prog Lipid Res. 1996, 35: 227-241. 10.1016/S0163-7827(96)00007-0.
Nakamura MT, Nara TY: Gene regulation of mammalian desaturases. Biochem Soc Trans. 2002, 30: 1076-1079.
Wang Y, Botolin D, Christian B, Busik J, Xu J, Jump DB: Tissue-specific, nutritional, and developmental regulation of rat fatty acid elongases. J Lipid Res. 2005, 46: 706-715. 10.1194/jlr.M400335-JLR200.
Leaver MJ, Villeneuve LA, Obach A, Jensen L, Bron JE, Tocher DR, Taggart JB: 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. 2008, 9: 299-10.1186/1471-2164-9-299.
Taggart JB, Bron JE, Martin SAM, Seear PJ, Høyheim B, Talbot R, Carmichael SN, Villeneuve LAN, Sweeney GE, Houlihan DF, Secombes CJ, Tocher DR, Teale AJ: A description of the origins, designs and performance of the TRAITS-SGP Atlantic salmon Salmo salar L. cDNA microarray. J Fish Biol. 2008, 72: 2071-2094. 10.1111/j.1095-8649.2008.01876.x.
Morais S, Monroig O, Zheng X, Leaver MJ, Tocher DR: Highly unsaturated fatty acid synthesis in Atlantic salmon: characterization of ELOVL5- and ELOVL2-like elongases. Mar Biotechnol. 2009, 1: 627-639.
van Schothorst EM, Flachs P, Franssen-van Hal NL, Kuda O, Bunschoten A, Molthoff J, Vink C, Hooiveld GJ, Kopecky J, Keijer J: Induction of lipid oxidation by polyunsaturated fatty acids of marine origin in small intestine of mice fed a high-fat diet. BMC Genomics. 2009, 10: 110-10.1186/1471-2164-10-110.
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.
Rosen ED, Walkey CJ, Puigserver P, Spiegelman BM: Transcriptional regulation of adipogenesis. Genes Dev. 2000, 14: 1293-1307.
Schmid B, Rippmann JF, Tadayyon M, Hamilton BS: Inhibition of fatty acid synthase prevents preadipocyte differentiation. Biochem Biophys Res Commun. 2005, 328: 1073-1082. 10.1016/j.bbrc.2005.01.067.
Davidson MH: Mechanisms for the hypotriglyceridemic effect of marine omega-3 fatty acids. Am J Cardiol. 2006, 98: 27i-33i.
Mori T, Kondo H, Hase T, Tokimitsu I, Murase T: Dietary fish oil upregulates intestinal lipid metabolism and reduces body weight gain in C57BL/6J mice. J Nutr. 2007, 137: 2629-2634.
Madsen L, Frøyland L, Kryvi H, Auwerx J, Staels B, Berge RK: The mitochondrion is the principal target for nutritional and pharmacological control of plasma triglyceride. Lipids. 1999, 34 (Suppl): S167-
Kjær MA, Vegusdal A, Gjøen T, Rustan AC, Todorcević M, Ruyter B: Effect of rapeseed oil and dietary n-3 fatty acids on triacylglycerol synthesis and secretion in Atlantic salmon hepatocytes. Biochim Biophys Acta. 2008, 178: 112-122.
Jacobs MN, Covaci A, Schepens P: Investigation of selected persistent organic pollutants in farmed Atlantic salmon (Salmo salar), salmon aquaculture feed, and fish oil components of the feed. Environ Sci Technol. 2002, 36: 2797-2805. 10.1021/es011287i.
Bell JG, McGhee F, Dick JR, Tocher DR: Dioxin and dioxin-like polychlorinated biphenyls (PCBs) in Scottish farmed salmon (Salmo salar): effects of replacement of dietary marine fish oil with vegetable oils. Aquaculture. 2005, 243: 305-314. 10.1016/j.aquaculture.2004.10.016.
Takahashi M, Tsuboyama-Kasaoka N, Nakatani T, Ishii M, Tsutsumi S, Aburatani H, Ezaki O: Fish oil feeding alters liver gene expressions to defend against PPARα activation and ROS production. Am J Physiol Gastrointest Liver Physiol. 2002, 282: G338-G348.
Saera-Vila A, Benedito-Palos L, Sitjà-Bobadilla A, Nácher-Mestre J, Serrano R, Kaushik S, Pérez-Sánchez J: Assessment of the health and antioxidant trade-off in gilthead sea bream (Sparus aurata L.) fed alternative diets with low levels of contaminants. Aquaculture. 2009, 296: 87-95. 10.1016/j.aquaculture.2009.07.028.
Moore MJ, Mitrofanov IV, Valentini SS, Volkov VV, Kurbskiy AV, Zhimbey EN, Eglinton LB, Stegeman JJ: Cytochrome P4501A expression, chemical contaminants and histopathology in roach, goby and sturgeon and chemical contaminants in sediments from the Caspian Sea, Lake Balkhash and the Ily River Delta, Kazakhstan. Mar Pollut Bull. 2003, 46: 107-119. 10.1016/S0025-326X(02)00325-9.
Lee RF, Anderson JW: Significance of cytochrome P450 system responses and levels of bile fluorescent aromatic compounds in marine wildlife following oil spills. Mar Pollut Bull. 2005, 50: 705-723. 10.1016/j.marpolbul.2005.04.036.
Bacanskas LR, Whitaker J, Di Giulio RT: Oxidative stress in two populations of killifish (Fundulus heteroclitus) with differing contaminant exposure histories. Mar Environ Res. 2004, 58: 597-601. 10.1016/j.marenvres.2004.03.048.
Galar-Martínez M, Gomez-Olivan LM, Amaya-Chavez A, Razo-Estrada C, Garcia-Medina S: Oxidative stress induced on Cyprinus carpio by contaminants present in the water and sediment of Madín Reservoir. J Environ Sci Health, Part A: Tox Hazard Subst Environ Eng. 2010, 45: 155-160. 10.1080/10934520903425780.
Wood ZA, Schröder E, Robin Harris J, Poole LB: Structure, mechanism and regulation of peroxiredoxins. Trends Biochem Sci. 2003, 28: 32-40. 10.1016/S0968-0004(02)00003-8.
Tolosano E, Altruda F: Hemopexin: structure, function, and regulation. DNA Cell Biol. 2002, 21: 297-306. 10.1089/104454902753759717.
Burk RF, Hill KE, Motley AK: Selenoprotein metabolism and function: evidence for more than one function for selenoprotein P. J Nutr. 2003, 133 (5 Suppl 1): 1517S-1520S.
Kang YJ: Metallothionein redox cycle and function. Exp Biol Med. 2006, 231: 1459-1467.
Berntssen MH, Julshamn K, Lundebye AK: Chemical contaminants in aquafeeds and Atlantic salmon (Salmo salar) following the use of traditional- versus alternative feed ingredients. Chemosphere. 2010, 78: 637-646. 10.1016/j.chemosphere.2009.12.021.
Moret S, Conte LS: Polycyclic aromatic hydrocarbons in edible fats and oils: occurrence and analytical methods. J Chromatogr A. 2000, 882: 245-253. 10.1016/S0021-9673(00)00079-0.
Shimada T: Xenobiotic-metabolizing enzymes involved in activation and detoxification of carcinogenic polycyclic aromatic hydrocarbons. Drug Metab Pharmacokinet. 2006, 21: 257-276. 10.2133/dmpk.21.257.
Hellou J, Leonard J, Anstey C: Dietary exposure of finfish to aromatic contaminants and tissue distribution. Arch Environ Contam Toxicol. 2002, 42: 470-476. 10.1007/s00244-001-0038-x.
Yuen BB, Wong CK, Woo NY, Au DW: Induction and recovery of morphofunctional changes in the intestine of juvenile carnivorous fish (Epinephelus coioides) upon exposure to foodborne benzo[a]pyrene. Aquat Toxicol. 2007, 82: 181-194. 10.1016/j.aquatox.2007.02.010.
Zhang Y, Lu Y, Zhou H, Lee M, Liu Z, Hassel BA, Hamburger AW: Alterations in cell growth and signaling in ErbB3 binding protein-1 (Ebp1) deficient mice. BMC Cell Biol. 2008, 9: 69-10.1186/1471-2121-9-69.
Celis JE, Madsen P, Celis A, Nielsen HV, Gesser B: Cyclin (PCNA, auxiliary protein of DNA polymerase delta) is a central component of the pathway(s) leading to DNA replication and cell division. FEBS Lett. 1987, 220: 1-7. 10.1016/0014-5793(87)80865-7.
Liu Z, Ahn JY, Liu X, Ye K: Ebp1 isoforms distinctively regulate cell survival and differentiation. Proc Natl Acad Sci USA. 2006, 103: 10917-10922. 10.1073/pnas.0602923103.
Sweet LI, Passino-Reader DR, Meier PG, Omann GM: Xenobiotic-induced apoptosis: significance and potential application as a general biomarker of response. Biomarkers. 1999, 4: 237-253. 10.1080/135475099230778.
Ramachandran A, Madesh M, Balasubramanian KA: Apoptosis in the intestinal epithelium: its relevance in normal and pathophysiological conditions. J Gastroenterol Hepatol. 2000, 15: 109-120. 10.1046/j.1440-1746.2000.02059.x.
Paclik D, Lohse K, Wiedenmann B, Dignass AU, Sturm A: Galectin-2 and −4, but not galectin-1, promote intestinal epithelial wound healing in vitro through a TGF-beta-independent mechanism. Inflamm Bowel Dis. 2008, 14: 1366-1372. 10.1002/ibd.20499.
Loser K, Sturm A, Voskort M, Kupas V, Balkow S, Auriemma M, Sternemann C, Dignass AU, Luger TA, Beissert S: Galectin-2 suppresses contact allergy by inducing apoptosis in activated CD8+ T cells. J Immunol. 2009, 182: 5419-5429. 10.4049/jimmunol.0802308.
Tsukamoto K, Miura F, Fujito NT, Yoshizaki G, Nonaka M: Long-Lived Dichotomous Lineages of the Proteasome Subunit Beta Type 8 (PSMB8) Gene Surviving More than 500 Million Years as Alleles or Paralogs. Mol Biol Evol. 2012, May 2. [Epub ahead of print]
Nandi D, Tahiliani P, Kumar A, Chandu D: The ubiquitin-proteasome system. J Biosci. 2006, 31: 137-155. 10.1007/BF02705243.
Widmer CC, Pereira CP, Gehrig P, Vallelian F, Schoedon G, Buehler PW, Schaer DJ: Hemoglobin can attenuate hydrogen peroxide-induced oxidative stress by acting as an antioxidative peroxidase. Antioxid Redox Signal. 2010, 12: 185-198. 10.1089/ars.2009.2826.
Shimizu S, Matsuoka Y, Shinohara Y, Yoneda Y, Tsujimoto Y: Essential role of voltage-dependent anion channel in various forms of apoptosis in mammalian cells. J Cell Biol. 2001, 152: 237-250. 10.1083/jcb.152.2.237.
Jeon YJ, Kim DH, Jung H, Chung SJ, Chi SW, Cho S, Lee SC, Park BC, Park SG, Bae KH: Annexin A4 interacts with the NF-kappaB p50 subunit and modulates NF-kappaB transcriptional activity in a Ca2+−dependent manner. Cell Mol Life Sci. 2010, 67: 2271-2281. 10.1007/s00018-010-0331-9.
Parag HA, Raboy B, Kulka RG: Effect of heat shock on protein degradation in mammalian cells: involvement of the ubiquitin system. EMBO J. 1987, 6: 55-61.
Raboy B, Sharon G, Parag HA, Shochat Y, Kulka RG: Effect of stress on protein degradation: role of the ubiquitin system. Acta Biol Hung. 1991, 42: 3-20.
Pirkkala L, Alastalo TP, Zuo X, Benjamin IJ, Sistonen L: Disruption of heat shock factor 1 reveals an essential role in the ubiquitin proteolytic pathway. Mol Cell Biol. 2000, 20: 2670-2675. 10.1128/MCB.20.8.2670-2675.2000.
Di Lullo GA, Sweeney SM, Korkko J, Ala-Kokko L, San Antonio JD: Mapping the ligand-binding sites and disease-associated mutations on the most abundant protein in the human, type I collagen. J Biol Chem. 2002, 277: 4223-4231. 10.1074/jbc.M110709200.
Lussier C, Sodek J, Beaulieu JF: Expression of SPARC/osteonectin/BM4O in the human gut: predominance in the stroma of the remodeling distal intestine. J Cell Biochem. 2001, 81: 463-476. 10.1002/1097-4644(20010601)81:3<463::AID-JCB1060>3.0.CO;2-L.
Layland J, Solaro RJ, Shah AM: Regulation of cardiac contractile function by troponin I phosphorylation. Cardiovasc Res. 2005, 66: 12-21. 10.1016/j.cardiores.2004.12.022.
Olivares-Hernández R, Usaite R, Nielsen J: Integrative analysis using proteome and transcriptome data from yeast to unravel regulatory patterns at post-transcriptional level. Biotechnol Bioeng. 2010, 107: 865-875. 10.1002/bit.22868.
National Research Council (NRC): Nutrient requirements of fish and shrimp. 2011, National Academy Press, Washington
Folch J, Lees M, Sloane-Stanley GH: A simple method for the isolation and purification of total lipides from animal tissues. J Biol Chem. 1957, 226: 497-509.
Christie WW: Lipid analysis. 2003, The Oily Press, Bridgewater, UK, 3
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 FC, 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. Nat Genet. 2001, 29: 365-371. 10.1038/ng1201-365.
Pfaffl MW, Horgan GW, Dempfle L: Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002, 30: e36-10.1093/nar/30.9.e36.
Zheng X, Tocher DR, Dickson CA, Bell JG, Teale AJ: Highly unsaturated fatty acid synthesis in vertebrates: new insights with the cloning and characterization of a Δ6 desaturase of Atlantic salmon. Lipids. 2005, 40: 13-24. 10.1007/s11745-005-1355-7.
Torstensen BE, Nanton DA, Olsvik PA, Sundvold H, Stubhaug I: Gene expression of fatty acid-binding proteins, fatty acid transport proteins (cd36 and FATP) and β-oxidation-related genes in Atlantic salmon (Salmo salar L.) fed fish oil or vegetable oil. Aquacult Nutr. 2009, 15: 440-451. 10.1111/j.1365-2095.2008.00609.x.
Rees CB, Li W: Development and application of a real-time quantitative PCR assay for determining CYP1A transcripts in three genera of salmonids. Aquat Toxicol. 2004, 66: 357-368. 10.1016/j.aquatox.2003.10.004.
Kleveland EJ, Ruyter B, Vegusdal A, Sundvold H, Berge RK, Gjøen T: Effects of 3-thia fatty acids on expression of some lipid related genes in Atlantic salmon (Salmo salar L.). Comp Biochem Physiol B Biochem Mol Biol. 2006, 145: 239-248. 10.1016/j.cbpb.2006.07.012.
Olsvik PA, Kristensen T, Waagbø R, Rosseland BO, Tollefsen KE, Baeverfjord G, Berntssen MHG: mRNA expression of antioxidant enzymes (SOD, CAT and GSH-Px) and lipid peroxidative stress in liver of Atlantic salmon (Salmo salar) exposed to hyperoxic water during smoltification. Comp Biochem Physiol C Toxicol Pharmacol. 2005, 141: 314-323. 10.1016/j.cbpc.2005.07.009.
Takle H, McLeod A, Andersen O: Cloning and characterization of the executioner caspases 3, 6, 7 and Hsp70 in hyperthermic Atlantic salmon (Salmo salar) embryos. Comp Biochem Physiol B Biochem Mol Biol. 2006, 144: 188-198. 10.1016/j.cbpb.2006.02.006.
This study was funded by the EU FP6 IP “AQUAMAX” (Sustainable Aquafeeds to Maximise the Health Benefits of Farmed Fish for Consumers; 016249–2). SM was supported by a Marie Curie Intra European Fellowship (FP7-PEOPLE-2007-2-1-IEF, Proposal N219667) and by AQUAGENOME Mobility and Resource Exchange Grants in 2008. The authors acknowledge Landcatch Ltd for provision of the Lean and Fat smolts and thank staff of Marine Harvest Ltd, Ardnish FTU, for assistance with fish husbandry and sample collection. Technical assistance from Jacquie Ireland in microarray hybridizations is deeply appreciated. We also acknowledge Mahaut de Vareilles for performing some of the protein identifications by mass spectrometry, and discussions on the use of DeCyder V7.0.
The authors declare that they have no competing interests.
SM performed laboratory work and data analysis, with assistance from TS and OC in proteomics; DRG was responsible for family selection; JBT and JEB supported the microarray analysis; PR supported the proteomic analysis; SM wrote the first draft of the manuscript, followed by contributions from remaining authors; SM, JGB and DRT planned and coordinated the research; JGB and DRT were project leaders. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Intestine transcripts corresponding to the top 100 most significant features exhibiting a significant diet × family interaction. (DOC 127 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Morais, S., Silva, T., Cordeiro, O. et al. Effects of genotype and dietary fish oil replacement with vegetable oil on the intestinal transcriptome and proteome of Atlantic salmon (Salmo salar). BMC Genomics 13, 448 (2012). https://doi.org/10.1186/1471-2164-13-448