![]() |
|
|
4 Department of Pediatrics and 5 Ilyssa Center for Molecular and Cellular Endocrinology, Johns Hopkins University School of Medicine, Baltimore, MD 21287 and 6 Gene Expression and Genomics Unit, National Institute on Aging, NIH, Baltimore, MD 21224
* To whom correspondence should be addressed. E-mail: rsmiller{at}jhmi.edu.
| ABSTRACT |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Given the high prevalence of obesity, there has been interest in identifying factors associated with successful weight reduction. Most human studies have demonstrated that obese individuals typically return to pretreatment weight within 5 y of successful weight reduction, with a substantial proportion exceeding pretreatment weight (6–8). Expansion of the adipocyte population and secretion of small molecules from adipocytes may be involved in the failure to maintain weight loss; however, many factors responsible for rebound weight gain after weight reduction have yet to be identified.
Adipocyte-secreted factors influence diverse processes such as appetite and energy balance, lipid metabolism, insulin sensitivity, and inflammation (9,10) and visceral adiposity has been shown to be more directly related to cardiovascular risk than subcutaneous fat (11,12). We hypothesized that the visceral fat depot may be differentially responsive to diet composition and may produce factors not yet identified that may contribute to difficulty in maintaining weight loss. Using an oligonucleotide microarray, we performed a comprehensive evaluation of adipocyte gene expression in obesity as a function of diet composition (high fat vs. high carbohydrate). We then characterized gene expression in adipocytes derived from weight-reduced mice, identifying genes that may be involved in rebound weight gain following weight reduction.
| Materials and Methods |
|---|
|
|
|---|
We recorded body weight on arrival and at weekly intervals in the morning. Food intake was recorded twice weekly. All mice were given commercially prepared diets (Research Diets; Supplemental Table 1). After 1 wk of acclimatization, mice were divided into 4 groups: mice fed a low-fat diet (LF;7 10% energy as fat; catalog no. D12450B); mice fed a high-fat diet (HF; 60% energy as fat; no. D12492); mice fed a high-carbohydrate diet (HC) consisting of LF supplemented with a liquid formula, prepared by dissolving the powdered LF in water to a concentration of 103.5 g/L, thus providing 1674 kJ/L; and mice fed a HF/restricted diet (HFR) provided in sufficient quantity to match total body weight to the HC group. For wk 1, mice in the HFR group were given a quantity of food calculated, based in part on pilot studies, to provide the same energy we estimated would be consumed by the mice fed the HC diet. We weighed mice weekly and adjusted the amount of food provided to match total body weight in the HFR and HC groups. Fresh formula was changed daily and provided to mice in the HC group in place of drinking water to increase total energy intake. After 13 wk, HF-fed mice were subdivided into 2 groups. One-half continued to consume the HF ad libitum and the other half was switched to a LF diet and restricted to 2 g per mouse per day to induce weight reduction (WR).
We used a total of 50 mice for this experiment (10 mice per group). Microarray and RT-PCR data are reported on 42 mice. Data on the remaining 8 mice are not available due to inadequately prepared labeled cRNA or degraded RNA.
Insulin tolerance testing and fasting insulin concentrations.
Insulin tolerance testing was performed at wk 17. After food deprivation overnight, mice were administered 0.75 units/kg insulin (recombinant human regular insulin; Eli Lilly) intraperitoneally (i.p.). Venous tail blood was obtained for glucose determination by glucometer (Accu-Chek Advantage; Roche Diagnostics) at baseline and at 20, 40, and 60 min after insulin injection. The following week, following overnight food deprivation, mice were anesthetized with Avertin (2, 2, 2 tribromoethanol in tert-amyl alcohol; Sigma-Aldrich) 0.5 mg/g i.p.
100 µL blood was obtained by mandibular bleed, centrifuged at 4000 x g; 5 min at 4°C and serum was stored at –80°C for subsequent insulin concentration determined by RIA (Linco Research).
Adipocyte isolation and RNA purification. After 18 wk, mice were killed by CO2 asphyxiation and gonadal fat pads were harvested and weighed. Adipocytes were purified from adipose tissue according to established methods (15). Gonadal fat pads were weighed and transferred to vials containing 6 mL Krebs ringer bicarbonate HEPES buffer with 5% bovine fraction V albumin (Millipore), 200 nmol/L adenosine (Sigma), and 10 mg collagenase type 1 powder (Worthington Biochemical). Tissue was minced with sharp scissors and samples were incubated in a shaking 37°C water bath for 1 h. Samples were passed through a 500-µm Nitex mesh (Sefar Filtration) into 50-mL conical tubes, centrifuged at 200 x g for 10 min at room temperature and washed in Krebs ringer bicarbonate HEPES/albumin/adenosine (KRBH) buffer followed by 2 cycles of centrifugation at 200 x g for 1 min at room temperature and washing with KRBH buffer.
Immediately after the final wash, total RNA was extracted from adipocytes using Trizol reagent (Invitrogen) according to the manufacturer's protocol. For each sample, 200 µL of adipocytes was homogenized in 1 mL of Trizol and remaining cells were frozen at –80°C. The quality of total RNA samples was assessed using an Agilent BioAnalyzer (Agilent Technologies) and was equivalent across all samples analyzed.
Oligonucleotide microarray. A 0.5-µg aliquot of total RNA from each sample was labeled using the Illumina TotalPrep RNA Amplification kit (Ambion) as previously described (16). A total of 0.85 µg of biotin-labeled cRNA was hybridized for 16 h to Illumina's Sentrix MouseRef-8 Expression BeadChips (Illumina). The arrays were washed, blocked, and then hybridized. Biotinylated cRNA was detected with streptavidin-Cy3 and quantitated using Illumina's BeadStation 500GX Genetic Analysis Systems scanner. Image processing and data extraction were performed using BeadStudio ver. 15 (Illumina) according to previously described methods (17).
RT-PCR confirmation. RT-PCR was performed with the iScript One-Step RT-PCR kit with SYBR Green (Bio-Rad) using the MyiQ Single Color Real-Time PCR Detection system with an iCycler thermocycler (Bio-Rad). RNA for RT-PCR was obtained from the same samples used for the microarray. Reactions were optimized for each primer set (Supplemental Table 2; 18) and PCR products were separated on a 1% agarose gel to verify product size. For every reaction, a standard curve was generated using serial dilutions of mouse adipocyte RNA. Primer concentration was 200 nmol/L for all reactions. All reactions included a dissociation curve analysis at the end of the amplification to confirm a single product at the expected melting temperature. Relative quantitation of gene expression was performed using the threshold cycle. Differences in mRNA expression were calculated after normalizing to 18S expression.
Data analysis. Values for weights, glucose, and insulin concentrations are expressed as means ± SD. Results of RT-PCR are expressed as relative differences in gene expression. Statistical analyses were performed using a t test for RT-PCR data. Repeated-measures ANOVA or 1-way ANOVA with Student-Newman-Keuls post test was performed for all other comparisons. Statistical analyses were performed using GraphPad InStat version 3.0 for Windows NT (GraphPad Software). Differences were considered significant at P < 0.05.
Microarray data were analyzed using DIANE 6.0, a spreadsheet-based microarray analysis program. Raw microarray data were subjected to Z normalization and tested for significant changes as previously described (19). Genes were determined to be differentially expressed after calculating the Z-ratio and false discovery rate (20). Individual genes with P-value < 0.01, Z-ratio > 2, and false discovery rate < 0.3 were considered significantly changed. Ingenuity pathways analysis (IPA) (Ingenuity Systems) was used to identify biological networks with the greatest number of differentially expressed genes. For each network or pathway, probability scores were calculated using the right-tailed Fisher's exact test. A full discussion of the statistical methodologies employed can be found in Calvano et al. (21)
| Results |
|---|
|
|
|---|
|
|
Following insulin administration, whole blood glucose concentrations were lower at 20 min in LF compared with HFR, HF, HC, and WR (Fig. 1B; P < 0.01), suggesting greater insulin sensitivity in LF than in all other groups. This physiologic difference between LF and WR indicates there was incomplete recovery from prior obesity. Whole blood glucose concentrations 60 min after insulin administration were lower in LF, WR, and HFR (P < 0.05) compared with HF. The latter comparison demonstrates there was an effect of degree of obesity on insulin sensitivity in these mice.
Gene expression. We identified genes differentially expressed in HF vs. LF, HFR vs. HF, HFR vs. LF, HFR vs. HC, WR vs. LF, and WR vs. HF (Table 2; Supplemental Table 3; Supplemental Summary of Analysis). IPA demonstrated altered gene expression in a broad range of cellular functions (Supplemental Table 4). Interestingly, the most highly involved pathway involved PPAR signaling, which was downregulated in WR compared with LF.
|
|
We generated a list of genes that were differentially expressed in both WR vs. LF and HF vs. LF. We found that there were 268 genes upregulated and 130 genes downregulated in HF that were similarly differentially expressed in WR. Thus, for a large number of genes, altered expression in diet-induced obesity persisted even after weight reduction had been achieved and maintained in steady state. We used IPA to identify altered pathways (Supplemental Table 4) and to sort genes by location to compile a list of gene products that localize to the extracellular compartment. This resulted in identifying 59 genes differentially expressed in HF and WR (compared with LF) that code for extracellular proteins (Table 3).
|
|
| Discussion |
|---|
|
|
|---|
Several studies have identified differences in physiology and metabolism in high-fat vs. high-carbohydrate–induced obesity in rodents (3,24,25). Using a comprehensive microarray, we discovered that differences in diet composition affect the expression of a large number of genes. The clustering of gene expression data demonstrates that diet and pattern of adipocyte gene expression are associated. There are many signals originating from adipocytes that play a role in total body insulin sensitivity and hypothalamic control of appetite and the patterns of adipocyte gene expression we have identified may play a role in these complex biological processes.
Our data confirm that changes in degree of obesity are associated with changes in physiologic measures such as fasting insulin concentrations and insulin sensitivity. Not surprisingly, we found that as obesity increases, the altered expression of a number of genes that occurs in moderate obesity (HFR) persists into severe obesity (HF). However, we also determined that there are a large number of genes whose expression levels in moderately obese mice are not different from those in lean mice but are altered only in severe obesity. PCA and hierarchical clustering also demonstrate distinct segregation of HF and HFR mice by adipocyte gene expression. Thus, increasing obesity is associated with qualitative changes in gene expression, with unique patterns of gene expression in HFR and HF.
By merging WR and HF gene expression lists, we were able to identify factors that could play a role in failure to maintain weight loss after energy restriction. By the end of the study, WR mice had maintained weight in a steady state equivalent to the weight of lean controls for 2 wk. Despite near normalization of measured physiologic variables, WR mice displayed a pattern of adipocyte gene expression markedly different from that of either the obese, HF group or the lean, LF mice. PCA and hierarchical clustering revealed broad differences in adipocyte gene expression among the WR, HF, and LF mice. It remains unclear why a number of obesity-related genes continue to have altered expression following weight reduction. A number of different epigenetic phenomena may be involved, including DNA methylation or DNA/protein acetylation (26).
A number of genes we identified as persistently altered following weight reduction have previously been described as playing a role in appetite and adipose tissue development, thus confirming the physiological importance of our findings. Interestingly, a number of factors in this list are cytokines that have been implicated in cachexia related to illness and inflammation. These findings are consistent with previous studies by Clement et al. (27), which showed differential expression of 100 inflammation-related transcripts in adipose tissue of obese individuals that had been on a very low-energy diet for 1 mo. A number of factors that remain upregulated following weight loss are known to act locally in obesity. For example, colony-stimulating factor 1 has been shown to be upregulated in adipose tissue of human subjects who gain weight with overfeeding (28). Matrix metalloproteinases 2, 12, and 14 are strongly induced in adipose tissue in obesity (29) and matrix metalloproteinase 2 is elevated in obese mice (30).
As many signals related to appetite and satiety originate from adipose tissue, we hypothesized that genes differentially expressed in HF and remaining differentially expressed after diet could play a role in the rebound weight gain that often follows dieting and weight loss. Two proteins in Table 3 have been shown to affect appetite in rodents. Tenascin c (TNC) is an extracellular matrix protein expressed in the central nervous system. TNC-knockout mice were found to have lower neuropeptide Y (NPY) mRNA concentrations in limbic structures of the brain (31). As NPY is known to be orexigenic, one could speculate that increased circulating concentrations of TNC could result in increased NPY concentrations and increased appetite. Additional studies are required to determine whether these factors are consistently elevated postobesity and if they act on the hypothalamus to alter appetite.
To our knowledge, this study is unique in demonstrating that in mice of equivalent body weight, differences in diet composition lead to broad differences in adipocyte gene expression. Altered adipocyte gene expression as a consequence of diet composition may be adaptive to a diet high in saturated fat but may also be responsible for increasing the risk of diabetes and other obesity-related disorders. Although other investigators have shown changes in gene expression in adipocytes and adipose tissue in obesity (32–34), we have demonstrated differences in metabolism and adipocyte gene expression between formerly obese mice and weight-matched controls. Furthermore, this study has demonstrated that there are genes with altered expression in obesity that do not correct following weight reduction. Although a small number of adipocyte-secreted factors are known to be involved in appetite and energy balance, our findings reveal a number of adipocyte-secreted factors with persistently altered expression following weight reduction that may play a role in rebound weight gain following weight loss.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
2 Author disclosures: R. Miller, K. Becker, V. Prabhu, and D. Cooke, no conflicts of interest. ![]()
3 Supplemental Tables 1–4, Supplemental Figures 1 and 2, and a Summary of Analysis are available with the online posting of this paper at jn.nutrition.org. ![]()
7 Abbreviations used: HC, high-carbohydrate/-energy; HF, high-fat; HFR, high-fat–restricted; i.p., intraperitoneally; IPA, ingenuity pathways analysis; KRBH, krebs ringer bicarbonate HEPES/albumin/adenosine; LF, low-fat; WR, weight-reduced; NPY, neuropeptide Y; PCA, principal components analysis; TNC, tenascin c; TNF, tumor necrosis factor. ![]()
Manuscript received 17 September 2007. Initial review completed 17 October 2007. Revision accepted 21 March 2008.
| LITERATURE CITED |
|---|
|
|
|---|
1. Cowie CC, Rust KF, Byrd-Holt DD, Eberhardt MS, Flegal KM, Engelgau MM, Saydah SH, Williams DE, Geiss LS, et al. Prevalence of diabetes and impaired fasting glucose in adults in the U.S. population: National Health and Nutrition Examination Survey 1999–2002. Diabetes Care. 2006;29:1263–8.
2. Ogden CL, Carroll MD, Curtin LR, McDowell MA, Tabak CJ, Flegal KM. Prevalence of overweight and obesity in the United States, 1999–2004. JAMA. 2006;295:1549–55.
3. Dourmashkin JT, Chang GQ, Gayles EC, Hill JO, Fried SK, Julien C, Leibowitz SF. Different forms of obesity as a function of diet composition. Int J Obes (Lond). 2005;29:1368–78.[Medline]
4. Stubbs RJ, Johnstone AM, Mazlan N, Mbaiwa SE, Ferris S. Effect of altering the variety of sensorially distinct foods, of the same macronutrient content, on food intake and body weight in men. Eur J Clin Nutr. 2001;55:19–28.[Medline]
5. Willett WC. Is dietary fat a major determinant of body fat? Am J Clin Nutr. 1998;67:S556–62.[Abstract]
6. Douketis JD, Macie C, Thabane L, Williamson DF. Systematic review of long-term weight loss studies in obese adults: clinical significance and applicability to clinical practice. Int J Obes (Lond). 2005;29:1153–67.[CrossRef][Medline]
7. Miller WC. How effective are traditional dietary and exercise interventions for weight loss? Med Sci Sports Exerc. 1999;31:1129–34.[Medline]
8. Wadden TA. Treatment of obesity by moderate and severe caloric restriction. Results of clinical research trials. Ann Intern Med. 1993;119:688–93.
9. Rajala MW, Scherer PE. Minireview: The adipocyte–at the crossroads of energy homeostasis, inflammation, and atherosclerosis. Endocrinology. 2003;144:3765–73.
10. Trayhurn P. Endocrine and signalling role of adipose tissue: new perspectives on fat. Acta Physiol Scand. 2005;184:285–93.[CrossRef][Medline]
11. Arner P. Not all fat is alike. Lancet. 1998;351:1301–2.[CrossRef][Medline]
12. Pouliot MC, Despres JP, Lemieux S, Moorjani S, Bouchard C, Tremblay A, Nadeau A, Lupien PJ. Waist circumference and abdominal sagittal diameter: best simple anthropometric indexes of abdominal visceral adipose tissue accumulation and related cardiovascular risk in men and women. Am J Cardiol. 1994;73:460–8.[CrossRef][Medline]
13. Collins S, Martin TL, Surwit RS, Robidoux J. Genetic vulnerability to diet-induced obesity in the C57BL/6J mouse: physiological and molecular characteristics. Physiol Behav. 2004;81:243–8.[CrossRef][Medline]
14. Surwit RS, Feinglos MN, Rodin J, Sutherland A, Petro AE, Opara EC, Kuhn CM, Rebuffe-Scrive M. Differential effects of fat and sucrose on the development of obesity and diabetes in C57BL/6J and A/J mice. Metabolism. 1995;44:645–51.[CrossRef][Medline]
15. Smith U, Sjostrom L, Bjornstorp P. Comparison of two methods for determining human adipose cell size. J Lipid Res. 1972;13:822–4.[Abstract]
16. Cheadle C, Becker KG, Cho-Chung YS, Nesterova M, Watkins T, Wood W III, Prabhu V, Barnes KC. A rapid method for microarray cross platform comparisons using gene expression signatures. Mol Cell Probes. 2007;21:35–46.[CrossRef][Medline]
17. Gunderson KL, Kruglyak S, Graige MS, Garcia F, Kermani BG, Zhao C, Che D, Dickinson T, Wickham E, et al. Decoding randomly ordered DNA arrays. Genome Res. 2004;14:870–7.
18. Rozen S, Skaletsky H. Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000;132:365–86.[Medline]
19. Cheadle C, Vawter MP, Freed WJ, Becker KG. Analysis of microarray data using Z score transformation. J Mol Diagn. 2003;5:73–81.
20. Hochberg Y, Benjamini Y. More powerful procedures for multiple significance testing. Stat Med. 1990;9:811–8.[Medline]
21. Calvano SE, Xiao W, Richards DR, Felciano RM, Baker HV, Cho RJ, Chen RO, Brownstein BH, Cobb JP, et al. A network-based analysis of systemic inflammation in humans. Nature. 2005;437:1032–7.[CrossRef][Medline]
22. Zhou YP, Ling ZC, Grill VE. Inhibitory effects of fatty acids on glucose-regulated B-cell function: association with increased islet triglyceride stores and altered effect of fatty acid oxidation on glucose metabolism. Metabolism. 1996;45:981–6.[CrossRef][Medline]
23. Roche E, Maestre I, Martin F, Fuentes E, Casero J, Reig JA, Soria B. Nutrient toxicity in pancreatic beta-cell dysfunction. J Physiol Biochem. 2000;56:119–28.[Medline]
24. Kahn BB, Pedersen O. Suppression of GLUT4 expression in skeletal muscle of rats that are obese from high fat feeding but not from high carbohydrate feeding or genetic obesity. Endocrinology. 1993;132:13–22.
25. Kanarek RB, Marks-Kaufman R. Developmental aspects of sucrose-induced obesity in rats. Physiol Behav. 1979;23:881–5.[CrossRef][Medline]
26. Feinberg AP, Cui H, Ohlsson R. DNA methylation and genomic imprinting: insights from cancer into epigenetic mechanisms. Semin Cancer Biol. 2002;12:389–98.[CrossRef][Medline]
27. Clement K, Viguerie N, Poitou C, Carette C, Pelloux V, Curat CA, Sicard A, Rome S, Benis A, et al. Weight loss regulates inflammation-related genes in white adipose tissue of obese subjects. FASEB J. 2004;18:1657–69.
28. Levine JA, Jensen MD, Eberhardt NL, O'Brien T. Adipocyte macrophage colony-stimulating factor is a mediator of adipose tissue growth. J Clin Invest. 1998;101:1557–64.[Medline]
29. Chavey C, Mari B, Monthouel MN, Bonnafous S, Anglard P, Van Obberghen E, Tartare-Deckert S. Matrix metalloproteinases are differentially expressed in adipose tissue during obesity and modulate adipocyte differentiation. J Biol Chem. 2003;278:11888–96.
30. Ahima RS, Qi Y, Singhal NS, Jackson MB, Scherer PE. Brain adipocytokine action and metabolic regulation. Diabetes. 2006;55 Suppl 2:S145–54.
31. Fukamauchi F, Aihara O, Kusakabe M. Reduced mRNA expression of neuropeptide Y in the limbic system of tenascin gene disrupted mouse brain. Neuropeptides. 1998;32:265–8.[CrossRef][Medline]
32. Baranova A, Collantes R, Gowder SJ, Elariny H, Schlauch K, Younoszai A, King S, Randhawa M, Pusulury S, et al. Obesity-related differential gene expression in the visceral adipose tissue. Obes Surg. 2005;15:758–65.[CrossRef][Medline]
33. Moraes RC, Blondet A, Birkenkamp-Demtroeder K, Tirard J, Orntoft TF, Gertler A, Durand P, Naville D, Begeot M. Study of the alteration of gene expression in adipose tissue of diet-induced obese mice by microarray and reverse transcription-polymerase chain reaction analyses. Endocrinology. 2003;144:4773–82.
34. Urs S, Smith C, Campbell B, Saxton AM, Taylor J, Zhang B, Snoddy J, Jones VB, Moustaid-Moussa N. Gene expression profiling in human preadipocytes and adipocytes by microarray analysis. J Nutr. 2004;134:762–70.
This article has been cited by other articles:
![]() |
G. J. Hausman, M. V. Dodson, K. Ajuwon, M. Azain, K. M. Barnes, L. L. Guan, Z. Jiang, S. P. Poulos, R. D. Sainz, S. Smith, et al. BOARD-INVITED REVIEW: The biology and regulation of preadipocytes and adipocytes in meat animals J Anim Sci, April 1, 2009; 87(4): 1218 - 1246. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Nairizi, P. She, T. C. Vary, and C. J. Lynch Leucine Supplementation of Drinking Water Does Not Alter Susceptibility to Diet-Induced Obesity in Mice J. Nutr., April 1, 2009; 139(4): 715 - 719. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||