Besides being a source of high-quality protein and essential micronutrients for humans, fish are unique sources of n-3 long-chain PUFA (LC-PUFA), namely EPA (20 : 5n-3) and DHA (22 : 6n-3), which were proven to be beneficial for human health( Reference Tocher 1 ). Driven by the crescent awareness of the health beneficial effects of n-3 LC-PUFA in a range of human pathologies (including CVD, inflammatory and neurological diseases) the global fish consumption is rising( Reference Tocher 1 ) and an increasing proportion of this fish is now being supplied by aquaculture( 2 ).
Fishmeal (FM) and fish oil (FO) have been widely used as main ingredients in aquafeeds for carnivorous fish species. However, environmental sustainability and economical issues, related to the limited availability of fisheries resources and escalating costs, forced the aquafeed industry to search for alternative and eco-friendly ingredients such as plant feedstuffs and vegetable oils (VO). However, the unbalanced amino acid and fatty acid (FA) profiles, relatively high amounts of carbohydrates and the presence of antinutritional factors may limit their use in aquafeeds, especially for carnivorous fish species that are metabolically adapted to diets rich in LC-PUFA and almost devoid of carbohydrates( Reference Tocher 1 , Reference Enes, Panserat and Kaushik 3 – Reference Stone 5 ).
It is well known that modifications of dietary macronutrients (e.g. lipid source, carbohydrates, etc.) may have marked effects on tissue FA composition and lipid deposition of fish( Reference Betancor, Sprague and Usher 6 – Reference Kjaer, Vegusdal and Gjoen 9 ). For instance, replacement of FO by VO generally leads to a decrease of n-3 LC-PUFA and an increase of 18C PUFA precursors, linoleic acid (18 : 2n-6) and α-linolenic acid (18 : 3n-3), in the fillet( Reference Castro, Corraze and Pérez-Jiménez 7 , Reference Kjaer, Vegusdal and Gjoen 9 – Reference Morais, Edvardsen and Tocher 11 ). Increased tissue lipid deposition with dietary incorporation of VO( Reference Betancor, Sprague and Usher 6 , Reference Jordal, Lie and Torstensen 8 , Reference Kjaer, Vegusdal and Gjoen 9 ) or carbohydrates( Reference Castro, Corraze and Pérez-Jiménez 7 , Reference Brauge, Medale and Corraze 12 ) was also reported. Such effects may occur because of different mechanisms that are recognised to regulate the quantity and quality of the fish lipid depots, such as modification of lipogenesis, β-oxidation, tissue lipid uptake and transport or FA desaturation and elongation processes( Reference Sheridan 13 , Reference Tocher 14 ).
Although not always consistent, several studies in fish reported that dietary lipid source and carbohydrates regulate gene expression, key transcription factors and/or activity of enzymes involved in lipogenesis, β-oxidation and lipid uptake. Accordingly, it was observed that replacing dietary FO by VO decreased( Reference Menoyo, Izquierdo and Robaina 15 , Reference Panserat, Kolditz and Richard 16 ) or increased( Reference Morais, Pratoomyot and Taggart 17 – Reference Peng, Xu and Mai 19 ) gene expression or activity of lipogenic enzymes (such as glucose-6-phosphate dehydrogenase (G6PD), malic enzyme (ME), fatty acid synthase (FAS)). Similarly, dietary VO was reported to regulate gene expression or activities of lipolytic enzymes either positively( Reference Liang, Ogata and Oku 20 , Reference Stubhaug, Froyland and Torstensen 21 ) or negatively( Reference Peng, Xu and Mai 19 , Reference Liang, Ogata and Oku 20 , Reference Jordal, Torstensen and Tsoi 22 , Reference Turchini, Mentasti and Froyland 23 ) (including carnitine palmitoyltransferase (CPT) I and II, lipoprotein lipase) in a tissue-specific manner. Dietary carbohydrate or glucose administration was also reported to enhance lipogenesis( Reference Kamalam, Médale and Kaushik 24 ). However, regarding β-oxidation, data are contradictory, as either stimulation( Reference Kamalam, Médale and Larroquet 25 ) or inhibition( Reference Jin, Médale and Kamalam 26 ) effects were reported. Regarding LC-PUFA biosynthesis, strong induction of desaturase (fatty acyl desaturase (FAD) Δ5 and Δ6 desaturases, FADS1 and FADS2, respectively) and elongase expression by VO administration was reported in freshwater fish and salmonids( Reference Leaver, Villeneuve and Obach 10 , Reference Seiliez, Panserat and Kaushik 27 ), but in marine fish such induction of gene expression is not so clear( Reference Morais, Edvardsen and Tocher 11 , Reference González-Rovira, Mourente and Zheng 28 – Reference Vagner and Santigosa 30 ). Desaturases and elongases were also shown to be up-regulated by dietary carbohydrates in salmonids( Reference Kamalam, Médale and Larroquet 25 , Reference Seiliez, Panserat and Kaushik 27 ), but in marine fish such an effect has never been demonstrated( Reference Castro, Corraze and Pérez-Jiménez 7 ).
Besides the reported effects on tissue FA composition and lipid deposition of fish, the use of plant feedstuffs and VO can also compromise fish physiological functions and, ultimately, fish health( Reference Montero and Izquierdo 31 , Reference Torstensen and Tocher 32 ). For example, there is an increased awareness of potential effects due to reduced dietary phospholipids (PL) or cholesterol in plant feedstuff-based diets. Besides their important roles in membrane structure, cholesterol and PL also have important functional roles. Cholesterol is a precursor of physiologically active compounds such as bile acids, vitamin D, adrenal corticoids and sex hormones, and PL are precursors of eicosanoids, diacylglycerol, inositol phosphates and platelet-activating factors( 33 , Reference Tocher, Bendiksen and Campbell 34 ). In a number of fish species, dietary VO was reported to decrease plasma PL( Reference Castro, Corraze and Pérez-Jiménez 7 , Reference Richard, Kaushik and Larroquet 35 , Reference Richard, Mourente and Kaushik 36 ), cholesterol and LDL-cholesterol( Reference Jordal, Lie and Torstensen 8 , Reference Richard, Kaushik and Larroquet 35 , Reference Morais, Pratoomyot and Torstensen 37 ), and to modulate the activity or expression of genes involved in PL synthesis( Reference Caballero, Gallardo and Robaina 38 ) and cholesterol synthesis and absorption( Reference Castro, Corraze and Pérez-Jiménez 7 , Reference Leaver, Villeneuve and Obach 10 , Reference Liland, Espe and Rosenlund 39 ). Recently, it was also demonstrated that dietary carbohydrates regulate plasma cholesterol and PL concentrations and the transcription of proteins involved in cholesterol metabolism( Reference Castro, Corraze and Pérez-Jiménez 7 ).
Therefore, to provide adequate background for successful use of plant feedstuffs in aquafeeds, all aspects related to fish physiological functions and the nutritional quality of the final product must be better understood. For that purpose, we assessed the long-term effect of dietary lipid source, carbohydrate content and interactions between both on growth performance, tissue composition, liver and intestine enzymatic activity and expression of genes related with lipid metabolism (lipogenesis, β-oxidation, FA bioconversion, cholesterol and PL metabolism) in gilthead sea bream juveniles. In this study, the selection of the relevant genes involved in lipid metabolism was performed taking advantage of the recent advances in the molecular and functional characterisation of a number of new gilthead sea bream sequences related to FA, PL and cholesterol metabolism( Reference Agaba, Tocher and Zheng 40 – Reference Seiliez, Panserat and Corraze 47 ). In addition, as an increasing number of studies have reported that glucose metabolism was distinctly regulated by different dietary lipid sources( Reference Castro, Corraze and Pérez-Jiménez 7 , Reference Morais, Pratoomyot and Taggart 17 , Reference Castro, Corraze and Panserat 48 – Reference Menoyo, Diez and Lopez-Bote 50 ), we also investigated the effects of these dietary manipulations on mechanisms involved in glucose utilisation/metabolism (such as glycolysis and gluconeogenesis pathways).
Methods
Experimental diets
Four diets differing in carbohydrate content (0 and 20 % gelatinised starch, diets CH– and CH+, respectively) and lipid source (diets FO or VO) were formulated (Table 1). The increase in carbohydrate content in CH+ diets was achieved by decreasing protein, which was kept well above the requirements of the species( Reference Oliva-Teles 51 ). The major lipid source of FO diets was cod liver oil. In VO diets, 100 % of the cod liver oil was replaced by a VO blend composed of 20 % rapeseed, 50 % linseed and 30 % palm oils. FM was added as a major dietary protein source to isolate the impacts of dietary VO and to avoid the interference of dietary plant protein on lipid metabolism, especially on cholesterol metabolism.
FO, fish oil; VO, blend of vegetable oils; CH, carbohydrate; CH content, 0 % (CH–) or 20 % (CH+) gelatinised maize starch.
* Steam-dried low-temperature fishmeal (Superprime; Inproquisa) (CP: 74·6 % DM; crude lipid: 10·1 % DM).
† C-Gel Instant-12018; Cerestar.
‡ Labchem; Laborspirit Lda.
§ 30 % palm oil (Colmi), 50 % linseed oil (Sociedade Portuense de Drogas) and 20 % rapeseed oil (Huilerie Emile Noël S.A.S).
|| Vitamins (mg/kg diet): retinyl acetate, 18 000 IU/kg diet (6·19mg); cholecalciferol, 2000 IU/kg diet (0·04mg); α-tocopheryl acetate, 35; sodium menadione bisulphate, 10; thiamin-HCl, 15; riboflavin, 25; calcium pantothenate, 50; nicotinic acid, 200; pyridoxine HCl, 5; folic acid, 10; cyanocobalamin, 0·02; biotin, 1·5; ascorbic acid, 50; inositol, 400 (premix).
¶ Minerals (mg/kg diet): cobalt sulphate, 1·91; copper sulphate, 19·6; iron sulphate, 200; sodium fluoride, 2·21; potassium iodide, 0·78; magnesium oxide, 830; manganese oxide, 26; sodiumselenite, 0·66; zinc oxide, 37·5; dibasic calcium phosphate, 5·93 (g/kg diet); potassium chloride, 1·15 (g/kg diet); sodium chloride, 0·40 (g/kg diet) (premix).
** Aquacube (guar gum, polymethyl carbamide, manioc starch blend, hydrate calcium sulphate) (Agil).
Animals, experimental conditions and sampling
The experiment was directed by trained scientists (following FELASA category C recommendations) and conducted according to the European Union Directive (2010/63/EU) on the protection of animals for scientific purposes. The study was performed at the Marine Zoological Station, University of Porto, Portugal, in a thermoregulated recirculation water system equipped with twelve fibreglass cylindrical tanks of 300 litres water capacity, and supplied with continuous flow of filtered seawater. After 2 weeks of adaptation to the experimental conditions, twelve groups of twenty-one gilthead sea bream (Sparus aurata) juveniles with an initial body weight of 70·8 (sd 0·03) g were established and the experimental diets were randomly assigned to triplicate groups of these fish.
The growth trial lasted for 81 d, and during this period fish were hand-fed twice a day, 6 d/week, to apparent visual satiety. At the end of the trial, fish were unfed for 1 d to empty the gut content and then bulk-weighed after mild anesthesia with 0·3 ml/l methylethanol. At the beginning and at the end of the growth trial, fifteen fish from the stock population and three fish from each tank were, respectively, sampled, pooled and frozen until whole-body composition analysis. During the trial, salinity averaged 34·7 (sd 0·8) g/l, dissolved O2 was kept near saturation and water temperature was regulated to 24·0 (sd 0·5)°C.
To eliminate handling stress, after the growth trial fish, continued to be fed for one more week, and then 18 h after the last meal (the previous day afternoon meal) nine fish from each tank were randomly sampled for blood, liver, intestine and muscle collection. Blood was collected from the caudal vein using heparinised syringes and centrifuged at 2500 g for 10 min and the recovered plasma was kept at –20°C until analysis. Thereafter, fish were killed with a sharp blow to the head, and whole body, viscera and liver were weighed for determination of hepatosomatic (HSI) and viscerosomatic (VSI) indexes. Liver, intestine and muscle sections were frozen in liquid N2 and then stored at –80°C until biochemical, enzymatic and molecular analyses.
Diets, whole fish, liver, muscle and plasma analysis
Chemical analysis of the diets, whole fish, liver and muscle was conducted according to AOAC( 52 ) and by following the procedures given below: DM after drying at 105°C until constant weight; ash by incineration in a muffle furnace at 450°C for 16 h; protein content (N×6·25) by the Kjeldahl method after acid digestion using a Kjeltec digestion and distillation unit (models 1015 and 1026, Tecator Systems; Höganäs); and lipid by petroleum diethyl ether extraction (Soxtec HT System; Höganäs). Starch was determined according to Beutler( Reference Beutler 53 ) and gross energy by direct combustion in an adiabatic bomb calorimeter (PARR model 6200; Parr Instruments).
Hepatic and muscular glycogen contents were determined as described by Roehrig & Allred( Reference Roehrig and Allred 54 ), and lipids were determined according to the method of Folch et al.( Reference Folch, Lees and Sloane Stanley 55 ). FA methyl esters were prepared by acid transmethylation of total lipids using boron trifluoride in methanol (14 %), as described by Shantha & Ackman( Reference Shantha and Ackman 56 ), and analysed by GC (Varian 3900; Varian; for details see Castro et al.( Reference Castro, Corraze and Panserat 48 )). Total cholesterol in diets, liver and muscle was assayed on total lipid extract by the Liebermann–Burchard method( Reference Stadtman 57 ). Plasma metabolites were analysed using commercial kits from Spinreact: glucose (ref: 1001191), TAG (ref: 1001312), total cholesterol (ref: 1001090) and PL (ref: 1001140).
Enzymatic activity assays
The activity of key lipogenesis enzymes was determined in the liver (three fish per tank). For that purpose, liver was homogenised (dilution 1:4) in ice-cold buffer (100 mm-Tris-HCl, 0·1 mm-EDTA and 0·1 % triton X-100 (v/v), pH 7·8). All procedures were performed on ice. Homogenates were centrifuged at 30 000 g for 30 min at 4°C. After centrifugation, the resultant supernatant was collected and aliquots were stored at –80°C until analysis. All enzyme activities were measured at 37°C, by monitoring the changes in absorbance of NADPH at 340 nm in a Multiskan GO microplate reader (model 5111 9200; Thermo Scientific), using 6·22 mm/cm as the millimolar extinction coefficient for NADPH. The optimal substrate and protein concentrations for measurement of each enzyme activity were established by preliminary assays. G6PD (EC 1.1.1.49), ME (EC 1.1.1.40) and FAS (EC 2.3.1.38) activities were determined as previously described by Castro et al.( Reference Castro, Corraze and Pérez-Jiménez 7 ).
Gene expression
Analyses of mRNA levels were performed on liver and intestine samples (two fish per tank). Tissues for RNA analyses were homogenised in 2-ml tubes containing Trizol reagent (Invitrogen) using rapid vibration (liver: 2×10 s, with an interval of 10 s, at 5000 rpm; intestine: 3×10 s, with 10-s intervals, at 6500 rpm) in Precellys®24 (Bertin Technologies). Extraction of total RNA was then performed according to the manufacturer’s recommendations. RNA quality and quantity were assessed by gel electrophoresis and spectrophotometry (NanoDrop ND1000; NanoDrop Labtech). Complementary DNA (cDNA) synthesis was performed with 1 µg of the resulting total RNA using SuperScript III RNaseH-Reverse Transcriptase kit (Invitrogen) and random primers (Promega). Real-time quantitative PCR (q-PCR) analyses were performed in a total volume of 6 µl (detailed information of the reaction mix in Castro et al.( Reference Castro, Corraze and Pérez-Jiménez 7 )) using LightCycler® 480 II apparatus (Roche Diagnostics) to assess the gene expression levels. Primers were either obtained in the literature (Mininni et al.( Reference Mininni, Milan and Ferraresso 44 ); Pérez-Sánchez et al.( Reference Pérez-Sánchez, Borrel and Bermejo-Nogales 45 ); Sánchez-Gurmaches et al.( Reference Sánchez-Gurmaches, Cruz-Garcia and Ibarz 46 ); Enes et al.( Reference Enes, Panserat and Kaushik 59 ); Diez et al.( Reference Diez, Menoyo and Perez-Benavente 74 )) or designed from gilthead sea bream-expressed sequence tag sequences available on the SIGENAE database (http://www.sigenae.org) using the Primer3 software( Reference Rozen and Skaletsky 58 ) (online Supplementary Table S1). For gene targets that had not been previously validated, primers were tested on a pool of cDNA and amplified products were systematically sequenced. The PCR protocol followed the conditions described previously by Castro et al.( Reference Castro, Corraze and Pérez-Jiménez 7 ). Each PCR run included duplicates of reverse transcription for each sample and negative controls (RT-free samples, RNA-free samples). PCR run for reference gene included quadruplicates for each sample (duplicates of reverse transcription and PCR amplification, respectively) and negative controls. Quantification of the target gene transcripts in the liver and intestine was done using β-actin gene expression as reference, as previously used in gilthead sea bream by Pérez-Sánchez et al.( Reference Pérez-Sánchez, Borrel and Bermejo-Nogales 45 ), and that was stably expressed in the present study (data not shown). Relative quantification of the target gene transcript with the β-actin reference gene transcript was performed using the mathematical model described by Pfaffl( Reference Pfaffl 60 ). The relative expression ratio (R) of a target gene was calculated on the basis of real-time PCR efficiency (E) and the threshold cycle (Ct) deviation (ΔCT) of the unknown sample compared with a control sample and expressed in comparison with the β-actin reference gene:
Efficiency of q-PCR was measured by the slope of a standard curve using serial dilutions of cDNA. The fish fed the FOCH – diet was used as the control group.
Statistical analysis
Data were checked for normality and homogeneity of variances, and they were normalised when appropriate. Statistical evaluation of data was carried out by a 2×2 factorial arrangement of treatments in a completely randomised experimental design (two-way ANOVA) with carbohydrate content and lipid source as fixed factors. The significance level of 0·05 was used for rejection of the null hypothesis. In cases in which interaction was significant, one-way ANOVA was performed for each factor. All statistical analyses were conducted using the SPSS 22.0 software package (IBM Corp.) for Windows.
Results
Dietary fatty acid composition
The four diets presented small differences in the proportions of total SFA, whereas MUFA were higher in FO diets, and n-3 an n-6 PUFA were higher in VO diets (Table 2). Within MUFA, high levels of oleic acid (18 : 1n-9) were recorded in VO diets, whereas the opposite occurred for palmitoleic acid (16 : 1n-7), eicosenoic acid (20 : 1n-9) and erucic acid (22 : 1n-9). Among n-3 PUFA, VO diets were particularly rich in linolenic acid (18 : 3n-3) and poor in EPA and DHA. The proportion of total n-6 PUFA was strongly higher in VO diets mainly because of linoleic acid (18 : 2n-6) levels.
FO, fish oil; VO, blend of vegetable oils; CH, carbohydrate; CH content, 0 % (CH–) or 20 % (CH+) gelatinised maize starch; n-3 LC-PUFA, n-3 long-chain PUFA; unsat. index, unsaturation index=sum (fatty acid percentage)×(number of double bonds).
Growth performance and feed utilisation
Fish promptly accepted the experimental diets, and no mortality was recorded during the trial. Dietary treatments had no effects on fish growth performance or feed utilisation (Table 3). Feed intake (g/kg average body weight per d) was similar among diets. N retention, expressed per unit weight gain, was not affected by diet composition. However, protein efficiency ratio was higher with the CH+ diets, which had a lower protein content. Dietary carbohydrate intake increased lipid retention per unit weight gain only when fish were fed the VO diet (Table 3). In fish fed CH– diets, lipid retention per unit weight gain was lower in the VO group. In addition, within the VO group, lipid retention per unit weight gain was lower in fish fed the CH– diet (carbohydrate and lipid source interaction).
LS, lipid source; FO, fish oil; VO, blend of vegetable oils; CH, carbohydrates; CH content, 0 % (CH–) or 20 % (CH+) gelatinised maize starch; IBW, initial body weight; FBW, final body weight; ABW, average body weight; FE, feed efficiency; PER, protein efficiency ratio.
a,b,A,B If interaction was significant, one-way ANOVA was performed for each factor, and means in the same line with different capital and small letters indicate significant differences (P<0·05) between the two tested LP and two CH levels, respectively; means with no letters are not significantly different (P>0·05).
* Significant differences at P<0·05 two-way ANOVA.
† DGI: ((FBW1/3−IBW1/3)/time (d))×100.
‡ ABW: (IBW+FBW)/2.
§ FE: wet weight gain/dry feed intake.
|| PER: wet weight gain/crude protein intake.
¶ N retention=((FBW×carcass N content)−(IBW×carcass N content))/(ABW×the number of days).
** Lipid retention=((FBW×carcass lipid content)−(IBW×carcass lipid content))/(ABW×the number of days).
Whole-body, liver and muscle composition
At the end of the trial, only whole-body lipid and DM contents were affected by dietary treatments. Under a VO-based diet regimen, whole-body lipid content was higher in CH+ than in CH– groups. In addition, whole-body lipid was lower in VO than in FO groups only when fish were fed no-carbohydrate diets (CH– diets) (carbohydrate and lipid source interaction).Whole-body DM content was higher in fish fed the CH+ diets (Table 4).
LS, lipid source; FO, fish oil; VO, blend of vegetable oils; CH, carbohydrates; CH content, 0 % (CH–) or 20 % (CH+) gelatinised maize starch.
a,b,A,B If interaction was significant, one-way ANOVA was performed for each factor, and means in the same line with different capital and small letters indicate significant differences (P<0·05) between the two tested LS and two CH levels, respectively; means with no letters are not significantly different (P>0·05).
* n 3 For whole-body composition; n 6 for lipids and cholesterol; n 9 for glycogen; n 18 for HSI and VSI.
† Significant differences at P<0·05 two-way ANOVA.
‡ Initial body composition on the fish: DM 28·96 %; protein 16·15 %; lipid 7·87 %; ash 5·84 %.
§ HSI: (liver weight/body weight)×100.
║ VSI: (viscera weight/body weight)×100.
Higher HSI and VSI were observed in the CH+ groups, but no effect of dietary lipid source was noticed (Table 4).
Lipid content in muscle was higher in the CH+ group, whereas in the liver an increase in lipid content with carbohydrate intake was only evident when fish were fed VO-based diets (carbohydrate and lipid source interaction). Carbohydrate intake also increased the glycogen content in liver, but in muscle a similar effect was only recorded when fish were fed FO-based diets (carbohydrate and lipid source interaction). Dietary lipid source induced no changes on hepatic and muscular lipid content, and glycogen content in both tissues increased in the VO group when fish were fed CH– diets (carbohydrate and lipid source interaction). Hepatic cholesterol content was lower in CH+ groups, but no differences were recorded in muscle cholesterol content. Dietary lipid source did not affect hepatic and muscular cholesterol content.
Liver and muscle fatty acid profiles
Muscle and liver FA composition were affected by diet composition and resembled the FA composition of the dietary lipid sources (Tables 5 and 6).
LS, lipid source; FO, fish oil; VO, blend of vegetable oils; CH, carbohydrates; CH content, 0 % (CH–) or 20 % (CH+) gelatinised maize starch; n-3 LC-PUFA, n-3 long-chain PUFA; unsat. index, unsaturation index=sum (fatty acid percentage)×(number of double bonds).
a,b,A,B If interaction was significant, one-way ANOVA was performed for each factor, and means in the same line with different capital and small letters indicate significant differences (P<0·05) between the two tested LS and two CH levels, respectively; means with no letters are not significantly different (P>0·05).
* FA≥0·02 %, <0·02 % was not considered in the table.
† Significant differences at P<0·05 two-way ANOVA.
LS, lipid source; FO, fish oil; VO, blend of vegetable oils; CH, carbohydrates; CH content, 0 % (CH–) or 20 % (CH+) gelatinised maize starch; n-3 LC-PUFA, n-3 long-chain PUFA; unsat. index, unsaturation index=sum (fatty acid percentage)×(number of double bonds).
a,b,A,B If interaction was significant, one-way ANOVA was performed for each factor, and means in the same line with different capital and small letters indicate significant differences (P<0·05) between the two tested LS and two CH levels, respectively; means with no letters are not significantly different (P>0·05).
* FA≥0·02 %, <0·02 % was not considered in the table.
† Significant differences at P<0·05 two-way ANOVA.
Except for muscle SFA content that was similar among CH+ and CH– diets, liver and muscle FA profiles of fish fed the CH+ diets were characterised by a higher proportion of SFA (particularly 16 : 0) and MUFA (particularly 18 : 1n-9) and a lower proportion of n-3 PUFA and n-6 PUFA. In the liver, replacing FO by VO resulted in higher proportions of SFA (mainly 16 : 0) and n-6 PUFA (mainly 18 : 2n-6) and a lower proportion of n-3 PUFA (mainly 22 : 6n-3, 20 : 5n-3, 22 : 5n-3) and MUFA (mainly 16 : 1n-7, 20 : 1n-9, 22 : 1n-9).
The muscle of fish fed the VO diets presented a lower proportion of SFA (mainly 14 : 0) and MUFA (mainly 16 : 1n-7, 20 : 1n-9, 22 : 1n-9) and an increased proportion of n-6 (mainly 18 : 2n-6) and n-3 (mainly 18 : 3n-3) PUFA.
Plasma metabolites and enzyme activity
After 18 h of feeding diets, CH+ plasma glucose concentration was lower than in fish fed CH– diets (Table 7).
LS, lipid source; FO, fish oil; VO, blend of vegetable oils; CH, carbohydrates; CH content, 0 % (CH–) or 20 % (CH+) gelatinised maize starch; GLU, glucose; CHOL, total cholesterol; PL, phospholipids; G6PD, glucose-6-phosphate dehydrogenase; ME, malic enzyme; FAS, fatty acid synthase.
a,b,A,B If interaction was significant, one-way ANOVA was performed for each factor, and means in the same line with different capital and small letters indicate significant differences (P<0·05) between the two tested LS and two CH levels, respectively; means with no letters are not significantly different (P>0·05).
* Significant differences at P<0·05 two-way ANOVA.
Plasma cholesterol was also lower in fish fed the CH+ diet when combined with VO (carbohydrate and lipid source interaction). On the contrary, plasma TAG was higher in fish fed the CH+ diet, but only in fish fed FO (carbohydrate and lipid source interaction). Plasma PL and cholesterol were lower in fish fed the VO diets. Among fish fed the CH+ diets, plasma TAG was lower in the VO group than in the FO group (carbohydrate and lipid source interaction).
The CH+ diet promoted an increase in FAS, G6PD and ME activities. G6PD activity was also responsive to dietary lipid source, being higher in VO diets (Table 7).
Gene expression
Hepatic transcript levels of glucokinase (GK), the first glycolytic enzyme, was higher in CH+ than in CH– groups, whereas 6-phosphofructokinase (PFK) transcript levels were not affected by dietary carbohydrate (Fig. 1). Hepatic transcript levels of phosphoenolpyruvate carboxykinase (PEPCK), the first key enzyme involved in gluconeogenesis, were lower in CH+ than in CH– groups, whereas the opposite occurred for hepatic glucose 6-phosphatase (G6Pase) mRNA levels, the enzyme involved in the last step of gluconeogenesis. Among glycolytic and gluconeogenic enzymes, only hepatic G6Pase transcript levels were up-regulated by dietary VO. At the intestinal level, no transcriptional regulation by diet composition of PFK and of PEPCK was observed (Fig. 1).
Hepatic and intestinal transcript levels of key enzymes involved in lipogenesis (FAS), β-oxidation (carnitine palmitoyltransferase 1A (CPT1A)) and PL synthesis (glycerol-3-phosphate acyltransferase (GPAT)) were not affected by diet composition (Fig. 2).
With the exception of lanosterol 14 α-demethylase, cytochrome P450 51 (CYP51A1) transcript levels in the intestine, which were higher in CH+ groups, no variation in hepatic or intestine transcript levels of proteins involved in cholesterol efflux (ATP-binding cassette G5 (ABCG5)), and cholesterol synthesis (7-dehydrocholesterol reductase (DHCR7); CYP51A1) and catabolism (liver X receptor α (LXRα)) were observed (Fig. 3).
The expression of genes encoding key proteins involved in the LC-PUFA-biosynthesis pathway (∆6 fatty acyl desaturase-liver isoform (FADS2); elongase 5 (elovl5)) were nutritionally regulated both in the liver and in the intestine (Fig. 4). In the liver, FADS2 transcript levels were upregulated by dietary carbohydrate and VO. Elovl5 expression in the liver was down-regulated in CH+ groups, but in the intestine down-regulation was only evident in fish fed FO-based diets (carbohydrate and lipid source interaction). In addition, intestinal elovl5 mRNA levels increased in VO groups when fed a carbohydrate-rich diet (carbohydrate and lipid source interaction).
In the liver transcription of PPARα and PPARγ was down-regulated in CH+ groups while in the intestine down-regulation of PPARβ transcript levels in CH+ groups was only observed when fish were fed the VO diet (carbohydrate and lipid source interaction) (Fig. 4). Additionally, intestine from the CH– groups exhibited an up-regulation of PPARβ transcripts when fish were fed the VO diet (carbohydrate and lipid source interaction).
Discussion
Effects of dietary carbohydrates
Although gilthead sea bream is a carnivorous species (trophic level 3·3–3·5 according to Fish Base), it tolerates up to 20 % dietary starch without detrimental effects in growth performance and feed efficiency. These results support previous evidences( Reference Couto, Enes and Peres 61 , Reference Enes, Panserat and Kaushik 62 ) and indicate the possibility to reduce aquafeed costs and alleviate the overexploitation of fisheries marine resources through the use of the carbohydrate component in gilthead sea bream diets.
The low glycaemia values at 18 h after feeding are in accordance with results of Peres et al. ( Reference Peres, Gonçalves and Oliva-Teles 63 ) that observed in a glucose tolerance test that seabream was able to restore glucose levels within 12 h after receiving an overdose of glucose.
Hepatic GK is a key player in blood glucose homoeostasis by catalysing the phosphorylation of glucose and providing the first substrate for glycolysis, glycogenesis and the pentose phosphate pathway( Reference Engelking 64 ). An induction of hepatic GK transcripts by dietary carbohydrates was recorded in the present study, as in a previous study in this species( Reference Panserat, Medale and Blin 65 ). On the contrary, dietary carbohydrates did not affect transcription levels of hepatic PFK, another key glycolytic enzyme.
Hepatic transcript levels of PEPCK, the first key enzyme of gluconeogenesis, were down-regulated by dietary carbohydrate, but the transcript levels of G6Pase, another key enzyme of gluconeogenesis, was up-regulated by dietary carbohydrates. Although a similar response of hepatic G6Pase by dietary carbohydrate at the transcriptional level was previously reported in carnivorous rainbow trout( Reference Kamalam, Médale and Kaushik 24 ), present data apparently contradict results of Panserat et al.( Reference Panserat, Plagnes-Juan and Kaushik 66 ) in this species, who observed that G6Pase and FBPase were down-regulated by dietary carbohydrates at the transcriptional level, whereas PEPCK was not affected. At the enzymic activity level, carbohydrates in gilthead sea bream diets induced minor effects( Reference Enes, Panserat and Kaushik 62 , Reference Fernández, Miquel and Córdoba 67 ) or even an increase( Reference Couto, Enes and Peres 61 ) in the hepatic activities of gluconeogenesis enzymes. Recent studies pointed out that the intestine has important functions in glucose homoeostasis( Reference Kamalam, Panserat and Aguirre 68 , Reference Polakof, Álvarez and Soengas 69 ). However, present results do not support that assumption, as expression of glycolytic (PFK) and gluconeogenesis (PEPCK) enzymes in this tissue did not respond to dietary carbohydrate.
In parallel with the hepatic up-regulation of GK transcripts, increased HSI and VSI were recorded in fish fed carbohydrate-rich diets, which may indicate that the hepatic glucose pool was directed towards glycogen and/or FA synthesis. Indeed, in liver, an increased deposition of glycogen and a higher lipogenic potential, indicated by FAS, G6PD and ME activities, were recorded in CH+ groups. Increased HSI may have unwanted physiological effects; therefore, histomorphological analysis of liver should be considered in future studies to discard any histopathological damage induced by diet.
Despite the fact that lipogenesis was nutritionally regulated by dietary carbohydrates, the contribution to the overall lipid deposition in liver or whole body of fish fed FO-based diet seemed to be minor, as lipids in liver and whole body increased with carbohydrate intake only under a VO-based diet regimen.
Recently, it was suggested that dietary carbohydrates have a role in cholesterol biosynthesis by inducing (at least at a molecular level) the capacity to produce it( Reference Castro, Corraze and Pérez-Jiménez 7 ). In the present study, it was observed that gilthead sea bream fed carbohydrate-rich diets exhibited reduced liver and plasma (only in VOCH+ group) cholesterol content. Furthermore, an increased expression at the intestinal level of one gene encoding for an enzyme involved in the last steps of cholesterol biosynthesis, CYP51A1, was also observed. We hypothesise that the depletion of the cholesterol pool was promoted by a lower dietary input of cholesterol, and the up-regulation of CYP51A1 may reflect an increased synthesis of sterols by the intestinal cells in response to the low dietary supply. Metabolic adjustments in de novo cholesterol biosynthesis in response to the dietary load of cholesterol were previously described in Atlantic salmon( Reference Kortner, Björkhem and Krasnov 70 ).
In this work, we demonstrated a transcriptional up-regulation by dietary carbohydrate of liver FADS2, a key desaturase involved in LC-PUFA synthesis. With the up-regulation of FADS2 gene, it was expected to observe increased transcript levels of elovl5, a fatty acyl elongase involved in LC-PUFA biosynthesis. However, and somehow surprisingly, the reverse pattern was observed both in the intestine (only in FOCH+ group) and in the liver of gilthead sea bream fed CH+ diets. In the present state of knowledge, we have no clear explanation for these conflicting results, and therefore further studies are needed to better understand the apparent atypical molecular regulation of some enzymes of the LC-PUFA synthesis pathway by dietary carbohydrate in gilthead sea bream. Irrespective of the regulation of FADS2 expression, increased liver or muscle n-3 LC-PUFA content of gilthead sea bream fed CH+ diets were not observed. On the contrary, reduced n-3 LC-PUFA content in CH+ groups was recorded, which could be related to an increase in SFA derived by lipogenesis from carbohydrates, as previously reported in other species such as rainbow trout and European sea bass( Reference Castro, Corraze and Pérez-Jiménez 7 , Reference Álvarez, Lopez-Bote and Diez 71 ).
PPAR are a family of nuclear receptors that have three isoforms in mammals known as PPARα, PPARβ and PPARγ that play key roles in regulation of lipid metabolism( Reference Viana Abranches, Esteves de Oliveira and Bressan 72 ). It is believed that PPAR control metabolic pathways that support FA β-oxidation (especially PPARα), tissue lipid deposition and lipogenesis (especially PPARγ), and the overall lipid homoeostasis (especially PPARβ), mainly by means of the action of its ligands (such as FA and their derivatives, and also glucose)( Reference Viana Abranches, Esteves de Oliveira and Bressan 72 , Reference Leaver, Bautista and Bjornsson 73 ). In this study, despite the down-regulation of PPARα,β gene expression in the liver and of PPARβ in the intestine (only in VO diets) by dietary carbohydrates, gene expression of CPT1, a marker of mitochondrial FA β-oxidation, was unaffected both in the liver and intestine of fish fed CH+ diets. In previous studies, a synchronised regulation between piscine PPARα,β gene expression and the activity or expression of other enzymes involved in FA β-oxidation, such as acyl-CoA oxidase and L-3-hydroxyacyl-CoA dehydrogenase, was described( Reference Diez, Menoyo and Perez-Benavente 74 , Reference Du, Demizieux and Degrace 75 ). However, present results are not completely unexpected, as in a number of studies the expression of the CPT1 gene was not nutritionally regulated or directly linked to the FA β-oxidation capacity( Reference Kamalam, Médale and Kaushik 24 , Reference Morais, Pratoomyot and Torstensen 37 , Reference Kennedy, Leavera and Campbell 76 ). Furthermore, FA β-oxidation and lipogenesis are two pathways usually regulated in opposite directions( Reference Bonacic, Estevez and Bellot 77 , Reference Zheng, Luo and Zhuo 78 ). Therefore, as lipogenic potential increased in liver of fish fed CH+ diets, it is possible that FA β-oxidation might have been repressed. Additional data on the activity or expression of the other FA β-oxidation related enzymes would be necessary to confirm this hypothesis.
In the above discussion, differences in metabolic responses were assumed to be related to differences in dietary carbohydrate (0 v. 17–18 %). However, a potential effect of dietary protein cannot be discarded, particularly in relation to regulation of gluconeogenesis and lipogenesis, as there is increasing evidence that protein and amino acids also modulate these pathways( Reference Dai, Panserat and Plagnes-Juan 79 – Reference Lansard, Panserat and Plagnes-Juan 82 ).
Effects of dietary lipid source
In the present study, gilthead sea bream performed as well as in previous studies with gilthead sea bream juveniles in which dietary FO was partially (70 %) replaced by VO in FM( Reference Izquierdo, Montero and Robaina 83 , Reference Montero, Mathlouthi and Tort 84 ) or plant protein( Reference Benedito-Palos, Saera-Vila and Calduch-Giner 85 , Reference Fountoulaki, Vasilaki and Hurtado 86 )-based diets. However, plasma cholesterol and PL concentrations were lower in fish fed the VO diets. Replacing FO with VO usually increases dietary phytosterol and reduces cholesterol content. Such dietary modifications have been reported to induce a decrease in plasma cholesterol and LDL-cholesterol both in humans( Reference Brufau, Canela and Rafecas 87 , Reference MacKay and Jones 88 ) and in fish( Reference Richard, Kaushik and Larroquet 35 – Reference Morais, Pratoomyot and Torstensen 37 , Reference Gilman, Leusch and Breckenridge 89 ).
Phytosterols are structurally similar to cholesterol and may induce a relative cholesterol deficiency by mechanisms that interfere with cholesterol absorption such as competition for space in mixed micelles or competition with cholesterol transporters (such as ABC transporters)( Reference Brufau, Canela and Rafecas 87 , Reference Ostlund 90 ). For instance, in Atlantic salmon, dietary inclusion of soyabean meal and soya saponins decreased the expression of ABCG5 in the distal intestine( Reference Kortner, Gu and Krogdahl 91 ) and in the liver( Reference Gu, Kortner and Penn 92 ). In the present study, however, no repression of ABCG5 expression related to diet composition was observed at hepatic or intestinal levels. Furthermore, no induction in the expression of the selected genes (CYP51A1 and DHCR7) involved in the cholesterol biosynthesis was recorded in the intestine and liver of fish fed the VO diets. This differs from other studies in which hepatic or intestinal gene expression of key enzymes involved in cholesterol biosynthesis pathways (such as 3-hydroxy-3-methylglutaryl-CoA reductase, isopentenyl-diphosphate Δisomerase, CYP51A1 or DHCR7) increased in response to the decreased plasma cholesterol concentration in fish fed plant-based( Reference Kortner, Gu and Krogdahl 91 – Reference Geay, Ferraresso and Zambonino-Infante 93 ) or VO( Reference Castro, Corraze and Pérez-Jiménez 7 , Reference Leaver, Villeneuve and Obach 10 , Reference Liland, Espe and Rosenlund 39 ) diets. Whether the lack of activation of cholesterol biosynthesis at the molecular level as a response to the reduced plasma cholesterol concentration that was observed in the present study is due to variations in the dietary cholesterol supply or due to species specific sensitivities to alternative ingredients, such as VO, remains to be clarified.
An intestinal and hepatic induction of expression or increased activity of enzymes involved in PL biosynthesis has been reported in carnivorous fish fed VO( Reference Caballero, Gallardo and Robaina 38 ) or plant feedstuff-based diets( Reference Geay, Ferraresso and Zambonino-Infante 93 ). However, the biosynthesis rate seems to be insufficient to avoid an accumulation of lipid droplets in the enterocytes( Reference Caballero, Izquierdo and Kjorsvik 94 , Reference Olsen, Dragnes and Myklebust 95 ). This condition is usually associated with a PL deficit that promotes an impairment in lipoprotein assembly and export from enterocytes and reduces PL concentration in plasma( Reference Castro, Corraze and Pérez-Jiménez 7 , Reference Richard, Kaushik and Larroquet 35 , Reference Richard, Mourente and Kaushik 36 ).
In the present study, despite the lower plasma PL concentration in VO groups, a nutritional regulation of transcript levels of GPAT was not observed in the intestine or the in liver of gilthead sea bream fed the experimental diets. GPAT is involved in the first steps of glycerol-3-phosphate pathway, which according to Caballero et al.( Reference Caballero, Gallardo and Robaina 38 ) is the main pathway implicated in intestinal PL synthesis in gilthead sea bream. Therefore, further investigation on the nutritional regulation of the enzymes involved in downstream steps of the PL biosynthesis pathways is required to understand PL biosynthesis regulation in gilthead sea bream.
Dietary lipid source is known to regulate lipogenesis and FA bioconversion pathways, which may affect tissue lipid deposition and FA composition. Although in some studies increased hepatic lipogenic enzymes activity and lipid content have been reported in fish fed diets in which dietary FO was replaced with VO( Reference Betancor, Sprague and Usher 6 , Reference Peng, Xu and Mai 19 ), in other studies such effects were not demonstrated( Reference Castro, Corraze and Pérez-Jiménez 7 , Reference Richard, Kaushik and Larroquet 35 , Reference Richard, Mourente and Kaushik 36 , Reference Torstensen, Froyland and Lie 96 ). In the present study, dietary lipid source did not affect hepatic and muscular lipid content, and a decrease in whole-body lipid content was reported within CH– groups when FO was substituted by VO. On the other hand, liver and muscle FA composition was strongly influenced by dietary replacement of FO by VO, as previously reported in gilthead sea bream( Reference Benedito-Palos, Bermejo-Nogales and Karampatos 97 ) and in other species( Reference Castro, Corraze and Pérez-Jiménez 7 , Reference Jordal, Lie and Torstensen 8 , Reference Tocher, Zheng and Schlechtriem 29 , Reference Benitez-Dorta, Caballero and Izquierdo 98 ). According to the predicted effects of dietary oils in muscle FA composition of 1-year-old gilthead sea bream with different nutritional backgrounds( Reference Benedito-Palos, Bermejo-Nogales and Karampatos 97 ), in the present study liver and muscle lipids increased the content of 18 : 1n-9, 18 : 2n-6 and 18 : 3n-3, whereas total lipid reduced the content of 20 : 5n-3 and 22 : 6n-3 when FO was replaced by VO. Thus, despite the hepatic up-regulation of FADS2 in fish fed the VO-based diets, the low proportion of LC-PUFA or their intermediary products such as 18 : 3n-6 and 18 : 4n-3 in VO groups suggest a limited FADS2 activity. Maximal FADS2 efficiency is thought to be modulated by the levels of substrate and product availability( Reference Vagner and Santigosa 30 ). In this sense, in the present study the limited accumulation of FA intermediates could partly have been caused by an inadequate dietary supply of substrate and/or product availability. Accordingly, in a previous study with this species using a diet completely devoid of LC-PUFA and containing olive oil as the sole lipid source, Seiliez et al.( Reference Seiliez, Panserat and Corraze 47 ) observed higher transcript levels of FADS2 and increased accumulation of FA intermediates (18 : 2n-9, 20 : 2n-9 and 18 : 3n-6) in fish fed that diet comparatively to fish fed an LC-PUFA-rich diet, which suggested the existence of FADS2 activity.
In this study, the low efficacy in induction of FADS2 expression and of nutritional regulation of elovl5 transcripts by VO observed, along with lack of ∆5 desaturase activity observed in vitro ( Reference Tocher and Ghioni 99 ), may explain the low capacity of conversion of C18 PUFA into LC-PUFA at an appreciable rate in gilthead sea bream.
Dietary carbohydrate content and lipid source interaction
The interactions between dietary carbohydrate content and lipid source that were recorded on the whole-body, liver and plasma lipid content suggest that the overall effect of starch intake on lipid deposition was disturbed by the change in dietary lipid source. Indeed, an increase in lipid retention with dietary starch intake was more evident when fish were fed the VO-based diet. This assumption is supported by the higher, whole-body and liver lipid content in the VOCH+ group than in the VOCH– group and can be related with the coupled increase in the lipogenic potential of dietary carbohydrates (higher liver FAS, ME and G6PD activities) and VO (higher liver G6PD activity). On the other hand, under an FO-based diet regimen, the unaltered lipid retention, liver and whole-body lipid content, but elevated TAG levels following starch intake, suggest a higher lipid mobilisation, transport and/or utilisation in the FOCH+ group than in the FOCH– group. No molecular markers for lipid uptake were assessed in this study in liver and intestine, and the one evaluated here and related to catabolism (CPT1A) does not help to understand or clarify this hypothesis because transcriptional regulation of these proteins by nutrients was not recorded.
At the same time, we also noticed that the stimulatory effect of dietary starch on muscle glycogen deposition was attenuated and became non-significant under the VO-based dietary regimen, probably related with the fact that dietary VO per se, but not in combination with carbohydrates, had an inductor effect on glycogenesis, as suggested by the higher muscle and liver glycogen levels in the VOCH– group than in the FOCH– group.
Furthermore, it was observed that intake of starch when coupled with dietary VO seemed to enhance the hypocholesterolemic effects of dietary VO as plasma cholesterol concentration decreased with carbohydrate intake only when fish were fed VO-based diet. It is important to note here that, within the molecular actors involved in cholesterol-related processes (CYP51A1, DHCR7, ABCG5, LXR α) assessed in liver and intestine in the present study, no molecular difference induced by a nutritional interaction of nutrients that could reflect this phenotype was found. Therefore, the physiological or metabolic mechanisms underlying this finding remain to be demonstrated.
Conclusion
To our knowledge, this is the first study in a marine fish species reporting a transcriptional induction of FADS2 by dietary carbohydrate, in addition to VO. Although the n-3 LC-PUFA biosynthesis pathway was insufficient to compensate for the reduced dietary EPA and DHA in VO-based diets, this study provides new perspectives on the use of nutritional strategies for inducing LC-PUFA biosynthesis in marine fish species.
Furthermore, change in dietary lipid source seemed to modify the overall effect of starch intake on mechanisms involved in cholesterol body pools, lipid and glycogen body allocation. Considering that these metabolites greatly influence the fish quality, the present findings highlight the critical need to assess the potential effects between dietary nutrients on metabolic-related processes involved in tissue metabolites deposition, especially in the context of alternative aquafeeds rich in VO and carbohydrates.
Acknowledgements
The authors express their thanks to P. Correia and M. Cluzeaud for their technical assistance.
This work was partially supported by national funds through Foundation for Science and Technology (FCT) – under the project ‘PEst-C/MAR/LA0015/2011’ and co-financed by the European Regional Development Fund through the COMPETE Operational Competitiveness Programme. C. C. was supported by a grant (SFRH/BD/76297/2011) from FCT, Portugal, and A. D. was supported by the National Council of Technological and Scientific Development (CNPq), Brazil.
C. C. carried out the main experimental work and wrote the draft of the manuscript under the direction of the project designer and leaders A. O.-T., G. C. and S. P.; A. D. assisted with the biochemical analyses; and L. L. performed the fatty acid analyses. All authors contributed to and approved the manuscript.
The authors declare that there are no conflicts of interest.
Supplementary material
For supplementary material/s referred to in this article, please visit http://dx.doi.org/10.1017/S000711451600163X