|Home | About | Journals | Submit | Contact Us | Français|
Conceived and designed the experiments: SND DJF AKS CS TH BJC VV VS GM. Performed the experiments: SND DJF AKS CS RH TH VLV JVS GM. Analyzed the data: SND DJF AKS CS RH TH VLV JVS VS GM. Contributed reagents/materials/analysis tools: VS GM. Wrote the paper: SND DJF AKS CS RH TH VLV BJC VV JVS VS GM. Collected the samples: DJF. Collected the samples: VV.
In obesity, impaired adipose tissue function may promote secondary disease through ectopic lipid accumulation and excess release of adipokines, resulting in systemic low-grade inflammation, insulin resistance and organ dysfunction. However, several of the genes regulating adipose tissue function in obesity are yet to be identified.
In order to identify novel candidate genes that may regulate adipose tissue function, we analyzed global gene expression in abdominal subcutaneous adipose tissue before and one year after bariatric surgery (biliopancreatic diversion with duodenal switch, BPD/DS) (n=16). Adipose tissue from lean healthy individuals was also analyzed (n=13). Two different microarray platforms (AB 1700 and Illumina) were used to measure the differential gene expression, and the results were further validated by qPCR. Surgery reduced BMI from 53.3 to 33.1 kg/m2. The majority of differentially expressed genes were down-regulated after profound fat loss, including transcription factors involved in stress response, inflammation, and immune cell function (e.g., FOS, JUN, ETS, C/EBPB, C/EBPD). Interestingly, a distinct set of genes was up-regulated after fat loss, including homeobox transcription factors (IRX3, IRX5, HOXA5, HOXA9, HOXB5, HOXC6, EMX2, PRRX1) and extracellular matrix structural proteins (COL1A1, COL1A2, COL3A1, COL5A1, COL6A3).
The data demonstrate a marked switch of transcription factors in adipose tissue after profound fat loss, providing new molecular insight into a dichotomy between stress response and metabolically favorable tissue development. Our findings implicate homeobox transcription factors as important regulators of adipose tissue function.
Obesity has reached epidemic proportions and is associated with increased risk of type 2 diabetes, cardiovascular diseases, several forms of cancer, and other diseases . As a potent endocrine organ as well as the body's primary energy storage reserve, adipose tissue plays key roles in systemic metabolism . The pathogenic potential of adipose tissue is conferred by changes in morphology and cellular functions that involve local inflammation, aberrant hormonal signaling and adipokine secretion, and altered lipid storage dynamics , . While adipocyte turnover appears to be tightly regulated , obesity may result from a combination of adipocyte hypertrophy and hyperplasia with varying contributions by these processes in different individuals . Recent studies have indicated that hypertrophy may promote inflammation and dysfunctional adipose tissue via hypoxia and aberrant extracellular matrix remodeling , , . Moreover, hypertrophy was found to be associated with a reduced adipocyte turnover . It has been proposed that limited recruitment and development of adipocytes, and thereby reduced ability to expand the adipose tissue during energy surplus, may promote ectopic lipid accumulation and adverse effects on other organs , , .
Previous studies have demonstrated a marked increase in the expression of pro-inflammatory genes in adipose tissue in obesity , . While isolated adipocytes show an increased inflammatory gene expression profile , increased infiltration and activation of non-fat cells, such as macrophages, may propagate local inflammation in adipose tissue ,  and promote aberrant extracellular matrix remodeling . Analyses of adipose tissue before and three months after gastric bypass surgery revealed a reduction in macrophages and altered expression of genes regulating inflammation and extracellular matrix , . Moreover, adipose tissue gene expression is responsive to dietary intervention. Energy-restricted diets have been shown to improve the inflammatory profile  and to down-regulate genes involved in extracellular matrix  and production of polyunsaturated fatty acids . It has also been demonstrated that dietary modification alters gene expression in adipose tissue irrespective of weight loss . Interestingly, differences in adipose tissue gene regulation were observed depending on individual responses to caloric restriction .
Due to their profound fat loss, morbidly obese patients undergoing bariatric surgery represent a powerful model for studying changes in adipose tissue biology. Among the different forms of bariatric surgery, the biliopancreatic diversion with duodenal switch (BPD/DS) induces the most rapid and pronounced fat loss , and loss of excess fat is especially rapid within the first year of surgery . Hence, we have analyzed global gene expression in subcutaneous abdominal adipose tissue (SAT) from 16 morbidly obese subjects before and one year after BPD/DS. To obtain a reference of normal expression values, we also included adipose tissue from 13 lean healthy controls undergoing inguinal hernia repair. The aim of the study was to identify novel candidate genes that may regulate adipose tissue function.
Each enrolled subject gave their written consent after being informed about the study and legal rights. The study was approved by the Western Norway Regional Committee for Medical Research Ethics (REK).
Anthropometric and biochemical data were recorded for patients before and one year after bariatric surgery (biliopancreatic diversion with duodenal switch) between 2003 and 2007. Biopsies of subcutaneous adipose tissue were collected from 16 patients (12 women) per-operatively and one year after surgery, and from 13 lean healthy control subjects who underwent inguinal hernia repair (6 women). Each biopsy was gently pressed flat in a sealable polyethylene bag, immediately transferred into liquid nitrogen, and stored at −80°C. Bariatric patients with Norwegian Caucasian origin and morbid obesity (BMI>45) were included. Average BMI was 53.3 kg/m2 (range 47 to 61) and average age was 39.3±10.9 (range 23 to 62 years). Fasting glucose, triglycerides, total-, HDL- and LDL cholesterol, insulin, insulin C-peptide, and CRP were measured in serum, and HbA1C was measured in EDTA blood (Table 1). Significance values of anthropometric and biochemical data were calculated using paired-samples T-test. Seven of the bariatric patients were diabetic before surgery, and one patient continued on a low dose of insulin one year after surgery (Table S1). Inclusion criteria for the healthy lean control subjects were BMI<27 and no history of disorders or drug abuse. The control subjects had an average BMI of 23.0±2.48 (range 18.5 to 26.8) and an average age of 47.6±17.1 (range 20 to 77 years). All the measured fasting serum values of the control subjects (glucose, triglycerides, total-, HDL- and LDL cholesterol, insulin, insulin C-peptide) were within the normal range, except three CRP values (6, 8 and 18 mg/L). Homeostasis model assessment (HOMA) index, a measure of insulin resistance where a value of 1 is optimal, was calculated as fasting serum insulin (mIU/L)×glucose (mmol/L)/22.5.
To avoid systematic bias due to batch/lot and processing variability, bariatric surgery samples and control samples were balanced and randomized in each of the following steps: adipose tissue homogenization, RNA extraction, labeling, and microarray hybridization. Pre- and post-surgery pairs of adipose tissue samples were handled simultaneously.
Approximately 200–300 mg of frozen adipose tissue was sliced on dry ice and transferred to 2 ml microcentrifuge tubes with round bottom (Eppendorf). One ml lysing buffer (Qiazol, Qiagen) and a 5 mm metal bead (Millipore) were added and immediately followed by homogenization in a TissueLyser (Qiagen), with shaking three times for 2 minutes at a frequency of 25 Hz. The lysates were stored at −80°C. RNA from the adipose tissue was extracted using the RNeasy Lipid Tissue Midi Kit (Qiagen) according to the manufacturer protocol, and samples were treated with the RNase-Free DNase Set (Qiagen). Amount and quality of the extracted RNA were measured by the NanoDrop® ND-1000 spectrophotometer (NanoDrop Technologies, USA) and the Agilent 2100 Bioanalyzer (Agilent Technologies, USA).
A pilot microarray study of adipose tissue from 9 of the patients before and after surgery was performed using the Applied Biosystems (AB) 1700 Expression Array system (Text S1). The main microarray experiment was performed using the Illumina iScan, which is based upon fluorescence detection of biotin-labeled cRNA. Using the Illumina TotalPrep RNA Amplification Kit (version 280508, Applied Biosystems/Ambion, USA), 300ng of total RNA from each sample was reversely transcribed, amplified and Biotin-16-UTP–labeled. The amount (15–52 µg) and quality of labeled cRNA were measured using both the NanoDrop spectrophotometer and Agilent 2100 Bioanalyzer. 750ng of biotin-labeled cRNA was hybridized to the HumanRef-8 v.3 Illumina Sentrix BeadChip according to manufacturer's instructions. The HumanRef-8 v.3 BeadChip targets approximately 24,500 annotated RefSeq transcripts and covers 18,631 unique curated genes.
Bead summary data was imported into GenomeStudio to remove control probes and to produce a text file containing the signal and detection p-values per probe for all samples. The text file was imported into J-Express Pro v.2.7, and signal intensity values were quantile normalized  and log transformed (base 2). Correspondence Analysis (CA)  and hierarchical clustering with Pearson Correlation as a distance measure were performed to look for global trends in the data. In the CA plot, the microarray data for genes and samples are projected onto a two-dimensional plane defined by the first and second principal components. The first principal component (along the x-axis) explains most of the total chi square, the second principal component explains second most of the total chi square. Samples that are close together in the plot have more similarity than samples further apart. The plots revealed sample Pre_11 as an outlier. Technical measures showed that this sample had lower yield from the labeling step than the other samples, and most of the signals were higher, including the P95, P05 and negative control signals. Therefore, the samples of patient 11 were removed from the dataset before further analysis.
Significance Analysis of Microarrays (SAM)  was used to look for differentially expressed genes. To obtain manageable datasets, differentially expressed genes were defined by fold change ≥1.5 and q-value=0. Protein ANalysis THrough Evolutionary Relationships (PANTHER) (dated July 15th 2009, www.pantherdb.org) was then used to search for over-represented functional categories among the differentially expressed genes. PANTHER describes genes in terms of biological processes, molecular functions, and pathways in hierarchies up to three levels. Each entry has an equivalent in Gene Ontology (www.geneontology.org). The Bonferroni correction for multiple testing was used in the calculation of PANTHER p-values.
The microarray data are MIAME compliant and are available in ArrayExpress with accession numbers E-TABM-862 (Illumina data) and E-TABM-864 (AB 1700 data).
Ten genes were selected for validation by quantitative real-time RT-PCR (qPCR). The analysis was performed on a subset of the samples corresponding to the AB 1700 microarray analysis (pre- versus post-surgery, n=9). cDNA synthesis was performed using the Transcriptor First Strand cDNA Synthesis Kit (Roche) (1 µg RNA per sample), followed by quantitative real-time RT-PCR with the LightCycler480 Probes Master kit and the LightCycler480 rapid thermal cycler system (Roche Applied Science). Target gene expression was quantified relative to the constitutive control gene TATA-binding protein (TBP) using specific Universal ProbeLibrary (UPL) probes and target-specific primers (Table S5). Target genes were amplified in duplex with TBP using the UPL Human TBP Gene Assay (Roche) with the exception of HOXA5, HOXB5, and IRX5 as a duplex run could not be performed for these genes. Instead, calculations were made with mean TBP values from five of the duplex runs. Efficiency of the transcripts was measured using standard curves with 15 dilutions of an adipose tissue sample with concentrated cDNA. Negative controls were prepared by replacing the mRNA template with PCR-grade H2O. RT-PCR was performed according to the optimized protocol (Roche Applied Science). Fold change in mRNA expression was calculated pair-wise for each patient using the crossing point (CP) for each sample (triplicate) and the amplification efficiency for the transcript. The mRNA expressions were found to be non-linear and are presented as geometric mean with 95% confidence intervals (CI).
Genes whose promoters contain potential binding sites for homeobox transcription factors were identified using the Gene2Promoter software provided by Genomatix (http://www.genomatix.de/). Genes with a differential expression after bariatric surgery (152 up-regulated and 469 down-regulated) were input for the analysis. The software retrieved human gene promoters for each gene, and the data were filtered in order to extract genes with one or more binding sites for selected homeobox transcription factors.
The average age of the bariatric patients before surgery and of the lean controls was 39.3±10.9 and 44.0±17.5 years, respectively. Average BMI of the bariatric patients was reduced by 20.2 kg/m2 one year after surgery (Table 1). Average BMI of the lean control subjects was 23.0 kg/m2. As expected, there was an improvement in glucose homeostasis and cardiovascular risk profile (Table 1). Six of the 16 bariatric patients had clear indications of manifest diabetes before surgery, assessed by fasting glucose values and diabetes history. Four of the patients were taking exogenous insulin prior to the surgery, and only one continued on a low dose after surgery (Table S1). Comparison of the average homeostasis model assessment (HOMA) index of the patients before and after surgery (8.6±6.41 versus 1.5±1.33, p-value=0.0003) indicated a marked reversal of insulin resistance. Average HOMA index of the control subjects was 1.3±0.65. Before surgery there was sign of a low-grade inflammation evaluated by C-reactive protein (CRP), and this was reversed after fat loss.
Microarray analysis revealed a substantial and highly consistent change in adipose tissue gene expression after surgery. Correspondence Analysis , which projects the differences in global gene expression in a two-dimensional plot, showed a clear-cut separation of the samples from the three groups (pre-surgery, post-surgery, and controls) (Figure 1). The first principal component, representing the largest component variance in the dataset (10.69%), showed that all the pre-operative samples grouped together, clearly separated from the post-operative and control samples. One exception was the post-operative sample of patient 17, which grouped together with the pre-operative samples. The second principal component, representing the second largest component variance in the dataset (6.36%), revealed that the control samples grouped together separately from the post-surgery samples.
Fat loss overall resulted in a down-regulation of genes in adipose tissue. Among the 621 differentially expressed genes using Illumina microarrays, 469 genes were down-regulated and 152 were up-regulated after surgery (Dataset S1). Strikingly, ninety-six of the 100 most differentially expressed genes showed a reduced expression after fat loss. Examples of strongly down-regulated genes included IL6 (Entrez Gene 3569), IL8 (Entrez Gene 3576), IL1B (Entrez Gene 3553), CCL2/MCP1 (Entrez Gene 6347), SOCS3 (Entrez Gene 9021), HIF1A (Entrez Gene 3091), PTGS2/COX-2 (Entrez Gene 5743) and NAMPT/VISFATIN (Entrez Gene 10135) (Table S2), which have previously been described in the context of obesity. Well-known genes involved in lipid metabolism were among the most differentially expressed single genes post-operatively, including LDLR (Entrez Gene 3949) and CH25H (Entrez Gene 9023) (down-regulated) and APOE (Entrez Gene 348), APOC1 (Entrez Gene 341) and SREBF1/SREBP-1c (Entrez Gene 6720) (up-regulated) (Table S2).
To test whether the use of medication or the presence of diabetes might affect the results, we analyzed the data excluding subjects on medication or with diabetes (n=7). We found similar fold change values as when including all subjects, showing no clear effects of diabetes or medication on the global gene expression (data not shown).
Next, we searched for over-represented functional categories among the most differentially expressed genes using Protein ANalysis THrough Evolutionary Relationships (PANTHER) (dated July 15th 2009, www.pantherdb.org). We found that the PANTHER Biological Process functional categories Immunity and defense and Signal transduction were strongly over-represented among down-regulated genes post-surgery (Figure 2). Over-represented sub-categories in the Immunity and defense category included Cytokine/chemokine mediated immunity, Macrophage mediated immunity, Granulocyte mediated immunity, and Stress response. Signal transduction sub-categories included Cell surface mediated signal transduction, Cytokine and chemokine mediated signaling pathway, and Cell communication. Interestingly, the analysis of up-regulated genes identified Developmental processes as the only PANTHER Biological Process category that was significantly over-represented post-operatively (Figure 2).
Over-represented PANTHER Molecular Function categories post-operatively included Signalling molecule, Receptor, Transcription factor, and Extracellular matrix (Figure 3). Genes belonging to Signalling molecule, Receptor, and Transcription factor were strongly over-expressed pre-surgery compared to post-operation and controls. Interestingly, our data revealed a switch in the expression of transcription factors, from a high expression of stress response related factors such as FOS (Entrez Gene 2353), JUN (Entrez Gene 3725), CEBPB (Entrez Gene 1051), CEBPD (Entrez Gene 1052), ETS1 (Entrez Gene 2113) and ETS2 (Entrez Gene 2114) in the morbidly obese state, to an increased expression of homeobox genes after surgery (Table 2 and and3,3, Table S3). The expression levels of homeobox genes (IRX3 (Entrez Gene 79191), IRX5 (Entrez Gene 10265), HOXA5 (Entrez Gene 3202), HOXA9 (Entrez Gene 3205), HOXB5 (Entrez Gene 3215), HOXC6 (Entrez Gene 3223), EMX2 (Entrez Gene 2018)) after surgery approximated the expression levels in healthy controls, except PRRX1 (Entrez Gene 5396) which was expressed at a higher level in post-surgery samples than in controls (Table 3). Finally, we found considerably higher expression levels of extracellular matrix genes (e.g. COL1A1 (Entrez Gene 1277), COL1A2 (Entrez Gene 1278), COL3A1 (Entrez Gene 1281), COL5A1 (Entrez Gene 1289), COL6A3 (Entrez Gene 1293)) post-surgery, whereas the expression levels were more similar in control and pre-surgery biopsies (Table 4, Table S3).
To validate the results, we compared the Illumina microarray analysis to the AB 1700 microarray analysis (n=9). The AB 1700 data confirmed the results (Figure S1, Dataset S2). Significance Analysis of Microarrays (SAM) identified 710 genes that were differentially expressed, of which 499 were down-regulated and 211 were up-regulated after surgery. As in the Illumina analysis, close to all of the top 100 differentially expressed genes were down-regulated after surgery (three up-regulated). Moreover, we also performed qPCR on the eight most significantly up-regulated homeobox transcription factors, in addition to COL1A1 (most up-regulated gene) and COL6A3 (recently connected with obesity). The qPCR analysis showed highly consistent results to the microarray analyses (Table 5).
Based on the strong down-regulation, and since surgery normalized the biochemical parameters, we asked whether the gene expression levels after surgery might be similar to lean controls despite a higher BMI (33 versus 23). Comparing pre-surgery samples to controls we found 861 differentially expressed genes (244 more expressed and 617 less expressed in control samples). PANTHER functional analysis showed a high degree of similarity to the pre- versus post-surgery comparison (Figure 2 and and3).3). Moreover, the comparison of post-surgery samples to controls revealed 254 differentially expressed genes (78 more expressed and 176 less expressed in control patient samples) (Dataset S3). Functional analysis showed a distinct over-representation post-operatively of the Immunity and defense sub-categories B-cell, T-cell, and natural killer cell mediated immunity, and also the category Proteolysis (Figure 2), suggesting that these processes were not normalized after surgery. Moreover, the expression of several genes involved in metabolic processes was higher in the lean controls, including genes categorized to Fatty acid metabolism, Amino acid metabolism, Carbohydrate metabolism, and Coenzyme and prosthetic group metabolism (Figure 2). Of particular note, genes encoding oxidoreductases, including dehydrogenases, were expressed at higher levels in controls than in post-surgery subjects (Figure 3, Table 6).
To further analyze the potential functions of homeobox transcription factors in adipose tissue, we performed promoter analysis of the most differentially expressed genes pre- versus post-surgery. Between eight and 48% of the genes contained one or more homeobox binding sites, depending on the homeobox transcription factor (Table 7). Many individual genes contained binding sites for several of the homeobox transcription factors (Table S6). We also examined genes with homeobox binding sites that were enriched in specific functional categories. The analysis revealed a high number of genes with homeobox binding sites in the PANTHER Molecular Function categories Transcription factor, Extracellular matrix, Signaling molecule, Select regulatory molecule and Receptor (Table 8, Table S7 and S8).
In this study we have shown that profound fat loss substantially alters the global gene expression profile of adipose tissue. To our knowledge, this is the first microarray study on adipose tissue obtained from patients treated with BPD/DS, the most effective bariatric surgical procedure for fat loss . The prospective design and paired analysis reduces potential effects of genetic individuality. Our analysis revealed a marked switch of transcription factor usage, from a high expression of factors such as FOS, JUN, CEBPB, CEBPD, ETS1 and ETS2 in the morbidly obese state, to an increased expression of homeobox genes after profound fat loss. This switch of transcription factor profiles was associated with a strong reduction in inflammation and an increase in collagen expression.
Improved treatment of obesity-related diseases requires a better characterization of molecular mechanisms that determine adipose tissue function. Transcription factors mediate cellular responses to environmental and physiologic stimuli and intracellular signaling cascades. The homeobox genes constitute a special group of highly conserved transcription factors, characterized by a common DNA binding motif, the homeodomain, usually comprising 60 amino acids. Transcription factors encoded by the homeobox genes are master regulators of developmental patterning and segment diversification along the antero-posterior axis in every animal . As targets of endocrine stimuli, they also regulate a host of processes related to normal and abnormal tissue development and function in adults . However, the downstream homeobox target genes are largely unknown . Moreover, little is known about the functions of homeobox genes in relation to obesity, though it has been suggested that they affect processes that are integral to adipose tissue function, such as adipogenesis , . A total of 235 homeobox genes were recently proposed to be functional, categorized into eleven classes subdivided in 102 families .
Global gene expression studies comparing subcutaneous and visceral adipose tissue have revealed depot-specific expression of homeobox genes , . These observations have raised the possibility that differential homeobox gene expression contributes to an increased disease risk associated with visceral obesity. The present study shows that the expression of several homeobox genes is responsive to fat loss, shedding new light on a potentially important role of homeobox transcription in adipose tissue. However, it is unclear at this time whether and to what extent homeobox transcription factors may promote fat loss, improve adipose tissue function, and/or simply respond to changes in adipose tissue mass. In the latter scenario, homeobox transcription factors may mediate the induction of factors involved in adipose tissue remodeling. A role for homeobox genes in tissue regeneration and repair has been observed . Conceivably, the suggested role for homeobox genes in fibroblast differentiation  and adipogenesis ,  may also point towards an involvement in adipose tissue development and adipocyte turnover. Functional studies of homeobox genes in preadipocytes as well as in earlier adipocyte precursor cells may lead to important new knowledge about adipose tissue biology. Our promoter analysis points towards a role for homeobox transcription factors in regulating a host of other transcription factors and a variety of molecular functions. Moreover, our study highlights the possibility that homeobox transcription factors may mitigate inflammation, or that a low expression during obesity may predispose to inflammation because of a reduced capability of tissue repair. Such a protective role for homeobox genes was demonstrated in the lungs of Hoxa5 knock-out mice, which were characterized by a perpetuated infiltration of activated macrophages .
The extracellular matrix is integral to tissue biology, and alteration of its functions has been implicated in a host of chronic conditions including obesity , . In the present study of adipose tissue we found a strong over-expression in morbid obesity of several fibrogenic factors, including macrophage-secreted factors , ETS transcription factors , hypoxia-inducible factor (HIF)-1A , plasminogen activator inhibitor (PAI)-1/SERPINE1 (Entrez Gene 5054) , tenascin C (TNC) (Entrez Gene 3371) , serum/glucocorticoid kinase (SGK) (Entrez Gene 6446) , and connective tissue growth factor (CTGF/CCN2) (Entrez Gene 1490) . It has been proposed that the inflammation and dysfunction of adipose tissue in obesity results from hypoxia-mediated fibrosis, which may promote adipocyte necrosis during continued adipocyte expansion . A possible link between adipose tissue collagen and metabolic health was demonstrated in collagen VI knock-out mice, which were protected from metabolic derangements despite a massive increase in fat mass . In addition, a recent study found a positive correlation between adipose tissue COL6A3 mRNA and body-mass index in humans and suggested that increased collagen VI expression may promote adipose tissue inflammation . Of note, although we found a significant mean increase in COL6A3 expression after fat loss, some patients showed a higher expression in obesity, suggesting some inter-individual variation. Nonetheless, our finding that collagen gene expression was reduced in obesity relative to normal controls and after profound fat loss appears to conflict with these studies, and with the proposed negative impact of adipose tissue collagen on metabolic health. It should be noted that a previous study found no reversal of adipose tissue fibrosis in morbidly obese patients three months after gastric bypass surgery , suggesting that metabolic health can improve despite adipose tissue fibrosis being irreversible. It is conceivable that adipose tissue fibrosis impairs tissue function only during energy surplus, which is unlikely in the post-operative and malabsorptive state. Moreover, excessive collagen deposition may mainly occur in the presence of chronic cellular stress such as hypoxia and inflammatory responses , . Contrary to pro-inflammatory transcription factors, there is evidence from mice suggesting that homeobox transcription factors promote normal tissue repair, despite a substantial increase in collagen synthesis and deposition . Taken together, it is possible that a switch from pro-inflammatory transcription factors to homeobox transcription factors enables a metabolically favorable remodeling of adipose tissue after profound fat loss.
The amount of down-regulated genes and the degree of down-regulation were not as high in a previous study of differential gene expression after Roux-en-Y gastric bypass , which may be explained by the longer duration of our study (one year versus three months). Importantly, we did not employ a very low calorie diet (VLCD) prior to surgery, and post-surgery biopsies were obtained when fat loss began to plateau. We believe that the profound fat loss of our subjects represents close to the largest change in adipose tissue biology that can be obtained with current obesity treatment. The marked down-regulation in the pilot study prompted us to determine whether gene expression may be reduced to a below-normal level after such profound fat loss. By also including lean controls in the main study, we found that most genes returned to a normal level post-surgery despite a higher BMI than in lean controls (33.1 versus 23.0 kg/m2). In this regard, it is interesting to note that insulin sensitivity measures in weight-stable bariatric patients one year post-surgery were closer to lean subjects than to weight-matched controls . Also of note we found a normalization of local markers of insulin resistance in adipose tissue, such as SOCS3 and IRS1 , . Furthermore, the inflammatory gene expression in adipose tissue of obese subjects after a 28-day VLCD was more similar to lean than obese subjects . On the other hand, in the present study many genes involved in metabolic processes did not reach the same expression level post-surgery as in lean controls, suggesting an effect of the higher BMI or other factors in the post-surgery state. This observation encourages further study of the identified metabolic genes in the post-operative state, and weight-matched controls should be used.
We found some variation between the pilot and main study. The higher number of patients in the main study yielded a higher statistical significance, revealing over-representation of PANTHER categories such as Extracellular matrix structural protein, which in the pilot study did not meet the strict significance criterion (p-value<0.01 with Bonferroni correction). In the PANTHER analysis we used fold change ≥1.5 and q-value=0 as inclusion criteria for differentially expressed genes. Several of the same homeobox genes were differentially expressed with fold change ≥1.5 in the pilot study as in the main study, though only IRX3, IRX5, HOXC9 (Entrez Gene 3225), HOXC10 (Entrez Gene 3226) and NANOG (Entrez Gene 79923) also met the q-value=0 criterion (Table S4). Moreover, variation in fold change values in the two analyses may be ascribed to the different detection systems. The chemiluminescence based AB expression array system is more sensitive and covers a broader range of signal intensities than the fluorescence based Illumina system. The combination of two microarray platforms strengthens our database.
Alterations in adipose tissue gene expression after bariatric surgery may partly be due to changes in the cellular composition of the tissue. In this regard, it is noteworthy that we obtained biopsies by surgical excision. Surgical biopsy was demonstrated to be superior to the needle biopsy technique, as the latter incompletely extracts the stromal vascular fraction . Stromal vascular cells include non-adipose cell types such as preadipocytes, macrophages and lymphocytes that are integral to adipose tissue biology and whose functions are altered in obesity , , , , , , . The strong down-regulation of gene expression in our study may partly reflect a reduced number and/or activity of these cell types, in particular immune cells. On the other hand, the genes with an increased expression after surgery are more likely to be actively up-regulated in cells whose phenotypes are altered by profound fat loss, including adipocytes, preadipocytes, and fibroblasts.
In summary, a switch from stress response to increased expression of homeobox transcription factors may represent a mechanism whereby adipose tissue function improves after fat loss and is maintained in normal tissue. Our results encourage further study of the potential roles for homeobox genes in regulating adipose tissue function.
Medication in 16 patients before and one year after bariatric surgery.
(0.01 MB PDF)
Expression values of selected genes previously implicated in obesity, and that were differentially expressed in adipose tissue one year after bariatric surgery.
(0.01 MB PDF)
Description of differentially expressed genes in adipose tissue before versus after bariatric surgery, belonging to the PANTHER functional categories Other transcription factor, Developmental Processes, Homeobox transcription factor and Extracellular Matrix.
(0.02 MB PDF)
Differentially expressed homeobox transcription factors in adipose tissue.
(0.02 MB PDF)
Forward and reverse primers and UPL probes used for qPCR.
(0.01 MB PDF)
Down-regulated genes in adipose tissue after bariatric surgery that contain one or more binding sites for one or more homeobox transcription factors.
(0.03 MB PDF)
Up-regulated genes in adipose tissue after bariatric surgery that contain one or more binding sites for one or more homeobox transcription factors.
(0.01 MB PDF)
PANTHER categories enriched with differentially expressed genes that contain one or more binding sites for one or more homeobox transcription factors.
(0.02 MB PDF)
Materials and methods (AB 1700 Expression Array system)
(0.03 MB DOC)
Functional categorization of differentially expressed genes in adipose tissue after fat loss (AB 1700, n=9). PANTHER was used to search for over-represented functional categories among the most differentially expressed genes (q-value=0, fold change at least 1.5). The color intensity displays the statistical significance (−log p-value) of over- and under-represented PANTHER functional categories. Red color signifies an over-representation of genes mapping to a certain term, blue color an under-representation and white a representation as expected based on the overall distribution on the array. A p-value<0.01 was used as inclusion criterion for categories, with Bonferroni correction for multiple testing. Numbers presented in the table indicate the percentage of genes within a gene set that map to the given category, e.g. 18% of the 499 down-regulated genes map to the biological process ‘Immunity and defense’. The first column states the overall distribution of a term among all human NCBI genes (25,431), e.g. 5% of the genes are expected to map to ‘Immunity and defense’, hence this category is significantly over-represented among the down-regulated genes. Of note, unlike the Illumina data, the data showed an up-regulation of genes involved in B-cell and antibody-mediated immunity (e.g. immunoglobulins). However, the majority of these genes had records that were discontinued in the Entrez Gene database or were listed as hypothetical proteins. Ref, reference (based on all human NCBI genes); Pre, pre-surgery biopsies; Post, post-surgery biopsies; Ctr, lean controls; Arrow up, up-regulated/more expressed genes; Arrow down, down-regulated/less expressed genes (e.g. arrow up in Ctr/Post signifies higher expression in Ctr).
(0.48 MB TIF)
Illumina microarray analysis of adipose tissue before versus one year after bariatric surgery.
(7.15 MB XLS)
AB 1700 microarray analysis of adipose tissue before versus one year after bariatric surgery.
(8.84 MB XLS)
Illumina microarray analysis of adipose tissue comparing morbidly obese and lean subjects.
(7.22 MB XLS)
We thank the staff at Haraldsplass Deaconess Hospital, Bergen, Norway, particularly Drs. Inge Glambek, Barbara Auras Jaatun and Torodd Bruland. We also thank Anne Merete Sellevold, Carol Cook and Solveig Mjelstad Angelskår for expert technical assistance. We are grateful to Drs. Inge Jonassen, Bjarte Dysvik, Kjell Petersen and others at the Norwegian Bioinformatics Platform and the Norwegian Microarray Consortium (national FUGE technology platforms, www.fuge.no).
Competing Interests: The authors have declared that no competing interests exist.
Funding: Funding for this project was provided by the Research Council of Norway (RCN)(http://www.forskningsradet.no), The Norwegian Cancer Society (http://www.kreftforeningen.no), Samarbeidsorganet Helse Vest RHF (http://www.helse-bergen.no/forskning/samarbeidsorganet), Meltzerfondet (http://meltzerfondet.no) and Programstyret for ernaering at the University of Bergen (http://www.uib.no/mofa). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.