Introduction

One of the core symptoms of major depressive disorder is anhedonia (diminished interest or pleasure).1 This suggests that dysfunction of brain reward pathways contributes to the pathophysiology of the disease. Dysfunction of the nucleus accumbens (NAc), an important relay station in the brain’s response to rewarding stimuli, has been consistently implicated in depression.2, 3, 4, 5 Chronic stress is known to induce anhedonia in animal models of depression,6 and its impact on the fine structure of the NAc has been described at the synaptic and molecular levels.7

The NAc, a part of the ventral striatum, is involved in motivation, reward, learning, problem-solving and motor function, all of which are affected in clinical depression.8 The NAc, a part of the ventral striatum, consists of two functionally distinct divisions: the core and the shell. Whereas the core is associated with goal-directed behavior, instrumental learning and motivation,9, 10 the shell is thought to integrate motivational valence and novelty.11 The majority of the NAc is composed of medium spiny γ-aminobutyric acid (GABA)ergic neurons, but the structure also contains smaller populations of interneurons that express GABA or acetylcholine.12 Importantly, the NAc is part of a complex network receiving glutamatergic, dopaminergic, serotonergic, noradrenergic and cholinergic afferents that originate in brain regions sensitive to stress-induced depression.13, 14, 15 Previous studies have reported that stress-induced dendritic and synaptic plasticity in the NAc are modulated by brain-derived neurotrophic factor (BDNF)16 and the neuropeptide melanocortin.17

In the present study, we searched for the structural and molecular correlates of anhedonic behavior in the NAc of rats that had been subjected to a chronic mild stress (CMS) protocol. In addition, we examined whether the stress-induced changes in plasticity within the NAc are influenced by two commonly used antidepressants, fluoxetine and imipramine. Our results show that stress-induced anhedonic behavior is associated with changes in the morphology of dendrites of medium spiny neurons of the NAc with concomitant changes in the levels of expression of various neurotrophins, cell adhesion molecules and synaptic proteins in the NAc. Antidepressant treatment with either fluoxetine or imipramine reversed these stress-induced morphological and neurochemical alterations while restoring hedonic behavior.

Materials and methods

Animals

Male Wistar rats (Charles-River Laboratories, Barcelona, Spain), weighing 300–400 g and aged 3 months, were used in this study. Animals were housed (three per cage) under standard laboratory conditions (12 h light: 12 h dark cycle, at 22 °C, relative humidity of 55%; free access to food and water). Animals were assigned to one of the two main treatment groups (control, n=14 and CMS, n=42). All procedures were carried out in accordance with the Portuguese national authority for animal experimentation, Direção Geral de Veterinária (ID: DGV9457) and in accordance with the guidelines for the care and handling of laboratory animals in the Directive 2010/63/EU of the European Parliament and of the Council.

Drugs

The drugs used were fluoxetine (10 mg kg−1, Kemprotec, Middlesbrough, UK) and imipramine (10 mg kg−1, Sigma-Aldrich, St Louis, MO, USA). Compounds were dissolved in 5% DMSO in 0.9% saline and administered intraperitoneally (1 ml kg−1) to animals, daily at 20:00 hours, during the last 2 weeks of the CMS protocol. The doses of antidepressant drugs administered daily were chosen based on their therapeutic effects as described in the literature.18

CMS

A slightly modified version of an unpredictable CMS protocol was used.6 It consisted of chronic exposure to unpredictable mild stressors (confinement to a restricted space for 1 h, placement in a tilted cage (30°) for 3 h, housing on damp bedding for 8 h, overnight illumination, food deprivation for 18 h followed by exposure to inaccessible food for 1 h, water deprivation for 18 h followed by exposure to an empty bottle for 1 h and reversed light/dark cycle for 48 h every 7 days) over 6 weeks. During the last 2 weeks of CMS, animals were given daily injections with saline (n=14), fluoxetine (n=14) or imipramine (n=14).

Open-field test

Locomotor activity was investigated using the open-field test in a room brightly illuminated by white light. Briefly, rats were placed in the center of an arena (43.2 × 43.2 cm2, transparent acrylic walls and white floor, MedAssociates, St Albans, VT, USA) and instant position was monitored online over a period of 5 min with the aid of two 16-beam infrared arrays. Total distances travelled and average speeds were used as measures of locomotor activity.

Sucrose-preference test

Anhedonia was assessed weekly during exposure to CMS using the sucrose preference test. Briefly, animals were allowed to habituate to the sucrose solution 1 week before the CMS protocol to establish baseline preference levels. To test sucrose preference, animals that were food- and water-deprived for 18 h were presented with two pre-weighed bottles containing 1% sucrose solution or tap water for a period of 1 h. Sucrose preference was calculated according to the formula: sucrose preference=(sucrose intake/(sucrose intake+water intake)) × 100, as previously described.19 Anhedonia was defined as a reduction in sucrose preference relative to baseline levels.

Forced swimming test (FST)

Learned helplessness was evaluated in the FST on the last day of exposure to CMS. Twenty four hours after a pre-test session (10 min), rats were placed in cylinders filled with water (25 °C; depth 30 cm) for a period of 5 min. Test sessions were assessed using a camera connected to a video tracking system (Viewpoint, Lyon, France); the system automatically calculated immobility time and latency to immobility. Learned helplessness behavior was defined as an increase in time of immobility and a decrease in latency to immobility.

Structural analysis

At the end of the experimental procedures, animals were killed under anesthesia. The total number of neurons in the core and shell divisions of the NAc was estimated using the optical fractionator20 in the right hemisphere of four animals from each experimental group. Brains were postfixed in 4% paraformaldehyde for 48 h and embedded in glycolmethacrylate. Volumes were determined by Cavalieri’s principle using the StereoInvestigator software (Microbrightfield, Willinston, ND, USA). For the tri-dimensional morphometric analysis, six animals from each treatment group were transcardially perfused with 0.9% saline and processed.21 Briefly, brains were immersed in Golgi–Cox solution22 for 21 days, transferred to a 30% sucrose solution and cut on a vibratome. Coronal sections (200 μm thick) were collected in 6% sucrose and blotted dry onto gelatin-coated microscope slides. They were subsequently alkalinized in 18.7% ammonia, developed in Dektol (Kodak, Rochester, NY, USA), fixed in Kodak Rapid Fix, dehydrated and xylene-cleared before coverslipping. Dendritic arborization and spine numbers and shape were analyzed in the core and shell divisions of the NAc. For each selected neuron, all branches of the dendritic tree were reconstructed at × 600 (oil) magnification using a motorized microscope (Axioplan 2, Carl Zeiss, Thornwood, NY, USA) and Neurolucida software (Microbrightfield). A 3D analysis of the reconstructed neurons was performed using the NeuroExplorer software (Microbrightfield). Forty neurons were studied for each animal and measurements from individual neurons from each animal were averaged. Several aspects of the dendritic morphology were examined: (i) total dendritic length was compared across the experimental groups; (ii) dendritic spine density (number of spines/dendritic length) was determined in branches that were either parallel or at acute angles to the coronal surface of the section; and (iii) spine morphology, specifically by classifying spines in the selected segments as mushroom-shaped, thin, wide and ramified spines, with the proportion of each type of spine being calculated for each neuron.23

Gene-expression measurements by quantitative real time PCR (qRT-PCR)

Levels of Bdnf, Ncam1 (neural cell adhesion molecule) and Syn1 (synapsin 1) mRNA were determined by qRT-PCR in NAc derived from four animals of each treatment group. Tissue from the NAc was dissected under the microscope. Briefly, 2 μg of total RNA were reverse transcribed using oligo-dt primers of the Superscript First-strand Synthesis system for RT-PCR (Invitrogen, Carlsbad, CA, USA). Hypoxanthine guanine phosphoribosyl transferase (Hprt) was used as internal standard for normalization. Oligonucleotide primers for Bdnf (sense 5′-GCGGCAGATAAAAAGACTGC-3′, antisense 5′-GCAGCCTTCCTTCGTGTAAC-3′), Ncam1 (sense 5′-AAAGGATGGGGAACCCATAG-3′, antisense 5′-TAGGTGATTTTGGGCTTTGC-3′), Syn1 (sense 5′-CACCGACTGGGCAAAATACT-3′, antisense 5′-TCCGAACTTCCATGTCC-3′) and Hprt (sense 5′-GCAGACTTTGCTTTCCTTGG-3′, anti-sense 5′-TCCACTTTCGCTGATGACAC-3′) were designed using the Primer3 software (Simgene, San Francisco, CA, USA) on the basis of the GenBank sequences NM_031521, NM_019133 and NM_012583, respectively. qRT-PCR was performed on a CFX 96TM real-time system instrument (Bio-Rad Laboratories, Hercules, CA, USA), with the QuantiTect SYBR Green RT-PCR reagent kit (Qiagen, Hamburg, Germany) according to the manufacturer’s instructions, using equal amounts of RNA from each of the samples. Product fluorescence was detected at the end of the elongation cycle. All melting curves exhibited a single sharp peak at the expected temperature.

Statistical analysis

After confirmation of homogeneity, appropriate statistical tests were applied to the data. Repeated-measures analysis of variance were used to analyze the results of the sucrose consumption test (SCT). Two-factor analysis of variance was used to evaluate further behavioral, structural and gene-expression data. Differences between groups were then determined by Tukey’s honestly significant difference test post-hoc analysis. Statistical significance was accepted for P<0.05. Results are expressed as mean±s.e.m.

Results

Behavioral results

Anhedonia was assessed weekly in the SCT, performed during the 6 weeks of the experimental protocol. Animals exposed to CMS revealed a significant decrease in sucrose preference when compared with control animals (F1,26=89.53, P<0.001) (Figure 1a), indicating an anhedonic phenotype. The administration of antidepressants in the last 2 weeks of the CMS protocol significantly reversed this behavioral phenotype (F2,52=4.875, P=0.011). As previously shown,24 the tricyclic antidepressant imipramine elicited an earlier onset of antidepressant action, reversing anhedonic behavior within the first week of treatment (P<0.001), whereas the response to the selective serotonin reuptake inhibitor fluoxetine was only evident after 2 weeks (P<0.001) (Figure 1a). Consistent with the SCT data, learned helplessness in the FST revealed that animals exposed to CMS showed increased immobility time (F1,26=37.568, P<0.001) (Figure 1b) and decreased latency to immobility (F1,26=32.894, P<0.001) (Figure 1c). Antidepressant treatment proved to be a significant factor in recovery from learned helplessness (F2,52=27.962, P<0.001); when compared with stressed rats that did not receive drug treatments, immobility time was significantly reduced by fluoxetine (P<0.001) and imipramine (P<0.001) and the latency time to immobility was significantly increased by both drugs (P<0.001). As previously described, the concordance between the anhedonic behavior in the SCT and the learned helplessness in the FST further strengthens the consistency of the CMS protocol as an animal model of depression.24 No significant effects of CMS or treatment with fluoxetine and imipramine were observed in locomotor activity evaluated in the open-field test (Supplementary Figure S1).

Figure 1
figure 1

Behavioral effects of chronic mild stress (CMS) on hedonic and learned helplessness behavior. (a) Sucrose preference in the sucrose-consumption tests performed during CMS. (b) Immobility time and (c) latency to immobility time in the forced swimming test performed after CMS and drug administration. *P<0.005; **P<0.001. Data presented as mean±s.e.m.

Structural analysis

Volumetric measurements of the core and shell divisions of the NAc revealed a tendency for increased volume in both the subregions in animals exposed to CMS without antidepressant treatment that, however, failed to reach statistical significance (shell: P=0.612; core: P=0.339) (Figures 2c and d). Stereological estimation of the total number of cells in the core and shell divisions of the NAc also did not reveal significant differences between experimental groups (shell: P=0.770; core: P=0.666) (Figures 2e and f).

Figure 2
figure 2

Stereological volumetric measurements and estimation of total number of neurons. (a) Paxinos atlas scheme, figure 56 (Bregma 1.8). (b) Low-power micrograph of a glycolmethacrylate-embedded brain section stained with Giemsa. Volumetric measurements of the core (c) and shell (d) divisions of the nucleus accumbens (NAc). Stereological estimation of the total number of cells in the core (e) and shell (f) divisions of the NAc. AC, anterior comissure; CMS, chronic mild stress; LV, lateral ventricle.

On the other hand, 3D morphometric analysis of Golgi-impregnated neurons revealed that exposure to CMS induced dendritic hypertrophy in medium spiny neurons, with significant increase in their total dendritic length in the core (F1,10=7.195, P=0.023) and shell (F1,10=8.148, P=0.017) divisions of the NAc (Figures 3g and h). Importantly, this hypertrophic effect of stress was reversed after administration of the antidepressants fluoxetine (core: P<0.001; shell: P=0.021) and imipramine (core: P=0.001; shell: P=0.041).

Figure 3
figure 3

3D morphometric analysis of Golgi-impregnated neurons using computer- assisted reconstructions of nucleus accumbens (NAc) medium spiny neurons. Representative neurons of different experimental groups in the core (a, c, e), and shell (b, d, f) divisions of the NAc. Cells are depicted in the XY orthogonal plan. Total dendritic length of neurons in the core (g) and shell (h) divisions of the NAc. Asterisk represents the comparison between the control and chronic mild stress (CMS) groups. Double asterisk represents the effects of antidepressant treatment. *P<0.03, **P<0.03. Data represented as mean±s.e.m.

Furthermore, exposure to CMS resulted in a significant increase in spine density in the core division of the NAc (F1,10=7.903, P=0.018) (Figure 4b). This gain in synaptic contacts was reversed by administration of fluoxetine (P=0.016) and imipramine (P=0.002). No significant changes were observed in spine densities in the shell division of the NAc (Figure 4c). None of the treatments led to significant differences in the relative abundance of mature (mushroom) and immature (thin) spines in the core and shell divisions of the NAC (Figures 4d and e).

Figure 4
figure 4

Spine density and morphological classification. Representative dendritic segments for spine analysis from each experimental group (a). Spine densities of nucleus accumbens (NAc) medium spiny neurons in the core (b) and shell (c) divisions of the NAc. Morphological classification of dendritic spines in the core (d) and shell (e) subregions of NAc. Asterisk represents the comparison between the control and chronic mild stress groups. Double asterisk represents the effects of antidepressant treatment. *P<0.02, **P<0.02. Data represented as mean±s.e.m.

Gene-expression studies

The mRNA expression of Bdnf, Ncam1 and Syn1 was measured by qRT-PCR. The expression levels of Bdnf were significantly increased in the NAc (F1,6=8.794, P=0.025) of animals exposed to CMS. These changes were reversed after the administration of fluoxetine (P=0.003) and imipramine (P=0.021) (Figure 5a). Additionally, the CMS group revealed significantly increased Ncam1 expression in the NAc (F1,6= 10.767, P=0.017) that was reversed by fluoxetine (P=0.001) and imipramine (P=0.027) (Figure 5b). Similarly, the expression levels of Syn1 were significantly increased in the NAc of rats exposed to CMS (F1,6=6.771, P=0.041) and reversed by administration of both fluoxetine (P=0.039) and imipramine (P=0.039) (Figure 5c).

Figure 5
figure 5

mRNA expression levels of (a) Bdnf (brain-derived growth factor), (b) Ncam1 (neural cell adhesion molecule 1) and (c) Syn1 (synapsin 1) in the nucleus accumbens measured by quantitative real-time PCR. All the genes analyzed showed an increased expression level in response to chronic mild stress (CMS), this effect being abrogated by the antidepressant treatment. Asterisk represents the comparison between control and CMS groups. Double asterisk represents the effects of antidepressant treatment. *P<0.05, **P<0.05. Data represented as mean±s.e.m. Hprt, hypoxanthine guanine phosphoribosyl transferase.

Discussion

The brain regions involved in the genesis of the psychopathological phenomena associated with depression are not completely identified. Many studies have suggested a role for the NAc in depression as anhedonia is a cardinal symptom in this condition, and the NAc is implicated in the regulation of motivation and reward processing.1, 2, 3, 4, 5 The results of the present study add important new support to this view. Specifically, they show that exposure to unpredictable CMS, an established mean to induce depression-like behavior in rodents, increases anhedonic behavior (reduced sucrose preference) in association with changes in the fine structure of the NAc. Although not affecting volumes or neuronal numbers in either the core or the shell of the NAc, chronic stress caused hypertrophy of the dendritic trees (increased dendritic lengths) and increased spine densities in the core division of the NAc. Importantly, the latter morphological adaptations were reversed by treatment with the tricyclic antidepressant imipramine as well as by the serotonin reuptake inhibitor fluoxetine; contemporaneously, these drugs attenuated the behavioral deficits seen in stress-exposed animals. Together, these observations add to the evidence that neuroplastic changes in the NAc contribute to the pathophysiology of depression and its pharmacologically induced recovery and point to the role of the NAc in the regulation of (an)hedonia.

The morphological findings in the NAc are interesting in that they contrast to the effects of chronic stress in the hippocampal formation and prefrontal cortex14, 25 that have been the focus of most earlier studies aimed at identifying the neuroanatomical and neurochemical substrates of depression; in both areas, chronic stress leads to significant dendritic atrophy.20 On the other hand, stress has been shown previously to induce dendritic hypertrophy in other corticolimbic structures implicated in the regulation of mood and affect, namely, the amygdala,26 bed nucleus of stria terminals,27 dorsolateral striatum and orbitofrontal cortex,28 and to increase dendritic spine density and functional excitatory synaptic plasticity in the NAc.7, 29 These observations reflect the complex neural networks that underlie the cluster of depressive symptoms.13

Stress-induced increases in dendritic arborization in the extended amygdala were previously ascribed to increases in BDNF expression.30 Importantly, the stress-induced increases in spine density and synaptic plasticity in the NAc are dependent on the activity of the inhibitor of KB kinase, a downstream target of BDNF.31 Here, we also observed an upregulation in BDNF mRNA levels that was concomitant with dendritic hypertrophy in the NAc of stressed rats; moreover, both events were reverted by the antidepressants fluoxetine and imipramine. Here it is important to note that, although the recent work by Malenka and colleagues dissociates NAc-mediated anhedonic behavior from behavioral despair,17 learned helplessness evaluated in the FST can be triggered by intra-VTA (ventral tegmental area) BDNF and blocked by overexpression of a dominant-negative mutant form of the BDNF receptor, TrkB.32 In addition, increased BDNF signalling in the VTA-NAc pathway has also been implicated in the display of other depression-like behaviors (decreased sexual behavior, decreased social interaction and a general decrease in locomotor activity), most of which are reversible by antidepressant treatment. Further, mice with selective knockdown of BDNF in the mesolimbic pathway fail to develop depression-like behavior after exposure to social defeat stress.16 Thus, it would appear that an imbalance in BDNF signalling in the NAc contributes to dysregulated mood by altering hedonic tonus. Interestingly, recent evidences suggest a similar modulatory effect of BDNF in addictive behavior, where the suppression of BDNF in the mesolimbic pathways seems to increase dopaminergic neuron excitability and promote reward.33

An interesting outcome of the present work is the confirmation that some of the same molecular pathways implicated in depression and antidepressant action in the hippocampus and in the prefrontal cortex also influence these phenomena at the level of the NAc, although very different effects are seen in these various regions. The stress-induced changes in behavior are not only accompanied by increased BDNF expression but also by alterations in the expression profiles of other genes implicated in neuroplasticity. We report here that chronic stress upregulates the levels of both NCAM and SYN1, changes that are commensurate with the observed increases in dendritic length and synaptic contacts of medium spiny neurons of the NAc. Given that this region is remarkably interconnected with the hippocampus and the prefrontal cortex, it remains to be demonstrated whether the changes herein described are causal or a mere consequence of the neurodegenerative effects triggered by stress in these ‘cortical’ regions. Interestingly, however, the orbitofrontal cortex, an area also highly connected with the NAc, displays hypertrophic changes that are concordant with the present findings, which may suggest that these changes may be circuit specific.28

In conclusion, the present results provide a clear demonstration of structural changes in the NAc, a central component of the mesolimbic pathway, in stressed animals displaying anhedonic behavior. Of notice, these structural changes are reversible after efficient antidepressant treatment and may relate with the altered expression of molecules involved in neuroplasticity. These observations pave the way for future studies addressing whether primary abnormalities in the reward circuitry are involved in the etiology of depression and whether antidepressants operate through direct actions in the mesolimbic pathway or indirectly through other brain regions.