Rat models for low and high adaptive response to exercise differ for stress‐related memory and anxiety

Abstract Physical exercise and fitness may serve as resilience factors to stress exposure. However, the extreme range in human exercise performance suggests that genetic variation for exercise capacity could be a confounding feature to understanding the connection between exercise and stress exposure. To test this idea, we use laboratory rat models selectively bred for a low and high gain in aerobic running capacity in response to training to examine whether an inherent capacity to respond to physical exercise reflects how stress changes neurobiological functioning and regulates fear‐associated memory processing. Utilization of this contrasting rat model system of low and high responders has the potential to guide the interpretation of the reported association with exercise involvement and the reduction of stress‐induced anxiety disorders. Our data show that aerobic fitness may be linked to the ability to regulate fear‐associated memories. We also show that acquired exercise capacity may play a key role in regulating responses to an acute stressor. Exercise sensitivity plays a significant role in the activation of the plasticity‐associated molecule extracellular signal‐regulated kinase, changes in stress hormone activity, and anatomical modifications to the noradrenergic locus coeruleus. These data identify a unique operational mechanism that may serve as translational targets for lessening symptoms of stress and anxiety.

2 of 10 | VANDERHEYDEN Et Al. psychological stress, exercise is generally accepted as being healthy for the mind and body by increasing neurogenesis (Praag et al., 1999) and improving learning and memory (Ahmadiasl et al., 2003). In fact, exercise may be a useful therapy for anxiety and stress disorders Blumenthal et al., 2007).
Physical fitness is hypothesized to be a resilience factor in the acquisition of stress-related disorders Zen et al., 2012). However, the neural and molecular mechanisms leading to this type of resilience are poorly understood. To assess the molecular and anatomical regulation of fear-associated memory processing, we examined the response to fear conditioning, fear extinction, and fear extinction recall in strains of rats that differ in their adaptive response to aerobic treadmill training. We developed high response trainer (HRT) and low response trainer (LRT) rats by divergent selective breeding for gain in exercise capacity after performing the same amount of treadmill training . Compared to HRT rats, LRT rats fail to increase VO 2max , show blunted cardiac muscle remodeling response, have altered signal transduction for JNK and p38 mitogen-activated protein kinase, reduced mitochondrial biogenesis regulating factors in skeletal muscle (PGC1-α, NRF1, and TFAM), and lower hippocampal neurogenesis in response to aerobic exercise training suggesting impaired plasticity (Lessard et al., 2013;Marton et al., 2014;Nokia et al., 2016;Wisloff et al., 2015). It also appears that exercise training differently modified the results on a passive avoidance test used to evaluate learning and memory. The HRT rats had a decline in recall after gaining an improvement in exercise capacity whereas the LRT rats were not altered (Marton et al., 2016). Thus, the motivation for this study was to discover whether this animal model system of contrasting sensitivity to exercise training may be uniquely suited for studying the genotype-to-phenotype connection between physical exercise, stress, and the processing of fear-associated memories.
In addition to behavioral responses to fear conditioning, we also examined the neuro-anatomical regulation of plasticity-related molecules in response to an acute 2-hour physical restraint stressor in LRT and HRT animals. Extracellular signal-regulated kinase (ERK), glucocorticoid receptor (GR), tyrosine hydroxylase (TH), and brain-derived neurotrophic factor (BDNF) have been identified to be regulated by both exercise (Holmes, 2014;Rothman & Mattson, 2013;Widegren et al., 2001) and stress (Sciolino & Holmes, 2012;Takei et al., 2011). We performed these experiments in order to assess whether key molecular contributions to the fear memory processing centers of the brain differ between LRT and HRT rats. We chose to assess changes in brain areas previously identified to be contributing to aberrant stress responses and fear memory processing abnormalities (reviewed in Liberzon & Sripada, 2008;Pitman et al., 2012;Sripada et al., 2012;Vanderheyden et al., 2014).
We hypothesize that inherited variation for aerobic gain in response to physical training may play a role in mediating plasticity-related responses that affect stress, memory, and anxiety pathways. Due to the differential activation of plasticity-related molecules and their importance in mediating gene transcription, the data herein suggest that fear-associated memory impairments may be due to altered epigenetic mechanisms in a circuit dependent manner or by structural changes in the noradrenergic system. Further experimentation on these animals, however, is required to determine the specific epigenetic mechanisms mediating the interaction between physical fitness and stress exposure. Last, we show that the HRT and LRT rats are a unique animal model system that can be used to understand the relationships between the processing of fear-associated memories, and the molecular and anatomical interactions of stress and physical exercise.

| General animal procedures
All animal procedures were carried out in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals and with approval from the University of Michigan Committee on the Use and Care of Laboratory Animals (UCUCA). The LRT and HRT selected lines are developed by Koch and Britton and currently maintained in a specific pathogen-free facility at the University of Toledo . Starting with a founder population of genetically heterogeneous rats (N/ NIH stock), the two-way artificial selection was performed across several generations based on the magnitude of change in running capacity after completing 8 weeks of

Highlights
• Rats selectively bred for low and high response to physical exercise training show a differential ability to extinguish or recall fear-associated memories. • High and low responders to exercise training show the differential activation of the stress hormone system. • Plasticity-associated molecules are differentially regulated by stress in these selectively bred animals. • Rats bred to be more sensitive to exercise training results in structural anomalies in the noradrenergic locus coeruleus. standardized aerobic treadmill training. Animals are phenotyped for the response to training starting at 3 months of age and completed by 6 months of age. For this study, a cohort of 20-month-old males from generation 20 of selection was housed individually in sound-attenuating boxes on a 12:12 hr, light:dark schedule at constant temperature (23°C) and humidity (20%) and given ad libitum access to rat chow and water.

| Fear conditioning, fear extinction, and extinction recall
Fear conditioning experiments were conducted using eight High Response to Training (HRT) and eight Low Response to Training (LRT) mature male rats. Fear conditioning, extinction, and extinction recall were performed as previously published (experiment 3 from Knox et al. [2012]). All fear conditioning, extinction, and extinction recall experiments were performed in four identical Coulbourn Instruments Rat Test Cages (12"W × 10"D × 12"H) (Whitehall, PA) containing a Shock Floor with 18 current-carrying metal bars, a wall-mounted speaker and in-chamber lighting. Test cages were housed in wooden sound-attenuating boxes. Tones were delivered via speakers mounted in the housing of the test cages and controlled by FreezeFrame data acquisition software (Coulbourn Instruments). Shocks were delivered through precision animal shockers (Coulbourn Instruments) also controlled by FreezeFrame software. Ceiling mounted cameras recorded behavior for analysis and FreezeFrame was used to assess freezing levels.
As previously published (Knox et al., 2012), two unique contexts were created using two different sets of olfactory and visual cues. Context A consists of 50 ml of 1% acetic acid solution placed in a small dish above the test cage and standard lighting which illuminates the chamber walls of the Rat Test Cages. Context B consists of 50 ml of a 1% ammonium hydroxide solution placed in a small dish above the test cage along with patterned paper placed on the chamber walls to alter the visual context. Other labs, as well as our own, have used these specific methods repeatedly in the past and found no evidence of increased stress with these concentrations of acetic acid or ammonium hydroxide, that is, they show no effects on behavior or HPA response (Knox et al., 2012).
Fear conditioned animals were exposed to five, 1 mA, 1-s foot-shocks paired with the cessation of a 10 s 80 dB tone in Context A. The first tone was presented 180 s after the animal was placed in the test cage and the subsequent tones occurred with a 60-s inter-tone interval. Sixty seconds after the last tone, animals were removed to their home cages. Fear extinction was conducted 24 hr after fear conditioning and was performed in the distinctly different Context B. Fear extinction consisted of 180 s acclimation to the new context and presentation of 30 ten-second tones without the paired foot-shock with each tone followed by a 60-s inter-tone interval. Extinction recall was assessed 24 hr after extinction and consisted of the animals being placed back into the same fear extinction context (Context B) for 180 s acclimation followed by 10 tones (60-s inter-tone interval), again without foot shock. The percent time spent immobile (freezing) within each 70-s long block (the 10-s tone and 60-s inter-tone interval combined) was assessed by setting threshold values of movement (number of pixels which moved between frames) via FreezeFrame Software. Threshold values were verified to be similar to experimenter-confirmed immobility times.

| Nociceptive testing
Rats (n = 8 HRT and 8 LRT) were conditioned to handling and to the IITC Model 336 T Paw Stimulator Analgesia Meter with a temperature-controlled glass floor (IITC Model 400 Heated Base; IITC Life Science, CA, USA) for one week prior to nociceptive testing. Rats were allowed 30 min to habituate to the experiment chambers before nociceptive testing began. Nociceptive measures started at 2:00 PM and were conducted by an investigator blinded to the rat strain. Thermal hyperalgesia was assessed using the Hargreaves' paw withdrawal latency (PWL) method (Hargreaves et al., 1988). As described in detail previously (Vanini et al., 2014), to test for thermal nociception a light beam was aimed directly at the plantar surface of a hind paw. The light was then switched from idle to the active intensity and a timer was simultaneously started at the onset of the stimulus. An automatic cut-off time for the thermal stimulus was set at 15 s to prevent tissue damage. Upon withdrawal of the paw away from the thermal source, the light beam and the timer were stopped and the latency to withdrawal in seconds was recorded. Four withdrawal latency values were obtained from each paw by alternating the stimulus between the hind paws. A minimum of 30 s was allowed between thermal stimuli. No significant difference in PWL was observed between the left and right paw. The mean PWL from each rat was averaged to compare nociceptive thresholds between HRT and LRT rats.

| Acute stress
Animals (4 LRT and 4 HRT) were restrained in custom-built Plexiglas restraining devices for 2 hr beginning at ZT0. At ZT 2, animals were briefly anesthetized in a glass bell jar containing vaporized isoflurane. Biological tissues (blood, brain) were collected in under 2 min and quickly frozen and stored in a −80 freezer for later processing.

| Western blot
Frozen brains were sectioned in 100 μm on a −20°C, Leica (Buffalo Grove, IL) cryostat and specific brain areas (amygdala, hippocampus, and locus coeruleus [LC]) were dissected into cold micro-centrifuge tubes using a metal 2 mm circular tissue biopsy punch. Tissues from individual animals were homogenized in 300 μl of Homogenization Buffer (200 mM Tris, 10 mM EDTA, 10 mM Na 4 P 2 O 7 , 100 mM PMSF, 1x Sigma Phosphatase inhibitor II, III [Sigma Aldrich, St Louis, MO]). Total protein was quantified using BC Assay (Bio-Rad, Hercules, CA) and all samples were equalized to the same concentration in Homogenization Buffer. Samples were combined with Laemmli Sample Buffer (Bio-Rad), heated to 100°C for 5 min and then centrifuged at max speed for 3 min and loaded (35 μl) on a gel (4-15% TGX [Bio-Rad, Hercules, CA]). Gel was electrophoresed at 100v for 60 min and then transferred to PVDF membrane at 4°C on ice at 100v for 1.5 hr. Blots were probed with mouse anti-ppERK, 1:1000 (Sigma Aldrich), rabbit anti-total ERK antibodies 1:1000 (Sigma Aldrich), rabbit anti-GR, 1:1000 (Santa Cruz Biotechnology, Dallas, TX), rabbit anti-BDNF antibodies, 1:1000 (Santa Cruz Biotechnology). Blots were visualized using a Li-cor infrared detector and quantified using ImageJ software (NIH).

| Stress hormone assay
Trunk blood was collected rapidly after decapitation and allowed to coagulate at room temperature for 30 min before being centrifuged at 1,000 g for 10 min. Serum was collected into clean micro-centrifuge tubes and rat stress hormone levels were detected using the MILLIPLEX Rat Stress Hormone Magnetic Bead Panel (Millipore, St. Charles, MO.) according to the manufacturer's directions.

| Immunostaining
Animals were sacrificed at ZT2, immediately following a 2-hr physical restraint stressor. Four percent Paraformaldehyde fixed and frozen brain tissue was sectioned in 30 µm increments through the LC on a Leica (Buffalo Grove, IL) cryostat, stored in 1x PBS with sodium azide, and later processed for immunohistochemistry. Briefly, all sections were washed 3x for 10 min in 1x PBS, permeabilized in PBS-Triton (0.5%), blocked in 5% normal goat serum (in PBS-Triton 0.5%), and incubated in primary antibody overnight (24 hr) at 4°C. Sections were then washed in PBS and incubated in secondary antibody (24 hr) at 4°C. Tissue was mounted on poly-lysine coated microscope glass and coverslipped under Prolong Gold (Molecular Probes). The following antibodies were used: mouse anti-TH (ImmunoStar Hudson, WI), rabbit anti-GR (Santa Cruz, Dallas, TX), and Alexa 488 conjugated anti-mouse and Alexa 594 conjugated anti-rabbit IgG (Molecular Probes) at 1:1000. Fluorescent images were collected on an Olympus fluorescent microscope provided by the University of Michigan Microscopy and Imaging core facility (University of Michigan). Immuno-positive GR cells were counted using the ImageJ binary thresholding algorithm (San Diego Plugin). The average number of cells was compared between the LRT and HRT animals. The average values for each group were then evaluated using an independent sample t-test.

| Statistical analysis
Statistical comparison of time spent freezing on the fear conditioning, extinction, and extinction recall tasks was made between HRT and LRT selected rat strains using a 2-way (time × group) ANOVA. Average freezing levels between groups were analyzed independently on fear conditioning, extinction, and extinction recall days.
Western Blot data were quantified using ImageJ software (NIH), normalized to control protein levels and statistical comparisons were made using GraphPad Prism Software.
One-tailed and two-tailed Mann-Whitney tests were performed on Western Blot, pain sensitization PWL, and hormone assay data where appropriate.

| Fear extinction and recall is elevated in HRT animals
In order to test for differences in fear-associated memory between LRT and HRT rats, we performed a 3-day fear conditioning study (detailed in Methods). ANOVA of % freezing during the fear conditioning session (Day 1, Figure 1a) revealed a significant main effect of time (F (5,100) = 87.3, p < 0.0001). There was no significant effect of genotype (F (1,20) = 0.011, p = 0.92) or genotype × time interaction (F (5,100) = 0.567, p = 0.73) (Figure 1a). ANOVA of cued freezing during fear extinction (Day2, Figure 1b animals (Figure 1b,c, respectively) suggests that HRT rats have diminished the processing of fear-associated memories compared to LRT rats.

| Thermal nociception is not different between HRT and LRT rats
To test whether diminished fear extinction and recall may have been due to heightened pain sensitivity in HRT animals, we exposed eight HRT and eight LRT rats to thermal nociception testing as described in the methods section. A one-tailed Mann-Whitney test revealed no significant differences in PWL to nociceptive thermal stimuli between HRT (M = 7.43 ± 0.83 s) and LRT (M = 7.08 ± 1.16 s) (p = 0.25) rats ( Figure 2).

| HRT and LRT rats display a differential hormonal response to acute stress
The exaggerated freezing responses in HRT rats on extinction and fear recall testing days suggest that HRT rats may have an altered hormonal response to stress. To test this hypothesis, we subjected four HRT and four LRT animals to 2 hr of physical restraint stress and immediately anesthetized, decapitated, and collected trunk blood for serum stress hormone detection. Serum ACTH, Corticosterone, and Melatonin levels were measured. Serum ACTH levels were significantly elevated in HRT rats (M = 3.17 ± 2.0 pg/ml) compared to LRT rats (M = 0.76 ± 0.19 pg/ml) (p = 0.0001) (Figure 3a). Serum Corticosterone (p = 0.21) and Melatonin (p = 0.94) levels were not significantly different between HRT and LRT rats (Figure 3b,c).

| Hippocampal ERK activation
Following blood collection for the assessment of hormonal changes in response to stress, we also collected brain tissue of neural circuits involved in fear memory processing (hippocampus, amygdala, and LC) and using Western blot, assessed changes in the plasticity-related molecules pERK, GR, and BDNF. Two-tailed Mann-Whitney tests revealed significant differences in hippocampal pERK activation between LRT (M = 1.0 ± 0.144) and HRT (M = 1.43 ± 0.059) F I G U R E 1 Extinction and Extinction Recall Deficits. On 3 successive days, we performed fear conditioning, extinction, and extinction recall. Freezing values (y-axis) are plotted for the LRT (black triangles) and HRT (grey circles) rats against tones delivered (x-axis). Fear conditioning consisted of one 3-min baseline period followed by the presentation of 5, 1 mA foot shock presented at the cessation of a 10-s 80 dB tone, in Context A (60-s inter-tone interval) (a). Fear extinction followed 24 hr after fear conditioning. All animals received 30 tones in Context B, on the same inter-tone interval (15 tones shown). HRT rats froze significantly more than their LRT comparators (b). Fear extinction recall followed 48 hr after fear conditioning. All animals received 10 tones in the now extinguished Context B, on the same 60-s inter-tone interval. Again, HRT rats froze significantly more than their LRT comparators (c)  (Figure 4)a-c. No statistically significant differences were found in GR levels in the hippocampus, amygdala, or LC between LRT and HRT rats (p = 0.06, 0.11, and 0.65, respectively) ( Figure 5). No statistically significant differences were found in BDNF levels in the hippocampus, amygdala, or LC between LRT and HRT rats (p = 0.99, 0.20, and 0.49, respectively) ( Figure 6).

| Differences in molecular and structural characteristics of the Locus Coeruleus
Based on the data shown here and previously published research on stress, we hypothesized that the noradrenergic LC may differentially undergo molecular (Li et al., 2011;Weiss et al., 1994) or structural changes (Bracha et al., 2005) between LRT and HRT rats in response to stress. Immunohistochemistry for TH, which specifically stains the noradrenergic cells of the LC, revealed a differential anatomical localization of TH containing cells between the LRT and HRT rats. The dorsal LC lacked staining of TH in the HRT animals that are present in LRT animals and normally present in wild-type animals (Circles in Figure 7a,b), while the dorsal section of the LC looked relatively similar between groups. Using Western Blot, we further quantified the change in TH between the LRT and HRT animals. A one-tailed Mann-Whitney test revealed significant differences in TH levels between HRT (M = 1.00 ± 0.04) and LRT (M = 1.15 ± 0.12) (p = 0.028) rats (Figure 7c). We additionally confirmed that this change in TH was not due to differences in cell number within the LC, as the response to physical restraint resulted in F I G U R E 3 Stress hormone activation. The hormones ACTH, corticosterone, and melatonin were quantified after a brief 2-hr physical restraint stress. Compared to LRT rats, HRT rats displayed a significant increase in ACTH (a). A non-significant reduction in Corticosterone levels was found in HRT rats relative to the LRT rats (b). Melatonin activation was not different between LRT and HRT rats (c) F I G U R E 5 Glucocorticoid receptor. Western blot was performed to quantify GR levels following 2 hr of physical restraint. Western blots for GR from hippocampus, amygdala, and LC tissues (a). Quantification was performed using ImageJ software. No statistically significant differences were found for GR between LRT and HRT rats (b) F I G U R E 6 Brain-derived neurotrophic factor. Western blot was performed to quantify BDNF levels following 2 hr of physical restraint. Western blots for BDNF from hippocampus, amygdala, and LC tissues (a). No statistically significant differences in BDNF levels were found between LRT and HRT rats (b)

| DISCUSSION
Physical fitness has been hypothesized to be a resilience factor to stress exposure. Additionally, physical activity has the potential to assist in the treatment of stress-related and anxiety disorders ). Yet, a mechanistic understanding of the association between exercise activity and stress is lacking. The response to exercise training is noted as having wide heterogeneity such that some individuals improve exercise capacity more than others (e.g., higher rise in maximal oxygen consumption) (Vollaard et al., 2009), largely due to unidentified genotype-environment interactions (Bouchard et al., 2012). This work is the first of its kind to identify molecular and anatomical pathways that overlap in their function for stress regulation and adaptation to exercise. In this work, we have shown that, compared to LRT rats, HRT rats show impaired fear-associated memory processing following fear conditioning. HRT rats have been artificially selected to show an enhanced ability to respond to aerobic training. However, the selection process consisting of repeated exposure to forced treadmill exercise at increasing speed and duration, designed to select animals' responsiveness to running training, may have additionally co-selected traits for heightened arousal and stress sensitivity. Exercise training by these methods modulates the dose-response phenomenon that balances the health-promoting effects of exercise with the exercise-induced production of oxidative stress markers (Radak et al., 2017) Indeed, aberrant ACTH signaling arises as a consequence of exposure to stress (Knox et al., 2012). Further, similar to the HRT animals shown here, ERK activation has been shown to be increased in the amygdala when presented with extreme stress (Liu et al., 2010;Xiao et al., 2011). Therefore, the selection process for exercise responsiveness has generated a contrasting animal system that has molecular and behavioral characteristics similar to animals more (HRT) and less responsive (LRT) to extreme stress.
The HRT and LRT rats showed no differences in nociception ( Figure 2). Therefore, the behavioral learning differences detected during fear extinction and fear extinction recall cannot be attributed to an increased perception of pain during fear conditioning. More than likely, these fear-associated memory variances are due to molecular signaling or anatomical remodeling abnormalities. Indeed, pERK activation was greater in HRT rats in the hippocampus and the amygdala following 2 hours of physical restraint stress (Figure 4), and the LC showed significant histochemical changes which may contribute to behavioral differences in response to this fear conditioning paradigm.
Increased noradrenergic signaling contributes to features of hyper-arousal such as those seen in HRT rats with reduced abilities to extinguish fear-associated memories (Geracioti et al., 2001;Mellman et al., 1995;Southwick et al., 1993Southwick et al., , 1999Strawn & Geracioti, 2008). To examine the contribution of the noradrenergic LC in mediating this behavior, we exposed animals to a 2-hr physical restraint stress and assessed differences in GR and TH, a rate-limiting step for catecholamine synthesis. The structural abnormalities in TH localization in the LC revealed a reduced number of LC cells and suggest that some of the behavioral impairments may be due, in part, to structural abnormalities in the noradrenergic LC. Loss of noradrenergic cells is also seen in human populations exposed to extreme levels of physical or psychological stress (Bracha et al., 2005). These data further connect the stress response system to exercise physiology.
Acute stress can provide useful information on how the brain responds to over-activity of some molecular systems and underactivity of others (summarized in Figure 8). This study provides a framework by which the molecular response F I G U R E 7 HRT rats have reduced tyrosine hydroxylase levels. Immuno-histochemical analysis revealed that TH positive cells in the LC were significantly reduced in the HRT rats (b) compared to the LRT (a) rats. Histochemical analysis was verified using Western blot, confirming that TH levels were quantifiably different between LRT and HRT rats (c). This reduction in TH was not due to a reduction in cell numbers in the LC, the number of GR containing cells in the LC remained the same between LRT and HRT animals (d). to acute stress can be examined within the context of exercise. These data, showing that ERK activation is unique between HRT and LRT animals, identify a unique molecular pathway underlying stress and exercise neurobiology. Future work may focus on the contribution of ERK activation in physical exercise as a potential intervention for stress-induced maladaptive behaviors such as fear-memory impairments and hyper-arousal. Additional examination of the BDNF pathway is also warranted given that exercise is sufficient to activate BDNF in the brains of LRT animals exposed to exercise, but not in HRT animals (Marton et al., 2016). As a first report, the data presented here show that rats selectively bred for high adaptation for exercise training have a higher level of molecular and structural brain adaptations in response to stress and a higher level of retention for fear-related memory compared to those rats selected for a lessened adaptation. Further work is required to assess if the HRTspecific fear-associated memory impairments are the result of increased threat avoidance behavior or increased threat vigilance (Waters & Kershaw, 2015) and to understand the mechanisms regulating impaired fear-associated memory processing in these animals in a more biologically comprehensive cross-sectional and longitudinal studies. This study was uniquely conducted in aged animals which lend itself to further efforts needed to assess the association between stress disorders and impaired cognitive performance (Yaffe et al., 2010).
Last, exercise and stress activate many similar anatomical and molecular pathways. For example, both exercise (Jacks et al., 2002) and stress (McEwen, 2000) have been shown to increase cortisol response via the ACTH stimulus. Further complicating the story is the evidence that exercise type, intensity, and duration along with genetic predisposition for high capacity for aerobic training is critical to the beneficial response in the brain for adult hippocampal neurogenesis (Nokia et al., 2016). There is also the idea that the practice of extreme sports with high adrenaline and cortisol release such as whitewater rafting, rock climbing, and skydiving may in turn help to replace bad memories previously linked with high adrenaline and cortisol levels in combat veterans (xsports4vets.org). Indeed, physical activity is prescribed to help reduce symptoms of post-traumatic stress disorder (PTSD) and promote recovery and resilience. However, a non-linear and U-shaped association showing veterans reporting no weekly exercise and daily exercise having the highest prevalence of probable PTSD is consistent with endpoints of the exercise-associated hormesis curve (Adams et al., 2020;Radak et al., 2008). The importance of exercise capacity as a biomarker in the onset of stress-related disorders is currently unknown and the mechanisms of overlapping stress hormone molecular pathways are still poorly understood (McEwen, 2015). It is our contention, therefore, that the HRT and LRT rats we developed may provide a differential genetic substrate to aid in our understanding of the complex biological process by which stress and physical exercise interact.

ACKNOWLEDGMENTS
This work was supported in part by the Office of Research Infrastructure Programs grant (P40OD021331) to LGK & SLB from the National Institutes of Health. The rat models of low and high response to exercise are maintained as an international resource at The University of Toledo, Toledo, OH. Contact Lauren.Koch2@UToledo.Edu for LRT and HRT rats available for study.