Journal of Applied Physiology Information on EB 2010
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Appl Physiol 96: 1943-1953, 2004. First published February 20, 2004; doi:10.1152/japplphysiol.00886.2003
8750-7587/04 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow All Versions of this Article:
96/5/1943    most recent
00886.2003v1
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (23)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Sonna, L. A.
Right arrow Articles by Lilly, C. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Sonna, L. A.
Right arrow Articles by Lilly, C. M.

Exertional heat injury and gene expression changes: a DNA microarray analysis study

Larry A. Sonna,1,2 C. Bruce Wenger,{dagger} Scott Flinn,3 Holly K. Sheldon,2 Michael N. Sawka,1 and Craig M. Lilly2

1Thermal and Mountain Medicine Division, United States Army Research Institute of Environmental Medicine, Natick 01760; 2Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital/Harvard Medical School, Boston, Massachusetts 02115; and 3Directorate for Primary Care, Naval Medical Center San Diego, San Diego, California 92134

Submitted 19 August 2003 ; accepted in final form 9 January 2004


    ABSTRACT
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
This study examined gene expression changes associated with exertional heat injury (EHI) in vivo and compared these changes to in vitro heat shock responses previously reported by our laboratory. Peripheral blood mononuclear cell (PBMC) RNA was obtained from four male Marine recruits (ages 17-19 yr) who presented with symptoms consistent with EHI, core temperatures ranging from 39.3 to 42.5°C, and elevations in serum enzymes such as creatine kinase. Controls were age- and gender-matched Marines from whom samples were obtained before and several days after an intense field-training exercise in the heat ("The Crucible"). Expression analysis was performed on Affymetrix arrays (containing ~12,600 sequences) from pooled samples obtained at three times for EHI group (at presentation, 2-3 h after cooling, and 1-2 days later) and compared with control values (average signals from two chips representing pre- and post-Crucible samples). After post hoc filtering, the analysis identified 361 transcripts that had twofold or greater increases in expression at one or more of the time points assayed and 331 transcripts that had twofold or greater decreases in expression. The affected transcripts included sequences previously shown to be heat-shock responsive in PBMCs in vitro (including both heat shock proteins and non-heat shock proteins), a number of sequences whose changes in expression had not previously been noted as a result of in vitro heat shock in PBMCs (including several interferon-induced sequences), and several nonspecific stress response genes (including ubiquitin C and dual-specificity phosphatase-1). We conclude that EHI produces a broad stress response that is detectable in PBMCs and that heat stress per se can only account for some of the observed changes in transcript expression. The molecular evidence from these patients is thus consistent with the hypothesis that EHI can result from cumulative effects of multiple adverse interacting stimuli.

heat stroke; exercise; peripheral blood mononuclear cells; genomics


EXERTIONAL HEAT ILLNESS encompasses a spectrum of clinical conditions, of which the most severe are exertional heat injury (EHI) and exertional heatstroke (EHS). EHI is defined as severe illness associated with high body temperatures in the setting of physical exertion, typically involving exhaustion or collapse, which is accompanied by evidence of tissue and organ damage (11). EHS is a more severe form of heat illness defined clinically by the presence of significant concurrent central nervous system dysfunction (11). EHI and EHS account for many medical emergencies each year in athlete, worker, and military populations (25) and have recently received considerable attention because of a number of high-profile deaths among athletes.

The pathophysiological mechanisms that produce EHI and EHS are under active investigation (2). Experimental approaches have included investigations of animal mortality (10, 14), tissue damage (13, 17), and cellular factors [such as heat shock protein (HSP) expression] that affect tissue susceptibility to thermal injury (9, 16, 20). Animal models and observations in humans have identified mechanisms such as heat-induced translocation of lipopolysaccharide from the bowel lumen into the circulation (2, 12, 17, 22), reticuloendothelial system deactivation (4), and heat-induced changes in pro- and antiinflammatory cytokines (2) as contributing to the pathogenesis of EHI and EHS. Additionally, human studies have examined the possibility that an inflammatory response induced by previous muscle injury leads to an accentuated hyperthermic response during exercise-induced heat loads (19). The picture that emerges from these mechanistic studies is that EHI and EHS share a common pathophysiological basis (3, 11) and represent systemic manifestations of derangements that occur at the level of cells and tissues, in which the immune system is highly involved (2). Identification of cellular pathways involved in EHI and EHS is a logical and important extension of this body of work. Until recently, this has been difficult because of limitations in our ability to simultaneously assay large numbers of mediator molecules in the context of thermal stress.

At the cellular level, it is generally accepted that thermal stress leads to increases in expression of HSPs and that the expression of these proteins correlates closely with the acquisition of thermotolerance (16, 18, 21). However, it is also increasingly apparent that the cellular response to heat shock involves more than HSPs (16, 28). For example, a recent DNA microarray study of rats exposed to whole body hyperthermia found substantial changes in liver expression of both HSPs and genes belonging to other functional classes (32). Similarly, isolated human peripheral blood mononuclear cells (PBMCs) subjected to in vitro heat shock demonstrated extensive changes in expression of both HSPs and genes not traditionally considered as HSPs (29).

Knowledge of the molecular processes involved in EHI and EHS in humans has lagged behind our understanding of systemic processes, in part because it is difficult to obtain samples of human tissues thought to be of high potential interest, such as intestine, liver, skeletal muscle, and heart. Because leukocytes are easily obtained from human volunteers and show changes in HSP expression after heat shock in vitro (6, 7, 23, 24, 29) and after exercise in vivo (6-8, 23, 24), and because blood is exposed to all body compartments and thus all body temperatures (which are typically substantially elevated during EHI and EHS), PBMCs likely represent an informative cell type to identify molecular changes in EHI and EHS patients. Accordingly, we thought it likely that individuals suffering from EHI would demonstrate substantial changes in PBMC gene expression and that at least some of these changes would be similar to those found in human in vitro models of heat shock. This study, therefore, examined in vivo changes in PBMC gene expression associated with EHI and compared them to in vitro heat shock responses in PBMCs that were previously reported by our laboratory (29). We hypothesized that gene expression would be similar to heat shock but that many other genes associated with known systemic processes, particularly those related to inflammatory responses, would be expressed with EHI.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Subjects and controls. The four EHI cases were male Marine recruits undergoing basic training at Parris Island, who presented to the field corpsman (medic) with evidence of heat illness and were evacuated to the Branch Medical Clinic Parris Island. The criteria for EHI included hyperthermia from physical exercise, with signs and symptoms such as collapse during exertion in the heat, malaise, disorientation, or ataxia, accompanied by laboratory evidence of organ dysfunction [such as elevated serum creatine kinase (CK), transaminases, or lactate dehydrogenase], a level of consciousness that returned to normal after cooling, and the absence of seizure activity. All subjects had active cooling started in the field (icing). On arrival to the clinic, they underwent further active cooling with ice sheets, typically followed by showering. In addition, all subjects received intravenous hydration. Blood samples were drawn on presentation, after cooling [mean ± SD, 3:02 ± 0:24 (h:min) after presentation], and again on a follow-up visit 1-2 days later. Several blood samples were drawn on case 4 in the days after the initial scheduled follow-up. All cases were presented between July and September 2000. All four subjects had been in basic training at Parris Island long enough (range, 9-20 days from start of training to injury) to have become heat acclimatized.

Control subjects were three Marine recruits matched for age, gender, and ethnic group, who volunteered to donate blood samples in August 2000, several days before and several days after an intense field-training exercise ("The Crucible," the capstone event of US Marine Basic Training, which represents 54 h of vigorous field training with average energy expenditures of >6,000 kcal/day for men). Because this event was preceded by weeks of basic combat training in the summer heat, it is assumed that the controls were also heat acclimatized.

All Marine Corps recruits are medically screened before enlistment for an extensive battery of disqualifying conditions, including a variety of chronic liver diseases, infectious diseases such as chronic hepatitis, autoimmune conditions such as systemic lupus erythematosus that can produce renal and hepatic dysfunction, autoimmune conditions capable of producing rhabdomyolysis, alcohol dependence, and renal disease. Also disqualifying for enlistment are a known predisposition to heat illness, a prior history of malignant hyperthermia, recurrent episodes of heat injury requiring medical attention, and evidence of residual injury from a prior heat illness.

Blood collection. Blood was collected from all individuals in cell-preparation tubes (Fisher Scientific). Mononuclear cells were separated per the manufacturer's instructions, suspended in ~1 ml of RNALater (Ambion), and frozen at -80°C until RNA processing occurred.

RNA extraction. Samples were thawed, and the cells were pelleted by centrifugation. RNA was extracted by using the RNeasy Mini Kits (Qiagen). RNA yield was estimated by absorbance spectrophotometry, and the samples were stored at -80°C.

DNA microarray analysis. Due to limiting sample volumes and RNA yields, equal amounts of RNA from the samples obtained from the four cases were pooled. Similarly, equal amounts of RNA from the samples obtained from the three control subjects were pooled. DNA microarray analysis was performed by using Affymetrix U95Av2 gene chips as described previously (27, 29). Analyses using Affymetrix arrays (containing ~12,600 sequences) were performed from pooled samples obtained at each of the three time points for EHI group (1 chip each for samples obtained at presentation, after cooling, and at follow-up) and compared with control values (1 chip each for samples obtained before and after the Crucible) as described in Statistical analysis.

Expression measurement by RT-PCR. RT-PCR on the pooled samples was performed using standard techniques, as described previously (27, 29). The primers used for these reactions are listed (see Table 4).


View this table:
[in this window]
[in a new window]
 
Table 4. PCR Primers used

 

Statistical analysis. The output signal from the Affymetrix U95Av2 array was preprocessed using MAS 5.0 software (Affymetrix). For each sequence on the array, the mean and standard deviation of the software-reported signals were calculated from the two pooled control chips. From these values, 98.33% population intervals were computed by adding or subtracting 2.39 standard deviations from the mean. Corresponding sequences from the heat-injury cases were considered to have a statistically significant deviation from the controls if their signals fell outside of these population intervals. Population intervals of 98.33% were chosen rather than 95% intervals to correct for the three multiple comparisons being performed; i.e., a P value of 0.05/3 = 0.0167 was considered statistically significant. Other statistical analyses were performed as noted throughout the manuscript using SigmaStat 2.03 for Windows and taking a P value of <=0.05 as statistically significant.

Where noted in this manuscript, post hoc filtering of significantly affected sequences was performed by two criteria. First, sequences had to show a twofold or greater change in expression. Second, sequences had to be detected by the chip-reading software as "present" or "marginal" in both control chips (for downregulated sequences) or in the chip(s) corresponding to the time point(s) of interest (for upregulated sequences). This post hoc filtering was performed to maximize the comparability of the results of this experiment with the analysis performed in our laboratory's previous in vitro experiment (29) involving five sets of normal human PBMCs subjected to heat shock in serum-free media. In this previous study, preprocessing was performed using MAS 4.0 software, and the post hoc filters used were a geometric mean change in expression of twofold or greater and a present or marginal call in five of five experiments in the controls (for downregulated sequences) or in the heat-shocked cells (for upregulated sequences).


    RESULTS
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
EHI and control subjects. Table 1 provides the demographic characteristics of the EHI and control subjects. Tables 2 and 3 provide the clinical and laboratory characteristics of the EHI patients. Values outside the range of normal for the clinical laboratory are highlighted in bold font in Table 3.


View this table:
[in this window]
[in a new window]
 
Table 1. Demographic characteristics of cases and controls

 

View this table:
[in this window]
[in a new window]
 
Table 2. Clinical characteristics of subjects experiencing exertional heat illness

 

View this table:
[in this window]
[in a new window]
 
Table 3. Clinical laboratory characteristics of subjects experiencing exertional heat illness

 

All four EHI cases had at least one recorded body temperature in excess of 102.2°F (39°C) with their highest core temperatures ranging from 102.8°F (39.3°C) to 108.5°F (42.5°C) (as recorded either in the field or at the branch medical clinic). All subjects received initial cooling in the field with ice sheets and up to 14 min of additional treatment with ice sheets in the cool room of the acute care treatment area of the Branch Medical Clinic (Table 2).

Each of four EHI cases had prodromal symptoms in the days leading up to the acute event. The precipitating event involved running in three of the four cases and participation in an obstacle course in the fourth. No subject had neurological impairment more severe than mild confusion. All subjects had at least one blood sample showing elevation of unfractionated CK. With the exception of the fourth case, none of the cases had evidence of urinary heme pigment on the day of injury as judged by urine dipstick test. Case 4 had both a positive urinary dipstick test for heme pigment and a minor degree of microscopic hematuria on presentation (with 10-25 erythrocytes per high power field on urinalysis). Case 1 had evidence of mononucleosis (based on elevated serum Epstein Barr virus viral capsid IgM levels), case 3 had pharyngeal/tonsillar erythema and exudates on exam, and the chest radiograph from case 4 was read by the treating practitioners as having an infiltrate in the left lower lobe that was felt to be suspicious for pneumonia. All cases displayed elevations of serum CK above normal, although only case 4 displayed an elevation of CK above 3,000, suspicious for clinically significant rhabdomyolysis. All four cases received intravenous hydration and showed a resulting decrease in serum creatinine. There were varying degrees of leukocytosis, with a relative shift in the distribution of cell types toward predominance of granulocytes (>80%) in the blood samples taken after cooling (Table 3). In summary, each of our case subjects met the definition of EHI, although a broad range of disease severity was represented.

Number of sequences affected by EHI. Of the ~12,600 sequences present on the U95Av2 chip, 3,605 and 4,478 sequences were expressed as "present" or "marginal" in the pooled controls (drawn before and after the field-training exercise, respectively). In the cases, 4,040 sequences were similarly expressed at presentation to the acute care treatment room, 4,056 were expressed after cooling, and 3,634 were expressed at the time of follow-up. These numbers are consistent with the number of sequences consistently expressed in previous in vitro experiments with human PBMCs (~3,700) that used an earlier version of the U95 array (29).

Although the absolute number of sequences expressed during EHI did not differ significantly from those expressed at baseline, the number of sequences that showed significantly different expression between cases and controls was large. As shown in Table 5, ~1,800 sequences showed a statistically significant increase in expression at each of the time points examined in the subjects with EHI, and ~2,000 sequences showed a significant decrease in expression. However, most of these changes in expression were relatively small, with only ~800 sequences showing increases of twofold or greater at each of the time points and a similar number showing decreases (Table 5). Even fewer sequences met our presence/absence post hoc filter criteria; in total, we estimate that 361 sequences were upregulated at one or more time points and that 331 sequences were downregulated.


View this table:
[in this window]
[in a new window]
 
Table 5. Number of sequences affected by exertional heat illness

 

The sequences that were significantly affected by EHI changed as a function of time. As illustrated by the Venn diagrams in Fig. 1 (center of each diagram), only a minority of sequences were significantly upregulated (or downregulated) at all three time points. A substantial number sequences showed a time-dependent response; for example, HSPs were typically most highly upregulated on the day of injury and had diminished or returned to control levels by the time of follow-up (Table 6).



View larger version (71K):
[in this window]
[in a new window]
 
Fig. 1. Differentially expressed sequences grouped by time after exposure for signifi-cantly unregulated (left) and downregulated (right) sequences.

 

View this table:
[in this window]
[in a new window]
 
Table 6. HSPs affected by exertional heat illness

 

We classified all sequences that were significantly affected and met our post hoc filter criteria into broad functional classes. For purposes of classification and illustration, each sequence was assigned to only one primary functional class. Of the 361 upregulated sequences, about three-fourths fell into one of the following classes: immune function, 85 sequences (24%); cell growth, proliferation, differentiation, and apoptosis, 37 sequences (10%); metabolism and redox control (including heme oxygenase-1), 33 sequences (9%); transcription, 29 sequences (8%); unknown, 29 sequences (8%); protein degradation (including ubiquitins, proteases, and antiproteases), 23 sequences (6%); HSPs, chaperonins, and cochaperonins, 21 sequences (6%); and signal transduction, 22 sequences (6%). Of the 331 downregulated sequences, about two-thirds fell into one of the following classes: unknown, 58 sequences (18%); immune function, 42 sequences (13%); cell growth, proliferation, differentiation, and apoptosis, 44 sequences (13%); transcription, 41 sequences (12%); signal transduction, 32 sequences (10%); and metabolism and redox control, 18 sequences (5%).

Heat shock response. EHI produced a strong, time-dependent heat shock response in PBMCs. As shown in Table 6, a significant change in expression was found in at least one representative member of most major families of known human HSP. As might be expected from a heat-induced process, the number of significantly upregulated HSP sequences was greatest at presentation and diminished over time. At the time of follow-up, only eight HSP sequences were still significantly upregulated, and, of these, only two showed a change in expression that was twofold or more greater than controls. The finding of a widespread increase in HSP expression is highly congruent with our laboratory's previous microarray study of the effect of in vitro heat shock on gene expression in PBMCs (29), and as noted in Table 6, almost all of the HSP sequences upregulated by EHI in vivo had previously been found to be upregulated in our in vitro heat shock experiments.

Control sequences. Although EHI produced extensive changes in gene expression, these changes in expression were not universal. Table 7 shows the effect of EHI on a number of control sequences, including {beta}-actin, GAPDH, cyclophilin A, and the five ribosomal protein sequences most highly expressed in the control subjects. Although some of these sequences showed changes in expression that were significantly different from the controls, none showed an absolute change of twofold or greater, the cutoff used in our post hoc filtering process.


View this table:
[in this window]
[in a new window]
 
Table 7. Control sequences

 

Effect of EHI on non-HSPs. The non-HSPs most strongly (>=5.0-fold after rounding, at one or more time points) affected by EHI are listed in Tables 8 and 9. Importantly, of the 36 upregulated sequences listed in Table 8, at least one-fourth are known to be interferon inducible. Interestingly, in comparing these results to previous in vitro data (29), only two of these non-HSPs had previously been found to be significantly up-regulated after an in vitro heat shock delivered to normal PBMCs (43°C x 20 min followed by 2 h 40 min of recovery at 37°C): the phosphatidylserine receptor (upregulated 3.1-fold in vitro) and pyruvate carboxylase (upregulated 1.8-fold in vitro but did not meet our post hoc filter criteria in the previous study). Furthermore, seven of these upregulated genes were actually downregulated in our previous in vitro experiment: granzyme B (down 0.49-fold in vitro), adrenomedullin (down 0.47-fold in vitro), complement component 3a receptor 1 (down 0.40-fold in vitro), myxovirus resistance A (down 0.23-fold in vitro), endothelial cell growth factor 1 (down 0.38-fold in vitro), interferon inducible protein 60 (down 0.12-fold in vitro), and interferon-induced protein with tetratricopeptide 2 (down 0.22-fold in vitro but did not meet the previous study's post hoc filter criteria). Among the differences between the two experiments, it is noteworthy that the in vitro experiment was performed in serum-free media.


View this table:
[in this window]
[in a new window]
 
Table 8. Non-HSP genes most strongly upregulated by EHI

 

View this table:
[in this window]
[in a new window]
 
Table 9. Non-HSP genes most strongly downregulated by EHI

 

A sequence-encoding interferon-{gamma} (GenBank no. X13274 [GenBank] ) was significantly upregulated by EHI, with changes of 3.7-fold at presentation, 3.0-fold after cooling, and 7.6-fold at followup. However, this sequence showed a low level of expression and did not meet our post hoc presence/absence filter criteria. In our laboratory's previous in vitro experiment (29), this sequence was upregulated 4.3-fold, but the variance from experiment to experiment was high enough to render the change not statistically significant (increased in 4 experiments, decreased in 1 experiment; 95% confidence interval, 0.66 to 28). Thus, both in vitro and in vivo, our microarray evidence is equivocal but favors the possibility that heat stress leads to increased expression of interferon-{gamma} mRNA.

Two interferon-{alpha} sequences on the Affymetrix array (GenBank V00541 [GenBank] and X02956 [GenBank] , both encoding interferon-{alpha}5) showed statistically significant increases of twofold or greater at presentation (2.4-fold and 2.1-fold, respectively) but not after cooling or at recovery. However, neither of these sequences met our post hoc criteria for inclusion, and neither had been upregulated in our laboratory's previous in vitro study (29).

Among the 35 most strongly downregulated sequences, only 5 were previously found to be downregulated as a result of in vitro PBMC heat shock (Table 9). It is noteworthy that the time point at which maximum downregulation occurred was most commonly at follow-up (1-2 days after initial presentation), which is substantially later than the time point examined in the in vitro experiment (2 h 40 min after heat shock).

Confirmatory RT-PCR. RT-PCR was performed on a select number of genes identified by the microarray as being strongly upregulated by heat shock (Fig. 2). The source material for this PCR came from control subject B (lanes 1 and 2) and case 2 (lanes 3-5). As illustrated, the PCR confirmed the presence of an increased expression of HSP70B', phosphatidylserine receptor, interferon-induced proteins 27 and 60, and AB000115 [GenBank] (an open reading frame on chromosome 1) in the EHI case. By contrast, there was no clear effect on cyclophilin A, as predicted by the microarray results.



View larger version (68K):
[in this window]
[in a new window]
 
Fig. 2. Confirmatory RT-PCR. The controls are from subject B, taken before (lane 1) and after (lane 2) participation in a field-training exercise. The heat-injury samples are from subject 2, taken at presentation (lane 3), after cooling (lane 4), and at follow-up (lane 5). Sequences studied were heat shock protein (HSP) 70B'; interferon-induced protein 60 (IFI-60); interferon-{alpha}-inducible protein 27 (IFI-27); phosphatidylserine receptor (PSR); chromosome 1 open reading frame 29 (AB000115 [GenBank] ); and cyclophilin A. The ladder bands occur in 100-bp increments. The arrows point to the bright 500-bp band in each ladder.

 

Comparison to previous in vitro heat shock in PBMCs. Previous work in our laboratory (27) has suggested that human cells display a relatively small nonspecific response to environmental stress in addition to stress- and cell type-specific responses. From the data in Tables 6, 8, and 9, it became apparent that the PBMCs of the subjects with EHI had evidence of not only a conventional heat shock response (defined as increased expression of HSPs and other genes previously found to be upregulated by heat shock in vitro) but also of a response that differed from the previously observed in vitro effects of heat shock on PBMCs. To better estimate the extent to which our in vivo responses could be accounted for by thermal stress alone, we compared the list of 361 upregulated and 331 downregulated genes to the list of all genes signifi-cantly affected by heat shock in our in vitro PBMC experiment (29). Similarly, to estimate the magnitude of the stress nonspecific response (previously estimated at ~10-15% of all sequences affected), we compared the changes in expression noted here to those found in a database of cells exposed to a nonthermal stress (HepG2 cells exposed to hypoxia in vitro for 24 h) (27). To reduce the biasing effect of post hoc filtering, we included all sequences from the in vitro experiments that showed a statistically significant change in expression, independent of the magnitude of change or the presence/absence calls.

Of the 144 sequences upregulated at presentation, 35 (24%) had previously been found to be significantly upregulated by in vitro heat shock. By contrast, only 13 (9%) were similarly affected by hypoxic exposure in HepG2 cells (P < 0.001 by {chi}2 analysis). Of the 237 sequences upregulated after cooling, 41 (17%) were similarly affected by in vitro heat shock, which was comparable to the number of sequences affected by hypoxia [34 (14%), P = 0.45]. When we compared the 179 sequences that were upregulated at recovery, we found that only 17 (9%) were similarly affected by in vitro heat shock and 15 (8%) by hypoxia (P = 0.85). Although these numerical findings would suggest an overlap between the three studies of ~10-15%, when the actual sequences affected were compared, only eight sequences representing seven unique genes were significantly upregulated in all three studies at one or more time points. These included the well-established stress genes ubiquitin C, dual-specificity phosphatase-1, and colligin 2 (HSP47), as well as the phosphatidylserine receptor, calcium-binding protein S100A2, pyruvate carboxylase, and the collagen synthesis enzyme lysine hydroxylase.

Among the downregulated sequences, 11-18% have previously been found to be significantly downregulated by in vitro heat shock in PBMCs or by hypoxic exposure in HepG2 cells. Neither of the two in vitro experiments showed a significantly greater degree of overlap than the other with the effects of EHI.


    DISCUSSION
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Our most important finding is that the PBMC gene expression response to EHI is far more extensive than previously recognized. This response appears to include a time-dependent series of changes that has at least three components: 1) a heat shock response that involves elements such as HSPs, as described in PBMCs and other cell lines in vitro; 2) a response that includes a number of genes not previously found to be upregulated by heat shock in our in vitro experiments in PBMCs and that includes a substantial number of interferoninducible genes; and 3) a small nonspecific stress response that is shared by other cell lines and stressors.

Although limited by several technical considerations (discussed below), the findings reported here are biologically plausible within the context of what is known about thermal stress in vitro and in vivo. For example, a time-dependent heat shock response was observed, a phenomenon that has been noted in vivo in contexts such as intense exercise (7, 8). It is also noteworthy that a sequence corresponding to c-myc was among the most highly downregulated sequences at presentation, because in vitro experiments have found that decreased expression of this gene is essential to the process of cellular recovery from severe thermal stress (31). Additionally, the congruence between the findings of the present study and our previous in vitro study in PBMCs was greatest for the samples obtained on the day of presentation (17-24%), as would be expected given the time point used in the in vitro study (2 h 40 min after heat shock). Importantly, the observed differences between the in vivo and in vitro responses suggest that, as might be anticipated from a multifactorial process, EHI also involves pathways other than those affected by in vitro heat shock in serum-free media. The molecular data reported here likely reflect not only the direct cellular effects of thermal stress on PBMCs but also the responses of PBMCs to the humoral or other signals detectable by cells traveling through tissues reacting to thermal stress.

It is noteworthy that many of the sequences most highly increased by EHI in this report are interferon-inducible genes (Table 8). Interferons-{alpha} and -{gamma} have been variably reported to be elevated in exertional heat illness (reviewed in Ref. 2), and these immune modulators are known to be capable of producing flu-like syndromes, at least when administered at pharmacological doses (1, 15). As is commonly observed in EHI (5, 26), all of our cases reported feeling ill in the days leading up to the acute events, and at least one of the case subjects had firm evidence of an interferon-associated viral process (mononucleosis) on clinical evaluation. The finding of increased expression of interferon-induced genes at the time of presentation may simply be a molecular correlate of antecedent viral infections in our cases but may also identify candidate mediators that might account for the apparent association of EHI with prodromal symptoms suggestive of viral illness. Alternatively, at least some of the increases in expression of these sequences might reflect heat-induced increases in plasma levels of interferons. Further research will be required to establish whether the pathways associated with these sequences participate in the pathophysiology of EHI or whether their increased expression merely serves as markers that identify individuals at increased risk of EHI. Importantly, our microarray evidence to date suggests that the interferon-induced sequences (Table 8) are not directly upregulated by heat itself, because their expression was not significantly affected in our laboratory's previous in vitro heat shock experiment, and, indeed, at least three were found to be downregulated in vitro (29).

A comparison of the present work to our laboratory's previous in vitro work suggests that, as is generally believed, there is a cell type and stressor nonspecific component to the response to stress at the level of gene expression. This nonspecific component includes well-known stress genes such as ubiquitin C, HSP47/colligin 2, and dual-specificity phosphatase-1. The limited microarray evidence our laboratory has produced to date (including the present work) suggests that the number of genes involved in this nonspecific component may be a smaller fraction of the genomic response than previously thought, amounting to no more than ~10-15% of the total number of genes affected by any given environmental stress (27).

Our work confirms and extends, in important respects, the findings made in a mouse model of thermal stress reported recently by Zhang et al. (32). Among other findings, these authors reported that young mice exposed to two sequential thermal stresses had significant changes in liver gene expression in four principal categories: stress response genes (including HSPs); cell growth, death, and signaling-related genes; antioxidant enzymes and drug metabolism enzymes; and DNA/RNA/protein repair-related genes. Although we found genes related to immune function to be the largest single category affected by EHI (as might be expected from the fact that we studied PBMCs rather than other cell types), we also found a large number of genes to be affected by EHI that are functionally thought to be involved in the cell stress response, signal transduction, metabolism/redox control, and cell growth, proliferation, differentiation, and apoptosis (Tables 8 and 9). This suggests that a number of the functional responses to thermal stress at the level of RNA expression are conserved across species and cell types and adds support to the use of rodent models for the study of heat-related illnesses.

There are several important limitations to this study. First, as an observational rather than an experimental study, the general applicability of its conclusions are limited by the low number of subjects and the presence of significant intersubject variation in severity of illness, prodromal symptoms, clinical course, and timing of sample collection. Accordingly, we are unable to determine the extent to which particular individual sequences might influence severity of disease from the present data. Nonetheless, although the degree of severity was broad, all subjects met clinical criteria for EHI (11). Second, although our controls were matched for gender, ethnic origin, and both season and year of exposure, they were not drawn in parallel with the cases during the events that precipitated EHI. It is, therefore, possible that some of the gene-expression changes identified in our subjects suffering from EHI reflect the effects of variables, such as different types or levels of physical activity, and are not pathophysiologically related to EHI itself. Indeed, increases in leukocyte HSP expression have previously been described in individuals undergoing intense physical activity (6-8, 24), although it can be difficult in these studies to separate the effects of exercise from the effects of exercise-induced increases in body temperature. Third, the low number of replicates (only 2 control chips) coupled with the high number of sequences tested (~12,600) makes it likely that some of the changes in expression we report as significant represent false-positive results. We have attempted to reduce this by application of strict post hoc filter criteria. Fourth, it is possible that some of the changes in gene expression noted reflect a shift in the relative distributions of different subpopulations of circulating PBMCs rather than changes in expression within the same cells. Although this is unlikely to be the case for molecules such as HSPs (which are almost universally expressed by cells in response to thermal stress), it could account for some of the apparent changes in expression of genes known to be highly expressed by specific cell types [e.g., granzyme B, which is selectively expressed by cytolytic T cells and natural killer cells (30)]. Finally, we do not have measurements of plasma interferon levels.

In summary, this is the first large-scale survey of gene expression changes related to EHI or EHS. Our findings suggest that EHI produces gene expression changes that are far more extensive than previously realized. The gene expression signatures that occur during EHI include a classic heat shock response, a small nonspecific cell stress response, and a response that cannot be fully explained by our previous in vitro work. This last component involves a significant number of interferon-inducible sequences and might, therefore, account for the known association of EHI with prodromal flulike symptoms. Finally, PBMCs provide an easily obtainable tissue to study gene expression with EHI and EHS.


    GRANTS
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
This study was sponsored in part by the Navy Bureau of Medicine and Surgery Clinical Investigation Program No. P00-L-H00000-025:A and by National Heart, Lung, and Blood Institute Grant HL-64104.


    ACKNOWLEDGMENTS
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
The authors thank Michael Cullivan for assistance in sample preparation and Dr. Richard Pratt, whose laboratory performed the gene chip hybridizations and first-pass MAS software analysis. The authors also thank Dr. Jacques Reifman, Dr. Sorin Draghichi, and the US Army Medical Research and Material Command Bioinformatics Cell for valuable discussions pertaining to the statistical approach used in this manuscript.

The views, opinions, and findings contained in this publication are those of the authors and should not be construed as an official US Department of the Army or Department of the Navy position, policy, or decision, unless so designated by other documentation.

For protection of test subjects, the investigators adhered to the protections of 45 CFR 46.

Approved for public release; distribution unlimited.


    FOOTNOTES
 

Address for reprint requests and other correspondence: L. A. Sonna, Thermal and Mountain Medicine Div., USARIEM, Natick, MA 01760 (E-mail: larry.sonna{at}na.amedd.army.mil).

The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

{dagger} Deceased 22 November 2002. Back


    REFERENCES
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 

  1. Bonnem EM and Oldham RK. Gamma-interferon: physiology and speculation on its role in medicine. J Biol Response Mod 6: 275-301, 1987.
  2. Bouchama A and Knochel JP. Heat stroke. N Engl J Med 346: 1978-1988, 2002.
  3. Costrini AM, Pitt HA, Gustafson AB, and Uddin DE. Cardiovascular and metabolic manifestations of heat stroke and severe heat exhaustion. Am J Med 66: 296-302, 1979.
  4. DuBose DA, McCreary J, Sowders L, and Goode L. Relationship between rat heat stress mortality and alterations in reticuloendothelial carbon clearance function. Aviat Space Environ Med 54: 1090-1095, 1983.
  5. Epstein Y, Moran DS, Shapiro Y, Sohar E, and Shemer J. Exertional heat stroke: a case series. Med Sci Sports Exerc 31: 224-228, 1999.
  6. Fehrenbach E, Niess AM, Schlotz E, Passek F, Dickhuth HH, and Northoff H. Transcriptional and translational regulation of heat shock proteins in leukocytes of endurance runners. J Appl Physiol 89: 704-710, 2000.
  7. Fehrenbach E, Niess AM, Veith R, Dickhuth HH, and Northoff H. Changes of HSP72 expression in leukocytes are associated with adaptation to exercise under conditions of high environmental temperature. J Leukoc Biol 69: 747-754, 2001.
  8. Fehrenbach E, Passek F, Niess AM, Pohla H, Weinstock C, Dickhuth HH, and Northoff H. HSP expression in human leukocytes is modulated by endurance exercise. Med Sci Sports Exerc 32: 592-600, 2000.
  9. Flanagan SW, Ryan AJ, Gisolfi CV, and Moseley PL. Tissue-specific HSP70 response in animals undergoing heat stress. Am J Physiol Regul Integr Comp Physiol 268: R28-R32, 1995.
  10. Fruth JM and Gisolfi CV. Work-heat tolerance in endurance-trained rats. J Appl Physiol 54: 249-253, 1983.
  11. Gardner JW and Kark JA. Clinical diagnosis, management, and surveillance of exertional heat illness. In: Medical Aspects of Harsh Environments, edited by Pandolf KB and Burr RE. Washington, DC: Office of the Surgeon General, Department of the Army, USA/Borden Institute, 2001, vol. 1, p. 231-279.
  12. Hall DM, Buettner GR, Oberley LW, Xu L, Matthes RD, and Gisolfi CV. Mechanisms of circulatory and intestinal barrier dysfunction during whole body hyperthermia. Am J Physiol Heart Circ Physiol 280: H509-H521, 2001.
  13. Hall DM, Xu L, Drake VJ, Oberley LW, Oberley TD, Moseley PL, and Kregel KC. Aging reduces adaptive capacity and stress protein expression in the liver after heat stress. J Appl Physiol 89: 749-759, 2000.
  14. Hubbard RW, Bowers WD, Matthew WT, Curtis FC, Criss RE, Sheldon GM, and Ratteree JW. Rat model of acute heat stroke mortality. J Appl Physiol 42: 809-816, 1977.
  15. Koeller JM. Biologic response modifiers: the interferon alfa experience. Am J Hosp Pharm 46: S11-S15, 1989.
  16. Kregel KC. Heat shock proteins: modifying factors in physiological stress responses and acquired thermotolerance. J Appl Physiol 92: 2177-2186, 2002.
  17. Lambert GP, Gisolfi CV, Berg DJ, Moseley PL, Oberley LW, and Kregel KC. Selected contribution: Hyperthermia-induced intestinal permeability and the role of oxidative and nitrosative stress. J Appl Physiol 92: 1750-1761, 2002.
  18. Lindquist S. The heat-shock response. Annu Rev Biochem 55: 1151-1191, 1986.
  19. Montain SJ, Latzka WA, and Sawka MN. Impact of muscle injury and accompanying inflammatory response on thermoregulation during exercise in the heat. J Appl Physiol 89: 1123-1130, 2000.
  20. Moseley PL. Exercise, heat, and thermotolerance: molecular mechanisms. In: Exercise, Heat, and Thermoregulation, edited by Gisolfi CV, Lamb DR, and Nadel ER. Dubuque, IA: William C. Brown, 1993, p. 305-334.
  21. Parsell DA and Lindquist S. The function of heat-shock proteins in stress tolerance: degradation and reactivation of damaged proteins. Annu Rev Genet 27: 437-496, 1993.
  22. Ryan AJ. Heat stroke and endotoxin: sensititzation or tolerance to endotoxins? In: Exercise, Heat, and Thermoregulation, edited by Gisolfi CV, Lamb DR, and Nadel ER. Dubuque, IA: Brown, 1993, p. 335-385.
  23. Ryan AJ, Gisolfi CV, and Moseley PL. Synthesis of 70K stress protein by human leukocytes: effect of exercise in the heat. J Appl Physiol 70: 466-471, 1991.
  24. Schneider EM, Niess AM, Lorenz I, Northoff H, and Fehrenbach E. Inducible HSP70 expression analysis after heat and physical exercise: transcriptional, protein expression, and subcellular localization. Ann NY Acad Sci 973: 8-12, 2002.
  25. Shapiro Y and Seidman DS. Field and clinical observations of exertional heat stroke patients. Med Sci Sports Exerc 22: 6-14, 1990.
  26. Shibolet S, Coll R, Gilat T, and Sohar E. Heatstroke: its clinical picture and mechanism in 36 cases. Q J Med 36: 525-548, 1967.
  27. Sonna LA, Cullivan ML, Sheldon HK, Pratt RE, and Lilly CM. Effect of hypoxia on gene expression by human hepatocytes (HepG2). Physiol Genomics 12: 195-207, 2003.
  28. Sonna LA, Fujita J, Gaffin SL, and Lilly CM. Invited review: Effects of heat and cold stress on mammalian gene expression. J Appl Physiol 92: 1725-1742, 2002.
  29. Sonna LA, Gaffin SL, Pratt RE, Cullivan ML, Angel KC, and Lilly CM. Effect of acute heat shock on gene expression by human peripheral blood mononuclear cells. J Appl Physiol 92: 2208-2220, 2002.
  30. Trapani JA. Granzymes: a family of lymphocyte granule serine proteases. Genome Biol 2: 3014, 2001.
  31. Wennborg A, Classon M, Klein G, and von Gabain A. Downregulation of c-myc expression after heat shock in human B-cell lines is independent of 5' mRNA sequences. Biol Chem Hoppe Seyler 376: 671-680, 1995.
  32. Zhang HJ, Drake VJ, Morrison JP, Oberley LW, and Kregel KC. Selected Contribution: Differential expression of stress-related genes with aging and hyperthermia. J Appl Physiol 92: 1762-1769, 2002.



This article has been cited by other articles:


Home page
Am. J. Physiol. Regul. Integr. Comp. Physiol.Home page
J. P. McClung, J. D. Hasday, J.-r. He, S. J. Montain, S. N. Cheuvront, M. N. Sawka, and I. S. Singh
Exercise-heat acclimation in humans alters baseline levels and ex vivo heat inducibility of HSP72 and HSP90 in peripheral blood mononuclear cells
Am J Physiol Regulatory Integrative Comp Physiol, January 1, 2008; 294(1): R185 - R191.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
J. T. Selsby, S. Rother, S. Tsuda, O. Pracash, J. Quindry, and S. L. Dodd
Intermittent hyperthermia enhances skeletal muscle regrowth and attenuates oxidative damage following reloading
J Appl Physiol, April 1, 2007; 102(4): 1702 - 1707.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
E. Fehrenbach
Multifarious microarray-based gene expression patterns in response to exercise
J Appl Physiol, January 1, 2007; 102(1): 7 - 8.
[Full Text] [PDF]


Home page
Physiol. GenomicsHome page
L. A. Sonna, M. M. Kuhlmeier, H. C. Carter, J. D. Hasday, C. M. Lilly, and K. D. Fairchild
Effect of moderate hypothermia on gene expression by THP-1 cells: a DNA microarray study
Physiol Genomics, September 14, 2006; 26(1): 91 - 98.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
D. S. Moran, L. Eli-Berchoer, Y. Heled, L. Mendel, M. Schocina, and M. Horowitz
Heat intolerance: does gene transcription contribute?
J Appl Physiol, April 1, 2006; 100(4): 1370 - 1376.
[Abstract] [Full Text] [PDF]


Home page
Physiol. GenomicsHome page
D. Zieker, E. Fehrenbach, J. Dietzsch, J. Fliegner, M. Waidmann, K. Nieselt, P. Gebicke-Haerter, R. Spanagel, P. Simon, A. M. Niess, et al.
cDNA microarray analysis reveals novel candidate genes expressed in human peripheral blood following exhaustive exercise
Physiol Genomics, November 17, 2005; 23(3): 287 - 294.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Regul. Integr. Comp. Physiol.Home page
J. T. Selsby and S. L. Dodd
Heat treatment reduces oxidative stress and protects muscle mass during immobilization
Am J Physiol Regulatory Integrative Comp Physiol, July 1, 2005; 289(1): R134 - R139.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Regul. Integr. Comp. Physiol.Home page
M. B. Reid
Response of the ubiquitin-proteasome pathway to changes in muscle activity
Am J Physiol Regulatory Integrative Comp Physiol, June 1, 2005; 288(6): R1423 - R1431.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
C. M. Lilly, H. Tateno, T. Oguma, E. Israel, and L. A. Sonna
Effects of Allergen Challenge on Airway Epithelial Cell Gene Expression
Am. J. Respir. Crit. Care Med., March 15, 2005; 171(6): 579 - 586.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow All Versions of this Article:
96/5/1943    most recent
00886.2003v1
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (23)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Sonna, L. A.
Right arrow Articles by Lilly, C. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Sonna, L. A.
Right arrow Articles by Lilly, C. M.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2004 by the American Physiological Society.