EE alleviated avoidance behaviors impaired by MD in male BALB/c mice
The experimental timeline is shown in Fig. 1a. The hippocampal regions of CA1–3 and DG sampled for microscopy were shown in Fig. 1b. The PAT and SOT schematic drawings are shown in Fig. 1c and 1d. In the PAT, because the mice prefer darkness, the mice entered the darkrooms immediately after the sliding doors were open; because the grids in the darkrooms were applied an electric current, the mice returned to the lit rooms after receiving the foot shocks. Thus, the longer latency to enter the darkroom represents the better avoidance. A two-way nonparametric ANOVA of the data at 0 h showed that EE increased the spending time in the lit room before entering the darkroom (F (1, 44) = 6.966, P = 0.011, Figs. 1e and f), suggesting that EE improved immediate avoidance. MD did not affect the spending time in the lit room (F (1, 44) = 0.052, P = 0.820, Fig. 1f). The results at 24 h showed that MD decreased, but EE increased the spending time in the lit room before entering the darkroom (MD, F (1, 44) = 7.676, P = 0.008; EE, F (1, 44) = 22.696, P < 0.001; Fig. 1g) and MD + EE mice spent less time in the lit room before entering the darkroom than ND + EE mice (P < 0.001), suggesting that MD impaired the mid-term avoidance memory and the effect could not be improved by EE in the PAT. A two-way nonparametric ANOVA of the data at 1 w showed that EE increased the spending time in the lit room before entering the darkroom (F (1, 44) = 6.892, P = 0.012, Fig. 1h), suggesting that EE improved long-term avoidance memory. A Kruskal-Wallis post hoc test showed that MD + EE mice spent less time in the lit room than ND + EE mice (P < 0.001, Fig. 1h), suggesting that EE improved the long-term avoidance memory, but the combination of MD and EE could not.
In the SOT, once the switch was open, the mice stepped on the platform immediately to avoid the electric shock from the grid with a 24V current. Thus, if ignoring the ability to jump, the longer latency to step on the platform represents the worse avoidance. A two-way nonparametric ANOVA of the data at 0 h showed that MD increased and EE decreased the latency to step on the platform (EE, F (1, 44) = 12.835, P = 0.002; MD, F (1, 44) = 15.488, P < 0.001; Fig. 1i), suggesting that EE improved immediate avoidance, but MD impaired. A Kruskal-Wallis post hoc test showed that MD + EE mice spent less time stepping on the platform than MD + NE mice (P < 0.001, Fig. 1i), suggesting that EE could restore the avoidance impairment of the SOT induced by MD. At 24 h, two-way nonparametric ANOVA showed that MD and EE had interaction effects (F (1, 44) = 9.072, P = 0.004, Fig. 1j), the post hoc test suggested that MD impaired the mid-term avoidance memory of NE mice (P = 0.003) and EE could alleviate that (P = 0.027, Fig. 1j). A two-way nonparametric ANOVA of the data at 1 w showed that MD increased the latency to step on the platform (F (1, 44) = 9.314, P = 0.004, Fig. 1k) and had interaction with EE (F (1, 44) = 5.186, P = 0.028). The Kruskal-Wallis test showed that MD + NE mice spent more time stepping on the platform than ND + NE mice (P = 0.001, Fig. 1k), suggesting that MD impaired long-term avoidance memory. The Kruskal-Wallis test also showed that MD + EE mice spent similar time stepping on the platform as ND + EE mice (P = 0.625, Fig. 1k), suggesting that EE could alleviate long-term avoidance memory impairment induced by MD. Together, these results of SOT showed that MD impaired the avoidance behaviors and EE improved.
EE increased cells in the CA1–3 of the dHIP and the interaction of MD and EE normalized in the CA1 but not CA3
DAPI is a fluorescent probe that forms a DAPI-nucleic acid complex by attaching DNA (Kapuscinski 1995). We performed DAPI staining to evaluate the nuclei containing DNA in the dHIP regions (Figs. 2a–d). Because most eukaryotes have only one nucleus, the area% of DAPI staining can indicate the cell number in the mouse dHIP. A two-way ANOVA of the data in the CA1 (Fig. 2a) showed that, probably, EE (F (1, 20) = 3.542, P = 0.074) increased nuclei in the CA1 and meanwhile had interaction with MD (F (1, 20) = 7.068, P = 0.015). A Tukey’s test showed that compared to ND + NE, the treatments of ND + EE (P = 0.021) and, probably, MD + NE (P = 0.083) increased cells in the CA1 (Fig. 2a). But, MD + NE mice had similar cells in the CA1 as MD + EE mice (P = 0.946). These results showed that EE increased cells in the CA1 but the interaction of MD and EE normalized the cells in the CA1 as the control mice. The two-way ANOVA of the DAPI staining in the CA2 showed that EE probably increased cells in the CA2 (F (1, 20) = 4.203, P = 0.053, Fig. 2b), MD did not affect cell numbers in the CA2, but the interaction of MD and EE normalized the cell number (F (1, 20) = 7.208, P = 0.014, Fig. 2b). The two-way ANOVA of the DAPI staining in the CA3 showed that EE increased cells in the CA3 (F (1, 20) = 5.776, P = 0.026, Fig. 2c), MD did not affect cell numbers in the CA3, and EE had not interaction with MD (F (1, 20) = 0.912, P = 0.352, Fig. 2c). The two-way ANOVA of the DAPI staining in the DG showed that neither MD nor EE affect cell numbers in the DG (MD, F (1, 20) = 0.031, P = 0.861; EE, F (1, 20) = 1.293, P = 0.269, Fig. 2c). Together, these results suggested that EE increased in the CA1–3 (Figs. 2e1–4), and the interaction of MD and EE normalized the cells in the CA1 but not CA3.
EE could restore the Nissl bodies increased by MD in the CA1 but not the CA3 and DG
The Nissl body constructs endoplasmic reticulum and fine granules (Palay and Palade 1955), representing the biosynthetic function of the neuron. We performed Nissl staining of Cresyl violet to evaluate the neuronal biosynthetic function in the regions of the dHIP (Figs. 3a–d). We calculated the integrated density/µm2 of Cresyl violet staining to compare the difference in Nissl bodies levels. The two-way ANOVA of the data in the CA1 revealed that EE (F (1, 20) = 5.049, P = 0.036, Fig. 3a) and its interaction with MD (F (1, 20) = 6.034, P = 0.023, Fig. 3a) decreased Nissl bodies. The Tukey’s test showed that MD + NE mice had higher Nissl bodies levels in the CA1 than ND + NE mice (P = 0.035, Fig. 3a) and MD + EE mice (P = 0.016, Fig. 3a), suggesting that EE could reduce the Nissl bodies increased by MD in the CA1. The two-way ANOVA of the data in the CA2 showed that neither MD nor EE affected Nissl bodies levels in the CA2 (Fig. 3b). The two-way ANOVA of the data in the CA3 showed that EE (F (1, 20) = 4.640, P = 0.043, Fig. 3c) and, probably, MD (F (1, 18) = 3.602, P = 0.072, Fig. 3c) increased Nissl bodies levels in the CA3. The Tukey’s test showed that ND + EE mice had higher Nissl bodies levels in the CA2 than ND + NE mice (P = 0.01, Fig. 3c) and MD + EE mice (P = 0.015), but MD + EE mice had similar Nissl bodies levels as MD + NE mice (P = 0.958), suggesting that EE could not alter the Nissl bodies levels in the CA3 of the MD mice. The data of the DG showed that MD increased Nissl bodies levels (F (1, 20) = 6.960, P = 0.016, Fig. 3d), and the MD + EE mice had similar Nissl bodies levels as MD + NE mice (P = 0.346, Fig. 3d), suggesting that MD increased Nissl bodies in the DG and EE could not restore that. Together, these results showed that for Nissl bodies, MD probably increased in the CA1, CA3, and DG, EE increased in the CA3 and decreased in the CA1, and EE could restore the Nissl bodies altered by MD in the CA1 but not the CA3 and DG (Figs. 3e1–4). These results revealed that experiences such as MD and EE might affect social behaviors via altering the neuronal protein synthesis in the dHIP.
MD decreased synaptic connections in the CA2–3 and DG, and EE restored that in the CA3 but not CA2 and DG
To observe the synaptic connections in the dHIP, we performed immunofluorescence staining of SYP and PSD95. SYP is a presynaptic membrane marker (Wiedenmann and Franke 1985), and PSD95 is a marker of the postsynaptic membrane (Kim et al. 2007). We calculated the Pearson’s R values for SYP and PSD95 colocalization to determine the synaptic connection levels. A two-way ANOVA of the data in the CA1 showed that neither MD nor EE altered the synaptic connections (MD, F (1, 20) = 3.522., P = 0.075, Fig. 4a). In the CA2, whether EE or not, MD downregulated the colocalization of PSD95 and SYP (F (1, 19) = 9.757, P = 0.005, Fig. 4b), suggesting that EE could not restore the synaptic connection’s downregulation induced by MD in the CA2. A two-way nonparametric ANOVA of the data in the CA3 showed that EE increased but MD decreased the synaptic connections in the CA3 (EE, F (1, 20) = 16.364, P = 0.001; MD, F (1, 20) = 9.205, P = 0.007; Fig. 4c). The Kruskal-Wallis test showed that MD + EE mice had higher synaptic connection levels than MD + NE mice (P = 0.021, Fig. 4c), which had lower synaptic connection levels than ND + NE mice (P = 0.005, Fig. 4c), suggesting that EE could restore the synaptic connection downregulation induced by MD in the CA3. The two-way ANOVA of the data from the DG region showed that MD reduced synaptic connection levels (F (1, 19) = 9.248, P = 0.007, Fig. 4d), and EE did not affect (F (1, 19) = 0.005, P = 0.942, Fig. 4d). Tukey’s test showed that MD + EE mice probably had lower synaptic connection levels than ND + EE mice (P = 0.0501, Fig. 4d), suggesting that MD affected the effect of EE to increase synaptic connections in the DG. Together, the abovementioned results showed MD decreased synaptic connections in the CA2–3 and DG, EE increased synaptic connection levels in the CA3, and thus, EE restored synaptic connection downregulation induced by MD in the CA3 but not CA2 and DG (Figs. 4e1–4).
EE restored the dendritic branches alteration induced by MD in the CA1 but not CA3
We conducted Golgi-Cox staining to exhibit neuron morphology in the dHIP. Using Sholl analysis, we evaluated the dendritic branches number of the pyramidal neurons in the CA1, CA2, and CA3 and granule cells in the DG. The intersections numbers in Sholl analysis represented the dendritic branches complexity of the neurons. A two-way nonparametric ANOVA showed that MD probably increased but EE decreased dendritic branches of pyramidal neurons in the CA1 (EE, F (1, 68) = 44.095, P < 0.001; MD, F (1, 68) = 3.306, P = 0.073; Fig. 5a). The Kruskal-Wallis test showed that MD + NE mice had more dendritic branches than ND + NE mice (P = 0.008, Fig. 5a) and MD + EE mice (P = 0.001), suggesting that EE could restore the dendritic branches alteration induced by MD in the CA1. A two-way nonparametric ANOVA of the total intersections of CA2 pyramidal neurons showed that neither MD nor EE altered branches of pyramidal neurons in the CA2 (MD, F (1, 68) = 0.002, P = 0.969; EE, F (1, 68) = 1.653, P = 0.203; Fig. 5b). A two-way ANOVA of the total intersections of CA3 pyramidal neurons showed that EE increased but MD decreased the branches of the CA3 pyramidal neurons (EE, F (1, 68) = 8.712, P < 0.004; MD, F (1, 68) = 20.80, P < 0.001, Fig. 5c). The Tukey’s test showed that ND +EE mice had more dendritic branches than ND + NE mice (P = 0.007) and MD + EE mice (P < 0.001), suggesting that EE increased the branches of pyramidal neurons in the CA3 but could not restore the effect of MD on the branches of pyramidal neurons in the CA3. A two-way nonparametric ANOVA of total intersections of granular cells in the DG showed that neither MD nor EE altered the dendritic branches of granular cells in the CA3 (EE, F (1, 68) = 0.327, P = 0.569; MD, F (1, 68) = 2.473, P = 0.120, Fig. 5d). Together, these abovementioned results indicated that MD decreased dendritic branches in the CA3 and, probably, increased dendritic branches in the CA1, and oppositely, EE increased dendritic branches in the CA3 and decreased dendritic branches in the CA1. However, EE restored the dendritic branches alteration induced by MD in the CA1 but not CA3 (Figs. 5e1–4).
MD increased CRHR1 levels in the CA1 and decreased the CRHR1 levels in the CA2. The effects of EE on the CRHR1 levels depend on whether MD or not
Next, we checked CRHR1 and CRH levels in the dHIP regions using immunohistochemistry to find out the mechanism for the neuronal morphology alteration. We calculated the integrated density/µm2 of immunochemistry staining to evaluate the protein levels among the groups. The analysis of the CA1’s data showed that EE probably decreased CRH protein levels (F (1, 20) = 3.278, P = 0.085, Fig. 6a) and MD elevated CRHR1 protein levels in the CA1 (F (1, 20) = 5.387, P = 0.031, Fig. 6b). The post hoc tests showed that MD + NE mice had higher levels of CRH than ND + NE mice (P = 0.027, Fig. 6a) and lower levels of CRH than MD + EE mice (P = 0.0018, Fig. 6a), suggesting that EE could decrease the CRH levels increased by MD. However, EE could not affect the CRHR1 levels increased by MD in the CA1 (F (1, 20) = 1.327, P = 0.255, Fig. 6b). The analysis of the CA2’s data showed MD did not affect CRH levels (F (1, 20) = 1.721, P = 0.204, Fig. 6c) but decreased the CRHR1 levels (F (1, 18) = 9.502, P = 0.0064, Fig. 6d), and EE did not affect both levels of CRH and CRHR1, and meanwhile, the interaction effect of MD and EE was to decrease the levels of CRH (F (1, 20) = 25.31, P < 0.001, Fig. 6c) and CRHR1 (F (1, 20) = 15.996, P = 0.001, Fig. 6d) in the CA2. The post hoc tests showed that compared to the ND + NE mice, both MD + NE (P = 0.001) and ND + EE mice (P = 0.015, Fig. 6c) had higher CRH levels in the CA2. Besides CRH, ND + EE mice also had higher CRHR1 levels in the CA2 (P = 0.045, Fig. 6d). MD + EE mice had lower CRH levels than MD + NE mice (P = 0.006, Fig. 6c) and, probably, CRHR1 levels (P = 0.070, Fig. 6d) than ND + EE mice, suggesting that EE decreased the CRH levels increased by MD and MD could affect the effect of EE to increase the CRHR1 levels in the CA2. A two-way ANOVA of the CA3’s data showed the interaction effect of MD and EE was to decrease the levels of CRH (F (1, 20) = 13.83, P = 0.001, Fig. 6e) and CRHR1 (F (1, 20) = 40.74, P < 0.001, Fig. 6f) in the CA3. The Tukey’s test showed that MD + NE mice had higher CRHR1 levels and, probably, CRH levels than ND + NE mice (P = 0.001, Fig. 6f; P = 0.099, Fig. 6e) and MD + EE mice (P = 0.001, Fig. 6f; P = 0.042, Fig. 6e), suggesting that MD increased CRH-CRHR1 signaling of the CA3 and EE could restore. The Tukey’s test also showed that ND + EE mice had higher CRH levels than MD + EE mice (P = 0.050, Fig. 6e), and higher CRHR1 levels than ND + NE mice (P = 0.003, Fig. 6f) and MD + EE mice (P = 0.001, Fig. 6f), suggesting that MD could affect the effects of EE to increase CRH and CRHR1 levels in the CA3. A two-way ANOVA of the DG’s data showed that neither MD nor EE altered the CRH and CRHR1 levels, but they had the interaction effect on the levels of CRH (F (1, 20) = 7.593, P = 0.012, Fig. 6g) and CRHR1 levels (F (1, 20) = 6.533, P = 0.019, Fig. 6h). Tukey’s tests showed that MD + EE mice tended to have lower CRHR1 levels than MD + NE mice (P = 0.060, Fig. 6h), suggesting that EE was beneficial for decreasing CRHR1 levels in the MD mice. Taken together, for CRHR1 levels, MD increased in the CA1 and CA3 but decreased in the CA2. EE increased CRH levels in the CA1. For MD mice, the EE treatment decreased CRH levels of the CA1–3 but only decreased CRHR1 of the CA3 region.
To explore the effects of CRH-CRHR1 signaling to modulate the neuron morphology in the dHIP, we performed the nonparametric Spearman correlations between CRHR1 levels and the tested parameters of neuron morphology. Without considering the effects of the groups, the results showed that CRHR1 protein levels were positively correlated with CRH protein levels in all tested regions of the dHIP (CA1, n = 22 mice, r = 0.581, P = 0.005, Fig. s2a1; CA2, n = 22 mice, r = 0.429, P = 0.047, Fig. s2a2; CA3, n = 24 mice, r = 0.628, P = 0.001, Fig. s2a3; DG, n = 24 mice, r = 0.534, P = 0.007, Fig. s2a4). The CRHR1 levels were positively correlated with %area of DAPI staining in the CA1 (n = 23 mice, r = 0.433, P = 0.039, Fig. s2b1) and CA2 (n = 21 mice, r = 0.407, P = 0.067, Fig. s2b2) but not CA3 (Fig. s2b3) and DG (Fig. s2b4). The CRHR1 levels were positively correlated with integrated density/µm2 of cresyl violet-staining in the CA3 (n = 24 mice, r = 0.567, P = 0.004, Fig. s2c3) but not CA1 (Fig. s2c1), CA2 (Fig. s2c2), and DG (Fig. s2c4). The CRHR1 levels were correlated with the Pearson’s R of SYP and PSD95 colocalization positively in the CA1 (n = 23 mice, r = 0.426, P = 0.043, Fig. s2d1), CA2 (n = 20 mice, r = 0.458, P = 0.042, Fig. s2d2) , and DG (n = 20 mice, r = 0.515, P = 0.020, Fig. s2d4) , and negatively in the CA3 (n = 20 mice, r = −0.438, P = 0.053, Fig. s2d3). Together, the results suggested that CRHR1 levels were positively correlated with CRH protein levels in all regions of the dHIP. The results also suggested that CRHR1 levels were associated with synaptic connection levels in all regions of the dHIP (positively in the CA1, CA2, and DG, negatively in the CA3), Nissl body levels in the CA3, and cell nuclei levels in the CA1 and CA2.
MD elevated the OTR level in the CA1 region of the dHIP
Accumulating evidence shows that reciprocal regulation of the CRH and OT systems balances stress (Dabrowska et al. 2011). We checked the OTR level associated with stress, energy metabolism, and social attachment (Onaka et al. 2012). A two-way ANOVA of the data in the CA1 showed that MD and EE had interaction effects (F (1, 20) = 10.76, P = 0.004, Fig. 7a). The Tukey’s test showed that MD + NE mice had higher OTR levels in the CA1 than ND + NE mice (P = 0.008, Fig. 7a); the MD + EE mice had lower OTR levels than MD + NE mice (P = 0.041, Fig. 7a), suggesting the EE could decrease OTR levels increased by MD in the CA1. A two-way ANOVA of the data of the CA2 showed that neither MD nor EE altered the OTR levels in the CA2, but they probably interacted with each other (F (1, 20) = 3.443, P = 0.078, Fig. 7b). In the CA3 and DG, the two-way nonparametric ANOVA showed that neither MD nor EE affected the OTR levels, but MD and EE probably had interaction effects (CA3, F (1, 20) = 6.961, P = 0.016, Fig. 7c; DG, F (1, 20) = 3.137, P = 0.092, Fig. 7d). The above results showed EE restored the OTR levels increased by MD in the CA1. Both MD and EE did not affect the OTR levels of the CA2–3 and DG.
To ensure the OTR level in the dHIP, we checked the OTR level in the lateral entorhinal area (ENTL). Fig. s3a shows the relationship between the ENTL and HIP. Layers Ⅲ and Ⅴ of the ENTL are closely associated with the CA1 in the HIP, and Layer Ⅱ is linked to the CA3 and DG (Deng et al. 2010). Fig. s3b shows the representative images of OTR immunochemistry. A two-way ANOVA of the total OTR levels in the ENTL showed that neither MD nor EE affected the total OTR levels, but MD and EE had interaction effects (F (1, 20) = 9.302, P = 0.006, Fig. s3c). However, Tukey’s test did not show any differences. Next, we analyzed the OTR levels in Layer Ⅱ, Ⅲ–Ⅳ, and Ⅴ–Ⅵ of the ENTL (Fig. s3d–f). The results of Layer Ⅱ showed that neither MD nor EE probably affect OTR level (F (1, 20) = 3.180, P = 0.090, Fig. s3d). The two-way ANOVA of Layer Ⅲ–Ⅳ also showed that neither MD nor EE affect the OTR levels (Fig. s3e). The results of Layer Ⅴ–Ⅵ showed that MD probably increased the OTR levels (F (1, 19) = 3.951, P = 0.062, Fig. s3f), and meanwhile, MD and EE interacted with each other (F (1, 19) = 6.344, P = 0.021, Fig. s3f). The Tukey’s test showed that MD + NE mice had higher OTR levels in the Layer Ⅴ–Ⅵ than ND + NE mice (P = 0.027, Fig. s3f) and, probably, MD + EE mice (P = 0.066, Fig. s3f), suggesting that EE could decrease OTR levels in the Layer Ⅴ–Ⅵ increased by MD. Thus, the results reassured that MD elevated the OTR level in the CA1 region, which EE could alleviate.
Next, we performed the nonparametric Spearman correlation analysis to show the relationships between CRHR1 levels and OTR levels and the parameters of neuronal morphology. The result showed that the CRHR1 levels correlated with OTR levels in the CA1 (n = 20 mice, r = 0.483, P = 0.031, Fig. s4a1) and CA3 (n = 21 mice, r = 0.436, P = 0.048, Fig. s4a3), but not DG (Fig. s4a4) and CA2 (Fig. s4a2). The OTR levels were probably correlated with the %area of DAPI in the CA1 (n = 21 mice, r = 0.400, P = 0.072, Fig. s4b1), CA2 (n = 21 mice, r = 0.487, P = 0.025, Fig. s4b2), and DG (n = 20 mice, r = 0.550, P = 0.012, Fig. s4b4), but not CA3 (Fig. s4b3). The OTR levels were correlated with integrated density/µm2 of cresyl violet-staining positively in the CA1 (n = 22 mice, r = 0.654, P = 0.001, Fig. s4c1) and, probably, negatively in the CA2 (n = 21 mice, r = −0.395, P = 0.077, Fig. s4c2), but not CA3 (Fig. s4c3) and DG (Fig. s4c4). The OTR levels were correlated with the Pearson’s R of SYP and PSD95 colocalization positively in the CA2 (n = 20 mice, r = 0.462, P = 0.040, Fig. s4d2) and DG (n = 20 mice, r = 0.484, P = 0.031, Fig. s4d4), negatively in the CA3 (n = 21 mice, r = −0.579, P = 0.006, Fig. s4d3). In the CA1 of the ND + NE mice, the OTR levels and Pearson’s R of SYP and PSD95 colocalizations were negatively correlated (n = 5 mice, r = −0.975, P = 0.033, Fig. s4d1). Together, these results indicated that the OTR levels positively correlated with CRHR1 levels in the CA1 and CA3, cells number in the CA2, Nissl bodies in the CA1, and synaptic connections in CA2 and DG (Fig. 8).
Further, to explore the balances of CRHR1 and OTR levels, we compared the ratios of CRHR1 to OTR among groups in the tested dHIP regions. A two-way ANOVA of the CA1 showed that MD and EE had interaction effects of elevating the OTR levels of the CA1 (F(1, 20) = 6.142, P = 0.022, Fig. 7e1). However, Tukey’s test did not show any differences between groups. The two-way nonparametric ANOVAs showed that neither MD nor EE affect the ratios of CRHR1 to OTR in the CA2 and CA3 (Fig. 7e2 and Fig. 7e3). A two-way ANOVA of the DG’s data showed that EE decreased the ratios of CRHR1 to OTR (F (1, 20) = 4.87, P = 0.039, Fig. 7e4), but no difference was found in the following Tukey’s test, suggesting that EE disrupted the balance of the CRHR1 and OTR in the DG. Together, EE alone decreased the ratios of CRHR1 levels to OTR levels in the DG, MD and EE had interaction to elevate the ratios in the CA1.