Altered expression of genes involved in lipid metabolism in obese subjects with unfavourable phenotype

Abstract

Obesity (BMI ≥30 kg/m²) increases the risk of developing lifestyle-related diseases. A subgroup of obese individuals has been described as “metabolically healthy, but obese” (MHO). In contrast to at-risk obese (ARO), the MHO phenotype is defined by a favourable lipid profile and a normal or only slightly affected insulin sensitivity, despite the same amount of body fat. The objective was to characterize the metabolic phenotype of MHO subjects. We screened a variety of genes involved in lipid metabolism and inflammation in peripheral blood mononuclear cells (PBMC). Obese subjects (men and women; 18–70 years) with BMI ≥30 kg/m² were characterized as MHO (n = 9) or as ARO (n = 10). In addition, eleven healthy, normal weight subjects characterized as healthy by the same criteria as described for the MHO subjects were included. We found that with similar weight, total fat mass and fat mass distribution, the ARO subjects have increased plasma levels of gamma-glutamyl transpeptidase and free fatty acids. This group also has altered expression levels of a number of genes linked to lipid metabolism in PBMC with reduced gene expression levels of uncoupling protein 2, hormone-sensitive lipase and peroxisome proliferator-activated receptor δ compared with MHO subjects. The present metabolic differences between subgroups of obese subjects may contribute to explain some of the underlying mechanisms causing the increased risk of disease among ARO subjects compared with MHO subjects.

Introduction

A subgroup of obese individuals has been described as “metabolically healthy, but obese” (MHO) (Karelis et al. 2004a, b; Sims 2001). In contrast to at-risk obese (ARO) subjects, the MHO phenotype is defined by a favourable lipid profile and a normal to high level of insulin sensitivity, despite the same amount of body fat (Karelis et al. 2004a, b; Wildman 2009). Whether the expression of genes involved in lipid metabolism including uptake, oxidation and de novo lipogenesis differs between MHO and ARO subjects remains to be elucidated. To study human hepatic gene expression in vivo is challenging. Peripheral blood mononuclear cells (PBMC) include monocytes and lymphocytes and are easily obtainable in humans. Expression pattern of genes involved in lipid metabolism is reflected in human PBMC (Bouwens et al. 2007, 2008, 2010; Marx et al. 2001, 2002). Additionally, since liver, adipose tissue and PBMC evolutionary derive from the same body compartment (Hotamisligil 2006), PBMC may be a good model system for reflecting important metabolic changes in the liver. In order to understand the underlying molecular mechanisms determining why a subgroup of obese subjects manages to keep a healthy metabolic profile compared with other obese subjects with similar BMI, one approach would be to study the gene expression in PBMC. The aim of the present study was therefore to characterize the metabolic phenotype of MHO subjects compared with ARO subjects by screening a variety of genes involved in lipid metabolism in PBMC.

Methods

Study population

Obese subjects (men and women; 18–70 years) with BMI ≥30 kg/m² were included in this study and were characterized as MHO (n = 9) when three out of the following five criteria were fulfilled (HOMA_ir index ≤1.95; triacylglycerol (TAG) ≤1.7 mmol/L; total cholesterol ≤5.2 mmol/L; LDL cholesterol ≤2.6 mmol/L and HDL cholesterol ≥1.3 mmol/L) or as ARO subjects (n = 10) when four out of the following five criteria were fulfilled (HOMA_ir index ≥1.95; TAG ≥1.7 mmol/L; total cholesterol ≥5.2 mmol/L; LDL cholesterol ≥2.6 mmol/L and HDL ≤1.3 mmol/L). The criteria used in the present study are based on the National Cholesterol Education Program’s Adult Treatment Panel III report (ATP III) for lipid profiles as previously described by Karelis et al. (2004a, b). Furthermore, eleven healthy, normal weight subjects characterized as healthy when four out of five of the same criteria described for the MHO subjects were included. Exclusion criteria were T2DM; kidney, liver, gall bladder, coronary, endocrine or chronic rheumatic disease; malign cancer the last 5 years; hypertension (≥160/100); pregnancy and lactation. Regular use of anti-inflammatory, lipid-lowering and antihypertensive medications was not permitted. The study was approved by the Regional Committee of Medical Ethics (approval no. 6.2008.1368) and by the Norwegian Social Science Data Services (approval no. 19667). Written informed consent for participation was obtained from each participant, and the study complied with the Declaration of Helsinki.

Measurements of body composition

Subjects wore light clothing and no shoes. Two trained persons performed all measurements. Height was measured by a wall-mounted stadiometer to the nearest 0.1 cm. Weight (to the nearest 0.1 kg), BMI and body composition were estimated by the Tanita BC-418 bioelectric impedance analysis (BIA), a hand-to-foot system, according to the manufacturer’s manual (Tanita Corporation, Tokyo, Japan). One kilogram was subtracted from the measured value to correct for light clothing. Waist circumference and hip circumference were measured to the nearest millimetre with a standard, non-stretch tape. Waist circumference was measured at the point midway between the iliac crest and the lower rib margin, while hip circumference was measured at the maximum circumference of the buttocks posterior and the symphysis anterior. All subjects were measured while standing in a relaxed position and with normal respiration. Age, height and body type (all were routinely classified as “standard” body type) were entered before measuring. BIA measurements were carried out at a frequency of 50 kHz.

Blood sampling

The day prior to blood sampling the subjects were told to refrain from alcohol consumption and vigorous physical activity. Venous blood samples were drawn after an overnight fast (≥12 h). Serum was obtained from silica gel tubes [Becton–Dickinson (BD) vacutainer] and kept at room temperature for at least 30 min, until centrifugation (1,500g, 12 min). Serum was kept at room temperature and immediately prepared for subsequent analysis of routine laboratory analyses or aliquoted and stored at −80 °C until further analyses. Plasma was obtained from EDTA tubes (BD vacutainer), immediately placed on ice and centrifuged within 10 min (1500g, 4 °C, 10 min). Plasma samples were aliquoted and stored at −80 °C until further analyses.

Routine laboratory analysis

Fasting serum high-sensitivity C-reactive protein (hsCRP), total cholesterol, LDL cholesterol, HDL cholesterol, TAG, glucose, C-peptide, insulin, HbA1c, gamma-glutamyl transpeptidase (γGT), alkaline phosphatase (ALP), alanine aminotransferase (ALAT) and aspartate aminotransferase (ASAT) were measured by standard methods at a routine laboratory (Fürst Medical Laboratory, Norway). Free fatty acids (FFA) were measured by standard methods at Oslo University Hospital, Norway.

PBMC and RNA isolation and RT-qPCR analysis

After blood collection, PBMC were isolated using BD Vacutainer Cell Preparation tubes according to the manufacturer’s instructions (Becton, Dickinson and Company, NJ 07417, USA). Pellets were frozen and stored at −80 °C prior to RNA isolation.

Total RNA was isolated from all PBMC samples using RNeasy mini kit (Qiagen, Hilden, Germany) with lysis buffer–added β-mercaptoethanol according to the manufacturer’s instructions and stored at −80 °C. RNA quantity and quality measurements were performed using ND 1000 Spectrophotometer (Saveen Werner, Malmö, Sweden) and Agilent Bioanalyser (Agilent Technologies, Santa Clara, CA, USA), respectively. All samples had RNA integrity number (RIN) >8. Four hundred nanogram of RNA was reverse transcribed by high-capacity RNA-to-cDNA kit (cat# 4387406, Applied Biosystems, Foster City, CA, USA). Reverse transcription real-time polymerase chain reaction (RT-qPCR) was performed on an ABI PRISM 7900HT or 7500HT (single assays) Sequence Detector System (Applied Biosystems). Custom-designed TaqMan Low-Density array (LDA) cards (Applied Biosystems) and Taqman or Power SYBR green (designed in Primer Express software version 1) (Applied Biosystems) single assays were used for RT-qPCR amplification of the target genes, all shown in the supplementary table. Due to errors in five of the wells on the LDA card, ATP-binding cassette (ABC) A1, ABCG1, carnitine palmitoyltransferase 1A (CPT1A), CD36 molecule (CD36) and acyl-coenzyme A oxidase 1 (ACOX1) were not included in the analysis. TATA-binding protein (TBP) was selected as housekeeping gene based on running a TaqMan Human Endogenous Control Plate-test (data not shown). Ten microlitre of cDNA (200 ng) was added to water and Taqman master mix to a total volume of 100 μl before applying to the LDA card and centrifuged according to the manufacturer’s protocol. For the single assays, cDNA was diluted 1:5. The relative mRNA level for each transcript was calculated by the ΔΔcycle threshold (Ct) method (Livak and Schmittgen 2001). Briefly, the Ct values of each target gene were subtracted from the Ct values of the housekeeping gene TBP (=ΔCt). Target gene ΔΔCt of MHO or ARO individuals was calculated as ΔCt minus median ΔCt of normal weight individuals. The fold change in mRNA expression was calculated as 2^−ΔΔCt.

Enzyme immunoassay (EIA)

Serum fetuin-A, adiponectin, CD40 ligand (CD40L), leptin, resistin and IL-6 were measured by EIA from R&D Systems (Minneapolis, MN, USA) according the manufacturer’s instructions.

Statistics

As the number of subjects in the study was relatively small, the significance was assessed by nonparametric analyses and correction for multiple testing was performed with Kruskal–Wallis test, followed by Mann–Whitney U test at significant values. All values are presented as median (25–75 percentiles). All analyses were performed using SPSS for Windows (version 19.0). Probability values (asymptotic) were considered statistically significant at a value of ≤0.05.

Results

Characteristics of the participants

The age of the subjects were 49 years (42–63 years; MHO [n = 9]), 52 years (43–59 years; ARO [n = 10]) and 47 years (42–54 years; normal weight [n = 11]) with a BMI of 33 (30–37 kg/m²), 32 (30–34 kg/m²) and 23 (21–24 kg/m²), respectively. There were no significant differences between the two obese groups in weight or BMI (Table 1). As expected, the parameters included in the inclusion criteria were significantly different between the ARO subjects compared with the MHO subjects (Table 1). hsCRP was significantly higher in both MHO subjects (P = 0.033) and ARO subjects (P = 0.020) compared with healthy, normal weight subjects. There was no significant difference in the level of hsCRP between the two obese groups.

Table 1 Characteristics of the population (median with 25–75 percentiles)

Anthropometric measurements and adipokines in serum

Total fat mass, trunk fat mass, fat-free mass, hip and waist circumference were measured, and waist/hip ratio was calculated in all subjects (Fig. 1). Anthropometric measures concerning fat distribution did not differ between the two obese groups. Both MHO and ARO subjects had larger total fat mass, trunk fat mass, hip circumference, waist circumference, waist/hip ratio (P < 0.001 for all parameters for both groups) and fat-free mass (P = 0.053 and P = 0.003) than normal weight subjects. In addition, the circulating levels of leptin, adiponectin and resistin were measured. There was no difference in any of these parameters between the two obese groups. As expected, normal weight subjects had lower levels of leptin compared with both obese groups, while adiponectin levels were higher only compared with ARO subjects (Fig. 2).

Fig. 1

Anthropometric measurements in normal weight (n = 11), ARO (n = 10) and MHO subjects (n = 9). Total fat mass (a), trunk fat mass (b), fat-free mass (c), hip circumference (d), waist circumference (e) and waist/hip ratio (f). ARO at-risk obese, MHO metabolically healthy but obese. Data are shown as median, and bars indicate (25–75) percentiles

Fig. 2

Plasma levels of leptin (a), adiponectin (b) and resistin (c) in normal weight (n = 11), ARO (n = 10) and MHO (n = 9) subjects. ARO at-risk obese, MHO metabolically healthy but obese. Data are shown as median, and bars indicate (25–75) percentiles

Circulating markers of hepatic fat accumulation and inflammation

We measured established markers of liver dysfunction and liver fat accumulation such as γGT, ASAT, ALAT, ALP as well as fetuin-A. γGT levels were significantly higher in ARO subjects compared with MHO (P = 0.008) and normal weight subjects (P < 0.001). In addition, MHO subjects had increased γGT concentration levels compared with normal weight subjects (P = 0.043) (Fig. 3). There were no significant differences in ALAT, ASAT, ALP or fetuin-A between any of the groups (Fig. 3). We also analysed the circulating level of the inflammatory markers IL6 and CD40L. There was no significant difference between the two obese groups; however, both obese groups had significantly elevated serum levels of IL6 and CD40L compared with normal weight subjects (data not shown).

Fig. 3

Plasma levels of γGT (a), ALP (b), ASAT (c), ALAT (d) and fetuin-A (e) in normal weight (n = 11), ARO (n = 10) and MHO subjects (n = 9). ARO at-risk obese, MHO metabolically healthy but obese. Data are shown as median, and bars indicate (25–75) percentiles

Free fatty acids in plasma

The availability of FFA from circulation may influence the progression of fat accumulation in the liver (Mendez-Sanchez et al. 2007). ARO subjects had significant higher plasma levels of FFA (P = 0.020) compared with MHO subjects (Fig. 4).

Fig. 4

Plasma levels of FFA in ARO (n = 10) and MHO subjects (n = 9). ARO at-risk obese, FFA free fatty acids, MHO metabolically healthy but obese. Data are shown as median, and bars indicate (25–75) percentiles

PBMC gene expression

The increase in plasma FFA and γGT in ARO subjects may point to development of fat accumulation and lipid disturbances in the liver. Expression and regulation of genes involved in lipid metabolism in PBMC have previously been shown to reflect hepatic changes (Bouwens et al. 2007, 2008). The expression levels of selected genes involved in lipid uptake, transport, lipolysis, de novo lipogenesis and fatty acid oxidation were analysed (Table 2). Additionally, we analysed the expression of the transcription factors PPARδ and LXRα which have been described as fatty acid sensors regulating genes involved in lipid metabolism (Sanderson et al. 2009; Strable and Ntambi 2010). Whereas there were no differences in PBMC gene expression levels of UCP2, HSL and PPARδ between normal weight and MHO subjects, ARO subjects had significantly lower expression levels of these genes compared with both normal weight subjects (P = 0.041, P = 0.035 and P = 0.05, respectively) and MHO subjects (P = 0.003, P = 0.003 and P = 0.006, respectively) (Fig. 5). Minor differences between normal weight subjects and the two obese groups were observed for some genes; however, for the majority of genes, we found no differences in the expression level between the groups (Table 2).

Table 2 PBMC mRNA expression levels in normal weight, at-risk obese (ARO) and metabolically healthy but obese (MHO) subjects (median with 25–75 percentiles)

Fig. 5

PBMC mRNA levels of UCP2 (a), HSL (b) and PPARδ (c) in normal weight (n = 11), ARO (n = 10) and MHO subjects (n = 9). The figures are shown as fold change (median) from normal weight subjects, and the bars indicate (25–75) percentiles. ARO at-risk obese, HSL hormone-sensitive lipase, MHO metabolically healthy but obese, PBMC peripheral blood mononuclear cells, PPARδ peroxisome proliferator-activated receptor delta, UCP2 uncoupling protein 2

Discussion

In the present study, we have screened a variety of markers associated with fat accumulation in the liver and lipid metabolism in order to characterize differences in metabolic pathways between MHO and ARO subjects. We show that despite similar weight, total fat mass and fat mass distribution, ARO subjects have increased plasma levels of γGT and FFA, and reduced PBMC gene expression level of UCP2, HSL and PPARδ compared with MHO subjects. These metabolic differences may contribute to explain some of the underlying mechanisms causing increased risk of disease among ARO subjects compared with MHO subjects.

MHO subjects are characterized with a favourable lipid profile and normal insulin sensitivity despite an excessive fat mass. The mechanistic explanation as to how a subgroup of obese subjects manages to obtain these metabolically favourable conditions is poorly understood. Elevated levels of circulating γGT and ALAT, without fatty liver or hepatic dysfunction, in healthy and asymptomatic subjects may predict future development of metabolic disease (Vozarova et al. 2002; Wannamethee et al. 2005). We found that ARO subjects had higher levels of the liver marker γGT than both MHO and normal weight individuals. This is in line with reports from a previous study where the level of ALAT, even within the normal range, was shown to be significantly associated with obesity phenotypes (Mojiminiyi et al. 2010). We were unable to distinguish between subcutaneous and visceral fat; however, there were no differences between MHO and ARO subjects regarding trunk fat mass or waist circumference. Even so, the level of the different fat depots may still differ between these two obese groups as shown by Brochu and co-workers who found that the MHO subjects had up to 49 % less visceral fat mass (Brochu et al. 2001).

In the present study, we found genes related to β-oxidation to be changed. PBMC gene expression level of UCP2 in ARO subjects was significantly decreased compared with both normal weight and MHO subjects. UCP2 has previously been associated with obesity (Dalgaard 2011; Dalgaard and Pedersen 2001). UCP2 is important for uncoupling respiration and promotes mitochondrial fatty acid oxidation (Fisler and Warden 2006), although the physiological significance of its uncoupling activity is discussed (Pecqueur et al. 2009).

We found that ARO subjects had significant lower expression of PPARδ than the MHO subjects. PPARδ has been suggested to be a fatty acid sensor, and recent studies have shown that activation of PPARδ may lead to beneficial effects on lipoprotein and glucose metabolism (Ooi et al. 2011). Activation with the synthetic PPARδ agonist GW501516 decreases plasma TAG and plasma FFA in dyslipidaemic subjects with central obesity (Ooi et al. 2011). PPARδ appears to sense and respond to elevated FFA by up-regulating genes such as Lipin 2 (Sanderson et al. 2009), and in mice models PPARδ agonists have previously been shown to protect against weight gain, steatosis, reduced HDL levels and elevated TAG levels, while PPARδ null mice were more susceptible to weight gain (Barish et al. 2006). This collectively suggests that activation of PPARδ protects against obesity and/or dyslipidaemia. One interpretation of our results may be that ARO subjects have an impaired lipid metabolism with a slower fatty acid turnover due to impaired fatty acid oxidation. In the present study, ARO subjects had higher plasma levels of FFA than MHO subjects. If fatty acid oxidation and lipid storage is impaired, one might expect that FFA will remain in plasma, leading to increased levels.

A tight regulation of hepatic TAG synthesis, hydrolysis and β-oxidation is required to prevent lipid accumulation. HSL was initially observed in adipose tissue, but observations in mice suggest that HSL is also expressed in liver (Mulder et al. 2003; Reid et al. 2008; Sekiya et al. 2008). Hepatic overexpression of HSL (in cells) is shown to reduce hepatic steatosis by promoting β-oxidation, and it has therefore been suggested that HSL may be a therapeutic target for the treatment of fatty liver in human subjects (Reid et al. 2008).

Our data are in line with the previous observation from Taskinen and co-workers (Taskinen et al. 2011), who showed that hypertriglyceridaemia in obese men with similar BMI and waist circumference was caused by increased hepatic secretion of VLDL (induced by increased liver fat) in combination with severely impaired clearance of triglyceride-rich VLDL₁ particles. The liver fat content seemed to be the driving force for the overproduction of VLDL₁-triglycerides and apoB. In the present study, we extend and support these novel findings by showing reduced PBMC gene expression, potentially reflecting liver expression, of UCP-2 and PPARδ in ARO subjects compared with MHO subjects, further supporting such a notion.

Furthermore, also in line with our data, in a recent study, the specific PPARδ agonist GW501516 showed potent hypolipidaemic actions in a two-week randomized controlled study with six healthy moderately overweight subjects in each intervention arm (Riserus et al. 2008). The fasting plasma TAG concentration was reduced by 30 %, and there was a concomitantly lowering of the fasting non-esterified fatty acid concentrations which were followed by an increase in the oxidation of fatty acids measured as increased gene expression of CPT1b in skeletal muscle. This mechanism was shown to be unique for the PPARδ agonist and not seen with the PPARα agonist. In the present study, we show that ARO subjects have reduced PBMC gene expression of PPARδ, increased level of non-esterified fatty acids in circulation and decreased PBMC expression UCP-2, which promotes mitochondrial fatty acid oxidation.

The strength of the study, even though being descriptive in nature, is that we have thorough clinical description, extensive anthropometric measurements, gene expression analysis and plasma markers in the two obese groups. The limitations of the study are the relative small number of subjects in each group, the lack of measurements of FFA in the normal weight subjects and the skewed gender distribution. When recalculating the data, using only the male population in each group, the same pattern was found between MHO and ARO subjects and the level of significance was maintained for most of the markers (γGT, FFA, UCP2, PPARδ and HSL).

In conclusion, ARO subjects had increased plasma levels of γGT and FFA, and reduced PBMC gene expression level of UCP2, HSL and PPARδ compared with MHO subjects. These metabolic differences may contribute to explain some of the underlying mechanisms responsible for the increased risk of disease among ARO subjects compared with MHO subjects. However, our results need to be confirmed in larger studies which also include clinical examination of liver fat content. A better understanding of the metabolic phenotypes of obese subjects may have important implications for future therapeutics.