1 Division of Endocrinology, Metabolism, and Lipid Research, Department of Internal Medicine, Washington University School of Medicine, St. Louis, Missouri
2 Genome Sequencing Center, Washington University School of Medicine, St. Louis, Missouri
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Pancreatic islet ß-cells can be regulated by multiple stimuli, including nutrients and growth factors. ß-Cell proliferation and function are controlled by plasma glucose concentration and by growth factors acting via multiple intracellular signaling pathways (1). Changes in gene expression that result from the activation of these signaling pathways are likely responsible for the adaptation of ß-cells to physiological and pathological states. However, large gaps in our knowledge currently exist regarding the changes in gene expression and the molecular mechanisms mediating these ß-cell responses to nutrients and growth factors.
Some genes likely to be involved in chronic glucose regulation of islet ß-cell mass or function have been identified (26). We have focused on early signaling events initiated by glucose treatment of insulinoma cells that result in rapid transient activation of a number of immediate early genes (IEGs). These include Egr1, Egr2, c-fos, and c-jun, known to respond to growth factor stimulation in a number of other tissues (7). The signaling pathways for induction of IEGs exhibit considerable stimulus and tissue specificity and in general involve activation of kinase/phosphatase cascades (8). Initial glucose-mediated signaling can represent the first step in elucidating long-term changes in gene expression and islet physiology. These signaling pathways, limited to the initial kinase/phosphatase cascades, are critical for understanding how the ß-cell responds to its environment. The events occurring from the time the stimulus reaches the ß-cell until the signal is transmitted to the nucleus to activate or repress transcription of a particular set of genes may be crucial in understanding the defects in islet growth in diabetic subjects or the adverse consequences of glucose toxicity on ß-cell function.
IEGs are often transcription factors that in turn activate expression of downstream target genes, thus generating distinct biological responses by inducing specific long-term programs of gene expression. In pancreatic ß-cells, activation of expression of these IEGs was shown to depend on depolarization activation of voltage-gated Ca2+ channels and subsequent influx of extracellular Ca2+. This resulted in activation of Ca2+-regulated kinases, including calmodulin-dependent kinase IV and protein kinase A, leading to phosphorylation and activation of several transcription factors (cAMP-responsive element binding protein, serum response factor, and Elk-1) (9,10). The results of these experiments defined the rapid glucose-signaling pathways for a small number of IEGs whose transcription is rapidly activated by glucose, but the results now pose additional questions addressed by the current study.
Animal models perfused with glucose for 45 days, or transgenic animals overexpressing a particular gene (11), result in more readily measured physiological changes, yet the sequence of molecular events leading to these physiological changes are difficult to discern. This result highlights the desirability of beginning to dissect these mechanisms using other models. Thus, we designed experiments using insulinoma cells to elucidate early transcriptional responses to islet growth factors. The results of the present experiments extend our knowledge of IEGs regulated by glucose, by KCl-induced depolarization, and by insulin through use of high-resolution custom cDNA microarrays that contain clones from the Endocrine Pancreas Consortium (EPCon: http://www.cbil.upenn.edu/ EPConDB). The arrays used for these experiments contain up to 9,700 cDNAs with >3,000 novel clones not currently available on commercial arrays (1214). The results of this work suggest that glucose activation of IEGs is mediated primarily via depolarization and that glucose and insulin activate an overlapping set of genes. Further, both of these growth stimuli appear to activate transcription through a phosphatidylinositol (PI) 3-kinase-dependent pathway. These results further illustrate how monitoring expression gene profiles can serve to elucidate important ß-cell biological functions.
![]() |
RESEARCH DESIGN AND METHODS |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
Microarray construction.
After the sequencing of cDNA libraries derived from various mouse pancreas tissues gathered by the EPCon (12,13), expressed sequence tags (ESTs) were evaluated for similarity to the existing nonredundant GeneBank entries (http://www.ncbi.nlm.nih.gov/Genbank/index.html) and for redundancy within the EST clone set using standard BLAST analyses (17).
A nonredundant set of clones was selected and PCR amplified in preparation for the microarray slide production. PCR products were purified using Millipore 96-well plates, DNA quantity was determined using Picogreen (Molecular Probes) reagents, and, finally, each of the PCR products was normalized with Millipore purified water to a standard concentration of 400 ng/µl. Before microarray spotting, all PCR products were adjusted to a final DNA concentration of 200 ng/µl in a spotting buffer of 3x sodium chloride-sodium citrate (SSC) and 0.75 mol/l betaine. PCR products were spotted onto epoxy-coated slides (MWG Biotech) using an ArrayMaker2 arrayer (designed by P. Browns laboratory; Stanford University). Spotted slides were incubated 1214 h at 40% humidity in a 42°C oven. After incubation, the slides were cross-linked using the UV Stratalinker 2400 (Stratagene) at 700 µJ x100 and processed by the following steps: 1) gently shaking the slides in a 0.2% SDS bath for 2 min; 2) three room temperature water bath washes, each at 1 min; 3) a 50°C water bath for 20 min; 4) a 95°C water bath for 2 min; and 5) spinning the slides dry. The slides were stored in a desiccant cabinet for future hybridization experiments.
The different clones on the array were identified through BLASTn similarity matches using the NCBI BLAST (http://www.ncbi.nlm.nih.gov/blast) and WU-BLAST version 2.0 (17) against the nonredundant subsets of the public Mouse databases and RefSeq (http://www.ncbi.nlm.nih.gov/RefSeq), and only matches with an E value below at least 1 x 1050 and a score superior to 200 were considered as a match. Annotation was further confirmed comparing the results to the annotation provided on the EPCon website (http://www.cbil.upenn.edu/EPConDB) for the PancChip 5.0 and manually curating discrepancies. Further annotation (official names and symbols, gene ontology [GO] functions) was gathered from Source (http://source.stanford.edu) and in some cases imported directly from LocusLink (http://www.ncbi.nlm.nih.gov/LocusLink).
Hybridization.
Two hybridizations were carried out in a sequential manner. The primary hybridization was performed by adding 38 µl of sample to the microarray under a supported glass coverslip (Erie Scientific) at 50°C for 1620 h at high humidity in the dark. Before the secondary hybridization, slides were gently submerged into 2x SSC, 0.2% SDS for 12 min; transferred to 2x SSC for 12 min; transferred to 0.2x SSC for 12 min; transferred to 95% ethanol for 2 min; and then spun dry by centrifugation. Secondary hybridizations were carried out using the complimentary capture reagents provided in the 3DNA Array 50 kit (Genisphere).
Data acquisition.
Slides were scanned on a Perkin Elmer ScanArray Express HT scanner to detect Cy3 and Cy5 fluorescence. Laser power was kept constant for Cy3/Cy5 scans, and photo-multiplier tube was varied for each experiment based on background fluorescence. Gridding and analysis of images was performed using QuantArray (Perkin Elmer). Each spot was defined on a pixel-by-pixel basis using a modified Mann-Whitney statistical test.
Pairing scheme.
To compare each of the conditions with every other condition, samples were paired as indicated in Fig. 1B. For each pair, transcripts from one condition were labeled with either Cy3 or Cy5 and hybridized with the Cy5- or Cy3-labeled transcripts from the partner condition. Another pair of RNAs was also labeled inverting dyes ("dye-flip" hybridization). This sample pairing and dye flip required a total of 12 microarrays per experiment (each arrow in Fig. 1B represents one slide of microarray with a pair of RNAs). As a result, each sample for an experiment was hybridized six times: against three other conditions twice, once with Cy3, and the other with Cy5. RNA from the first experiment was hybridized to two sets of microarray, a 5,700 clone set and then to a 9,700 clone microarray consisting of a superset of the 5,700 clone set; therefore, in each case, 12 slides were used. Therefore, a total of 24 microarray slides was used for these analyses.
Normalization.
In the microarray studies performed herein, the probes that gave a signal intensity 2.0-fold above the corresponding background intensity in both Cy5 and Cy3 channels were chosen for calculation of the fold change. The intensity of each spot was first adjusted by subtracting background intensity. Then, the log ratio of Cy5 and Cy3 channel intensity of each spot was calculated, and the median of the log ratio from all probes was subtracted from each log ratio (the global normalization procedure [18]). This is based on the assumption that the expression levels of most genes are unchanged, and thus log fold change should center at zero. Indeed, data from all hybridizations indicated that the expression of the majority of genes (
95%, data not shown) were unchanged.
The log ratios from six pairs of dye-flip hybridizations were used to estimate probe-specific measurement variation, usually derived from self-vs.-self experiments (19). Probe-specific measurement variation was subtracted from log fold change to calculate the final normalized log fold change (probe-specific normalization).
Statistical analysis.
From our pairing scheme, two direct and two indirect ratios were obtained when two conditions were compared. When glucose-stimulated and -unstimulated samples were compared, for example, two ratios were acquired from dye-flip hybridizations directly comparing those two samples. Two additional ratios were calculated indirectly from eight hybridizations involving all glucose, KCl, insulin, and unstimulated samples. These four ratios were used to calculate average fold change as well as standard error. The 95% CI was calculated with attention to the change of distribution when converting log fold change to fold change (20). Criteria were set to assess changes in gene expression that included a fold change of 1.3 or more compared with "unstimulated." This criterion was selected after observing that the mean variance of all probes derived from a single self-to-self hybridization from the first 5,700 and 9,700 probe sets was 0.12 (Fig. 1C). Because we had four measurements, that mean variance was further reduced to 0.06 (= 0.12/4). With this small variance, the false-positive rate is 0.0032%. With similar multiple hybridization and cDNA microarray, a fold change of 1.3 was used to identify significant gene expression change in other studies (2123). Additionally, significant fold changes required that the 95% CI of the fold change, two standard errors away from the mean fold change, excluded 1. Thus, any fold changes meeting these criteria could be considered significant with P < 0.05.
Hierarchical clustering.
Hierarchical clusters were performed using the Genesis software version 1.3 (Institute for Biomedical Engineering, Graz University of Technology, Graz, Austria) (24). Heat maps generated by hierarchical clustering were created using the average linkage clustering (Euclidian distances) directly from the logarithmic ratios (base e), and all color ranges for the "heat maps" were adjusted to a maximum ratio of 1.
Quantitative real-time RT-PCR.
Total RNA was isolated, and 1 µg was used to prepare cDNA, primed with random hexamers, and reverse transcribed with Superscript II (Invitrogen) according to manufactures protocol. Quantitative RT-PCR (qRT-PCR) was performed by monitoring in real time the increase in fluorescence of the SYBR Green dye (ABI) as described (25,26) using the ABI 7000 sequence detection system (Applied Biosystems). For comparison of transcript levels between samples, a standard curve of cycle thresholds for serial dilutions of a cDNA sample was established and then used to calculate the relative abundance of each gene. Values were then normalized to the relative amounts of 18S ribosomal RNA, which were obtained from a similar standard curve. All PCRs were performed at least in replicates of four. Standard error of the quantity of transcript normalized to the amount of 18S ribosomal RNA was calculated from a formula with consideration of error propagation. When gene expression levels of two conditions were compared, the ratio was expressed with standard error calculated from the same formula. Specificity of each primer pair was confirmed by melting curve analysis and agarose gel electrophoresis of PCR products. Sequences of primers used in this study are included in an online appendix at http://diabetes.diabetesjournals.org.
![]() |
RESULTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Gene expression in response to glucose.
Stimulation by glucose resulted in a significant increase in mRNA levels for 72 of the 9,700 cDNAs on the arrays. The 50 genes most increased by glucose treatment are shown in Table 1. The complete list of regulated genes is shown in Table 1S in the online appendix. Glucose-stimulated genes of diverse GO functions, described more fully below, included two immediate early genes previously identified as glucose responsive in islet ß-cells, early growth response 1 (Egr1) (10,28), and an inhibitor of DNA binding 1 (Idb1) (29). Several other immediate early genes found to respond to mitogenic stimuli in other tissues (immediate early response 2 [Ier2] and 3 [Ier3] and inhibitor of DNA binding 2 [Idb2]) were also shown to respond to glucose treatment. Interestingly, among the glucose-regulated genes, the overwhelming majority was not previously known to be glucose responsive. Glucose treatment also resulted in significant downregulation of 21 genes (Table 1), of which only one, Chop (Ddit3), was previously described as glucose regulated (3). Chop has been incriminated in cell cycle arrest and apoptosis (30). At least 13 of the glucose-responsive genes do not yet have a match in LocusLink or UniGene.
|
|
|
Examination of the genes whose expression was significantly affected by insulin treatment revealed a surprising result, i.e., there appeared to be considerable overlap between the genes regulated by insulin and by glucose, as can be readily seen in Table 1, where glucose- and insulin-activated genes are compared (bold). The expression profiles after the two treatments are also graphically assessed in Fig. 2B by a hierarchical cluster analysis. Immediately apparent is that 8 of the top 10 genes activated by glucose were not at all stimulated by insulin treatment. On the other hand, most of the remainder of the glucose-regulated genes was expressed to the same extent by insulin treatment.
Validation of gene expression profiles by quantitative RT-PCR (qRT-PCR) and by biological replicates.
The observation by microarray analysis that genes appeared to be activated by both glucose/depolarization and insulin required validation by independent means. This was performed by qRT-PCR on RNA samples from a subset of genes that were shown to be activated on the microarrays (experiment 1). As can be seen in Table 3, in general, the fold changes noted with qRT-PCR were larger than those observed with microarray, yet there was a highly significant correlation (R2 = 0.90, data not shown) between the levels of gene expression measured by the two methods for experiment 1.
|
Assessment of the PI 3-kinase pathway in depolarization and insulin-mediated gene expression.
The observation of a common set of genes activated by glucose and insulin was consistent with glucose/depolarization and insulin inducing gene transcription via completely separate pathways that target the same genes or that the two stimuli activate pathways that converge on a common signaling pathway that activates the same genes. To differentiate between the two possibilities, gene expression was evaluated in another set of experiments using a pharmacological inhibitor (LY294002) for PI 3-kinase, known to be involved in insulin signaling (32). KCl and insulin treatment were compared in this experiment, because KCl induced depolarization and glucose stimulated considerably common sets of genes, but the fold changes by KCl were more robust (Table 1). After 45 min of KCl or insulin treatment, the cells were harvested and analysis was performed in identical fashion to the previous experiments. As shown in Fig. 3, the genes were represented in log fold changes to allow a representation of the transcripts on a same scale for both up- and downregulation. The addition of LY294002 resulted in almost complete elimination of gene regulation by each stimulus. These results suggest that glucose, acting through depolarization, and insulin share common signaling pathways leading to gene expression and that this PI 3-kinase-dependent pathway seems to play a major role.
|
|
![]() |
DISCUSSION |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
A number of biologically relevant observations were made by the results of these experiments. 1) Over 90 potential candidates either up- or downregulated by glucose were uncovered in this study. Many of these genes are currently unknown ESTs and are not yet listed in UniGene or AllGenes (DoTS assemblies). 2) By comparing gene expression profiles after glucose- and KCl-induced depolarization, there was considerable overlap of genes activated by both agents, consistent with the hypothesis that glucose activation of rapidly responding genes is mediated predominantly via depolarization. These results are in accord with those of earlier studies showing that expression of several early response genes depended on depolarization (10,29). This premise was further confirmed by the results of the study of the measurements of gene expression in the presence of the depolarization inhibitor diazoxide in a subset of these genes. 3) The current results showed that exogenous insulin treatment at 100 nmol/l activated transcription of a number of newly described early response genes in insulinoma cells. The most notable finding, however, was that glucose and insulin treatments each activated a common subset of genes among several thousands on the arrays.
Validation of microarray gene expression results.
The validity of the observed microarray results ideally requires both technical validation and biological replication. Technical validation can be seen in that several clones corresponding to the same gene were present on our microarray set, for example, splicing factor arginine/serine rich (Sfrs5) and inhibitor of DNA binding 2 (Idb2) (Table 1). We also used qRT-PCR for validation of the microarray results and observed excellent correlation between microarray and qRT-PCR results testing the same samples (Table 3). Further, biological replication was achieved by two additional independent experiments in which expression of a subset of genes was measured by qRT-PCR (Table 3).
Limitations of the study.
Although a short 45-min treatment is suitable for evaluating initial transcriptional responses, the long-term changes that may result from these early changes remain to be determined. When MIN6 cells were incubated in high glucose for 4 h, preliminary results indicated that the majority of genes regulated at 45 min returned to baseline. This finding was consistent with characteristics of early response genes (8). On the other hand, new sets of genes not previously regulated at 45 min were now found to be regulated at 4 h (M.O., C.C.-M., Y.Z., E.B.-M., M.A.P., unpublished data). This may suggest cascades of gene expression in which genes responding early to stimuli may initiate or suppress expression of other genes.
The maximum concentrations of stimuli for glucose, KCl, and insulin were used. One may raise an issue regarding potential effect of osmotic stress as a glucose stimulus. In previous publications (10,29) and unpublished observations (E.B.-M. and M.A.P.), we have shown failure of nonmetabolizable glucose analogs to activate transcription, that the effects of glucose were blocked by diazoxide, and that KCl in the presence of a PI 3-kinase inhibitor showed no activation of genes (Fig. 3).
Another limitation of the current study was that the levels of expression of genes activated by glucose and insulin were for the most part in the 1.3- to 1.5-fold range. However, multiple determinations of fold change by microarray resulted in small variances of measurements, as seen in Fig. 1C, and several of the glucose and insulin response genes were validated by other means (Table 3). Although we believe that in general these 3050% changes in gene transcriptional rates are significant, rather than speculating on the importance of these modest fold changes in gene expression, we found that the main conclusion of the study is that common signaling pathways are similarly activated by two distinctive important ß-cell stimuli: glucose and insulin.
Not all ß-cell early response genes such as c-myc and c-fos were observed in the present study (7,34). These genes were not included in the original EPCon arrays, presumably because of the low abundance of those transcripts. Despite those omissions, our current studies are confirming and extending the previous findings, such as Egr1, and highlighting the use of microarray to assess similarities of treatments (e.g., glucose and insulin) and to identify signaling pathways activated by those treatments.
The microarray results with insulinoma cells are limited to ß-cells adapted to culture, yet isolated islets are a mixture of several cell types and are further complicated by this issue. Isolation of pure ß-cells from isolated islets raises additional concerns regarding the conditions of preparation, so every method has its limitations. The results of our experiments in insulinoma cells can now be used to measure responses to primary islets in culture, both rodent and human. Importantly, our microarray data now identify dozens of new early response genes, some known and some only described in the databases as ESTs.
Depolarization is the major component of glucose-regulated early transcriptional activation.
Remarkably, >90% of the effects of glucose on gene transcription were mimicked by KCl-induced depolarization (Fig. 2B). Although Ca2+ influx after glucose-induced depolarization is a critical mechanism resulting in insulin secretion, the role of depolarization in glucose-mediated gene expression has not been fully evaluated. Only few genes were shown to be regulated via depolarization in previous studies of early response genes, but the results of this study demonstrated that the majority of the glucose-regulated genes are activated through depolarization. Considering the results of previous studies performed on limited numbers of genes (10,29), this is likely to result in activation of Ca2+-regulated kinases leading to phosphorylation and activation of transcription factors. Although signaling events leading to transcriptional activation or suppression after depolarization remain to be further defined, the current studies indicated that the PI 3-kinase pathway seems to be involved (Fig. 3). It remains to be determined whether PI 3-kinase activation occurs via an autocrine/paracrine effect of insulin.
Activation of a common set of genes by glucose and insulin.
The observation that glucose and insulin induce an overlapping set of genes is interesting in light of recent data that provide evidence for an important role of insulin in the regulation of gene transcription in ß-cells. Insulin stimulation of ß-cells in vitro results in transcriptional induction of the preproinsulin, liver-type pyruvate kinase, and acetyl-CoA carboxylase I genes by activation of the insulin receptor substrate/PI 3-kinase pathway (3537). Mouse models (3840) clearly indicated the importance of the roles of the insulin receptor (ßIRKO mouse) and the downstream targets of the insulin-signaling pathway (IRS2) in islet ß-cell survival and function, but its precise mechanisms are still elusive.
Regarding the GO functions of the genes commonly regulated by glucose and insulin, there was a high proportion of regulated genes that belonged to catalytic and binding activity classes, and there were lower proportions in the other categories. The differences in the number of genes in these GO function categories differed slightly from the distribution of genes on the entire array (data not shown), although the array as a whole is not representative of the MIN-6 cells, because it comprises genes expressed in all cells of the pancreas. Further, it would be speculative to suggest that these categories of biological functions represent true changes at the cellular level considering the small numbers of genes in each GO category.
The broadly overlapping set of genes common to both glucose and insulin treatment in the current studies suggested at least three possible explanations. One is that these two islet ß-cell stimuli act through independent pathways that finally converge on the same gene promoters (Fig. 5A). Another is that the two stimuli act on two different pathways that converge on a common pathway that ends by activating the same gene promoters (Fig. 5B). This latter alternative seems the most likely because the PI 3-kinase inhibitor almost completely eliminated gene expression induced by both stimuli (Fig. 3). A third possibility is that the major effect of glucose/depolarization on ß-cell gene transcription is through an autocrine/paracrine effect of glucose-stimulated insulin secretion (Fig. 5C) (4143). It has been difficult to test the autocrine/paracrine effects of insulin, because MIN6 cells are immersed in a large amount of insulin and anti-insulin receptor antibody studies have been inconclusive (M.O., E.B.-M., M.A.P., unpublished data). To test the third hypothesis, experiments are currently underway to silence insulin receptor gene expression with small interfering RNA (44) in MIN6 cells.
|
![]() |
ACKNOWLEDGMENTS |
---|
We gratefully acknowledge the D. Melton lab (Harvard University) and K. Kaestner and C. Stoeckert labs (University of Pennsylvania), as well as Ellen Ostlund, Jessica Murray, Sandy Clifton, Hiroshi Inoue, Chris Sawyer, Mike Heinz, Elaine Mardis, and other members of the Genome Sequencing Center for their work with EPCon and microarrays. We would also like to thank Gary Stormo for helpful suggestions and Gary Skolnick for preparation of the manuscript.
![]() |
FOOTNOTES |
---|
Additional information for this article can be found in an online appendix at http://diabetes.diabetesjournals.org.
Address correspondence and reprint requests to M. Alan Permutt, MD, Division of Endocrinology, Metabolism,Lipid Research, Washington University School of Medicine, 660 S. Euclid Ave., Campus Box 8127, St. Louis, MO 63110. E-mail: apermutt{at}im.wustl.edu
Received for publication March 5, 2004 and accepted in revised form March 19, 2004
DMEM, Dulbeccos modified Eagles medium; EPCon, Endocrine Pancreas Consortium; EST, expressed sequence tag; FBS, fetal bovine serum; GO, gene ontology; IEG, immediate early gene; KATP channel, ATP-sensitive potassium channel; PI, phosphatidylinositol; qRT-PCR, quantitative RT-PCR; SSC, sodium chloride-sodium citrate
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|