Introduction

Serum lipid and lipoprotein abnormalities are inherent components of the insulin resistance syndrome [1, 2]. Insulin resistance in the liver is characterised by the inability of insulin to inhibit hepatic very-low-density lipoprotein (VLDL) production [3]. This leads to elevation of serum triacylglycerols (TGs) [4, 5] and a low HDL-cholesterol concentration [6].

The relationship between fasting serum TGs and insulin resistance measured using the hyperinsulinaemic–euglycaemic clamp technique in 1,308 non-diabetic individuals was weak [7]. However, TGs consist of multiple different fatty acids and are found not only in VLDL but also in other lipoprotein particles [8]. A recent cross-sectional study indicated that the fatty acid composition of TG determined using GC is differently related to surrogate markers of insulin resistance [9]. The proportion of saturated TG 16:0 fatty acid was positively associated with fasting serum insulin concentrations and that of essential 18:3 n-6 fatty acid negatively associated [9]. This raises the possibility that the specific increases in serum TGs containing saturated fatty acids and decreases in the TGs containing polyunsaturated fatty acids may characterise insulin resistance and thus weaken the association between insulin resistance and total serum TG concentrations. The relationships between individual serum TGs and insulin resistance have not been previously studied.

The type of fat in the diet is an important determinant of the fatty acid composition in serum [10, 11]. Dietary intervention studies have shown that an increased intake of saturated fat is associated with an increase in insulin resistance [12], whereas an increased intake of polyunsaturated fat has been shown to improve insulin sensitivity [13]. The majority of serum fatty acids are carried by lipoproteins [8], some of which are considered as pro-atherogenic while others as anti-atherogenic. Long-term prospective studies show decreased essential and increased saturated fatty acid concentrations in total serum to predict the metabolic syndrome [14] and cardiovascular disease [15]. However, it is unknown how fatty acid composition of individual lipoproteins relates to insulin sensitivity.

Lipoprotein particles have traditionally been analysed by measuring their total protein, phospholipid, cholesteryl ester and TG content [8]. Modern lipidomics techniques such as ultra-performance liquid chromatography coupled to MS (UPLC/MS) afford analysis of multiple abundant lipids (e.g. TG, cholesteryl esters, phosphatidylcholines) as well as bioactive lipid species (e.g. ceramides, plasmalogens and lysophosphatidylcholines) in parallel [16]. Although the bioactive lipid species represent a minor component of the lipid milieu, several lines of evidence suggest that they may play a role in, for example, inflammation [17] and antioxidative defence [18], which might contribute to atherogenity or cardioprotective properties of different lipoprotein particles. The relative abundance of reactive lipid species in different lipoproteins has not been previously examined in any species.

In the present study, we performed a comprehensive lipidomics analysis to determine concentrations of individual TGs and fatty acid composition in serum and major lipoprotein particles. We then analysed how changes in different TGs and fatty acid composition are related to features of insulin resistance and abdominal obesity. We also determined whether and in which lipoprotein particles bioactive lipids are present.

Methods

Participants and study design

A total of 16 non-diabetic individuals were recruited by newspaper advertisements based on the following inclusion criteria: (1) age 18–60 years; (2) no known acute or chronic disease based on history, physical examination and standard laboratory tests (blood counts, serum creatinine, thyroid-stimulating hormone and electrolyte concentrations) and ECG; and (3) alcohol consumption <20 g/day. The study protocol was approved by the ethics committee of the Helsinki University Central Hospital, and each participant provided written informed consent. Blood samples were taken after an overnight fast for measurement of fasting plasma glucose, fasting serum insulin, HbA1c, serum lipids and lipoproteins as previously described [19, 20]. Waist circumference was measured midway between spina iliaca superior and the lower rib margin, and hip circumference at the level of the greater trochanters [19]. BP was measured in the sitting position after 10–15 min of rest using a random-zero sphygmomanometer (Erka, Germany). None of the participants used oral contraception or hormonal replacement. One woman was postmenopausal. Clinical characteristics of study participants are listed in Table 1.

Table 1 Clinical characteristics of the study participants
Full size table

Lipoprotein fractionation

The lipoprotein fractions were separated by sequential flotation in an ultracentrifuge [21] (Beckman L-60; Beckman, Palo Alto, CA, USA) at 4°C by using a Ti 50.4 Beckman rotor. Fresh serum was overlayered with 0.16 mol/l NaCl, 1 mmol/l EDTA (density = 1.006 g/ml). The VLDL, intermediate-density lipoprotein (IDL), LDL and HDL fractions 2 and 3 (HDL2 and HDL3) were fractionated using a KBr/NaCl gradient as previously described [20]. The normalisation based on total protein content was used to estimate the lipid content of an individual lipoprotein particle. Protein concentrations in lipoprotein fractions were measured by a modified Lowry assay method [22].

Lipidomic analyses

UPLC/MS-based analysis was used for lipid profiling. The applied platform affords broad screening of multiple lipid classes, including TGs, cholesteryl esters, and major phospholipids, from total lipid extracts within a single sample run and has been previously described in detail [16]. The platform is not optimal for detection of some negatively charged phospholipids such as phosphatidylserines and phosphatidylinositols. The method does not cover glycosphingolipids or the low molecular mass reactive lipids such as eicosanoids and NEFAs.

In short, 10 µl of an internal standard mixture and 10 µl NaCl (0.9%) were added to 10 µl of the sample. Lipids were extracted from the samples with 100 µl chloroform:methanol (2:1, vol./vol.) solvent and the sample was homogenised with a glass rod. Lipid extracts were analysed on a Q-ToF Premier mass spectrometer (Waters, Milford, MA, USA) combined with an Acquity UPLC (Waters). The lipid profiling was carried out using positive ion mode. The data were collected at mass range of m/z 300–1,200 with scan duration of 0.2 s.

The obtained data were converted into netCDF file format using Dbridge software from MassLynx (Waters). The converted data were processed using MZmine software version 0.60 [23]. All reported lipids were identified with tandem mass spectrometry (UPLC/MS/MS) in positive and if necessary negative ion mode. The lipid concentrations were calibrated with class-specific internal standards.

Fatty acid analyses

Analyses of serum and lipoprotein fatty acids were performed by adding 20 μl of internal standard (heptadecanetrienoate TG C17:0 + NEFA C17:0, 500 mg/l), 200 µl choloroform:methanol (2:1, vol./vol.) solvent and 20 µl 0.15 mol/l NaCl to 50 µl of the sample. The samples were prepared for GC analysis as previously described [24]. GC was carried out using an HP 5890 GC-FID System that was directly connected to a FID detector (Agilent Technologies, Waldbronn, Germany).

Other measurements

Serum unbound insulin concentrations were measured with a time-resolved immunofluorometric assay (PerkinElmer, Waltham, MA, USA). Inter-assay CV was <4% in the concentration range 40–720 pmol/l. Cross-reaction with proinsulin and C-peptide was ≤0.1% [25]. HOMA for insulin resistance (HOMA-IR) was calculated using the formula of Matthews et al. [26].

Statistical methods

Data are shown as means ± SD. Correlation analyses were performed using Spearman’s non-parametric rank correlation coefficient. The weighted correlation coefficient for each serum TG was calculated by multiplying the abundance (%) of each serum TG by the Spearman’s correlation coefficient. Five serum TGs with the greatest positive (representing 29.4% of all TGs) and negative (36.9%) weighted correlation coefficients are shown in Fig. 1. The results remained essentially unchanged when men were excluded from the analyses, and when the concentrations of lipids (μmol/ml) were used instead of the protein-normalised values. Calculations were made using GraphPad Prism version 4.00 for Windows (GraphPad Software, San Diego, CA, USA), SysStat Statistical Package (SysStat version 10; SysStat, Evanston, IL, USA) and SPSS 14.0 for Windows (SPSS, Chicago, IL, USA). p < 0.05 was considered statistically significant.

Fig. 1
figure 1

Relationships between TG concentrations in different lipoprotein fractions and HOMA-IR. Spearman rank correlations were used to correlate the concentrations of each individual TGs (μmol/ml) within each lipoprotein fraction with HOMA-IR. The colour code corresponds to Spearman’s correlation coefficient (r Spearman) values as shown. The lipids which were not detected within lipoproteins are marked in grey. Cross-hatching: p < 0.05. The five serum TGs with the greatest positive (representing 29.4% of all serum TGs) and negative (36.9%) weighted correlation coefficients are marked in red and blue, respectively

Full size image

Results

Serum and lipoprotein TGs

With the lipidomic approach, we identified 45 different TGs in serum. The most abundant serum TG molecules were TGs containing one palmitic (16:0) and two oleic (18:1) acids (TG[18:1/16:0/18:1], 19.9% of total) and TG(18:2/16:0/18:1) (18.3% of total). The Spearman’s correlation coefficient for each individual serum TG and HOMA-IR were calculated.

Serum TGs with the greatest positive weighted correlation coefficient were significantly positively related to HOMA-IR when measured in VLDL, IDL and LDL, but not in HDL2 or HDL3 (Fig. 1). For the serum TGs that were negatively correlated with HOMA-IR, no clear trends were observed for the respective associations in different lipoprotein (Fig. 1). The results remain essentially unchanged if HOMA-IR was replaced by waist circumference (data not shown).

Serum and lipoprotein fatty acid composition

Total serum esterified palmitic (16:0), palmitoleic (16:1 n-7), and oleic (18:1 n-9) acids were positively and linoleic acid (18:2 n-6) negatively related to HOMA-IR. Analyses on total esterified fatty acids within lipoproteins revealed that palmitic acid was positively related to HOMA-IR when measured in VLDL, IDL and LDL, but not in HDL2 or HDL3. Linoleic acid was negatively related to HOMA-IR in all lipoproteins (Table 2). The distribution of esterified fatty acids across lipoproteins is shown in Electronic supplementary material (ESM) Fig. 1 and given in ESM Table 1.

Table 2 Relationships between esterified fatty acids and HOMA-IR of all the lipid species, expressed as molar % in serum and its different lipoprotein fractions
Full size table

Bioactive lipid species

An overview of the lipid composition of lipoprotein fractions is shown in Fig. 2. Ceramides were found only in VLDL and LDL, whereas ethanolamine plasmalogens were found only in LDL and HDL2. Lysophosphatidylcholines and ether-linked phosphatidylcholines were present in all lipoproteins with the greatest abundance in HDL2, HDL3 and LDL. Concentrations of bioactive lipids in different lipoproteins are shown in Fig. 3 and given in ESM Table 2.

Fig. 2
figure 2

Lipidomic characterisation of different lipoprotein fractions, based on per cent of the total lipid content within each fraction. Major (a) and minor (b) lipid classes are shown. Cer, ceramide; ChoE, cholesteryl ester; PC, phosphatidylcholine; PE, phosphatidylethanolamine; SM, sphingomyelin

Full size image
Fig. 3
figure 3

Absolute (μmol/ml) concentrations of ceramides (Cer), lysophosphatidylcholines (lysoPC), ethanolamine plasmalogens (PE(p)) and ether-linked phosphatidylcholines (PC(e)) in lipoproteins

Full size image

Discussion

Analyses on individual TGs in the present study revealed that serum TG molecules containing saturated and monounsaturated fatty acids, such as TG(16:0/16:0/18:1) and TG(16:0/18:1/18:0), correlated positively, whereas those containing essential fatty acids, such as TG(18:1/18:2/18:2), negatively with HOMA-IR and waist circumference. These TGs correlated positively with features of insulin resistance only within VLDL, IDL and LDL, but not within HDL. Fatty acid analyses revealed that saturated 16:0 fatty acid (palmitate) was positively related to HOMA-IR only when measured in VLDL, IDL and LDL, whereas unsaturated 18:2 n-6 fatty acid (linoleate) was inversely related to HOMA-IR in all lipoproteins. In addition, we found that lipotoxic ceramides are present within VLDL and LDL in human serum.

We found that saturated, monounsaturated and polyunsaturated fatty acids containing TGs within the VLDL–IDL–LDL axis and HDLs were differently related to insulin resistance. This might explain why the association between total serum TGs and insulin resistance is weak [7].

Saturated fatty acids containing TGs correlated positively with features of insulin resistance in VLDL, which is produced by the liver, and in IDL and LDL, which are derived from VLDL. The fatty acids in VLDL TGs have been shown to closely reflect sources of intrahepatocellular TGs [27]. Hepatic insulin resistance is associated with increased liver fat content even independently of obesity [28]. The fatty liver overproduces large TG-rich VLDL particles [29, 30]. In individuals with a fatty liver, the contribution of TGs derived from de novo lipogenesis to both hepatic and VLDL TGs is significantly increased [27, 31]. Accelerated hepatic de novo lipogenesis leads to increased incorporation of palmitate (16:0), stearate (18:0), and oleate (18:1 n-9) into VLDL particles [32]. This in turn leads to decreased relative amounts of linoleate (18:2 n-6) in VLDL TGs [32]. In the present study, the significant positive relationships between features of insulin resistance and VLDL TGs containing saturated and monounsaturated fatty acids, such as TG(16:0/16:0/18:1) and TG(16:0/18:1/18:0), could therefore reflect increased hepatic de novo lipogenesis.

The analyses on esterified fatty acids of all lipids revealed that linoleic fatty acid decreases with increasing HOMA-IR and waist circumference in all lipoproteins. On the other hand, changes in saturated and monounsaturated fatty acids with respect to insulin resistance were lipoprotein-specific. Epidemiological prospective studies have shown that high proportions of esterified palmitic and stearic, and low proportions of esterified linoleic fatty acids in serum are risk factors for type 2 diabetes [33–35], myocardial infarction [36, 37], the metabolic syndrome [38] and overall and cardiovascular mortality [39]. The relationship between the amount of polyunsaturated fat in the diet and the respective proportions of serum esterified polyunsaturated fatty acids is strong [40, 10]. The same relationships for saturated fatty acids are weaker [10, 41], and those between monounsaturated fatty acids non-existent [10, 41]. These data would suggest that endogenous sources of saturated and monounsaturated fatty acids significantly contribute to serum and lipoprotein fatty acid concentrations, and that exogenous sources of fatty acids influence fatty acid composition in a similar manner in all lipoproteins.

Saturated esterified fatty acids were positively correlated with features of insulin resistance only in VLDL, IDL and LDL, and monounsaturated fatty acids only in HDL. The proportion of saturated fatty acids in serum lipid esters have been shown to be positively related to insulin sensitivity in some [42–44] but not all [45] studies. A recent study of 294 healthy 63-year-old Swedish men showed that correlation coefficients between HOMA-IR and plasma saturated fatty acids in phospholipids is 0.29 [44]. Both weak, but significant [43, 45], and insignificant [44] relationships between serum palmitoleic fatty acid and HOMA-IR have been reported. The specific pattern of the changes in lipoprotein fatty acid composition may explain the divergent results of the associations between features of insulin resistance and total serum fatty acid composition. Our results suggest that analyses of lipoprotein-specific rather than total serum fatty acid composition may more accurately reflect changes of fatty acids characteristic to insulin resistance.

Of the minor lipid sub-classes, proinflammatory and lipotoxic ceramides [17] were found within LDL and VLDL. These data in humans are novel and consistent with the ability of hepatocytes to synthesise ceramides de novo and incorporate into them into VLDL [46]. The absence of ceramides in IDL in the present study may be due to undetectably low concentrations of ceramides in IDL. Ethanolamine plasmalogens, known as endogenous antioxidants [18], were found only within LDL and HDL2 fractions. This may contribute to atherogenecity of low HDL2 [47].

Our results indicate that serum concentrations of specific TGs, such as TG(16:0/16:0/18:1) or TG(16:0/18:1/18:0), may be more precise markers of insulin resistance than total serum TG concentrations. Analyses of fatty acid composition in lipoproteins rather than total serum may help to define the role of saturated, monounsaturated and polyunsaturated fatty acids in the pathogenesis of the insulin resistance syndrome.