Adipose tissue gene expression analysis reveals changes in inflammatory, mitochondrial respiratory and lipid metabolic pathways in obese insulin-resistant subjects
BMC Medical Genomics volume 5, Article number: 9 (2012)
To get insight into molecular mechanisms underlying insulin resistance, we compared acute in vivo effects of insulin on adipose tissue transcriptional profiles between obese insulin-resistant and lean insulin-sensitive women.
Subcutaneous adipose tissue biopsies were obtained before and after 3 and 6 hours of intravenously maintained euglycemic hyperinsulinemia from 9 insulin-resistant and 11 insulin-sensitive females. Gene expression was measured using Affymetrix HG U133 Plus 2 microarrays and qRT-PCR. Microarray data and pathway analyses were performed with Chipster v1.4.2 and by using in-house developed nonparametric pathway analysis software.
The most prominent difference in gene expression of the insulin-resistant group during hyperinsulinemia was reduced transcription of nuclear genes involved in mitochondrial respiration (mitochondrial respiratory chain, GO:0001934). Inflammatory pathways with complement components (inflammatory response, GO:0006954) and cytokines (chemotaxis, GO:0042330) were strongly up-regulated in insulin-resistant as compared to insulin-sensitive subjects both before and during hyperinsulinemia. Furthermore, differences were observed in genes contributing to fatty acid, cholesterol and triglyceride metabolism (FATP2, ELOVL6, PNPLA3, SREBF1) and in genes involved in regulating lipolysis (ANGPTL4) between the insulin-resistant and -sensitive subjects especially during hyperinsulinemia.
The major finding of this study was lower expression of mitochondrial respiratory pathway and defective induction of lipid metabolism pathways by insulin in insulin-resistant subjects. Moreover, the study reveals several novel genes whose aberrant regulation is associated with the obese insulin-resistant phenotype.
Instead of merely being a storage depot, adipose tissue is now recognized as an important regulator of energy homeostasis and a central player in the development of insulin resistance [1, 2]. The development of insulin resistance is closely associated with obesity . However, it is not known why all obese subjects are not insulin-resistant and insulin resistance can also be observed in lean subjects, highlighting the fact that the molecular mechanisms underlying insulin resistance are still largely unknown.
A number of transcriptome studies have reported differences in adipose tissue mRNA expression at fasting state between obese and lean subjects (reviewed in ). Most of these investigations have been performed after an overnight fast [4–14] and they compare differences between obese and non-obese subjects [5, 6, 10–14], predominantly in males or mixed groups consisting of both males and females.
A number of studies employing different set-ups, including global transcriptome analyses, have implicated low-grade adipose tissue inflammation as one of the central factors involved in the development of insulin resistance . As an example, analysis of subcutaneous adipose tissue transcript profiles in monozygotic twins revealed elevated expression of several pro-inflammatory cytokines and complement components in the in the obese, insulin-resistant co-twins as compared to the non-obese ones .
Using human whole genome microarrays, Dahlman et al.  observed down-regulation of mitochondrial electron transport chain genes in both visceral and subcutaneous adipose tissue of type 2 diabetic subjects during fasting. A study on monozygotic twins discordant for BMI reported decreased branched-chain amino acid catabolism and reduced mitochondrial DNA copy number in the adipose tissue of obese co-twins, further emphasizing a role for mitochondrial energy and amino acid metabolism in obesity and insulin resistance .
Baranova et al.  reported mRNA expression patterns in visceral adipose tissue of morbidly obese patients compared with non-obese controls. The differentially expressed genes included those related to lipid and glucose metabolism, membrane transport, and promotion of the cell cycle. Using isolated adipocytes, Yang et al.  discovered through comparison of gene expression between non-diabetic insulin-resistant first degree relatives of T2D patients and healthy controls, that insulin resistance is associated with impaired adipogenesis.
While previous studies have identified a number of factors involved in adipogenesis, adipose tissue lipid and energy metabolism, as well as in inflammatory processes to be abnormally expressed, the alterations could be unrelated to insulin action and merely reflect the impact of factors other than insulin on gene expression in obese as compared to non-obese subjects. To relate the abnormalities to impaired insulin response, comparison of differences in acute regulation of the transcriptome by insulin between obese and non-obese subjects would seem of interest. To our knowledge there are no such data available. Therefore the present study was undertaken to compare transcriptional differences in subcutaneous adipose tissue derived from obese insulin-resistant and lean insulin-sensitive women before and during euglycemic hyperinsulinemia.
The nature and potential risks of the study were explained to all subjects prior to obtaining their written informed consent. The study was carried out in accordance with the principles of the declaration of Helsinki. The protocol was approved by the ethics committee of the Helsinki University Central Hospital.
A total of 20 non-diabetic Caucasian women were recruited based on the following inclusion criteria: 1) age 18-60 years; 2) no known acute or chronic disease other than obesity based on history, physical examination and standard laboratory tests (blood counts, serum creatinine, thyroid stimulating hormone, electrolyte concentrations, and electrocardiogram); and 3) body mass index (BMI) < 40 kg/m2. Other exclusion criteria included pregnancy or treatment with drugs that may alter glucose tolerance. Data on expression of some individual genes from these women have previously been published [17, 18]. Whole body insulin sensitivity was measured in each subject using the euglycemic hyperinsulinemic clamp technique, as described . The duration of the insulin infusion was 6 hours and insulin infusion rate 1 mU/kg·min.
Metabolic studies were performed after an overnight fast. Two 18-gauge catheters (Venflon; Viggo-Spectramed, Helsingborg, Sweden) were inserted, one in an antecubital vein for infusion of insulin (Insulin Actrapid; Novo Nordisk, Denmark) and glucose, and another retrogradely in a heated hand vein to obtain arterialized venous blood for measurement of glucose concentrations every 5 min and serum free insulin concentration every 30 min. The rate of the continuous insulin infusion was 1 mU kg - 1 min - 1 for 6 h. Normoglycemia was maintained by adjusting the rate of a 20% glucose infusion based on plasma glucose measurements from arterialized venous blood every 5 min. Wholebody insulin sensitivity was determined from the glucose infusion rate required to maintain normoglycemia between 30 and 360 min .
Blood samples were taken after an overnight fast for measurement of plasma glucose, serum insulin, C-peptide, serum triglyceride and total and HDL cholesterol concentrations, as described .
Adipose tissue biopsies, total RNA preparation and microarray processing
Needle aspiration biopsies of abdominal subcutaneous fat were taken under local intracutaneous anesthesia at baseline (0 h) and after 3 and 6 hours of hyperinsulinemia from the left and right lower abdominal region  and frozen in liquid nitrogen until analysis. Total RNA was prepared from frozen fat tissue (on the average 250 mg) using the RNeasy Lipid Tissue Mini Kit (Qiagen) according to the manufacturer's protocol. Quality of RNA was analyzed using the 2100 Bioanalyzer platform (Agilent Technologies). For microarray gene expression analysis, the five most insulin-sensitive and the five most insulin-resistant patients were selected based on their whole body insulin sensitivity (M-value). Two micrograms of total RNA from time point 0 h (fasting) and 3 h (hyperinsulinemia) were processed according to the conventional Affymetrix eukaryotic RNA labeling protocols (Affymetrix, Santa Clara, CA). Hybridization, staining and washing of the Affymetrix U133 Plus 2.0 chips were performed using the Affymetrix Fluidics Station 450 and Hybridization Oven 640 under standard conditions.
Analysis of the microarray expression data was done using the R/Bioconductor through a graphical user interface, Chipster (v1.4.3, CSC, Finland, http://chipster.csc.fi/). The microarray data has been submitted to GEO database (accession number, GSE26637). Pre-processing of the expression data was performed using GC-RMA algorithm and re-annotated probe set libraries by the University of Michigan http://brainarray.mbni.med.umich.edu. Expression data were filtered before pathway analysis by flags. All probes with status 'present' or 'marginal' in at least 3 out of 5 chips in one of the four condition groups (fasted insulin-sensitive, fasted insulin-resistant, hyperinsulinemic insulin-sensitive, hyperinsulinemic insulin-resistant) were accepted for analysis. Probes passing these criteria were then divided into up and down-regulated probes by their fold change in expression. Analysis was done separately for fasting samples (fasted insulin-resistant/fasted insulin sensitive fold change) and for hyperinsulinemia samples (hyperinsulinemic insulin-resistant/hyperinsulinemic insulin-sensitive fold change). Probe lists were ranked according to statistical significance (Student's t-test) and submitted to the pathway analysis software. The top 10 up and down-regulated pathways are reported. Similarly, we analyzed which biological pathways show the strongest response in gene expression after systemic stimulation by insulin. Analyses were performed separately for the insulin-resistant and the insulin-sensitive groups (hyperinsulinemic insulin-sensitive/fasted insulin-sensitive fold change and hyperinsulinemic insulin-resistant/fasted insulin-resistant fold change). After defining the most insulin responsive pathways, we selected two of them, and compared their expression between the study groups during hyperinsulinemia. Pathway analysis was performed using an in-house developed non-parametric pathway analysis software to resolve which gene sets from GO classification are enriched in a list of genes as described recently . As a feature of the dendrogram-like structure of the GO classifications, gene sets get progressively larger and less descriptive when moving down the tree. Therefore, in order to concentrate on the biologically more meaningful pathway collections, an arbitrary cut-off of a maximum of 150 genes per pathway was chosen as the largest reported gene list. In addition, all pathways referring to cellular localization were removed.
List of differentially expressed genes
Expression data was filtered with stringent criteria to reveal only the most significant differentially expressed probes from the dataset. Every probe with flag status 'present' was accepted and the data were further filtered by standard deviation (95% filtered out). Statistical analysis was performed with 2-way analysis of variance (ANOVA), using R-statistics, and taking into account the repeated measures from the same subjects (Additional file 1). All probes with p < 0.05 after Benjamini-Hochberg multiple testing correction  were accepted. Two supplement tables were generated: one for significantly regulated probes between the groups (insulin-resistant vs. insulin-sensitive) and a second one for significantly regulated probes by insulin treatment (fasting vs. hyperinsulinemia).
Quantitative RT-PCR (qPCR) was used to measure the relative abundance of interesting transcripts found from the microarray analysis. Samples from 20 women, including the 10 samples used for the microarray analysis, at three time points (0 h, 3 h and 6 h) were analyzed. Reverse transcription was performed using M-MLV reverse transcriptase (Invitrogen) using oligo(dT)12-18 primers. Power SYBR® Green PCR Master Mix and ABI 7900 HT system (Applied biosystems) were used for quantitative PCR analysis with gene specific primers (Additional file 2: Table S1Q1) and 1 ng of template in a final reaction volume of 10 μl. Genorm algorithm  was used to validate three different control genes for normalization. Data were normalized to geometric mean of two control genes, beta actin (ACTB) and large ribosomal protein 0 (RPLP0).
Clinical and biochemical characteristics in the insulin-sensitive and the insulin-resistant groups were compared using the Mann-Whitney U test and 2-way ANOVA. The calculations were performed with SPSS statistics 18.0 for Windows (SPSS, Chicago, IL). All data are shown as mean ± standard error of mean.
Clinical characteristics of the study subjects
The clinical characteristics of the study subjects are summarized in Table 1. The characteristics are reported separately for subjects in the microarray study and for the subjects in the qRT-PCR experiment. Markers of insulin resistance, serum fasting insulin and C-peptide concentrations were higher, and HDL-cholesterol concentration lower in the insulin-resistant group. During insulin infusion, serum insulin concentrations were similar between the lean insulin-sensitive and the obese insulin-resistant groups (68 ± 6 vs. 84 ± 10 mU/l, respectively; NS). By definition, whole-body insulin sensitivity was 131% higher in the insulin-sensitive than the insulin-resistant group in the microarray experiment (n = 10) (9.7 ± 0.5 vs. 4.2 ± 0.4 mg kg - 1 min - 1, p < 0.009) and 107% higher in the entire study group (n = 20) (8.7 ± 0.4 vs. 4.2 ± 0.4 mg kg 1 min 1, p < 0.0001). Moreover, all indicators of obesity (body weight, BMI, whole-body-fat, fat mass, and waist-to-hip ratio were significantly higher in the insulin-resistant subjects than in the sensitive ones.
Inflammatory pathways are up-regulated in the insulin-resistant women
Activation of inflammatory pathways was a predominant phenotype in the insulin-resistant group after overnight fast (Table 2). The inflammatory gene expression pattern consisted mainly of two aspects of the immune response: genes involved in complement activation (inflammatory response, GO:0006954, p = 2.0E-10) and those involved in cellular chemotactic activity (taxis/chemotaxis, GO:0006935, p = 1.8E-07). The inflammatory response pathway included complement 1q subcomponent, B chain (C1QB), A chain (C1QA) and C chain (C1QC). In this pathway, the gene most strongly elevated in the insulin-resistant group was LDL-associated phospholipase A2 (PLA2G7 or LDL-PLA2) (Figure 1A).
The second aspect of the inflammatory phenotype involved a number of chemokines and adhesion molecules in the taxis/chemotaxis pathways. The most prominently elevated genes in the insulin-resistant subjects were the spleen tyrosine kinase (SYK), chemokine CCL13, and beta 2 integrin (ITGB2) involved in leukocyte adhesion (Figure 1B). Up-regulation of inflammatory pathways similar to those in the fasted situation was observed between the groups during hyperinsulinemia: The most divergent pathways here included inflammatory response (GO:0006954, p = 3.7E-07) and regulation of immune system process (GO:0002682, p = 1.0E-06) (Table 3). Expression of marcrophage marker CD36 was 14.2% increased in the insulin-resistant group.
ANOVA analysis of single genes expressed differentially between the groups revealed additional genes of interest such as osteopontin (SPP1), matrix metalloproteinase 9 (MMP9), branched chain amino-acid transaminase 1 (BCAT1) and EGF-like protein 6 (EGFL6). These were expressed at higher levels in the insulin-resistant than the insulin-sensitive group (Additional file 3: Table S2).
Monocarboxylic acid metabolism is down-regulated in the insulin-resistant women during fasting
Expression of genes involved in monocarboxylic acid metabolism (GO:0032787, p = 5.6E-06), most of which represent the fatty acid metabolism pathway (GO:0006631, p = 5.8E-05), was reduced in the insulin-resistant group as compared with the insulin-sensitive group after an overnight fast (Table 4). Fatty acid transporter 2 (SLC27A2 or FATP2) was the most differentially expressed gene in this pathway (Figure 2). In addition, decreased mRNA expression of the oxidation reduction (GO:0055114) and cellular alcohol metabolic (GO:0006066) pathways was detected in the insulin-resistant group (Table 4).
Defective mRNA expression of nuclear mitochondrial genes in insulin-resistant women during hyperinsulinemia
Reduced mRNA expression of mitochondrial pathways was a predominant phenotype in the insulin-resistant group during hyperinsulinemia (Table 5). Especially genes in the mitochondrial respiratory chain (GO: 0005746, p = 9.7E-06) pathway were expressed at markedly reduced levels in the insulin-resistant group as compared with the insulin-sensitive group. The genes were mostly those involved in mitochondrial complex I function (Figure 3). In addition to genes encoding mitochondrial complexes I, III, and IV, the differentially expressed genes in this pathway included fatty acid desaturase 1 (FADS1) and peroxisome proliferator-activated receptor gamma, coactivator 1 alpha (PPARGC1) transcribed at reduced levels in the insulin-resistant women.
Induction of lipid metabolism pathways by insulin in the two study groups
We next investigated the pathways most prominently induced by insulin (Tables 6 &7). The ten most up-regulated pathways were predominantly those involved in lipid metabolism, in both the insulin-resistant and the insulin-sensitive groups. More detailed analysis revealed that the lipid metabolic process (GO:0006629) and sterol metabolic process (GO:0016125) pathways, despite that fact that they were induced by insulin in both groups, were expressed at lower levels in the insulin-resistant as compared with the insulin-sensitive group during hyperinsulinemia (Figure 4A &4B). Fatty acid binding protein 2 (FATP2 or SLC27A2), fatty acid elongation factor 6 (ELOVL6), and patatin-like phospholipase domain containing 3 (PNPLA3) were among the most divergently expressed genes in the GO:0006629 pathway. The sterol metabolic process pathway contained several genes with well established roles in cholesterol metabolism, including 3-hydroxy-3-methylglutaryl-Coenzyme A reductase (HMGCR) and sterol regulatory element binding transcription factor 1 (SREBF1) (Figure 4B).
ANOVA analysis of single insulin-regulated genes revealed a number of additional interesting lipid metabolism-related transcripts whose levels differed between the insulin-sensitive and -resistant groups (Additional file 4: Table S3Q4). One of these, angiopoietin-like protein 4 (ANGPTL4 or FIAF), is currently recognized as an important regulator of lipoprotein lipase function, the expression of which is strongly induced in mouse adipose tissue during fasting . This transcript was elevated during hyperinsulinemia in the insulin-resistant as compared to the insulin-sensitive group.
Quantitative RT-PCR analysis of mRNA levels of selected genes altered by insulin or insulin resistance status
To verify some of the most interesting findings of the microarray study, we performed qPCR quantification of mRNAs for 8 genes, for 11 insulin-sensitive and 9 insulin-resistant women, under fasting and at the 3 h and 6 h time points of hyperinsulinemia. In all cases, the microarray findings were confirmed. FATP2 mRNA remained at lower levels in the insulin-resistant group during fasting (p < 0.01) and the levels were not affected by hyperinsulinemia in either group (Figure 5A). ELOVL6 was strongly induced by insulin only in the insulin-sensitive group (p < 0.001) resulting in strikingly decreased mRNA levels in the insulin-resistant compared to insulin-sensitive group (p < 0.01) at the 6 hour time point (Figure 5B). ANGPTL4 was strongly down-regulated by insulin in the insulin-sensitive group (p = < 0.001) whereas this regulation was almost absent in the insulin-resistant group (p = NS). As a consequence we observed higher ANGPTL4 expression in the insulin-resistant compared to the insulin-sensitive group after 6 hours of hyperinsulinemia (Figure 5C). Furthermore, significantly higher mRNA levels of ANGPTL4 in the insulin-sensitive group were observed during fasting in analysis of the 5 most insulin-resistant and 5 most insulin-sensitive persons analyzed in the microarray experiment (p < 0.05, data not shown). PNPLA3 expression was lower in the insulin-resistant compared to the insulin-sensitive group during hyperinsulinemia (p < 0.001 at 3 h and p < 0.001 at 6 h) and the expression was increased by insulin in both groups at 6 h time point and in the insulin-sensitive group already at 3 h time point (p < 0.001) (Figure 5D). Expression of SREBF1 and its target HMGCR displayed significantly lower mRNA levels in the insulin-resistant group after 3 h (p < 0.05 and p < 0.001) and 6 h (p < 0.05 for both genes) of hyperinsulinemia (Figure 5E &5F). The expression of both SREBF1 (p < 0.001 at 6 h for insulin-sensitive, p < 0.01 for insulin-resistant) and HMGCR (p < 0.001 at 6 h for the insulin-sensitive, p < 0.05 for the insulin-resistant) in both study groups was up-regulated by insulin, although the activation was clearly weaker in the insulin-resistant group. Expression of three nuclear mitochondrial genes, NDUFS4, UQCR and COX4I2, was verified by RT-PCR (Additional file 5). All have lower expression in insulin-resistant as compared to insulin-sensitive subjects, consistently at every time point (p < 0.05).
In the present study we compared the acute in vivo effects of insulin on adipose tissue transcriptional profiles of obese insulin-resistant and lean insulin-sensitive women, to gain insight into the molecular mechanisms underlying insulin resistance. The most striking difference in gene expression of the insulin-resistant group during hyperinsulinemia was reduced transcription of genes involved in mitochondrial respiration (mitochondrial respiratory chain, GO:0001934). Inflammatory pathways with complement components (inflammatory response, GO:0006954) and cytokines (chemotaxis, GO:0042330) were strongly up-regulated in insulin-resistant as compared to insulin-sensitive subjects both before and during hyperinsulinemia. Furthermore, differences were observed in genes contributing to fatty acid, cholesterol and triglyceride metabolism and in genes involved in regulating lipolysis, between the insulin-resistant and -sensitive subjects especially during hyperinsulinemia.
Low-grade adipose tissue inflammation has been postulated as one of the central factors in the development of insulin resistance. Increased inflammation was a predominant phenotype in the insulin-resistant vs. insulin-sensitive subjects in the present study. Many findings of the up-regulated genes in the insulin-resistant group, such as C1Q peptides, MMP9 and SPP1, are consistent with previous reports on adipose tissue gene expression in the transcriptomics analysis in the fasting state in obese mice [24, 25], and in monozygotic twins discordant for BMI . While these data characterize changes in adipose tissue transcriptome in humans and are unique as such, the role of complement components and inflammatory mediators in the development of insulin resistance cannot be determined in a cross-sectional study in humans. Since many of the inflammatory genes observed in this study are active in cells of the monocyte/macrophage lineage, our observation may reflect increased inflammatory cell content in adipose tissue of the insulin-resistant women. This is in line with the observation that the adipose tissue of lean subjects usually consists of approximately 5-10% of macrophages, whereas in obese patients, adipose tissue macrophage content can be as high as 50% of the total cell number . However, there is increasing evidence that some of the genes we found up-regulated are expressed also by adipocytes, such as complement components, PLA2G7 [27, 28], SYK  and ITGB2. Adipocyte and macrophage trancriptomes are similar and they can become even more similar after macrophages engulf lipids  from the dying adipocytes observed in obesity .
Alterations in fatty acid handling and release characterize adipose tissue in obesity. In the present study one of the most notable changes in adipose tissue transcriptome of the insulin-resistant women during fasting was down-regulation of the monocarboxylic acid metabolism pathway, FATP2 being the most down-regulated gene in this pathway. There are no previous reports describing FATP2 expression in human adipose tissue. Moreover, FATP2 is not expressed in mouse 3T3-L1 adipocytes  and in murine tissues it is expressed most strongly in liver and kidney cortex (data on adipose tissue not available in this study) . It seems unlikely that the FATP2 signal in our study would originate from macrophages, since its transcription is strongly down-regulated in the obese insulin-resistant group where macrophage numbers are elevated . Interestingly, a recent report in rat peripheral blood mononuclear cells suggested FATP2 to be an early marker of obesity .
ELOVL6 is a fatty acid elongation factor specific for long chain fatty acids . It which was recently shown to be regulated by SREBP1c . In the present study we made a novel observation of ELOVL6 down-regulation in the obese insulin resistant as compared to the lean insulin-sensitive group in hyperinsulinemia. Interestingly, liver deficiency of ELOVL6 significantly ameliorates insulin resistance in mice by modifying hepatic fatty acid composition . PNPLA3 is expressed mainly in liver and adipose tissue and its genetic variants associate in multiple studies with increased hepatic fat. PNPLA3 has been shown to have a lipolytic and weak lipogenic effect in vivo but its precise role in vivo is unclear . Its expression in adipose tissue has been reported to be similar [38, 39] or increased  in obese compared to lean patients at fasting state. Both insulin and glucose stimulate adipose tissue PNPLA3 expression . The dramatically lower transcript levels of ELOVL6 as well as PNPLA3 in hyperinsulinemic insulin-resistant subjects might have an impact to adipocyte lipid composition and could further decline adipose tissue function in insulin resistance.
Mitochondrial pathways, especially genes involved in mitochondrial respiration, have been shown to be down-regulated in muscle and adipose tissues of insulin-resistant and type 2 diabetic subjects in fasting state [41, 42]. Acquired obesity and poor physical fitness are known to impair the expression of genes of oxidative phosphorylation (41). However, their response to insulin in insulin resistant subjects has not been reported before. Mootha et al. (41) presented evidence that a number of genes involved in oxidative phosphorylation in skeletal muscle are subject to regulation by PGC-1α encoded by PPARGC1, and are down-regulated in type 2 diabetes. Therefore, the down-regulation of PPARGC1 we observed in the subcutaneous fat of insulin resistant subjects may provide one mechanistic explanation to the reduced expression of mitochondrial respiratory chain genes. Interestingly, decreased expression of mitochondrial pathways is the most prominent finding during hyperinsulinemia reflecting a regulatory defect that may further aggravate the pathogenesis of insulin resistance. Since regular practice of physical exercise is known to improve insulin sensitivity and reduce body weight [43, 44], it is possible that differences in physical activity could have amplified the observed differences in gene expression between the insulin sensitive and insulin resistant subjects.
Activity of lipoprotein lipase (LPL) is an important determinant in the development of obesity in mouse models. As a general rule, high fat diet-induced adipogenesis is aggravated by stimulated LPL activity (e.g. by adipose tissue-specific over expression of LPL or deficiency for APOCIII), and attenuated by inhibited LPL activity . A physiologically important LPL inhibitor, ANGPTL4, is strongly down-regulated by insulin in mouse 3T3-L1 cells . In the present study we found insulin to significantly decrease ANGPTL4 expression in human adipose tissue. This finding is novel as was the finding of impaired acute insulin regulation of ANGPTL4 in obese insulin resistant as compared to lean insulin sensitive subjects. In mice, ANGPTL4 over expression inhibits LPL, which in turn slows down adipogenesis and increases plasma triglyceride concentrations. When ANGPTL4 is deleted the reverse phenotype arises [23, 45, 47]. Interestingly, serum ANGPTL4 levels are inversely correlated with plasma glucose concentrations and the serum levels are significantly lower in type 2 diabetic patients than healthy subjects . ANGPTL4 expression is also stimulated by insulin sensitizing drugs thiazolidinediones via PPARy mediated mechanism suggesting that changes in ANGPTL level could have a role in the development of insulin-resistance .
The insulin-resistant group displayed defective induction of important insulin-induced lipid metabolism pathways. While several genes in these pathways are important players in cholesterol metabolism and have a well established role in cholesterol synthesis in liver, little is known of their importance in adipose tissue. Adipose tissue contains the body's largest pool of free cholesterol and cholesterol imbalance is recognized as a characteristic of enlarged adipocytes in obesity . Therefore, the differentially expressed genes in the sterol metabolic pathway in our data could play an important part in the insulin-resistant phenotype. The molecular mechanism underlying the decreased expression of lipid biosynthesis pathways in the insulin-resistant group could involve SREBP-1c since several of the genes shown in Figure 4 are targets of SREBP-1c including HMGCR and ELOVL6 [35, 51]. SREBP-1c is regulated mainly at the transcriptional level and is a major mediator of fatty acid and triglyceride synthesis induction by insulin in liver and adipose tissue [51–54]. Down-regulation of SREBP-1c was observed in the adipose tissue of leptin deficient mice (ob/ob) , in insulin-resistant humans during fasting and in response to insulin , and in adipose tissue of patients with type 2 diabetes , similar to our study. In contrast to the latter study, we detected no difference in SREBF1 expression between the groups in the basal state. The small sample size of the study could have contributed to failure to detect such a difference and also other small scale changes in gene expression. The defective response of SREBF1 to insulin in the insulin-resistant group could explain down-regulation of ELOVL6 and HMGCR mRNAs in this group. Reduced SREBF1 transcription could be a homeostatic reaction of enlarged adipocytes to prevent further lipid synthesis. On the other hand, the reduced SREBF1 expression could also decrease lipogenesis in smaller adipocytes and hinder adipocyte differentiation via PPARG down-regulation, leading to increased insulin resistance .
In conclusion, our data demonstrate that insulin-resistant subjects have a marked decrease in the expression of numerous genes involved in lipid and mitochondrial metabolism in human adipose tissue. This study also replicates previous findings of increased inflammation in adipose tissue of insulin-resistant human subjects. Insulin resistance is a complex phenotype characterized by disturbances in the transcriptional networks controlling immune response, lipid metabolism and mitochondrial function. While this study uncovers several novel aberrations in the obese and insulin-resistant human subjects' adipose tissue, the ultimate trigger of the transcriptional changes leading to adipose tissue dysfunction and insulin resistance remains to be identified.
Guilherme A, Virbasius JV, Puri V, Czech MP: Adipocyte dysfunctions linking obesity to insulin resistance and type 2 diabetes. Nat Rev Mol Cell Biol. 2008, 9 (5): 367-377. 10.1038/nrm2391.
Kahn SE, Hull RL, Utzschneider KM: Mechanisms linking obesity to insulin resistance and type 2 diabetes. Nature. 2006, 444 (7121): 840-846. 10.1038/nature05482.
Sun G: Application of DNA microarrays in the study of human obesity and type 2 diabetes. OMICS. 2007, 11 (1): 25-40. 10.1089/omi.2006.0003.
Baranova A, Collantes R, Gowder SJ, Elariny H, Schlauch K, Younoszai A, King S, Randhawa M, Pusulury S, Alsheddi T, et al: Obesity-related differential gene expression in the visceral adipose tissue. Obes Surg. 2005, 15 (6): 758-765. 10.1381/0960892054222876.
Clement K, Viguerie N, Poitou C, Carette C, Pelloux V, Curat CA, Sicard A, Rome S, Benis A, Zucker JD, et al: Weight loss regulates inflammation-related genes in white adipose tissue of obese subjects. FASEB J. 2004, 18 (14): 1657-1669. 10.1096/fj.04-2204com.
Dahlman I, Linder K, Arvidsson Nordstrom E, Andersson I, Liden J, Verdich C, Sorensen TI, Arner P: Changes in adipose tissue gene expression with energy-restricted diets in obese women. Am J Clin Nutr. 2005, 81 (6): 1275-1285.
Gabrielsson BG, Johansson JM, Lonn M, Jernas M, Olbers T, Peltonen M, Larsson I, Lonn L, Sjostrom L, Carlsson B, et al: High expression of complement components in omental adipose tissue in obese men. Obes Res. 2003, 11 (6): 699-708. 10.1038/oby.2003.100.
Hardy OT, Perugini RA, Nicoloro SM, Gallagher-Dorval K, Puri V, Straubhaar J, Czech MP: Body mass index-independent inflammation in omental adipose tissue associated with insulin resistance in morbid obesity. Surg Obes Relat Dis.
Henegar C, Tordjman J, Achard V, Lacasa D, Cremer I, Guerre-Millo M, Poitou C, Basdevant A, Stich V, Viguerie N, et al: Adipose tissue transcriptomic signature highlights the pathological relevance of extracellular matrix in human obesity. Genome Biol. 2008, 9 (1): R14-10.1186/gb-2008-9-1-r14.
Lee YH, Nair S, Rousseau E, Allison DB, Page GP, Tataranni PA, Bogardus C, Permana PA: Microarray profiling of isolated abdominal subcutaneous adipocytes from obese vs non-obese Pima Indians: increased expression of inflammation-related genes. Diabetologia. 2005, 48 (9): 1776-1783. 10.1007/s00125-005-1867-3.
Marrades MP, Milagro FI, Martinez JA, Moreno-Aliaga MJ: Differential expression of aquaporin 7 in adipose tissue of lean and obese high fat consumers. Biochem Biophys Res Commun. 2006, 339 (3): 785-789. 10.1016/j.bbrc.2005.11.080.
Pietilainen KH, Naukkarinen J, Rissanen A, Saharinen J, Ellonen P, Keranen H, Suomalainen A, Gotz A, Suortti T, Yki-Jarvinen H, et al: Global transcript profiles of fat in monozygotic twins discordant for BMI: pathways behind acquired obesity. PLoS Med. 2008, 5 (3): e51-10.1371/journal.pmed.0050051.
Vohl MC, Sladek R, Robitaille J, Gurd S, Marceau P, Richard D, Hudson TJ, Tchernof A: A survey of genes differentially expressed in subcutaneous and visceral adipose tissue in men. Obes Res. 2004, 12 (8): 1217-1222. 10.1038/oby.2004.153.
Yang X, Jansson PA, Nagaev I, Jack MM, Carvalho E, Sunnerhagen KS, Cam MC, Cushman SW, Smith U: Evidence of impaired adipogenesis in insulin resistance. Biochem Biophys Res Commun. 2004, 317 (4): 1045-1051. 10.1016/j.bbrc.2004.03.152.
Hotamisligil GS, Erbay E: Nutrient sensing and inflammation in metabolic diseases. Nat Rev Immunol. 2008, 8 (12): 923-934. 10.1038/nri2449.
Dahlman I, Forsgren M, Sjogren A, Nordstrom EA, Kaaman M, Naslund E, Attersand A, Arner P: Downregulation of electron transport chain genes in visceral adipose tissue in type 2 diabetes independent of obesity and possibly involving tumor necrosis factor-alpha. Diabetes. 2006, 55 (6): 1792-1799. 10.2337/db05-1421.
Westerbacka J, Corner A, Kannisto K, Kolak M, Makkonen J, Korsheninnikova E, Nyman T, Hamsten A, Fisher RM, Yki-Jarvinen H: Acute in vivo effects of insulin on gene expression in adipose tissue in insulin-resistant and insulin-sensitive subjects. Diabetologia. 2006, 49 (1): 132-140. 10.1007/s00125-005-0075-5.
Westerbacka J, Corner A, Kolak M, Makkonen J, Turpeinen U, Hamsten A, Fisher RM, Yki-Jarvinen H: Insulin regulation of MCP-1 in human adipose tissue of obese and lean women. Am J Physiol Endocrinol Metab. 2008, 294 (5): E841-E845. 10.1152/ajpendo.00653.2006.
DeFronzo RA, Tobin JD, Andres R: Glucose clamp technique: a method for quantifying insulin secretion and resistance. Am J Physiol. 1979, 237 (3): E214-E223.
Yki-Jarvinen H, Nikkila EA, Kubo K, Foley JE: Assay of glucose transport in human fat cells obtained by needle biopsy. Diabetologia. 1986, 29 (5): 287-290. 10.1007/BF00452064.
Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: a Practical and Powerful Approach to Multiple Testing. J Roy Statistical Society B. 1995, 289-300. 57
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3 (7): RESEARCH0034.
Mandard S, Zandbergen F, van Straten E, Wahli W, Kuipers F, Muller M, Kersten S: The fasting-induced adipose factor/angiopoietin-like protein 4 is physically associated with lipoproteins and governs plasma lipid levels and adiposity. J Biol Chem. 2006, 281 (2): 934-944.
Kiefer FW, Zeyda M, Gollinger K, Pfau B, Neuhofer A, Weichhart T, Saemann MD, Geyeregger R, Schlederer M, Kenner L, et al: Neutralization of osteopontin inhibits obesity-induced inflammation and insulin resistance. Diabetes. 59 (4): 935-946.
Zhang J, Wright W, Bernlohr DA, Cushman SW, Chen X: Alterations of the classic pathway of complement in adipose tissue of obesity and insulin resistance. Am J Physiol Endocrinol Metab. 2007, 292 (5): E1433-E1440. 10.1152/ajpendo.00664.2006.
Weisberg SP, McCann D, Desai M, Rosenbaum M, Leibel RL, Ferrante AW: Obesity is associated with macrophage accumulation in adipose tissue. J Clin Invest. 2003, 112 (12): 1796-1808.
Gountopoulou A, Leondaritis G, Galanopoulou D, Mavri-Vavayanni M: TNFalpha is a potent inducer of platelet-activating factor synthesis in adipocytes but not in preadipocytes. Differential regulation by PI3K. Cytokine. 2008, 41 (2): 174-181. 10.1016/j.cyto.2007.11.008.
Nair S, Lee YH, Rousseau E, Cam M, Tataranni PA, Baier LJ, Bogardus C, Permana PA: Increased expression of inflammation-related genes in cultured preadipocytes/stromal vascular cells from obese compared with non-obese Pima Indians. Diabetologia. 2005, 48 (9): 1784-1788. 10.1007/s00125-005-1868-2.
Wang H, Malbon CC: G(s)alpha repression of adipogenesis via Syk. J Biol Chem. 1999, 274 (45): 32159-32166. 10.1074/jbc.274.45.32159.
Cinti S, Mitchell G, Barbatelli G, Murano I, Ceresi E, Faloia E, Wang S, Fortier M, Greenberg AS, Obin MS: Adipocyte death defines macrophage localization and function in adipose tissue of obese mice and humans. J Lipid Res. 2005, 46 (11): 2347-2355. 10.1194/jlr.M500294-JLR200.
Sandoval A, Chokshi A, Jesch ED, Black PN, Dirusso CC: Identification and characterization of small compound inhibitors of human FATP2. Biochem Pharmacol. 79 (7): 990-999.
Hirsch D, Stahl A, Lodish HF: A family of fatty acid transporters conserved from mycobacterium to man. Proc Natl Acad Sci USA. 1998, 95 (15): 8625-8629. 10.1073/pnas.95.15.8625.
Caimari A, Oliver P, Rodenburg W, Keijer J, Palou A: Slc27a2 expression in peripheral blood mononuclear cells as a molecular marker for overweight development. Int J Obes (Lond). 34 (5): 831-839.
Matsuzaka T, Shimano H, Yahagi N, Yoshikawa T, Amemiya-Kudo M, Hasty AH, Okazaki H, Tamura Y, Iizuka Y, Ohashi K, et al: Cloning and characterization of a mammalian fatty acyl-CoA elongase as a lipogenic enzyme regulated by SREBPs. J Lipid Res. 2002, 43 (6): 911-920.
Kumadaki S, Matsuzaka T, Kato T, Yahagi N, Yamamoto T, Okada S, Kobayashi K, Takahashi A, Yatoh S, Suzuki H, et al: Mouse Elovl-6 promoter is an SREBP target. Biochem Biophys Res Commun. 2008, 368 (2): 261-266. 10.1016/j.bbrc.2008.01.075.
Matsuzaka T, Shimano H, Yahagi N, Kato T, Atsumi A, Yamamoto T, Inoue N, Ishikawa M, Okada S, Ishigaki N, et al: Crucial role of a long-chain fatty acid elongase, Elovl6, in obesity-induced insulin resistance. Nat Med. 2007, 13 (10): 1193-1202. 10.1038/nm1662.
Romeo S, Huang-Doran I, Baroni MG, Kotronen A: Unravelling the pathogenesis of fatty liver disease: patatin-like phospholipase domain-containing 3 protein. Curr Opin Lipidol. 21 (3): 247-252.
Liu YM, Moldes M, Bastard JP, Bruckert E, Viguerie N, Hainque B, Basdevant A, Langin D, Pairault J, Clement K: Adiponutrin: A new gene regulated by energy balance in human adipose tissue. J Clin Endocrinol Metab. 2004, 89 (6): 2684-2689. 10.1210/jc.2003-031978.
Moldes M, Beauregard G, Faraj M, Peretti N, Ducluzeau PH, Laville M, Rabasa-Lhoret R, Vidal H, Clement K: Adiponutrin gene is regulated by insulin and glucose in human adipose tissue. Eur J Endocrinol. 2006, 155 (3): 461-468. 10.1530/eje.1.02229.
Johansson LE, Hoffstedt J, Parikh H, Carlsson E, Wabitsch M, Bondeson AG, Hedenbro J, Tornqvist H, Groop L, Ridderstrale M: Variation in the adiponutrin gene influences its expression and associates with obesity. Diabetes. 2006, 55 (3): 826-833. 10.2337/diabetes.55.03.06.db05-1075.
Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, Puigserver P, Carlsson E, Ridderstrale M, Laurila E, et al: PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003, 34 (3): 267-273. 10.1038/ng1180.
Mustelin L, Pietilainen KH, Rissanen A, Sovijarvi AR, Piirila P, Naukkarinen J, Peltonen L, Kaprio J, Yki-Jarvinen H: Acquired obesity and poor physical fitness impair expression of genes of mitochondrial oxidative phosphorylation in monozygotic twins discordant for obesity. Am J Physiol Endocrinol Metab. 2008, 295 (1): E148-E154. 10.1152/ajpendo.00580.2007.
Li G, Zhang P, Wang J, Gregg EW, Yang W, Gong Q, Li H, Jiang Y, An Y, Shuai Y, et al: The long-term effect of lifestyle interventions to prevent diabetes in the China Da Qing Diabetes Prevention Study: a 20-year follow-up study. Lancet. 2008, 371 (9626): 1783-1789. 10.1016/S0140-6736(08)60766-7.
Lindstrom J, Ilanne-Parikka P, Peltonen M, Aunola S, Eriksson JG, Hemio K, Hamalainen H, Harkonen P, Keinanen-Kiukaanniemi S, Laakso M, et al: Sustained reduction in the incidence of type 2 diabetes by lifestyle intervention: follow-up of the Finnish Diabetes Prevention Study. Lancet. 2006, 368 (9548): 1673-1679. 10.1016/S0140-6736(06)69701-8.
Voshol PJ, Rensen PC, van Dijk KW, Romijn JA, Havekes LM: Effect of plasma triglyceride metabolism on lipid storage in adipose tissue: Studies using genetically engineered mouse models. Biochim Biophys Acta. 2009
Yamada T, Ozaki N, Kato Y, Miura Y, Oiso Y: Insulin downregulates angiopoietin-like protein 4 mRNA in 3 T3-L1 adipocytes. Biochem Biophys Res Commun. 2006, 347 (4): 1138-1144. 10.1016/j.bbrc.2006.07.032.
Okubo M, Horinishi A, Saito M, Ebara T, Endo Y, Kaku K, Murase T, Eto M: A novel complex deletion-insertion mutation mediated by Alu repetitive elements leads to lipoprotein lipase deficiency. Mol Genet Metab. 2007, 92 (3): 229-233. 10.1016/j.ymgme.2007.06.018.
Xu A, Lam MC, Chan KW, Wang Y, Zhang J, Hoo RL, Xu JY, Chen B, Chow WS, Tso AW, et al: Angiopoietin-like protein 4 decreases blood glucose and improves glucose tolerance but induces hyperlipidemia and hepatic steatosis in mice. Proc Natl Acad Sci USA. 2005, 102 (17): 6086-6091. 10.1073/pnas.0408452102.
Yoon JC, Chickering TW, Rosen ED, Dussault B, Qin Y, Soukas A, Friedman JM, Holmes WE, Spiegelman BM: Peroxisome proliferator-activated receptor gamma target gene encoding a novel angiopoietin-related protein associated with adipose differentiation. Mol Cell Biol. 2000, 20 (14): 5343-5349. 10.1128/MCB.20.14.5343-5349.2000.
Yu BL, Zhao SP, Hu JR: Cholesterol imbalance in adipocytes: a possible mechanism of adipocytes dysfunction in obesity. Obes Rev. 2009
Le Lay S, Lefrere I, Trautwein C, Dugail I, Krief S: Insulin and sterol-regulatory element-binding protein-1c (SREBP-1C) regulation of gene expression in 3T3-L1 adipocytes. Identification of CCAAT/enhancer-binding protein beta as an SREBP-1C target. J Biol Chem. 2002, 277 (38): 35625-35634. 10.1074/jbc.M203913200.
Nadeau KJ, Leitner JW, Gurerich I, Draznin B: Insulin regulation of sterol regulatory element-binding protein-1 expression in L-6 muscle cells and 3 T3 L1 adipocytes. J Biol Chem. 2004, 279 (33): 34380-34387. 10.1074/jbc.M403596200.
Repa JJ, Liang G, Ou J, Bashmakov Y, Lobaccaro JM, Shimomura I, Shan B, Brown MS, Goldstein JL, Mangelsdorf DJ: Regulation of mouse sterol regulatory element-binding protein-1c gene (SREBP-1c) by oxysterol receptors, LXRalpha and LXRbeta. Genes Dev. 2000, 14 (22): 2819-2830. 10.1101/gad.844900.
Shimano H, Horton JD, Hammer RE, Shimomura I, Brown MS, Goldstein JL: Overproduction of cholesterol and fatty acids causes massive liver enlargement in transgenic mice expressing truncated SREBP-1a. J Clin Invest. 1996, 98 (7): 1575-1584. 10.1172/JCI118951.
Yahagi N, Shimano H, Hasty AH, Matsuzaka T, Ide T, Yoshikawa T, Amemiya-Kudo M, Tomita S, Okazaki H, Tamura Y, et al: Absence of sterol regulatory element-binding protein-1 (SREBP-1) ameliorates fatty livers but not obesity or insulin resistance in Lep(ob)/Lep(ob) mice. J Biol Chem. 2002, 277 (22): 19353-19357. 10.1074/jbc.M201584200.
Ducluzeau PH, Perretti N, Laville M, Andreelli F, Vega N, Riou JP, Vidal H: Regulation by insulin of gene expression in human skeletal muscle and adipose tissue. Evidence for specific defects in type 2 diabetes. Diabetes. 2001, 50 (5): 1134-1142. 10.2337/diabetes.50.5.1134.
Sewter C, Berger D, Considine RV, Medina G, Rochford J, Ciaraldi T, Henry R, Dohm L, Flier JS, O'Rahilly S, et al: Human obesity and type 2 diabetes are associated with alterations in SREBP1 isoform expression that are reproduced ex vivo by tumor necrosis factor-alpha. Diabetes. 2002, 51 (4): 1035-1041. 10.2337/diabetes.51.4.1035.
Kim JB, Spiegelman BM: ADD1/SREBP1 promotes adipocyte differentiation and gene expression linked to fatty acid metabolism. Genes Dev. 1996, 10 (9): 1096-1107. 10.1101/gad.10.9.1096.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/5/9/prepub
We gratefully acknowledge Ms Lea Puhakka, Anne Vikman, Kaisa Makkonen, Katja Tuominen and Mia Urjansson from the University of Helsinki for excellent technical assistance. Main investigators contributing to the original clinical study, Jukka Westerbacka and Janne Makkonen from the University of Helsinki, are also greatly acknowledged. The study was supported by grants from the Finnish Foundation for Cardiovascular Research (M.J., J.S.), Academy of Finland (H. Y-J.), the Academy of Finland Centre of Excellence in Complex Disease Genetics (J.S., P-P.L., J.N.), the Liv och Hälsa Foundation (V.M.O.), the Sigrid Juselius foundation (H.Y-J), EVO foundations (H. Y-J.), Paulo foundation (J.S.), Ida Montin foundation (J.S.), Jalmari ja Rauha Ahokas foundation (J.S) and The Diabetes Research Foundation of Finland (J.S.).
We are indebted to late Professor Leena Peltonen for her enthusiasm, expertise and support in initiating, supervising and funding the genome wide expression analyses in this study.
The authors declare that they have no competing interests.
JS researched data and drafted the manuscript. PPL, JN, IS and SR researched data and edited the manuscript and contributed to discussion. MJ and VMO edited the manuscript and contributed to discussion. HJ designed the clinical study, edited the manuscript and contributed to discussion. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 3: 2-way ANOVA of differentially expressed genes associated with insulin resistance status. All probes with p-value < 0.05 after Benjamini-Hochberg correction are reported in ascending order. Non-corrected p-values are reported. Fold change (FC) and direction of regulation in insulin-resistant compared to insulin-sensitive group are reported in fasting state and during hyperinsulinemia. (DOC 200 KB)
Additional file 4: 2-way ANOVA of differentially expressed genes associated with hyperinsulinemia. All probes with p-value < 0.05 after Benjamini-Hochberg correction are reported in ascending order. Non-corrected p-values are reported. Fold change (FC) and direction of regulation in insulin-resistant compared to insulin-sensitive group are reported in fasting state and during hyperinsulinemia. (DOC 150 KB)
Additional file 5: Relative mRNA levels of selected genes in mitochondrial respiratory chain pathway. Expression levels were analyzed in insulin-resistant (black bar) and insulin-sensitive subjects (white bar) at 0 h, 3 h, and 6 h during euglycaemic hyperinsulinaemia. All genes are differentially expressed between insulin-resistant and insulin-sensitive groups: p(group) < 0.05, p(group*insulin) = NS, p(insulin) = NS for all genes via 2-way ANOVA. (PPT 146 KB)
About this article
Cite this article
Soronen, J., Laurila, PP., Naukkarinen, J. et al. Adipose tissue gene expression analysis reveals changes in inflammatory, mitochondrial respiratory and lipid metabolic pathways in obese insulin-resistant subjects. BMC Med Genomics 5, 9 (2012). https://doi.org/10.1186/1755-8794-5-9