Introduction

In the past, there have been numerous studies vastly covering the effects of physical exercise (PE) and environmental enrichment (EE) on affective and cognitive behaviors, locomotion and spatial learning, and memory at various age points, but a comparative study over the life span has never been undertaken. PE, which broadly includes forced running on a treadmill or voluntary running on wheels in cages for rodents, has been shown to improve locomotion (Gentile et al., 1987), enhance spatial memory acquisition, retention, and retrieval (Van Praag et al., 2005; Chen et al., 2006; Van der Borght et al., 2007; Ahmadiasl et al., 2003; Albeck et al., 2006; Alomari et al., 2013; Anderson et al., 2000), augment cognitive abilities (Radák et al., 2001; Nichol et al., 2009; Griesbach et al., 2009; Ma et al., 2017), and alleviate anxiety- and depressive-like behaviors (Binder et al., 2004; Duman et al., 2008; Salam et al., 2009; Helmich et al., 2010; Huang et al., 2017). Physical exercise also has been shown to alter molecular biology, such as monoamine metabolism (Ahmadiasl et al., 2003; Samorajski et al., 1987) and neurotrophin levels in the hippocampus (Alomari et al., 2013; Griesbach et al., 2009; Fang et al., 2013; Ferris et al., 2007; Gómez-Pinilla et al., 2002). Evidence suggests that improvement in hippocampal plasticity and neurogenesis is a prerequisite for change in brain functions in response to PE (Van Praag et al., 2005; Chen et al., 2006; Nichol et al., 2009; Ma et al., 2017).

In rodents, the external environment can be enriched with a variety of toys and novel objects, and this has been shown to elicit a considerable influence on neuronal plasticity and brain functions during aging (Singhal et al., 2014). Such rodents exhibited higher locomotion (Gentile et al., 1987; Van de Weerd et al., 1994), improved visuospatial attention performances and recognition memory (Falkenberg et al., 1992; Harati et al., 2011), enhanced cognition (Kobayashi et al., 2002), and decreased anxiety- and depressive-like behaviors (Richter et al., 2013; Brenes et al., 2008) in preclinical trials. Also, EE has been shown to enhance neurogenesis and neuronal survival (Birch et al., 2013), increase synaptic plasticity and synaptogenesis (Birch et al., 2013; Nakamura et al., 1999), and improve neurotrophin (Falkenberg et al., 1992; Birch et al., 2013), monoamine (Brenes et al., 2008; Rasmuson et al., 1998), and glucocorticoid (Olsson et al., 1995) metabolism in the hippocampus of rodents.

Many studies in rodents investigated the effects of either PE or EE, and not their combination, on behavior, memory, neurobiology, and underlying molecular biology (Gentile et al., 1987; Van Praag et al., 2005; Radák et al., 2001; Samorajski et al., 1987; Van de Weerd et al., 1994; Falkenberg et al., 1992; Rasmuson et al., 1998; Hoffmann et al., 1987; Bowling et al., 1993). Most studies on EE in recent time, however, also have utilized running wheels as a tool of enrichment (Duffy et al., 2001; Morley-Fletcher et al., 2003; Leggio et al., 2005; Segovia et al., 2006; Cui et al., 2006; Yang et al., 2006; Benaroya-Milshtein et al., 2004; Novkovic et al., 2015). The rationale behind the inclusion of running wheels in the EE paradigms is still unclear, and this may have been a source of bias in such behavioral studies on EE. Nevertheless, like PE and EE, the combination of PE and EE (PE+EE) has been shown to improve cognition, spatial learning and memory (Kobayashi et al., 2002; Leggio et al., 2005; Segovia et al., 2006; Frick et al., 2003; Frick & Fernandez, 2003; Hüttenrauch et al., 2016), and alleviate anxiety- and depressive-like behaviors (Cui et al., 2006; Yang et al., 2006; Benaroya-Milshtein et al., 2004; Hüttenrauch et al., 2016; Sampedro-Piquero et al., 2013). Neurobiological and molecular changes in the brain, for example, enhanced hippocampal neurogenesis (Segovia et al., 2006) and dendritic growth (Leggio et al., 2005), improved neurotransmission (Segovia et al., 2006) and HPA axis reactivity (Morley-Fletcher et al., 2003; Laviola et al., 2004), and increased BDNF levels in the brain (Novkovic et al., 2015) also have been reported after PE+EE.

The abovementioned functional change in the brain after PE, EE, and PE+EE may involve change in gene expression, particularly in the hippocampus, because the effects of EE and PE on the brain results primarily from a change in hippocampal integrity (Van der Borght et al., 2009; Pham et al., 1999; Gobbo & O'Mara, 2004). Indeed, it has been reported that external environment influences gene expression and gene-gene interactions, thereby suggesting that modulation of the external environment could be used as a tool to modify gene expression in the brain (Izzotti, 2009). Independent studies have shown that after PE, EE, and PE+EE, genes associated with immune functions, aminergic pathways, neurotrophin metabolism, HPA axis activity, transcription factors, and mitochondrial function express differentially in the brain, thereby improving neural functions and biology (Falkenberg et al., 1992; Rasmuson et al., 1998; Olsson et al., 1995; Hüttenrauch et al., 2016; Keyvani et al., 2004; Mlynarik et al., 2004; Hullinger et al., n.d.; Buschler & Manahan-Vaughan, 2017; Torasdotter et al., 1998; Dahlqvist et al., 1999; Zajac et al., 2010; Tong et al., 2001; Rampon et al., 2000; Kaliman et al., 2011). Any change in gene expression, therefore, can alter underlying neurobiology and may manifest into observable behavioral changes. Indeed, enhanced hippocampal neurogenesis and improved affective-like behaviors after EE have been shown to be the outcomes of an increase in the expression of Ngf and Bdnf (Falkenberg et al., 1992; Torasdotter et al., 1998), and 5-ht1a receptor (Rasmuson et al., 1998) mRNA in the hippocampus of rats. Similarly, PE has been shown to regulate the Sirt1 and Tnf activity, thereby playing a key role in cellular energy metabolism, gene transcription, immune response, and, consequently, in cell survival (Kaliman et al., 2011). The differential effects of PE, EE, and PE+EE on hippocampal gene expression across life-span, however, remain unexplored.

The changes in behavior, memory, neurobiology, cellular and molecular functions, and hippocampal gene expression after PE, EE, and PE+EE may differ. Perhaps, the age and sex of an animal also influence behavioral and neurobiological response after PE, EE, and PE+EE as has been reported previously (Harati et al., 2011; Frick et al., 2003; Frick & Fernandez, 2003; Kim et al., 2004; Lin et al., 2011). It also seems probable that when PE and EE are combined, this could potentiate the individual effects of both treatments. These research questions remain unanswered in the absence of a direct comparison between the three treatments. This study, therefore, examined the differential effects of short-term PE, EE, and the combination of PE and EE on the locomotor activity, affective-like behaviors, and spatial and retention memory at three time-points over the lifespan of healthy male and female C57BL/6 mice. The change in the expression of 43 genes of interest in the hippocampus, which were selected after extensive literature review and included 12 immune (e.g., Il1b, Tnf, Aif1, Gfap), 9 monoamine (e.g., Htr1aa, Grin2a, Tph2), 4 growth factors (e.g., Bdnf, Ngf, Igf1), 5 transcription factor (e.g., Creb1, Smad2, Foxo3), 3 HPA axis (Nr3c1, Crhr1, Crh), and 10 mitochondrial (e.g. Sirt1, Gpx1, Prkaa1) function-related genes also has been investigated in response to the abovementioned treatments. We tested the hypothesis that short-term PE and EE of 4 weeks will reverse the aging-associated impairment in brain functions, and hippocampal genes expression in early (4 months), middle (9 months), and late-middle age (14 months) male and female C57BL/6 mice compared with age-matched control mice and that the combination of short-term PE and EE will have additive beneficial effects than PE and EE alone. Together, this allowed us to investigate the potential underlying mechanisms whereby these three treatments can enhance neuronal function and ensure healthy aging.

Methods

Animals

Wild-type (C57BL/6) mice (n = 158; 80 males and 78 females), parental substrain Nhsd (derived from a colony from the National Institutes of Health), were bred in-house and housed in same-sex groups of four to five in individually ventilated cages (IVCs) in the Laboratory Animal Services facility at the University of Adelaide. In time-period of the experiments, the C57BL/6Nhsd line has inbred between 9-13 generations. The strain of the mice is important as it has been reported that C57BL/6 mice are more reactive and alert to the external environmental stimuli (Van de Weerd et al., 1994). All mice were given ad libitum access to standard laboratory food and water and were maintained under controlled conditions of temperature (21±1 °C), humidity (55 %) and a 12:12-hour dark-light cycle. This study was approved by the University of Adelaide Animal Ethics Committee (M-216-12), and all guidelines as prescribed for handling the experimental animals were followed during the study.

Experimental design

Once of the desired age (3, 8, and 13 months), mice that showed no signs of injury and sickness were distributed equally into three treatment (Physical exercise (PE), Environmental Enrichment (EE), and PE+EE) and control groups (Figure 1). Mice were randomly paired (males and females paired separately) and transferred to open top cages on the morning of Monday, week 1 (2 mice per cage unless fighting necessitated separation). Control mice received no intervention and were kept in cages with the following dimensions: 48.5 cm x 15.5 cm x 12 cm. PE, EE, and PE +EE mice were kept in plexiglass cages with dimensions: 37 cm x 20.5 cm x 13.5 cm, because these had more breadth and depth to provide extra space for running wheel and the objects associated with EE. Mice were inspected daily but handled only once a week while transferring them to clean cages every Friday morning, starting on the Friday of week 1. Friday was selected for change of cages to allow mice to acclimatize to the fresh bedding for 2 days before the Home Cage test is performed after 4 weeks of treatment. At the same time, mice were weighed on a digital weighing scale. Mice were monitored for dominancy throughout the experiments, and those found to be dominant were segregated to prevent dominance effects on behavior and gene expression. Only the EE mice received nesting material (paper shreds) during the experiments.

Figure 1.
figure 1

Methodology used during research. Four weeks of treatments are considered equivalent to 1 month of treatments for the sake of simplicity. IVC = individually ventilated cage; PE = physical exercise; EE = environmental enrichment; M = month; W = week.

The groups assigned as PE and PE+EE were provided with a running wheel for 4 weeks before behavioral testing, i.e., starting when the mice were 3, 8, and 13 months of age. Wheels were made of metal and had the following dimensions: 12 cm in diameter and 5.5 cm wide. Running wheels remained in cages throughout the behavioral testing period. Wheel revolutions were counted using an automated counter on the wheels every Friday during change of cages. At the same time, wheels were cleaned using F10SC veterinary disinfectant. Final readings of the wheel revolutions were taken the day after the behavioral testing ended. The total number of revolutions was then multiplied by the circumference of the wheels, i.e., 37.7 cm, to obtain the total distance traveled. Because the mice were kept in pairs, the total distance traveled for one cage reflects the distance traveled by a pair of mice. Distance traveled was then converted into kilometers for statistical analysis.

The groups assigned as EE and PE+EE were provided with a variety of nontoxic objects (toys, nonchewable novel plastic, and wooden objects, two types of houses — house on the cage floor and house suspended from the open-top roof, colored plastic balls of different texture, hanging toys, ladder, and tunnels) and extra bedding as per previously published protocols (Leggio et al., 2005; Spires et al., 2004; Jankowsky et al., 2005) during the same time period (see Table 1 for the EE protocol that we followed over 8 weeks). The objects and bedding were changed once every week to maintain novelty (during change of cages on Fridays to minimize handling stress to mice). The objects remained in the cages throughout the behavioral testing period.

Table 1. Enrichment protocol over 8 weeks

A priori power analysis of the experimental design using the program G*Power indicated a 98% chance of detecting a large effect size and a 62% chance of detecting a medium effect size (defined by Cohen, 1992, as 0.4 and 0.25 of a population standard deviation between the means respectively (Cohen, 1992)) between the groups significant at the 5% level (ANOVA: fixed effects, main effects, and interactions). We observed that the power of our analyses would be 0.8 at the moderate effect size of 0.3 with n = 158, which we deemed acceptable for research.

The first week consisted of only 5 days, because the EE objects were changed every Friday morning during change of cages.

Behavioral analysis

After 4 weeks of exposure to PE and/or EE, mice underwent a behavioral test battery following established procedures in the laboratory starting at 4, 9, and 14 months of age (4 weeks of treatments are considered equivalent to 1 month of treatments for the sake of simplicity), as per the below schedule:

Week 1: Home cage (Monday), open field (Tuesday)

Week 2: Elevated zero maze (Monday)

Week 3: Barnes maze (Monday to Friday)

Week 4: Forced swim test (Monday)

The behavioral testing was done in order of the least to most stressful tests. All trials were recorded by a ceiling-mounted camera and analyzed using Any-maze software ver. 4.70 from Stoelting, USA. To remove the olfactory traces, F10SC veterinary disinfectant was used to clean the testing area during and between behavioral tests.

Locomotor and Exploratory Behaviour

  1. 1.

    Home cage activity: The general locomotor activity and exploratory behavior of mice were analyzed in home cages with 2-day-old bedding under basal nonstressful conditions according to previously published protocols (Baune et al., 2008). The total distance covered over 5 minutes was recorded as a measure of locomotor activity. Running wheels were removed during the Home cage test, because wheels occupied significant space in the cage and may have affected the home cage activity results. However, mice assigned as EE and PE+EE were left with two random EE objects to keep the environment enriched during the Home cage test.

  1. 2.

    Open Field: Locomotor activity and exploratory behavior of mice under more anxiety-inducing conditions were measured in the Open Field where mice were individually tested in a brightly lit (approximately 600 lumens) square arena of 40 cm x 40 cm, with clear 35-cm high walls, according to previously published protocols (Baune et al., 2008; Hart et al., 2010). The floor was divided into inner (26.6 cm x 26.6 cm) and outer (13.4 cm x 13.4 cm) zones. Total distance traveled in Open Field over 5 minutes was measured as an indication of locomotor activity and exploratory behavior.

Anxiety-like Behaviour

  1. 1.

    Elevated zero maze: The elevated zero maze is a round maze, 105 cm in diameter, with a 5-cm wide platform 40 cm above the floor and divided into four equal quadrants (Braun et al., 2011). Two quadrants have 15-cm high walls (closed), and two are open in alternate order. The EZM was placed in a brightly lit (approximately 600 lumens) area during the behavioral testing. The mouse was placed on the open quadrant and allowed to explore for 5 minutes according to published protocols (Shepherd et al., 1994). The time spent by the mouse in closed and open quadrants and the number of head dips were recorded as a measure of anxiety in mice. Increased time spent in closed as opposed to open quadrants, and lower number of head dips is considered indices of anxiety-like behaviors (Wahlsten, 2010).

  1. 2.

    Open Field: Time spent by mice in the inner and outer zones of the Open Field was recorded as a measure of anxiety-like behavior. Less time spent in the inner zone is an index of anxiety-like behavior (Baune et al., 2008; Hart et al., 2010).

Spatial learning and retention memory

  1. 1.

    Barnes Maze: Change in spatial learning and retention memory was measured using the Barnes Maze, a circular grey platform 91 cm in diameter and elevated 90 cm above the ground, with 20 holes on the perimeter: 1 with a hidden escape box and the rest with false boxes. Because the false boxes are too small to enter and looked the same as the target escape box to the mouse, they removed visual cues that might be observed through an open hole. Barnes Maze procedures were performed according to published protocols (Baune et al., 2008; McAfoose & Baune, 2009; McLay et al., 1998). Time to find the location of the escape box over 4 days of training was recorded. The mice were placed in the centre of the maze and were allowed to locate the escape box for 3 minutes. Three trials separated by 15 minutes were performed for each mouse on a day. Latency to find the location of the escape box was used to assess the spatial learning and visual memory of the mouse in the new environment (O’leary & Brown, 2012).

On day 5, the position of escape box was changed from the original training position (in NW quadrant) to the probe trial position (in NE quadrant – at 90 degrees). Time to find the location of the escape box in NW quadrant was recorded to assess the retention of spatial memory of the mouse in the new environment (O’leary & Brown, 2012).

We also calculated the entry errors, i.e., searches of any holes that did not contain the hidden escape box beneath it following the protocol published by Nithianantharajah et al. (Nithianantharajah et al., 2008). This included nose pokes and head deflections over a hole. Increased number of entry errors were considered indices of spatial learning and memory impairment.

Depressive-like Behaviour

  1. 1.

    Forced swim test: Duration of immobility (when mice were floating with no movement of limbs in any direction) over 6 minutes in a 4-L circular cylinder, 20 cm in diameter and 45-cm high and filled two thirds with water (at 23-24 °C) was recorded using a ceiling-mounted camera as a measure of despair and depression-like behavior (Wahlsten, 2010; Petit-Demouliere et al., 2005; Porsolt et al., 1977).

Gene expression analysis

Mice were lethally injected intraperitoneally with pentobarbital (60 mg/kg IP) a day after forced swim test. Brains were extracted, dissected and hippocampal tissues were stored in RNAlater (Ambion, Life Technologies) at −80 °C until further processing.

Quantification of the levels of mRNA of 43 genes of interest (GOI) across various cytokines, monoamines, neurotrophins, and other genes was performed using TaqMan assays (Life Technologies, ThermoFisher, Australia) in the high-throughput qPCR system BioMarkHD™ (Fluidigm Inc., USA). Briefly, total RNA was extracted from HC tissues stored in RNALater using PureLink RNA mini extraction kit (Ambion) following the manufacturer’s instructions. Total RNA was then subjected to reverse transcription using the SuperScript III first-strand cDNA synthesis system (Invitrogen, Australia) according to the manufacturer’s instructions. The samples were then prepared for High-throughput qPCR in BioMark HD™ (Fluidigm Inc.) using a single 14-cycle pre-amplification consisting of 20 ng of each cDNA sample mixed with pooled TaqMan assays (47 assays, including 4 endogenous reference genes, listed in Table 2) and PreAmp Master Mix (Fluidigm Inc.). Real-time quantitative PCR was performed for each TaqMan assay for each sample in a 96.96 dynamic array nanofluidic chip (Fluidigm Inc.). A total of 47 x 96 (Assays x Samples in duplicates) PCR reactions were performed. Cycle threshold (Ct) values were generated by Fluidigm real-time PCR analysis software (Fluidigm Inc.).

Table 2. Genes of Interest quantified using high-throughput qPCR

Statistical Analysis

In GraphPad Prism version 7.02 (GraphPad Software Inc.), all data outliers were removed using the ROUT method, and normality of data distribution was determined by visual inspection of histograms. Further analyses were performed in SPSS statistics version 25.0 (IBM Corporation). The three-way interactions between age, treatments, and sex for body weight, distance traveled on running wheels, locomotion, anxiety- and depressive-like behaviors, and spatial learning and memory were determined using three-way ANOVAs. Following this, the analyses of simple two-way interaction effects between age and treatments for males and females, simple main effects of the treatments, and simple comparisons between the groups were conducted. If the three-way interaction was significant, then two-way ANCOVA was used to control for sex differences while determining the effects of age and treatment on the dependent variable. If the three-way interaction was nonsignificant, the two-way interaction effect of age and treatment on the dependent variable was determined using a two-way ANOVA. Post-hoc analysis was conducted using Tukey’s multiple comparison test. Results were presented as mean ± standard error (SE). Differences were considered statistically significant when p < 0.05.

Differentially expressed genes were identified by analysis of Ct data measurements taken from the BioMarkHD arrays in R. Briefly, Input expression values (Ct) of house-keeping genes (B2m, Gapdh, Gusb, and Hprt) were compared across all samples to identify outlier samples. Delta Ct values were normalized against the geometric mean of all sample expression values. Linear mixed-effects models were then used to compare normalized expression between control and treatment groups and adjusted for multiple comparisons using the R package multcomp (Hothorn et al.,2008). Genes were identified as differentially expressed between control and EE groups when the adjusted p values were <0.05 with Z scores indicating the direction of expression change. GeneMANIA was used to identify the molecular pathways modulated by the differentially expressed genes at a false discovery rate (FDR) < 0.05 (Warde-Farley et al., 2010). We were not able to analyze the sex differences for gene expression due to the relatively smaller sample size of gene expression data compared with behavioral data. This also limits us from performing the correlation analysis between behavioral and gene expression data.

Results

Animal body weight

Body weights of mice at the start and end of experiments were statistically analyzed using a two-way and three-way ANOVA, respectively. Body weights were measured in grams and are presented as a mean ± standard error unless stated otherwise. No statistically significant interaction between age and sex was observed for body weights at the start of the experiment (F (2, 152) = 1.001, p = 0.370). There was, however, a statistically significant main effect of age on starting body weights (F (2, 152) = 110.741, p < 0.001). The simple pairwise comparison revealed that the mean body weight of 4-month mice was significantly lower than 9- and 14-month-old mice at the start of experiments (21.3 ± 0.3 vs. 26.2 ± 0.3 and 26.7 ± 0.3; p < 0.001). The analysis of body weights taken at the end of experiments using a three-way ANOVA showed a nonsignificant interaction between age, treatment, and sex (F (6, 134) = 1.808, p = 0.102). A two-way ANOVA also revealed a nonsignificant interaction between age and treatment (F (6, 146) = 0.837, p = 0.543), as well as the main effect of treatment (F (3, 146) = 0.904, p = 0.441), but the main effect of age was significant (F (2, 146) = 37.438, p < 0.001). Body weights of four months mice were significantly lower than both 9- and 14-month-old mice in all experimental groups (Table 3; p < 0.001).

Table 3. Mouse body weights at the start and end of the experiment

Distance traveled on wheels during the experimental period

Three-way ANOVA found a statistically significant three-way interaction between age, treatment, and sex for distance traveled on running wheels (F (2, 62) = 3.96, p = 0.024). However, the simple two-way interactions between age and treatment for both males and females were nonsignificant (F (2, 62) = 1.666, p = 0.197, F (2, 62) = 2.832, p = 0.067). Hence, a two-way ANOVA was conducted to analyze the two-way interaction effect of age and treatment on distance traveled on running wheels, but the result was still nonsignificant (F (2, 68) = 0.459, p = 0.634; Fig. 2).

Figure 2.
figure 2

Two-way ANOVA for distance traveled on wheels during the experimental period by separate cohorts of C57BL/6 mice exposed to PE and PE+EE conditions. All data represented as estimated marginal means (EMM), n = 11-14 per group.

Behavior

Locomotor activity

  • Home Cage: Distance traveled

Baseline locomotor activity was measured in the home-cage following 1-month exposure to either control conditions, PE, EE, or PE+EE at 4, 9, and 14 months of age. Distance traveled was measured in meters and is presented as mean ± standard error, unless stated otherwise. Three-way ANOVA found a significant three-way interaction between age, treatment, and sex (F(6, 133) = 2.581, p = 0.021). Further analysis revealed a significant simple two-way interaction between age and treatment for both males (F(6, 133) = 2.320, p = 0.037) and females (F(6, 133) = 4.018, p = 0.001). There was a significant simple main effect of treatment for males at 14 months (F(3, 133) = 6.901, p < 0.001) and females at 4 (F(3, 133) = 5.604, p = 0.001), 9 (F(3, 133) = 4.198, p = 0.007), and 14 (F(3, 133) = 4.163, p = 0.007) months. Hence, simple pairwise comparisons were conducted for males and females with a Bonferroni adjustment applied (Figure 3A). The 14-month male PE+EE mice were significantly more mobile in home cage than 14-month male control mice (14.5 ± 1.2 vs. 8.2 ± 1.2; p = 0.001) and PE mice (14.5 ± 1.2 vs. 8.6 ± 1.1; p = 0.002). Conversely, PE and EE treatments, and not PE+EE, significantly improved the baseline locomotion compared to age matched controls in 14-month female mice (14.5 ± 1.2 and 15.7 ± 1.2 vs. 9.4 ± 1.2; p = 0.013 and 0.001 respectively). EE also improved the baseline locomotion of female mice significantly at 9 month compared to age matched controls (18.1 ± 1.2 vs. 12.4 ± 1.2; p = 0.003) and 4-month EE mice (18.1 ± 1.2 vs. and 12.2 ± 1.2; p = 0.001). Female PE mice at 4 months, however, showed significantly reduced baseline locomtion than 4-month female control, EE, and PE+EE mice (7.7 ± 1.2 vs. 12.2 ± 0.9, 12.2 ± 1.2, and 12.5 ± 1.2; p = 0.016, 0.036, and 0.021 respectively), as well as 9- and 14-month female PE mice (7.7 ± 1.2 vs. 15.2 ± 1 and 14.5 ± 1.2; p < 0.001 for both).

Figure 3.
figure 3

Locomotor activity in Home Cage. (A) Simple pairwise comparisons for distance traveled in the home cage at 4, 9, and 14 months by separate cohorts of male and female C57BL/6 mice exposed to PE, EE, PE+EE or control conditions. All data represented as mean ± SE, n = 6-10 per group. (B) Two-way ANCOVA with sex as a covariate at the following value: 1.4904. All data represented as estimated marginal means (EMM), n = 11-19 per group. *Significant difference between a treatment and age-matched control or between two treatments at one age point. #Significant difference between the matched treatments over two age points. *,#p < 0.05, ** p < 0.01, ***,###p < 0.001.

Following above analysis, we controlled baseline locomotion for sex differences using two-way ANCOVA and found a significant two-way interaction between age and treatment (F (6, 144) = 3.384, p = 0.004; Figure 3B). An analysis of simple main effects for age and treatment with statistical significance receiving a Bonferroni adjustment and being accepted at the p < 0.0125 level for age and p < 0.0167 for treatment was conducted. EE mice were significantly more mobile than control mice at 9 months (16.0 ± 0.876 vs. 11.7 ± 0.785; p = 0.003) and 14 months (13.8 ± 0.876 vs. 8.8 ± 0.876; p < 0.001), and PE+EE mice showed significantly higher locomotion than controls at 14 months (14.0 ± 0.915 vs. 8.8 ± .876; p < 0.001). Also, EE and PE mice at 9 months showed significantly higher locomotion than 4-month treatment-matched cohorts (16.0 ± 0.876 vs. 11.9 ± 0.876; p = 0.004 and 12.7 ± 0.784 vs. 9.1 ± 0.915; p = 0.009).

  • Open Field: Distance traveled

Mice were again tested for locomotor activity in the open field test 1 day after the home cage. Distance traveled was measured in meters and are presented as mean ± standard error unless stated otherwise. A three-way ANOVA found nonsignificant interaction between age, treatment, and sex for locomotion in the open field (F (6, 132) = 1.694, p = 0.127). The simple two-way interactions between age and treatment were also nonsignificant for both male and female mice (F (6, 132) = 1.344, p = 0.242 and F (6, 132) = 1.209, p = 0.306). Hence, no corrections were made for sex differences, and further analysis was performed using a two-way ANOVA. No significant main effects of interaction and age were noted (F (6,144) = 0.782, p = 0.586 and F (2,144) = 0.245, p = 0.783); however, the main effect of treatment was found to be significant (F (3,144) = 5.807, p = 0.001; Figure 4). At 14 months, PE+EE-treated mice traveled significantly less distance than controls (10.5 ± 1.2 vs. 14.7 ± 0.9 m; p = 0.0216). At 9 months, both PE and PE+EE treated mice showed significantly lower locomotion than EE treated mice (12.7 ± 1.1 vs. 15.8 ± 0.8 m; p < 0.05; p = 0.0228 and 10.6 ± 0.7 vs. 15.8 ± 0.8 m; p = 0.0002).

Figure 4.
figure 4

Locomotor activity in Open Field. Two-way ANOVA for distance traveled in the open field at 4, 9, and 14 months by separate cohorts of C57BL/6 mice exposed to PE, EE, PE+EE or control conditions. All data represented as estimated marginal means (EMM), n = 11-19 per group. *Significant difference between a treatment and age-matched control or between two treatments at one age point. * p < 0.05, *** p < 0.001.

Anxiety-like behavior

  • EZM: Time in open arms

Time spent in the open quadrants of EZM is a measure of anxiety-like behavior, with less anxious mice spending more time in the open quadrants. Time was measured in seconds and is presented as a mean ± standard error unless stated otherwise. A three-way ANOVA showed statistical significance for the three-way interaction between age, treatment, and sex for time in open arms (F (6, 118) = 2.477, p = 0.027). Further analysis showed a statistically significant simple two-way interaction between age and treatment for males (F(6, 118) = 3.165, p = 0.006), but not for females (F(6, 118) = 1.599, p = 0.153). There was a significant simple main effect of treatment for males at 9 months (F(3, 118) = 3.680, p = 0.014) and females at 4 months (F(3, 118) = 5.211, p = 0.002). Hence, simple pairwise comparisons were run for males and females with a Bonferroni adjustment applied (Figure 5A). Nine-month-old male EE mice were found to spend significantly more time in the open arms of EZM than 4- and 14-month-old male EE mice (64.3 ± 8 vs. 36.8 ± 8 and 14.6 ± 7.3; p = 0.049 and < 0.001 respectively), and 9-month-old male PE mice (64.3 ± 8 vs. 30.7 ± 6.8; p = 0.010). The time in open arms was significantly low for 14-month male PE+EE mice compared to 4- and 9-month-old male PE+EE mice (9.3 ± 7.3 vs. 42.7 ± 8 and 48.3 ± 8; p = 0.007 and 0.001 respectively). For females, 4-month-old control mice spent significantly more time in the open arms of EZM than age-matched PE (52 ± 6.3 vs. 9 ± 8.9; p = 0.001) and 14-month-old control mice (52 ± 6.3 vs. 23.9 ± 7.3; p = 0.013).

Figure 5.
figure 5

Anxiety-like behavior in EZM. (A) Simple pairwise comparisons for time in open arms in the EZM at 4, 9, and 14 months by separate cohorts of male and female C57BL/6 mice exposed to PE, EE, PE+EE or control conditions. All data represented as mean ± SE, n = 6-10 per group. (B) Two-way ANCOVA for time in open arms in the EZM with sex as a covariate at the following value: 1.4904. All data represented as estimated marginal means (EMM), n = 11-19 per group. *Significant difference between a treatment and age-matched control or between two treatments at one age point. #Significant difference between the matched treatments over two age points. *,#p < 0.05, **,##p < 0.01, ###p < 0.001.

Using a two-way ANCOVA, we then controlled time in open arms in the EZM for sex differences and found a nonsignificant two-way interaction effect between age and treatment (F (6, 129) = 2.152; p = 0.052). An analysis of the simple main effects for age and treatment with statistical significance receiving a Bonferroni adjustment and being accepted at the p < 0.0125 level for age and p < 0.0167 for treatments showed non-significant main effect of treatment (F (3, 129) = 1.736, p = 0.163) but significant main effect of age (F (2, 129) = 16.965, p < 0.0001; Figure 5B). The EE mice at 9 months, and PE+EE mice at 4 and 9 months spent significantly more time in open arms than 14-month EE (56.7 ± 5.6 vs. 19.7 ± 5.6; p = 0.049) and PE+EE (38.0 ± 5.6 and 45.2 ± 5.6; p = 0.011 and <0.001) mice respectively. Also, 9-month-old EE mice were found to spend significantly more time in the open arms than 9-month-old PE mice (56.7 ± 5.6 vs. 31.6 ± 5.1; p = 0.007).

  • EZM: Number of Head dips

The number of head dips along the open quadrants in EZM is a measure of anxiety-like behavior, with less anxious mice showing more number of head dips. Data for head dips is presented as a mean ± standard error unless stated otherwise. The three-way ANOVA showed a nonsignificant three-way interaction between age, treatment, and sex for head dips in the elevated zero maze (F (6, 132) = 1.830, p = 0.098). Two-way interactions between age and treatment for both males and females were also nonsignificant (F (6, 132) = 2.083, p = 0.059 and F (6, 132) = .332, p = 0.919). Hence, no corrections were made for sex differences. Follow-up analysis with a two-way ANOVA showed that the main effect of interaction between age and treatments on the number of head dips in EZM was nonsignificant (F (6,144) = 0.453; p = 0.842), but the main effects of treatment and age were significant (F (3,144) = 5.155; p = 0.002 and F (2,144) = 7.367; p = 0.001; Figure 6). At 4 months, EE mice showed a significantly higher number of head dips than PE+EE mice (33.1 ± 2.6 vs. 19.8 ± 1.7; p = 0.0099). Also, the 4-month controls and EE treated mice showed a significantly higher number of head dips than 14-month controls and EE treated mice respectively (25.2 ± 1.7 vs. 15.9 ± 3.4; p = 0.0429 and 33.1 ± 2.6 vs. 21.1 ± 4.4; p = 0.0143).

Figure 6.
figure 6

Anxiety-like behavior in EZM. Two-way ANOVA for number of head dips in EZM. All data represented as estimated marginal means (EMM), n = 11-19 per group. *Significant difference between a treatment and age-matched control or between two treatments at one age point. #Significant difference between the matched treatments over two age points. #p < 0.05, ** p < 0.01.

  • Open Field: Time spent in the inner zone

Time spent in the inner zone of the open field is a measure of anxiety-like behavior, with more anxious mice spending less time in the inner zone. Time was measured in seconds and is presented as mean ± standard error unless stated otherwise. The three-way ANOVA showed a nonsignificant three-way interaction between age, treatment and sex (F (6, 134) = 0.368, p = 0.152). Similarly, the two-way interactions between age and treatment for both males and females were also non-significant (F (6, 134) = 1.195, p = 0.313 and F (6, 134) = 0.202, p = 0.976). Hence, no corrections were made for sex differences. Following this, a two-way ANOVA was utilized to examine the effects of age and treatments on time spent in the inner zone of open field and found that the main effects of interaction and treatment were not significant (F (6,146) = 1.048; p = 0.397 and F (3,146) = 2.540; p = 0.059), but the main effect of age was significant (F (2,146) = 8.805; p < 0.001; Figure 7). Post-hoc analysis using Tukey’s multiple comparison test showed that 9-month-old control mice spent significantly more time in the inner zone than 4-month control mice (80.2 ± 16.2 vs. 57.9 ± 8.1; p = 0.0178). Similarly, 9-month PE mice showed significantly higher time in the inner zone than 4- and 14-month PE mice (127.9 ± 15.9 vs. 75.5 ± 15.5 and 76.0 ± 9.3; p = 0.0124 and 0.0124 respectively).

Figure 7.
figure 7

Anxiety-like behavior in Open Field. Two-way ANOVA for time in the inner zone of Open Field. All data represented as estimated marginal means (EMM), n = 11-19 per group. #Significant difference between the matched treatments over two age points. #p < 0.05.

Depressive-like behavior

  • FST: Immobility time

Mice exposed to PE, EE, PE+EE, or control conditions were tested for depressive-like behavior in the Forced swim test. Time was measured in seconds and is presented as a mean ± standard error unless stated otherwise. The three-way ANOVA showed a nonsignificant three-way interaction between age, treatment, and sex for immobility time in FST (F (6, 131) = 0.627, p = 0.708). The two-way interactions between age and treatment for immobility time in FST also were nonsignificant for both male and female mice (F (6, 131) = 1.005, p = 0.425 and F (6, 131) = 0.658, p = 0.684). As such, no corrections were made for sex. Further analysis with the two-way ANOVA to determine the effects of age and treatments on immobility time in FST showed non-significant main effects of interaction and age (F (6, 143) = 1.052, p = 0.395 and F (2, 143) = 0.816, p = 0.444), but the main effect of treatment was significant (F (3,143) = 5.673; p = 0.001; Figure 8). Post-hoc analysis using Tukey’s multiple comparison test showed that at 4 months PE mice were immobile for a significantly longer time than both controls and EE treated mice (82.6 ± 17.0 vs. 29.7 ± 6.9 and 22.1 ± 6.9; p = 0.0171 and 0.0171).

Figure 8.
figure 8

Depressive-like behavior in FST. Two-way ANOVA for immobility time at 4, 9, and 14 months in separate cohorts of animals exposed to PE, EE, PE+EE or control conditions. All data represented mean ± SEM, n = 12-14 per group. *Significant difference between a treatment and age-matched control or between two treatments at one age point. * p < 0.05.

Spatial learning and memory

  • Spatial learning over 4 days of training in Barnes maze

Spatial learning at 4, 9, and 14 months in separate cohorts of mice exposed to PE, EE, PE+EE, or control conditions was tested by measuring the latency to find the escape box and entry errors over 4 days of training in Barnes maze. Latency to escape was measured in seconds and is presented as a mean ± standard error unless stated otherwise. Data was analyzed using three-way ANOVA, which showed a three-way interaction between age, treatment, and sex for spatial learning over 4 days of training (F (6, 536) = 2.600, p = 0.017). Post-hoc analysis using Tukey’s multiple comparison test showed significant differences between groups at 4 and 14 months for females and 14 months for males (Table 4; Figure 9).

Table 4 Multiple comparisons for the effects of EE, PE, and PE+EE on latency to escape over 4 days of training in Barnes maze in male and female mice at 4, 9, and 14 months
Figure 9.
figure 9

Spatial learning. Latency to find the escape box (LtoE) over 4 days of training at 4 months (A), 9 months (B), and 14 months (C) in separate cohorts of male and female animals exposed to PE, EE, PE+EE, or control conditions in Barnes Maze. All data represented mean ± SEM, n = 5-10 per group. α: Control vs. PE, β: PE vs. EE, γ: PE vs. PE+EE, μ: EE vs. PE+EE. α,β,γ,μ: p < 0.05, ββ,γγ: p < 0.01.

Figure 10.
figure 10

Spatial learning. Entry errors (EnEr) over four days of training at 4 months (A), 9 months (B), and 14 months (C) in separate cohorts of male and female animals exposed to PE, EE, PE+EE, or control conditions in Barnes Maze. All data represented mean ± SEM, n = 5-10 per group. @: Control vs. PE+EE, β: PE vs. EE, γ: PE vs. PE+EE. @,β,γ: p < 0.05.

Three-way ANOVA also showed a significant three-way interaction between age, treatment, and sex for entry errors over 4 days of training in Barnes maze (F (6, 536) = 4.391, p < 0.001). However, post-hoc analysis with Tukey’s multiple comparison tests showed that none of the male and female groups differed from controls on all 4 days of training, except for female 4-month PE+EE mice, which showed significantly lower entry errors than 4-month female control mice on day 3 (9.6 ± 2.8 vs. 33.4 ± 8.8; p = 0.0149). Also, 9- and 14-month male PE mice showed significantly higher entry errors than age-matched male PE+EE (43.7 ± 4.1 vs. 20.8 ± 4.7; p = 0.0334) mice on day 1 and EE mice on day 2 (49.0 ± 9.3 vs. 18.3 ± 5.5; p = 0.0013) respectively (Fig 10).

The sex differences for latency to escape and entry errors over 4 days of training in Barnes maze were controlled using a two-way ANCOVA. We observed that the two-way interaction was significant for latency to escape (F (6, 583) = 3.191, p = 0.004) but nonsignificant for entry errors (F (6, 583) = 0.504, p = 0.805). Post-hoc analysis using Tukey’s multiple comparisons for latency to escape showed significant differences between groups at 4 and 14 months (Table 5; Figure 11).

Table 5 Differentially expressed genes between control, PE, EE, and PE+EE mice at 9 months of age (listed in order of increasing adjusted p values)
Figure 11.
figure 11

Spatial learning. Latency to find the escape box (LtoE) over 4 days of training in Barnes Maze at 4 months (A), 9 months (B), and 14 months (C) in separate cohorts of animals exposed to PE, EE, PE+EE, or control conditions with sex as covariate at the following value: 1.4937. All data represented mean ± SEM, n = 11-19 per group. α: Control vs. PE, β: PE vs. EE, γ: PE vs. PE+EE. α,β,γ: p < 0.05, βββ: p < 0.001.

Table 5. Multiple comparisons for the effects of EE, PE, and PE+EE on latency to escape over 4 days of training in Barnes maze in mice at 4, 9, and 14 months (with sex as a covariate at the following value: 1.4937).

Months

Day

Tukey's multiple comparisons test

Mean comparison

Adjusted p

4

1

Control vs. PE

46.8 ± 5.7 vs. 80.3 ± 8.2

0.0498

3

Control vs. PE

28.6 ± 5.9 vs. 59.7 ± 10.4

0.0282

PE vs. EE

59.7 ± 10.4 vs. 35.7 ± 8.0

0.0368

9

No significant differences between any of the groups were noted.

14

1

PE vs. EE

120.8 ± 9.7 vs. 82.3 ± 11.6

0.0255

PE vs. PE+EE

120.8 ± 9.7 vs. 81.1 ± 11.2

0.0192

2

PE vs. EE

81.8 ± 9.5 vs. 41.0 ± 5.4

0.0182

PE vs. PE+EE

81.8 ± 9.5 vs. 45.8 ± 7.5

0.0369

3

Control vs. PE

57.1 ± 15.3 vs. 95.0 ± 16.8

0.0234

PE vs. EE

95.0 ± 16.8 vs. 37.0 ± 10.4

0.0002

PE vs. PE+EE

95.0 ± 16.8 vs. 56.9 ± 14.0

0.0228

4

Control vs. PE

50.5 ± 13.2 vs. 85.0 ± 15.9

0.0490

PE vs. EE

85.0 ± 15.9 vs. 31.4 ± 8.4

0.0006

PE vs. PE+EE

85.0 ± 15.9 vs. 46.6 ± 13.7

0.0216

  1. All data represented as mean ± SEM, n = 11-19 per group.
  • Retention of spatial learning in Barnes maze probe trial on day 5

After 4 days of training on Barnes maze, mice were tested for retention of spatial memory in a probe trial on day 5. Three-way ANOVA showed nonsignificant three-way interaction between age, treatment and sex for both latency to escape (F (6, 126) = 0.270, p = 0.950) and entry errors (F (6, 134) = 0.739, p = 0.619). The two-way interactions between age and treatment for latency to escape and entry errors during the probe trial were also non significant for both male (F (6, 126) = 0.446, p = 0.846 and F (6, 134) = 1.627, p = 0.144) and female (F (6, 126) = 1.431, p = 0.208 and F (6, 134) = 0.935, p = 0.472) mice. Therefore, the data were analyzed again without controlling for sex differences using a two-way ANOVA.

No significant main effects of interaction and treatment were noted for escape latency, but the main effect of age was significant (F (2,139) = 23.11; p < 0.0001; Figure 12A). Post-hoc analysis showed that at 14-month controls and PE mice had significantly higher latency to find the NW quadrant location of the escape box than EE mice (110.5 ± 21.6 and 107.5 ± 20.4 vs. 21.0 ± 8.3; p = 0.0391 and 0.0406). Similarly, 14-month control, PE and PE+EE mice showed significantly higher latency compared with 4- and 9-month control (110.5 ± 21.6 vs. 3.7 ± 0.7 and 4.1 ± 1.0; p = 0.0002 and < 0.0001), 9-month PE (107.5 ± 20.4 vs. 6.2 ± 1.0; p < 0.0001), and 9-month PE+EE (73.4 ± 21.1 vs 7.7 ± 1.4; p = 0.0140) mice, respectively. Also, 4-month PE mice showed significantly higher latency to find the NW quadrant location of the escape box than 9-month PE mice (68.9 ± 24.3 vs. 6.2 ± 1.0; p = 0.0163).

Figure 12.
figure 12

Retention of spatial learning (A) Escape latency to find the escape box and (B) entry errors during probe trial at 4, 9, and 14 months in separate cohorts of C57BL/6 mice exposed to PE, EE, PE+EE or control conditions. All data represented as (A) mean ± SEM, and (B) estimated marginal means (EMM), n = 11-19 per group. *Significant difference between a treatment and age-matched control or two different treatments at one age point. #Significant difference between the matched treatments over two age points. *,#p < 0.05, ##p < 0.01, ###p < 0.001, ####p < 0.0001.

Using two-way ANOVA, we also noted a significant main effect of age for entry errors (F (2,146) = 8.365; p = 0.0004; Figure 12B). Again, this was because of the 14-month-old EE and PE mice showing significantly higher number of entry errors than 4- and 9-month EE (69.8 ± 10.6 vs. 40.8 ± 11.8 and 30.7 ± 7.3; p = 0.0396 and 0.0054) and PE (63.9 ± 13.8 vs. 26.5 ± 6.2 and 29.3 ± 5.5; p = 0.0085 and 0.0085) cohorts.

Changes in the expression of genes of interest in the hippocampus of C57BL/6 mice exposed to short-term PE, EE or PE+EE

The differential effects of PE, EE, and PE+EE on hippocampal gene expression were investigated with high-throughput qPCR (43 target genes; Table 2). No effects of any of the treatments on hippocampal gene expression, with an adjusted p < 0.05 for significance, were observed at 4 and 14 months of age. However, at 9 months, several genes were differentially expressed amongst the groups. Three genes, i.e., Creb1, Smad2, and Il1r1 differentially expressed during PE vs. control, PE vs. EE, PE+EE vs. control, and PE+EE vs. EE comparisons. Similarly, genes Sirt1 and Il12a were found to be differentially expressed on comparing PE with control, and PE with EE groups. The gene Prkaa2 was found to be differentially regulated during PE vs. control, and PE vs. EE group comparisons. Other genes that were differentially expressed during single comparisons only were Actb, Ifng, Htr1a, Tnfrsf1a, and Foxo3 (PE vs. control), and Htr2a, Cs, and Prkaa1 (PE+EE vs. EE). None of the genes differentially expressed during the remaining two comparisons, i.e., EE vs. C and PE vs. PE+EE groups. The differentially expressed genes between the mice groups at 9 months of age are listed in Table 5.

Discussion

In our study, none of the treatments, i.e., PE, EE, and PE+EE, affected body weight significantly during experiments. Similarly, distance traveled on running wheel during the experiment was not significantly different between the PE and PE+EE mice at all three age points, i.e., early (4 months), middle (9 months), and late-middle (14 months) age. While short-term EE improved locomotor activity, reversed the age-related decline in spatial memory and reduced anxiety-like behavior, short-term PE had no effects on spatial memory and affective-like behaviors in middle and late-middle age mice. In early age mice, short-term PE reduced spatial learning and increased depressive-like behavior, particularly affecting female mice. These results on the effects of PE are interesting, knowing that short-term PE alone and in combination with EE modified the expression of many hippocampal genes of interest in middle-age mice. EE alone showed no effects on the expression of genes of interest in the hippocampus of mice at all ages. Main findings of our research are discussed below pointwise.

Short-term PE, EE, and combined PE and EE do not affect the mean body weight of C57BL/6 mice.

Early evidence suggests that PE reduces body weight by mobilizing fats, both within rodents (Crews 3rd et al., 1969; Brown & Van Huss, 1973; Levin & Dunn-Meynell, 2004) and humans (Donnelly et al., 2000; Messier et al., 2004). The lack of significant changes in body weight of both male and female mice within our study may be attributable to the short duration of the treatment protocols or the young baseline age of mice with minimal body fat to begin with. Our reasoning is consistent with a previous study that reported a decrease in body weight only in response to continuous exercise and not intermittent exercise (Donnelly et al., 2000). Similarly, EE of 6 weeks that included a running wheel in 10-week-old mice resulted in no change in body weight (Benaroya-Milshtein et al., 2004). We did not investigate the body composition of mice or food intake in the study, and recommend this be investigated in similar future studies to get more clarity on the effects of external interventions on body weight. We also theorize that a lengthier exposure to voluntary PE on mice older than late-middle age may show observable changes in body weights.

EE improves baseline locomotion from middle age onwards. PE+EE, although beneficial at late-middle age, reduces locomotion under stress

Short-term EE under non-threatening environmental conditions of the home cage induced increased locomotion at 9 and 14 months of age, whereas PE+EE resulted in increased baseline locomotion at 14 months only. The effects of EE on baseline locomotion is primarily driven by a significant increase in the baseline locomotion of female mice and not male mice at 9 and 14 months. On the contrary, the effects of PE+EE at 14 months were driven by the enhanced baseline locomotion of male mice. Furthermore, 4-month-old female PE mice showed significantly reduced baseline locomotion; however, this seems to be compensated by the baseline locomotion of the 4-month male PE mice, thereby resulting in an overall non-significant difference between 4-month PE and age-matched control PE mice.

Under more threatening environmental conditions in the open field, PE+EE mice were less mobile at 14 months of age compared with controls. This may be due to the greater anxiety-inducing environment of the open field compared with the home cage. While not significantly different to age-matched controls, EE mice outperformed PE and PE+EE mice in exploring the open field at 9 months of age. We observed no sex differences in open field locomotor activity. Our results for locomotor activity comply with the previously published findings. For example, voluntary wheel running of 6 weeks reduced the locomotor activity of 8-weeks-old spontaneously hypertensive male rats in the Open Field test (Hoffmann et al., 1987). EE without running wheel for 39 days has been shown to enhance locomotion in young male Sprague-Dawley rats in the Open Field test (Bowling et al., 1993). Similarly, EE without running wheel for 11 weeks enhanced play behavior in home conditions in newly weaned male mice (Marashi et al., 2004). Overall, our result suggests that EE primarily enhances the locomotor activity of female C57BL/6 mice, but the effects of PE+EE on locomotor activity depends on the presence or absence of stressful conditions. These results contrast with the previously reported findings in male rodents; however, this can be attributed to the effect of different rodent species and strain used (Van de Weerd et al., 1994). Further research into differentiating the effects of PE, EE, and PE+EE on locomotion after a longer duration of treatment in C57BL/6 mice is required.

EE, PE, and PE+EE did not affect anxiety-like behavior compared to age-matched controls at all ages but had an effect across the lifespan

While none of the treatments were different to controls at all ages for anxiety-like behavior, we found age effects within treatments. Short-term EE and PE mice appeared to exhibit less anxiety-like behavior at 9 months compared with their 4- and 14-month cohorts as evident by significantly more time spent in the open quadrants of EZM and inner zone of Open Field respectively. However, sex influenced anxiety-like behavior. We found that EE mice were different due to the significantly less anxiety-like behavior of 9-month male EE mice compared with 4- and 14-month male EE mice. Similarly, PE was found to be anxiogenic for females compared with age-matched controls at 4 months. Indeed, PE, when added to EE at 4 months, increased anxiety, as evident by significantly fewer head dips by 4-month PE+EE mice compared with 4-month EE mice in the EZM. Furthermore, the observed significant higher anxiety-like behavior of 14-month PE+EE mice compared with 9-month counterparts, as evident by differences between the treatment matched age cohorts for time in open arms of EZM, was primarily due to higher anxiety-like behavior in male PE+EE mice. The 14-month controls and EE mice showed significantly less head dips in the EZM than their cohorts at 4 months, suggesting significantly higher anxiety-like behavior at late-middle age compared to early age cohorts in untreated and EE mice during our experiments. Together, our findings point toward EE being anxiolytic at young and middle age, PE being anxiogenic for females at young age and EE when combined with PE may induce anxiety at late-middle age in male C57BL/6 mice.

PE and EE, as well as the combination of PE and EE, are established anxiolytic interventions (Benaroya-Milshtein et al., 2004; Sampedro-Piquero et al., 2013; Fulk et al., 2004; Petruzzello et al., 1991; Görtz et al., 2008). Indeed, both PE and EE seem to positively affect the ventral hippocampus, a region involved in anxiety response (Görtz et al., 2008; Pietropaolo et al., 2014). Short-term PE+EE has shown to reduce anxiety-like behavior at a young age in female transgenic models of AD (Görtz et al., 2008; Pietropaolo et al., 2014), male C3H/eB mice (Benaroya-Milshtein et al., 2004), as well as in healthy male Wistar rats (Sampedro-Piquero et al., 2013). This contradicts our above-reported findings. Similarly, our results also challenge the previously reported findings on PE. For example, voluntary wheel running of 4 weeks when provided to 12- to 14-week-old single house C57BL/6 mice reduced anxiety-like behavior and impulsiveness (Binder et al., 2004). Similarly, 8 weeks old male C57BL/6 mice, which underwent voluntary wheel running of 2 weeks showed reduced anxiety-like behavior (Salam et al., 2009). The discrepancy in our findings for PE and PE+EE effects at a young age could be attributed to differences in cages, running wheel dimensions, or handling time during the experiment and may require further validation. However, there also is evidence to suggest that physical exercise is anxiogenic in young male C57BL/6 mice (Fuss et al., 2010) and has no effects on anxiety-like behavior in young male Sprague–Dawley rats (Kim et al., 2004). In the former study, however, mice were single-housed, and social isolation is an established etiology for the onset of psychological stress and subsequent anxiety-like behavior in rodents (Hart et al., 2010).

Short-term PE induces depressive-like behavior at an early age

EE mice showed the least depressive-like behavior among all groups at 4, 9, and 14 months of age, although results were not significant when compared to controls. EE has been shown to attenuate depressive-like behavior associated with prenatal chronic stress (Yang et al., 2006) but not with early life stress (Cui et al., 2006). We house mice in controlled environmental conditions where chances of stress were minimal. Hence, it is possible that the non-significant results for EE and PE+EE in FST could be due to the short duration of the treatment, and hence longer treatment with EE could be more effective in the treatment of depressive-like behavior. Like EE, PE+EE mice were not significantly different to control mice for depressive-like behavior at all ages. PE, however, acted as an inducer of depressive-like behavior at 4 months of age, as evidenced by significantly higher immobility time of 4-month PE mice compared with both 4-month-old controls and EE mice in the FST. However, research has shown that moderate PE alleviates depressive-like behavior in rodents (Duman et al., 2008; Patki et al., 2014; Carek et al., 2011). Indeed, our finding contradicts another similar research where voluntary wheel running of 3-4 weeks produced an antidepressant-like behavioral response in the FST in 8- to 10-week-old wild-type mice (Duman et al., 2008). No sex differences were observed for depressive-like behavior. The longer immobility time in FST could represent a faster adaptation of young PE mice to an acute stressor, thereby switching from active to passive behavior to conserve energy for better survival (Molendijk & de Kloet, 2015). However, in the absence of similar adaptation by PE mice at middle and late-middle age, the longer immobility time of young PE mice in FST during our experiments is indicative of depressive-like behavior and not coping strategy, since we find no reason for PE mice losing their capability to adapt to an acute stressor after young age. The discrepancy in our findings for PE in FST from past studies could be due to differences in cages, running wheel dimensions, or handling time during the experiment and may require further validation.

Short-term PE adversely affects spatial learning of female mice at early and late-middle age. Short-term EE improves retention of spatial learning at late-middle age

The observed impaired spatial learning in 4- and 14-month-old PE mice was primarily driven by 4- and 14-month female PE, showing significantly higher latency to escape than age-matched female controls. At 14 months, PE mice also showed significantly lower spatial learning than EE and PE+EE mice, with contribution from both PE sexes but more significantly from female PE mice. The analysis of entry errors also shows a trend for male and female PE mice low on spatial learning. Overall, this suggests that PE at a young age and late-middle age may affect spatial learning of female C57BL/6 mice adversely.

Furthermore, analyzing the escape latency of controls in the probe trial of the Barnes maze revealed that while retention of spatial learning improved from 4 to 9 month, it deteriorated after that. Interestingly, EE reversed the age-related decline in spatial memory at 14 months, but without any significant sex differences. EE mice, however, also showed a significantly higher number of entry errors in locating the escape box at 14 months compared with 4- and 9-month EE mice. This low escape latency with a high number of entry errors could indicate higher exploration activity and enhanced cognition in searching the escape box with visual input. These findings on EE that we report comply with the previously published literature. For example, short-term EE without running wheels has been shown to improve spatial memory of the young male Wistar and Sprague-Dawley rats (Falkenberg et al., 1992; Birch et al., 2013). Similarly, short-term EE with running wheel has been shown to improve spatial memory by promoting hippocampal neurogenesis, dendritic arborization, and spine density in the brain of Wistar rats (Leggio et al., 2005; Segovia et al., 2006). There is also evidence suggesting that EE from weaning to 4 months of age in rats accelerate the acquisition of spatial information in naïve male Wistar rats (Van Waas & Soffié, 1996). Because spatial memory is the function of the dorsal hippocampus (Sampedro-Piquero et al., 2013; Lee & Kesner, 2004), it seems that there is a substantial improvement in the integrity of dorsal hippocampus after EE.

Unlike EE, PE showed detrimental effects on the spatial learning and retention of spatial memory on a Barnes maze test. However, PE has widely been accepted as a modulator of brain functions and shown to improve spatial memory in young, middle age and old age male rodents (Albeck et al., 2006; Anderson et al., 2000; Radák et al., 2001; Nichol et al., 2009; Singhal et al., 2014; Van Praag et al., 1999; Cassilhas et al., 2012; Uysal et al., 2005). Notably, all of these studies except for two utilized only male animals for research. Among the two, one balanced research for gender (Nichol et al., 2009), and the other used only female animals (Anderson et al., 2000) but both reported improvement in the spatial memory after PE. Lack of research on female mice in past studies may be a reason for the discrepancy that we report for spatial learning and memory after PE in this study. The species and strain of rodents also could have a significant influence on the rate of spatial memory acquisition. Moreover, PE in short-term may have served as a stressor at an early age for female C57BL/6 mice. At late-middle age, the adverse effects of PE were seen on both male and female mice, which could be attributed to the altered physiological response of the body to acute physical strain after middle age. This hypothesis, however, needs validation through future research. Other factors, such as the design of our study, cage dimensions, the presence of one wheel per cage, handling frequency, and mouse-mouse interactions, may have played a role as well.

PE, PE+EE elicit differential gene expression at 9 months, but EE has no effect on gene expression at all ages

There is evidence to suggest that the change in behavior and memory is a measure of hippocampal integrity (Sampedro-Piquero et al., 2013; Scoville & Milner, 1957; Schacter et al., 1996; Engin & Treit, 2007; Campbell & MacQueen, 2004). Both PE and EE have a considerable influence on the expression of hippocampal genes associated with immune functions, aminergic pathways, neurotrophin metabolism, HPA axis, transcription factors, and mitochondrial function, and this can manifest into known anti-inflammatory, cell-proliferative and neuroprotective effects of the two treatments (Falkenberg et al., 1992; Rasmuson et al., 1998; Olsson et al., 1995; Keyvani et al., 2004; Mlynarik et al., 2004; Hullinger et al., n.d.; Buschler & Manahan-Vaughan, 2017; Torasdotter et al., 1998; Dahlqvist et al., 1999; Zajac et al., 2010; Tong et al., 2001; Rampon et al., 2000; Kaliman et al., 2011). In this context, we analyzed the differential expression of hippocampal genes of interest in response to short-term PE, EE, and combined PE and EE. We found no significant differential expression after all treatments at 4 and 14 months of age; however, at 9 months, a set of genes expressed differentially after PE and PE+EE treatments, whereas EE showed no effect.

The absence of treatment effects on the expression of hippocampal genes during early and late-middle age was unexpected. However, this could be attributed to genes being more responsive to the external environment during middle age than young and late-middle age. In middle age mice, we observed 11, 3, 6, and 7 differentially expressed genes during PE vs. controls, PE+EE vs. controls, PE vs. EE, and PE+EE vs. EE comparisons, respectively. None of the genes differentially expressed during EE vs. C and PE vs. PE+EE comparisons (see Table 5 for the list of differentially expressed genes between the mice groups at 9 months of age). Surprisingly, the differential expression of hippocampal genes after PE at 9 months did not manifest into significant behavioral changes. Mice were tested on a behavioral battery immediately after 4 weeks of treatments. It is possible that the effects of differentially expressed genes on behavior could have appeared at a later time point. Moreover, PE of longer duration may be required for observable behavioral changes that could be correlated to differentially expressed genes. Also, the genes that differentially expressed after PE at 9 months elicit neuroprotective effects, which could have allowed 9-month old PE mice to adapt to behavior similar to the age-matched control mice. All of these hypotheses need to be tested for validity in future research.

We performed pathway analysis for the differentially expressed genes using GeneMANIA to understand the molecular effects PE might have had on mice. Genes that were differentially expressed during simple comparisons between groups were mapped in GeneMANIA (Figure 13) and molecular processes statistically significant at FDR < 0.01 were noted (Table 6).

Figure 13
figure 13

Gene interaction and function prediction network obtained using GeneMANIA for (A) PE vs. C, (B) PE+EE vs. C, (C) PE vs. EE, and (D) PE+EE vs. EE. The colors represent significantly modulated molecular pathways.

Table 6. Differentially expressed genes and the associated statistically significant key molecular processes as observed during group comparisons

Three of the genes, i.e., Creb1, Smad2, and Il1r1 were upregulated in all groups with PE when compared to groups without PE (Table 6), suggestive of the effects of PE and not EE on these genes. Creb1 encodes for the transcription factor CREB (cyclic adenosine monophosphate response element-binding protein), which is involved in the regulation of genes important for neuronal regulation and survival of proteins such as brain-derived neurotrophic factor. The expression of Creb1 has been reported to increase in the hippocampus, temporal cortex, and nucleus accumbens of animal models of major depressive disorder and related disorders after chronic treatment with antidepressant drugs (Zubenko et al., 2003; Young et al., 1998; Blendy, 2006). Smad2 mediates TGF-β signals controlling a wide variety of cellular processes, including cell differentiation, proliferation, apoptosis, developmental fate specification, and differentiation of hippocampal granule neurons (Lu et al., 2005; Wang et al., 2011; Bottner et al., 2000). Similarly, IL-1R1 signaling mediates innate and acquired immune responses (Dinarello, 1988), associated with protection or aggravation of the disease states in the CNS, infections, and cancers (Durrant et al., 2014; Wolf et al., 1994). Knocking out these genes have shown to result in apoptotic cell death and affects brain functions and neurobiology adversely (Wang et al., 2011; Castro-Jorge et al., 2017; Becher et al., 2002). The upregulation of these three genes in all groups after PE at nine months suggests the neuroprotective role of PE, which has also been reported earlier (Campos et al., 2016).

Two of the genes, i.e., Sirt1 and Il12a, were upregulated in PE mice compared to control and EE mice. These genes were not differentially expressed in PE+EE mice. This suggests that the addition of EE to PE reduces the expression of these two genes. Sirt1 gene is a mitochondrial gene that encodes for SIRT1 protein, which influences energy metabolism and is important for neuronal plasticity, cognitive functions, as well as protection against aging-associated neuronal degeneration and cognitive decline (Ng et al., 2015; Paraíso et al., 2013). PE has been shown to enhance Sirt1 activity (Kaliman et al., 2011). Il12a gene encodes cytokine IL 12a, which act as a growth factor for activated T and NK cells, enhance the lytic activity of NK cells, and stimulate the production of IFN-γ by resting peripheral blood mononuclear cells. Upregulation of these two genes after PE is again suggestive of the neuroprotective role of PE at nine months of age.

Other genes that were upregulated in response to PE or PE+EE when compared to control or EE mice are Foxo3, Cs, Actb, Prkaa1, Prkaa2, Htr1a, and Htr2a. Foxo3 activates genes that preserve cellular quiescence, prevents premature differentiation, and controls oxygen metabolism during aging (Salih & Brunet, 2008). Foxo3 maintains the neural stem cells (NSC) pool by preventing the premature proliferation of progenitors and preserving the relative quiescence of NSC. Indeed Foxo3 KO mice were found to have fewer NSC in vivo than wild-type counterparts (Renault et al., 2009). Interestingly, a study has shown that Foxo3 deficient mice are significantly less immobile in FST (Polter et al., 2009), suggesting that Foxo3 regulates the behavioral manifestation of depression. A review validated this further by reporting that Forkhead box O transcription factor influences the development of major depressive disorder by regulating neuronal morphology, synaptogenesis, and adult hippocampal neurogenesis (Wang et al., 2015). We, however, observed no significant change in depressive-like behavior after PE at middle age. This is because none of the groups, including controls, showed the development of depressive-like behavior at middle age during our experiments.

Mitochondrial genes Cs, Actb, Prkaa1, and Prkaa2 were also upregulated in the hippocampus of 9-month-old PE and PE+EE mice when compared to control and EE mice. A study has indeed shown that PE stimulates the expression of mitochondrial genes in the hippocampus (Tong et al., 2001), thus suggesting that there can be a possible link between mitochondrial function and behavioral changes. Moreover, during aging, there is an increase in dysfunctional brain mitochondria, increased oxidation products, and diminished mitochondrial functional activity, resulting in long-term behavioral changes. Moderate exercise has been reported to reverse mitochondrial dysfunction and lower content of oxidation products and improve related neurological performances (Navarro & Boveris, 2007). Cs encode citrate synthase, an enzyme of the tricarboxylic acid cycle that plays a key role in regulating energy generation of mitochondrial respiration. Actb encodes for protein beta (β)-actin that forms the structural framework inside developing neurons. PE has been previously reported to upregulate the Cs (Huang et al.,2006) and Actb expression (Cappelli et al., 2008). Prkaa1 and Prkaa2 genes encode proteins that are the catalytic subunit of the AMP-activated protein kinase (AMPK). AMPK maintains energy homeostasis within cells. Enhanced expression of Prkaa1 and Prkaa 2 gene has been reported after short-term PE in mice (King-Himmelreich et al., 2016). Together, the upregulation of these mitochondrial genes after PE suggests that short-term PE play a key role in cellular energy metabolism and neurogenesis during middle age.

The upregulated Htr1a and Htr2a genes encode serotonin 1A and 1B receptors which mediate inhibitory neurotransmission. Inactivation of Htr1a and Htr2a genes in mice results in behavior consistent with an increased anxiety and stress response. For example, allelic variation in 5-HT1A and 5-HT2A receptor expression has been shown to be associated with anxiety- and depression-related personality traits (López-Figueroa et al., 2004). Studies have shown that PE alters levels of serotonin in the brain of rodents (Brown et al., 1979). Upregulation of Htr1a and Htr2a genes, therefore, suggests improved serotonin neurotransmission in the hippocampus after short-term PE at middle age.

The downregulated immune genes Ifng and Tnfrsf1a encode IFN-γ and members of the TNF receptor family of membrane-bound and soluble proteins that interact with the ligand TNF-α. Both IFN-γ and TNF-α regulate complex injury promoting proinflammatory and apoptotic and protective cell proliferation and survival roles in the CNS depending on dose, disease phase, and cell development stage (Ottum et al., 2015; Álvarez et al., 2011). Moreover, evidence suggests that IFN-γ is involved in modulation of CNS plasticity and the activity of related brain monoamine systems, thereby also regulates memory function and anxiety and depressive states (Litteljohn et al., 2014; Campos et al., 2014). Similarly, TNF-α has been shown to play a role in the pathophysiology of affective-like behavior (Liu et al., 2012; Haji et al., 2012). The downregulation of Ifng and Tnfrsf1a again suggests that PE at nine months is neuroprotective.

Given the effects of PE on gene expression, as mentioned above, the observed changes in gene expression appear primarily due to the effects of PE and not EE during PE+EE vs. C comparison. It is, however, interesting to find other genes, for example, Prkaa2, Ifng, and Htr1a that were differentially expressed after PE lost significance after PE+EE. Also, finding the same genes being upregulated in 9-month PE mice when compared to age-matched controls and EE mice suggests that the EE might have no or minimal effects on the hippocampal gene expression.

A short note on the effects of short-term EE on the hippocampal gene expression

EE did not alter the expression of any GOI. This contrasts with previously published literature. For example, a study reported that EE upregulates the Htr5B gene expression and reverses the depression-like state caused by early life seizure experiences (Koh et al., 2007). Likewise, EE in the absence of exercise showed anxiolytic effects in 5-HT1AR KO mice (Rogers et al., 2016). EE has also been shown to modulate the expression of 5-HT receptors in the hippocampus of R6/1 transgenic mice enriched for four weeks, from 8 to 12 weeks of age (Pang et al., 2008). However, the authors also noted gender effects, and the Htr2a gene differentially expressed only in females. Hence, the observed change in the behavior despite no GOI being expressed differentially after EE could be attributed to gender effects that we could not measure due to low sample size for gene expression analysis. While our behavioral results suggest that EE elicits age-dependent effects on anxiety-like behavior and could be anti-depressive in long-term, longer exposure of EE may be required to elicit significant changes in hippocampal gene expression. We also propose that study confounders such as the age of mice, and differences in the delivery of enrichment such as box dimensions, number/type of toys and how frequently they were changed may have a role to play for the reported change in behavior. Since aging is a known factor for the differential expression of genes in the hippocampus, it is possible that the effects of short-term EE on hippocampal gene expression could be more significant in mice older than late-middle age. Also, we analyzed 43 genes only, which leaves scope for future studies to investigate if other genes are responsible for the change in behavior after EE that we have reported in this study.

Past studies that have reported gene modulation effects of EE also used running wheels to enrich the environment of rodents (Zajac et al., 2010; Braun et al., 2011). Moreover, the absence of EE effects on hippocampal gene expression could be due to factors such as differences in the delivery of enrichment such as box dimensions, number/type of toys (i.e., quantity, quality and complexity of EE) and how frequently they were changed (i.e., the novelty of EE). A study indeed reported that EE, in the absence of running wheels, has limited transcriptomic effects genome-wide (Hart et al., 2010). More research with long-term EE using a wide variety of objects may be required to understand the effects of EE on hippocampal gene expression fully. Also, if the inclusion of running wheels affects gene expression results in EE studies needs to be investigated further, especially when the presence of running wheel alone has been shown to modulate hippocampal gene expression in C57BL/6 mice (Petit-Demouliere et al., 2005; Porsolt et al., 1977).

Limitations of the study

Several study confounders such as the age of mice, controlled environmental conditions, design of our study, cage dimensions, handling frequency, mouse-mouse interactions, number of wheels per cage, differences in the delivery of enrichment such as box dimensions, number/type of toys and how frequently they were changed may account for findings in our study. Indeed, the paradigms of PE and EE varied greatly across studies in the past. For example, while some researchers provided continuous PE and EE (Segovia et al., 2006; Marashi et al., 2004), others provided these treatments intermittently to the rodents (Huang et al., 2017; Samorajski et al., 1987). Cage dimensions also varied from small one layered cages (Nakamura et al., 1999) to large multilayered cages (Marashi et al., 2004), and some studies used one (Berardi et al., 2007) and others two wheels per cage (Segovia et al., 2006). Few researchers handled rodents every day (Albeck et al., 2006) while others handled rodents once a week (Leggio et al., 2005). The presence of a dominant/aggressive mouse can also affect the behavior of the subordinate mouse/mice in a cage (Singhal et al., 2014). Likewise, the presence or absence of running wheels and enrichment objects during home cage test can also affect the results on locomotion. Future studies investigating the longer duration of PE in older transgenic rodent strains relevant to specific disease models and EE with a wide variety of objects in the presence and absence of running wheels may be required to validate our findings. Also, the sex differences analysis for gene expression results could not be performed during our study due to the relatively smaller sample size of gene expression data. Hence, this remains to be a subject for future research.

Conclusions

Our research work is one of the few to investigate differential effects of PE, EE, and their combination on behavior, memory, and hippocampal gene expression over the lifespan. We observed that short-term EE alleviated age-related cognitive decline and anxiety-like behavior, but without altering hippocampal gene expression. Conversely, PE was found to be detrimental at a young age for both cognitive- and affective-like behaviors and showed no effects at middle and late middle age despite altering hippocampal gene expression. Hence, EE seems to provide a better opportunity to reduce anxiety and depression symptoms in the short-term, and PE might not benefit brain functions as much as short-term EE. Perhaps, PE of longer duration may be required to elicit observable changes in behavior and memory. We also found it unusual that changes in the expression of genes after PE were not manifested in observable behavioral alterations. We suggest further research with a longer duration of EE and PE and their combination to validate our findings. The gene expression results are in general agreement with studies suggesting that PE is a major neurogenic and neurotrophic stimulus (Bechara & Kelly, 2013), and modulator of gene expression (Rampon et al., 2000). This, in turn, also suggests that EE in the absence of exercise improves behavior and enhances memory probably via mechanisms that are independent of BDNF upregulation and neurogenesis in the dentate gyrus. It has been reported that PE and EE enhance hippocampal neurogenesis via dissociable pathways: while EE increases the likelihood of survival of new cells, PE increases the level of proliferation of progenitor cells (Olson et al., 2006). Such a protective effect of EE on new cells may explain the beneficial effects of EE that we have reported in this work.