|
|
||||||||
Department of Genomics, Wyeth Research, Cambridge, Massachusetts 02140
| ABSTRACT |
|---|
|
|
|---|
C2C12 cells; oligonucleotide array; functional category enrichment
| INTRODUCTION |
|---|
|
|
|---|
The myogenic program consists of two temporally separated processes: myoblast proliferation and differentiation. Proliferating mononucleate myoblasts expressing MyoD and Myf5 are committed to the muscle lineage and will continue to proliferate in the presence of mitogens under high-serum conditions in vitro. Upon serum deprivation, myoblasts activate transcription of myogenin and undergo irreversible cell cycle arrest following transcription of the cyclin-dependent kinase (Cdk) inhibitor, p21, and dephosphorylation of pRb (1, 40). Skeletal muscle differentiation then proceeds through the induction of muscle-specific gene expression and fusion of myoblasts into myotubes (1, 13, 14, 29, 41).
Here, we have used high-density oligonucleotide arrays to investigate transcriptional changes occurring during myoblast proliferation and differentiation in C2C12 cells, a well-characterized in vitro model of mouse skeletal muscle cell differentiation (1, 6). DNA microarray technology enables highly parallel analysis of mRNA expression patterns. The use of DNA microarrays to identify genes involved in skeletal muscle function and pathology has been limited to studies of mouse and monkey muscle aging and the effects of caloric restriction, as well as alveolar rhabdomyosarcomas (19, 20, 22, 42). We have identified sets of genes whose transcripts are differentially regulated between proliferating and differentiating C2C12 cells and have clustered these genes according to expression profiles. Systematic evaluation of the distribution of genes by biological function has revealed many clusters statistically enriched in functions relevant to skeletal muscle growth and differentiation.
| MATERIALS AND METHODS |
|---|
|
|
|---|
80% confluence), cells were rinsed briefly with PBS then grown in differentiation medium (DM; DMEM plus 1% bovine serum albumin) to initiate myogenic differentiation. Cells were harvested by lysis with Trizol Reagent (Invitrogen, Carlsbad, CA) 48 h and 72 h after plating in GM [GM and day 0 (d0) time points, respectively] and 24 h and 96 h after switching to DM [day 1 (d1) and day 4 (d4) time points, respectively]. Total RNA was isolated by Trizol extraction, and mRNA was isolated using the PolyATract mRNA Isolation System IV (Promega, Madison, WI).
Preparation of labeled targets and high-density oligonucleotide array hybridization.
Poly(A)+ mRNA samples were prepared independently from two total RNA samples at each time point. One microgram of each poly(A)+ mRNA sample was used to generate amplified, biotin-labeled cRNA. Ten micrograms of fragmented, biotin-labeled cRNA was mixed with 11 internal control labeled spike-in transcripts and hybridized in duplicate to oligonucleotide arrays containing
11,000 distinct murine genes and expressed sequence tags (ESTs). cRNA labeling, microarray hybridization, staining, and visualization was performed as described by Hill et al. (15).
Data analysis.
Absolute decision calls ("present," "absent," or "marginal") for each gene on the arrays were determined by the GeneChip Software (Affymetrix, Santa Clara, CA). Transcript levels, indicated as gene frequency, were quantified as described by Hill et al. (15). There were 5,303 genes we called "present," with a frequency higher than the sensitivity of detection in a minimum of 4 of 16 arrays (4 replicates at 4 time points), and these were selected for further data analysis.
The experimental design has both a treatment variable and a time variable, where treatment refers to the two conditions in which the cells were grown: GM or DM. Cells grown under each treatment condition were harvested at two time points (treatment intervals). Frequency values for 5,303 genes were subjected to one-way nested analysis of variance (ANOVA) tests. ANOVA P values calculated for each gene were adjusted with the Bonferroni method of multiple comparison to control for type I error rate (33). Genes with Bonferroni-adjusted significance levels less than 1% (n = 644) were considered to have statistically significant changes in gene expression between treatments and/or between treatment intervals. Mean frequencies for each gene at each time point were calculated from the individual frequency values over four replicates. For low-frequency genes (maximum mean frequencies less than 40 for all time points), additional criteria were imposed to increase confidence that the estimated changes reflect real differences in expression levels. To qualify, low-abundance mRNAs had to exhibit an arithmetic difference greater than 10 and at least a 2.5-fold difference between the maximum and minimum mean frequencies for all time points. Fifteen genes with low mean frequencies did not pass these additional criteria. For cluster analysis, mean expression frequencies of 629 genes with differential expression were scaled so that each gene had a mean value of zero and a variance of one, and the scaled expression profiles were clustered using a self-organizing map (SOM) algorithm (36).
Enrichment of functional categories within expression clusters.
BLAST searches were performed with EST sequences against GenBank and Affymetrix Hu6800 gene sequences to identify gene hits and provide annotation. We classified 544 distinct genes into functional categories using a hierarchical classification scheme modified from Cho et al. (8) with 167 categories. Assignments were made to one or more categories using gene annotation from GeneExpress 2000 (Gene Logic, Gaithersburg, MD) and publicly available databases. Of 433 distinct annotated genes, 382 were assigned to functional categories, with 193 being assigned more than one category. We designated 111 sequences as ESTs. Since the classification scheme is hierarchical, with several categories having multiple subcategories, fold enrichment and P values were calculated for each level of the hierarchy independently. For each category in each cluster, a fold enrichment value was calculated based on the frequency of the category in the cluster compared with the frequency of the category in the entire set of 544 transcripts. P values were calculated using the hypergeometric probability distribution (38) and represent the probability of observing fold enrichment of the category among a random sampling of genes.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
Cell adhesion and extracellular matrix genes.
Cluster 1,3 is significantly enriched in genes with extracellular matrix (ECM) and cell adhesion functions, and cluster 1,4 is significantly enriched in cell adhesion genes (Table 2). In both clusters, genes are often classified as having both ECM and cell adhesion functions.
Multiple components of the myocyte basal lamina [type IV collagens (Col4a1, Col4a2), nidogen 2 (Nid2), and ß -dystroglycan (Dag1)] have maximal transcript levels 24 h after serum withdrawal (Fig. 3A). Five additional collagen transcripts (Col6a1, Col6a2, Col8a1, Col3a1, Col5a2) and other ECM protein-encoding transcripts, including matrix metalloprotease 2, matrix-associated glycoprotein 2 (Magp2), matrillin 2, syndecan 2, and matrix Gla protein (Mglap), are also upregulated during early myoblast differentiation (Table 2). Several genes with cell adhesion functions [decorin (Dcn), Cd81, tetraspanin-3 (Tm4sf8), Vcam1, Ncam, and M-cadherin (Cdh15)] are expressed at highest levels in terminally differentiated myocytes (Fig. 3).
|
Genes involved in cellular metabolism.
Cluster 2,1 is enriched in genes with roles in carbohydrate metabolism (Tables 1 and 2), and other key glycolytic enzymes (Pfka and Pygm) are observed in the related cluster 1,1. Cluster 1,1 is enriched in amino acid transport genes, Slc1a6 (Eeat4) and Mdu1 (4F2hc), which are members of the transport systems XAG- and L,y+, respectively (Tables 1 and 2).
Cluster 3,1 is significantly enriched in transcripts encoding mitochondrial transport proteins (Timm9, Timm13a, Timm23, Slc25a5), which are coregulated with genes required for oxidative respiration [cytochrome c (Cycs) and ATP synthase (Atp5g1)].
Genes involved in cell cycle regulation, DNA replication, mRNA transcription, and chromatin modification.
Clusters in rows 3 and 4 of the cluster matrix (Fig. 2) are enriched in genes with functions necessary for actively proliferating cells: mRNA transcription, DNA replication, and cell cycle control and mitosis. Cluster 3,2 is significantly enriched in genes functioning in cell cycle control and mitosis as well as the subcategory chromosome segregation. In addition, 20 cell cycle control genes are present in clusters 3,3 and 4,3 (Tables 1 and 2; and Supplemental Data1
published online at the Physiological Genomics web site). Such coordinately regulated genes include cyclins (Ccna2, Ccnb1-rs1), protein kinases (Cdc2a, Cdk4, Cks1, Cks2), protein phosphatases (Cdc25c, Cdkn3, Ppp1r7), and chromosome segregation genes (e.g., Bub3, Ran, Rangap1, Ranbp1) (Fig. 4A). The cyclin-dependent kinase inhibitors, p21Cip1 (Cdkn1a; cluster 1,1) and p57Kip2 (Cdkn1c; cluster 1,4), as well as other genes functioning in negative growth control (Bin1, Gadd45a, Gadd45b, and Gadd45g), are appropriately expressed in an opposite pattern to cell cycle promoting genes (Figs. 2 and 4B).
|
DNA replication genes are enriched in cluster 4,3, and additional replication genes are found in rows 3 and 4 of the cluster matrix, including minichromosome-deficient family members, ribonuclease reductase M1 genes, and proliferating cell nuclear antigen A (Pcna) (Tables 1 and 2; Fig. 4A; and Supplemental Data).
Other proliferation and differentiation genes.
A number of transcription factors and signaling molecules with regulatory roles in cellular proliferation and differentiation were differentially regulated during C2C12 myogenesis (Fig. 5). Transcription factors increasing during muscle differentiation include MyoD1, Six1, Stat3, Hes6, Klf4, and Osf2 (Fig. 5A). Members of the Id gene family of differentiation inhibitors, Idb1, Idb2, and Idb3, were all expressed at higher levels in proliferating myoblasts and decrease in expression in differentiating cells (Fig. 5B). Expression of the secreted growth and differentiation factor, Igf2, is increased greater than 10-fold upon differentiation in C2C12 cells, as are transcripts for several IGF binding proteins (Igfbp2, Igfbp5, and Nov) (Fig. 5C). Members of the TGFß gene family of secreted signaling molecules, activin (Inhbb) and TGFß3 (Tgfb3) as well as the bone morphogenetic protein (BMP) inhibitor, follistatin (Fst), increase in expression at day 0 and peak in differentiating muscle cells (Fig. 5C). Other signaling pathways important in cellular differentiation, the Notch and Wnt pathways, have members [Notch3 and Sfrp-2 (Sdf5), respectively] that are expressed at highest levels in differentiated C2C12 cells.
|
| DISCUSSION |
|---|
|
|
|---|
Clustering and functional classification of differentially regulated genes and subsequent statistical analysis of the distribution of functional categories identified both expected and unexpected results. Enrichment for genes involved in muscle contraction, cytoskeletal organization, cell adhesion, cell-matrix interactions, cellular metabolism, cell cycle regulation, DNA replication, and mRNA transcription were anticipated.
Transcripts encoding muscle contractile apparatus proteins, proteins regulating muscle contraction, cell-cell interactions, and cell-matrix interactions increase in expression during the differentiation process (Fig. 2; Table 2). Cell-matrix and cell-cell interactions mediate terminal differentiation of myoblasts including cell fusion and assembly of the muscle contractile apparatus (26, 27). Thus, besides upregulating the expression of genes responsible for muscle contraction, differentiating myocytes promote muscle integrity and function by expressing genes encoding structural and regulatory components of the myocyte extracellular matrix and plasma membrane.
Muscle cells utilize glucose, fatty acid, ketone bodies, and amino acids as fuel sources and store glucose mainly as glycogen. In addition, muscle is the major storage tissue of amino acids. Our study identified several genes that are upregulated in differentiating cells that function in carbohydrate metabolism (Table 2). Elevated expression of amino acid transporters in differentiated myotubes (Table 2) may be due to an increase in amino acid storage requirements or may be a cellular response to low nutrient levels induced by serum starvation.
Cell cycle progression requires ATP for many cellular processes including protein synthesis, degradation, and cytoskeletal rearrangements. Most cellular ATP is synthesized through mitochondrial oxidative phosphorylation. We found that transcripts encoding proteins required for oxidative respiration were coregulated with transcripts encoding mitochondrial transport proteins (Table 2; Fig. 2). In addition to metabolic requirements, actively proliferating cells require DNA replication, mRNA transcription, and cell division. Several clusters with genes that are downregulated upon the initiation of myogenic differentiation are enriched in such functional categories (Table 2).
Surprisingly, our study revealed several clusters that are enriched in genes with immune functions (Table 2). Of all genes with immunity roles, the majority function during the inflammatory process. Inflammation in response to muscle injury or disease is intimately associated with muscle regeneration (18). In muscle tissue, inflammation involves chemokine and cytokine production by damaged fibers, leukocyte and macrophage migration, and activation of muscle satellite cell proliferation and differentiation, followed by phagocytic apoptosis (4, 18, 31). Upregulation of chemokine and cytokine mRNA expression during early or late differentiation may act as a trigger to promote muscle regeneration through leukocyte and macrophage recruitment, bone-marrow-derived stem cell recruitment (10), or may be part of a cellular response to serum deprivation.
Many of the genes that we have identified as differentially regulated transcripts by DNA microarray expression profiling have been characterized previously in C2C12 cells or other skeletal muscle cell lines by mRNA transcription methods and/or by functional analysis. Those characterized by both expression and functional analyses by others include Igf2 (reviewed in Ref. 11), p21Cip1 (Cdkn1a) (14), Bin1 (25), caveolin 3 (Cav3) (12, 37), M-cadherin (Cdh15) (21, 43), FasL (31), decorin (Dcn) (30), Idb1 (17), Idb3 (3, 7), and DNA methyltransferase (Dnmt1) (24, 35). In all of these cases, the RNA expression profile we have detected agrees with published reports.
In addition to identifying genes with previously characterized functions in muscle cell differentiation, our results highlight genes for which a potential role in skeletal muscle differentiation has not been defined. Several of these genes function in osteoblasts: the transcription factor, Osf2 (cluster 1,2) (9) and the osteoblast matrix proteins, osteomodulin/osteoadherin (Omd; cluster 1,1) (34) and osteoglycin (Ogn; cluster 1,3) (32). Osf2 (Cbfa1) mRNA expression is increased by BMP treatment of C2C12 cells, which causes a transformation from myogenic to osteoblast differentiation (23). The induction of Osf2 transcripts in differentiating C2C12 myoblasts that was observed in our study has not been observed by others. This discrepancy may be due to different culturing conditions, specifically the amount of serum in the medium. The presence of osteoblast extracellular matrix components in differentiating muscle cells may be indicative of the pluripotent state of C2C12 cells, which have the capacity to differentiate into osteoblasts when challenged with the appropriate cues. Alternatively, these extracellular matrix proteins may perform a generalized ECM function in differentiated cells.
A comparison of the genes identified in our study to those that have been functionally characterized in muscle cell lines and other expression profiling studies can provide clues as to the roles these genes are playing during the C2C12 differentiation process. Two relevant large-scale expression profiling studies examining the serum response in human fibroblasts (16) and the human fibroblast cell cycle (8) have identified genes that are also regulated in our study. In addition to the induction of cell cycle and proliferation genes [e.g., p57Kip2 (Cdkn1c), Pcna, Rrm1, Rrm2] and immediate-early transcription factors (Id2, Id3), Iyer and colleagues (16) observed a striking upregulation of genes involved in the wound healing process during the serum response. These genes are involved in a variety of cellular processes: inflammation/angiogenesis [Scya2 (MCP1)], tissue remodeling (Plaur, Fmod, Col3a1), signal transduction (Sgk) and chromosome segregation (Ranbp1). Fibroblast cell cycle regulated genes that are also differentially regulated during C2C12 differentiation include those involved in cell cycle control (e.g., Ccna2, Ccnb1, Cdc2, Cdkn1a), proliferation and differentiation (Notch3, Inhbb, Gadd45a, Gro1), and cell-cell adhesion/ extracellular matrix genes (Matn2, Nid2, Plaur, Hmmr, Tm4sf1) (8). Thus the identification of common genes regulated by serum addition and/or by cell cycle progression suggests that several of the transcriptional changes we have identified in the C2C12 time course may be due to mitotic synchronization, a physiological response to serum withdrawal (8).
Importantly, our study revealed 111 distinct ESTs with no appreciable similarity to known genes. Although ESTs were not significantly enriched as a functional category in any gene cluster, it is interesting to note that 50% of all ESTs are coordinately regulated in clusters 4,1 through 4,4, representing
30% of the unique genes in each of the four clusters (Fig. 2; and Supplemental Data). Our functional category enrichment data can be used to predict gene function and would suggest that many of these uncharacterized genes may function in mRNA transcription, chromatin modification, DNA replication, and cell cycle control and mitosis.
This study provides a useful database of genes that are differentially expressed throughout myogenic proliferation and differentiation. As we have reported functional categories that are statistically enriched in expression clusters and have highlighted only a portion of genes in this study, a more detailed analysis of all the data provided (see Supplemental Data) will allow others to investigate particular genes, pathways, or cellular processes of interest. Functional analysis of genes identified in this study as well as characterization of ESTs in in vitro and in vivo models, will likely lead to the identification of genes required for cell cycle progression and withdrawal and skeletal muscle differentiation.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
Address for reprint requests and other correspondence: C. P. Miller, Dept. of Genomics, Wyeth Research, 35 Cambridge Park Drive, Cambridge, MA 02140 (E-mail: cmiller{at}wyeth.com).
10.1152/physiolgenomics.00011.2002.
1 Supplementary materials (Supplemental Tables 1 and 2 and Supplemental Fig. 1) to this article are available online at http://physiolgenomics.physiology.org/cgi/content/full/10/2/103/DC1. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
S. J. Ho Sui, D. L. Fulton, D. J. Arenillas, A. T. Kwon, and W. W. Wasserman oPOSSUM: integrated tools for analysis of regulatory motif over-representation Nucleic Acids Res., July 13, 2007; 35(suppl_2): W245 - W252. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Seidel, N. Rozwadowska, K. Tomczak, and M. Kurpisz Myoblast preparation for transplantation into injured myocardium Eur. Heart J. Suppl., December 1, 2006; 8(suppl_H): H8 - H15. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Nedachi and M. Kanzaki Regulation of glucose transporters by insulin and extracellular glucose in C2C12 myotubes Am J Physiol Endocrinol Metab, October 1, 2006; 291(4): E817 - E828. [Abstract] [Full Text] [PDF] |
||||
![]() |
I.-H. Park and J. Chen Mammalian Target of Rapamycin (mTOR) Signaling Is Required for a Late-stage Fusion Process during Skeletal Myotube Maturation J. Biol. Chem., September 9, 2005; 280(36): 32009 - 32017. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Y. Kostrominova, D. E. Dow, R. G. Dennis, R. A. Miller, and J. A. Faulkner Comparison of gene expression of 2-mo denervated, 2-mo stimulated-denervated, and control rat skeletal muscles Physiol Genomics, July 14, 2005; 22(2): 227 - 243. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Blais, M. Tsikitis, D. Acosta-Alvear, R. Sharan, Y. Kluger, and B. D. Dynlacht An initial blueprint for myogenic differentiation Genes & Dev., March 1, 2005; 19(5): 553 - 569. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Paris, C. Virtanen, Z. Lu, and M. Takahashi Identification of MEF2-regulated genes during muscle differentiation Physiol Genomics, December 15, 2004; 20(1): 143 - 151. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. C. Goetsch, T. J. Hawke, T. D. Gallardo, J. A. Richardson, and D. J. Garry Transcriptional profiling and regulation of the extracellular matrix during muscle regeneration Physiol Genomics, August 15, 2003; 14(3): 261 - 271. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |