- Research article
- Open Access
Patterns of gene expression associated with recovery and injury in heat-stressed rats
BMC Genomics volume 15, Article number: 1058 (2014)
The in vivo gene response associated with hyperthermia is poorly understood. Here, we perform a global, multiorgan characterization of the gene response to heat stress using an in vivo conscious rat model.
We heated rats until implanted thermal probes indicated a maximal core temperature of 41.8°C (Tc,Max). We then compared transcriptomic profiles of liver, lung, kidney, and heart tissues harvested from groups of experimental animals at Tc,Max, 24 hours, and 48 hours after heat stress to time-matched controls kept at an ambient temperature. Cardiac histopathology at 48 hours supported persistent cardiac injury in three out of six animals. Microarray analysis identified 78 differentially expressed genes common to all four organs at Tc,Max. Self-organizing maps identified gene-specific signatures corresponding to protein-folding disorders in heat-stressed rats with histopathological evidence of cardiac injury at 48 hours. Quantitative proteomics analysis by iTRAQ (isobaric tag for relative and absolute quantitation) demonstrated that differential protein expression most closely matched the transcriptomic profile in heat-injured animals at 48 hours. Calculation of protein supersaturation scores supported an increased propensity of proteins to aggregate for proteins that were found to be changing in abundance at 24 hours and in animals with cardiac injury at 48 hours, suggesting a mechanistic association between protein misfolding and the heat-stress response.
Pathway analyses at both the transcript and protein levels supported catastrophic deficits in energetics and cellular metabolism and activation of the unfolded protein response in heat-stressed rats with histopathological evidence of persistent heat injury, providing the basis for a systems-level physiological model of heat illness and recovery.
Heat illness, a continuum of disorders caused by hyperthermia, includes heat cramps, heat exhaustion, heat injury, and heat stroke [1–3]. Hyperthermia occurs when the body’s core temperature exceeds its hypothalamic set point as a result of overwhelming external or internal sources of heat or the effects of drugs or disease which impair heat dissipation [1–3]. Without adequate fluid replacement, dehydration and loss of blood volume lead to circulatory collapse, severe hyperthermia, and eventually heat stroke, defined as neurological dysfunction associated with a body temperature higher than 40°C (104°F) in humans [3–5].
Sudden death related to heat stroke occurs primarily in individuals with pre-existing cardiac abnormalities, organ failure secondary to rhabdomyolysis, or acute onset of organ failure . Systemic inflammatory response syndrome (SIRS) is considered to be the primary cause of organ dysfunction related to heat stroke [3, 7, 8], resulting in severe encephalopathy, rhabdomyolysis, acute renal failure, acute respiratory distress syndrome, myocardial injury, hepatocellular injury, intestinal ischemia, pancreatic injury, and hemorrhagic complications . SIRS may also be complicated by endotoxemia and/or disseminated intravascular coagulation through leakage of the gut epithelium  and/or damage to the vascular endothelium , respectively (particularly exertional heat stroke). Thus, multiorgan dysfunction syndrome is the primary reason for heat-related morbidity and mortality.
At the cellular and molecular levels, severe hyperthermia causes acidosis, hypoxia, and cellular fatigue. Thermal stress leads to an increased expression of heat shock proteins (HSPs) associated with thermotolerance [10–13]. The cellular response extends beyond HSPs to include networks of gene changes in organs and peripheral blood mononuclear cells after whole-body hyperthermia  or heat stroke . Moran and colleagues  describe at least three components of a time-dependent cascade of events related to the heat-stress response: 1) the involvement of HSPs, 2) a response that involves interferon-inducible genes (e.g., cytokines), and 3) a small, nonspecific stress response shared by other cell lines and stressors . Heat stress–induced molecular perturbations are also a likely contributing factor [10, 15–19] to altered heat intolerance .
Recently, we described an anatomically accurate, three-dimensional rat model of thermoregulation during heat stress , in which rats were placed at 37°C until their core temperature reached 41.8°C (Tc,Max); they were then allowed to recover. Here, we perform a global characterization of the heat-stress–induced gene response in liver, lung, kidney, and heart tissues at Tc,Max, and during recovery at 24 hours and 48 hours. Previous studies examined transcriptional alterations related to heat stress in various mammalian models and tissues [10, 11, 15, 17–19, 21, 22]. In this study, we use a whole transcript, gene array approach, which avoids issues related to 3’ bias, genes which lack poly-A tails, and confounders related to isoform variation . We present a pathway analysis to establish the regulatory framework of heat stress in our conscious rat model. Within this framework, we examine the consensus heat-stress response (cHSR) found in all four organs, and discuss genes not previously linked to in vivo heat stress. Finally, we identify changes in gene and protein expression in the heart which are anchored to histological evidence of organ injury, and we propose that mitochondrial dysfunction and injury pathways are involved in the observed cardiomyopathy.
Results and discussion
To induce the pathophysiological effects of heat stress, we placed conscious rats in an incubator at 37°C and monitored their core temperatures using remote telemetry implants (Figure 1A; see also Rakesh et al.). Animals were euthanized when their core temperature reached 41.8°C (Tc,Max), or they were returned to their normal housing environment and allowed to recover for up to 48 hours prior to sacrifice. We characterized the response of the animals to heat stress using a standard panel of hematological and clinical chemistry assays, as well as a histopathological evaluation of the liver, heart, kidney, and lung. In the context of these well-established assays, we also analyzed global changes in gene expression in the four tissues using microarrays and performed a “shotgun” proteomic analysis of cardiac tissue using mass spectrometry. This approach allowed us to link changes in gene and protein expression to tissue injury, the animals’ physiological state, and their overall health (Figure 1B).
Complete blood count and concomitant decrease in body weight were significantly different between heat-stressed animals and time-matched control animals only at Tc,Max (Additional file 1: Table S1). Further, blood urea nitrogen (BUN), a measure of kidney function, was significantly higher in heat-stressed animals, consistent with dehydration and/or kidney function impairment (Additional file 1: Table S1). As expected, no animals demonstrated significant heat-related histopathological changes at Tc,Max (Additional file 2: Table S2). The incidence and severity of lesions, such as cardiomyopathy and chronic progressive nephropathy, were generally minimal and consistent with the normal background lesions found in rats at 2–3 months of age (Additional file 2: Table S2).
Animals sacrificed at 24 hours demonstrated background lesions which were minimal in severity and not considered related to heat exposure (Additional file 3: Table S3). In one animal, we found diffuse, acute tubular necrosis in the kidney (Additional file 4, Figure S1A, B) and acute necrosis in the liver, mostly affecting hepatocytes in the periportal regions (Additional file 4: Figure S1C, D). Similar lesions did not occur in any other animal at any other sacrifice time point.
Although we found essentially no evidence of physiological or histopathological changes related to heat stress at 24 hours, we observed substantial cardiac injury in three out of the six animals sacrificed at 48 hours (Additional file 5: Table S4). All treated animals (and one control) had degeneration of cardiomyocytes with subacute inflammation consisting of mostly macrophages. This type of lesion is common in Fisher 344 rats (Figure 2A and B, time-matched controls), but the pattern of injury and severity of the degeneration and inflammation observed at the 48-hour time point suggested treatment-related tissue injury rather than an incidental, strain-related cardiomyopathy (Figures 2C and D). Extensive tissue damage often leads to kidney proteinosis similar to that observed in the animals most severely affected by cardiac inflammation (Additional file 5: Table S4). Thus, while the physiological parameters at Tc,Max suggest that the animals had undergone heat-related dehydration, the development of tissue injury appeared to be a delayed process. No injury was evident at 24 hours, but there was clear, persistent, heat-induced cardiac damage in some animals at 48 hours.
The consensus heat-stress response: pathway analysis of the transcriptome at Tc,Max
In order to elucidate the biological events related to heat stress at the molecular level, after heat exposure we identified genes which were differentially expressed in the heart, lung, liver, and kidney using Partek® Genomics Suite software (version 6.3 Copyright © 2008 Partek Inc., St. Louis, MO, US) (Additional file 6: Table S5), and we performed a functional enrichment analysis of the genes using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) [24, 25]. Gene set analyses were conducted in Partek to determine the top enriched gene ontology (GO) terms for each condition (Additional file 7: Table S6). At Tc,Max, 467 unique biological processes (BPs), cellular components (CCs), and molecular functions (MFs) were enriched across all organs, of which only 25 were common to all four organs, and 156 were common at least to two organs. Using DAVID, we also determined enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (Additional file 8: Table S7). The differentially expressed genes (DEGs) from all four organs were enriched for functional terms related to protein folding and/or response to unfolded protein, regulation of apoptotic process, and the response to lipopolysaccharide (LPS)/cytokine in the gene set enrichment analysis. There was also a set of stress and transcriptional responses associated with the mitogen-associated protein kinase (MAPK) KEGG pathway (Table 1).
At Tc,Max, we identified 78 differentially expressed transcripts in all four organs (Additional file 6: Table S5), comprising a cHSR. We examined the cHSR further with DAVID tools and found that 51 nonredundant transcript identifiers (IDs) were identified in 11 functional annotation clusters (Additional file 9: Table S8). The most significantly enriched terms in the dataset were then subcategorized as “protein folding” (PF) and “regulation of apoptosis” (RA) (Table 2). Twenty-seven transcripts were not functionally clustered by DAVID (Additional file 6: Table S5), but literature supports the roles related to PF [26–34] and/or RA [35, 36] in 14 of the 27 (Additional file 6: Table S5). The remaining 13 transcripts had unknown functions or were related to inflammation or other biological processes [37–41] (Additional file 6: Table S5).
The unfolded protein response is activated in response to the accumulation of unfolded or misfolded proteins in the lumen of the endoplasmic reticulum (ER). Once activated, global protein translation is halted and chaperones (e.g., HSPs) are produced to manage protein folding. If the unfolded protein response fails, a cellular alarm response occurs, leading to autophagy and recovery [42–48] or apoptosis  and/or cytolysis . At the molecular level, imbalance of the proteostasis network results in a cascade of responses to cellular stressors, involving molecular chaperones, clearance mechanisms, detoxifying enzymes, and an accumulation of misfolded aggregation-prone and proteotoxic species, which can lead to cell-specific dysfunction [50, 51]. The heat shock factor 1 (Hsf1) cycle can be initiated by proteotoxic stress, ultimately leading to the expression of heat shock proteins [10, 15], however, less is known about how heat stress, proteostasis networks, apoptosis, and autophagy contribute to organ-specific injury. In the heart, aberrant folding and unfolded protein response signaling lead to a number of pathologies, including myocardial ischemia, cardiac hypertrophy, and heart failure [52, 53]. Thus, the enriched pathways in protein folding and unfolded protein response at Tc,Max may point toward a global imbalance of the proteostasis network .
The consensus heat stress response: differential gene expression in protein folding and regulation of apoptosis pathways at Tc,Max
The environmentally induced heat shock transcripts and heat-stress–related transcripts enriched in PF and RA displayed the largest fold changes in heat-stressed animals relative to controls at Tc,Max (Table 2). Interestingly, our three-dimensional thermoregulatory model for the rat predicted that the liver temperature in heat-stressed animals exceeds that of other organs . This prediction is consistent with our findings that the largest changes in temperature-dependent HSPs occurred in the liver.
A wide variety of physiological and environmental stimuli can induce expression of HSPs, including physical, chemical, and biological insults (e.g., inflammation) . Genes encoding HSPs are divided into five superfamilies based on molecular weight: the small HSPs , the Hsp40 (DnaJ) family, the Hsp70 (Hspa) family, and the Hsp90 (HspC) and Hsp110 (HspH) families . As molecular chaperones, HSPs play central roles in 1) protein folding, aggregation, transport, and/or stabilization; 2) cytoprotection, by preventing post-mitochondrial apoptosis in caspase-dependent pathways (e.g., Hspb1, Hspa1b, and Hsp90) and/or caspase-independent pathways (e.g., Hsp70) ; and 3) cytoprotection, by preventing apoptosis at the premitochondrial level through inhibition of stress-activated kinases (i.e., Jnk1) in an adenosine-triphosphate (ATP)-independent manner (e.g., Hsp70) . Although a recent report demonstrates a heat-stress–induced increase in rat liver expression of the inducible HSP70 protein (HSP72) but not the constitutive forms (HSP70 or HSP73) , our study demonstrates large increases in the products of both constitutive and inducible HSP genes at Tc,Max (Table 2), including Hspb1, Hspa1b, Hspe1, Serpinh1, Hspd1, Hsp90aa1, Hsph1, and several members of the DnaJ (Hsp40) chaperone family, transcripts homologous to genes coding for HSPs [i.e., LOC498996 (similar to DnaJ homolog subfamily A, member 4), LOC500476 (similar to Hsp8), and RGD1561150 (similar to genes encoding CPN10-like protein)], and heat-stress–related genes (e.g., Slc5a3, SNORD14, Ahsa2, Chordc1, and Banp). SNORD14 encodes a small nucleolar RNA which is transcribed from intronic regions of the Hsp70 gene in mice . Heat-stress–induced expression of SNORD14 may facilitate the repair process involving ribosomal RNA after heat injury . The environmentally inducible HSP gene, Hspa1b, was expressed at very high levels in all tissues (ranging from 24- to 209-fold), although one copy of the gene did not functionally cluster with the PF response in DAVID (Additional file 6: Table S5). Hspa1b codes for an ATP-dependent chaperone regulating intracellular anti-apoptosis, and is considered to be the most universal stress-induced HSP [54, 60, 61]. In baboons, expression of Hsp72 is proportional to severity of heat stress .
Some of the differentially expressed HSPs at Tc,Max are also found in the extracellular environment (including Hspb1 and Hspa1b), where they are thought to have immunogenic functions (Table 2) . Putative mechanisms for HSP release into biofluids include cell death, active non-classical secretory pathways, lysosome-endosome pathways, release by secretory-like granules, or insertion into the lipid bilayers of export vesicles. The protein product of extracellular Hsp72 binds to CD14, toll-like receptor (TLR) 4, and CD40 on dendritic cells, monocytes, and macrophages [63–65] and induces expression of NF-κB and chemokines (MIP-2 or CXCL2) by a TLR 2/4-dependent mechanism [66, 67]. It remains to be determined whether products of the differentially expressed HSP genes in this study are actually released into the extracellular space in our heat stress model.
The apoptosis pathway was also significantly enriched in the KEGG analysis. In the heart, liver, and kidney, activating transcription factor 3 (Atf3) was significantly induced (Table 2). The protein product of Atf3 is a member of the ATF/cyclic AMP response element-binding (ATF/CREB) family of transcription factors [68, 69]. The product of Atf3 is a critical negative regulator of pro-inflammatory cytokine production in response to TLR activation . Atf3 knockout mice exhibit no developmental abnormalities  unless challenged with LPS or another stress stimulus, resulting in an overt susceptibility to endotoxic shock and death . Zfand2a, also known as arsenite-inducible RNA-associated protein (AIRAP), was also up-regulated in all four organs in heated animals (Table 2). In rats, Zfand2a is closely associated with the febrile response, a heat-associated physiology. Although early experiments in rodent fibroblasts suggested that Zfand2a was induced by arsenite stress but not heat shock, later experiments in human cell lines demonstrated that Zfand2a was induced by heat stress in an HSF-1–dependent manner and closely paralleled Hsp70 expression after heat stress. While we can associate most of the members of the cHSR with heat stress, a few genes are associated with other stressors and some are uncharacterized (see Additional file 6: Table S5). Up-regulated apoptosis-related genes included Uspl1, MGC95208, and Dusp8[72, 73]. How the products of these genes contribute to the consensus response will require further investigation.
No transcripts were down-regulated in the cHSR, although several transcripts were down-regulated in at least three organs and were near the threshold for the fourth organ. For example, Cd180 was down-regulated in the liver, kidney, and lung, but was just below the two-fold significance threshold in the heart (i.e., 1.99-fold down-regulated; see Additional file 6: Table S5). Cd180 codes for a cell surface marker in the TLR family that forms the receptor RP105/MD1 and negatively regulates the expression of TLR4 [74, 75]. Reduction of the Cd180 expression by heat stress may promote TRL signaling and inflammation.
Transcriptomic profile at 24 and 48 hours
Our results at Tc,Max provide evidence of genes and functional regulatory pathways involved with the initial in vivo heat-stress response. Most of the DEGs at Tc,Max in the cHSR returned to control levels by 24 hours after heat stress (Additional file 6: Table S5). These results are consistent with the blood chemistries, suggesting that animals achieved recovery of most organ function by 24 hours. At Tc,Max and 24 hours, heat-stressed animals showed limited histopathological effects, although one animal showed evidence of liver histopathology at 24 hours (Additional file 3: Table S3 and Additional file 4: Figure S1). However, at 48 hours in three out of six animals, both histopathology and the pattern of gene expression indicated evidence of persistent, moderately severe damage to heart tissue.
We explored the patterns of gene expression in our experimental animals using self-organizing maps, and observed that the DEGs in the hearts of three out of six heat-stressed animals at 48 hours differed from both controls and the other heat-stressed animals (Figure 3 and Additional file 8: Table S7; see Additional file 8: Table S7 and Additional file 10: Figure S2 for a self-organizing map analysis of the single animal with liver injury at 24 hours). Strikingly, these same animals also displayed cardioinflammation in the histological analysis. Node A (Figure 3) contains 1,351 down-regulated transcripts and Node B contains 887 up-regulated transcripts (see Additional file 8: Table S7). We next performed a functional analysis of the DEGs in Nodes A and B using DAVID, because these nodes were evident in the three animals with the highest histopathological scores. The down-regulated genes were associated with highly enriched KEGG pathways, including oxidative phosphorylation (enrichment score of 6.7, 63 genes), cardiac muscle contraction (enrichment score of 6.2, 33 genes), and pathways associated with chronic proteotoxic stress [Parkinson’s disease (enrichment score of 6.3 fold, 61 genes), Huntington’s disease (enrichment score of 4.8, 63 genes), and Alzheimer’s disease (enrichment score of 4.6, 62 genes)]. The oxidative phosphorylation pathway is central to Parkinson’s Disease, Huntington’s Disease, Alzheimer’s Disease, and Cardiac Muscle Contraction, suggesting that mitochondrial dysfunction and disruption in energetics are key features in animals with heat-stress-induced heart injury. Mitochondrial dysfunction, disruption of ER and membrane trafficking, alterations in protein folding and clearance, and activation of inflammation are all thought to be underlying causes of protein-folding disorders in neurons and other organ systems [76–79]. Primary and secondary aggregation of peptides and proteins disrupt normal protein homeostasis, leading to pathological states involving aberrant protein aggregation [76–79]. A recent study shows a strong correlation between cytological evidence of aggregated proteins and the functional pathways we observed in our clustering analysis (i.e., Oxidative Phosphorylation, Parkinson’s Disease, Huntington’s Disease, Alzheimer’s Disease, and cardiac muscle contraction ). Further, supersaturated protein complexes in protein aggregates are associated with highly enriched Ribosome and Proteasome pathways. The Ribosome pathway is enriched in genes upregulated in node C (+2.8 fold), cell cycle (+5.7 fold), and DNA replication (+6.3 fold) (Additional file 8: Table S7). Thus, the three animals with cardiac injury in response to heat stress show genomic signatures consistent with neurodegenerative diseases at the pathway level, suggesting that heat-injured animals have transcriptomic evidence of a greater propensity of protein to aggregate.
Proteomic analysis of heart tissue at Tc,Max, 24 hours, and 48 hours
To determine whether the pathway analysis was complementary at the protein level, proteins were extracted from heart tissue, analyzed by iTRAQ (Additional file 11: Table S9; data includes proteins with at least a 1.3-fold change in expression when compared to time-matched controls), and compared with gene expression at all three time points. Proteins and gene expression data were largely discordant at Tc,Max (Figure 4). This time point had significant enrichment in genes encoding proteins involved in the heat-stress response, especially the HSPs (Figure 4, genes in bold). However, there were too few differentially expressed proteins at Tc,Max to derive functional clusters (Additional file 11: Table S9). Although HSP70 and the gene encoding this protein were the most upregulated gene/protein pair at Tc,Max, HSPB1 was the only other heat shock protein modulated in the protein data at this time point, and its expression decreased while the transcript encoding the protein product increased. Proteins associated with serine proteinase activity were upregulated, including α-1 macroglobulin precursor and the acute phase protein α-1-inhibitor 3 precursor (MUG2, A1I3), but the transcripts for these products were below the fold-change threshold. Both proteins are serum proteins and changes may reflect changes blood contamination not necessarily changes in the heart tissue itself. Further, several mitochondrial proteins were modulated, including mitochondrial precursor hydroxyacyl-coenzyme A dehydrogenase and mitochondrial import inner membrane translocase subunit, but the transcript expression did not correlate with the protein expression.
At 24 hours, 88 proteins were significantly altered (Additional file 11: Table S9), but few protein/gene expression pairs were in concordance (HSPA1B and its corresponding transcript; Figure 5). By 24 hours, multiple HSPs were expressed above the fold-change threshold (i.e., HSPA5, HSP90B1, HSPB1, HSPA7, and HSP90AA1; bold in Figure 5). However, all of the transcripts encoding the HSPs which had been up-regulated at Tc,Max were either down-regulated (HSPA5) or below the fold-change threshold, with the exception of HSPB1 (Figure 5). HSP90AA1, HSP90B1, and calreticulin precursor (CALR) clustered to the functional KEGG pathway Antigen Processing and Presentation (Figure 5, inset), although their respective transcripts did not reach the fold-change threshold (Figure 5). The enrichment of the Complement and Coagulation cascades KEGG pathway was driven by up-regulation of the proteins KNG1, KNG2, KNG1L1, C4, F2, and FGB, although there was not a complete correlation between transcript/protein expression level and direction (Figure 5).
Thus, at 24 hours environmentally inducible HSPs remained elevated, despite a lack of blood chemistry or histopathologic evidence of injury. The HSPs elevated at 24 hours are associated with cardiomyocyte injury and immune-mediated functions (e.g., Antigen Processing and Presentation). HSPA5 is an ER-resident chaperone associated with hypoxia in cardiomyocytes [81, 82]. HSPA5 retargets to mitochondria and contributes to the mitochondrial unfolded protein response . Extracellular HSPs associated with Antigen Processing and Presentation have been associated with immune cell recruitment to cardiomyocytes, a potential mechanism of inflammatory injury .
At 48 hours, we subdivided the animals based on histopathologic evidence of injury (Figure 6) or recovery (Figure 7). Five hundred and twenty-one proteins were significantly altered in animals with histopathological evidence of heart injury (Figure 6 and Additional file 11: Table S9), whereas only 71 were modulated in heat-stressed animals with no histopathological evidence of injury (Figure 7 and Additional file 11: Table S9). The up-regulated proteins were enriched for pathways associated with protein folding (Ribosome; Figure 6, inset) and, by extension, propensity to aggregate . Pathways associated with immune dysregulation were enriched (Complement and Coagulation, including concordant transcript and protein expression of KNG1, KNG1L1, KNG2, C3, C4, CFD, F2, FGB, FGG, SERPINC1, SERPIND1, and SERPINH1; and Antigen Processing and Presentation, including concordant protein expression of HSPA1B, HSP90AA1, HSPA5, CALR, HSPA8, and HSPAB1). HSP expression correlated with transcript expression only in the heat-injured animals (red squares in Figure 6), not the heat-stressed animals (Figure 7). Down-regulation of pathways enriched for Oxidative Phosphorylation, citrate cycle, and fatty acid metabolism suggested an energy crisis with massive mitochondrial dysfunction, supporting the observed necrosis in the cardiac histopathology. In addition, KEGG enrichment analysis of down-regulated proteins identified pathways associated with chronic proteotoxic injury (Parkinson’s Disease, Huntington’s Disease, Alzheimer’s Disease; Figure 6, inset) . Protein/gene pairs which showed concordant regulation in 48 hour heat-injured animals were diagrammed by subcellular organelle (Figure 8 and Additional file 12: Table S10 and Additional file 13: Table S11). Most of the protein/gene pairs associated with mitochondria were associated with the Oxidative Phosphorylation, with nearly 100% overlap with protein/gene pairs in the KEGG pathways of Parkinson’s Disease, Huntington’s Disease, and Alzheimer’s Disease. Likewise, the heat shock proteins and proteins/gene pairs associated with proteasomal degradation were localized to the proteasome, including biomolecules in the immunomodulatory pathway Antigen Processing and Presentation. The Cardiomyocyte Contraction pathway contained protein/gene pairs which were downregulated in conjunction with the sarcomere.
Protein aggregation and supersaturation by subcellular localization in heart tissue
Supersaturation (i.e., high cellular concentration of proteins with a high propensity to aggregate) of proteins has been associated with some of the KEGG pathways we found modulated by heat stress. Many of the KEGG pathways we identified as enriched in heat-injured animals at 48 hours have been associated with neurodegenerative protein misfolding disorders, including Oxidative Phosphorylation (and the related pathways Alzheimer’s Disease, Parkinson’s Disease, and Huntington’s Disease), Cardiomyocyte Contractility, and Proteasome . In addition, proteins in these biochemical pathways have a combined high propensity to aggregate and elevated concentrations in proteotoxic diseases characterized with protein misfolding . Because heat induces protein misfolding and dysregulation, we investigated whether proteins affected by heat stress and heat injury had higher propensities to aggregate as a function of subcellular organelle localization (Figure 9 and Additional file 12: Table S10).
The supersaturation scores for proteins modulated after heating at 24 hours and at 48 hours in heat-injured animals (Figure 9; p<0.05), but not at Tc,Max, was higher than for proteins whose abundance was unchanged at these time points. Average supersaturation score of proteins modulate in animals without evidence of cardiac injury at 48 hours was also higher, but was not significantly different from the average for the unchanged proteins (in part due to a smaller number of proteins changing in abundance at 24 hours vs 48 hours). These data indicate that: (1) preferential modulations in abundance of proteins with high supersaturation scores occur as early as 24 hours; and (2) proteins prone to aggregation are more likely to be regulated in cardiac injured rats than in animals that do not show histopathological evidence of injury at 48 hours. This suggests that protein aggregation is contributing to changes in protein abundance and perturbation of the biological processes associated with these proteins. Further, aggregation as a result of heating is possibly a leading mechanism in the heat injury observed at 48 hours.
Here, we present a global characterization of gene and protein response to heat stress in a conscious rat model that provides new insight into the regulatory processes associated with heat stress and recovery. 1) Transcriptional analyses of the heart, lung, liver, and kidney at Tc,Max revealed a consensus heat stress-induced gene response that was enriched for genes related to protein folding and the regulation of apoptosis. 2) We used self-organizing maps to group genes into clusters of similar expression profiles in animals with evidence of similar histopathologies. In the hearts of three animals at 48 hours, significantly enriched Oxidative Phosphorylation, Parkinson’s Disease, Huntington’s Disease, Antigen Processing and Presentation, and Cardiac Muscle Contractility KEGG pathways were consistent with protein-folding disorders, mitochondrial dysfunction, and perturbation in cellular energetics. 3) Transcriptional responses were corroborated by the results of proteomic analysis in heat-injured cardiac tissue at 48 hours, with pathway analysis, protein aggregation scores, and subcellular localization supporting up-regulation of proteins and genes in Ribosome, Antigen Processing and Presentation, Regulation of Actin Cytoskeleton, Complement and Coagulation Cascade, and Pentose Phosphate Pathway, and downregulation of Hypertrophic/Dilated Cardiomyopathy, Fatty Acid Metabolism, and Oxidative Phosphorylation (Figure 8). Proteins involved in Antigen Processing and Presentation emerge as early as 24 hours. Disruption in these mitochondrial mediators and deficits in pathways related to energetics precede large changes in mediators of oxidative phosphorylation and protein folding disorders by 48 hours in heat-injured cardiac tissue. 5) Supersaturation scores indicated an increased propensity to aggregate in proteins changing in abundance after heat stress at 24 hours and 48 hours. Increased propensity to aggregate at 48 hours in animals with histopathological evidence of heat injury was more pronounced than in animals without cardiac injury. Disruption in protein folding and the resulting aggregation is likely to contribute to the persistent subcellular indicators of energy crisis and regulation of apoptosis observed in heat-injured animals.
The Fischer 344 strain was chosen because they do not gain body weight as quickly as Sprague–Dawley, which can have a large impact on thermal dynamics. It is possible that naturally occurring lesions in the heart were exacerbated by heat stress, however, only those with severe lesions were correlated with animals that demonstrated both gene and protein alterations (i.e. uninjured animals, with or without background lesions, were indistinguishable at the molecular level). It is important to note that many humans have pre-existing conditions leading to heat susceptibility, and it is possible to extrapolate our findings to these populations. Given that most heat stroke deaths occur in the elderly, this would be particularly relevant to those populations, although not exclusive. Future studies should attempt to corroborate the findings of our study using rat strains less susceptible to the congenital cardiomyopathy phenotype. Further, intermediate time points between Tc,Max, 24 hours, and 48 hours should be investigated to more accurately chronicle the initiation and progression of organ injury in the heart. However, our model may prove useful in identifying and understanding mediators of heart injury after heat stress in susceptible individuals.
This work in heat-stressed animals provides the basis for an assessment of global gene and protein responses in future heat-stroke–induced organ injury experiments. Top-down integration of heat-induced changes in metabolomics, proteomics, and transcriptomics (to include small RNA) will provide the foundation for a computationally based linkage to human three-dimensional thermoregulation models. Anchoring these systemic stress responses to the physiological model of heat stroke will provide further insight into predicting risk, severity, and timing of organ injury in response to hyperthermia. Such a model will potentially accelerate the development of strategies to improve prevention, classification, and treatment of heat-related illnesses.
Animal model and tissue collection
In vivo rat experiments were performed at the US Army Research Institute of Environmental Medicine. The Institutional Animal Care and Use Committee approved all experimental procedures, which complied with the American Physiological Society’s guiding principles for research involving animals and the Guide for Care and Use of Laboratory Animals. As previously described , we utilized male Fischer 344 rats (F344; n = 36; Charles River Laboratories, Stone Ridge, NY) weighing 234–336 g (~2–3 months old). Briefly, rats were housed under standard laboratory conditions (22°C, 12:12 hours light:dark cycle, lights on at 6:00 AM) in an Association for Assessment and Accreditation of Laboratory Animal Care–accredited facility. We provided chow (Harlan Teklad, LM-485, Madison, WI) and water ad libitum. Rats were implanted with TL11M2-C50-PXT PhysioTel® Multiplus Transmitters (Data Sciences International, St. Paul, MN) to measure core temperature (Tc; ±0.25°C), heart rate (HR; bpm, beats per minute), and mean arterial pressure (MAP; ±3 mm Hg).
We conducted all experiments in conscious, free-moving animals, as previously described . We used 36 rats with remote telemetry implants for this study (Figure 1). Briefly, we placed rats in the heat-stress arm of the experiment into a floor-standing incubator (Thermo Scientific, Ashville, NC) set at room temperature (22°C) 24 hours prior to initiation of heat stress experiments. Rats in the control arm (n=18) were never introduced to the incubator environment and were kept at a normal housing temperature. Once Tc ≤ 37.3°C was measured (Figure 1), we weighed the control rats and returned them to their original cages for the remainder of the experiment. Rats in the heat-stress arm were heated in the incubator at 37.0 ± 0.2°C until reaching a Tc of 41.8°C (Tc,Max). Six heat-stressed rats were weighed at Tc,Max and euthanized as previously described . The remaining rats were placed in a new cage at the ambient housing temperature (22.0 ± 0.2°C) and euthanized at 24 hours (n=6) or 48 hours (n=6) (Figure 1). Six time-matched controls were euthanized at times corresponding to Tc,Max, 24 hours, and 48 hours in the heat-stress arm. Time points were chosen to maximize the likelihood of observing changes in gene expression and to reflect the clinical observation that recovery from heat exhaustion usually takes 24 to 48 hours following mild injury. We provided control and experimental animals with food and water ad libitum throughout recovery.
Physiological and hematological parameters
We measured body weight (BW) on a top-loading balance with an accuracy of ±0.1 g. The difference between pre- and post-heat BW was used as an indicator of percent dehydration using the following calculation: [(post–heat-stress BW – pre–heat-stress BW)/pre–heat-stress BW] × 100. We did not take into account changes in body weight due to urine or feces loss. Rats were deeply anesthetized under isoflurane anesthesia, and a thoracotomy was performed, followed by a cardiac puncture to collect blood for hematologic analysis. We determined hemoglobin (hgb) and hematocrit (hct) using a VetScan HM5 Hematology System (Abaxis, Union City, CA), and BUN (mg/dL) using a handheld iStat clinical analyzer (Abbott Laboratories, Abbott Park, IL) and 95 μL of the collected blood. We placed the remaining blood in heparinized tubes and immediately stored it on ice. Plasma was isolated via centrifugation (4°C; 5 minutes, ~2,000 rpm), placed in aliquots, and stored at -80°C until assays were performed. We determined plasma alanine aminotransferase (ALT) (U/L) and aspartate aminotransferase (AST) levels (U/L) in duplicate using the Poly-Chem® System (Polymedco Inc., Cortlandt Manor, NY) per manufacturer specifications. Precision of the ALT and AST assays were valid to a minimum of ±10 U/L and ±12.1 U/L, respectively. Statistical significance among groups was determined by Kruskal-Wallis analysis of variance with post-hoc Tukey’s honest significant differences test relative to control, with an α of 0.05 denoting statistical significance.
Fixed tissues (heart, liver, lung, and kidney) were mounted and hematoxylin and eosin (H&E)-stained (IHC World, Woodstock, MD). A total of 24 rats were separated into three groups: Tc,Max, 24 hours recovery, and 48 hours recovery. Each group comprised six heat-stressed and two time-matched control animals. We performed serial sections to produce 20 unstained slides per tissue, with three to five sections per slide. Experimental Pathology Laboratories® (Sterling, VA) conducted the histopathological evaluation. The magnitude of inflammatory or degenerative lesions was graded on a scale of 1 to 5, with grade 1 being minimal and grade 5 being severe (See Additional file 1: Table S1, Additional file 2: Table S2, Additional file 3: Table S3).
Isolation of RNA
Frozen tissues (heart, lung, liver, and kidney) stored at -80°C in cryovials were placed on dry ice, cut into aliquots with a sterile scalpel on a pre-chilled titanium block, and placed in new tubes pre-chilled on dry ice. The block was washed with RNase Zap (Ambion, Life Technologies, Grand Island, NY) between samples. While working on dry ice, whole tissues (~25 mg) from the lung, liver, and kidney were cut from each sample and placed in 700 μl of QIAzol Lysis Reagent. Samples were homogenized with a TissueLyzer LT (Qiagen, Valencia, CA) for 5 minutes at 25 Hz twice, allowed to sit at room temperature for 5 minutes, then placed at -80°C.
We homogenized heart tissue with a 6750 Freezer/Mill with three sample microvials (Spex SamplePrep, Metuchen, NJ). Approximately 1 cm3 of heart tissue was cut, dipped into liquid nitrogen, and then placed into pre-assembled and pre-chilled microvials containing 200 μL of QIAzol. An impactor was added to each microvial, securely capped, and placed in the mill. The Freezer/Mill settings were as follows: T1 run time = 3 min, T2 intermission = 1 min, T3 pre-cool = 5 min, rate = 15, and cycles = 2. After completion, we added 500 μL of QIAzol to each vial, allowed the vial to sit at room temperature for 5 minutes, transferred the sample to a microfuge tube, and then placed it at -80°C.
The following day, we allowed all samples to thaw and isolated the RNA using the miRNeasy 96 kit (Qiagen), according to manufacturer’s instructions. The quality and quantity of RNA samples were evaluated with a 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA), using the Agilent RNA 6000 Nano Reagents, and a multiwell Nanodrop 8000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA).
Affymetrix gene array
We prepared cDNA from total tissue RNA using the Ambion Whole Transcript (WT) Expression Kit (Ambion) according to manufacturer’s recommendations, fragmenting the cDNA using the GeneChip WT Terminal Labeling Kit (Affymetrix, Santa Clara, CA). Following manufacturer’s instructions, we prepared labeled cDNA, trays, and arrays with the GeneTitan Hybridization, Wash and Stain Kit for WT array plates (Affymetrix). Samples were then applied to the Rat Gene 1.1 ST 16 Array Plate or 24 Array Plate and placed in the GeneTitan System (Affymetrix), according manufacturer’s recommendations. Of the 144 arrays, 2 liver, 1 lung, and 1 heart array did not pass the GeneTitan scanning quality control. Thus, 4 arrays (2.7%) were excluded from further analysis.
Gene array data analysis
For each organ, we normalized arrays using an extension of the PLIER (Probe Logarithmic Intensity Error) algorithm, called the iterPLIER procedure, in the Affymetrix Expression Console . The iterPLIER (exon level) procedure discards feature sets which perform poorly, as described by Qu et al.. We imported the resulting CHP files into Partek. Affymetrix library files included all available reference files related to RaGene-1_1-st-v1.na32.rn4. To address signal variance following the iterPLIER analysis, we applied a variance-stabilizing transformation (PLIER+16) to raw intensity (Guide to PLIER Estimation) . Standard deviation for each transcript ID were determined across all samples for each organ, and we removed those lower than the median and retained the remaining transcript IDs for determination of differentially expressed genes and as background for our pathway analysis of each organ. The low-variance criterion from Bourgon et al. was implemented by computing and sorting the expression variance of each gene over the complete condition set, then removing the low-variance genes (lower half). To determine DEGs, a two-way analysis of variance (ANOVA) was conducted for each organ, and gene lists were generated by applying a two-fold change cutoff with the Benjamini-Hochberg False Discovery Rate (FDR <0.05) [87, 88]. Microarray data described in this study have been deposited in the Gene Expression Omnibus database with accession number GSE56740.
Gene Set enrichment analysis
DEGs were analyzed for enriched GO terms using Partek. We considered only terms with <500 members and removed terms with <4 members identified. Additionally, only terms with enrichment scores above two-fold with p-values below 0.05 were reported. Gene lists were examined with the gene set analysis tools in Partek software and DAVID [24, 25]. Enriched KEGG pathways and FACs (functional annotation clusters) were determined using default settings for each organ. Transcript IDs from functionally associated clusters with an enrichment score >2.0 were consolidated and redundant genes removed. Gene lists were consolidated and sub-categorized to account for redundant annotation terms comprising FACs. For representation purposes, unique transcripts were consolidated from similar FACs.
ITRAQ proteomic analysis
We dried 100 μg of protein in a centrifugal evaporator and prepared the sample for mass spectrometry according to the iTRAQ reagent protocol (SCIEX, Framingham, MA). Samples were reduced in 50 mM tris-(2-carboxyethyl) phosphine, cysteines were blocked in 200 mM methanethiosulfonate, and proteins were digested in 5 μg trypsin overnight. We fractionated the samples using a Polysulfomethyl A column (200 × 4.6 mm, 5 μm, 1000 Å, PolyLC, Columbia, MD) at a flow rate of 0.95 μL/min. An Agilent 1100 high performance liquid chromatography (HPLC) equipped with an autosampler with an expanded injection loop and needle seat, diode array detector, fraction collector and Chemstation data system was used for the separation. The column compartment was set to 35°C, and 1 mL of the sample was injected. The solvents used for the separation were solvent A (distilled water), solvent B (500 mM ammonium formate, and 3% formic acid in water), and solvent C (acetonitrile). The initial conditions were 73% A: 2% B: 25% C. A linear gradient was performed to 75% B: 25% C in 20 minutes, with a stop time of 30 min. One-minute fractions were collected starting at 4 minutes and continuing until 19 minutes, producing 15 fractions. Samples were dried and reconstituted in 80 μL of 0.1% formic acid. Following reconstitution, fractions 1, 13, 14, and 15 were combined into one sample, dried, and reconstituted again in 80 μL of 0.1% formic acid.
We separated the peptides using a Thermo Scientific Proxeon EASY nano-LC prior to direct injection into a Thermo Orbitrap Velos mass spectrometer (ThermoScientific Corp, West Palm Beach, FL). The peptides were trapped using an Agilent Zorbax C-18 5 × 0.3 mm, 5 μm particle size column. The analytical separation was performed using a Waters NanoAcquity UPLC column BEH C18 100 μm × 100 mm, 1.7 μm particle. Ten microliters of the sample were loaded on the column and flushed with 60 μL of 0.1% formic acid in water. The solvents used were 0.1% formic acid in water (Buffer A) and 0.1% formic acid in acetonitrile (Buffer B) were used. The gradient profile was as follows: initially 3% B, 45% B at 90 minutes, 95% B at 100 minutes, 95% B at 115 minutes, 3% B at 116 min, and the analysis stopped at 120 minutes. The flow rate was 0.30 μL/min. Precursor selection was done in a Thermo Orbitrap using 100,000 resolution. The top 10 peptides, based on intensity, were selected for fragmentation. Sample sets were analyzed with a single scan, using stepped higher-energy collisional dissociation from 40 and 50 V. Species with single or unidentified charge states were excluded from precursor ion selection. Dynamic exclusion was used with a 20-second window.
We processed mass spectral data using Thermo Proteome Discoverer 1.3 (ThermoFisher Scientific) and searched for samples using the Sequest-generated randomized rat database based on all Rattus norvegicus sequences found in the National Center for Biotechnology Information (NCBI) Reference Sequence Database (4-12-2012), with the mass tolerance values set at 10 ppm for precursor ions and 0.8 Da for fragment ions . Other parameters included fixed-modification methylthio (cysteine) and iTRAQ™ 8-plex (N-terminus/lysine). The variable modifications allowed were oxidation (methionine), phosphylation (serine/tyrosine/threonine), acetylation (histidine/serine/tyrosine/threonine), modification by N-acetylglucosamine (asparginine/serine/threonine), modification by hexose-linked N-acetylglucosamine (serine/threonine) and iTRAQ™ 8-plex (tyrosine). Scans having peptide identifications with a 1% FDR were used for quantification. Reporter ions were quantified with a mass tolerance window of 20 ppm using the most confident centroid. Reporter ion intensities were corrected using the Discoverer 1.3 default method of iTRAQ™ 8-plex mass tags by Applied Biosystems optimized for ThermoFisher Scientific.
Quantitative iTRAQ data analysis
To quantitate peptide and protein identifications from MS/MS spectra, we used Scaffold Q+ (version Scaffold 4.2.1, Proteome Software Inc., Portland, OR). Peptide identifications of greater than 95.0% probability by the Scaffold Local FDR algorithm were accepted, and protein identifications of greater than 99.0% probability with at least two identified peptides and probabilities were assigned by the Protein Prophet algorithm . To satisfy the principles of parsimony, we grouped into clusters the proteins containing similar peptides not differentiable by MS/MS analysis alone, correcting channels as described in i-Tracker . Acquired intensities in the experiment were globally normalized across all acquisition runs, and individual quantitative samples were normalized within each acquisition run. We normalized intensities for each peptide identification within the assigned protein, and normalized the reference channels to produce a 1:1 fold change. All normalization calculations were performed using medians to multiplicatively normalize data. Differentially expressed proteins were determined using a ±1.3-fold threshold as described by Ahn et al..
Calculation of supersaturation scores
Supersaturation scores were calculated based on the methods of Ciryam et al.. Briefly, supersaturation scores are the sum of propensity to aggregate and the log of the protein concentration normalized to the median score for all proteins detected. The propensity to aggregate for unfolded proteins (Zagg) and proteins in their native state (ZSCagg) were kindly calculated by Tartaglia and Venduscolo for all proteins identified by our MS analysis . Protein abundance was estimated based on the sum of the total ion intensity for the top three most abundant peptide ions, which has been shown to correlate well with protein abundance . Precursor ion intensities determined by extracted ion chromatograms in Proteome Discoverer were summed across all fractions and charge states. For peptides where multiple precursor intensities were reported within a fraction, only the highest intensity was used. Only proteins with three identified peptides were used in this analysis. Total ion intensity for all samples in a single iTRAQ experiment was converted to sample ion intensities using the protein ratio data from Scaffold. Ion intensity were converted to protein concentrations using the known concentration of human homologs MB, HBB, CKM, MYL2, GAPDH, FABP3, and LDHB . Average supersaturation scores were calculated for proteins whose abundance changed (>±1.3 fold) or didn’t change (<±1.1) between an experimental condition and the time matched control. A one way ANOVA with post hoc Bonferroni correction Statistics was conducted in Graphpad Prism (v5). Data is plotted as mean +/- standard error of the mean.
Arsenite-inducible RNA-associated protein
Analysis of variance
- Atf3 :
Activating transcription factor 3
Blood urea nitrogen
Consensus heat stress response
Cyclic AMP response element-binding
Database of Annotation, Visualization and Integrated Discovery
Differentially expressed genes
Functional annotation clusters
False discovery rate
Hydroxyacyl-coenzyme A dehydrogenase
- Hsf1 :
Heat shock factor 1
Heat shock protein
Isobaric tag for relative and absolute quantitation
Kyoto Encyclopedia of Genes and Genomes
Mitogen-associated protein kinase
National Center for Biotechnology Information
Probe logarithmic intensity error
Regulation of apoptosis
Systemic inflammatory response syndrome
Becker JA, Stewart LK: Heat-related illness. Am Fam Physician. 2011, 83 (11): 1325-1330.
Bouchama A, Knochel JP: Heat stroke. N Engl J Med. 2002, 346 (25): 1978-1988. 10.1056/NEJMra011089.
Epstein Y, Roberts WO: The pathopysiology of heat stroke: an integrative view of the final common pathway. Scand J Med Sci Sports. 2011, 21 (6): 742-748. 10.1111/j.1600-0838.2011.01333.x.
Wenger RH: Cellular adaptation to hypoxia: O2-sensing protein hydroxylases, hypoxia-inducible transcription factors, and O2-regulated gene expression. FASEB J. 2002, 16 (10): 1151-1162. 10.1096/fj.01-0944rev.
Winkenwerder M, Sawka MN: Disorders due to heat and cold. Cecil Medicine. Edited by: Goldman L, Ausiello D. 2007, Philadelphia: Elsevier Science, 763-767. 23
Lanner JT, Georgiou DK, Dagnino-Acosta A, Ainbinder A, Cheng Q, Joshi AD, Chen ZW, Yarotskyy V, Oakes JM, Lee CS, Monroe TO, Sabtillan A, Dong K, Goodyear L, Ismailov II, Rodney GG, Dirksen RT, Hamilton SL: AICAR prevents heat-induced sudden death in RyR1 mutant mice independent of AMPK activation. Nat Med. 2012, 18 (2): 244-251. 10.1038/nm.2598.
Leon LR, Helwig BG: Heat stroke: role of the systemic inflammatory response. J Appl Physiol Respir Environ Exerc Physiol. 2010, 109 (6): 1980-1988.
Leon LR, Helwig BG: Role of endotoxin and cytokines in the systemic inflammatory response to heat injury. Front Biosci. 2010, 2: 916-938.
Grogan H, Hopkins PM: Heat stroke: implications for critical care and anaesthesia. Br J Anaesth. 2002, 88 (5): 700-707. 10.1093/bja/88.5.700.
Kregel KC: Heat shock proteins: modifying factors in physiological stress responses and acquired thermotolerance. J Appl Physiol Respir Environ Exerc Physiol. 2002, 92 (5): 2177-2186.
Moran DS, Eli-Berchoer L, Heled Y, Mendel L, Schocina M, Horowitz M: Heat intolerance: does gene transcription contribute?. J Appl Physiol Respir Environ Exerc Physiol. 2006, 100 (4): 1370-1376.
Parsell DA, Lindquist S: The function of heat-shock proteins in stress tolerance: degradation and reactivation of damaged proteins. Annu Rev Genet. 1993, 27: 437-496. 10.1146/annurev.ge.27.120193.002253.
Parsell DA, Taulien J, Lindquist S: The role of heat-shock proteins in thermotolerance. Philos Trans R Soc Lond B Biol Sci. 1993, 339 (1289): 279-285. 10.1098/rstb.1993.0026. discussion 285–276
Zhang HJ, Drake VJ, Morrison JP, Oberley LW, Kregel KC: Selected contribution: differential expression of stress-related genes with aging and hyperthermia. J Appl Physiol Respir Environ Exerc Physiol. 2002, 92 (4): 1762-1769. discussion 1749
Sonna LA, Wenger CB, Flinn S, Sheldon HK, Sawka MN, Lilly CM: Exertional heat injury and gene expression changes: a DNA microarray analysis study. J Appl Physiol Respir Environ Exerc Physiol. 2004, 96 (5): 1943-1953.
Sonna LA, Glueck SB, Jeunemaitre X: Exercise, genetics, and blood pressure: focus on "physical exercise and blood pressure with reference to the angiotensinogen M235T polymorphism" and on "angiotensinogen M235T polymorphism associates with exercise hemodynamics in postmenopausal women". Physiol Genomics. 2002, 10 (2): 45-47.
Sonna LA, Gaffin SL, Pratt RE, Cullivan ML, Angel KC, Lilly CM: Effect of acute heat shock on gene expression by human peripheral blood mononuclear cells. J Appl Physiol Respir Environ Exerc Physiol. 2002, 92 (5): 2208-2220.
Sonna LA, Fujita J, Gaffin SL, Lilly CM: Invited review: effects of heat and cold stress on mammalian gene expression. J Appl Physiol Respir Environ Exerc Physiol. 2002, 92 (4): 1725-1742.
Horowitz M, Eli-Berchoer L, Wapinski I, Friedman N, Kodesh E: Stress-related genomic responses during the course of heat acclimation and its association with ischemic-reperfusion cross-tolerance. J Appl Physiol Respir Environ Exerc Physiol. 2004, 97 (4): 1496-1507.
Rakesh V, Stallings JD, Helwig BG, Leon LR, Jackson DA, Reifman J: A 3-D mathematical model to identify organ-specific risks in rats during thermal stress. J Appl Physiol. 2013, 5 (12): 1822-1837.
Horowitz M: Genomics and proteomics of heat acclimation. Front Biosci. 2010, 2: 1068-1080.
Islam A, Deuster PA, Devaney JM, Ghimbovschi S, Chen Y: An exploration of heat tolerance in mice utilizing mRNA and microRNA expression analysis. PLoS One. 2013, 8 (8): e72258-10.1371/journal.pone.0072258.
Ha K, Coulombe-Huntington J, Majewski J: Comparison of affymetrix gene array with the exon array shows potential application for detection of transcript isoform variation. BMC Genomics. 2009, 10: 519-10.1186/1471-2164-10-519.
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4 (5): 3-10.1186/gb-2003-4-5-p3.
Huang J, Wang H, Xie X, Gao H, Guo G: Developmental changes in DNA methylation of pollen mother cells of David lily during meiotic prophase I. Mol Biol. 2010, 44 (5): 853-858.
Leverette RD, Andrews MT, Maxwell ES: Mouse U14 snRNA is a processed intron of the cognate hsc70 heat shock pre-messenger RNA. Cell. 1992, 71 (7): 1215-1221. 10.1016/S0092-8674(05)80069-8.
Chen MS, Goswami PC, Laszlo A: Differential accumulation of U14 snoRNA and hsc70 mRNA in Chinese hamster cells after exposure to various stress conditions. Cell Stress Chaperones. 2002, 7 (1): 65-72. 10.1379/1466-1268(2002)007<0065:DAOUSA>2.0.CO;2.
Brocker C, Thompson DC, Vasiliou V: The role of hyperosmotic stress in inflammation and disease. Biomol Concepts. 2012, 3 (4): 345-364.
Lotz GP, Lin H, Harst A, Obermann WM: Aha1 binds to the middle domain of Hsp90, contributes to client protein activation, and stimulates the ATPase activity of the molecular chaperone. J Biol Chem. 2003, 278 (19): 17228-17235. 10.1074/jbc.M212761200.
Gerstner JR, Landry CF: The zinc-binding protein chordc1 undergoes complex diurnal changes in mRNA expression during mouse brain development. Neurochem Res. 2007, 32 (2): 241-250. 10.1007/s11064-006-9271-z.
Gano JJ, Simon JA: A proteomic investigation of ligand-dependent HSP90 complexes reveals CHORDC1 as a novel ADP-dependent HSP90-interacting protein. Mol Cell Proteomics. 2010, 9 (2): 255-270. 10.1074/mcp.M900261-MCP200.
Pavithra L, Mukherjee S, Sreenath K, Kar S, Sakaguchi K, Roy S, Chattopadhyay S: SMAR1 forms a ternary complex with p53-MDM2 and negatively regulates p53-mediated transcription. J Mol Biol. 2009, 388 (4): 691-702. 10.1016/j.jmb.2009.03.033.
Pavithra L, Singh S, Sreenath K, Chattopadhyay S: Tumor suppressor SMAR1 downregulates Cytokeratin 8 expression by displacing p53 from its cognate site. Int J Biochem Cell Biol. 2009, 41 (4): 862-871. 10.1016/j.biocel.2008.08.038.
Pavithra L, Sreenath K, Singh S, Chattopadhyay S: Heat-shock protein 70 binds to a novel sequence in 5' UTR of tumor suppressor SMAR1 and regulates its mRNA stability upon Prostaglandin A2 treatment. FEBS Lett. 2010, 584 (6): 1187-1192. 10.1016/j.febslet.2010.02.025.
Schultz CR, Golembieski WA, King DA, Brown SL, Brodie C, Rempel SA: Inhibition of HSP27 alone or in combination with pAKT inhibition as therapeutic approaches to target SPARC-induced glioma cell survival. Mol Cancer. 2012, 11: 20-10.1186/1476-4598-11-20.
Oehrl W, Cotsiki M, Panayotou G: Differential regulation of M3/6 (DUSP8) signaling complexes in response to arsenite-induced oxidative stress. Cell Signal. 2013, 25 (2): 429-438. 10.1016/j.cellsig.2012.11.010.
Bauer H, Willert J, Koschorz B, Herrmann BG: The t complex-encoded GTPase-activating protein Tagap1 acts as a transmission ratio distorter in mice. Nat Genet. 2005, 37 (9): 969-973. 10.1038/ng1617.
Kuchta K, Knizewski L, Wyrwicz LS, Rychlewski L, Ginalski K: Comprehensive classification of nucleotidyltransferase fold proteins: identification of novel families and their representatives in human. Nucleic Acids Res. 2009, 37 (22): 7701-7714. 10.1093/nar/gkp854.
Xu X, Zeng W, Popov S, Berman DM, Davignon I, Yu K, Yowe D, Offermanns S, Muallem S, Wilkie TM: RGS proteins determine signaling specificity of Gq-coupled receptors. J Biol Chem. 1999, 274 (6): 3549-3556. 10.1074/jbc.274.6.3549.
Rother S, Clausing E, Kieser A, Strasser K: Swt1, a novel yeast protein, functions in transcription. J Biol Chem. 2006, 281 (48): 36518-36525. 10.1074/jbc.M607510200.
Zhao W, Wang L, Zhang M, Wang P, Yuan C, Qi J, Meng H, Gao C: Tripartite motif-containing protein 38 negatively regulates TLR3/4- and RIG-I-mediated IFN-beta production and antiviral response by targeting NAP1. J Immunol. 2012, 188 (11): 5311-5318. 10.4049/jimmunol.1103506.
Bernales S, McDonald KL, Walter P: Autophagy counterbalances endoplasmic reticulum expansion during the unfolded protein response. PLoS Biol. 2006, 4 (12): e423-10.1371/journal.pbio.0040423.
Fujita E, Kouroku Y, Isoai A, Kumagai H, Misutani A, Matsuda C, Hayashi YK, Momoi T: Two endoplasmic reticulum-associated degradation (ERAD) systems for the novel variant of the mutant dysferlin: ubiquitin/proteasome ERAD(I) and autophagy/lysosome ERAD(II). Hum Mol Genet. 2007, 16 (6): 618-629.
Hoyer-Hansen M, Jaattela M: Connecting endoplasmic reticulum stress to autophagy by unfolded protein response and calcium. Cell Death Differ. 2007, 14 (9): 1576-1582. 10.1038/sj.cdd.4402200.
Kamimoto T, Shoji S, Hidvegi T, Mizushima N, Umebayashi K, Perlmutter DH, Yoshimori T: Intracellular inclusions containing mutant alpha1-antitrypsin Z are propagated in the absence of autophagic activity. J Biol Chem. 2006, 281 (7): 4467-4476. 10.1074/jbc.M509409200.
Kouroku Y, Fujita E, Tanida I, Ueno T, Isoai A, Kumagai H, Ogawa S, Kaufman RJ, Kominami E, Momoi T: ER stress (PERK/eIF2alpha phosphorylation) mediates the polyglutamine-induced LC3 conversion, an essential step for autophagy formation. Cell Death Differ. 2007, 14 (2): 230-239. 10.1038/sj.cdd.4401984.
Ogata M, Hino S, Saito A, Morikawa K, Kondo S, Kanemoto S, Murakami T, Taniguchi M, Tanii I, Yoshinaga K, Shiosaka S, Hammarback JA, Urano F, Imaizumi K: Autophagy is activated for cell survival after endoplasmic reticulum stress. Mol Cell Biol. 2006, 26 (24): 9220-9231. 10.1128/MCB.01453-06.
Yorimitsu T, Nair U, Yang Z, Klionsky DJ: Endoplasmic reticulum stress triggers autophagy. J Biol Chem. 2006, 281 (40): 30299-30304. 10.1074/jbc.M607007200.
Chakrabarti A, Chen AW, Varner JD: A review of the mammalian unfolded protein response. Biotechnol Bioeng. 2011, 108 (12): 2777-2793. 10.1002/bit.23282.
Morimoto RI: The heat shock response: systems biology of proteotoxic stress in aging and disease. Cold Spring Harb Symp Quant Biol. 2011, 76: 91-99. 10.1101/sqb.2012.76.010637.
Morimoto RI, Driessen AJ, Hegde RS, Langer T: The life of proteins: the good, the mostly good and the ugly. Nat Struct Mol Biol. 2011, 18 (1): 1-4. 10.1038/nsmb0111-1.
Fonseca SG, Burcin M, Gromada J, Urano F: Endoplasmic reticulum stress in beta-cells and development of diabetes. Curr Opin Pharmacol. 2009, 9 (6): 763-770. 10.1016/j.coph.2009.07.003.
Glembotski CC: Endoplasmic reticulum stress in the heart. Circ Res. 2007, 101 (10): 975-984. 10.1161/CIRCRESAHA.107.161273.
Seigneuric R, Mjahed H, Gobbo J, Joly AL, Berthenet K, Shirley S, Garrido C: Heat shock proteins as danger signals for cancer detection. Front Biol. 2011, 1: 37-
Kampinga HH, Hageman J, Vos MJ, Kubota H, Tanguay RM, Bruford EA, Cheetham ME, Chen B, Hightower LE: Guidelines for the nomenclature of the human heat shock proteins. Cell Stress Chaperones. 2009, 14 (1): 105-111. 10.1007/s12192-008-0068-7.
Macario AJ, Conway de Macario E: The pathology of cellular anti-stress mechanisms: a new frontier. Stress. 2004, 7 (4): 243-249. 10.1080/10253890400019706.
Gallucci S, Matzinger P: Danger signals: SOS to the immune system. Curr Opin Immunol. 2001, 13 (1): 114-119. 10.1016/S0952-7915(00)00191-6.
Park HS, Lee JS, Huh SH, Seo JS, Choi EJ: Hsp72 functions as a natural inhibitory protein of c-Jun N-terminal kinase. EMBO J. 2001, 20 (3): 446-456. 10.1093/emboj/20.3.446.
Cvoro A, Korac A, Matic G: Intracellular localization of constitutive and inducible heat shock protein 70 in rat liver after in vivo heat stress. Mol Cell Biochem. 2004, 265 (1–2): 27-35.
Garrido C, Brunet M, Didelot C, Zermati Y, Schmitt E, Kroemer G: Heat shock proteins 27 and 70: anti-apoptotic proteins with tumorigenic properties. Cell Cycle. 2006, 5 (22): 2592-2601. 10.4161/cc.5.22.3448.
Didelot C, Lanneau D, Brunet M, Joly AL, De Thonel A, Chiosis G, Garrido C: Anti-cancer therapeutic approaches based on intracellular and extracellular heat shock proteins. Curr Med Chem. 2007, 14 (27): 2839-2847. 10.2174/092986707782360079.
Dehbi M, Baturcam E, Eldali A, Ahmed M, Kwaasi A, Chishti MA, Bouchama A: Hsp-72, a candidate prognostic indicator of heatstroke. Cell Stress Chaperones. 2010, 15 (5): 593-603. 10.1007/s12192-010-0172-3.
Beere HM: Death versus survival: functional interaction between the apoptotic and stress-inducible heat shock protein pathways. J Clin Invest. 2005, 115 (10): 2633-2639. 10.1172/JCI26471.
Johnson JD, Fleshner M: Releasing signals, secretory pathways, and immune function of endogenous extracellular heat shock protein 72. J Leukoc Biol. 2006, 79 (3): 425-434.
Wang ZZ, Wang CL, Wu TC, Pan HN, Wang SK, Jiang JD: Autoantibody response to heat shock protein 70 in patients with heatstroke. Am J Med. 2001, 111 (8): 654-657. 10.1016/S0002-9343(01)00974-3.
Galloway E, Shin T, Huber N, Eismann T, Kuboki S, Schuster R, Blanchard J, Wong HR, Lentsch AB: Activation of hepatocytes by extracellular heat shock protein 72. Am J Physiol Cell Physiol. 2008, 295 (2): C514-520. 10.1152/ajpcell.00032.2008.
Salari S, Seibert T, Chen YX, Hu T, Shi C, Zhao X, Cuerrier CM, Raizman JE, O'Brien ER: Extracellular HSP27 acts as a signaling molecule to activate NF-kappaB in macrophages. Cell Stress Chaperones. 2013, 18 (1): 53-63. 10.1007/s12192-012-0356-0.
Thompson MR, Xu D, Williams BR: ATF3 transcription factor and its emerging roles in immunity and cancer. J Mol Med. 2009, 87 (11): 1053-1060. 10.1007/s00109-009-0520-x.
Lu D, Wolfgang CD, Hai T: Activating transcription factor 3, a stress-inducible gene, suppresses Ras-stimulated tumorigenesis. J Biol Chem. 2006, 281 (15): 10473-10481. 10.1074/jbc.M509278200.
Gilchrist M, Thorsson V, Li B, Rust AG, Korb M, Roach JC, Kennedy K, Hai T, Bolouri H, Aderem A: Systems biology approaches identify ATF3 as a negative regulator of Toll-like receptor 4. Nature. 2006, 441 (7090): 173-178. 10.1038/nature04768.
Hartman MG, Lu D, Kim ML, Kociba GJ, Shukri T, Buteau J, Wang X, Frankel WL, Guttridge D, Prentki M, Grey ST, Ron D, Hai T: Role for activating transcription factor 3 in stress-induced beta-cell apoptosis. Mol Cell Biol. 2004, 24 (13): 5721-5732. 10.1128/MCB.24.13.5721-5732.2004.
Horowitz M: Heat acclimation and cross-tolerance against novel stressors: genomic-physiological linkage. Prog Brain Res. 2007, 162: 373-392.
Horowitz M, Robinson SD: Heat shock proteins and the heat shock response during hyperthermia and its modulation by altered physiological conditions. Prog Brain Res. 2007, 162: 433-446.
Karper JC, de Jager SC, Ewing MM, de Vries MR, Bot I, van Santbrink PJ, Redeker A, Mallat Z, Binder CJ, Arens R, Jukema JW, Kuiper J, Quax PH: An unexpected intriguing effect of toll-like receptor regulator RP105 (CD180) on atherosclerosis formation with alterations on B-cell activation. Arterioscler Thromb Vasc Biol. 2013, 33 (12): 2810-2817. 10.1161/ATVBAHA.113.301882.
Karper JC, Ewing MM, de Vries MR, de Jager SC, Peters EA, de Boer HC, van Zonneveld AJ, Kuiper J, Huizinga EG, Brondijk TH, Kuiper J, Quax PH: TLR accessory molecule RP105 (CD180) is involved in post-interventional vascular remodeling and soluble RP105 modulates neointima formation. PLoS One. 2013, 8 (7): e67923-10.1371/journal.pone.0067923.
Balch WE, Morimoto RI, Dillin A, Kelly JW: Adapting proteostasis for disease intervention. Science. 2008, 319 (5865): 916-919. 10.1126/science.1141448.
Dobson CM: Protein folding and misfolding. Nature. 2003, 426 (6968): 884-890. 10.1038/nature02261.
Querfurth HW, LaFerla FM: Alzheimer's disease. N Engl J Med. 2010, 362 (4): 329-344. 10.1056/NEJMra0909142.
Selkoe DJ: Alzheimer's disease. Cold Spring Harb Perspect Biol. 2011, 3: 7-
Ciryam P, Tartaglia GG, Morimoto RI, Dobson CM, Vendruscolo M: Widespread aggregation and neurodegenerative diseases are associated with supersaturated proteins. Cell Rep. 2013, 5 (3): 781-790. 10.1016/j.celrep.2013.09.043.
Pan YX, Lin L, Ren AJ, Pan XJ, Chen H, Tang CS, Yuan WJ: HSP70 and GRP78 induced by endothelin-1 pretreatment enhance tolerance to hypoxia in cultured neonatal rat cardiomyocytes. J Cardiovasc Pharmacol. 2004, 44 (Suppl 1): S117-120.
Pan YX, Ren AJ, Zheng J, Rong WF, Chen H, Yan XH, Wu C, Yuan WJ, Lin L: Delayed cytoprotection induced by hypoxic preconditioning in cultured neonatal rat cardiomyocytes: role of GRP78. Life Sci. 2007, 81 (13): 1042-1049. 10.1016/j.lfs.2007.08.015.
Sun FC, Wei S, Li CW, Chang YS, Chao CC, Lai YK: Localization of GRP78 to mitochondria under the unfolded protein response. Biochem J. 2006, 396 (1): 31-39. 10.1042/BJ20051916.
Davies B, Morris T: Physiological parameters in laboratory animals and humans. Pharm Res. 1993, 10 (7): 1093-1095. 10.1023/A:1018943613122.
Qu Y, He F, Chen Y: Different effects of the probe summarization algorithms PLIER and RMA on high-level analysis of Affymetrix exon arrays. BMC Bioinformatics. 2010, 11: 211-10.1186/1471-2105-11-211.
Bourgon R, Gentleman R, Huber W: Independent filtering increases detection power for high-throughput experiments. Proc Natl Acad Sci U S A. 2010, 107 (21): 9546-9551. 10.1073/pnas.0914005107.
Hochberg Y, Benjamini Y: More powerful procedures for multiple significance testing. Stat Med. 1990, 9 (7): 811-818. 10.1002/sim.4780090710.
Reiner A, Yekutieli D, Benjamini Y: Identifying differentially expressed genes using false discovery rate controlling procedures. Bioinformatics. 2003, 19 (3): 368-375. 10.1093/bioinformatics/btf877.
Eng JK, McCormack AL, Yates JR: An approach to correlate tandem mass spectral data of peptides with amino acid sequences in a protein database. J Am Soc Mass Spectrom. 1994, 5 (11): 976-989. 10.1016/1044-0305(94)80016-2.
Nesvizhskii AI, Keller A, Kolker E, Aebersold R: A statistical model for identifying proteins by tandem mass spectrometry. Anal Chem. 2003, 75 (17): 4646-4658. 10.1021/ac0341261.
Shadforth IP, Dunkley TP, Lilley KS, Bessant C: i-Tracker: for quantitative proteomics using iTRAQ. BMC Genomics. 2005, 6: 145-10.1186/1471-2164-6-145.
Ahn JM, Sung HJ, Yoon YH, Kim BG, Yang WS, Lee C, Park HM, Kim BJ, Kim BG, Lee SY, An HJ, Cho JY: Integrated glycoproteomics demonstrates fucosylated serum paraoxonase 1 alterations in small cell lung cancer. Mol Cell Proteomics. 2014, 13 (1): 30-48. 10.1074/mcp.M113.028621.
Silva JC, Gorenstein MV, Li GZ, Vissers JP, Geromanos SJ: Absolute quantification of proteins by LCMSE: a virtue of parallel MS acquisition. Mol Cell Proteomics. 2006, 5 (1): 144-156.
Scholten A, Heck JR: Determining protein concentrations of the human ventricular proteome. Heart Proteomics: Methods and Protocols. Volume 1005. Edited by: Vivanco F. 2013, New York: Springer Science+Business Media, 11-24.
We would like to thank Dr. Roy Vigneulle (Military Operational Medicine Research Program), CAPT Carroll D. Forcino (Director, Military Operational Medicine Research Program), COL Richard P. Duncan (former Commander, USACEHR), and COL Thomas C. Timmes (current Commander, USACEHR) for their programmatic support, encouragement, and insightful discussion. We also thank J. Ward, S. Dineen, M. Blaha, and R. Duran for technical support with the rat heat-stress experiments. We thank Gian G. Tartaglia and Michele Vendruscolo (University of Cambridge and Northwestern University) for batch calculation of the zagg and zSCagg used to calculate the σf scores reported in this paper. This research was supported in part by an appointment to the Research Participation Program at the US Army Center for Environmental Health Research administered by the Oak Ridge Institute for Science and Education through an interagency agreement between the U.S. Department of Energy and MRMC and USACEHR. The research was supported by the Military Operational Medicine Research Program, US Army Medical Research and Materiel Command, Fort Detrick, Maryland.
Research was conducted in compliance with the Animal Welfare Act, and other Federal statutes and regulations relating to animals and experiments involving animals and adheres to principles stated in the “Guide for Care and Use of Laboratory Animals” (NRC 2011) as prepared by the Committee on Care and Use of Laboratory Animals of the Institute of Laboratory Animal Resources, National Research Council in facilities that are fully accredited by the Association for Assessment and Accreditation of Laboratory Animal Care, International. The authors have declared that they have no competing interests. The views, opinions, assertions, and/or findings contained herein are those of the authors and should not be construed as official US Department of Defense or Department of the Army position, policy, or decision, unless so designated by other official documentation. Citations of commercial organizations or trade names in this report do not constitute an official Department of the Army endorsement or approval of the products or services of these organizations. This paper has been approved for public release with unlimited distribution.
The initiation and conception of this study involved JDS, VR, BGW, JAL, LRL, JR, and DAJ. LRL and BGW were responsible for all animal exposures and model development. CEB performed all protein and RNA preparation for transcriptomics and proteomics, and performed all transcriptomic experiments. JDS performed all gene array and pathway analyses. WED performed all mass spectrometry, and DLI analyzed the mass spectrometry data. JAL performed the aggregation score analysis. JDS (overall) was responsible for writing the manuscript, which was edited by VR, DLI, LRL, and JR. All authors reviewed the final manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 4: Figure S1: Figure illustrating the histopathological evidence of kidney and liver injury that occurred in one rat 24 hours after heat stress. (DOCX 1 MB)
Additional file 6: Table S5: Table of all differentially expressed genes in the heart, liver, lung, and kidney at Tc,Max, 24 hours, and 48 hours. (XLSX 214 KB)
Additional file 8: Table S7: Table listing DEGs mapping to KEGG pathways and genes of interest identified by self-organizing maps for both the heart (48 hours) and liver (24 hours). (XLSX 97 KB)
Additional file 9: Table S8: Table listing the functional annotation clusters identified by DAVID analysis. (XLSX 43 KB)
Additional file 10: Figure S2: Figure showing the self-organizing maps identifying gene nodes enriched in the liver at 24 hours. (DOCX 458 KB)
Additional file 11: Table S9: Table of differentially expressed proteins (DEP) identified by quantification of normalized intensities by iTRAQ proteomics. (XLSX 2 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Stallings, J.D., Ippolito, D.L., Rakesh, V. et al. Patterns of gene expression associated with recovery and injury in heat-stressed rats. BMC Genomics 15, 1058 (2014). https://doi.org/10.1186/1471-2164-15-1058
- Heat stress
- Systems biology
- Protein aggregation