Abstract
Adipose tissue plays a role in obesity-related cancers via increased production of inflammatory factors, steroid hormones, and altered adipokines. The impact of weight loss on adipose tissue gene expression may provide insights into pathways linking obesity with cancer risk. We conducted an ancillary study within a randomized trial of diet, exercise, or combined diet + exercise versus control among overweight/obese postmenopausal women. In 45 women, subcutaneous adipose tissue biopsies were conducted at baseline and after 6 months, and changes in adipose tissue gene expression were determined by microarray with an emphasis on prespecified candidate pathways as well as by unsupervised clustering of more than 37,000 transcripts (Illumina). Analyses were conducted first by randomization group and then by degree of weight change at 6-months in all women combined. At 6 months, diet, exercise, and diet + exercise participants lost a mean of 8.8, 2.5, and 7.9 kg (all P < 0.05 vs. no change in controls). There was no significant change in candidate gene expression by intervention group. In analysis by weight change category, greater weight loss was associated a decrease in 17β-hydroxysteroid dehydrogenase-1 (HSD17B1, Ptrend < 0.01) and leptin (LEP, Ptrend < 0.01) expression, and marginally significant increased expression of estrogen receptor-1 (ESR1, Ptrend = 0.08) and insulin-like growth factor–binding protein-3 (IGFBP3, Ptrend = 0.08). Unsupervised clustering revealed 83 transcripts with statistically significant changes. Multiple gene expression changes correlated with changes in associated serum biomarkers. Weight loss was associated with changes in adipose tissue gene expression after 6 months, particularly in two pathways postulated to link obesity and cancer, that is, steroid hormone metabolism and IGF signaling. Cancer Prev Res; 6(3); 217–31. ©2013 AACR.
Introduction
Overweight or obesity and low levels of physical activity are associated with an increased risk of several types of cancer, particularly colorectal and postmenopausal breast cancer (1–4). The metabolic and hormonal consequences of excess fat mass and a sedentary lifestyle are associated with disease promotion and progression. Possible mechanistic pathways, many of which are modulated by change in weight or body composition, appear to involve inflammatory factors, steroid hormones, insulin-like growth factors, or insulin resistance (5–8).
Biomarkers associated with these mechanisms have been investigated in peripheral blood, but there has been little evaluation directly at the level of adipose tissue in humans. Adipose tissue may play a role in several of these proposed mechanisms (9), as it is an endocrine organ integral to regulation of energy homeostasis, insulin sensitivity, and glucose tolerance (10, 11). In addition to adipocytes, white adipose tissue also contains pre-adipocytes, endothelial cells, fibroblasts, and leukocytes (12). In obesity, the number of adipose tissue leukocytes increases dramatically, which, together with a switch in their activation status, results in a state of low-grade chronic inflammation (13). The activated leukocytes, including macrophages and T cells, secrete inflammatory cytokines such as interleukin-1 (IL)-1, IL)-6, and TNFα (14), which likely mediate the inflammation. Inflammation in turn modulates the endocrine functions of the adipocytes, including secretion of adipokines (15).
Systemic inflammation and the changes in circulating adipokine concentrations are now thought to be a major cause of obesity-associated insulin resistance (13) and have also been linked to an increased risk for several types of cancer (16). Systemic inflammation [e.g., elevated serum C-reactive protein (CRP) or serum amyloid A (SAA)] has been associated with multiple cancer types, particularly colorectal cancer (17), whereas low serum concentrations of the anti-inflammatory and insulin-sensitizing adipokine adiponectin have been associated with an increased risk for breast cancer (18). Insulin resistance and the associated hyperinsulinemia, as well as type II diabetes, are associated with breast, colorectal, endometrial, and pancreatic cancer (16). Finally, increased adipose tissue estrogen synthesis plays a role in postmenopausal breast, as well as endometrial cancer (3, 19). In this context, adipose tissue expresses aromatase (cytochrome P450 19A1) and is an important site for extragonadal conversion of androgens to estrogens (20).
Research on the effect of changes in energy balance on human adipose tissue biology is limited and primarily focused on the short-term effects of very low-calorie diets (21–29). The effect of exercise and/or caloric restriction on gene expression in normal mouse mammary glands has been reported (30); together these experimental studies provide intriguing initial evidence for alterations in mRNA expression of inflammatory cytokines and adipokines, particularly IL-6 and leptin (21, 24, 28, 29). The research presented here adds human data on the effect of weight loss on adipose tissue biology. The purpose of this study was to examine, within a randomized controlled trial, the 6-month effects of caloric restriction, exercise, or both on subcutaneous adipose tissue biology in overweight, postmenopausal women. A secondary purpose was to examine the association of degree of weight loss with changes in adipose tissue biology.
Materials and Methods
This study, conducted from 2008 to 2009, was ancillary to a 12-month randomized, controlled trial that tested the effects of 3 lifestyle change interventions (caloric restriction weight loss diet, exercise, or both caloric restriction weight loss diet plus exercise) compared with a usual lifestyle control group on biomarkers of breast cancer risk among 439 women (31–34). This ancillary study examined the effect of the interventions from baseline to 6 months on adipose tissue gene expression, the associations of degree of weight loss with gene expression change, and correlations between adipose tissue gene expression and serum concentrations of specific adipokines, pro-inflammatory cytokines, and sex steroid hormones. All study procedures were reviewed and approved by the Fred Hutchinson Cancer Research Center (FHCRC) Institutional Review Board (Seattle, WA), and all participants provided informed consent.
Subject recruitment
Details of the parent trial have been published previously (31). In brief, participants were female, overweight or obese [body mass index (BMI) ≥25.0 or ≥23.0 kg/m2 if Asian-American], engaging in <100 min/wk of moderate or vigorous activity, and postmenopausal (ages 50–75 years). Specific exclusion criteria included: use of postmenopausal hormones within the past 3 months; history of breast cancer or other serious medical conditions; diagnosed diabetes, fasting blood glucose ≥126 mg/dL, or use of diabetes medications; alcohol intake >2 drinks/d; current smoking; current or planned participation in another structured weight loss program; and use of weight loss medications. Eligible women were randomized into 1 of the 4 study arms, as described previously (31).
During the yearlong period of recruitment for the ancillary study, 178 women attended an information session and 84 were eligible to participate in the parent trial. Of these, 55 expressed interest in participating in the ancillary study, 49 enrolled and consented to the adipose tissue biopsy at baseline and 6 months, and 47 provided both baseline and 6-month adipose tissue samples. Consistent with the primary endpoint analysis of the parent trial (34), 2 participants were removed because of subsequently confirmed high baseline serum estradiol levels (>42 pg/mL); thus, 45 participants were included in the analysis. Recruitment to the ancillary study is presented in Fig. 1.
Lifestyle change interventions
The parent trial interventions have been described in detail previously (31). The calorie-reduced, low-fat diet was modified from the Diabetes Prevention Program (35) and the Look AHEAD trial (36), with a goal of reducing total daily energy intake to 1,200 to 2,000 kcal/d based on baseline weight, daily energy intake from fat to <30%, and a 10% reduction in body weight by 6 months with maintenance thereafter. Participants met individually 2 to 4 times with a study dietitian and attended weekly dietitian-led group meetings in months 1 to 6.
The exercise intervention goal was ≥45 minutes of moderate-to-vigorous intensity exercise, 5 d/wk (≥225 min/wk). Participants attended ≥3 supervised sessions per week at our study facility and for remaining sessions exercised at home. Activities with ≥4 metabolic equivalents (MET; ref. 37) were counted toward the prescribed exercise target.
Women randomized to the diet + exercise intervention received both the dietary weight loss intervention and the aerobic exercise intervention. The diet sessions were run separately from those for the diet group, and participants were instructed not to discuss the diet intervention at the exercise facility. Women randomized to the control group were requested not to change their diet or exercise habits for the duration of the study.
Outcome measures
Baseline and 6-month study measures were obtained by trained study personnel blinded to randomization status. At baseline, anthropometric measures (height, weight, waist circumference) were obtained and body composition was measured using a dual energy x-ray absorptiometry (DEXA) scanner (GE Lunar). All participants wore pedometers (Accusplit) for 7 consecutive days, and maximal aerobic capacity (VO2max) was determined from a graded treadmill test with collection of expired gases. At 6 months, anthropometric measures and DEXA were repeated for the ancillary study, along with the use of pedometers. In addition, intervention women completed diet and/or activity logs, depending on assigned group, to document adherence (31).
Blood samples
At baseline and 6 months, 12-hour fasting blood samples (50 mL) were collected, processed within 1 hour and stored at −70°C. The following analytes were measured: estrone, estradiol, total testosterone, androstenedione, and sex hormone–binding globulin [Reproductive Endocrine Research Laboratory, University of Southern California, Los Angeles, CA; insulin, IL-6, and glucose (Northwest Lipid Research Laboratories, University of Washington, Seattle, WA]; high-sensitivity CRP and SAA (Department of Laboratory Medicine, University of Washington). The interassay coefficients were steroid hormones (8%–13%), sex hormone–binding globulin (SHBG; 5%–7%), insulin (4.5%), glucose (1.4%), CRP (4.7%), SAA (6.2%), and IL-6 (12.4%).
Subcutaneous adipose tissue biopsy and total RNA preparation
Details of the sampling procedure have been previously published (38). Briefly, following an overnight fast, an abdominal subcutaneous adipose tissue sample (≥∼350 mg) was obtained by needle aspiration under local anesthesia. Tissue was flash-frozen and then stored at −80°C. Total RNA was subsequently extracted using the Qiagen RNeasy lipid tissue mini kit (Qiagen), quantified using Ribogreen (Invitrogen) and integrity was assessed on a Bioanalyzer 2100 (Agilent).
BeadChip microarray assays and data preprocessing
One hundred nanograms of total RNA was transcribed using the Illumina TotalPrep RNA Amplification Kit (Ambion). Expression of more than 37,000 mRNA transcripts was assessed using the Human HT-12 v3 Expression Beadchips and the iScan System Reader (Illumina). Each chip was analyzed for assay performance by examination of internal quality control measurements. Data are publicly available in Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo/; Accession number GSE43471). The raw, non-normalized bead summary values output by BeadStudio (Illumina) were used for subsequent analysis. A final gene-level report was generated after background subtraction and quantile normalization. The data quality was assessed with the R package “arrayQualityMetrics” and no outliers were identified. Using Ingenuity Pathway Analysis (IPA; Ingenuity Systems), we examined overlaps with known biologic functions and canonical pathways.
Quantitative real-time PCR
Genes selected for quantitative real-time PCR (qRT-PCR) validation included 5 candidate genes that showed a statistically significant correlation with weight loss (LEP, IGFBP3, HSD17B1, SAA1, and ESR1), plus 4 candidate genes of high interest (ADIPOQ, CYP19A1, TNF, and IL-6) that were not significant in the microarray analysis of weight change. Individuals were included in the qRT-PCR validation if sufficient RNA was available from adipose tissue at both time points. A total of 20 individuals met this criteria, representing each weight change group: gained weight (n = 5), lost ≤ 5% baseline body weight (n = 3), lost 5%–10% baseline body weight (n = 6), and lost >10% baseline body weight (n = 6). Total RNA was reverse transcribed using SuperScript III (Invitrogen), and mRNA expression levels were determined using predeveloped TaqMan gene expression assays on a 7900HT real-time PCR system (Applied Biosystems). Reactions were carried out in triplicate using 2× Gene Expression Master Mix (Applied Biosystems). Levels of GUSB were used for normalization, chosen because of its stable expression in adipose tissue in our own results and those of others (39). The ΔΔCt method was used to relatively quantify mRNA levels across all samples tested.
Statistical analyses
We first investigated changes in gene expression by intervention compared with control. We then examined changes in gene expression by degree of weight change, combining all participants together, including controls. We chose the following weight change categories based on levels of clinically relevant weight loss and approximate quartiles based on our data: gained any weight, lost 0%–5%, lost 5%–10%, and lost >10% (40). We applied the generalized estimating equation (GEE) approach to account for repeated measurements on the same subjects. We conducted trend tests of weight change category on change of gene expression. A 3-tiered approach was used: first, we investigated effects on prespecified candidate genes; second, we investigated effects on candidate pathways and; third, we explored changes across all genes on the microarray. Baseline gene expression values and absolute changes in baseline gene expression were correlated to the serum values by Pearson correlation. Statistical analyses were conducted using SAS (version 9.2, SAS Institute Inc.).
Candidate gene analysis
The R limma Bioconductor package (41) was used to evaluate the significance of differences in change of expression between the control and intervention groups or by degree of weight change, pre- and postintervention. Eighty-two candidate genes were chosen before analysis based on literature reports and involvement with inflammation, immune function, or steroid/hormone metabolism (Supplementary Table S1). These specific genes were evaluated by intervention groups or by a trend test by weight change category for all women combined. P values were adjusted to account for 82 simultaneous comparisons, using the Benjamini and Hochberg method (42).
Gene set analysis
A gene set enrichment analysis (GSEA) approach was used to assess expression change differences between control and intervention or weight change category at the level of biochemical pathways (43, 44). Curated gene sets related to inflammation, immune function, and steroid/hormone metabolism were chosen from the Molecular Signatures Database (43). Fifteen gene sets were analyzed, all of which contained 10 or more of those genes (Supplementary Table S2). The gene sets together included 888 genes. As in the candidate gene analysis, the gene sets were evaluated by testing whether the change in expression of each of the intervention group samples was identical to the control. In addition, we assessed the significance of change in expression, postintervention minus preintervention, relative to weight change category.
Discovery gene analysis
To identify novel genes whose expression in adipose tissue was impacted by weight loss, all genes with measureable mRNA expression in our samples were evaluated for differences by weight change using the same statistical tools as for the candidate gene analysis. P values were adjusted for all 37,834 genes measured.
Results
Anthropometrics and main intervention effects
The 45 participants were randomly assigned to the diet intervention (n = 8), the exercise intervention (n = 14), the diet + exercise intervention (n = 16), or the control group (n = 7; Fig. 1). Unequal group sizes are attributable to this being a subset of consenting participants within a larger trial. Of these participants, 39 underwent DEXA scans at baseline and 6-months [diet (n = 7), exercise (n = 12), diet + exercise (n = 13), control (n = 7)] and 35 provided pedometer step counts at baseline and 6 months [diet (n = 6), exercise (n = 11), diet + exercise (n = 12), control (n = 6)].
The anthropometric characteristics of the subjects randomized to the 4 intervention groups were comparable at baseline (Table 1). The mean age of participants was 58.5 ± 4.6 years. Participants had a mean BMI of 31.3 ± 4.3 kg/m2, 47.6% ± 4.2% body fat, and VO2 max of 24.4 ± 4.3 mL/kg/min (consistent with poor aerobic fitness). Baseline values for serum levels of selected biomarkers of cancer risk are also presented in Table 1.
Baseline characteristics and fasting serum values of participants in the Adipose Biology Study (n = 45)
. | . | By intervention group . | |||
---|---|---|---|---|---|
. | All . | Control . | Diet . | Exercise . | Diet + exercise . |
. | Mean (SD) . | Mean (SD) . | Mean (SD) . | Mean (SD) . | Mean (SD) . |
. | n = 45 . | n = 7 . | n = 8 . | n = 14 . | n = 16 . |
Age, y | 58.5 (4.6) | 55.5 (2.9) | 54.5 (3.2) | 60.9 (5.3) | 59.7 (3.1) |
Race/ethnicity, n (%) | |||||
Non-Hispanic White | 43 (96) | 6 (86) | 8 (100) | 14 (100) | 15 (94) |
Othera | 2 (4) | 1 (14) | 0 (0) | 0 (0) | 1 (6) |
Weight, kg | 84.6 (12.6) | 92.8 (10.3) | 77.1 (10.0) | 85.3 (14.8) | 84.2 (11.2) |
BMI, kg/m2 | 31.3 (4.3) | 33.2 (3.8) | 28.9 (3.0) | 30.8 (3.5) | 32.1 (5.2) |
Body fat, % | 47.6 (4.2) | 49.0 (4.0) | 46.0 (3.4) | 47.1 (4.7) | 48.4 (4.2) |
Waist, cm | 96.4 (10.3) | 102.4 (4.4) | 91.3 (11.0) | 97.8 (12.8) | 95.0 (8.4) |
VO2max, mL/kg/min | 24.4 (4.3) | 23.3 (4.4) | 25.5 (5.7) | 23.2 (4.0) | 25.6 (3.7) |
Steps/wk | 5,450 (2,253) | 5,161 (2,565) | 4,955 (2,229) | 5,828 (1,815) | 5,518 (2,581) |
Fasting serum concentrations | |||||
Insulin, μIU/mL | 12.1 (6.3) | 15.8 (11.4) | 9.9 (5.4) | 13.1 (4.6) | 10.7 (4.7) |
Glucose, mg/dL | 95.8 (8.2) | 92.7 (5.4) | 92.8 (4.1) | 97.7 (10.8) | 97.0 (7.9) |
HOMA-IR | 2.9 (1.5) | 3.6 (2.4) | 2.3 (1.3) | 3.2 (1.2) | 2.6 (1.3) |
CRP, mg/L | 3.5 (3.0) | 2.9 (1.4) | 2.6 (1.9) | 4.5 (4.2) | 3.5 (2.6) |
SAA, mg/L | 6.8 (4.4) | 7.4 (3.6) | 5.6 (4.4) | 5.9 (2.8) | 8.0 (5.7) |
Estradiol, pg/mL | 12.3 (6.1) | 13.5 (8.4) | 9.7 (5.3) | 13.3 (6.0) | 12.1 (5.7) |
Estrone, pg/mL | 39.7 (17.0) | 40.9 (23.8) | 36.2 (20.3) | 43.3 (16.0) | 37.7 (13.5) |
Total testosterone, ng/dL | 32.2 (17.8) | 31.1 (24.1) | 29.8 (15.8) | 32.5 (15.3) | 33.5 (19.4) |
Androstenedione, ng/dL | 57.3 (24.0) | 51.9 (24.9) | 57.3 (28.6) | 59.3 (27.4) | 58.0 (19.8) |
Free estradiol, pg/mL | 0.33 (0.16) | 0.37 (0.21) | 0.26 (0.15) | 0.34 (0.16) | 0.33 (0.16) |
Free testosterone, pg/mL | 6.8 (3.6) | 6.9 (4.7) | 6.3 (3.1) | 6.4 (2.7) | 7.4 (4.4) |
SHBG, nmol/L | 38.5 (13.9) | 33.4 (12.0) | 39.2 (17.2) | 44.2 (16.4) | 35.5 (9.5) |
Leptin, ng/mL | 28.9 (11.3) | 32.6 (9.4) | 24.3 (11.0) | 27.4 (12.2) | 30.9 (11.4) |
Adiponectin, μg/mL | 13.8 (5.9) | 13.5 (5.8) | 12.2 (4.9) | 13.9 (6.3) | 14.7 (6.5) |
. | . | By intervention group . | |||
---|---|---|---|---|---|
. | All . | Control . | Diet . | Exercise . | Diet + exercise . |
. | Mean (SD) . | Mean (SD) . | Mean (SD) . | Mean (SD) . | Mean (SD) . |
. | n = 45 . | n = 7 . | n = 8 . | n = 14 . | n = 16 . |
Age, y | 58.5 (4.6) | 55.5 (2.9) | 54.5 (3.2) | 60.9 (5.3) | 59.7 (3.1) |
Race/ethnicity, n (%) | |||||
Non-Hispanic White | 43 (96) | 6 (86) | 8 (100) | 14 (100) | 15 (94) |
Othera | 2 (4) | 1 (14) | 0 (0) | 0 (0) | 1 (6) |
Weight, kg | 84.6 (12.6) | 92.8 (10.3) | 77.1 (10.0) | 85.3 (14.8) | 84.2 (11.2) |
BMI, kg/m2 | 31.3 (4.3) | 33.2 (3.8) | 28.9 (3.0) | 30.8 (3.5) | 32.1 (5.2) |
Body fat, % | 47.6 (4.2) | 49.0 (4.0) | 46.0 (3.4) | 47.1 (4.7) | 48.4 (4.2) |
Waist, cm | 96.4 (10.3) | 102.4 (4.4) | 91.3 (11.0) | 97.8 (12.8) | 95.0 (8.4) |
VO2max, mL/kg/min | 24.4 (4.3) | 23.3 (4.4) | 25.5 (5.7) | 23.2 (4.0) | 25.6 (3.7) |
Steps/wk | 5,450 (2,253) | 5,161 (2,565) | 4,955 (2,229) | 5,828 (1,815) | 5,518 (2,581) |
Fasting serum concentrations | |||||
Insulin, μIU/mL | 12.1 (6.3) | 15.8 (11.4) | 9.9 (5.4) | 13.1 (4.6) | 10.7 (4.7) |
Glucose, mg/dL | 95.8 (8.2) | 92.7 (5.4) | 92.8 (4.1) | 97.7 (10.8) | 97.0 (7.9) |
HOMA-IR | 2.9 (1.5) | 3.6 (2.4) | 2.3 (1.3) | 3.2 (1.2) | 2.6 (1.3) |
CRP, mg/L | 3.5 (3.0) | 2.9 (1.4) | 2.6 (1.9) | 4.5 (4.2) | 3.5 (2.6) |
SAA, mg/L | 6.8 (4.4) | 7.4 (3.6) | 5.6 (4.4) | 5.9 (2.8) | 8.0 (5.7) |
Estradiol, pg/mL | 12.3 (6.1) | 13.5 (8.4) | 9.7 (5.3) | 13.3 (6.0) | 12.1 (5.7) |
Estrone, pg/mL | 39.7 (17.0) | 40.9 (23.8) | 36.2 (20.3) | 43.3 (16.0) | 37.7 (13.5) |
Total testosterone, ng/dL | 32.2 (17.8) | 31.1 (24.1) | 29.8 (15.8) | 32.5 (15.3) | 33.5 (19.4) |
Androstenedione, ng/dL | 57.3 (24.0) | 51.9 (24.9) | 57.3 (28.6) | 59.3 (27.4) | 58.0 (19.8) |
Free estradiol, pg/mL | 0.33 (0.16) | 0.37 (0.21) | 0.26 (0.15) | 0.34 (0.16) | 0.33 (0.16) |
Free testosterone, pg/mL | 6.8 (3.6) | 6.9 (4.7) | 6.3 (3.1) | 6.4 (2.7) | 7.4 (4.4) |
SHBG, nmol/L | 38.5 (13.9) | 33.4 (12.0) | 39.2 (17.2) | 44.2 (16.4) | 35.5 (9.5) |
Leptin, ng/mL | 28.9 (11.3) | 32.6 (9.4) | 24.3 (11.0) | 27.4 (12.2) | 30.9 (11.4) |
Adiponectin, μg/mL | 13.8 (5.9) | 13.5 (5.8) | 12.2 (4.9) | 13.9 (6.3) | 14.7 (6.5) |
aHispanic white (n = 1) and “other” (n = 1).
Results for all participants from the parent trial have been previously reported (31–34). Participants in this ancillary study experienced mean weight changes as follows (P values in reference to control group): diet, −11.3% (P < 0.001); exercise, −3.0% (P = 0.03); and diet + exercise, −9.4% (P < 0.0001); controls, +1.0%. Compared with no change in controls, the mean percentage body fat loss was: diet, −12.6% (P = 0.02); exercise, −3.1% (P = 0.56); and diet + exercise, −13.2% (P < 0.01). Compared with the control group (−333 steps/wk), the diet group had no statistically significant change in number of steps per week (+1,071 steps/wk, P = 0.16), whereas both the exercise (+2,570 steps/wk, P < 0.0001) and diet + exercise (+4,182 steps/wk, P < 0.0001) had an increase in steps per week. Baseline and 6-month serum marker values are presented in Supplementary Table S3.
Changes in adipose tissue candidate genes and candidate gene sets: analyses by intervention arm
In the first candidate gene analysis, microarray data from 45 participants were used to compare changes in expression between baseline and 6 months, for 82 candidate genes, between intervention groups and the control group (Supplementary Fig. S1). The unadjusted analysis indicated that 9 genes had significantly altered expression in the diet intervention compared with control. Four genes had decreased expression: hydroxysteroid (17-beta) dehydrogenase 1 (HSD17B1; P = 0.01, adj. P = 0.29); leptin (LEP; P = 0.03/0.29), insulin receptor (INSR; P = 0.02/0.29), and telomeric repeat–binding factor 1 (TERF1; P = 0.03/0.29). Five genes had increased expression: telomeric repeat–binding factor 2 (TERF2; P = 0.03/0.29), telomeric repeat–binding factor 2, interacting protein (TERF2IP; P = 0.02/0.29); COX-2 (also known as prostaglandin-endoperoxide synthase 2, PTGS2; P = 0.004/0.29); IL1B (P = 0.03/0.29), and estrogen receptor (ESR1; P = 0.03/0.29). No significant change in gene expression was observed in the exercise-only group. In the diet + exercise intervention, expression of 3 genes was significantly decreased (HSD17B1: P = 0.02/0.45; LEP: P = 0.04/0.60; fatty acid–binding protein 1 FABP1: P = 0.02/0.45), and expression of 2 genes was significantly increased (ESR1: P = 0.02/0.45; TERF2IP: P = 0.01/P = 0.45).
For the initial gene set analysis, the microarray data were also used to compare changes in expression between baseline and 6-month visits for 15 predefined gene sets by intervention group (Supplementary Table S2). Differential expression was observed in 1 gene set in the diet intervention (jak-STAT signaling pathway: P = 0.004/0.06), and 2 gene sets in the diet + exercise intervention (IL-6 pathway: P = 0.0002/P = 0.003; and toll pathway: P = 0.03/0.20). However, when adjusted for multiple comparisons, only the IL-6 pathway in the diet + exercise intervention remained statistically significant.
Changes in adipose tissue candidate genes and candidate gene sets: analyses by weight change category
In the second candidate gene analysis, microarray data were used to evaluate the association of percentage change by weight change category between the baseline and 6-month visits on changes in gene expression among all study participants for significantly regulated genes (Fig. 2A) and all 82 candidate genes (Fig. 2B). In this analysis, the sample was split into approximate quartiles, gained weight (n = 10), lost ≤ 5% baseline body weight (n = 10), lost 5%–10% baseline body weight (n = 11), and lost >10% baseline body weight (n = 14). In sex hormone–related genes, a decrease in expression of HSD17B1 (P = 0.0002/0.008), STS (P = 0.02/0.20), and HSD17B10 (P = 0.045/0.37) was noted with greater weight loss, but only HSD17B1 remained statistically significant after adjustment for multiple comparisons. A trend for an increase in ESR1 (P = 0.004/0.08) with greater weight loss was noted.
Differential gene expression analyzed by weight change. A, heat map displaying significant gene expression changes over 6 months in individuals with no weight loss compared to those with <5%, 5>–10%, and >10% weight loss. Gradients indicate the level of gene expression change over time. B, heat map displaying the gene expression change of 82 candidate genes over 6 months in individuals with no weight loss compared to those with <5%, 5>–10%, and >10% weight loss. Gradients indicate the level of gene expression change over time.
Differential gene expression analyzed by weight change. A, heat map displaying significant gene expression changes over 6 months in individuals with no weight loss compared to those with <5%, 5>–10%, and >10% weight loss. Gradients indicate the level of gene expression change over time. B, heat map displaying the gene expression change of 82 candidate genes over 6 months in individuals with no weight loss compared to those with <5%, 5>–10%, and >10% weight loss. Gradients indicate the level of gene expression change over time.
In inflammation-related genes, the expression of ICAM4 (P = 0.01/0.13) and SAA1 (P = 0.01/0.13) decreased, whereas that of VCAM1 (P = 0.04/0.37) and CRP (P = 0.01/0.13) increased with weight loss in the unadjusted analysis; however, none remained statistically significant after adjustment for multiple comparisons.
Expression of LEP decreased with greater weight loss (P = 0.0003/0.002) and expression of IGFBP3 increased with greater weight loss (P = 0.003/0.08).
In the gene-set analysis, percent weight change was associated with 2 KEGG pathways: C21 steroid hormone metabolism (P = 0.005, adj. P = 0.07) and jak-STAT signaling (P = 0.042, adj. P = 0.30), but neither association was significant when adjusted for multiple comparisons (Supplementary Table S2).
To validate our findings, we carried out qRT-PCR gene expression analysis on a subset of 20 participants for 9 candidate genes. We observed statistically significant decreases in expression of LEP, SAA1, and HSD17B1 and increased expression of IGFBP3 and ESR with weight loss, which is consistent with our findings in the microarray analysis (Supplementary Fig. S2). No statistically significant changes were detected in ADIPOQ, CYP19A1, TNF, or IL-6, also consistent with our microarray results.
Changes in all analyzed adipose tissue genes: analyses by weight change
Unsupervised clustering of more than 37,000 transcripts by weight change category revealed 83 transcripts with statistically significant P values for trend (P < 0.05 after adjustment for multiple comparisons). This included 32 genes with decreased expression and 51 genes whose expression increased with greater weight loss (Supplementary Table S4). Of these, 31 genes had a ≥25% change in expression in the largest weight loss category (Table 2). Of note, LEP was the only one of our prespecified candidate genes that also was significant in the unsupervised analysis.
Change in adipose tissue gene expression (unsupervised) by weight change category (all women combined)
. | . | Weight change group . | ||||
---|---|---|---|---|---|---|
. | All . | Gained . | ≤5% loss . | 5%–10% loss . | >10% loss . | . |
. | n = 45 . | n = 10 . | n = 10 . | n = 11 . | n = 14 . | . |
Gene . | % change . | % change . | % change . | % change . | % change . | Adj. Ptrenda . |
BCYRN1 | −19% | 78% | 19% | −34% | −59% | 0.009 |
HSZFP36 | 2% | −17% | −9% | 0% | 30% | 0.009 |
ABCC6 | −13% | 4% | 1% | −4% | −35% | 0.013 |
PELI2 | −2% | −25% | −7% | −6% | 26% | 0.013 |
ZNF33B | 4% | −15% | −6% | 7% | 27% | 0.013 |
HIST1H4H | −22% | 8% | 9% | −16% | −53% | 0.014 |
CEBPZ | 4% | −16% | −2% | 2% | 27% | 0.014 |
C6 | 49% | −15% | 6% | 12% | 253% | 0.014 |
MRAP | −17% | 9% | −8% | −4% | −44% | 0.017 |
TMEM189 | −14% | 3% | −2% | −8% | −34% | 0.017 |
MGC13057 | −11% | 29% | −3% | −13% | −34% | 0.017 |
MAST4 | −3% | 41% | 6% | −4% | −30% | 0.017 |
LOC643496 | 0% | −24% | −11% | 2% | 32% | 0.017 |
HS.333084 | −30% | 42% | 15% | −21% | −74% | 0.024 |
HIST1H1C | −5% | 31% | 23% | −4% | −38% | 0.024 |
GPLD1 | −16% | 2% | 2% | −13% | −39% | 0.030 |
LEP | −19% | 0% | 1% | −13% | −44% | 0.031 |
BOAT | −4% | 32% | 2% | −3% | −26% | 0.031 |
HIST2H2AA3 | 0% | 47% | 11% | −6% | −26% | 0.031 |
LOC402560 | −3% | 19% | 16% | −2% | −26% | 0.031 |
FLJ10213 | 4% | −15% | −3% | 6% | 25% | 0.035 |
ATP6V0E2 | −7% | 14% | −4% | −1% | −26% | 0.035 |
ZNF818 | 2% | −23% | −7% | 8% | 26% | 0.035 |
ACVR1 | 2% | −14% | −14% | 1% | 30% | 0.036 |
CSTA | −5% | 34% | 5% | −7% | −30% | 0.038 |
LOC390349 | −7% | 4% | 28% | −12% | −29% | 0.038 |
PCDH18 | 1% | −29% | −10% | −1% | 43% | 0.038 |
C20ORF24 | −11% | 9% | −2% | −11% | −28% | 0.048 |
ZNF823 | 5% | −10% | −2% | 1% | 25% | 0.048 |
C9ORF116 | −3% | 26% | 3% | 5% | −27% | 0.049 |
LOC401115 | −4% | 51% | −1% | −13% | −28% | 0.049 |
. | . | Weight change group . | ||||
---|---|---|---|---|---|---|
. | All . | Gained . | ≤5% loss . | 5%–10% loss . | >10% loss . | . |
. | n = 45 . | n = 10 . | n = 10 . | n = 11 . | n = 14 . | . |
Gene . | % change . | % change . | % change . | % change . | % change . | Adj. Ptrenda . |
BCYRN1 | −19% | 78% | 19% | −34% | −59% | 0.009 |
HSZFP36 | 2% | −17% | −9% | 0% | 30% | 0.009 |
ABCC6 | −13% | 4% | 1% | −4% | −35% | 0.013 |
PELI2 | −2% | −25% | −7% | −6% | 26% | 0.013 |
ZNF33B | 4% | −15% | −6% | 7% | 27% | 0.013 |
HIST1H4H | −22% | 8% | 9% | −16% | −53% | 0.014 |
CEBPZ | 4% | −16% | −2% | 2% | 27% | 0.014 |
C6 | 49% | −15% | 6% | 12% | 253% | 0.014 |
MRAP | −17% | 9% | −8% | −4% | −44% | 0.017 |
TMEM189 | −14% | 3% | −2% | −8% | −34% | 0.017 |
MGC13057 | −11% | 29% | −3% | −13% | −34% | 0.017 |
MAST4 | −3% | 41% | 6% | −4% | −30% | 0.017 |
LOC643496 | 0% | −24% | −11% | 2% | 32% | 0.017 |
HS.333084 | −30% | 42% | 15% | −21% | −74% | 0.024 |
HIST1H1C | −5% | 31% | 23% | −4% | −38% | 0.024 |
GPLD1 | −16% | 2% | 2% | −13% | −39% | 0.030 |
LEP | −19% | 0% | 1% | −13% | −44% | 0.031 |
BOAT | −4% | 32% | 2% | −3% | −26% | 0.031 |
HIST2H2AA3 | 0% | 47% | 11% | −6% | −26% | 0.031 |
LOC402560 | −3% | 19% | 16% | −2% | −26% | 0.031 |
FLJ10213 | 4% | −15% | −3% | 6% | 25% | 0.035 |
ATP6V0E2 | −7% | 14% | −4% | −1% | −26% | 0.035 |
ZNF818 | 2% | −23% | −7% | 8% | 26% | 0.035 |
ACVR1 | 2% | −14% | −14% | 1% | 30% | 0.036 |
CSTA | −5% | 34% | 5% | −7% | −30% | 0.038 |
LOC390349 | −7% | 4% | 28% | −12% | −29% | 0.038 |
PCDH18 | 1% | −29% | −10% | −1% | 43% | 0.038 |
C20ORF24 | −11% | 9% | −2% | −11% | −28% | 0.048 |
ZNF823 | 5% | −10% | −2% | 1% | 25% | 0.048 |
C9ORF116 | −3% | 26% | 3% | 5% | −27% | 0.049 |
LOC401115 | −4% | 51% | −1% | −13% | −28% | 0.049 |
NOTE: Only transcripts with adjusted Ptrend < 0.05 and ± ≥25% change in expression in the >10% weight loss category are shown.
aAdjusted to account for 82 simultaneous comparisons, using the Benjamini and Hochberg method (41). A full table of all 83 genes is cited in Supplementary Table S4.
We computed the overlap between the 83 statistically significant genes and canonical pathways and GO (gene ontology) categories used GSEA and the Molecular Signatures Database (MSigDB; 43; Table 3). The results show overlap with signaling pathways, such as mTOR signaling (P = 0.011), and a strong overlap with gene sets related to mRNA metabolism and translation. Using IPA, we had similar results. In that analysis, mTOR signaling was the top canonical pathway (P = 0.002; data not shown).
Unsupervised clustering analysis canonical pathway and GO gene set hits from the Molecular Signatures Database
Gene set namea . | Database . | Description . | No. of genes in gene set (K) . | No. of genes in overlap (k) . | k/K . | P . |
---|---|---|---|---|---|---|
Signaling pathways | ||||||
mTOR pathway | Biocarta | mTOR signaling pathway | 23 | 2 | 0.087 | 0.011 |
ABC transporters | KEGG | ABC transporters | 44 | 2 | 0.0455 | 0.037 |
mRNA metabolism | ||||||
Deadenylation of mRNA | Reactome | Genes involved in deadenylation of mRNA | 22 | 2 | 0.0909 | 0.010 |
Metabolism of mRNA | Reactome | Genes involved in metabolism of mRNA | 46 | 3 | 0.0652 | 0.004 |
Metabolism of RNA | Reactome | Genes involved in metabolism of RNA | 96 | 4 | 0.0417 | 0.004 |
Transport of mature mRNA derived from an intron-containing transcript | Reactome | Genes involved in transport of mature mRNA derived from an intron-containing transcript | 51 | 2 | 0.0392 | 0.049 |
Translation | ||||||
Regulation of translational initiation | Gene Ontology | Genes annotated by the GO term GO:0006446. Any process that modulates the frequency, rate or extent of translational initiation. | 31 | 2 | 0.0645 | 0.019 |
Translation factor activity nucleic acid binding | Gene Ontology | Genes annotated by the GO term GO:0008135. Functions during translation by binding nucleic acids during polypeptide synthesis at the ribosome. | 39 | 2 | 0.0513 | 0.030 |
Translation initiation factor activity | Gene Ontology | Genes annotated by the GO term GO:0003743. Functions in the initiation of ribosome-mediated translation of mRNA into a polypeptide. | 24 | 2 | 0.0833 | 0.012 |
Translation regulator activity | Gene Ontology | Genes annotated by the GO term GO:0045182. Any substance involved in the initiation, activation, perpetuation, repression, or termination of polypeptide synthesis at the ribosome. | 41 | 2 | 0.0488 | 0.033 |
Translational initiation | Gene Ontology | Genes annotated by the GO term GO:0006413. The process preceding formation of the peptide bond between the first two amino acids of a protein. This includes the formation of a complex of the ribosome, mRNA, and an initiation complex that contains the first amino acyl-tRNA. | 39 | 2 | 0.0513 | 0.030 |
Cell life span | ||||||
Apoptotic execution phase | Reactome | Genes involved in apoptotic execution phase | 48 | 2 | 0.0417 | 0.044 |
Packaging of telomere ends | Reactome | Genes involved in packaging of telomere ends | 49 | 2 | 0.0408 | 0.045 |
Others | ||||||
ATP binding | Gene Ontology | Genes annotated by the GO term GO:0005524. Interacting selectively with ATP, adenosine 5′-triphosphate, a universally important coenzyme and enzyme regulator. | 151 | 4 | 0.0265 | 0.021 |
Gene set namea . | Database . | Description . | No. of genes in gene set (K) . | No. of genes in overlap (k) . | k/K . | P . |
---|---|---|---|---|---|---|
Signaling pathways | ||||||
mTOR pathway | Biocarta | mTOR signaling pathway | 23 | 2 | 0.087 | 0.011 |
ABC transporters | KEGG | ABC transporters | 44 | 2 | 0.0455 | 0.037 |
mRNA metabolism | ||||||
Deadenylation of mRNA | Reactome | Genes involved in deadenylation of mRNA | 22 | 2 | 0.0909 | 0.010 |
Metabolism of mRNA | Reactome | Genes involved in metabolism of mRNA | 46 | 3 | 0.0652 | 0.004 |
Metabolism of RNA | Reactome | Genes involved in metabolism of RNA | 96 | 4 | 0.0417 | 0.004 |
Transport of mature mRNA derived from an intron-containing transcript | Reactome | Genes involved in transport of mature mRNA derived from an intron-containing transcript | 51 | 2 | 0.0392 | 0.049 |
Translation | ||||||
Regulation of translational initiation | Gene Ontology | Genes annotated by the GO term GO:0006446. Any process that modulates the frequency, rate or extent of translational initiation. | 31 | 2 | 0.0645 | 0.019 |
Translation factor activity nucleic acid binding | Gene Ontology | Genes annotated by the GO term GO:0008135. Functions during translation by binding nucleic acids during polypeptide synthesis at the ribosome. | 39 | 2 | 0.0513 | 0.030 |
Translation initiation factor activity | Gene Ontology | Genes annotated by the GO term GO:0003743. Functions in the initiation of ribosome-mediated translation of mRNA into a polypeptide. | 24 | 2 | 0.0833 | 0.012 |
Translation regulator activity | Gene Ontology | Genes annotated by the GO term GO:0045182. Any substance involved in the initiation, activation, perpetuation, repression, or termination of polypeptide synthesis at the ribosome. | 41 | 2 | 0.0488 | 0.033 |
Translational initiation | Gene Ontology | Genes annotated by the GO term GO:0006413. The process preceding formation of the peptide bond between the first two amino acids of a protein. This includes the formation of a complex of the ribosome, mRNA, and an initiation complex that contains the first amino acyl-tRNA. | 39 | 2 | 0.0513 | 0.030 |
Cell life span | ||||||
Apoptotic execution phase | Reactome | Genes involved in apoptotic execution phase | 48 | 2 | 0.0417 | 0.044 |
Packaging of telomere ends | Reactome | Genes involved in packaging of telomere ends | 49 | 2 | 0.0408 | 0.045 |
Others | ||||||
ATP binding | Gene Ontology | Genes annotated by the GO term GO:0005524. Interacting selectively with ATP, adenosine 5′-triphosphate, a universally important coenzyme and enzyme regulator. | 151 | 4 | 0.0265 | 0.021 |
aOnly gene sets with P < 0.05 are shown.
Correlation of gene expression with serum blood markers
The baseline and end-of-study serum markers are presented in Supplementary Table S3. At baseline, there were statistically significant correlations between serum levels of sex hormones, inflammatory, and metabolic markers and adipokines, and gene expression for some sex steroid–related, inflammation-related, adipose tissue lipid metabolism genes, and other genes of interest (Table 4). For example, serum levels of estradiol were positively associated with the expression of HSD17B1 and HSD11B1 and negatively associated with HSD11B1L and HSD17B8. For adipokine/inflammation-related genes, serum levels of adiponectin were positively associated with the expression of ADIPOQ and RBP4 and negatively with the expression of LEP. Some associations were also observed between CRP levels and expression of SAA1 and LEP, but not CRP. Similarly, there were some associations about genes regulating lipid metabolism with serum glucose and insulin. As expected, a strong correlation was observed between serum leptin and the expression of LEP.
Associations (Pearson correlation coefficients) between baseline serum biomarker levels and baseline gene expression
Pathwaya . | Biomarkers . | |||
---|---|---|---|---|
Sex steroid–related genes | Estradiol, pg/mL | Testosterone, ng/dL | SHBG, nmol/L | Free estradiol, pg/mL |
ESR1 | −0.21b | −0.06 | 0.54 | −0.32 |
0.17c | 0.72 | 0.0001 | 0.03 | |
HSD17B11 | 0.34 | 0.04 | −0.09 | 0.34 |
0.02 | 0.79 | 0.54 | 0.02 | |
HSD11B1 | 0.33 | 0.08 | −0.17 | 0.37 |
0.03 | 0.59 | 0.27 | 0.013 | |
HSD11B1L | −0.38 | −0.18 | −0.09 | −0.34 |
0.01 | 0.24 | 0.56 | 0.02 | |
HSD17B8 | −0.31 | −0.08 | 0.13 | −0.29 |
0.04 | 0.58 | 0.39 | 0.049 | |
Inflammation-related genes | Leptin, ng/mL | Insulin, μU/mL | CRP, mg/L | Adiponectin, μg/mL |
ADIPOQ | 0.13 | 0.17 | −0.28 | 0.44 |
0.38 | 0.27 | 0.06 | 0.002 | |
IFNG | 0.32 | −0.09 | 0.15 | −0.09 |
0.034 | 0.54 | 0.33 | 0.56 | |
IL10 | −0.26 | −0.3 | 0.15 | −0.03 |
0.09 | 0.047 | 0.32 | 0.85 | |
SAA1 | 0.2 | −0.15 | 0.44 | 0.071 |
0.18 | 0.34 | 0.0023 | 0.65 | |
TLR4 | 0.27 | 0.39 | −0.01 | −0.13 |
0.08 | 0.009 | 0.97 | 0.41 | |
Adipose tissue lipid metabolism–related genes | Insulin, μU/mL | Glucose, mg/dL | Adiponectin, μg/mL | |
LIPE | −0.44 | −0.3 | 0.42 | NA |
0.003 | 0.047 | 0.004 | ||
ABCA1 | −0.04 | 0.33 | 0.01 | NA |
0.81 | 0.025 | 0.95 | ||
Other genes | Leptin, ng/mL | CRP, mg/L | Adiponectin, μg/mL | |
LEP | 0.46 | 0.33 | −0.37 | NA |
0.0016 | 0.026 | 0.013 | ||
RBP4 | −0.1 | −0.01 | 0.43 | NA |
0.51 | 0.97 | 0.0034 |
Pathwaya . | Biomarkers . | |||
---|---|---|---|---|
Sex steroid–related genes | Estradiol, pg/mL | Testosterone, ng/dL | SHBG, nmol/L | Free estradiol, pg/mL |
ESR1 | −0.21b | −0.06 | 0.54 | −0.32 |
0.17c | 0.72 | 0.0001 | 0.03 | |
HSD17B11 | 0.34 | 0.04 | −0.09 | 0.34 |
0.02 | 0.79 | 0.54 | 0.02 | |
HSD11B1 | 0.33 | 0.08 | −0.17 | 0.37 |
0.03 | 0.59 | 0.27 | 0.013 | |
HSD11B1L | −0.38 | −0.18 | −0.09 | −0.34 |
0.01 | 0.24 | 0.56 | 0.02 | |
HSD17B8 | −0.31 | −0.08 | 0.13 | −0.29 |
0.04 | 0.58 | 0.39 | 0.049 | |
Inflammation-related genes | Leptin, ng/mL | Insulin, μU/mL | CRP, mg/L | Adiponectin, μg/mL |
ADIPOQ | 0.13 | 0.17 | −0.28 | 0.44 |
0.38 | 0.27 | 0.06 | 0.002 | |
IFNG | 0.32 | −0.09 | 0.15 | −0.09 |
0.034 | 0.54 | 0.33 | 0.56 | |
IL10 | −0.26 | −0.3 | 0.15 | −0.03 |
0.09 | 0.047 | 0.32 | 0.85 | |
SAA1 | 0.2 | −0.15 | 0.44 | 0.071 |
0.18 | 0.34 | 0.0023 | 0.65 | |
TLR4 | 0.27 | 0.39 | −0.01 | −0.13 |
0.08 | 0.009 | 0.97 | 0.41 | |
Adipose tissue lipid metabolism–related genes | Insulin, μU/mL | Glucose, mg/dL | Adiponectin, μg/mL | |
LIPE | −0.44 | −0.3 | 0.42 | NA |
0.003 | 0.047 | 0.004 | ||
ABCA1 | −0.04 | 0.33 | 0.01 | NA |
0.81 | 0.025 | 0.95 | ||
Other genes | Leptin, ng/mL | CRP, mg/L | Adiponectin, μg/mL | |
LEP | 0.46 | 0.33 | −0.37 | NA |
0.0016 | 0.026 | 0.013 | ||
RBP4 | −0.1 | −0.01 | 0.43 | NA |
0.51 | 0.97 | 0.0034 |
aOnly genes with at least one statistically significant association (P < 0.05) are shown.
bCorrelation coefficient.
cP value.
The correlation between changes in serum levels and changes in gene expression is presented in Table 5. Again, here we observed some expected associations, for example, between changes in ESR1, SHBG, and HSD17B7 and corresponding changes in serum estradiol with fairly strong correlation coefficients (r = 0.33–0.47). A reduction in serum glucose was associated with a decrease in the expression of ADIPOQ and increases in IL1B, IL6, and IL8. For genes involved in adipose tissue lipid metabolism, decreases in serum leptin and insulin were associated with an increase in expression of ABCA1, As expected, changes in LEP were positively associated with serum leptin and glucose but inversely with adiponectin.
Associations (Pearson correlation coefficients) between changes in serum biomarker levels and changes in gene expression from baseline to 6 months
Pathwaya . | Biomarkers . | ||||
---|---|---|---|---|---|
Sex steroid–related genes | Estrone, pg/mL | Estradiol, pg/mL | Testosterone, ng/dL | SHBG, nmol/L | Free estradiol, pg/mL |
ESR1 | −0.07b | −0.45 | −0.01 | 0.47 | −0.48 |
0.67c | 0.0021 | 0.96 | 0.0011 | 0.0009 | |
HSD11B1L | −0.07 | −0.13 | −0.36 | −0.21 | −0.1 |
0.64 | 0.4 | 0.016 | 0.17 | 0.5 | |
SHBG | 0.31 | −0.06 | 0.2 | −0.23 | −0.04 |
0.036 | 0.68 | 0.2 | 0.12 | 0.78 | |
HSD17B7 | −0.02 | −0.33 | 0.08 | 0.16 | −0.34 |
0.91 | 0.028 | 0.59 | 0.3 | 0.023 | |
Inflammation-related genes | Glucose, mg/dL | ||||
ADIPOQ | 0.37 | NA | NA | NA | NA |
0.012 | |||||
IL1B | −0.33 | ||||
0.026 | |||||
IL6 | −0.38 | ||||
0.011 | |||||
IL8 | −0.39 | ||||
0.008 | |||||
Adipose tissue lipid metabolism–related genes | Leptin, ng/mL | Insulin, μU/mL | |||
ABCA1 | −0.35d | −0.35 | NA | NA | |
0.022 | 0.018 | ||||
Other genes | Leptin, ng/mL | Insulin, μU/mL | Glucose, mg/dL | Adiponectin, μg/mL | |
IGF1 | −0.04d | −0.3 | 0.2 | 0.07d | |
0.78 | 0.045 | 0.18 | 0.65 | ||
IGFBP3 | −0.22d | 0.02 | −0.16 | 0.49d | |
0.15 | 0.92 | 0.31 | 0.0011 | ||
LEP | 0.34d | 0.24 | 0.38 | −0.32d | |
0.028 | 0.11 | 0.01 | 0.039 | ||
RBP4 | 0.11d | 0.24 | 0.47 | −0.09d | |
0.5 | 0.11 | 0.0012 | 0.58 |
Pathwaya . | Biomarkers . | ||||
---|---|---|---|---|---|
Sex steroid–related genes | Estrone, pg/mL | Estradiol, pg/mL | Testosterone, ng/dL | SHBG, nmol/L | Free estradiol, pg/mL |
ESR1 | −0.07b | −0.45 | −0.01 | 0.47 | −0.48 |
0.67c | 0.0021 | 0.96 | 0.0011 | 0.0009 | |
HSD11B1L | −0.07 | −0.13 | −0.36 | −0.21 | −0.1 |
0.64 | 0.4 | 0.016 | 0.17 | 0.5 | |
SHBG | 0.31 | −0.06 | 0.2 | −0.23 | −0.04 |
0.036 | 0.68 | 0.2 | 0.12 | 0.78 | |
HSD17B7 | −0.02 | −0.33 | 0.08 | 0.16 | −0.34 |
0.91 | 0.028 | 0.59 | 0.3 | 0.023 | |
Inflammation-related genes | Glucose, mg/dL | ||||
ADIPOQ | 0.37 | NA | NA | NA | NA |
0.012 | |||||
IL1B | −0.33 | ||||
0.026 | |||||
IL6 | −0.38 | ||||
0.011 | |||||
IL8 | −0.39 | ||||
0.008 | |||||
Adipose tissue lipid metabolism–related genes | Leptin, ng/mL | Insulin, μU/mL | |||
ABCA1 | −0.35d | −0.35 | NA | NA | |
0.022 | 0.018 | ||||
Other genes | Leptin, ng/mL | Insulin, μU/mL | Glucose, mg/dL | Adiponectin, μg/mL | |
IGF1 | −0.04d | −0.3 | 0.2 | 0.07d | |
0.78 | 0.045 | 0.18 | 0.65 | ||
IGFBP3 | −0.22d | 0.02 | −0.16 | 0.49d | |
0.15 | 0.92 | 0.31 | 0.0011 | ||
LEP | 0.34d | 0.24 | 0.38 | −0.32d | |
0.028 | 0.11 | 0.01 | 0.039 | ||
RBP4 | 0.11d | 0.24 | 0.47 | −0.09d | |
0.5 | 0.11 | 0.0012 | 0.58 |
aOnly genes with at least one statistically significant association (P < 0.05) are shown.
bCorrelation coefficient.
cP value.
N = 45, except dN = 42.
Discussion
This study is the first to measure gene expression profiles in subcutaneous adipose tissue with an intervention study that examined the impact of dietary weight loss, exercise, and dietary weight loss plus exercise, compared with control. There was a change in gene expression in 3 genes, HSD17B1, LEP, and ESR1, in tissue obtained from women in both the diet and diet + exercise interventions. However, these effects by intervention arm did not remain statistically significant after the adjustment for multiple comparisons for 82 candidate genes and 15 candidate pathways. This finding may be due to the small sample size in this pilot study or the rather stringent adjustment. The only other intervention study to date to examine a change in gene expression with a dietary intervention in humans reported a decrease in LEP and ADIPOQ gene expression with a very low-calorie diet in 24 obese men (n = 18) and women (n = 6; ref. 45). LEP expression appears to be consistently altered with weight loss or caloric restriction, as a similar reduction in expression was noted in whole mouse mammary tissue with caloric restriction compared to ad libitum (30).
In gene set analyses, the IL-6 pathway was significantly modified by the diet and exercise intervention (adj. P = 0.003), whereas jak-STAT signaling was affected by diet alone with marginal significance (adj. P = 0.06). We observed more statistically significant signals in the analyses by weight change.
In an analysis combining women from all groups, a decrease in expression of genes involved in estrogen synthesis was noted with greater weight loss, with a statistically significant decrease in expression of HSD17B1 after accounting for multiple comparisons. This is the first report that illustrates the association of weight change with genes involved in synthesis of sex hormones in adipose tissue. HSD17B1 is a critical component of the estrogen metabolism pathway because it catalyzes the conversion of less-active estrone to estradiol.
We hypothesized that a reduction in body weight and fat mass would reduce adipose tissue inflammation, as supported by animal studies (46, 47), as well as human interventions including reduced or very low-calorie diets, or bariatric surgery (28, 29, 48). Surprisingly, we found no association between weight change and the expression of the vast majority of inflammatory cytokines and chemokines, including TNFα, IL-6, and MCP-1. Among the few inflammatory genes that showed a trend toward reduced expression from baseline to 6 months with weight loss, none reached statistical significance after adjustment for multiple comparisons, although some changes of substantial size with greater weight loss were observed, such as SAA1 (P = 0.01, adj. P = 0.13), where there was a reduction of expression of 38% among those who lost >10%. One potential explanation is that the weight loss achieved in this subset of participants was too small to have a sufficient effect on the inflammatory process. This is surprising given that the women in the diet and diet + exercise groups lost >12% of body fat and reductions in serum levels of inflammatory markers observed in all intervention groups (33). In a recent report from the CALERIE study, Tam and colleagues (49) also found no impact of a 25% weight reduction with diet or diet + exercise on genes relating to inflammation in subcutaneous adipose tissue among obese individuals. Our data, along with the findings of Tam and colleagues (49), strongly suggest that adipose tissue inflammation may be less impacted by changes in energy balance than has previously been hypothesized.
In the gene set analysis, the unsupervised clustering of more than 37,000 transcripts by percentage of weight change revealed 83 transcripts with statistically significant adjusted P values for trend. Using IPA to examine overlaps with known biologic functions and canonical pathways, several hits were noted in genes related to differentiation of adipocytes, the mTOR signaling pathway, and the IGF-1 signaling pathway. We also used GSEA to determine overlap between the genes that had statistically significant changes with weight loss and canonical pathways and GO categories. This analysis supported our IPA as we observed several hits for mRNA metabolism, the mTOR pathway, and others. Interestingly, LEP, the strongest signal in our candidate gene analysis, was the 42nd strongest hit in this clustering. This reaffirms the importance of leptin but also suggests that there are numerous other biologic pathways that are altered by weight change that require further investigation.
The ability to examine the association between serum biomarkers and gene expression is a novel aspect of our study. At baseline, serum levels of estradiol were associated with the expression of genes involved in synthesis of steroid hormones. Furthermore, serum levels of adiponectin, leptin, and SAA were positively associated with gene expression of ADIPOQ, LEP, and SAA1, respectively. Associations between serum insulin, glucose, and adiponectin and genes related to adipose tissue lipid metabolism were also observed. Furthermore, changes in serum levels with the intervention and changes in gene expression were examined. The interpretation of these findings is more difficult but does show the interplay between serum glucose and gene expression of several inflammatory genes, along with a negative association between serum leptin and insulin and gene expression of ABCA1, which is involved in adipose tissue lipid metabolism. While these findings suggest that serum levels of biomarkers may serve as a proxy for gene expression at the level of adipose tissue in some cases, more research is needed to understand the complex interplay of pathways at the level of the human adipose tissue.
Our study combined hypothesis-driven approaches at both the candidate gene and candidate pathway level with empirical strategies, using unsupervised clustering of more than 37,000 transcripts. Overall, the candidate gene/candidate pathway approach had a more than 10% success rate, even after adjustment of the level of multiple comparisons made (82 candidate genes). The unsupervised clustering (“-omic strategy”) revealed a total of 0.2% hits among all 37,000 transcripts. This approach highlights a number of new signals that require further follow-up. Notably, the very stringent level of multiple comparison adjustment for this approach would have missed significant changes in 8 of our candidate genes for which we had seen a statistically significant effect in the hypothesis-driven analysis. This suggests that the 2 approaches—hypothesis-driven and empirical—are complementary because the “–omic strategy” alone would result in a significant number of such false-negative results. When based on prior biologic plausibility, significant effects are to be expected when investigated via a targeted approach.
A major strength of this study is the investigation of effects of weight change directly in the target tissue. Here, a number of important processes (including estrogen biosynthesis and inflammatory processes) take place which influence overall metabolic states and cancer risk. We have shown that subcutaneous adipose tissue collection and gene expression analysis is feasible and has use in cancer prevention studies. We observed the strongest associations with gene expression changes when stratifying by weight loss with data from all intervention groups combined. An advantage of the study is the randomized design which included different weight loss groups (diet, exercise, diet + exercise) and a control group. In addition, the interventions used in this study are consistent with lifestyle interventions which could be adopted by the general population as a cancer prevention strategy (50), in contrast to the very-low-calorie diet (i.e., 450 kcal/d) examined in the study by Franck and colleagues (45).
There were several weaknesses in our study. First, the sample size for analyses by intervention arm was too small to reveal any statistically significant results after adjustment for multiple comparisons. Despite the pilot study nature of the investigation and small sample size, including a small number of participants in the diet and control groups, we were able to identify a number of prehypothesized and new candidates. We only assayed one adipose tissue biopsy sample from one anatomic site; subcutaneous abdominal fat. Gene expression could vary within the same fat depot and when compared with other fat depots, such as visceral, liver, or intramuscular fat. Visceral fat may be more biologic active for some analytes than adipose at other sites (51, 52). Furthermore, a target tissue of interest for breast cancer prevention is breast tissue and our study could not address the microenvironment in the breast itself. Therefore, the reported changes in gene expression in subcutaneous fat by intervention group or weight loss category may not be representative of changes in other biologically relevant tissues, such as breast or visceral fat. Our sample was limited to postmenopausal women ages 50 to 75 years, most of whom were non-Hispanic whites, and therefore our findings may not apply to men, women of different ages, or persons of other race/ethnic groups. Finally, we tested only one weight loss dietary plan and only one exercise program and cannot extrapolate with other methods of weight loss.
In conclusion, our study showed that diet- or exercise-induced weight loss results in measurable and sizable changes in adipose tissue gene expression, particularly in sex hormone steroid synthesis, leptin, and insulin signaling. Notable was the unexpected lack of effect on inflammatory pathways. The combination of a hypothesis-driven and empirical approach was particularly fruitful in that it supported hypothesized changes while also suggesting novel signals that require confirmatory study.
Disclosure of Potential Conflicts of Interest
A. McTiernan has ownership interest (including patents) in Merck and is a consultant/advisory board member of Metagenics. No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
Conception and design: K.L. Campbell, K.E. Foster-Schubert, K.W. Makar, M. Kratz, A. McTiernan, C.M. Ulrich
Development of methodology: K.L. Campbell, K.E. Foster-Schubert, K.W. Makar, M. Kratz, L.-Y. Kuan, A. McTiernan, C.M. Ulrich
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): K.L. Campbell, K.E. Foster-Schubert, K.W. Makar, D. Hagman, E.A. Schur, C. Abbenhardt, L.-Y. Kuan, A. McTiernan
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): K.L. Campbell, K.E. Foster-Schubert, K.W. Makar, M. Kratz, N. Habermann, M. Horton, L. Xiao, J. Davison, M. Morgan, C.-Y. Wang, C. Duggan, C.M. Ulrich
Writing, review, and/or revision of the manuscript: K.L. Campbell, K.E. Foster-Schubert, K.W. Makar, M. Kratz, D. Hagman, E.A. Schur, N. Habermann, M. Horton, L.-Y. Kuan, C.-Y. Wang, C. Duggan, A. McTiernan, C.M. Ulrich
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): K.L. Campbell, K.W. Makar, N. Habermann, M. Horton
Study supervision: K.L. Campbell, K.E. Foster-Schubert, C. Duggan, A. McTiernan, C.M. Ulrich
Grant Support
This study was supported by funding from NIH R21 CA131676-01(Principal Investigator: C.M. Ulrich) and NCI U54 CA116847 (Principal Investigator: K.E. Foster-Schubert).
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.