3 weeks of rest from CMS leads to a complete recovery of the behavioral phenotype
As we previously observed [5], 2 weeks of CMS induced the development of vulnerable and resilient phenotypes, as indicated by the significant reduction of sucrose intake after 1 (-56%, p < 0.001 vs No stress) and 2 (-68%, p < 0.001 vs No stress) weeks of CMS only in a subgroup of stressed animals that were, therefore, named vulnerable (CMS-vul) (Fig. 1B).
Vulnerable animals took 3 weeks to fully recover from the pathological phenotype (Fig. 1C). Indeed, we observed that starting from the first week after the end of CMS exposure, animals showed a statistically significant increase in the sucrose intake (1 week: +60%, p < 0.001 vs CMS-vul/2weeks; 2 weeks: +57%, p < 0.001 vs CMS-vul/2weeks) and, interestingly, 3 weeks of recovery from stress were enough to completely recover from the anhedonic phenotype (+ 148%, p < 0.001 vs CMS-vul/2weeks). Accordingly, at the last time point, sucrose intake of the CMS-vul + recovery group was not statistically different from the baseline, whereas, after 1 and 2 weeks of recovery, sucrose intake was significantly lower in comparison to the vulnerable group at the beginning of the experiment (week1 recovery: -38%, p < 0.001 vs CMS-vul/baseline; week2 recovery: -40%, p < 0.001 vs CMS-vul/baseline).
The IEGs Arc and Cfos are upregulated by ARS independently from the behavioral phenotype.
To map the neuronal responsiveness of brain areas involved in the acute stress response, we measured the expression levels of the immediate early genes Arc and Cfos in vHip, dHip, Amy and Pfc.
We observed that, in the vHip of non-stressed animals, both Arc and Cfos were increased by ARS (Arc: +33%, p < 0.01 vs No stress/Naive; Cfos: +404%, p < 0.001 vs No stress/Naive). Of note, regardless of the response to CMS in term of hedonic phenotype and of the period of rest, both the IEGs were upregulated by ARS in CMS-vul (Arc: +48%%, p < 0.001 vs CMS-vul/Naive; Cfos: +121%, p < 0.001 vs CMS-vul/Naive), in CMS-res (Arc: +34%, p < 0.05 vs CMS-res/Naive; Cfos: +154%, p < 0.001 vs CMS-res/Naive) and in CMS-vul + recovery (Arc: +56%, p < 0.001 vs CMS-vul + recovery/Naive; Cfos: +92%, p < 0.001 vs CMS-vul + recovery /Naive) (Fig. 2A-B).
Similarly, we observed the same effect due to ARS in all the experimental conditions in the dHip, Amy and Pfc (Table 2A-B-C).
These results indicate that the exposure to 2 weeks of CMS did not alter the ability of the brain regions considered to react with a novel challenging condition.
Table 2
Analysis of Arc and Cfos mRNA levels in the dorsal hippocampus (dHip) (A), the amygdala (Amy) (B) and the prefrontal cortex (Pfc) (C) of chronically stressed rats exposed to 1 hour of acute restraint stress (ARS).
IEGs | No stress | CMS VUL | CMS RES | CMS VUL-REC |
Naive | ARS | Naive | ARS | Naive | ARS | Naive | ARS |
A dHip | Arc | 100±8 | 152±3 *** | 94±6 | 151±8 ### | 114±10 | 159±11 §§§ | 127±10 * # | 162±10 ^^ |
Cfos | 100±16 | 470±34 *** | 202±40 | 475±80 ### | 189±15 | 422±25 §§§ | 207±20 | 437±57 ^^^ |
B Amy | Arc | 100±9 | 226±11 *** | 102±9 | 236±18 ### | 131±11 | 234±17 §§§ | 126±16 | 220±17 ^^^ |
Cfos | 100±19 | 451±43 *** | 93±12 | 444±29 ### | 161±14 | 473±21 §§§ | 186±23 *# | 419±38 ^^^ |
C Pfc | Arc | 100±18 | 245±35 *** | 103±14 | 251±15 ### | 209±25 *** | 321±7 §§ | 143±16 | 259±26 ^^ |
Cfos | 100±26 | 161±24 | 111±13 | 201±17 ## | 186±25 ** # | 298±34 §§ | 222±22 *** ## | 265±23 |
The data are the mean ± SEM: * p < 0.05, ** p < 0.01, *** p < 0.001 vs No stress/Naive; # p < 0.05, ##p < 0.01, ### p < 0.001 vs CMS-vul/Naive; §§p < 0.01 §§§p < 0.001 vs CMS-res/Naive; ^^p < 0.01, ^^^p < 0.001 VS CMS-vul + recovery/Naive) (two-way ANOVA, Fisher’s PLSD).
Corticosterone plasma levels are affected by chronic and acute stress exposure.
Given the fundamental role of the HPA axis activity in stress response [19], we measured corticosterone plasma levels and, as previously demonstrated [5], we observed that CMS elevated levels of circulating CORT selectively in vulnerable animals (+ 88%, p < 0.05 vs No stress/Naïve), whereas this modulation was observed neither in resilient nor in vulnerable + recovery groups. Moreover, in line with the notion that acute stressor leads to an enhancement of the HPA axis activity, we found that ARS increased levels of CORT in non-stressed animals (+ 100%, p < 0.05 vs No stress/Naïve) as well as in resilient rats (+ 76%, p > 0.05 vs No stress/Naïve) even if the augmentation in the latter group did not reach the statistical significance probably due to the complexity of the experimental settings (Fig. 3).
The results obtained highlight that the negative effects of chronic stress led to a disruption of the negative feedback that regulates the activity of the HPA axis specifically in vulnerable animals; at the same time, these results indicate that the ability of the HPA axis to properly respond to novel challenging condition, as it occurs in healthy subjects, is maintained also in animals resilient to CMS.
The ability of the vHip to deal with a novel challenge in terms of ERGs is preserved selectively in resilient animals.
Seen the effects on CORT levels, we decided to measure whether differences in its release could influence the transcription of the ERGs, also known as glucocorticoid responsive genes, namely Gadd45b, Sgk1, Dusp1 and Nr4a1, in a brain region-specific manner.
In vHip, we found that Gadd45b mRNA levels (Fig. 4A) were enhanced by ARS in non-stressed (+ 58%, p < 0.001 vs No stress/Naïve), in CMS-res (+ 20%, p < 0.05 vs CMS-res/Naïve), in CMS-vul + recovery (+ 40%, p < 0.001 vs CMS-vul + recovery /Naïve) groups but not in CMS-vul animals. Similarly, ARS induced a similar modulation of Sgk1 and Dusp1 (Fig. 4B-C respectively), with their mRNA levels being upregulated in unstressed (Sgk1: +84%, p < 0.001 vs No stress/Naïve; Dusp1: +61%, p < 0.001 vs No stress/Naïve), in resilient (Sgk1: +40%, p < 0.05 vs CMS-res/Naïve; Dusp1: +73%, p < 0.001 vs CMS-res/Naïve) as well as in recovery group (Sgk1: +42%, p < 0.05 vs CMS-vul + recovery/Naïve; Dusp1: +36%, p < 0.001 vs CMS-vul + recovery/Naïve), whereas this effect was completed blunted in vulnerable animals.
Differently, as shown in Fig. 4D, we observed that Nr4a1 gene expression was increased by ARS regardless of the behavioral phenotype (No stress: +45%, p < 0.001 vs No stress/Naïve; CMS-vul: +38%, p < 0.001 vs CMS-vul/Naïve; +30%, p < 0.001 vs CMS-res/Naïve; +40%, p < 0.001 vs CMS-vul/Naïve).
By contrast, we observed an overall upregulation of the ERGs expression due to the acute challenge in all the experimental settings in dHip, Amy and Pfc, as shown in Table 3A-B-C. Of note, in Pfc we found that ERGs were increased under basal conditions in resilient animals in comparison to non-stressed group (Table 3C).
These results suggest that the vHip seems to be the brain region mostly affected by CMS in the ability to modulate neuronal plasticity following a novel stressor.
Table 3
Analysis of Gadd45b, Sgk1, Dusp1 and Nr4a1 mRNA levels in the dorsal hippocampus (dHip) (A), the amygdala (Amy) (B) and the prefrontal cortex (Pfc) (C) of chronically stressed rats exposed to 1 hour of acute restraint stress (ARS).
ERGs | No stress | CMS VUL | CMS RES | CMS VUL-REC |
Naive | ARS | Naive | ARS | Naive | ARS | Naive | ARS |
A dHip | Gadd45b | 100±4 | 144±3 *** | 109±7 | 134±5 ## | 123±6 ** | 141±5 § | 104±6 § | 144±8 ^^^ |
Sgk1 | 100±2 | 162±12 *** | 97±5 | 129±11 # | 93±3 | 119±16 | 93±6 | 133±11 ^^ |
Dusp1 | 100±5 | 158±10 *** | 127±8 * | 145±9 | 118±3 | 151±7 §§ | 94±8 ## § | 136±5 ^^^ |
Nr4a1 | 100±5 | 135±6 ** | 111±6 | 133±7 # | 103±7 | 144±8 §§§ | 116±7 | 144±8 ^^ |
B Amy | Gadd45b | 100±7 | 170±8 *** | 113±10 | 166±5 ### | 128±7 * | 165±9 §§ | 128±7 * | 158±10 ^ |
Sgk1 | 100±6 | 152±12 *** | 103±6 | 124±8 # | 106±4 | 128±5 § | 97±6 | 112±5 |
Dusp1 | 100±5 | 168±9 *** | 92±6 | 166±10 ### | 108±7 | 150±9 §§§ | 114±9 | 136±7 |
Nr4a1 | 100±8 | 194±13 *** | 91±8 | 186±10 ### | 111±5 | 179±10 §§§ | 104±8 | 183±16 ^^^ |
C Pfc | Gadd45b | 100±10 | 200±22 *** | 96±14 | 235±23 ### | 144±10 *# | 257±26 §§§ | 139±12 | 226±15 ^^ |
Sgk1 | 100±4 | 126±11 ** | 92±3 | 123±6 ## | 111±7 | 104±7 | 103±4 | 112±9 |
Dusp1 | 100±7 | 149±14 * | 104±7 | 154±15 ### | 143±13 * | 150±13 | 148±5 * | 153±11 |
Nr4a1 | 100±13 | 271±20 *** | 141±7 | 263±28 ## | 185±27 * | 263±19 § | 197±13 * | 277±35 ^^ |
The data are the mean ± SEM: * p < 0.05, ** p < 0.01, *** p < 0.001 vs No stress/Naive; # p < 0.05, ## p < 0.01, ### p < 0.001 vs CMS-vul/Naive; § p < 0.05, §§ p < 0.01, §§§ p < 0.001 vs CMS-res/Naive; ^ p < 0.05, ^^ p < 0.01, ^^^ p < 0.001 VS CMS-vul + recovery/Naive) (two-way ANOVA, Fisher’s PLSD).
Z score activation of the ERGs highlights the role of vHip in the susceptibility to stress
We calculated the Z score activation in each brain region to get a combined overview of the modulation of the IEGs (Fig. 5A-B-C-D) and ERGs (Fig. 5E-F-G-H) following ARS.
We found a significant Z activation of the IEGs due to ARS, regardless of the behavioral phenotype, in all the brain regions examined (Fig. 5A-B-C-D). Moreover, in Pfc (Fig. 5D), we observed an effect of the stress exposure per sè, with an increased Z activation in resilient naïve animals in comparison to the non-stressed counterpart.
On the contrary, the z activation of the ERGs highlighted that 2 weeks of CMS affected the ability of the vHip to deal with ARS (Fig. 5E) selectively in vulnerable animals. By contrast, in dHip (Fig. 5F), Amy (Fig. 5G), and Pfc (Fig. 5H), Z activation indicated that, despite stress exposure, these brain areas were not impaired in their ability to mount a stress response.
Finally, as observed for the IEGs (Fig. 5D), the significant Z activation of the ERGs in the Pfc (Fig. 5H) of resilient animals suggests that the Pfc activates mechanisms of resilience to face the negative effects of the CMS procedure.
Z score activation of the ERGs in vHip significantly correlates with the anhedonic phenotype
Finally, we performed the Pearson correlation analysis to figure out the relationship between the behavioral phenotype and the Z activation of the ERGs in vHip, dHip, Amy, and Pfc. Interestingly, as revealed by Person product–moment correlation coefficient analysis, the sucrose intake positively correlated with Z activation of the ERGs specifically in vHip (R 2 = 0.368, p < 0.001) (Fig. 6A), but not in dHip (R2 = 0.056, p > 0.05) (Fig. 6B), Amy (R2 = 0.000, p > 0.05) (Fig. 6C) and in Pfc (R2 = 0.064, p > 0.05) (Fig. 6D), thus strengthening the role of the ventral subregion of the hippocampus in the susceptibility to develop stress-related disorders.