Adipose tissue stearoyl-CoA desaturase 1 index is increased and linoleic acid is decreased in obesity-prone rats fed a high-fat diet

Background Fatty acid (FA) composition and desaturase indices are associated with obesity and related metabolic conditions. However, it is unclear to what extent desaturase activity in different lipid fractions contribute to obesity susceptibility. Our aim was to test whether desaturase activity and FA composition are linked to an obese phenotype in rats that are either obesity prone (OP) or resistant (OR) on a high-fat diet (HFD). Methods Two groups of Sprague–Dawley rats were given ad libitum (AL-HFD) or calorically restricted (HFD-paired; pair fed to calories consumed by chow-fed rats) access to a HFD. The AL-HFD group was categorized into OP and OR sub-groups based on weight gain over 5 weeks. Five different lipid fractions were examined in OP and OR rats with regard to proportions of essential and very long-chain polyunsaturated FAs: linoleic acid (LA), alpha-linolenic acid, eicosapentaenoic acid, docosahexaenoic acid and the stearoyl-CoA desaturase 1 (SCD-1) product 16:1n-7. FA ratios were used to estimate activities of the delta-5-desaturase (20:4n-6/20:3n-6), delta-6-desaturase (18:3n-6/18:2n-6), stearoyl-CoA desaturase 1 (SCD-1; 16:1n-7/16:0, SCD-16 and 18:1n-9/18:0, SCD-18), de novo lipogenesis (16:0/18:2n-6) and FA elongation (18:0/16:0). Fasting insulin, glucose, adiponectin and leptin concentrations were measured in plasma. Results After AL-HFD access, OP rats had a significantly higher SCD-16 index and 16:1n-7 proportion, but a significantly lower LA proportion, in subcutaneous adipose tissue (SAT) triacylglycerols, as well as significantly higher insulin and leptin concentrations, compared with OR rats. No differences were found between the two phenotypes in liver (phospholipids; triacylglycerols) or plasma (cholesterol esters; phospholipids) lipid fractions or for plasma glucose or adiponectin concentrations. For the desaturase indices of the HFD-paired rats, the only significant differences compared with the OP or OR rats were higher SCD-16 and SCD-18 indices in SAT triacylglycerols in OP compared with HFD-paired rats. Conclusion The higher SCD-16 may reflect higher SCD-1 activity in SAT, which in combination with lower LA proportions may reflect higher insulin resistance and changes in SAT independent of other lipid fractions. Whether a lower SCD-16 index protects against diet-induced obesity is an interesting possibility that warrants further investigation.


Background
Obesity is a serious health concern that is increasing at an unprecedented rate worldwide. Excess intake of highenergy, highly palatable food may in turn result in a decreased metabolic response to increased fat intake, promoting a vicious cycle of weight gain [1]. However, even when presented with a high-energy diet, not all individuals become obese [2]. Thus, some are seemingly resistant to an obesity-inducing environment. Understanding the underlying mechanisms of this phenomenon is one of the central questions in the study of food intake regulation [3]. Just as humans differ in their susceptibility to diet-induced weight gain, individual animals within certain rodent strains also show different responses to excess energy intake [4][5][6]. On an ad libitum high-fat diet (HFD), the outbred Sprague-Dawley strain of rats can be divided into obesity-prone (OP) and obesity-resistant (OR) phenotypes to model human diet-induced obesity (DIO) [7,8]. This two-phenotype model can be useful for elucidating mechanisms that drive DIO, and to identify physiological and biochemical differences between the two groups.
Compared with their OR counterparts, OP rats show a more adverse metabolic risk profile as well as lower energy expenditure on a HFD [9,10], and a proadipogenic expression profile [11][12][13]. In addition, OP rats also have lower relative fat oxidation and fat oxidation enzyme expression [14][15][16][17]. However, the mechanisms behind these observed differences in DIO are still unclear, and further studies are needed to assess other potential metabolic and biochemical factors of importance.
Although OP/OR Sprague-Dawley rats represent key models for the study of differential susceptibility to DIO, the contribution of desaturase activities and fatty acid (FA) composition in various lipid fractions to these phenotypes has not been evaluated. Some differences in FA composition have been found between OP and OR animals, but the data are limited to a few FAs and lipid fractions [13], and not all desaturase indices have been studied.
Accumulating data suggest that the FA profile in plasma and metabolic tissues associates strongly with obesity and insulin resistance [18,19]. An altered FA metabolism and desaturase activity have been shown to be involved in body fat storage and oxidation [20,21]. Certain FAs and their ratios can also be utilized as important risk markers for various diseases and desaturase indices have been closely linked to several obesity-related conditions [22][23][24][25][26]. Especially, SCD-16 and D6D have been found to correlate positively with obesity and development of the metabolic syndrome [19,24,27]. The desaturase indices are FA product-to-precursor ratios that have been used to indirectly calculate the activities of the enzymes responsible for desaturating FAs [28]. By inserting double bonds into the FA carbon chains, desaturases convert saturated into unsaturated FAs, which in turn are substrates for various lipids and participate in signal transduction and food intake regulation [29,30]. The three desaturases, Δ5-(D5D), Δ6-(D6D) and, especially, the Δ9-desaturase (also known as stearoyl-CoA-desaturase or SCD-1), have been extensively studied (for review see Bjermo & Risérus [31]). SCD-1 is increased by saturated FAs and depressed by polyunsaturated FAs. Since high SCD-1 activity is associated with decreased fat oxidation and increased FA synthesis, and since an SCD-1 knockout leads to resistance to DIO, SCD-1 has been proposed as a potential target in the treatment of obesity [31,32].
Given the role for FA composition and desaturases in pathophysiology and metabolism we aimed to investigate these factors in response to HFD-induced weight gain. Potential differences between the OP and OR phenotypes in desaturase indices and biochemically important FAs, were investigated in five lipid fractions involved in FA handling and metabolism: subcutaneous adipose tissue triacylglycerols (SAT-TG), liver phospholipids (liver PL), liver triacylglycerols (liver TG), plasma cholesterol esters (PL-CE) and plasma phospholipids (PL-PL). In order to infer the impact of diet versus metabolic differences, we chose to study three groups: an ad libitum HFD group of rats divided into OP and OR, and a calorically restricted HFD group as control.

Results
We used two dietary groups divided into a total of three experimental groups. One group of rats had ad libitum access to a HFD (AL-HFD, n = 24); a HFD-paired group (n = 10) were fed the HFD but were calorically restricted to follow the weight trajectory of male Sprague-Dawley rats fed regular chow. The rats in the AL-HFD group were further divided into either obesity-prone (OP, n = 12) or obesityresistant (OP, n = 12), based on a median split of the gain in body weight (dBw) on the HFD after the 5-week experiment. There was no significant difference in initial body weight between the three groups.

Desaturase indices in OP, OR and HFD-paired rats
Comparing the desaturase indices of OP and OR, we found a significant difference for SCD-16 in SAT-TG, with OP rats having a higher index than OR rats (P < 0.01) (Table 2; Figure 2). Comparing OP and OR with HFDpaired rats, significant differences were detected in SAT-TG, where SCD-16 was significantly higher in OP compared with HFD-paired rats (P < 0.05), but no significant difference was seen between OR and HFD-paired rats. SAT-TG SCD-18 was significantly higher in OP compared to HFD-paired rats (P < 0.01), but no difference was detected between OR and HFD-paired rats. For all other tested desaturase indices in SAT-TG and the four other lipid fractions (PL-CE, PL-PL, liver PL and liver TG), no significant differences were found between OP, OR and HFD-paired rats.

De novo lipogenesis and elongation indices in OP, OR and HFD-paired rats
In the initial ANOVA there was a significant difference in the de novo lipogenesis (DNL) index of SAT-TG, which did not however remain significant in the post-hoc test.
In PL-CE, PL-PL, liver PL and liver TG there were no significant differences for the DNL indices.
The elongation index (EI) of SAT-TG was significantly lower in OP (P < 0.01) compared with HFD-paired rats. In PL-PL, the EI was significantly higher in both OP (P < 0.001) and OR (P < 0.01) compared with HFD-paired rats. In PL-CE, liver PL and liver TG, no significant differences between the three groups were found for the EI. Results for all the studied indices can be found in Table 2.

Fatty acid proportions in OP, OR and HFD-paired rats
In SAT-TG, 18:2n-6 (linoleic acid, LA) was found to be significantly lower in OP compared with OR (P < 0.01) ( Figure 3) and HFD-paired rats (P < 0.05), but there was no significant difference in LA proportions between OR and HFD-paired rats. SAT-TG 18:3n-3 (alpha-linolenic acid, ALA) was significantly higher in both OP and OR (P < 0.01) compared to HFD-paired rats. The product of SCD-1, 16:1n-7, was significantly higher in OP compared with OR (P < 0.01) and HFD-paired (P < 0.05) rats. The proportion of 22:6n-3 (docosahexaenoic acid, DHA) was significantly higher in OP compared with HFD-paired (P < 0.05), but not OR, rats. All FA comparisons can be found in Table 2.
In PL-CE, ALA was significantly lower in both OP and OR compared with HFD-paired rats (P < 0.05). While there were no differences in the 20:5n-3 (eicosapentaenoic acid, EPA) proportion between OP and HFD-paired rats, OR rats had a significantly lower proportion than HFDpaired (P < 0.05). In PL-PL, DHA was significantly higher in both OP and OR compared with HFD-paired rats (P < 0.05). In liver PL, the only significant difference was seen for DHA, of which OR had a higher proportion than HFD-paired rats (P < 0.01), which did not however differ significantly from OP rats. We were not able to statistically analyze ALA in liver PL due to a lack of data. Finally, no significant differences were found for liver TG or for any other of the studied FAs. Body weight (BW) and food intake in obesity-prone (OP), obesity-resistant (OR) and HFD-paired rats. The symbol # indicates a significant difference between OP and HFD-paired; † between OP and OR; ‡ between OR and HFD-paired. Level of significance: # / † / ‡ , P < 0.05; ## / † † / ‡ ‡ , P < 0.01; ### / † † † / ‡ ‡ ‡ , P < 0.001. All tests were done using One-Way ANOVAs, followed by Tukey's t-test when significant; or Kruskal Wallis test (for non-normally distributed data; indicated by a superscript "a" where performed), followed by Dunn's multiple comparison test when significant. Figure 1 Food intake was correlated with dBw only in OP rats. Average daily food intake was highly correlated with gain in body weight (dBw) in obesity-prone (OP; r = 0.87, P < 0.001) but not in obesity-resistant rats (OR; r = 0.34, P = 0.28). Data analyzed using Spearman's rank correlation. Glucometabolic and adipokine levels in OP, OR and HFD-paired rats Plasma insulin and leptin concentrations were significantly higher in OP compared with OR and HFD-paired rats, but there were no significant differences between OR and HFD-paired rats. No significant differences were observed for fasting plasma glucose values. Plasma adiponectin concentrations were significantly higher in OP compared to HFD-paired, but not OR, rats. All glucometabolic and adipokine values are shown in Table 1.

Discussion
In this study we examined the desaturase indices and FA composition of two different phenotypes, OP and OR, in five metabolically important lipid fractions in an outbred rat model suitable for modeling human DIO. We discovered a significant difference between OP and OR in the SCD-1 index SCD-16 and the proportion of linoleic acid in subcutaneous adipose tissue (see Figure 4 for overview), paralleled by hormonal changes indicating an insulin resistant state. Interestingly, in other lipid fractions, no significant differences between the two phenotypes were detected. This indicates that adipose tissue is affected early in the development of obesity, and alterations in the FA Fatty acid (FA) composition and FA enzyme indices in five lipid fractions in obesity-prone (OP) and obesity-resistant (OR) rats with ad libitum HFD access (AL-HFD), and energy-restricted HFD-paired rats (HFD-paired). FAs are presented as relative percentages of total FAs within each lipid fraction. The symbol # indicates a significant difference between OP and HFD-paired; † between OP and OR; ‡ between OR and HFD-paired. Level of significance: # / † / ‡ , P < 0.05; ## / † † / ‡ ‡ , P < 0.01; ### / † † † / ‡ ‡ ‡ , P < 0.001. N/A, not available. All tests were done using One-Way ANOVAs, followed by Tukey's t-test when significant; or Kruskal Wallis test (for nonnormally distributed data; indicated by a superscript "a" where performed), followed by Dunn's multiple comparison test when significant. Where data for only two groups was available, Student's t-test was used (indicated by a superscript "b" where performed) for normally distributed data. HFD-induced SCD-16 activity calculated from the subcutaneous adipose tissue triacylglycerol (SAT-TG) ratio 16:1n-7/16:0. Obesityprone (OP) rats had a significantly higher ratio than obesity-resistant (OR; P < 0.01) and HFD-paired (P < 0.05) rats, even after removing the lowest value (outlier) in OR. Scatter plot lines with error bars indicate means ± SD; data analyzed with One-way ANOVA and Tukey's pos-hoc test.
composition of the adipose tissue may precede or be independent of those in the liver or plasma, as supported by the observed opposing changes in the elongation index of subcutaneous adipose tissue triacylglycerols and plasma phospholipids.
Our results help to further elucidate the role that SCD-1 plays in driving susceptibility to diet-induced obesity, as it has not yet been determined whether differences in SCD-1 between OP and OR rats exist in the adipose tissue and whether these changes occur simultaneously across lipid fractions. It has previously been found that a HFD, especially a diet rich in SFA, decreases SCD expression in both rat liver and adipose tissue [33,34]. A HFD has also been shown to decrease the index SCD-16 in adipose tissue in a sex-independent manner [35]. Another study has instead shown that SFA are essential for increasing SCD-1 in mice [36], perhaps by counteracting the effect of PUFAs [37]. Human intervention studies have also reported increased SCD-16 on a SFA-rich diet [38]. It has also been found that hepatic SCD-1 levels are higher in OP versus OR rats [13] and OP versus OR mouse strains [39]. These diverging results may be due to differences in e.g. study design, diet, animal species and tissue investigated. Figure 3 The proportion of linoleic acid in SAT-TG decreased in OP rats fed a HFD. HFD-induced linoleic acid (18:2n-6) proportions in subcutaneous adipose tissue triacylglycerols (SAT-TG) as relative percentages of total fatty acids (FAs). Obesity-prone (OP) rats had a significantly lower proportion than obesity-resistant (OR; P < 0.01) and HFD-paired (P < 0.05) rats. Scatter plot lines with error bars indicate means ± SD; data analyzed with One-way ANOVA and Tukey's post-hoc test. The higher adipose tissue SCD-16 index in OP rats might reflect a genetically inherent difference in FA metabolism pathways that, at least at an early point in DIO, is only found in this specific tissue. In a previous study we found a genetic variant of the SCD-1 gene to be associated with the degree of abdominal obesity and insulin sensitivity in elderly men [40]. In DIO, OP rats may induce SCD-1 to protect them from excess SFA storage. In this case, increased SCD-1 could perhaps be used as a marker for increased fat preference and fat intake. Interestingly, a protective role for SCD-1 could be supported by our finding that the SCD-1 product palmitoleic acid (16:1n-7) was higher in the SAT-TG in OP compared with OR rats, as this FA has been reported to serve as an insulinsensitizing lipokine in the adipose tissue [41]. Alternatively, OR rats may downregulate SCD-1, to prevent adipose tissue storage of FAs and promote oxidation, which may be increased in OR rats [13,14]. If SCD-1 downregulation is protective against DIO, a higher activity could imply that such protective pathways are genetically dysregulated in OP animals. It has been shown that loss of SCD-1 both prevents obesity [19,42,43] and reduces inflammation in adipocytes [44]. A possible target in DIO could therefore be to selectively downregulate SCD-1 in adipose tissue, since a global knock-down of SCD-1 has been shown to lead to deleterious effects [4].
The non-uniform changes in SCD-16 activity across lipid fractions in this study, could be due to differences in species, strain and experimental duration: Li et al. used Wistar rats exposed to a 16-week HFDnotably longer than our 5-week treatment, which may therefore reflect earlier changeswhile Hu et al. fed obesityprone (C57BL/6) and obesity-resistant (FVB) mouse strains a HFD for 8 weeks [13,39]. Both of the aforementioned studies examined SCD-1 mRNA expression, whereas we studied the indices SCD-16 and SCD-18. In our study, the lack of difference in SCD-18, an index also employed to reflect SCD-1 enzyme activity, could be due to the high dietary content of the SCD-18 product oleic acid (OA; 18:1n-9), which could influence the SCD-18 index [31]. At least in our DIO model, the adipose tissue may be one of the first tissues affected by DIO. This was evident as an altered FA metabolism with up-regulated SCD-1 activity index and decreased linoleic acid proportions in OP rats. Changes in the liver, such as those observed in other studies, might then follow as insulin-resistance worsens. This further highlights the importance of the adipose tissue as a key metabolic tissue involved in the development of obesity.
We also found that OP rats had a significantly lower proportion of linoleic acid in SAT-TG, both compared to OR and the calorically restricted rats. Low proportions of LA in humans have previously been found in plasma of obese subjects [27]. In adipose tissue of older men, LA has recently been weakly but positively correlated with insulin sensitivity [45]. Intervention studies suggest that LA-rich diets may decrease SAT and improve insulin sensitivity [46], in line with a possible anti-diabetic effect of LA [47] that should be further investigated.
The higher SCD-1 in OP compared with OR rats could be due to the higher insulin levels seen in the OP rats, as insulin is known to regulate SCD-1 transcription [42,48]. Such interpretation is on the other hand complicated by the elevated leptin concentrations observed in OP rats, as leptin has been shown to decrease both hepatic and adipose tissue SCD-1 expression [49,50]. No difference in adiponectin levels were however found between OP and OR rats. We however did find that plasma adiponectin levels were higher in OP compared with HFD-paired rats, which may be due to an increased adiponectin resistance in OP compared with HFD-paired rats, with an at least initial compensatory higher adiponectin production [51][52][53][54].
OP but not OR rats showed a correlation between food intake and dBw. This strengthens the notion that OR rats have biochemical differences rendering them less efficient at converting food into fat tissue. This difference could perhaps simply be caused by a greater genetic heterogeneity within the OR group. It has however been discovered that on a HFD, OR rats seem to have a proportionately higher fat oxidation [14], and that OR in male but not female rats is associated with better food intake control [9]. Furthermore, a transcriptomic analysis previously found that OP rats have differences in metabolic pathways involving the Krebs cycle, increased ketone body production and a cholesterol transfer promoting fat storage [13]. Our finding that OR rats had a lower food efficiency than OP rats is supported by the findings reported by Chang et al. [14]. The two phenotypes seen in this study therefore seem to present two other major differences: a reduced food intake in OR rats, possibly due to a more optimal reward system [55], and a reduced food efficiency, perhaps linked to lower SCD-1 activity or other biochemical differences in the handling of FAs.
Since we had a control group that was calorically restricted but fed an identical HFD as the one fed ad libitum to the OP and OR groups, we could evaluate the effects of both ad libitum HFD intake and obesity phenotype on FA metabolism. The SCD-16 index of the SAT-TG was significantly higher in OP rats compared with the energy-restricted HFD-paired rats, whereas there was no difference between the OR and HFD-paired rats. Thus, even after moderate weight gain and ad libitum access to the HFD, the OR rats seem to be metabolically similar to the HFD-paired rats. However, for most of the desaturase indices and FAs no differences were observed when ad-libitum-fed rats were compared with the calorically restricted rats, suggesting that on a HFD, many of these parameters are weight-independent.
It should be noted that the diets used to induce DIO more than often differ between different studies [21], thus complicating comparisons between studies. The FA composition of the diet could mask potential differences in FA composition between OP and OR rats that would have been seen if the animals had been given a different diet. However, on the other hand, both the calorierestricted and the ad libitum groups received the same diet, making the results more attributable to metabolic differences. To elucidate what drives the OP and OR phenotypes, initial measurements of the FA composition would have been valuable. Such data could have provided additional predictive markers for susceptibility to DIO; some of which, e.g. plasma-derived FA, could be easily measurable in human subjects. E.g. SCD-1 is known to be correlated with markers of obesity [27], but it is still unknown whether these pathological changes are caused by SCD-1 per se, or if it is a marker of or perhaps even a compensatory mechanism for these ongoing processes [31]. Furthermore, the adipose tissue analyzed herein was derived from the inguinal fat depot, and does not necessarily reflect other adipose tissue depots, as reported by other researchers [35,56].

Conclusion
In this experiment we demonstrate that obesity-prone rats have a significantly higher SCD-16 index, significantly higher proportions of the SCD-1 product 16:1n-7, but a lower proportion of LA in adipose tissue. These changes are largely independent of weight gain and seem to appear before any changes in the plasma or liver tissue, and are paralleled by hormonal changes indicating insulin resistance with possible compensatory changes. Whether a lower SCD-16 index protects against DIO is an interesting possibility that warrants further investigation.

Methods
Thirty-four male outbred Sprague-Dawley rats (Scanbur B&K, Sollentuna, Sweden) were used in this study. At the start of the experiments the rats were 8 weeks old and weighed 352 ± 12 g (mean ± SD). The rats were housed one per cage, in standard macrolon cages (type IV), which had a wood chip bedding and a wooden house as enrichment. All rats had free access to water. Ambient temperature (21-22°C) and humidity (40-50%) were kept constant and a 12-h-light cycle with lights on at 07:00 was used. The rats were allowed one week of adaptation to the animal facility conditions before onset of the experiments. The rats were then randomly divided into two dietary groups. One group of rats were fed a HFD ad libitum (AL-HFD, n = 24); the control group (HFD-paired, n = 10) were fed the HFD but were calorically restricted to follow the weight trajectory of male Sprague-Dawley rats with ad libitum access to regular chow. The HFD (D12451, ResearchDiets) contained 19.79 kJ (4.73 kcal)/g (20% protein, 35% carbohydrates and 45% fat by energy; see Table 3 for composition). During the five weeks of access to the diets, food intake and body weight were measured daily during the first week and thereafter three times a week. After the five-week dietary intervention, the rats were fasted for three hours and were then killed by decapitation, at which point trunk blood was collected. Plasma was isolated by centrifugation and stored at −80°C until analyzed. Liver tissue and subcutaneous adipose tissue (SAT) from the inguinal region were isolated and stored at −80°C. Note that data from a subset of these rats have previously been analyzed and published [56].

Analysis of fatty acid composition
Gas chromatography was used as previously described by Boberg et al. [57] to analyze the FA composition of five different lipid fractions: subcutaneous adipose tissue triacylglycerols (SAT-TG), plasma cholesterol esters (PL-CE), plasma phospholipids (PL-PL), liver phospholipids (liver PL) and liver triacylglycerols (liver TG).
For the SAT and liver tissue, each tissue (5-25 mg of SAT; 40-50 mg of liver tissue) was extracted in 2.5 ml methanol and 5.0 ml chloroform (containing 0.005% butylated hydroxytoluene as an antioxidant). 7.5 ml 0.2 mol/l sodium dihydrogen phosphate (Na 2 H 2 PO 4 ) was added and the extract was left at +4°C over night. Both the liver and SAT chloroform phases were evaporated to dryness under nitrogen. The SAT lipid esters were transmethylated at 60°C overnight after addition of 2 ml of 5% H 2 SO 4 in methanol. For lipid extraction from plasma, 2.5 ml of methanol was added to 0.5 ml of plasma. After thorough mixing of the extract, 5.0 ml Chloroform (containing 0.005% butylated hydroxytoluene as an antioxidant) and 7.5 ml of 0.2 mM sodium dihydrogen phosphate (Na 2 H 2 PO 4 ) were added and were left at + 4°C over night. The chloroform phase was evaporated to dryness under nitrogen and the lipid residue was dissolved in 200-300 μl chloroform.
For the plasma and liver lipid fractions, the lipid esters (triacylglycerols and phospholipids from liver; cholesterol esters and phospholipids from plasma) were separated by thin-layer chromatography (TLC). The chloroform solution was applied to the adsorbent (Silica gel 60G, Merck) containing POPOP as a fluorescent agent. The TLC plates were eluted at room temperature with the solvent system petroleum benzine/diethyl ether/acetic acid (81:18:1, by vol). The lipid fractions were visualized in UV light and the spots were scraped off into vials and the lipid esters were transmethylated at 60°C overnight after addition of 2 ml 5% H 2 SO 4 in methanol.
For all five lipid fractions, the methyl esters were extracted into 3 ml of petroleum ether containing 0.005% butylated hydroxytolvene after addition of 1.5 ml distilled water. The phases were separated after thorough mixing and centrifugation at 1500 × g for 10 min. The petroleum ether phase was pipetted off and the solvent was evaporated under nitrogen. The methyl esters were then redissolved in Uvasol, grade hexane.
The FA methyl esters from the five different fractions were separated by gas chromatography (GC) on a 30-m glass capillary column coated with Thermo TR-FAME (Thermo Electron Corporation, USA), with helium gas as a carrier gas. An Agilent Technologies system consisting of model GC 6890N, autosampler 7683 and Agilent Chem-Station was used. The temperature was programmed to 150-260°C. The FAs were identified by comparing each peak's retention time with FA methyl ester standard Nu Check Prep (Elysian, MN, USA) and are expressed as relative percentages of total FAs. Due to low proportions, not all FAs were quantifiable in all samples.

Glucometabolic and adipokine analyses
Plasma insulin and glucose were analyzed with validated in-house methods, whereas plasma adiponectin and leptin were measured with RIA kits according to the manufacturer's instructions (Millipore AB, Solna, Sweden).
GraphPad Prism v 5.02 (GraphPad Software Inc., San Diego, CA) was used for all statistical calculations. Kolmogorov-Smirnov's test of normality was first used to determine whether data was normally distributed. The following parameters were found to be nonnormally distributed in OP rats: SAT-TG D6D; PL-CE DNL index and EI; liver PL SCD-18, D6D and D5D; and liver TG SCD-16, DNL index, EI, 16:1n-7, and LA. In OR rats the non-normally distributed parameters were: dBw; SAT-TG D6D; PL-CE LA; and PL-PL EI. In HFD-paired rats the non-normally distributed parameters were: initial body weight; plasma insulin concentrations; SAT-TG D6D and ALA; PL-PL ALA and EPA; and liver TG DHA. To detect differences between the three different groups, the investigated normally and non-normally distributed parameters were analyzed using, respectively, One-way ANOVA or Kruskal-Wallis ANOVA. For significant ANOVAs, we proceeded with post-hoc tests to find significant group differences: for One-way ANOVA, Tukey's test was used; for Kruskal-Wallis, Dunn's test was used. When only two groups were compared, the Student's t-test was used for normally distributed data. For correlations, non-normally distributed data was analyzed using Spearman's rank correlation. Values were considered significant at the P < 0.05 level and are presented as means ± SD.

Ethical statement
The Ethical Committee for Animal Experiments in Uppsala approved the experiments and procedures. Animal care procedures followed guidelines of EU (Convention ETS123 and Directive86/609/EEC) and Swedish (Animal Welfare Act SFS1998:56) legislation on animal experiments. experiment, analyzing the data, and writing the manuscript; Å. V. was responsible for analyzing the data, and writing the manuscript; U. R. was responsible for interpreting the data and writing the manuscript; H. B. S. was responsible for planning the study, interpreting the data and writing the manuscript. All authors read and approved the final manuscript.