Acute footshock stress induces anhedonic behavior in vulnerable animals that is preserved for up to two weeks
As previously reported [9], we found a significant reduction of sucrose intake in stressed animals (Mann Whitney test, p < 0.05), suggesting anhedonic behavior (Fig. 1B). However, looking at the distribution of single values, we observed that not all the rats displayed a reduction of sucrose intake compared to pre-stress condition. Thus, setting up an arbitrary cut-off of 75% of sucrose intake compared to baseline, we divided the animals into two groups: FS-vulnerable (FS-V), in which sucrose intake decreased by at least 25% from baseline, and FS-resilient (FS-R) (all the others) (Fig. 1C). Sucrose intake measured 24 h after FS remarkably decreased in FS-V rats compared to both control (CNT) and FS-R animals (Kruskal-Wallis p < 0.001; Dunn’s test: FS-V vs CNT p < 0.0001, FS-V vs FS-R p < 0.0001).
To assess whether different behavioral phenotypes were associated with changes in the stress response driven by the hypothalamic-pituitary-adrenal (HPA) axis at 24 h, we measured corticosterone serum levels (Fig. 1D), expression levels of mineralocorticoid (Fig. 1E) and glucocorticoid receptors (Fig. 1F), and phosphorylation of glucocorticoid receptors (Fig. 1G) in nuclear fractions of PFC from CNT, FS-R and FS-V rats sacrificed 24 h after stress. We found no significant changes among the three experimental groups (one-way ANOVA: MR F(2,26) = 0.4706, GR F(2,26) = 1.656, pGR F(2,26) = 1.679).
We also investigated how long after stress exposure FS-R and FS-V animals maintained a difference in behavioral phenotype. We measured sucrose intake at different time-points after stress exposure in a dedicated set of animals (Fig. 2A). Taking all the stressed rats as a single group, mixed-effects model revealed only a significant effect of time (F(4,85) = 4.228), but not of stress (F(1,23) = 1.770) or time x stress interaction (F(4,85) = 1.643) (Fig. 2B). However, when FS-R and FS-V animals were segregated 24 h after stress as above, we observed that FS-V rats showed a trend for reduced sucrose intake compared to both CNT and FS-R animals already 6 h after stress while the reduction was significant 24 h and 72 h after FS (significant effect of time F(4,81) = 6.762 and stress F(2,22) = 10.54; Tukey’s test: 6 h: FS-V vs CNT p = 0.0592, FS-V vs FS-R p = 0.0727; 24 h: FS-V vs CNT p < 0.01, FS-V vs FS-R p < 0.01; 72 h: FS-V vs CNT p < 0.05, FS-V vs FS-R p < 0.05) (Fig. 2C). Moreover, FS-V animals still displayed a reduction of sucrose intake compared to FS-R animals 1 and 2 weeks after stress (Tukey’s test, 1 week: FS-V vs FS-R p < 0.05; 2 weeks: FS-V vs FS-R p < 0.05).
Acute footshock stress induces selective changes of presynaptic glutamate release in the prefrontal cortex of FS-R and FS-V rats
We have previously demonstrated that acute FS induces a rapid enhancement of depolarization-dependent (but not basal) presynaptic glutamate release in the PFC of rats, which lasts for up to 24 h after stress [3, 5, 6, 9]. Here we asked whether resilience or vulnerability to the effects of acute stress is associated with alterations of glutamate release 24 h after stress. We selected this time point for further experiments for two main reasons: 1) we were interested in early determinants of acute stress resilience/vulnerability; 2) 24 h was the timepoint when FS-V rats were in the larger number (58% of stressed rats). The spontaneous glutamate release, here referred as basal release, selectively increased in the PFC of FS-V rats compared to controls (one-way ANOVA F(2,41) = 3.379 p < 0.05; Tukey’s test: FS-V vs. CNT p < 0.05) (Fig. 3A). Instead, 15 mM KCl depolarization-evoked glutamate release increased similarly in both FS-R and FS-V (one-way ANOVA F(2,31) = 9.882 p < 0.001; Tukey’s test: FS-R vs CNT p < 0.001, FS-V vs CNT p < 0.05) (Fig. 3B).
The measurement of synapsin I expression and phosphorylation at Ser9 in PFC synaptic membranes confirmed our previous evidence showing that Ser9 phosphorylation is a key mechanism for increasing depolarization-evoked glutamate release after acute FS exposure [5, 6]. Indeed, while total synapsin I expression was unchanged (one-way ANOVA F(2,33) = 0.1543 p > 0.05) (Fig. 3C), the levels of phospho-Ser9-synapsin I increased in both FS-R and FS-V at 24 h (one-way ANOVA F(2,32) = 5.354 p < 0.01; Tukey’s test: FS-R vs CNT p < 0.05, FS-V vs. CNT p < 0.05) (Fig. 3D).
Acute footshock stress induces changes in the expression of synaptic AMPA and NMDA receptors in the prefrontal cortex of FS-R and FS-V rats
To evaluate whether changes in presynaptic glutamate release were accompanied by postsynaptic alterations in AMPA and NMDA receptor levels and regulation, we measured protein expression and phosphorylation of main α-Ammino-3-idrossi-5-Metil-4-isossazol-Propionic Acid (AMPA) and N-Methyl-D-Aspartate (NMDA) receptor subunits in PFC from CNT, FS-R and FS-V animals 24 h after stress. We observed no differences in the total homogenate (Table 1). In synaptic membranes, we found no changes in GluA1 expression, phosphorylation at Ser845 and Ser831, GluA2 expression and phosphorylation at Ser880 (Table 2). On the other hand, while the expression levels of GluN1 (Fig. 3E) and GluN2B (Fig. 3G) were also unchanged (GluN1: one-way ANOVA F(2,18) = 1.695, p > 0.05; GluN2B: Kruskal-Wallis p > 0.05), GluN2A expression (Fig. 3F) and GluN2A/GluN2B ratio (Fig. 4H) were remarkably reduced in both FS-R and FS-V rats (GluNA2: one-way ANOVA F(2,16) = 9.008, p < 0.01; Tukey’s test: FS-R vs CNT p < 0.01, FS-V vs CNT p < 0.01; GluNA2/2B ratio: one-way ANOVA F(2,14) = 23.94, p < 0.0001; Tukey’s test: FS-R vs CNT p < 0.001, FS-V vs CNT p < 0.0001).
Table 1
Protein expression and phosphorylation of AMPA and NMDA receptor subunits in PFC homogenates.
| CNT | FS-R | FS-V | Statistics |
GluN1 | 100 ± 12.49 | 63.86 ± 10.32 | 71.78 ± 11.09 | ANOVA, F(2, 23) = 2.77, p = 0.084 |
GluN2A | 100 ± 11.65 | 66.09 ± 4.12 | 92.14 ± 13.95 | Kruskal-Wallis, p = 0.128 |
GluN2B | 100 ± 15.87 | 80.16 ± 8.81 | 87.80 ± 13.75 | ANOVA, F(2, 21) = 0.59, p = 0.564 |
GluN2A/GluN2B | 100 ± 9.46 | 77.68 ± 5.42 | 84.68 ± 8.80 | ANOVA, F(2, 22) = 1.98, p = 0.162 |
GluA1 | 100 ± 12.34 | 75.64 ± 2.94 | 88.98 ± 5.65 | ANOVA, F(2, 22) = 2.09, p = 0.148 |
pSer831-GluA1 | 100 ± 10.08 | 99.90 ± 14.30 | 92.01 ± 6.59 | ANOVA, F(2, 22) = 0.18, p = 0.838 |
pSer845-GluA1 | 100 ± 16.22 | 86.94 ± 6.93 | 70.42 ± 7.73 | Kruskal-Wallis, p = 0.262 |
GluA2 | 100 ± 7.33 | 86.81 ± 9.56 | 110.2 ± 9.58 | ANOVA, F(2, 24) = 1.69, p = 0.206 |
pSer880-GluA2 | 100 ± 10.84 | 91.16 ± 5.85 | 97.31 ± 7.32 | ANOVA, F(2, 24) = 0.55, p = 0.581 |
Data are reported as percentage compared to control (expression level of control sample is equal to 100) and as means ± SEM (N = 8–10). One-Way ANOVA or Kruskal-Wallis test. |
Table 2
Protein expression and phosphorylation of AMPA receptor subunits in PFC synaptic membranes.
| CNT | FS-R | FS-V | Statistics |
GluA1 | 100 ± 16.78 | 101.7 ± 18.75 | 75.25 ± 16.61 | ANOVA, F(2, 18) = 0.72, p = 0.498 |
pSer831-GluA1 | 100 ± 9.38 | 114.6 ± 20.29 | 90.78 ± 14.07 | ANOVA, F(2, 17) = 0.61, p = 0.556 |
pSer845-GluA1 | 100 ± 12.01 | 108.6 ± 18.70 | 92.98 ± 18.40 | ANOVA, F(2, 17) = 0.22, p = 0.806 |
GluA2 | 100 ± 16.06 | 54.79 ± 7.29 | 93.43 ± 10.27 | Kruskal-Wallis, p = 0.042 |
pSer880-GluA2 | 100 ± 13.67 | 98.47 ± 17.97 | 72.79 ± 15.75 | ANOVA, F(2, 17) = 0.87, p = 0.438 |
Data are reported as percentage compared to control (expression level of control sample is equal to 100) and as means ± SEM (N = 8–10). Statistics: One-Way ANOVA or Kruskal-Wallis test. |
Acute footshock stress induces changes in layers II/III dendrite morphology and neurotrophin levels in the prefrontal cortex of FS-R and FS-V rats
We have previously shown that acute FS induces a rapid retraction and simplification of prelimbic PFC layers II–III pyramidal neurons apical dendrites, which was measurable already after 24 h and sustained up to 14 days after stress [8, 9]. Here we found the total apical dendrite length (Fig. 4A) and number of dendritic branches (Fig. 4B) decreased in PFC layers II–III of both FS-R and FS-V rats (dendritic length: one-way ANOVA F(2,12) = 5.759, p < 0.05; Holm-Šídák test: FS-R vs CNT p < 0.05, FS-V vs CNT p < 0.05; N of branches: one-way ANOVA F(2,14) = 3.900, p < 0.05; Holm-Šídák test: FS-R vs CNT p < 0.05, FS-V vs CNT p < 0.05). However, Sholl analysis revealed a significant reduction in the number of intersections between 120 and 180 µm from soma only in FS-V animals compared to CNT (mixed-effects model, significant effect of distance F(14,189) = 20.58; Dunnett’s test: 120, 140, 160, 180 µm FS-V vs CNT p < 0.05) (Fig. 4C).
Local expression of neurotrophins has been linked to changes in dendrite development related to environmental factors, including exposure to stress [41–43]. Therefore, we measured the mRNA expression of the two main neurotrophins, Brain-Derived Neurotrophic Factor (BDNF) and Glial-Derived Neurotrophic Factor (GDNF). We found a significant reduction of total BDNF mRNA levels in the PFC of FS-V rats compared to CNT (Kruskal-Wallis p < 0.05; Dunn’s test: FS-V vs CNT p < 0.05) (Fig. 4D) and no changes of GDNF mRNA levels (Kruskal-Wallis p > 0.05) (Fig. 4E).
Acute footshock stress induces changes in glutamate release from glial perisynaptic processes in the prefrontal cortex of FS-R and FS-V rats
An increasing number of studies suggest that synapses require astrocyte signals for optimal performances [21]. Indeed, gliotransmitters released by astrocytes modulate presynaptic efficacy and postsynaptic responses [44, 45]. Although this neuron-astrocyte crosstalk mainly bases on glutamate signaling [46], possible changes in glutamate release from perisynaptic astrocyte processes after FS were never assessed before.
Here, for the first time we measured glutamate release from PFC glial perisynaptic processes (gliosomes) at different time points after acute FS. We did not found any difference in basal and 15 mM KCl depolarization-evoked glutamate release between control and FS animals immediately after stress (Student’s t-test; basal release: t = 0.5060, p > 0.05; depolarization-evoked release: t = 0.1749, p > 0.05) (Fig. 5A,B), as well as 6 h after stress (Student’s t-test; basal release: t = 0.9130, p > 0.05; depolarization evoked release: t = 0.5830, p > 0.05) (Fig. 5C,D) or 24 h after the start of stress (Student’s t-test; basal release: t = 0.07739, p > 0.05; depolarization evoked glutamate release: t = 1.501, p > 0.05) (Fig. 5E,F). However, when we separated rats into FS-R and FS-V, 24 h after stress, basal glutamate release from PFC gliosomes was unchanged (one-way ANOVA, F(2,37) = 0.5652, p > 0.05) (Fig. 5G) while depolarization-evoked glutamate release was significantly increased only in FS-V compared to both FS-R and control rats (one-way ANOVA, F(2,38) = 5.519, p < 0.05; Tukey’s test: FS-V vs CNT p < 0.05, FS-V vs FS-R p < 0.05) (Fig. 5H).
To understand the mechanisms underlying the increased glutamate release in PFC gliosomes of FS-V animals, we tested the involvement of calcium (Ca2+) and glutamate transporters. We found the basal glutamate release significantly increased in the absence of Ca2+. Still, it was strongly reduced by 10 µM DL-Tboa, a blocker of glutamate transporters [47], and unchanged in the presence of 10 µM KB-R7943, a blocker of the Na+/Ca2+ exchangers, [48] when working in a reverse mode, in both control and FS-V rats (two-way ANOVA, significant effect of treatment F(3,32) = 58.72, p < 0.0001 and of stress x treatment interaction F(3,32) = 3.750, p < 0.05; Tukey’s test: CNT: Ca2+-free vs physiological medium (PM) p < 0.01, DL-Tboa vs. PM p < 0.001; FS-V: Ca2+-free vs PM p < 0.001, DL-Tboa vs PM p < 0.01) (Fig. 5I). These experiments showed that spontaneous glutamate release from PFC gliosomes is mediated by glutamate transporters.
Notably, 10 µM DL-Tboa also abolished the depolarization-evoked glutamate release in both control and FS-V animals (in which glutamate release was confirmed to be higher compared to controls), suggesting a complete dependence on glutamate transporters (two-way ANOVA, significant effect of treatment F(3,32) = 27.22, p < 0.0001 and of stress x treatment interaction F(3,32) = 3.041, p < 0.05; Tukey’s test: FS-V vs PM p < 0.05; CNT: DL-Tboa vs PM p < 0.05; FS-V: DL-Tboa vs PM p < 0.0001). At the same time, no effects of Ca2+-free medium or KB-R7943 were measured (Fig. 5J).
Acute footshock stress induces changes in astrocyte proteins related to glutamate homeostasis in the prefrontal cortex of FS-R and FS-V rats
Protein expression levels of glutamate transporter 1 (GLT1), glutamine synthetase (GS), and cysteine-glutamate exchanger (xCt) were unchanged in PFC total homogenates (Table 3). In gliosomes from FS-R rats, GLT1 expression was increased compared to both CNT and FS-V rats (one-way ANOVA F(2,18) = 4.181 p < 0.05; Holm-Šídák test: FS-R vs CNT p < 0.05, FS-V vs FS-R p < 0.05) (Fig. 6A) and GS was reduced compared to CNT (Kruskal-Wallis p < 0.01; Dunn’s test: FS-R vs CNT p < 0.05) (Fig. 6B). xCt was increased compared to both CNT and FS-V animals (one-way ANOVA F(2,15) = 7.216 p < 0.01; Holm-Šídák test: FS-R vs CNT p < 0.05, FS-V vs FS-R p < 0.01) (Fig. 6C).
Table 3
– Expression of astrocytic proteins related to glutamate homeostasis in PFC homogenates.
| CNT | FS-R | FS-V | Statistics |
GLT1 | 100 ± 7.84 | 105.79 ± 12.14 | 84.38 ± 4.68 | ANOVA, F(2, 13) = 2.621, p = 0.110 |
GS | 100 ± 3.52 | 117.04 ± 9.34 | 98.28 ± 5.10 | Kruskal-Wallis, p = 0.0024 |
xCt | 100 ± 10.91 | 110.24 ± 6.52 | 113.52 ± 6.46 | ANOVA, F(2, 9) = 0.733, p = 0.506 |
Data are reported as percentage compared to control (expression level of control sample is equal to 100) and as means ± SEM (N = 4–6). One way ANOVA or Kruskal-Wallis test. |