Alteration in circulating metabolites during and after heat stress in the conscious rat: potential biomarkers of exposure and organ-specific injury
© Ippolito et al.; licensee BioMed Central. 2014
Received: 25 June 2014
Accepted: 11 December 2014
Published: 24 December 2014
Heat illness is a debilitating and potentially life-threatening condition. Limited data are available to identify individuals with heat illness at greatest risk for organ damage. We recently described the transcriptomic and proteomic responses to heat injury and recovery in multiple organs in an in vivo model of conscious rats heated to a maximum core temperature of 41.8°C (Tc,Max). In this study, we examined changes in plasma metabolic networks at Tc,Max, 24, or 48 hours after the heat stress stimulus.
Circulating metabolites were identified by gas chromatography/mass spectrometry and liquid chromatography/tandem mass spectrometry. Bioinformatics analysis of the metabolomic data corroborated proteomics and transcriptomics data in the tissue at the pathway level, supporting modulations in metabolic networks including cell death or catabolism (pyrimidine and purine degradation, acetylation, sulfation, redox alterations and glutathione metabolism, and the urea cycle/creatinine metabolism), energetics (stasis in glycolysis and tricarboxylic acid cycle, β-oxidation), cholesterol and nitric oxide metabolism, and bile acids. Hierarchical clustering identified 15 biochemicals that differentiated animals with histopathological evidence of cardiac injury at 48 hours from uninjured animals. The metabolic networks perturbed in the plasma corroborated the tissue proteomics and transcriptomics pathway data, supporting a model of irreversible cell death and decrements in energetics as key indicators of cardiac damage in response to heat stress.
Integrating plasma metabolomics with tissue proteomics and transcriptomics supports a diagnostic approach to assessing individual susceptibility to organ injury and predicting recovery after heat stress.
KeywordsHeat stress Metabolomics Systems biology Energetics Metabolic networks
Heat illness and heat stroke, the most severe form of heat illness, are life-threatening conditions characterized by elevations in core temperature (Tc) resulting from an inability to adequately dissipate excess body heat to the environment. The intrinsic nature of military operations (i.e., heavy physical activity in extreme environments) places military personnel at greater risk of developing heat-related illnesses. Sustained military operations in the Middle East have been accompanied by an increase in the number of heat stroke cases on the battlefield and in training. During the past two decades alone the incidence of heat stroke has increased over seven-fold . Moreover, the 30-year mortality rates from heart, kidney, and liver failure in US forces increases by 40% in individuals with a history of heat stroke ,. In 2013, US forces sustained over 2000 heat-related injuries requiring hospitalization, including 324 cases of heat stroke . The actual incidence is projected to be considerably higher when considering undocumented instances that never reach triage ,. In a five-year retrospective study, 10,319 cases of heat injury required medical resources to treat, including 1872 cases of heat stroke . The estimated cost to the military approaches $52 million per year (USD), assuming a cost of $6200 per day, an average hospital stay of 3.2 days, and 2620 cases of heat injuries and heat stroke per year . Other costs include duty days and salary lost during recovery (up to 5.5 months) and the loss in investment and training associated with service discharge for medical reasons . Taken together, these statistics indicate that there is an urgent need for earlier indicators of organ injury and susceptibility and molecular based indicators to improve return to duty decisions .
The physiological responses to an excess heat load include elevations in heart rate, a drop in mean arterial pressure, attenuated sweating rates, stupor, and collapse. The molecular level alterations and influence of the physiological events preceding and contributing to the systemic inflammatory response associated with heat stroke remain largely unknown. Heat stroke compromises tight junction integrity in the gut, resulting in leakage of bacteria into the circulation . If uncontrolled, the ensuing thermoregulatory and immune responses can progress to system inflammatory response syndrome (SIRS) and ultimately multi-organ failure and death -. Designing novel and effective treatment and detection strategies for heat stroke requires a better understanding of the physiological and molecular alterations that accompany heat stroke and characterize the SIRS event. However, in human studies 39.0 – 39.5°C is the highest ethically attainable maximal core temperature (Tc,Max), but these temperatures are insufficient to induce heat stroke. To overcome these ethical limitations, we recently described conscious rat and mouse models of heat illness ,,-. The models use abdominally implanted radiotelemetry units to track Tc while supporting physiological and behavioral adjustments to heat stress. In a recent study, we used our physiological model of heat stress to conduct an integrated systems biology evaluation of transcriptomic and proteomic changes in heart, liver, kidney, and lung after heat stress, heat injury, and recovery . We identified discriminatory gene and protein signatures in heat-injured cardiac tissue reflecting perturbations in oxidative phosphorylation, energy production, and inflammatory response.
These global changes in metabolic networks associated with energy production in target tissues suggest that more accessible biofluids (e.g., serum, plasma, urine, and saliva) are also likely to reflect changes in the physiological state of the heat-stressed and heat-injured organism. Metabolomic profiling of plasma in conjunction with proteomic and transcriptomic analysis has recently emerged as a powerful predictive tool reflecting the dynamic responses to genetic modification and physiological, pathophysioloigcal, and/or developmental stimuli. However, perturbation in metabolic networks has not been well studied in heat illness. Understanding the metabolic response to heat illness and subsequent organ injury provides unique insight into understanding how mammalian systems react to heat illness and recovery. Further, methods to interrogate metabolites in accessible biofluids have been developed to allow a global assessment of organism response to environmental stressors . One disadvantage of metabolomics is the likelihood of false positives given the metabolome’s exquisite sensitivity to subtle changes in physiology (e.g., food intake, changes in temperature, stress, etc.). Therefore, any study of the metabolome must discriminate pathological changes in metabolic networks from changes inherent in normal physiological functioning. Integrating metabolic networks at the level of molecular function is one approach to differentiate changes related to heat-stress and/or heat-injury from physiological variation in the unperturbed system. Bioinformatics methods can be used to identify an integrated panel of multiple physiological networks and biomarkers.
In this study, we evaluate the metabolomic profile in response to heat stress and heat injury, and compare the results with the proteomic and transcriptomic profiles described previously. We demonstrate plasma metabolomics profiles unique to the physiological conditions of heat stress, heat injury in cardiac tissue, and heat stress without injury. Metabolomics network analysis demonstrates perturbations in biological processes associated with energy usage and cell death, similar to the transcriptomics and proteomics analyses conducted previously. Further, integrated analysis of variance (ANOVA) analysis, random forest analysis, and hierarchical clustering analysis identify panels of biochemicals differentiating controls from heat-stressed and heat-injured animals.
Results and discussion
Thermoregulation and histopathology
Tc,Max was reached at 2–3 hours, as previously published by our laboratories in more detailed thermoregulation analyses conducted with these animals and other similar studies ,. Animals in the 24 and 48 hour cohorts were allowed to recover for the specified time prior to termination. In thermoregulation studies with these animals, temperatures return to baseline by 24–48 hours, but cardiac-injured animals exhibit hypothermia before recovery, as published in detail by our laboratory .
Liver, kidney, and heart were assessed by histopathology and scored for evidence of injury as described previously . The heart was the only tissue with treatment-related evidence for injury, with inflammation and cardiomyopathy in three out of six animals in the heated cohort as described in detail in a recent publication from our laboratory  (Additional file 1A). For the remainder of the study, animals in the 48 hour treatment group were sub-stratified into cohorts according to uninjured vs. injured.
The metabolomics heat stress response and recovery in plasma—biochemical networks
Overview of metabolic network changes
Tc,Max showed the greatest percentage of biochemicals altered in response to heat injury in virtually all categories. At Tc,Max, nearly 100% of identified biochemicals associated with cell death in a given pathway were elevated, suggesting an overall increase in biosynthesis to combat an overt exposure to redox stress (Figure 1). Biochemicals in all enriched metabolic pathways relating to cell death were up-regulated at Tc,Max, with progressively fewer biochemicals within each pathway elevated and/or reverting to down-regulation over time (Figure 1).
List of seven key biochemical differentially expressed in heat-stressed and control animals*
Identified biochemical (pathway)
Literature review (biochemical and/or pathway)
5,6-Dihydrouracil (pyrimidine metabolism)
Early signal of apoptosis; DNA damage from reactive oxygen species elevation; deficiency
3-Ureidopropionate (alanine/aspartate metabolism)
Increased reactive oxygen species; inhibition of mitochondrial energy metabolism; neurotoxic/excitotoxic
Ornithine (urea cycle, arginine metabolism)
Production of arginine and increase in autophagy, cell death, removal of excess NH4+, energy metabolism; slight renal dysfunction
Glutathione disulfide, oxidized (redox)
Apoptosis, DNA damage, cell proliferation, survival, differentiation, metabolism; redox stress and/or crisis due to elevated reactive oxygen species
Corticosterone (steroid/sterol metabolism)
Mitochondrial conversion of acetyl-CoA to cholesterol and conversion of cholesterol to corticosterone; adrenal cortex production of cholesterol to pregnenalone, and ultimately cortisol; involved in glucocorticoid activity and stress response
7-α-Hydroxy-3-oxo-4-cholestenoate [7-HOCA] (steroid/sterol metabolism)
Bile acid synthesis from cholesterol; CYP7A1 activity in the liver; bile acid synthesis
Pyridoxal (vitamin B metabolism)
decreased renal function; Anabolism; cofactor for reaction releasing glucose from glycogen
Redox stress and cell death
Production of reactive oxygen species (ROS) and oxidative stress occur when the cell’s protective mechanisms are saturated. ROS also play key signaling roles in cell proliferation, survival, disease and pathophysiology . Increased ROS may lead to apoptosis, DNA damage, cell proliferation, survival, differentiation, and disruption in metabolism . Heat exposure elevates free radicals . Free radicals and other ROS are neutralized by GSH in redox regulation cycles. Deficient and/or dysfunctional GSH can disrupt cell processes and lead to cell death . Oxidized GSH can be reduced back by glutathione reductase with NADPH (nicotinamide adenine dinucleotide phosphate) as an electron donor. GSH donates a reducing equivalent to unstable redox species. Under normal physiological conditions, more than 90% of the glutathione pool is in the GSH form and less than 10% is in the disulfide GSSG form. The ratio of GSH to GSSG is a metric of cellular toxicity, with elevated GSSG indicating oxidative stress . Taken together, these results suggest that early perturbations in metabolic networks elicit persistent redox stress out to 48 hours after exposure. These persistent metabolic changes in heat-stressed individuals may affect performance during subsequent episodes of heat injury.
Arginine metabolism, nitric oxide metabolism, and cell death
Elevated citrulline levels suggest a concomitant increase in nitric oxide, a product of inflammatory signaling (Figure 3B) . Increased urea may be associated with cell death, autophagy, and/or protein degradation. The combination of elevated urea and creatinine could also suggest renal dysfunction, consistent with the proteinosis observed in the renal histopathology as reported in the companion study. Elevated plasma creatinine at later time points could indicate continued use of skeletal muscle energy reserves and persistent renal dysfunction. Concomitantly decreased arginine and increased urea may indicate slower kinetics of arginine metabolism with increased urea production (Figure 3B). Increased urea may also be the result of an increase in NH4+, a byproduct of cellular death .
At 24 and 48 hours, citrulline levels were no longer elevated, suggesting that nitric oxide (NO) production had decreased. Interestingly, however, ornithine levels increased at 24 and 48 hours, possibly reflecting slower arginine metabolism kinetics in general . Taken together, these data suggest that a slight renal dysfunction was maintained during recovery, consistent with histopathologic evidence suggesting renal dysfunction.
Purine/Pyrimidine metabolism and cell death
Fourteen out of 18 identified biochemicals in the pyrimidine and purine metabolic pathways were significantly different at Tc,Max relative to control (Figure 1 and Additional file 2). Exceptions were adenosine monophosphate (AMP), uridine, cytosine, and dihydroorotate. Uracil, 5,6-dihydrouracil, 3-ureidopropionate, β-alanine, and an acetylated form of β-alanine, N-acetyl-β-alanine, were all significantly increased at Tc,Max compared to control rats, suggesting that a significant proportion of biochemicals involved in pyrimidine degradation was altered (Additional file 2). Significantly, nearly all the biochemicals immediately upstream of β-alanine metabolism in the pyrimidine metabolism reference Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway are modulated, suggesting that a significant proportion of this biochemical network is affected by heat stress at Tc,Max. In plant cells, alterations in pyrimidine nucleotide metabolism are considered an early signal of apoptosis and often are induced by an increase in endogenous NO . 3-ureidopropionate, the substrate of the enzyme dihydropyrimidinase, was elevated (Additional file 1) while the product, 5,6-dihydrouracil, was decreased (Additional file 2), suggesting a down-regulation in the level of dihydropyriminidase activity or 3-ureidopropionate. Further, carnosine is synthesized from L-histidine and β-alanine, with β-alanine as the rate-limiting precursor , and 3-ureidopropionate as the precursor to β-alanine production. The result of these alterations in pyrimidine degradation could suggest DNA damage by endogenous ROS over-riding cellular repair and protection mechanisms. Several altered biochemical pathways support elevated ROS.
Acetylation or sulfation and cell death
Of the N-acetylated biochemicals detected, 14 of 15 compounds were significantly elevated at Tc,Max (Additional file 3). While 8 of 15 biochemicals were significantly changed at 24 hours, 6 remained altered at 48 hours. Heat exposure appeared to change in the activity of acetyltransferases or deacetylases, consistent with cellular damage. This initial response probably reflects hypoxia and initiation of cell death . Acetyltransferases and deacetylases may allow an organism to adapt to recovery.
Eight sulfated compounds were detected and at Tc,Max, of which seven were significantly increased (Additional file 4). In addition, seven unnamed, putatively sulfated compounds also increased at Tc,Max (X-12182, X-12183, X-12184, X-12185, X-12230 and X-12307; Additional file 1). At 24 hours, all biochemicals which had been higher than control at Tc,Max were now significantly lower than controls, with a reversal of the trend after 48 hours (i.e., higher than controls; Additional file 4). Sulfotransferase enzymatic activity may have been induced or the activities of sulfatases may have reduced initially to facilitate clearance of cellular debris , then returned to control levels after recovery.
Amino acids and cell death
Of the 19 identified amino acids, seven were decreased and eight were increased relative to controls at Tc,Max, with a trend toward overall down-regulation at 24 hours and 48 hours. The exception was lysine, with a trend toward an increase at Tc,Max, with elevations at both 24 and 48 hours (Additional file 5). At Tc,Max, serine, threonine and alanine were lower than control. All three amino acids are synthesized from glycolytic intermediates. Although not definitive, these results may suggest a decrease in glycolysis . Moreover, glutamate, which is consumed in the production of these amino acids, was increased at Tc,Max. Whether these amino acid changes are indicative of energetic alterations or increased autophagy is unclear . However, it is interesting to note that, at Tc,Max, amino acids used as anaplerotic contributors to the TCA cycle (i.e., glutamate and the branched-chain amino acids leucine, isoleucine, and valine) were increased and amino acids potentially synthesized as a consequence of glycolysis were decreased (serine, alanine, and threonine) (Additional file 5). These results may indicate perturbations in glycolysis and an overall slowing of energy production.
TCA cycle intermediates and energetics
β-Oxidation, cholesterol synthesis, and bile acids
3-Hydroxybutyrate (BHBA) and carnitines were significantly increased at Tc,Max (Additional file 6), supporting an increase in fatty acid β-oxidation for energy production. Increased β-oxidation is also implied from a trend toward lower plasma levels of the medium-chain fatty acids (e.g., palmitoylcarnitine, stearoylcarnitine, and oleoylcarnitine) relative to controls. These fatty acids can enter the mitochondria without requiring transporter activity (Additional file 6). Further, hexanoylcarnitine, a long-chain fatty acid conjugated to carnitine to facilitate entry into the mitochondria , is increased in heated animals. Increased β-oxidation may increase levels of acetyl coenzyme A (acetyl CoA) (Additional file 1, Additional file 6A), but the profound down-regulation in glycolysis, TCA cycle, and amino acid metabolism suggests that the acetyl CoA is not being used for adenosine triphosphate (ATP) production at Tc,Max.
β-Oxidation may also be up-regulated to facilitate cholesterol synthesis and corticosterone production . At Tc,Max, corticosterone was significantly increased in heated animals relative to controls, and the increase persisted to 24–48 hours (Additional file 6). The increase in BHBA may reflect an increase in acetyl CoA for subsequent cholesterol metabolism and steroidogenesis, but functional mitochondria are required for this process. Alternatively, the elevated corticosterone could be released from the adrenal cortex during heat stress as part of a generalized stress response. Studies suggest prophylactic treatment with glucocorticoids protects against heat stress . Changes in corticosterone levels could gauge severity of perturbations during heat and recovery .
Bile acids and steroids are generated from cholesterol. Following heat stress, plasma levels of steroids increased with a concomitant decrease in bile acids, suggesting de novo steroidogenesis (Additional file 7). Bile acid concentration is usually low in the normal systemic circulation. With heat stress, levels of nearly all bile acids detected were even lower at Tc,Max (Additional file 7), possibly due to decreased reuptake within the small intestine and constricted intestinal circulation .
In addition to the broad categories of energetics and cell death, metabolites of Vitamin B6 are altered by heat stress (Additional file 8). Heat stress resulted in an accumulation of pyridoxate (a breakdown product secreted in the urine), potentially reflecting renal dysfunction. Pyridoxal levels were increased in heated animals relative to controls at Tc,Max, possibly due to increased catabolism. In contrast, pyridoxal and pyridoxate levels were lower than control at 24–48 hours, possibly due to increased anabolism (a process that requires B6 as a cofactor) and return to homeostasis (Additional file 8).
Data integration—random forest analysis to identify biochemicals which discriminate heat exposure from unheated controls
Identification of predictive indicators of persistent cardiac injury at 48 hours
Review of 15 metabolites differentiating animals with histopathological evidence of cardioinflammation from uninjured animals
Identified biochemical (pathway)
Literature review (biochemical and/or pathway)
Ornithine (phenylalanine and tyrosine metabolism)
inhibition of nitric oxide; relationship between nitric oxide modulation of the Frank-Starling response in heart; nitric oxide and nitric oxide synthase are sensitive to thermal stress in fish
3-(4-Hydroxyphenyl) propionate (phenylalanine and tyrosine metabolism)
Biological nitrification inhibition (in plants); phenylalanine and tyrosine concentrations are reduced after Hsp70 increase and heat stress (in yeast)
N1-Methyladenosine (purine metabolism, adenine containing)
N1-methyladenosine analogues are cardioprotective agents in ischemic reperfusion model; decreased infarction; purine metabolism associated with myocardial steatosis and down-regulation of adipose triglyceride In heart
Xylonate (nucleotide sugars, pentose metabolism)
Deficiency in pentose metabolism produces a protective effect through decreased cholesterol synthesis, superoxide production, and reductive stress
1-Stearoylglycerol (1-Monostearin) (monoacylglycerol)
Associated with increased lipid catabolism and remodeling mitochondrial oxidation to aerobic glycolysis (hepatocellular carcinoma)
Carnitine (carnitine metabolism)
Disrupted carnitine metabolism is associated with mitochondrial dysfunction and increased pulmonary flow (lamb model); cardioprotective by increasing heat shock protein synthesis in adriamycin-induced cardiomyopathy
Taurodeoxycholate (bile acid metabolism)
Bile acids exert a protective effect after ischemic injury in porcine hearts; cause endoplasmic reticulum mitochondrial stress; deoxycholate and taurodeoxycholate affect heart mitochondria by decreasing respiration, affecting membrane potential, inducing mitochondrial permeability transition, and altering mitochondrial bioenergetics; impaired cardiac mitochondrial function may cause cardiac alterations in cholestasis
Deoxycholate (bile acid metabolism)
Trigonelline (N’-Methylnicotinate) (nicotinate and nicotinamide metabolism)
Cardioprotective effects after isoproterenol induced myocardial dysfunction (reduction in Hsp27, αB-crystallin and calcium/calmodulin dependent kinase-II-δ)
Diisopropanolamine (xenobiotics - chemical metabolism)
Increases choline uptake without affecting phospholipid synthesis (Chinese hamster ovary cells)
X – 17502 (Unknown)
X – 12419 (unknown)
X – 15808 (unknown)
X – 12131 (unknown)
X – 15651 (unknown)
The results of the random forest analyses at each time point were most concordant with the two-way ANOVA analysis at 24 hours (two of the seven biomarkers identified in both analyses—7-HOCA and corticosterone). At Tc,Max, only corticosterone was common to both analyses, and at 48 hours, only pyroxidate was common to both analyses. It is important to note, however, that 88 of the 90 biochemicals that were excellent classifiers of exposure were also significantly altered according to the ANOVA analyses (as indicated by the red and green circles in Figures 6, 7, 8). Thus, rather than the manual and somewhat arbitrary selection of the top seven biochemical based on circulating levels, the random forest analysis pulls out a list of chemicals that in combination (whether up- or down-regulated) provide greater discriminatory power than a single biomarker, which may be a more powerful approach in populations that demonstrate significant individual variability in the metabolic response across time points.
Comparison of metabolic networks enriched in heat injury with complementary proteomics and transcriptomics pathway analysis
In this work, we present a plasma metabolomic profile of heat stress and heat injury in an in vivo rodent model of heat stress. The metabolomic profile in plasma is concordant with tissue proteomics and transcriptomic profiles indicating energy crisis and oxidative stress, suggesting that metabolic indicators in the plasma may provide surrogate markers for tissue injury in an accessible biofluid. Coupled with the earlier work in tissue transcriptomics and proteomics in the heat-stressed rodent model , the global metabolomic profiles identified the basis for future work for modeling the response to heat stress and heat injury . Integration of metabolomics, proteomics, and transcriptomics in a top-down manner will provide the foundation for network analysis and computational based experiments in a human 3D thermoregulation model . Anchoring these systemic stress responses to the physiological model of heat stroke will provide further insight into forecasting risk of disease, timing of disease onset, and intensity of disease at the organ level. As a result, we anticipate that such a model would accelerate the development of tools to improve disease prevention, classification, and ultimately treatment.
Materials and methods
In vivo rat experiments were performed at the US Army Research Institute of Environmental Medicine (USARIEM). The Institutional Animal Care and Use Committee approved all experimental procedures, which were performed in accordance with the American Physiological Society’s guiding principles for research involving animals and adhered to the high standard (best practice) of veterinary care as stipulated in the Guide for Care and Use of Laboratory Animals. 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. As previously described , male Fischer 344 (F344; n = 48; Charles River Laboratories, Stone Ridge, NY) rats weighing 234–336 g (~2–3 months old) were used. Briefly, rats were housed under standard laboratory conditions (22°C, 12:12 hour light:dark cycle, lights on at 6:00 AM) in an Association for Assessment and Accreditation of Laboratory Animal Care-accredited facility. A relatively cool housing temperature was chosen to support survival during heat-stress recovery , during which chow (Harlan Teklad, LM-485; Madison, WI) and water were provided 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, and mean arterial pressure (±3 mmHg). Physiological response and temperature regulation after heat stress has been published in recent studies by our laboratories ,,.
All experiments were conducted in conscious, free-moving animals, as previously described . Briefly, rats were placed in a floor-standing incubator (Thermo Scientific, Ashville, NC) set at room temperature (RT, 22°C) 24 hours prior to initiation of heat stress experiments. Non-heated rats were not introduced to the incubator environment. Heat-stress experimentation was initiated after Tc of the control and experimental rats reached values of approximately 37.3°C. The following day, food and water were removed from the animal cages; the incubator temperature was increased to 37.0 ± 0.2°C; and the rats were heated until a Tc of 41.8°C (Tc,Max) was reached, at which time they were removed from the incubator, weighed, placed in a new cage, and returned to normal housing temperature (22.0 ± 0.2°C). Animals in the 24 and 48 hour cohorts were allowed to recover for the specified time prior to termination. Time-matched control rats experienced the same experimental procedures as the heat-stressed rats, but remained at the normal housing temperature of 22.0°C throughout experimentation. Control and experimental animals were provided food and water ad libitum throughout recovery. Plasma, heart, liver, lung, and kidney were harvested at Tc,Max, 24, or 48 hours with Tc, mean arterial pressure, and heart rate monitored continuously throughout recovery.
At necropsy, tissues (heart, liver, lung, kidney) were fixed in paraformaldehyde, mounted, sectioned, and stained with hematoxylin and eosin (IHC World, Woodstock, MD). Twenty serial sections were cut per tissue, with three to five sections per slide. Inflammatory or degenerative lesions were graded on a scale of 1 to 5 (grade 1, minimal; grade 5, severe) by a board-certified pathologist (Experimental Pathology Laboratories® [Sterling, VA]). The results of the histopathological evaluation have been previously published .
Metabolomic analysis was conducted by Metabolon (Durham, NC). Briefly, frozen plasma samples (150 μl) were thawed, and extracts were prepared according to Metabolon’s standard protocol, which is designed to remove protein, dislodge small molecules bound to protein or physically trapped in the precipitated protein matrix and recover a wide range of chemically diverse metabolites. Samples were extracted and split into equal parts for analysis on the gas chromatography mass spectrometer and liquid chromatography mass spectrometer platforms. Proprietary software was used to match ions to an in-house library of standards for metabolite identification and for metabolite quantitation by peak area integration.
Two-way ANOVA with contrasts was used to analyze the data (See Additional file 1). For all analyses, missing values (if any) were inputted with the observed minimum for that particular compound (inputted values were added after block-normalization). The statistical analyses were performed on natural log-transformed data to reduce the effect of any potential outliers in the data. Two-way ANOVA and contrast comparisons were made between the means of each biochemical from the groups and were calculated using either or both of the statistical analysis software programs Array Studio (Omicsoft, Inc.) or ‘R’ (R Foundation for Statistical Computing, Vienna, Austria). Statistical cut-offs are typically used to detemine physiological significance in metabolomics studies. Conservative criteria of p < 0.05 and q < 0.1 are routinely used in metabolomic studies , allowing for the identification of significantly altered responses between groups with a false discovery rate of no more than 10%. Because we analyzed a constellation of metabolites in conjunction with biochemical pathways rather than single analytes, we considered p < 0.05 significant, regardless of the q value. Using this approach, we could be more inclusive of data that did not meet strict cut-off values, taking both p-value and number of metabolites changing in a given biochemical pathway.
Random forest, a supervised classification technique based on an ensemble of decision trees , was used to determine the predictive value of multiple biochemicals for both exposure and health effect. For a given decision tree, a random subset of the data with identifying true class information was selected to build the tree (“bootstrap sample” or “training set”), and then the remaining data, the OOB variables, were passed down the tree to obtain a class prediction for each sample, then repeated thousands of times to produce the forest. The final classification of each sample was determined by computing the class prediction frequency for the OOB variables over the whole forest. Class predictions were compared to the true classes, generating the “OOB error rate” as a measure of prediction accuracy. To determine which variables (biochemicals) made the largest contribution to the classification, a variable importance measure was computed, termed the mean decrease accuracy (MDA). The MDA was determined by randomly permuting a variable, running the observed values through the trees, and then reassessing the prediction accuracy. If a variable is not important, then this procedure had little change in the accuracy of the class prediction (permuting random noise gave random noise). By contrast, if a variable is important to the classification, the prediction accuracy will drop after such a permutation, which we record as the MDA. Thus, the random forest analysis provided an importance rank-ordering of biochemical, and the top 30 biochemicals were reported.
- acetyl CoA:
Acetyl coenzyme A
Analysis of variance
γ-glutamyl cysteine synthase
γ -glutamyl transferase
Glutathione disulfide, oxidized
Kyoto encyclopedia of genes and genomes
Mean decrease accuracy
Nicotinamide adenine dinucleotide phosphate
Reactive oxygenated species
System inflammatory response syndrome
US Army Center for Environmental Health Research
US Army Medical Research and Materiel Command
US Army Research Institute of Environmental Medicine
We would like to thank Dr. Roy Vigneulle (Military Operational Medicine Research Program), CAPT Carroll D. Forcino (Director, Military Operational Medicine Research Program), Dr. David A. Jackson (Director, Pulmonary Health, USACEHR), and LTC Thomas C. Timmes (current Commander, USACEHR) for their programmatic support, encouragement, and insightful discussion. We also thank Dr. Bryan Helwig, J. Ward, S. Dineen, M. Blaha, and R. Duran for access to biofluids and technical support with the rat heat stress experiments. We thank Dr. Andrea Eckhart and Dr. Rob Mohney from Metabolon, Inc., for the report and data analysis of the metabolomics data.
- Armed Forces Health Surveillance C: Update: heat injuries, active component, U.S. Armed Forces, 2013. Msmr. 2014, 21 (3): 10-13.Google Scholar
- Wallace RF, Kriebel D, Punnett L, Wegman DH, Wenger CB, Gardner JW, Gonzalez RR: The effects of continuous hot weather training on risk of exertional heat illness. Med Sci Sports Exerc. 2005, 37 (1): 84-90. 10.1249/01.MSS.0000150018.90213.AA.View ArticlePubMedGoogle Scholar
- Wallace RF, Kriebel D, Punnett L, Wegman DH, Wenger CB, Gardner JW, Kark JA: Risk factors for recruit exertional heat illness by gender and training period. Aviat Space Environ Med. 2006, 77 (4): 415-421.PubMedGoogle Scholar
- Armed Forces Health Surveillance Center: Update: heat injuries, active component, U.S. Armed Forces, 2011. MSMR 2012, 19(3):14–16.,
- Armed Forces Health Surveillance Center: Update: heat injuries, active component, U.S. Armed Forces, 2012. MSMR 2013, 20(3):17–20.,
- Merrill CT, Miller M, Steiner C: Hospital stays resulting from excessive heat and cold exposure due to weather conditions in U.S. Community Hospitals, 2005: statistical brief #55. In: Healthcare Cost and Utilization Project (HCUP) Statistical Briefs. Agency for Health Care Policy and Research (US): Rockville (MD); 2006.
- United States Army: Army Regulation 40–501: Standards of Medical Fitness. 2011Google Scholar
- Liu ZF, Sun XG, Tang J, Tang YQ, Tong HS, Wen Q, Liu YS, Su L: Intestinal inflammation and tissue injury in response to heat stress and cooling treatment in mice. Mol Med Rep. 2011, 4 (3): 437-443.PubMedGoogle Scholar
- Leon LR: The thermoregulatory consequences of heat stroke: are cytokines involved?. J Therm Biol. 2006, 31 (1–2): 67-81. 10.1016/j.jtherbio.2005.11.023.View ArticleGoogle Scholar
- Leon LR, Helwig BG: Heat stroke: role of the systemic inflammatory response. J Appl Physiol. 2010, 109 (6): 1980-1988. 10.1152/japplphysiol.00301.2010.View ArticlePubMedGoogle Scholar
- Leon LR, Helwig BG: Role of endotoxin and cytokines in the systemic inflammatory response to heat injury. Front Biosci. 2010, 2: 916-938. 10.2741/S111.View ArticleGoogle Scholar
- Zeller L, Novack V, Barski L, Jotkowitz A, Almog Y: Exertional heatstroke: clinical characteristics, diagnostic and therapeutic considerations. Eur J Intern Med. 2011, 22 (3): 296-299. 10.1016/j.ejim.2010.12.013.View ArticlePubMedGoogle Scholar
- Zhou FH, Song Q, Peng ZY, Pan L, Kang HJ, Tang S, Yue H, Liu H, Xie F: Effects of continuous venous-venous hemofiltration on heat stroke patients: a retrospective study. J Trauma-Injury Infect Crit Care. 2011, 71 (6): 1562-1568. 10.1097/TA.0b013e31822a71c2.View ArticleGoogle Scholar
- Grosman B, Shaik OS, Helwig BG, Leon LR, Doyle FJ: A physiological systems approach to modeling and resetting of mouse thermoregulation under heat stress. J Appl Physiol. 2011, 111 (3): 938-945. 10.1152/japplphysiol.00519.2010.View ArticlePubMedGoogle Scholar
- Leon LR, Dineen SM, Clarke DC: Early Activation of liver apoptotic signaling pathways during heat stroke recovery in mice. Faseb J 2013, 27:1201.7,
- Leon LR, Duran RM, Helwig BG: Complementing heat stroke: activation and amplification of the IL-6 and complement system during heat stroke recovery in F344 rats. Faseb J 2013, 27:1201.6.,
- Leon LR, Eustis HL, Urso ML: Skeletal muscle is a potential source of cytokines during heat stroke recovery in mice. Faseb J 2012, 26:1084.17.
- 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, 115 (12): 1822-1837. 10.1152/japplphysiol.00589.2013.View ArticlePubMedGoogle Scholar
- Stallings JD, Ippolito DL, Rakesh V, Baer CE, Dennis WE, Helwig BG, Jackson DA, Leon LR, Lewis JA, Reifman J: Patterns of gene expression associated with recovery and injury in heat-stressed rats. BMC Genomics. 2014, 15 (1): 1058-10.1186/1471-2164-15-1058.View ArticlePubMedPubMed CentralGoogle Scholar
- Malmendal A, Overgaard J, Bundy JG, Sorensen JG, Nielsen NC, Loeschcke V, Holmstrup M: Metabolomic profiling of heat stress: hardening and recovery of homeostasis in Drosophila. Am J Physiol Regul Integr Comp Physiol. 2006, 291 (1): R205-R212. 10.1152/ajpregu.00867.2005.View ArticlePubMedGoogle Scholar
- Quinn CM, Duran RM, Audet GN, Charkoudian N, Leon LR: Cardiovascular and thermoregulatory biomarkers of heat stroke severity in a conscious rat model. J Appl Physiol. 2014, 117 (9): 971-978. 10.1152/japplphysiol.00365.2014.View ArticlePubMedGoogle Scholar
- Stasolla C, Loukanina N, Yeung EC, Thorpe TA: Alterations in pyrimidine nucleotide metabolism as an early signal during the execution of programmed cell death in tobacco BY-2 cells. J Exp Bot. 2004, 55 (408): 2513-2522. 10.1093/jxb/erh259.View ArticlePubMedGoogle Scholar
- Piskur J, Schnackerz KD, Andersen G, Bjornberg O: Comparative genomics reveals novel biochemical pathways. Trends Genet. 2007, 23 (8): 369-372. 10.1016/j.tig.2007.05.007.View ArticlePubMedGoogle Scholar
- Lu Z, Zhang R, Diasio RB: Dihydropyrimidine dehydrogenase activity in human peripheral blood mononuclear cells and liver: population characteristics, newly identified deficient patients, and clinical implication in 5-fluorouracil chemotherapy. Cancer Res. 1993, 53 (22): 5433-5438.PubMedGoogle Scholar
- van Lenthe H, van Kuilenburg AB, Ito T, Bootsma AH, van Cruchten A, Wada Y, van Gennip AH: Defects in pyrimidine degradation identified by HPLC-electrospray tandem mass spectrometry of urine specimens or urine-soaked filter paper strips. Clin Chem. 2000, 46 (12): 1916-1922.PubMedGoogle Scholar
- Kolker S, Okun JG, Horster F, Assmann B, Ahlemeyer B, Kohlmuller D, Exner-Camps S, Mayatepek E, Krieglstein J, Hoffmann GF: 3-Ureidopropionate contributes to the neuropathology of 3-ureidopropionase deficiency and severe propionic aciduria: a hypothesis. J Neurosci Res. 2001, 66 (4): 666-673. 10.1002/jnr.10012.View ArticlePubMedGoogle Scholar
- Kellogg DL, Crandall CG, Liu Y, Charkoudian N, Johnson JM: Nitric oxide and cutaneous active vasodilation during heat stress in humans. J Appl Physiol. 1998, 85 (3): 824-829.PubMedGoogle Scholar
- Ray PD, Huang BW, Tsuji Y: Reactive Oxygen Species (ROS) homeostasis and redox regulation in cellular signaling. Cell Signal. 2012, 24 (5): 981-990. 10.1016/j.cellsig.2012.01.008.View ArticlePubMedPubMed CentralGoogle Scholar
- Hall DM, Baumgardner KR, Oberley TD, Gisolfi CV: Splanchnic tissues undergo hypoxic stress during whole body hyperthermia. Am J Physiol. 1999, 276 (5 Pt 1): G1195-G1203.PubMedGoogle Scholar
- Wu G, Fang YZ, Yang S, Lupton JR, Turner ND: Glutathione metabolism and its implications for health. J Nutr. 2004, 134 (3): 489-492.PubMedGoogle Scholar
- Lim CL, Wilson G, Brown L, Coombes JS, Mackinnon LT: Pre-existing inflammatory state compromises heat tolerance in rats exposed to heat stress. Am J Physiol Regul Integr Comp Physiol. 2007, 292 (1): R186-R194. 10.1152/ajpregu.00921.2005.View ArticlePubMedGoogle Scholar
- Meaney S, Babiker A, Lutjohann D, Diczfalusy U, Axelson M, Bjorkhem I: On the origin of the cholestenoic acids in human circulation. Steroids. 2003, 68 (7–8): 595-601. 10.1016/S0039-128X(03)00081-3.View ArticlePubMedGoogle Scholar
- Moccand C, Boycheva S, Surriabre P, Tambasco-Studart M, Raschke M, Kaufmann M, Fitzpatrick TB: The pseudoenzyme PDX1.2 boosts vitamin B6 biosynthesis under heat and oxidative stress in Arabidopsis. J Biol Chem. 2014, 289 (12): 8203-8216. 10.1074/jbc.M113.540526.View ArticlePubMedPubMed CentralGoogle Scholar
- Trachootham D, Lu W, Ogasawara MA, Nilsa RD, Huang P: Redox regulation of cell survival. Antioxid Redox Signal. 2008, 10 (8): 1343-1374. 10.1089/ars.2007.1957.View ArticlePubMedPubMed CentralGoogle Scholar
- Nelson DL, Lehninger AL, Cox MM: Lehninger Principles of Biochemistry. W. H. Freeman and Company New York, NY; 2008.
- Colleluori DM, Ash DE: Classical and slow-binding inhibitors of human type II arginase. Biochemistry. 2001, 40 (31): 9356-9362. 10.1021/bi010783g.View ArticlePubMedGoogle Scholar
- Derave W, Everaert I, Beeckman S, Baguet A: Muscle carnosine metabolism and beta-alanine supplementation in relation to exercise and training. Sports Med. 2010, 40 (3): 247-263. 10.2165/11530310-000000000-00000.View ArticlePubMedGoogle Scholar
- Evans DA: N-acetyltransferase. Pharmacol Ther. 1989, 42 (2): 157-234. 10.1016/0163-7258(89)90036-3.View ArticlePubMedGoogle Scholar
- Ezaki J, Matsumoto N, Takeda-Ezaki M, Komatsu M, Takahashi K, Hiraoka Y, Taka H, Fujimura T, Takehana K, Yoshida M, Iwata J, Tanida I, Furuya N, Zheng DM, Tada N, Tanaka K, Kominami E, Ueno T: Liver autophagy contributes to the maintenance of blood glucose and amino acid levels. Autophagy. 2011, 7 (7): 727-736. 10.4161/auto.7.7.15371.View ArticlePubMedPubMed CentralGoogle Scholar
- Gathiram P, Wells MT, Brockutne JG, Gaffin SL: Prophylactic corticosteroid increases survival in experimental heat-stroke in primates. Aviat Space Environ Med. 1988, 59 (4): 352-355.PubMedGoogle Scholar
- Hofmann AF, Molino G, Milanese M, Belforte G: Description and simulation of a physiological pharmacokinetic model for the metabolism and enterohepatic circulation of bile acids in man. cholic acid in healthy man. J Clin Invest. 1983, 71 (4): 1003-1022. 10.1172/JCI110828.View ArticlePubMedPubMed CentralGoogle Scholar
- Amelio D, Garofalo F, Capria C, Tota B, Imbrogno S: Effects of temperature on the nitric oxide-dependent modulation of the Frank-Starling mechanism: the fish heart as a case study. Comp Biochem Physiol A Mol Integr Physiol. 2013, 164 (2): 356-362. 10.1016/j.cbpa.2012.10.037.View ArticlePubMedGoogle Scholar
- Waagner D, Heckmann LH, Malmendal A, Nielsen NC, Holmstrup M, Bayley M: Hsp70 expression and metabolite composition in response to short-term thermal changes in Folsomia candida (Collembola). Comp Biochem Physiol A Mol Integr Physiol. 2010, 157 (2): 177-183. 10.1016/j.cbpa.2010.06.171.View ArticlePubMedGoogle Scholar
- Zakir HA, Subbarao GV, Pearse SJ, Gopalakrishnan S, Ito O, Ishikawa T, Kawano N, Nakahara K, Yoshihashi T, Ono H, Yoshida M: Detection, isolation and characterization of a root-exuded compound, methyl 3-(4-hydroxyphenyl) propionate, responsible for biological nitrification inhibition by sorghum (Sorghum bicolor). New Phytol. 2008, 180 (2): 442-451. 10.1111/j.1469-8137.2008.02576.x.View ArticlePubMedGoogle Scholar
- Inoue T, Kobayashi K, Inoguchi T, Sonoda N, Maeda Y, Hirata E, Fujimura Y, Miura D, Hirano K, Takayanagi R: Downregulation of adipose triglyceride lipase in the heart aggravates diabetic cardiomyopathy in db/db mice. Biochem Biophys Res Commun. 2013, 438 (1): 224-229. 10.1016/j.bbrc.2013.07.063.View ArticlePubMedGoogle Scholar
- Kasiganesan H, Wright GL, Chiacchio MA, Gumina G: Novel l-adenosine analogs as cardioprotective agents. Bioorg Med Chem. 2009, 17 (14): 5347-5352. 10.1016/j.bmc.2008.12.011.View ArticlePubMedGoogle Scholar
- Hecker PA, Leopold JA, Gupte SA, Recchia FA, Stanley WC: Impact of glucose-6-phosphate dehydrogenase deficiency on the pathophysiology of cardiovascular disease. Am J Physiol Heart Circ Physiol. 2013, 304 (4): H491-H500. 10.1152/ajpheart.00721.2012.View ArticlePubMedPubMed CentralGoogle Scholar
- Beyoglu D, Imbeaud S, Maurhofer O, Bioulac-Sage P, Zucman-Rossi J, Dufour JF, Idle JR: Tissue metabolomics of hepatocellular carcinoma: tumor energy metabolism and the role of transcriptomic classification. Hepatology. 2013, 58 (1): 229-238. 10.1002/hep.26350.View ArticlePubMedPubMed CentralGoogle Scholar
- Sharma S, Aramburo A, Rafikov R, Sun X, Kumar S, Oishi PE, Datar SA, Raff G, Xoinis K, Kalkan G, Fratz S, Fineman JR, Black SM: L-carnitine preserves endothelial function in a lamb model of increased pulmonary blood flow. Pediatr Res. 2013, 74 (1): 39-47. 10.1038/pr.2013.71.View ArticlePubMedPubMed CentralGoogle Scholar
- Strauss M, Anselmi G, Hermoso T, Tejero F: Carnitine promotes heat shock protein synthesis in adriamycin-induced cardiomyopathy in a neonatal rat experimental model. J Mol Cell Cardiol. 1998, 30 (11): 2319-2325. 10.1006/jmcc.1998.0793.View ArticlePubMedGoogle Scholar
- Sun X, Sharma S, Fratz S, Kumar S, Rafikov R, Aggarwal S, Rafikova O, Lu Q, Burns T, Dasarathy S, Wright J, Schreiber C, Radman M, Fineman JR, Black SM: Disruption of endothelial cell mitochondrial bioenergetics in lambs with increased pulmonary blood flow. Antioxid Redox Signal. 2013, 18 (14): 1739-1752. 10.1089/ars.2012.4806.View ArticlePubMedPubMed CentralGoogle Scholar
- Dai BH, Geng L, Wang Y, Sui CJ, Xie F, Shen RX, Shen WF, Yang JM: microRNA-199a-5p protects hepatocytes from bile acid-induced sustained endoplasmic reticulum stress. Cell Death Dis. 2013, 4: e604-10.1038/cddis.2013.134.View ArticlePubMedPubMed CentralGoogle Scholar
- Ejiri S, Eguchi Y, Kishida A, Ishigami F, Kurumi Y, Tani T, Kodama M: Cellular distribution of thrombomodulin as an early marker for warm ischemic liver injury in porcine liver transplantation: protective effect of prostaglandin I2 analogue and tauroursodeoxycholic acid. Transplantation. 2001, 71 (6): 721-726. 10.1097/00007890-200103270-00006.View ArticlePubMedGoogle Scholar
- Ferreira M, Coxito PM, Sardao VA, Palmeira CM, Oliveira PJ: Bile acids are toxic for isolated cardiac mitochondria: a possible cause for hepatic-derived cardiomyopathies?. Cardiovasc Toxicol. 2005, 5 (1): 63-73. 10.1385/CT:5:1:063.View ArticlePubMedGoogle Scholar
- Gao X, Fu L, Xiao M, Xu C, Sun L, Zhang T, Zheng F, Mei C: The nephroprotective effect of tauroursodeoxycholic acid on ischaemia/reperfusion-induced acute kidney injury by inhibiting endoplasmic reticulum stress. Basic Clin Pharmacol Toxicol. 2012, 111 (1): 14-23.PubMedGoogle Scholar
- Panda S, Biswas S, Kar A: Trigonelline isolated from fenugreek seed protects against isoproterenol-induced myocardial injury through down-regulation of Hsp27 and alphaB-crystallin. Nutrition. 2013, 29 (11–12): 1395-1403. 10.1016/j.nut.2013.05.006.View ArticlePubMedGoogle Scholar
- Stott WT, Kleinert KM: Effect of diisopropanolamine upon choline uptake and phospholipid synthesis in Chinese hamster ovary cells. Food Chem Toxicol. 2008, 46 (2): 761-766. 10.1016/j.fct.2007.10.001.View ArticlePubMedGoogle Scholar
- Khatri P, Sirota M, Butte AJ: Ten years of pathway analysis: current approaches and outstanding challenges. PLoS Comput Biol. 2012, 8 (2): e1002375-10.1371/journal.pcbi.1002375.View ArticlePubMedPubMed CentralGoogle Scholar
- Rakesh V, Stallings JD, Reifman J: A virtual rat for simulating environmental and exertional heat stress. J Appl Physiol 2014, 117(11):1278–1286.,
- Davies B, Morris T: Physiological parameters in laboratory animals and humans. Pharm Res. 1993, 10 (7): 1093-1095. 10.1023/A:1018943613122.View ArticlePubMedGoogle Scholar
- Bartel J, Krumsiek J, Theis FJ: Statistical methods for the analysis of high-throughput metabolomics data. Comput Struct Biotechnol J. 2013, 4: e201301009-PubMedPubMed CentralGoogle Scholar
- Breiman L: Using iterated bagging to debias regressions. Mach Learn. 2001, 45 (3): 261-277. 10.1023/A:1017934522171.View ArticleGoogle Scholar
- Goldstein BA, Hubbard AE, Cutler A, Barcellos LF: An application of random forests to a genome-wide association dataset: methodological considerations & new findings. BMC Genet. 2010, 11: 49-10.1186/1471-2156-11-49.View ArticlePubMedPubMed CentralGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.