DOI: 10.2337/diabetes.55.03.06.db05-1015 © 2006 by the American Diabetes Association
Necdin and E2F4 Are Modulated by Rosiglitazone Therapy in Diabetic Human Adipose and Muscle Tissue
1 Research Division, Joslin Diabetes Center, Boston, Massachusetts Address correspondence and reprint requests to Allison B. Goldfine, MD, Joslin Diabetes Center, One Joslin Place, Boston, MA 02215. E-mail: allison.goldfine{at}joslin.harvard.edu
Abbreviations:
CEBP-
To identify novel pathways mediating molecular mechanisms of thiazolidinediones (TZDs) in humans, we assessed gene expression in adipose and muscle tissue from six subjects with type 2 diabetes before and after 8 weeks of treatment with rosiglitazone. mRNA was analyzed using Total Gene Expression Analysis (TOGA), an automated restriction-based cDNA display method with quantitative analysis of PCR products. The expression of cell cycle regulatory transcription factors E2F4 and the MAGE protein necdin were similarly altered in all subjects after rosiglitazone treatment. E2F4 expression was decreased by 10-fold in muscle and 2.5-fold in adipose tissue; necdin was identified in adipose tissue only and increased 1.8-fold after TZD treatment. To determine whether changes were related to an effect of the drug or adipogenesis, we evaluated the impact of rosiglitazone and differentiation independently in 3T3-L1 adipocytes. While treatment of differentiated adipocytes with rosiglitazone did not alter E2F4 or necdin, expression of both genes was significantly altered during differentiation. Differentiation was associated with increased cytosolic localization of E2F4. Moreover, necdin overexpression potently inhibited adipocyte differentiation and cell cycle progression. These data suggest that changes in necdin and E2F4 expression after rosiglitazone exposure in humans are associated with altered adipocyte differentiation and may contribute to improved insulin sensitivity in humans treated with TZDs.
Thiazolidinediones (TZDs) effectively treat insulin resistance and type 2 diabetes. TZDs bind to peroxisome proliferator–activated receptor-
PPAR- To better understand the molecular mechanisms by which TZDs act in humans with type 2 diabetes, we sought to identify the genes differentially expressed in human adipose and skeletal muscle in vivo after exposure to the TZD rosiglitazone. We used Total Gene Expression Analysis (TOGA), an automated restriction-based cDNA display technique that offers advantages over other expression analysis methods. First, in contrast to high-density oligonucleotide arrays, this technology allows for the identification and unbiased evaluation of both known and novel mRNAs and splice variants. Second, TOGA is highly sensitive, detecting transcripts expressed at 1 in 106, whereas microarray sensitivity generally ranges from 1 in 105 to 1 in 104 transcripts. Therefore, TOGA analysis requires lower quantities of total RNA (20 ng), making this technique particularly suitable for small human biopsy samples. We have applied this technology for the first time to assess TZD-related responses in humans. We confirmed differential gene expression identified by TOGA using real-time quantitative PCR and performed functional analysis in cultured 3T3-L1 adipocytes.
The Joslin Diabetes Center Human Subject Committee approved the experimental protocol, and informed consent was obtained from all participants. For the study, six subjects with type 2 diabetes not previously treated with a TZD or insulin were recruited. All subjects had normal coagulation and liver function and had no other major systemic illness. Blood samples were obtained from the subjects in the fasted state for the analysis of insulin, glucose, and cholesterol levels and liver function. After subjects were given a local anesthetic with 1% lidocaine, percutaneous biopsies of abdominal subcutaneous adipose tissue and vastus lateralis muscle were performed using a triport cannula and Bergstrom needle, respectively. Samples were immediately frozen in liquid nitrogen. Subjects were then treated with 8 mg/day rosiglitazone (Avandia; GlaxoSmithKline) for 8 weeks, after which the adipose tissue and muscle biopsies were repeated. Compliance with the medication regimen, assessed by pill count, was >94%. (Although techniques are under development to assess TZD drug levels [13], they are not yet widely available for clinical or research use.) We chose the 8-week treatment period because it allowed us to identify gene expression changes that preceded major changes in glycemia, due to potentially confounding effects of glycemic changes on gene expression (14,15).
RNA isolation and TOGA analysis. TOGA analysis was performed (16) on mRNA isolated from adipose and muscle samples. cDNA was synthesized using a degenerate pool of biotinylated phasing primers that initiated synthesis at the poly(A) tail. cDNA was digested with MspI; 3' fragments were captured using streptavidin beads and released by cleavage with NotI, which recognized a site in the phasing primers. Captured fragments were modified at their 5' ends to harbor a start site for T3 polymerase in vitro transcription. RNA fragments were produced corresponding to the 3' portion of each starting mRNA, from its most 3' MspI recognition sequence to the beginning of its poly(A) tail, with each fragment flanked by linker tags of known sequence. cDNA was prepared from the RNA pool by reverse transcription and used in four separate PCRs, in which a 5' primer that extends by one of four possible nucleotides beyond the MspI site (N1 position) was paired with a universal 3' primer to generate an N1-specific double-stranded DNA template. Finally, 256 primers corresponding to all possible permutations of four nucleotides immediately adjacent to the MspI recognition site (N1, N2, N3, and N4) were matched with the appropriate N1 template and used in individual PCRs to produce 256 nonoverlapping product pools that were separated by electrophoresis. This assigned each PCR product an address consisting of an eight-nucleotide sequence (the four MspI recognition and adjacent four parsing nucleotides) and length, both of which are attributes of individual mRNAs. The fluorescent PCR product peak amplitudes correspond to initial concentrations of parent mRNAs. Data were queried to identify mRNAs whose concentrations differed among experimental samples. Each sample was analyzed in duplicate, independent reactions. cDNA clones were isolated from peaks of interest and sequenced.
Real-time quantitative RT-PCR.
Animal studies. Protocols were approved by the institutional animal care and use committee. ICR mice (Taconic Laboratories) were housed in an Office of Laboratory Animal Welfare–certified facility under standard light/dark cycle conditions. For the experimental group of mice, rosiglitazone was added to Purina 9F Chow at 3 mg · kg–1 · day–1 for 4 months (n = 7 per group).
Cultured cells.
The pRc/CMV vector (Invitrogen) subcloned with the 1.6-kb murine necdin fragment was ligated into the EcoR1-linearized pBABE retroviral construct. Then 3 µg of necdin-pBABE DNA or empty vector control DNA were transfected into
Protein analysis.
Immunocytochemistry.
Triglyceride assays.
Cell cycle analysis.
Statistical analysis.
Human studies. The subjects (four men and two women; five Caucasian and one African American) had a mean age of 49.5 ± 10.8 years, BMI 26.5 ± 3.3 kg/m2, HbA1c 6.9 ± 1.4%, and diabetes duration 1.9 ± 1.6 years. There was no significant change before and after 2 months of rosiglitazone treatment in weight (80.8 ± 14.5 vs. 83.9 ± 17.9 kg; P = 0.1) or fasting glucose (8.6 ± 2.3 vs. 7.6 ± 1.7 mmol/l; P = 0.1), insulin (91.7 ± 31.9 vs. 75.0 ± 33.3 pmol/l; P = 0.3), or C-peptide (0.9 ± 0.3 vs. 0.7 ± 0.3 nmol/l; P = 0.1) levels; all trends were of the expected direction and magnitude. No single subject appeared as a nonresponder across all metabolic variables.
Expression analysis. To assess the reproducibility of TOGA, a correlation coefficients matrix was generated from duplicate assays (Matlab, Natick, MA). Pearsons product moment correlation analysis demonstrated high reproducibility of replicates (r2 > 0.97) within adipose and muscle tissue. Despite the different sex and ethnic backgrounds of the subjects, the expression of most RNA species within adipose and muscle tissue did not significantly differ among subjects, as was demonstrated by the between-subject similarity scores (r2 = 0.84–0.98).
In four of six subjects, more peaks were altered in adipose than in muscle tissue, as might be predicted from the higher abundance of PPAR- The expression of one mRNA was decreased in all subjects after rosiglitazone in both muscle and fat tissue, with a mean reduction in muscle of 81 ± 9% (P = 0.001) and in adipose tissue of 49 ± 10% (P = 0.008) (Fig. 1A–C). Sequencing of the cloned PCR product identified this mRNA as transcription factor E2F4 (Genbank U15641). Decreased E2F4 mRNA expression was confirmed by real-time PCR in five of six subjects (due to limited tissue in one subject).
A second RNA species increased after rosiglitazone in adipose tissue of all subjects by 1.8 ± 0.2–fold (P = 0.008). Expression was low in muscle, with no net change after rosiglitazone ((Fig. 2). The product was identified as necdin (GenBank XM_007686.1). Expression differences were validated in adipose from two subjects due to limited mRNA availability. Thus we also analyzed necdin expression in mouse epididymal fat. Treatment of mice with rosiglitazone (3 mg · kg–1 · day–1) for 4 months increased necdin expression by 36% compared with controls (P = 0.003; data not shown).
Analysis in cultured adipocytes. The observed differential expression of E2F4 and necdin in humans after rosiglitazone treatment could be a direct result of the drug and/or PPAR- –regulated gene expression or an indirect result of alterations in the differentiation state or systemic changes in metabolic control or insulin sensitivity. To test the impact of rosiglitazone in the absence of complicating systemic metabolism, we evaluated the expression of E2F4 and necdin as a function of differentiation or drug treatment in 3T3-L1 adipocytes. Expression of E2F4 mRNA decreased significantly by day 2 of differentiation and remained low thereafter (75% reduction at day 2, P < 0.006; 56% at day 12, P = 0.03) (Fig. 3A). E2F4 protein expression in Triton-insoluble nuclear fractions also decreased with differentiation (92% reduction at day 2; P < 0.001) (Fig. 3B). Although the decrease in E2F4 protein expression was not prominent at day 6, it remained low thereafter (91% decrease at day 12; P < 0.001). We also observed the transient appearance of an additional 130-kDa band at day 2 that was not evident in subsequent days (data not shown); this could have represented nuclear translocation of a related E2F4-immunoreactive protein. In Triton-soluble cytosolic extracts, E2F4 expression increased during early differentiation (69% increase, days 0–4; P < 0.05) but declined by late differentiation (66% decrease, days 4–12; P < 0.01).
Because alterations in subcellular localization of E2F4 play an important role in cell differentiation (19), we evaluated E2F4 expression using immunocytochemistry (Fig. 3C). By day 2 of differentiation, nuclear E2F4 expression was reduced with relative preservation of cytosolic expression, thus suggesting cytosolic translocation. As differentiation progressed, expression decreased further in the nucleus and cytosol. To determine whether E2F4 was directly regulated by rosiglitazone, we treated preadipocytes and fully differentiated (day 9) adipocytes with rosiglitazone for 72 h. Although rosiglitazone did not alter mRNA expression of E2F4 in fully differentiated adipocytes, there was a small reduction in preadipocytes (7% decrease, P = 0.04; data not shown), perhaps related to the limited induction of differentiation with rosiglitazone alone. Similarly, the E2F4 protein level was modestly reduced in rosiglitazone-treated preadipocytes (18% decrease in Triton-insoluble fraction, P = 0.03; 14% decrease in Triton-soluble fraction, P = 0.06) (Fig. 4 and data not shown, respectively) but was not altered by rosiglitazone in fully differentiated cells. In contrast, E2F4 protein expression was markedly reduced as a function of the differentiation state alone (96% decrease; P < 0.001) (Fig. 4). These findings suggest a predominant role for the differentiation state rather than a direct effect of rosiglitazone in regulating E2F4 expression in adipocytes.
Like E2F proteins, necdin has been implicated in the regulation of cell cycle progression (20) and differentiation in cultured neuroblastoma cells (21). Therefore, we evaluated the effects of differentiation on necdin expression in 3T3-L1 adipocytes. Differentiation was associated with progressively increased expression of necdin mRNA (40% increase at day 12; P = 0.01) (Fig. 5A). Necdin protein levels increased throughout differentiation in both Triton-soluble (57-fold increase; P < 0.001) (Fig. 5B) and Triton-insoluble (2.3 fold increase, P = 0.001; data not shown) extracts. Increased expression of necdin was confirmed by immunocytochemistry (Fig. 5C).
We similarly evaluated direct effects of rosiglitazone on necdin expression. In preadipocytes, rosiglitazone increased necdin mRNA expression, albeit less robustly than differentiation-related effects (23%, P = 0.03; data not shown). Likewise, protein levels in cytosolic extracts increased dramatically with differentiation (14-fold increase; P < 0.001) but not with rosiglitazone (Fig. 6A and B).
E2F (22) and necdin (20) play important roles in cell cycle control, and E2F proteins regulate adipocyte differentiation. Necdin increases differentiation in neuroblastoma cells (21) but may play an inhibitory role in brown preadipocytes (23). To determine a potential functional role of necdin in white adipose tissue differentiation, we generated preadipocytes stably overexpressing necdin (17.5-fold increase in mRNA expression and 10-fold increase in protein levels in 3T3-L1/pBABE-necdin compared with 3T3-L1/pBABE cells). These cells did not differentiate after standard differentiation (Fig. 7A, panel b), even in the presence of 10 µmol/l rosiglitazone (Fig. 7A, panel d). In a similar finding, the accumulation of triglycerides during differentiation was decreased in 3T3-L1/pBABE-necdin cells (Fig. 7E). Although necdin overexpression did not affect E2F4 mRNA (data not shown), necdin reduced mRNA expression of both early (PPAR- and CAAAT/enhancer-binding protein- [CEBP- ]) and late (adipocyte fatty acid–binding protein) markers of adipocyte differentiation (Fig. 8A). Likewise, necdin overexpression reduced PPAR- protein levels in both cytosolic and nuclear fractions during early differentiation (data not shown). Given the potent inhibition of differentiation in necdin-expressing cells, we assessed the expression of preadipocyte factor 1 (Pref-1 [delta-like 1 homolog]), an important and potent inhibitor of adipocyte differentiation whose expression is normally reduced early after differentiation induction (24). Expression of Pref-1 increased 2.6-fold in 3T3-L1/pBABE-necdin cells (Fig. 8A) and remained substantially higher during differentiation in necdin-expressing cells, even in the presence of rosiglitazone.
Because necdin affects cell cycle progression in neuroblastoma cells (21), we analyzed the impact of necdin on cell cycle progression in 3T3-L1 preadipocytes. Necdin overexpression increased the percentage of cells in G1 (1.8-fold increase at 12 h, P < 0.001 vs. controls; 1.4-fold increase at 24 h, P < 0.001 vs. controls) and decreased the percentage of cells in the G2/M (3.3-fold decrease at 12 h, P < 0.001; 1.5-fold decrease at 24 h, P < 0.001) and S (2.7-fold decrease at 12 h, P = 0.002; 1.2-fold decrease at 24 h, P = 0.15) phases (Fig. 8B).
Despite the widespread use of TZDs in the treatment of type 2 diabetes, the cellular and molecular mechanisms by which these compounds improve insulin sensitivity remain incompletely understood. To identify genes that have altered expression in vivo in response to rosiglitazone and thus are of potential importance in the pathogenesis of insulin resistance and diabetes, we used the novel TOGA methodology to analyze the expression of both known and unidentified mRNA in human muscle and adipose tissue. Although studies of gene expression as a function of TZD in humans are ongoing (25–27), we are the first to report on the global differential expression profiles in these tissues from humans exposed to a TZD. Despite the inherent genetic variability in the human subjects in this study, the expression of most RNA species within adipose tissue and muscle did not differ significantly among the subjects, as demonstrated by their high similarity scores. Overall, rosiglitazone-induced expression changes were small and varied substantially among subjects. The expression of only 74 genes was altered by >1.5-fold in muscle or adipose tissue in at least two subjects after 8 weeks of rosiglitazone. Both the number and magnitude of changes in expression were substantially less than those seen in microarray studies of cells or rodents treated with TZDs (28,29). Several technical and methodological factors might account for these differences: 1) the detection of only mRNA containing Msp I restriction sites in our TOGA analysis, 2) between-study differences in drug dosage or duration of exposure, 3) greater intersubject genetic and environmental variability inherent to human studies, and 4) our stringent criteria for gene selection (requiring consistent per-gene expression changes in all subjects). In addition, it is interesting that despite potent changes in insulin sensitivity induced by TZDs, we and others (26,27) do not find insulin-signaling proteins to be highly regulated after TZD exposure. We focused our analysis on the two genes for which concordant differential regulation has been observed in all human subjects as a function of TZD: E2F4 and necdin. Although these two RNA species are unrelated members of the E2F and MAGE families, both are important regulators of cell proliferation, differentiation, and postmitotic survival. Rosiglitazone markedly suppressed E2F4 expression in muscle and adipose tissue but increased necdin expression in adipose tissue only. We observed similar effects of differentiation on increasing necdin and decreasing E2F4 expression in 3T3-L1 adipocytes.
Potential role of E2F4 in TZD effects.
Consistent with the major effects of E2F proteins on cell cycle control, E2F proteins also play an important role in the regulation of adipocyte differentiation (22). E2F1 induces PPAR- Our studies demonstrated reduced E2F4 expression after rosiglitazone in human adipose and muscle tissue in vivo and showed that E2F4 expression and nuclear localization are decreased by both rosiglitazone and differentiation in 3T3-L1 cells. The magnitude of change in E2F4 expression is much greater with differentiation than with rosiglitazone; this may reflect the relatively modest ability of rosiglitazone to independently induce differentiation in 3T3-L1 preadipocytes. Thus in vivo changes in E2F4 expression may be due to direct actions of rosiglitazone on E2F4 expression or be related to prodifferentiation effects of rosiglitazone.
Potential role of necdin in TZD effects. Necdin is of particular interest for metabolic diseases because it maps to the 15q11–q13 region associated with Prader-Willi syndrome, characterized by transient perinatal hypotonia, global developmental delay, hypogonadotropic hypogonadism, short stature, hyperphagia, and severe obesity (41). It is not surprising that insulin resistance, impaired glucose tolerance, and diabetes are also common (42). Genetic abnormalities in Prader-Willi syndrome include paternal deletions, maternal disomy, and imprinting mutations resulting in inactivation of paternally expressed genes in this region, including necdin. The disruption of necdin in mice results in phenotypes with some similarities to human Prader-Willi syndrome, including postnatal lethality (43–45). Necdin-null mice surviving the neonatal period demonstrate behavioral changes and reduced hypothalamic neuronal development but no major obesity phenotype. Given the potential relation between necdin and obesity, it is of interest that a 1352T/C promoter and silent 2311C/T coding region polymorphism have been identified in human necdin; however, these variants have not been linked to BMI in obese adolescents (46). Necdin was initially identified by subtraction cloning in embryonic carcinoma cells induced to neuronal differentiation by retinoic acid (37); necdin expression increased markedly during early differentiation. Overexpression of necdin in neuroblastoma cells both suppresses cell growth through inhibition of S phase entry and cell cycle arrest (20,47) and induces differentiation (21).
The effect of necdin on differentiation may be specific to cell type. In contrast to neuronal cells, where growth arrest is required before final differentiation, adipocyte differentiation occurs with initial growth arrest, cell cycle re-entry, and clonal expansion and then a second growth arrest followed by terminal differentiation (48). Tseng et al. (23) recently demonstrated increased expression of necdin in insulin receptor substrate 1–deficient brown preadipocytes, which are incapable of differentiation. Reduction of necdin expression (by siRNA or reconstitution of insulin receptor substrate 1) restored the differentiation capacity, indicating an inhibitory effect of necdin on differentiation. These inhibitory effects of necdin were accompanied by reductions in E2F4-stimulated PPAR-
Our data suggest the importance of necdin in the differentiation of white adipose tissue, a key tissue that is much more abundant than brown adipose tissue in adult humans and that is critical for the maintenance of metabolic homeostasis. Stable expression of necdin in 3T3-L1 adipocytes reduced expression of prodifferentiation regulators CEBP- In addition, we demonstrated that necdin can directly regulate cell cycle progression in white adipocytes. The overexpression of necdin inhibits cell cycle progression, as indicated by increased numbers of cells in the G1 phase and decreased numbers in G2/M and S phases. Inhibition of S phase entry by necdin overexpression has previously been demonstrated in stably transfected NIH-3T3 cells and N1E-115 cells, resulting in cell cycle arrest (20,21,47). However, in our studies, stable overexpression of necdin inhibited both cell cycle progression and differentiation of 3T3-L1 preadipocytes, thereby suggesting that 1) necdin effects on differentiation are independent of direct cell cycle regulation and are perhaps mediated via the E2F4 and transcriptional interactions described above or 2) necdin overexpression in preadipocytes also inhibits the initial growth arrest and mitotic clonal expansion required for terminal differentiation.
Potential coordinate roles of necdin and E2Fs in adipocyte differentiation. Our data in 3T3-L1 adipocytes support a model by which rosiglitazone, together with other prodifferentiation factors, rapidly decreases expression and nuclear localization of E2F4, thereby reducing E2F4 inhibition of cell cycle entry and stimulating early postinduction mitosis and differentiation initiation. At the same time, decreased E2F4 expression/function may derepress necdin, which, in turn, may induce growth arrest and limit further differentiation. Taken together, the effects of rosiglitazone to decrease E2F4 and increase necdin expression in humans likely reflects TZD-induced enhancement of differentiation (9). We speculate that the effects on necdin and E2F4 expression suggest one mechanism that might limit adipose accumulation with continued TZD exposure. In summary, our study identified two genes regulated by TZD therapy in humans; we determined that in vivo rosiglitazone exposure decreases E2F4 expression and increases necdin expression. Although both rosiglitazone and the induction of differentiation decreased E2F4 and increased necdin expression, the differentiation state had a greater effect than did TZD exposure, suggesting that changes in human tissue after rosiglitazone exposure are indirect. These data support the concept that one primary mechanism by which TZD therapy acts is via effects on cell cycle control and adipocyte differentiation, resulting in increased numbers of small, well-differentiated adipocytes that may be more insulin sensitive (9). Moreover, direct modulation of these protein families may provide new opportunities for the treatment of insulin resistance and type 2 diabetes.
This study was supported by National Institutes of Health Grants K23-DK-02795, RO1-DK-33201, RO1-DK-45935, P30-DK-36836, M01-RR-01032, K08-DK-02526, R01-DK-62948, and RO1-DK-060837 (Diabetes Genome Anatomy Project) and the Iacocca Foundation. We thank C. Ronald Kahn and Yu-Hua Tseng for their helpful discussion.
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. A.B.G. and S.C. contributed equally to this work. D.L. is currently affiliated with the La Jolla Institute for Allergy and Immunology, San Diego, California. Additional information for this article can be found in an online appendix at http://diabetes.diabetesjournals.org. Received for publication August 8, 2005 and accepted in revised form November 18, 2005
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||