|
|
||||||||
1Division of Diabetes, Department of Medicine, University of Texas Health Science Center at San Antonio, Texas; 2Cleveland Clinic Foundation, Endocrinology, Diabetes and Metabolism, Cleveland, Ohio; and 3Arizona State University, School of Life Sciences, Tempe, Arizona
Submitted 18 September 2007 ; accepted in final form 5 March 2008
| ABSTRACT |
|---|
|
|
|---|
, was significantly downregulated. SGK and CITED2 are located at chromosome 6q23, where we previously detected strong linkage to fasting plasma insulin concentrations. We independently validated the mRNA expression changes in an additional five subjects and closely paralleled the results observed in the original 12 subjects. A saline infusion in healthy, normal glucose-tolerant subjects without family history of diabetes demonstrated that the genes altered during the euglycemic hyperinsulinemic clamp were due to hyperinsulinemia and were unrelated to the biopsy procedure per se. The results of the present study demonstrate that insulin acutely regulates the levels of mRNAs involved in inflammation and transcription and identifies several candidate genes, including HES1 and BHLHB2, for further investigation. gene expression; muscle; insulin action; euglycemic hyperinsulinemic clamp; inflammation
Despite intensive investigation, the underlying cause of insulin resistance remains unknown. Many cellular defects have been shown to contribute to the pathophysiology of insulin resistance in type 2 diabetes, including impaired insulin signaling (4, 10), reduced insulin-stimulated glucose uptake, decreased hexokinase II expression and activity, diminished glycogen synthase activity, and decreased pyruvate dehydrogenase activity (8, 15, 26, 30, 39, 42). Impaired mitochondrial function has been associated with an increase in toxic lipid metabolites (fatty acyl-CoA, diacylglycerol, ceramide) in skeletal muscle of subjects with prediabetes and type 2 diabetes, and these toxic lipid metabolites have been implicated as a cause of the insulin resistance (1, 20, 23, 27, 29). We (36) previously reported that a physiological increase in plasma FFA concentration in normal healthy subjects resulted in changes in gene expression in skeletal muscle that are consistent with an inflammatory response. In addition, we also demonstrated that changes in the composition of the extracellular matrix are a general characteristic feature of insulin-resistant muscle (7).
Insulin has been shown to be a risk factor for the development of atherosclerotic cardiovascular disease (11). The effects of insulin on arterial tissues include proliferation of smooth muscle cells, enhanced cholesterol synthesis and low-density lipoprotein-receptor activity, increased formation and decreased regression of lipid plaques, stimulation of connective tissue synthesis, and stimulation of growth factors (11). In the present study, we set out to test the hypothesis that short-term exposure (4 h) to physiological hyperinsulinemia in normal, glucose-tolerant, healthy subjects without a family history of diabetes would induce a low-grade inflammatory response, independently of glycemic status, in skeletal muscle. By combining the euglycemic insulin clamp technique (13) with microarray profiling, we tested this hypothesis and, at the same time, identified novel pathways that are regulated by insulin in skeletal muscle.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Study design.
All subjects consumed a weight-maintaining diet containing 50% carbohydrate, 35% fat, and 15% protein over the 3 days prior to study and were instructed not to participate in any heavy exercise. After a 10-h overnight fast, all subjects reported to the General Clinical Research Center (GCRC) at 0800, and catheters were placed into an antecubital vein (for infusion of all test substances) and retrogradely into a vein on the dorsum of the hand (for blood withdrawal). The hand was placed in a heated box (65°C) to achieve arterialization of the venous blood. Twelve subjects received a 4-h euglycemic hyperinsulinemic (80 mU·m–2·min–1) clamp, which was performed with infusion of [3-3H]glucose to determine rates of glucose appearance and disposal, as previously described (10, 13). Briefly, 120 min prior to the start of the insulin clamp, a primed (25 µCi) continuous (0.25 µCi/min) infusion of [3-3H]glucose (DuPont-NEN, Boston, MA) was started and continued throughout the clamp to determine rates of total body glucose disposal and endogenous (primarily hepatic) glucose production (13). During the 30 min prior to the start of the insulin infusion and throughout the insulin clamp, plasma glucose concentration was determined every 5 min, and a variable glucose infusion was adjusted on the basis of the negative feedback principle to maintain the plasma glucose concentration at each subject's basal value. Plasma insulin and plasma tritiated glucose specific activity were measured every 10–15 min. Under euglycemic hyperinsulinemic conditions, the majority (
85%) of glucose disposal occurs in skeletal muscle (12). Biopsies of the vastus lateralis muscle were taken at –60, 30, and 240 min of the insulin infusion, as previously described (10). A control saline (0.2 ml·m–2·min–1) infusion was performed in three additional subjects, with a vastus lateralis muscle biopsy taken basally and at the end of the saline infusion. After completion of the microarray/RT-PCR determinations (see below), we studied an additional five subjects to validate our original observations. These five subjects received a 3-h euglycemic hyperinsulinemic (80 mU·m–2·min–1) clamp with tritiated glucose, as described above. Vastus lateralis muscle biopsies were performed at –60 min and after 180 min of insulin infusion. Muscle biopsy specimens (
200 mg) were frozen immediately in liquid nitrogen and stored until processed.
Muscle biopsy processing. Muscle biopsies from each subject were homogenized in RNAStat solution (Tel-Test, Friendswood, TX), using a Polytron homogenizer (Brinkmann Instruments Westbury, NY). Total RNA was purified with RNeasy and DNase I treatment (Qiagen, Chatsworth, CA). RNA isolated from the insulin infusion biopsies was used in the microarray and quantitative RT-PCR studies described below. The RNA isolated from the control saline infusion biopsies was used for the quantitative RT-PCR studies.
Microarray analysis: target preparation, hybridization, staining, scanning, and image analysis. RNA was prepared for hybridization to Affymetrix (Santa Clara, CA) HG-U133A arrays according to the manufacturer's instructions. A description of the protocol has been described previously (36). The probe arrays were scanned twice and the stored images aligned and analyzed using the GeneChip Operating System software (GCOS; Affymetrix). GeneChip quality assessments using GCOS demonstrated that the average percent present call across all arrays was 28 ± 1%. Average scaling values and noise levels were 6.8 ± 0.9 and 4.7 ± 0.2, respectively. bioB, bioC, bioD and cre were present on all arrays with the exception of bioB, which was present only 50% of the time, as expected. 3'/5' GAPDH and actin expression ratios were well within the acceptable range: 1.2 ± 0.1 and 0.8 ± 0.1, respectively.
Microarray data expression and analysis. The Affymetrix data acquisition programs in GCOS automatically generate a cell intensity file (CEL) from the stored images that contain a single intensity value for each probe cell on the array. The CEL data files have been deposited in the Gene Expression Omnibus (GEO) database as series GSE9105 [NCBI GEO] . The CEL files were imported into the R software package (http://www.r-project.org), and the probe level data were converted to expression measures using the Affy package from Bioconductor. Expression values for each mRNA were obtained by GC-Robust Multi-Array Analysis (GCRMA). CEL files were normalized together, and the expression values obtained were submitted to analysis with linear models of microarray data (LIMMA). To correct for multiple testing, P values were adjusted using the method of Benjamini and Hochberg by setting the false discovery rate (FDR) to 5% and a fold change of >1.5.
Functional classification of the significantly increased or decreased genes was performed using database for annotation, visualization and integrated discovery (DAVID) (http://david.niaid.nih.gov), Onto-Express (http://vortex.cs.wayne.edu/ontoexpress/), and expression analysis systemic explorer (EASE; downloaded from DAVID). DAVID uses clustering algorithms to classify related genes into functionally related groups. EASE discovers categories of genes that are overrepresented in the significant gene list compared with the total set of genes represented on the microarray. Such overrepresented categories represent biological "themes" of a given list of significant genes. The categories are sorted into three major gene ontology (GO) groups: 1) molecular function, 2) biological processes, and 3) cellular components. Onto-Express was used in the present study to translate the list of differentially expressed insulin-regulated genes into chromosome location.
Analysis of published microarray data. We compared the results of our microarray analysis with other previously published datasets by using the Web-based approach described by Parikh et al. (28). Briefly, we obtained normalized data files that were downloaded from either the Diabetes Genome Anatomy Project (DGAP; http://www.diabetesgenome.org) or the GEO database (http://www.ncbi.nlm.nih.gov/geo/). Using the Mann-Whitney U-test, a nonparametric statistical technique to identify mRNA expression changes, we compared the results in the DGAP and GEO databases with our list of mRNAs that were significantly altered following insulin infusion.
Quantitative real-time PCR.
Muscle expression of various genes was determined using a one-step quantitative real-time PCR (QRT-PCR) from total RNA (Applied Biosystems 7900HT). The quantity of mRNA for each gene of interest for each sample was normalized to that of 18S ribosomal RNA, using the comparative (2–
CT) method as described previously (36). Statistical comparisons were done using paired t-tests. Sequences of the primers and probes used for mRNA quantification by QRT-PCR can be supplied upon request.
Substrate and hormone determinations. Plasma glucose concentration was determined by the glucose oxidase method on a Beckman Glucose Analyzer (Beckman Instruments, Fullerton, CA). Plasma insulin was determined by radioimmunoassay (Diagnostic Products, Los Angeles, CA). Plasma tritiated glucose specific activity was determined on Somogyi precipitates of plasma as previously described (5).
Statistical analysis. Data were expressed as the means ± SE. Statistical significance of the difference between means was determined using paired or nonpaired Student's t-tests where appropriate. Correlations were analyzed using Pearson's product moment test. For all analyses, P < 0.05 was considered to be statistically significant. Statistical analysis of microarray data and QRT-PCR data were described above.
| RESULTS |
|---|
|
|
|---|
Analysis of microarray data. Microarray analysis was performed on the muscle biopsies obtained from the 12 subjects who received a 4-h euglycemic hyperinsulinemic clamp (Fig. 1). Using a false discovery rate of <5% and a fold change of >1.5, there were no individual probe sets with significantly increased or decreased expression after 30 min of insulin infusion compared with the basal state. In contrast, 121 probe sets were significantly altered (117 upregulated and 4 downregulated) after 240 min of insulin infusion compared with the basal state [Supplementary Tables 1 and 2, respectively (see online)]. The 117 upregulated probe sets represent in total 93 genes, since most genes are represented more than once due to the presence of more than one probe set on the HG-U133A GeneChip. The expression of 105 probe sets (96 increased and 9 decreased) was altered after 240 min compared with their expression 30 min after the start of the insulin infusion (Supplementary Tables 3 and 4, respectively). The 96 probe sets with significantly increased expression at 240 vs. 30 min of insulin infusion also were significantly altered at 240 min vs. the basal state. We identified five additional probe sets that were decreased at 240 vs. 30 min of the insulin infusion (Supplementary Table 4). Of note, one of these genes was thioredoxin-interacting protein (TXNIP), which has been implicated in insulin-stimulated peripheral glucose uptake in humans (28). The results of the present study indicate that most of the detectable changes occurred 30 min after the start of insulin.
|
|
The chromosomal location of the differentially expressed genes following insulin infusion was determined using the program Onto-Express. The location of the insulin-regulated genes was distributed among all chromosomes except Y, 13, and 7. Statistical analysis performed in Onto-Express yielded no chromosome that contained more insulin-regulated genes than expected (Supplementary Fig. 5).
Insulin signaling gene analysis. Because the microarray analysis identified significant increases in some genes involved in the insulin-signaling pathway including RRAD, KRAS, KLF10, and FOS (Supplementary Table 1), we looked at the microarray expression values of other insulin-signaling genes. Following 240 min of insulin infusion, the expression of 24 insulin-signaling genes increased (all P < 0.05, without correcting for multiple testing; Table 2). Many of these mRNA expression changes did not achieve statistical significance when corrected for multiple testing.
|
Correlations of microarray data with subject characteristics. To identify associations between insulin-regulated gene expression and insulin-stimulated glucose disposal, we performed correlation analysis using Pearson's product moment test. Insulin-stimulated whole body (primarily represents muscle) glucose disposal was strongly and inversely correlated with ADH1B, CYR61, and BHLHB2 mRNA expression (r = –0.74, –0.65, and –0.57, respectively). We also observed a negative correlation between basal EGP and serum/glucocorticoid-regulated kinase (SGK) mRNA expression (r = –0.74, P < 0.01), while DNAJB1, HBEGF, and TRIB1 mRNA gene expression was strongly correlated with fasting plasma insulin concentration (r = 0.74, 0.65, and 0.58, respectively). These and other correlations are depicted in Supplementary Table 7.
QRT-PCR analysis to confirm microarray data. To validate the microarray data, we quantitated the changes in mRNA levels of 20 genes (18 of which increased and 2 of which decreased) using QRT-PCR in the twelve subjects who received a 4-h euglycemic hyperinsulinemic clamp. These genes displayed wide differences in their response to insulin on the microarray (from –2.3 to 7.7-fold change). The correlation coefficient (r) between the genes quantitated by QRT-PCR and microarray analysis was 0.79 (P < 0.01) (Supplementary Fig. 6). Of note, the QRT-PCR fold changes were larger than the changes observed with microarray analysis. This is consistent with previous studies from our laboratory in which microarray analysis tends to underestimate the changes in gene expression (36). In addition, we independently validated the mRNA expression changes in an additional five subjects who received a 3-h euglycemic hyperinsulinemic clamp. mRNA expression in these 5 subjects, who were studied after the microarray/QRT-PCR analyses were completed in the first 12 subjects, closely paralleled that in the original 12 individuals (Table 3). We also observed significant interindividual variation in the basal mRNA levels in these normal, healthy subjects. This natural variation in mRNA expression most likely represents the complex interplay of intrinsic (i.e., age, sex, ethnicity, and genes), physiological, and environmental factors. This observation is consistent with other studies that have identified significant interindividual variation in basal mRNA expression in healthy subjects (18, 33, 43).
|
|
| DISCUSSION |
|---|
|
|
|---|
560 pmol/l) achieved during the insulin clamp was within the high physiological range, as described elsewhere (32), and was similar to that observed during an oral glucose/meal tolerance test, although we infused insulin for 4 h, which is longer than the usual tissue exposure to insulin following a meal in normal, healthy subjects. Therefore, we believe that the results in the present study have relevance to the effects of physiological hyperinsulinemia on muscle mRNA expression and may provide insights into the etiology of insulin resistance in type 2 diabetes and obesity (11). The uniqueness of the present study compared with previous studies (3, 15, 21, 25, 28, 37, 40, 44) is the inclusion of a control saline infusion group and validation of the results in a separate group of individuals. In our study, we felt that it was important to demonstrate that the observed changes in mRNA expression during the insulin infusion clamp were independent of the biopsy procedure per se, which could have led to increased mRNA expression of some genes involved in inflammation. The results of the saline infusion study clearly demonstrate that the genes altered during the euglycemic hyperinsulinemic clamp were due to insulin and independent of the muscle biopsy procedure.
The results of the present study demonstrate that a number of genes involved in inflammation were upregulated. Thrombomodulin, which is dysregulated during inflammation, had the largest fold change (7.69) from the microarray analysis. A number of inflammatory chemokines (CCL8, CXCL2, and CCL2) and cytokines (interleukin-1 receptor type I and interferon-
-inducible protein-16) were also upregulated by insulin. Chemokine, cc motif ligand 2 (CCL2) and chemokine, cc motif ligand 2 (CCL8), also known as monocyte chemoattractant protein-1 and -2, are members of the CC chemokine family and promote monocyte chemotaxis to sites of inflammation. These results are consistent with those of a previous study that identified a number of inflammatory genes (including THBD, CCL2, and CXCL2) that were upregulated following insulin treatment (21). Rome et al. (37) also reported that a number of immune-response genes were increased following a euglycemic hyperinsulinemic clamp, but those genes did not overlap with our study (37). QRT-PCR confirmed the increased mRNA expression of inflammation genes originally identified by microarray analysis. Furthermore, the control saline infusion demonstrated that the increased expression of inflammatory genes in response to the insulin infusion was due to the action of insulin on the muscle and was not secondary to the muscle biopsy per se. Taken collectively, these results provide evidence that in skeletal muscle insulin augments the expression of multiple genes that are involved in inflammation. This inflammation response may then provide positive feedback to cause/enhance insulin resistance.
Several candidate gene and microarray studies have identified genes whose expression was altered following an insulin infusion (3, 15, 21, 25, 28, 37, 40, 44). We examined some of these previously published microarray datasets (21, 28, 37) and found a number of genes that changed in parallel with our insulin-regulated genes. The majority of genes whose expression was altered by insulin coded for proteins involved in transcription regulator binding and activity. Other groups have reported similar findings concerning the effect of insulin on the mRNA expression of transcription factors (21, 37). Hansen et al. (21) described seven transcription genes (ATF3, BHLHB2, HES1, KLF10, JUNB, FOS, and FOSB) whose mRNA expression was upregulated and in common with those in the present study. Of particular note, we found that increased mRNA expression of HES1 and BHLHB2 was replicated in the majority of the datasets analyzed (21, 28, 37). Reduced insulin-mediated glucose disposal in skeletal muscle is a common metabolic feature of type 2 diabetes (4, 10, 14). In the present study, we observed an inverse relationship between BHLHB2 expression and whole body (primarily reflects muscle) glucose Rd. Moreover, we observed that the mRNA expression of HES1 and BHLHB2 was increased in muscle tissue of type 2 diabetic subjects (27, 29). Collectively, the HES1 and BHLHB2 mRNA expression results suggest that these two genes may be involved in the insulin-resistant phenotype observed in type 2 diabetes. However, the exact mechanism of action of the HES1 and BHLHB2 genes in skeletal muscle requires further investigation.
Oxidative stress has been implicated in the development of insulin resistance and microvascular/macrovascular complications in patients with type 2 diabetes (9). Metallothioneins (MTs) are cysteine-rich, low-molecular-weight intracellular proteins that have been shown to regulate the metabolism of metals and to be efficient scavengers of reactive oxygen species (ROS) (6, 38). In the present study, functional analysis using DAVID identified a group of MT mRNAs, whose expression is increased following insulin stimulation. Furthermore, analysis of the Parikh et al. (28) dataset identified a significant increase in expression of the mRNA for MT1H during a 2-h euglycemic hyperinsulinemic clamp. We also observed an increased mRNA expression of MT1H in muscle tissue of type 2 diabetic subjects (27), although opposite findings have been reported elsewhere (40). Another study found no differences in MT mRNA levels when type 2 diabetic groups were compared with control groups (38). Taken together, these results suggest that insulin activates the MT antioxidant defense system in normal, healthy controls, but further studies are needed to establish what role, if any, altered MT mRNA expression plays in the pathogenesis of human type 2 diabetes and its associated complications.
Insulin modulates intracellular metabolism by regulating pathways involved in protein, glycogen, and fat metabolism. In skeletal muscle, insulin stimulates glucose uptake, glycogen synthesis, and glucose oxidation (14). In the present study, we show that insulin increased the expression of a number of mRNAs in the phosphoinositol (PI) 3-kinase (metabolic) and MAP kinase (mitogenic) pathways (Table 2). In addition to the insulin receptor IRS2, regulatory/catalytic subunits of PI 3-kinase and hexokinase II, in particular, the ras-related associated with diabetes gene (RRAD) was increased following insulin stimulation and similar findings have been reported by other groups (15, 25). RRAD is overexpressed in the skeletal muscle of type 2 diabetes subjects (35, 40), and has been implicated in the pathogenesis of insulin resistance in type 2 diabetes. RRAD is involved in the classical MAP kinase pathway, which is an integral component of the insulin-signaling cascade and has been implicated in the development of atherosclerosis (22). We hypothesize that the increase in RRAD observed by previous investigators and in the present study may represent a response to the compensatory hyperinsulinemia that occurs as the β-cell attempts to offset the insulin resistance. RRAD may provide a link between insulin resistance, hyperinsulinemia, and atherosclerosis (11).
Insulin, as well as insulin-like growth factor, serum, and glucocorticoids, stimulate SGK activity by a mechanism requiring the participation and activation of PI 3-kinase (24, 31). We demonstrate in the present study that muscle mRNA expression of SGK is significantly increased by insulin in normal healthy subjects and that its expression was negatively correlated with basal EGP. To our knowledge, this result has not been replicated in any other microarray study. However, it has been shown elsewhere that SGK transcription is markedly enhanced in diabetic nephropathy (24). Interestingly, the SGK gene has been localized to chromosome 6q23, a major locus associated with fasting insulin concentrations and insulin resistance in Mexican Americans from the San Antonio Family Diabetes Study (16). These results suggest a potential role for SGK in insulin resistance in type 2 diabetes and other insulin resistant states.
cAMP-responsive element binding protein (CREB)/p300-interacting transactivator, with Glu/Asp-rich carboxy-terminal domain 2 (CITED2) also is located at 6q23, and in the present study was downregulated by insulin. CITED2 interacts with the CREB-binding protein (CBP)/p300, TFAP2, Lhx2, and nuclear receptors, such as peroxisome proliferator-activated receptor and estrogen receptor, to function as a transcriptional modulator (41). Another chromosomal region that has been linked with insulin resistance, type 2 diabetes and obesity is 10q24–25 (17). The peroxisome proliferative activated receptor-
coactivator-related 1 (PPRC1) is located at 10q24.32 and was found to be significantly upregulated by insulin in the present study. PPRC1 was first identified in database searches for sequences with similarities to peroxisome proliferator-activated receptor (PPAR)
coactivator-1 (PGC1) (2). PGC1 is a transcriptional regulator that activates (in conjunction with nuclear respiratory factor-1 and PPAR
), PPAR
, hepatocyte nuclear factor 4, and other transcription factors which are key metabolic regulators. PPRC1 is a growth-regulated coactivator that may coordinate the activities of multiple transcription factors required for cell growth (2). The exact role of the CITED2 and PPRC1 genes remains unclear, and further studies will be needed to understand their role in insulin signaling.
In summary, we have demonstrated that a sustained physiological increase in the plasma insulin concentration for as little as 4 h in nonexercised subjects alters the levels of multiple muscle mRNAs, including inflammatory factors, which may be involved in the development of insulin resistance in type 2 diabetes mellitus. Of note, several of these genes are localized to chromosomal regions that previously have been shown to be associated with type 2 diabetes, insulin resistance, obesity, and atherosclerosis.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. Granjon, M.-P. Gustin, J. Rieusset, E. Lefai, E. Meugnier, I. Guller, C. Cerutti, C. Paultre, E. Disse, R. Rabasa-Lhoret, et al. The microRNA Signature in Response to Insulin Reveals Its Implication in the Transcriptional Action of Insulin in Human Skeletal Muscle and the Role of a Sterol Regulatory Element-Binding Protein-1c/Myocyte Enhancer Factor 2C Pathway Diabetes, November 1, 2009; 58(11): 2555 - 2564. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |