Dynamic resolution of functionally related gene sets in response to acute heat stress
© Szustakowski et al; licensee BioMed Central Ltd. 2007
Received: 07 August 2006
Accepted: 05 June 2007
Published: 05 June 2007
Using a gene clustering strategy we determined intracellular pathway relationships within skeletal myotubes in response to an acute heat stress stimuli. Following heat shock, the transcriptome was analyzed by microarray in a temporal fashion to characterize the dynamic relationship of signaling pathways.
Bioinformatics analyses exposed coordination of functionally-related gene sets, depicting mechanism-based responses to heat shock. Protein turnover-related pathways were significantly affected including protein folding, pre-mRNA processing, mRNA splicing, proteolysis and proteasome-related pathways. Many responses were transient, tending to normalize within 24 hours.
In summary, we show that the transcriptional response to acute cell stress is largely transient and proteosome-centric.
Damaged skeletal muscle tissue retains the capacity to self-repair through activation, expansion and fusion of resident satellite, progenitor cells [1–3]. However, the signaling mechanisms required to trigger satellite cell activation remain unclear in this context. Emerging pathways involved in the injury-response profile of skeletal muscle include Notch , Akt/Foxo  and NFκB [6, 7], which appear to play distinct roles in regeneration, hypertrophy and atrophy. Although these signaling pathways are well established, there remains a dearth of molecular detail, including the transcriptional consequence of atrophic stimuli, and this aspect forms the focus of the current study.
We hypothesized that subtle variance in gene expression may underlie significant functional events within the cell. Consistently, previous studies have addressed the transcriptional profile of recovering muscle tissue following acute injury . However, such studies rely on the analysis of individual genes, which is limiting from both statistical and biological perspectives. Here, we performed a study with the tenet that mechanistic data may be more readily revealed in the context of pathway-related gene clustering [9, 10]. We employed Gene Set Enrichment Analysis (GSEA), a bioinformatics approach conceived to exploit functional regulation at the pathway level. The study focused on the dynamic expression profile of acutely stressed skeletal muscle cells, and the resultant data revealed major transcriptional adjustments in a wide range of functional categories. These findings may underscore a coordinated pathway-centric response, and reveal rational strategies for target discovery in atrophic muscle diseases.
Atrogin-1 is expressed specifically in skeletal and cardiac muscle, and has been shown to regulate muscle atrophy [19, 20]. In response to heat shock, Atrogin-1 (Fbxo32) was upregulated after 8 hours, and this level of expression was maintained for at least 15 hours (Figure 1D). This dynamic profile is consistent with previous reports of muscle injury response, whereby the expression of atrogin-1 was concomitant with atrophy and recovery . Atrogin-1 is regulated by Foxo transcription factors, via nuclear localization of Foxo in response to various stimuli that inhibit PI3K/Akt signaling and subsequently drive atrophy [5, 22].
Transient profiles of gene expression
Significantly upregulated genes at specific time points following heat shock
Hsp40 homolog – gi:12839599
Cystein rich protein 61
FBJ osteosarcoma oncogene
Bcl2-associated athanogene 3
Developmentally and sexually retarded with transient immune abnormalities
Procollagen, type VI, alpha 3
Placental growth factor
Procollagen, type V, alpha 1
Solute carrier family 19 (thiamine transproter), member 2
Procollagen, type I, alpha 1
AXIN1 up-regulated 1
Actin, beta, cytoplasmic
Moderate similarity to protein ref:NP_004410.2 (H.sapiens) phosphatase
Ribosomal protein S24
Heme oxygenase (decycling) 1
Low density lipoprotein receptor-related protein 1
Growth arrest and DNA-damage-inducible 45 gamma
Perlecan (heparan sulfate proteoglycan 2)
Myeloid differentiation primary response gene 116
Myosin heavy chain IX
v-maf musculoaponeurotic fibrosarcoma oncogene family, protein F (avian)
Early growth response 1
Dual specificity phosphatase 2
Expressed sequence – gi:2918507
Diphtheria toxin receptor
EST – gi:8772531
Transformed mouse 3T3 cell double minute 2
Interferon-related developmental regulator 1
Activity regulated cytoskeletal-associated protein
mRNA – gi:16877863
DnaJ (Hsp40) homolog, subfamily A, member 4
Acidic (leucine-rich) nuclear phosphoprotein 32 family, member A
Lysosomal acid lipase 1
Lysyl oxidase-like 2
ATPase, Na+/K+ transporting, alpha 1 polypeptide
Guanine nucleotide binding protein, alpha inhibiting 2
Monocyte to macrophage differentiation-associated
Protein tyrosine phosphatase 4a2
Purinergic receptor P2X, ligand-gated ion channel, 7
Histone 1, H4i
Similar to hypothetical protein MGC4368 – gi:15655075
Aldehyde dehydrogenase family 3, subfamily A1
Expressed sequence – gi:2519775
Pyridoxal (pyridoxine, vitamin B6) kinase
cDNA – gi:9515377
CDP-diacylglycerol synthase (phosphatidate cytidylyltransferase) 2
Protein kinase C-like 1
Related to uridine kinase – gi:12844097
mRNA – gi:4779190
Gene set enrichment analysis
Most significantly upregulated gene sets at 8 hours post heat shock
No. of probes in set
Protein metabolism and modification
mRNA transcription regulation
Protein complex assembly
Discussion and Conclusion
This study was conducted to explore changes in gene expression following acute injury of C2C12 myotubes, and to reveal regulation of pathways during initial recovery. 866 individual transcripts represented on the microarray varied more than 2 fold over the course of the study, and the data revealed significant transcriptional changes in a wide variety of functional classes following heat shock damage to the cells. We devised the current study with two distinct dimensions. Firstly, to trace the time-course of recovery over a 24 hour period, and secondly, to analyze in a pathway-centric manner in order to extract hypotheses of coordinated molecular mechanisms that occur during cellular stress and recovery. Two notable pathways significantly modulated during recovery were the BMP signaling pathway, and the extracellular matrix-protein mediated signaling pathway, both of which are known to function in myogenesis, either by directing lineage specification (BMP) or structural/regenerative signaling (ECM-M). Additionally, the muscle-specific E3 ligase gene, Atrogin-1 was upregulated significantly, consistent with its reported role in mediating muscle atrophy in response to catabolic signals , and this in turn was coincident with our observation of nuclear translocation of Foxo3a following heat shock [5, 22]. This reconstruction of a well documented atrophic pathway served to validate the heat shock treatment paradigm as a useful model of acute muscle stress.
The experimental approach carried out in this study demonstrated the value of characterizing cellular responses in a temporal fashion, as transcriptional regulation was commonly found to be non-static between contiguous time points within the first 24 hours of treatment, and often presented transient regulation that required a time-course dimension to clarify.
A notable trend that emerged from the data was an apparent increase in protein turnover apparatus-related gene sets. The most significantly upregulated gene sets included those for both transcriptional activity, translational activity and protein modification, as well as the proteolysis and proteosome pathways. This suggests that a primary response of skeletal myotubes to heat shock treatment is to renew the protein complement of the cell, presumably as a defense mechanism against proteomic damage, and potential aberrant consequences in the function or integrity of the cell. The data also show that much of this response is transient, and appears to normalize within 24 hours of the acute insult at the transcriptional level.
In summary, the data presented in this report demonstrate a dynamic resolution of gene expression in response to acute stress that implicates waves of pathway activity in the cell's metabolic recovery. Further validation of these data are required, however, in order to demonstrate clearly the practical nature of these pathways as part of the functional stress response profile.
Cell culture and treatment
C2C12 mouse skeletal myoblasts were cultured in DMEM high glucose with 10% FBS and 1% penicillin/streptomycin (Gibco), and maintained at 37°C and 5% CO2. Cells were grown to confluence and then differentiated for seven days in media containing 3% fetal bovine serum. Cells were heat shocked at 47°C for 30 minutes, and returned to the incubator at 37°C. After heat shock, cells were harvested at 1, 2, 4, 8, and 24 hours and processed to isolate the RNA using the RNeasy kit (Qiagen). Myotubes were isolated for expression profiling by serially exposing the culture to increasing concentrations of trypsin. This allowed for the multinucleated myotubes to become detached from the plate, while the mononucleate myoblasts remained adherent. This process ensured that we were profiling differentiated skeletal myotubes only, without contamination of myoblasts. 20 μg RNA per sample was microarrayed using the Mouse Affymetrix MOE430 Plus 2.0 chip. The heat-shock treatment conditions were determined not to be cytotoxic, as continued culture in normal growth conditions was sufficient for full recovery of the cells (data not shown).
Raw intensity files were imported into GeneSpring® v.6.2 (Silicon Genetics, US) and filtered such that genes with a raw expression level of less than 100 were excluded from any future analyses. The resulting gene list was filtered to identify genes with a fold change increase or decrease of 2 or greater. Statistically significant genes were found using one-way ANOVA with time as variable and no multiple testing correction.
where μi,jis the average expression value for gene i under condition j . Genes were then sorted according to their expression ratios. Next, the collection of available gene sets were projected onto the sorted gene list. In essence, this step applies a priori biological knowledge to the experimental data to identify functionally related genes that are expressed in a coordinated fashion. Gene sets were processed individually and used to partition the genes into two groups: those gene in a pathway, and all other genes measured on the Genechip. A two-tailed Wilcoxon rank-sum test was then calculated to determine if the genes in each gene set were enriched at either the top or bottom of the sorted list. The Wilcoxon p-values were then transformed into false discovery rate q-values [24, 25] as a means of multiple hypothesis testing correction.
We thank Qicheng Ma for help in the development of GSEA.
- Charge SB, Rudnicki MA: Cellular and molecular regulation of muscle regeneration. Physiol Rev 2004, 84: 209-238. 10.1152/physrev.00019.2003View ArticlePubMedGoogle Scholar
- Sherwood RI, Christensen JL, Conboy IM, Conboy MJ, Rando TA, Weissman IL, Wagers AJ: Isolation of adult mouse myogenic progenitors: functional heterogeneity of cells within and engrafting skeletal muscle. Cell 2004, 119: 543-554. 10.1016/j.cell.2004.10.021View ArticlePubMedGoogle Scholar
- Collins CA, Olsen I, Zammit PS, Heslop L, Petrie A, Partridge TA, Morgan JE: Stem cell function, self-renewal, and behavioral heterogeneity of cells from the adult muscle satellite cell niche. Cell 2005, 122: 289-301. 10.1016/j.cell.2005.05.010View ArticlePubMedGoogle Scholar
- Conboy IM, Rando TA: Aging, stem cells and tissue regeneration: lessons from muscle. Cell Cycle 2005, 4: 407-410.View ArticlePubMedGoogle Scholar
- Stitt TN, Drujan D, Clarke BA, Panaro F, Timofeyva Y, Kline WO, Gonzalez M, Yancopoulos GD, Glass DJ: The IGF-1/PI3K/Akt pathway prevents expression of muscle atrophy-induced ubiquitin ligases by inhibiting FOXO transcription factors. Mol Cell 2004, 14: 395-403. 10.1016/S1097-2765(04)00211-4View ArticlePubMedGoogle Scholar
- Hunter RB, Kandarian SC: Disruption of either the Nfkb1 or the Bcl3 gene inhibits skeletal muscle atrophy. J Clin Invest 2004, 114: 1504-1511. 10.1172/JCI200421696PubMed CentralView ArticlePubMedGoogle Scholar
- Cai D, Frantz JD, Tawa NE Jr., Melendez PA, Oh BC, Lidov HG, Hasselgren PO, Frontera WR, Lee J, Glass DJ, Shoelson SE: IKKbeta/NF-kappaB activation causes severe muscle wasting in mice. Cell 2004, 119: 285-298. 10.1016/j.cell.2004.09.027View ArticlePubMedGoogle Scholar
- Zhao P, Hoffman EP: Embryonic myogenesis pathways in muscle regeneration. Dev Dyn 2004, 229: 380-392. 10.1002/dvdy.10457View ArticlePubMedGoogle Scholar
- Patti ME, Butte AJ, Crunkhorn S, Cusi K, Berria R, Kashyap S, Miyazaki Y, Kohane I, Costello M, Saccone R, Landaker EJ, Goldfine AB, Mun E, DeFronzo R, Finlayson J, Kahn CR, Mandarino LJ: Coordinated reduction of genes of oxidative metabolism in humans with insulin resistance and diabetes: Potential role of PGC1 and NRF1. Proc Natl Acad Sci U S A 2003, 100: 8466-8471. 10.1073/pnas.1032913100PubMed CentralView ArticlePubMedGoogle Scholar
- Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, Puigserver P, Carlsson E, Ridderstrale M, Laurila E, Houstis N, Daly MJ, Patterson N, Mesirov JP, Golub TR, Tamayo P, Spiegelman B, Lander ES, Hirschhorn JN, Altshuler D, Groop LC: PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet 2003, 34: 267-273. 10.1038/ng1180View ArticlePubMedGoogle Scholar
- Mason MR, Lieberman AR, Anderson PN: Corticospinal neurons up-regulate a range of growth-associated genes following intracortical, but not spinal, axotomy. Eur J Neurosci 2003, 18: 789-802. 10.1046/j.1460-9568.2003.02809.xView ArticlePubMedGoogle Scholar
- Tsujino H, Kondo E, Fukuoka T, Dai Y, Tokunaga A, Miki K, Yonenobu K, Ochi T, Noguchi K: Activating transcription factor 3 (ATF3) induction by axotomy in sensory and motoneurons: A novel neuronal marker of nerve injury. Mol Cell Neurosci 2000, 15: 170-182. 10.1006/mcne.1999.0814View ArticlePubMedGoogle Scholar
- Nakagomi S, Suzuki Y, Namikawa K, Kiryu-Seo S, Kiyama H: Expression of the activating transcription factor 3 prevents c-Jun N-terminal kinase-induced neuronal death by promoting heat shock protein 27 expression and Akt activation. J Neurosci 2003, 23: 5187-5196.PubMedGoogle Scholar
- Sehl PD, Tai JT, Hillan KJ, Brown LA, Goddard A, Yang R, Jin H, Lowe DG: Application of cDNA microarrays in determining molecular phenotype in cardiac growth, development, and response to injury. Circulation 2000, 101: 1990-1999.View ArticlePubMedGoogle Scholar
- Muchowski PJ, Valdez MM, Clark JI: AlphaB-crystallin selectively targets intermediate filament proteins during thermal stress. Invest Ophthalmol Vis Sci 1999, 40: 951-958.PubMedGoogle Scholar
- Vilaboa NE, Garcia-Bermejo L, Perez C, de Blas E, Calle C, Aller P: Heat-shock and cadmium chloride increase the vimentin mRNA and protein levels in U-937 human promonocytic cells. J Cell Sci 1997, 110 ( Pt 2): 201-207.Google Scholar
- Oblander SA, Zhou Z, Galvez BG, Starcher B, Shannon JM, Durbeej M, Arroyo AG, Tryggvason K, Apte SS: Distinctive functions of membrane type 1 matrix-metalloprotease (MT1-MMP or MMP-14) in lung and submandibular gland development are independent of its role in pro-MMP-2 activation. Dev Biol 2005, 277: 255-269. 10.1016/j.ydbio.2004.09.033View ArticlePubMedGoogle Scholar
- Zahradka P, Harding G, Litchie B, Thomas S, Werner JP, Wilson DP, Yurkova N: Activation of MMP-2 in response to vascular injury is mediated by phosphatidylinositol 3-kinase-dependent expression of MT1-MMP. Am J Physiol Heart Circ Physiol 2004, 287: H2861-70. 10.1152/ajpheart.00230.2004View ArticlePubMedGoogle Scholar
- Bodine SC, Latres E, Baumhueter S, Lai VK, Nunez L, Clarke BA, Poueymirou WT, Panaro FJ, Na E, Dharmarajan K, Pan ZQ, Valenzuela DM, DeChiara TM, Stitt TN, Yancopoulos GD, Glass DJ: Identification of ubiquitin ligases required for skeletal muscle atrophy. Science 2001, 294: 1704-1708. 10.1126/science.1065874View ArticlePubMedGoogle Scholar
- Gomes MD, Lecker SH, Jagoe RT, Navon A, Goldberg AL: Atrogin-1, a muscle-specific F-box protein highly expressed during muscle atrophy. Proc Natl Acad Sci U S A 2001, 98: 14440-14445. 10.1073/pnas.251541198PubMed CentralView ArticlePubMedGoogle Scholar
- Sandri M, Sandri C, Gilbert A, Skurk C, Calabria E, Picard A, Walsh K, Schiaffino S, Lecker SH, Goldberg AL: Foxo transcription factors induce the atrophy-related ubiquitin ligase atrogin-1 and cause skeletal muscle atrophy. Cell 2004, 117: 399-412. 10.1016/S0092-8674(04)00400-3PubMed CentralView ArticlePubMedGoogle Scholar
- Kamei Y, Miura S, Suzuki M, Kai Y, Mizukami J, Taniguchi T, Mochida K, Hata T, Matsuda J, Aburatani H, Nishino I, Ezaki O: Skeletal muscle FOXO1 (FKHR) transgenic mice have less skeletal muscle mass, down-regulated Type I (slow twitch/red muscle) fiber genes, and impaired glycemic control. J Biol Chem 2004, 279: 41114-41123. 10.1074/jbc.M400674200View ArticlePubMedGoogle Scholar
- Katagiri T, Imada M, Yanai T, Suda T, Takahashi N, Kamijo R: Identification of a BMP-responsive element in Id1, the gene for inhibition of myogenesis. Genes Cells 2002, 7: 949-960. 10.1046/j.1365-2443.2002.00573.xView ArticlePubMedGoogle Scholar
- Storey JD, Tibshirani R: Statistical methods for identifying differentially expressed genes in DNA microarrays. Methods Mol Biol 2003, 224: 149-157.PubMedGoogle Scholar
- Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci U S A 2003, 100: 9440-9445. 10.1073/pnas.1530509100PubMed CentralView ArticlePubMedGoogle Scholar