- Research article
- Open Access
Dynamic resolution of functionally related gene sets in response to acute heat stress
BMC Molecular Biology volume 8, Article number: 46 (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.
In the present study, we exposed C2C12 skeletal myotubes to heat shock-induced damage in order to identify coordinate patterns of transcriptional regulation. To initially authenticate this cell-based injury model, we identified several stress-related biomarkers and regulatory patterns to confirm bona fide damage response in the cells. Following heat shock treatment, cells recovered in growth medium at 37°C, and gene expression profiles were generated at 0, 1, 2, 4, 8 and 24 hours post-treatment. ATF3 and c-jun are critically regulated in the neuronal response to damage and during degeneration. Furthermore, these genes may play an anti-apoptotic role in the survival response of neuronal cells [11–13]. Here, both ATF3 and c-jun transcript levels were upregulated in response to heat-shock, consistent with their reported coordinate role in cell damage response (Figure 1A). A transient peak in expression of both genes occurred at 8 hours and attenuated at 24 h post-damage induction, indicating a regulated damage-induced response. Vimentin, an intermediate filament protein that plays a role in wound healing  and response to thermal stress [15, 16], was augmented immediately following heat shock, and peaked at 2 hours. This rapid response to cell stress was consistent with the role of vimentin in protection against cytotoxic effects (Figure 1B). Also, the membrane type 1-matrix metalloproteases MMP14 and MMP2, which are functional partners during skeletal development  and stress response in the context of vascular injury , were robustly co-regulated at the transcriptional level in response to heat shock of C2C12 skeletal myotubes. Their expression level rose within 2–4 hours of insult, followed by a repression of expression levels at 8 and 24 hours (Figure 1C).
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
Of the approximately 46,000 probes arrayed on the gene chips, the most variable and dynamically regulated were characterized by basic clustering methodology (Figure 2). The most actively regulated group of genes immediately following treatment (ie. 1 hour post heat shock) unsurprisingly comprised a number of heat-shock proteins. hsp1a, hsp1b and hsp8 were among the most significantly upregulated transcripts at this time point (Figure 2A). Gene sets that clustered at 2, 4, 8 and 24 hours following heat shock are displayed in Figure 2 and are partially represented in Table 1. The 8 hour gene set was the largest group that displayed a significant upregulation, and is likely comprised of genes that were regulated downstream of immediate early genes. Therefore, this cluster may represent an important set of genes that functionally dictate the coordinated response of myogenic cells to stress or injury.
Gene set enrichment analysis
We employed a the gene set enrichment analysis strategy  to cluster functionally related genes in order to consider shifts in pathway activity. Individual gene sets were mined from various publicly and privately available data sources, and methodologies for the curation of functional pathways varied between sources. In all, a total of 789 gene sets made up the collection that was employed in this study, and the approach was employed to facilitate potential mechanistic insight to the myogenic cell response to injury . As a control, we looked at the gene set for stress response at the 8 hour time point (Figure 3A). This probe set was significantly upregulated (p < 1.0e-36), relative to the basal expression distribution (Figure 3–3A). Furthermore, a subsidence in the stress response could be observed between 8 and 24 hours post-heat shock treatment as the gene set underwent a significant downregulation (p = 1.4e-8) that may reflect a normalization of the cell's condition (Figure 3B). These data clearly demonstrate the statistical strength of pathway analysis in the context of the heat shock response.
To gain better insight into the dynamic nature of regulated functional pathways, we traced the expression pattern of specific gene sets within the time course, such as the extracellular matrix protein-mediated signaling (EMP-M) (Figure 4, see Additional file 1). The analysis revealed a pattern of regulation that reflected an initial increase in pathway activity at 1, 2, and 4 hours, followed by a time-dependent decrease in expression through 8 and 24 hours post treatment. Therefore, the overall response in the EMP-M pathway to heat shock reflected an augmentation immediately following treatment, followed by a subsidence and eventual decrease in functional pathway expression after 24 hours. The transient nature of the pathway response was clearly defined with this bioinformatics approach, and is presented in Figure 4F.
A similar analysis was carried out on the BMP signaling pathway (see Additional file 1). When traced as a function of time over the 24 hour period, the BMP signaling pathway was initially decreased at 1, 2 and 4 hours at the transcriptional level, and recovered through 8 and 24 hours post treatment (Figure 5). This result was consistent with the established role of BMP signaling in antagonizing myogenic specification during development, and in regenerating adult muscle tissue . Furthermore, along with the stress response pathway and the EMP-M pathway expression traces, it appears that much of the initial response to heat shock was reversed within 24 hours of recovery.
Analysis of the most active gene sets through the entirety of the time-course revealed a notable phenomenon, typified by the dynamic profile of the protein biosynthesis pathway gene set. In response to heat shock treatment, the genes in this pathway underwent a significant and immediate upregulation, which gradually declined over the next 24 hours (Figure 6). The trend for this particular functional pathway suggested that protein turnover was highly regulated in the response to heat shock-induced stress, and this hypothesis was further supported by the apparent upregulation of several related pathways involved in the expression and degradation of proteins.
After 8 hours recovery from heat shock, a selection of gene sets related to core gene expression machinery and protein turnover, were significantly upregulated (Table 2). This general increase in protein metabolism pathways as an immediate response to injury may reflect the clearance and replacement (turnover) of the protein complement of the cell.
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.
An in-house implementation of the Gene Set Enrichment Analysis method  was used to analyze microarray data. As input, GSEA requires microarray data from two conditions to be compared (e.g. myoblasts vs myotubes). Genes with expression levels below 100 on more than 75% of the Genechips were discarded. Then, the relative expression level under condition1 and condition2 was computed as an expression ratio r i
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.
Charge SB, Rudnicki MA: Cellular and molecular regulation of muscle regeneration. Physiol Rev 2004, 84: 209-238. 10.1152/physrev.00019.2003
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.021
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.010
Conboy IM, Rando TA: Aging, stem cells and tissue regeneration: lessons from muscle. Cell Cycle 2005, 4: 407-410.
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-4
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/JCI200421696
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.027
Zhao P, Hoffman EP: Embryonic myogenesis pathways in muscle regeneration. Dev Dyn 2004, 229: 380-392. 10.1002/dvdy.10457
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.1032913100
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/ng1180
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.x
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.0814
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.
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.
Muchowski PJ, Valdez MM, Clark JI: AlphaB-crystallin selectively targets intermediate filament proteins during thermal stress. Invest Ophthalmol Vis Sci 1999, 40: 951-958.
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.
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.033
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.2004
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.1065874
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.251541198
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-3
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.M400674200
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.x
Storey JD, Tibshirani R: Statistical methods for identifying differentially expressed genes in DNA microarrays. Methods Mol Biol 2003, 224: 149-157.
Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci U S A 2003, 100: 9440-9445. 10.1073/pnas.1530509100
We thank Qicheng Ma for help in the development of GSEA.
JDS developed the bioinformatics algorithms and analysis. PAK designed and performed the heat shock experiment. CAM performed statistical analysis. JHL and SJE participated in the design and coordination of the study. NRN and DMK conceived of the study and drafted the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
About this article
Cite this article
Szustakowski, J.D., Kosinski, P.A., Marrese, C.A. et al. Dynamic resolution of functionally related gene sets in response to acute heat stress. BMC Molecular Biol 8, 46 (2007). https://doi.org/10.1186/1471-2199-8-46
- Heat Shock Treatment
- Follow Heat Shock
- Skeletal Myotubes
- Acute Heat Stress
- Multiple Hypothesis Testing Correction