Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Molecular Profiling of the Lateral Habenula in a Rat Model of Depression

  • Trine Christensen,

    Affiliation Translational Neuropsychiatry Unit, Department of Clinical Medicine, Aarhus University, Risskov, Denmark

  • Line Jensen,

    Affiliation Translational Neuropsychiatry Unit, Department of Clinical Medicine, Aarhus University, Risskov, Denmark

  • Elena V. Bouzinova ,

    elena_bouz@hotmail.com

    Affiliation Translational Neuropsychiatry Unit, Department of Clinical Medicine, Aarhus University, Risskov, Denmark

  • Ove Wiborg

    Affiliation Translational Neuropsychiatry Unit, Department of Clinical Medicine, Aarhus University, Risskov, Denmark

Abstract

Objective

This study systematically investigated the effect of chronic mild stress and response to antidepressant treatment in the lateral habenula at the whole genome level.

Methods

Rat whole genome expression chips (Affymetrix) were used to detect gene expression regulations in the lateral habenula of rats subjected to chronic mild stress (mild stressors exchanged twice a day for 8 weeks). Some rats received antidepressant treatment during fifth to eights week of CMS. The lateral habenula gene expression profile was studied through the gene ontology and signal pathway analyses using bioinformatics. Real-time quantitative polymerase chain reaction (RT-PCR) was used to verify the microarray results and determine the expression of the Fcrla, Eif3k, Sec3l1, Ubr5, Abca8a, Ankrd49, Cyp2j10, Frs3, Syn2, and Znf503 genes in the lateral habenula tissue.

Results

In particular we found that stress and antidepressant treatment affected intracellular cascades like growth factor receptor signaling, G-protein-coupled receptor signaling, and Wnt signaling – processes involved in the neuroplastic changes observed during the progression of depression and antidepressant treatment.

Conclusion

The present study suggests an important role of the lateral habenula in the development of depression-like conditions and correlates to previous studies demonstrating a significant role of the lateral habenula in depressive-like conditions and antidepressant treatment.

Introduction

Depression is a major cause of disability affecting about 350 million people worldwide. However, little is known about the underlying pathology of the disease. The complexity of the disease is massive involving several neurological mechanisms and affecting numerous regions of the brain. The therapeutic intervention of depression-associated symptoms has evolved rapidly over the past two decades. Nevertheless, currently available antidepressants have significant limitations, including slow onset of action and low rates of response including treatment refraction. Despite substantial advances in understanding depression symptomatology, as well as the antidepressant machinery, the underlying molecular mechanisms still remain unclear [1].

The lateral habenula (LHb) is a small brain structure mediating behavioral responses to pain, stress, anxiety, sleep and reward and its dysfunction is associated with depression, schizophrenia and drug-induced psychosis [2], [3]. Accumulating lines of evidences indicate that the LHb, which innervates multiple brain regions and directly influences the serotonergic, noradrenergic and dopaminergic brain systems, exhibits hyperactivity during depressed states. Increased neuronal activity of the habenula has been observed in animal models of stress [4] and depression [5], [6], and also in human depression [7], [8]. Conversely, lesions of the LHb were shown to induce antidepressant-like effects in rats exposed to inescapable foot shock (learned helplessness) [9] and chronic mild stress (CMS) [10]. Similarly, neuronal inhibition of the LHb by deep brain stimulation has shown ameliorating effects in the CMS model of depression [11] while pharmacological inhibition of the LHb mediates an antidepressant effect in an animal model of treatment resistant depression [12]. Indeed, one study reports antidepressant effects of deep brain stimulation in a treatment resistant depressed patient [13].

The CMS model of depression was developed to mirror anhedonia – one of the core symptoms of depression. Anhedonia is characterized as a lack of interest in otherwise pleasurable events or stimuli, and in the CMS model of depression this stress-induced decrease in reward sensitivity is measured as a reduction in sucrose consumption/preference [14]. The model fulfills all standard validity criteria and has a face validity that includes several symptoms of depression [15]. The CMS model of depression has been validated in several studies linking anhedonic-like behavior to deficits in cognitive behavior [16], neuroendocrine and immunological dysfunctions [17], [18], alterations at a cellular and structural level [19], [20], and dysfunctions in serotoninergic and GABAergic neurotransmitter systems [21], [22]. Additionally, several studies indicate specific gene- and protein expression profiles associating with anhedonic-like behavior [23][27]. The validity of the CMS model of depression is further strengthen by the fact that the model in addition to anhedonia also mirrors the aspect of stress-resilience, as a fraction (about 20%) of the stress exposed rats does not reduce their sucrose intake. Furthermore, antidepressant treatment mediates reversing effects of the stress-induced decrease in sensitivity to reward in about 50% of the treated rats, resulting in a drug responder group and drug non-responder group, respectively [23][25], [27]. Hence the CMS model of depression not only simulates the chronic and episodic aspect of clinical depression, but it also models the fact that conventional antidepressant therapy is inefficient in a substantial fraction of treated patients [28], [29].

In the present study we investigated global transcriptomic profiles of CMS phenotypes to search for biomarkers/molecular mechanisms underlying vulnerability to depression (anhedonia) and the long-term action of antidepressant treatment. The aim was to clarify the role of the LHb in depression and antidepressant treatment at a molecular level. The results of the present study indicate that stress-induced anhedonia and antidepressant treatment primarily associated with changes in intracellular signaling mechanisms in the LHb.

Materials and Methods

Ethics statement

All procedures involving animals were accepted by the Danish National Committee Ethics in Animal Experimentation (2008/561–477) and approved by the Aarhus University review board.

Animals

Male Wistar rats were purchased from Taconic, Denmark. The initial experimental group contained 329 rats. 249 rats were submitted to CMS while the rest (80 rats) were used as unchallenged control rats. After four weeks, rats exposed to CMS were divided according to their hedonic status evaluated by decrease in sucrose intake. Thirty six rats were used in the present large-scale transcriptomic study, nine of which were stress-unchallenged controls treated with vehicle. Eighteen from 27 anhedonic-like CMS rats were exposed to escitalopram (5 mg/kg/day) treatment. Nine rats were injected with saline and used as treatment controls.

The animals were singly housed, except when grouping was applied as a stress parameter. Food and water was available ad libitum except when food and/or water deprivation was applied as a stress parameter. The standard 12-h light/dark cycle was only changed in course of stress regime.

During the CMS protocol all animals with symptoms of anhedonia demonstrated decreased weight gain and cognitive impairments as was shown repeatedly before [16], [18], [19].

Sucrose consumption test

The animals were first habituated to consume a palatable sucrose solution (1.5%). The habituation period lasted for five weeks. In this period, the sucrose test was made twice a week during the first three weeks and once a week during the last two weeks. Animals were food and water deprived 14 h prior to the test, which was a one hour period with free access to a bottle of the sucrose solution. During the stress period the sucrose consumption test was performed once a week.

Chronic mild stress protocol

On the basis of sucrose consumption in the two final baseline tests, the animals were divided into two matched groups and placed in separate rooms. One group was exposed to an initial four weeks of chronic mild stressors and the other one was left undisturbed. The stress procedure was performed according to a procedure optimized in our laboratory [19]. Briefly, rats were submitted to one period of intermittent illumination, stroboscopic light, grouping, food or water deprivation; two periods of soiled cage and no stress; and three periods of 45° cage tilting. During grouping, rats were housed in pairs with different partners, with the individual rat alternately being a resident or an intruder. All the stressors lasted from 10 to 14 hours.

Following exposure to stress, rats were characterized as being anhedonic-like (defined as a >30% within-subject decrease in sucrose intake, also denoted as an anhedonic index value of 0.70 or less when compared to baseline) or resilient (defined as a <10% within-subject decrease in sucrose intake). Rats not responding to either criterion were excluded from the experiment. Stress-resilient rats were not included in the present study. The proportion of anhedonic-like, stress-resilient and intermediate grouped rats is shown in a recent publication by Henningsen et al. [26].

After the initial four weeks of exposure to stress, the stress group was divided into two matched subgroups and subjected to chronic escitalopram and vehicle administration, respectively, for four weeks. Stress was continued during the entire period of treatment. During treatment, the drug-treated groups segregated into drug responders and drug non-responders. We have in a previous study shown that escitalopram-treated CMS animals follow a bimodal segregation into two subgroups of responders and non-responders based on readouts from all four weeks of treatment [19].

Drug or vehicle was administered intraperitoneally once a day in the morning. Escitalopram (Sequoia Research Products Limited, United Kingdom) was dissolved in saline and given at dosages of 5 mg/kg/day.

Following the exposure of antidepressant treatment rats were characterized as drug responders (defined as a≥10% within-subject increase in sucrose intake, if the readout is above the anhedonic index value of 0.70) or drug non-responders (defined as <10% increase in sucrose intake). Rats not responding to either criterion were excluded from the experiment. Individual CMS groups (n = 9) were divided into three subgroups (n = 3) (Table 1).

thumbnail
Table 1. Animal groups investigated in the present CMS study.

https://doi.org/10.1371/journal.pone.0080666.t001

Tissue processing

Rats were decapitated and brains were removed and immediately frozen on dry ice. Frozen brains were sectioned coronally (i.e. −2.12 to −4.16 mm relative to Bregma; [30] on a cryostat (CM3050S, Leica Microsystems, GmbH, Germany). 50 µm sections were collected throughout the lateral habenula. Immediately after sectioning, each tissue section was mounted on polyethylene napthalate (PEN) glass slides (Molecular Devices, USA) and stored at −80°C until further processing.

Tissue staining

Slides were removed from −80°C storage and placed in 96% ethanol for one minute. Tissue sections were then stained in 1% cresyl violet (CV) (dissolved in absolute ethanol) for 15–25 s. Slides were then dehydrated in graded alcohols; 30s in 96% ethanol, 30 s in absolute ethanol and finally one minute in a second ethanol (absolute) step. Tissue sections were air-dried for a couple of minutes prior to laser capture microdissection (LCM).

Laser capture microdissection

The LCM technique was used to isolate homogeneous LHb tissues devoid of contamination from the surrounding thalamic and medial habenular tissues. LCM was performed using the Veritas Microdissection Instrument model 704 (Molecular Devices, USA) with CapSure Macro caps (Molecular Devices, USA). The LHb was visualized in the microscope of the LCM instrument and captured by the “cut & capture” feature. The location of dissected area is demonstrated in Figure S1. The LCM procedure was recently described in a previous study [31]. Optimized settings were 80 mW pulse power, 3.500 µsec pulse duration, 45 µm laser spot diameter, and a UV laser power at 15. The LHb area was selected for capture using the 2× objective while the capture-part was performed on the 20× objective. Captured LHb tissue was removed from the caps by use of a pipette and placed in a 0.5 ml PCR tube containing 20 µl QIAzol buffer (Qiagen Inc., USA). The sample was centrifuged at 14.500 rpm for 30 s and stored at −80°C until RNA extraction.

RNA isolation and quality control

After tissue homogenization, total RNA was isolated from each subgroup by use of miRNeasy Mini Kit (Qiagen Inc., USA) according to the manufacturer's protocol. Prior to the microarray analyses, the isolated total RNA was assessed to confirm high quality with respect to integrity and purity. The RIN (RNA integrity number) value of our RNA samples, as measured by the 2100 BioAnalyzer (Agilent Technologies, USA) was 7 or higher, indicating high integrity of the total RNA used for downstream analyses. The quantity of RNA samples was measured by use of Nanodrop ND-1000 (Thermo Scientific, USA).

RNA amplification and microarray hybridization

In this study, gene expression profiles of rats exposed to CMS and antidepressant treatment (four different groups – table 1) were determined using whole genome-wide gene expression microarray analysis. All groups contained three biological replicates, therefore, a total of 12 microarrays were analyzed. Tissue from three individual rats was included in each biological replicate. Samples (15ng total RNA) were reverse transcribed and amplified using the Ribo-SPIA® Ovation™ Pico WTA system V2 (NuGEN, USA) according to the manufacturer's instructions. 5 μg cDNA was used for fragmentation and labeling by use of the FL-ovation™ cDNA biotin module (NuGEN, USA) according to the manufacturer's instructions. Labeled cDNA samples of 4 μg were hybridized to the Affymetrix GeneChip® High-throughput (HT) RG-230 Perfect Match (PM) 24-Array Plate for rat (Affymetrix, USA). Each array probes approximately 28.700 genes and contains more than 31.000 probes. Affymetrix Expression Console Software (Affymetrix, USA) was used for normalization (Probe Logarithmic Intensity Error Estimation – PLIER normalization strategy) of microarray data. Microarray results are available in the ArrayExpress database: accession E-MTAB-1857.

Verification of regulated genes in the lateral habenula by real-time qPCR

TaqMan® qPCR analysis was performed by AROS, ApS (Aarhus, Denmark) using Fluidigm® technology. All selected RNA specimens were the same as the ones used in the chip experiment. Primer sets and probes for 13 tested genes were obtained from Life Technologies (filial Denmark, Naerum, Denmark). NCBI gene classification and context sequences for the tested genes are shown in table S1. The standard curves for each gene amplification were used for quantification of gene expression. The average slope of curves was −3.31±0.07 demonstrating a high efficiency of all performed qPCR reactions.

Cluster Analysis

Cluster analysis was used to separate the four different groups based on between-groups similarities in global gene expression profiles. The clustering approach orders objects in a treelike structure and provides information about relations among groups. The cluster analyses were performed with Cluster software version 2.11 [32]. The first cluster analysis was based on all probes being expressed by the cDNA microarray chip. The second cluster analysis was based on the most regulated genes (13.976 genes) based on all four CMS groups calculated as standard deviation/mean (individual gene transcript intensity) > standard deviation/mean (all gene transcript intensities). Data were adjusted by log2 transformation, normalization and mean center of genes and arrays. The hierarchical clustering was performed by clustering of arrays followed by average linkage clustering. Subsequently, the cluster analysis was visualized by Tree View software version 1.60 (http://rana.lbl.gov/EisenSoftware.htm).

Statistical analyses

Overall stress effects and treatment effects, respectively, were investigated using Student's t-tests. Stress effects were estimated by comparing all stressed rats in the trial to all non-stressed rats (Figure 1A). Treatment effects were estimated by comparing all drug-treated rats (both responders and non-responders) to non-treated rats (Figure 1B). Data obtained from weekly sucrose tests were analyzed using one-way multivariate analysis of variance (MANOVA). Antidepressant effects were investigated separately, with stress and treatment, respectively, as between subject factors and time (weeks) as within subject factor. According to the statistically significant interactions revealed by MANOVA, additional analyses were performed by one-way analysis of variance (ANOVA) to reveal time-specific differences caused by stress- or treatment. Bonferroni corrections were performed in order to account for multiple ANOVAs being run. Significant ANOVAs were followed up with groupwise Bonferroni post-hoc tests. The statistical level of significance was set at p<0.05.

thumbnail
Figure 1. Effects of stress and treatment, respectively.

(A) Four weeks exposure to chronic mild stress (CMS rats, n = 249) resulted in a significantly decreased sucrose consumption when compared to unchallenged controls (n = 80). (B) Four weeks exposure to antidepressant treatment with escitalopram resulted in a significantly increased sucrose consumption when comparing all escitalopram-treated CMS rats (responders and non-responders, respectively, n = 18) with vehicle-treated anhedonic-like CMS-rats (CMS-V, n = 9). (C) Four weeks of treatment with escitalopram resulted in a significant increase in sucrose consumption in a subgroup of drug-treated rats, i.e. rats responding to escitalopram (Esc-R, n = 9) treatment. A subgroup of the treated rats did not increase in sucrose intake (CMS escitalopram non-responders, Esc-NR, n = 9). Statistical differences were set at *p<0.05 and **p<0.001. Data are shown as mean sucrose intake, indexed to baseline (before stress) and stress baseline (before treatment), respectively + standard error of mean.

https://doi.org/10.1371/journal.pone.0080666.g001

Differential changes in gene expression were investigated by Student's t-test to reveal significant differences. However, due to a high content of low regulations (with maximum regulations of 30–40%), significant genes (p<5%) were investigated collectively according to biological pathways.

Changes in gene expression detected by qPCR where analyzed in XLSTAT by one-way ANOVA using two-sided Dunnett's post-hoc comparison to identify significant differences between groups.

Functional pathway analysis of differentially regulated genes

Genes that were significantly altered by CMS and antidepressant treatment were functionally classified into cellular pathways using the Protein ANalysis THrough Evolutionary Relationships (PANTHER, www.pantherdb.org) system. This system classifies genes on the basis of their functions taking published experimental evidences and otherwise evolutionary relationships to predict functions. To assess the statistical enrichment or over-representation of these cellular pathways from our dataset relative to the global set of rat genes, a binomial statistics within the PANTHER system was applied. The statistical level of significance in PANTHER was set at p<0.05. Genes investigated for functional pathway analyses had to fulfill a certain criterion; the gene had to be significantly (p<5%) regulated between relevant groups (Table 2). Pools of significant genes were thus submitted to functional pathway analysis in PANTHER in order to look for pathways implicated in stress and antidepressant responses. Several pathways were shown to be implicated in stress and antidepressant responses. In order to limit the amount of pathways in the present study, pathways were only of interest if they were implicated in both the induction and recovery of depression. As illustrated in table 3, pathways have to be significantly changed when comparing U-V to CMS-V, but also when comparing CMS-V to Esc-R. The two distinct pathway regulation patterns in table 3, indicates whether pathways are affected when comparing the unchallenged healthy state (U-V) to a recovered state (Esc-R). Examples of genes within the individual pathways, affected by stress and treatment (Esc-R), respectively, are listed in table 4.

thumbnail
Table 2. Functional pathway analysis by the Protein Analysis through Evolutionary Relationships (PANTHER) system.

https://doi.org/10.1371/journal.pone.0080666.t002

thumbnail
Table 3. Biological pathways implicated in depression/recovery (Figure 5A–B).

https://doi.org/10.1371/journal.pone.0080666.t003

Similarly, numerous pathways are implicated in treatment responses. The amount of pathways affected by antidepressant treatment (CMS-V versus Esc-R) is reduced by selecting pathways with concomitant changes in escitalopram responding rats (Esc-R) versus escitalopram non-responding rats (Esc-NR), thus to focus on pathways related to treatment response, but with diverging effects in drug-responders (EscR) and drug-non-responders (Esc-NR), respectively. The two regulation patterns in table 5 indicate whether pathways are affected when comparing the anhedonic-like state (CMS-V) with rats not responding to antidepressant treatment with escitalopram (Esc-NR). Examples of genes within the individual pathways related to treatment response (Esc-R and Esc-NR respectively) are listed in table 6.

thumbnail
Table 6. Pathway categories implicated in treatment response.

https://doi.org/10.1371/journal.pone.0080666.t006

Results

Sucrose consumption

The sucrose consumption test was applied to assess stress-induced anhedonic-like behavior and antidepressant-induced recovery. Four weeks exposure to CMS resulted in significantly decreased sucrose consumption (p<0.001) when comparing all stress exposed rats (n = 249, anhedonic-like, intermediate and stress-resilient rat groups) to unchallenged controls (n = 80), (Figure 1A). Similarly, four weeks of antidepressant treatment with escitalopram resulted in significantly increased sucrose consumption (p<0.05) when comparing all escitalopram-treated anhedonic-like rats (both responders and non-responders, n = 18) with vehicle-treated CMS-rats (n = 9) (Figure 1B).

Figure 1C shows the sucrose intake (indexed to stress baseline) for the three experimental groups; CMS vehicle (CMS-V), CMS escitalopram responders (Esc-R), and CMS escitalopram non-responders (Esc-NR). The groups were defined according to changes in sucrose intake in response to exposure to chronic mild stressors and chronic antidepressant treatment with escitalopram. Antidepressant treatment proceeded for four weeks concurrently with stress exposure. As treatment progressed, chronic treatment with escitalopram reversed the decrease in sucrose intake in about 50% of the treated rats (responders).

MANOVA analysis on escitalopram treatment effects revealed significant differences over time (F8,42 = 3.96, P = 0.001), when comparing CMS-V (n = 9) to treatment responders (n = 9) and treatment non-responders (n = 9). Significant effects of drug treatment was present after one week of escitalopram administration (not shown) and retained in the entire period of treatment (F2,24 = 14.11, PANOVA <0.001, Ppost hoc <0.001) Significant segregation into drug responders and drug non-responders was achieved after three weeks of escitalopram treatment (not shown) (F2,24 = 14.11, PANOVA <0.001, Ppost hoc  = 0.001).

Transcriptomic analysis

The laser-capture microdissected LHb RNA was subjected to microarray analysis to identify transcriptomic changes related to depression and antidepressant treatment. 19.010 genes were up regulated after CMS while 12.129 genes were down regulated after CMS (Table 7). Similarly, 15.334 genes were up regulated after chronic treatment with escitalopram, while 15.805 genes were down regulated. However, this was also the case for rats responding to escitalopram treatment (Esc-R). On the contrary, escitalopram non-responders (Esc-NR) demonstrated treatment-associated up regulations of 10.441 genes while 20.698 genes were down regulated after antidepressant treatment. Ratios of up- and down regulations were similar also when focusing at the most regulated genes as illustrated in table 7.

thumbnail
Table 7. Gene regulations after CMS and antidepressant treatment.

https://doi.org/10.1371/journal.pone.0080666.t007

Hierarchical clustering analysis

The TreeView cluster analysis in Figure 2A was based on all probes on the microarray chip representing about 31.000 genes while the cluster analysis in figure 2B was based on the most regulated genes of the microarray data (i.e. standard deviation for a specific gene transcript normalized with mean expression level for the same gene is expected to be higher than the standard deviation normalized with the mean expression level for all genes). The global clustering analysis based on all genes indicated a clear differentiation between groups receiving antidepressant treatment or no treatment. However, clustering according to hedonic status was not evident from this cluster analysis. When including only the most regulated genes among the four groups in the analysis, treatment-associated effects appears, thus segregating rats responding and not responding to antidepressant treatment. Escitalopram non-responders are thus markedly diverging from all other groups in line with the results from the transcriptomic analyses (Table 7).

thumbnail
Figure 2. Hierarchical clustering analysis.

(A) Based on all probes represented on the cDNA microarray chip showing a clear differentiation between groups receiving antidepressant treatment (Esc-R + Esc-NR) or no treatment (U-V, CMS-V). (B) Based on most regulated genes (approx. 14.000 genes – defined as gene transcripts having a standard deviation/mean (individual gene transcript intensity) > standard deviation/mean (all gene transcript intensities)), showing a clear divergence of escitalopram non-responders (Esc-NR). U-V, unchallenged control rats; CMS-V, CMS vehicle; Esc-R, CMS escitalopram responders; Esc-NR, CMS escitalopram non-responders.

https://doi.org/10.1371/journal.pone.0080666.g002

Venn diagrams

First, we indentified genes involved in both the development of the anhedonia-like condition and antidepressant treatment-associated recovery. Four independent comparisons were performed to approach gene differentiation. Venn diagrams illustrated in Figure 3 demonstrate the models for comparison of genes involved in aforementioned processes. There are 88 genes differentiating CMS and stress unchallenged groups; 49 of them are identified in the NCBI database (Figure 3A and S2A, left panels). Comparison between CMS and treatment responder groups revealed changes in 138 genes, 77 of which are identified in the NCBI database (Figure 3A and S2A, right panels). Intersection area of Venn diagram in Figure 3A contains 6 genes and 3 of them are identified in NCBI database and up-regulated in both comparisons. Antidepressant treatment responders and non-responders differ in 565 genes, where 380 were identified in the NCBI database. A list of 46 genes from this comparison is shown in Figure S2B. Two comparisons describing differences in gene expression specific to response to antidepressant treatment revealed 13 common genes (Figure 3B). Ten of them are identified. Fcrla, Eif3k, Sec3l1, and Ubr5 were up-regulated in comparison to both CMS and escitalopram non-responders groups; levels of Abca8a, Ankrd49, Cyp2j10, Frs3, Syn2, Znf503 were lower in escitalopram responders compared to CMS and non-responders groups.

thumbnail
Figure 3. Venn diagrams for significant gene regulations associated to anhedonia and treatment recovery

. (A) Associated gene differences in comparisons between vehicle-treated anhedonic-like (CMS-V) versus unchallenged control groups (left panel) and CMS-V versus treatment responders (right panel); (B) associated gene differences according to response to treatment between escitalopram responders versus CMS-V (left panel) and versus treatment non-responders (right panel). Results in parentheses indicate amount of genes identified in the NCBI database. All identified genes present in the intersection areas are listed by gene symbols. U-V, vehicle treated unchallenged control rats; Esc-R, CMS escitalopram responders; Esc-NR, CMS escitalopram non-responders. * – significant changes by qPCR for left side comparisons; + – significant changes by qPCR for right side comparisons.

https://doi.org/10.1371/journal.pone.0080666.g003

Real-time qPCR

Secondly, 13 genes from intersection areas of both Venn diagrams (Figure 3) were used in qPCR experiments to confirm changes in gene expression discovered by microarray analysis. Results presented in Figure 4A demonstrate up-regulation of three genes in CMS rats compared to both unchallenged controls and escitalopram responders and thus confirming findings discovered in the microarray analysis. Figure 4B demonstrates qPCR revealed changes in expression in treatment responders versus anhedonic-like rats (CMS-V) and treatment non-responders rats, respectively. Interestingly, qPCR quantification resulted in relatively large deviations in all genes in treatment non-responders group; therefore the 95% of significance in comparison between groups was achieved only for Cyp2j10 and Eif3k genes, but the directions in most gene regulations were the same as in results of microarray analysis. Groupwise relationships according to the response to stress and response to antidepressant treatment are shown in Figure S3. CMS led to up-regulation of all tested genes according to the stress-unchallenged controls and positive response to antidepressant treatment was associated with gene expression recovery (Figure S3A). Comparison according to the response to treatment revealed a strong up-regulation of tested genes in both CMS and treatment non-responders groups with confirmed significant changes for most of the genes in the CMS group (Figure S3B).

thumbnail
Figure 4. Expression levels of the NCBI identified genes confirmed by qPCR.

(A) Genes changed in comparisons between anhedonic-like versus unchallenged controls and anhedonic-like versus treatment responders groups. (B) Genes regulated in comparisons between treatment responders versus anhedonic-like and treatment responders versus treatment non-responders groups. U-V, vehicle treated unchallenged control rats; CMS-V, vehicle treated anhedonic-like rats; Esc-R, CMS escitalopram responders; Esc-NR, CMS escitalopram non-responders. Red crosses indicate arithmetical mean values; *, **, *** – p<0.05, 0.01, 0.001, respectively, by Dunnett's ANOVA post-test: A – in comparison with CMS-V; B – in comparison with Esc-R.

https://doi.org/10.1371/journal.pone.0080666.g004

Biological pathway analysis

Due to relatively small stress- and treatment-associated gene regulations at single gene levels, significantly regulated genes (p<5%) were analyzed according to implications in biological pathways. As illustrated in table 3 and table 5, most pathways associated with recovery and treatment response belong to the two categories of pathways; growth factor receptor signaling and G protein–coupled receptor signaling. Furthermore, stress and escitalopram treatment affect additional intracellular pathways like integrin signaling, hedgehog signaling, and chemokine and cytokine signaling. The genes listed in table 4 and table 6 constitutes only a fraction of the genes represented in the pathways. Genes listed in the tables are randomly chosen from the entire set of genes in the individual pathways.

Discussion

In the present study we used the validated rat CMS model of depression to mirror symptoms of human depression and to investigate molecular correlates of depression and antidepressant treatment in the lateral habenula. In this study there was a general effect of four weeks of mild stress exposure on sucrose intake with a mean decrease of 29% (Figure 1A). As previously reported, rats segregate into two groups when exposed to CMS, one group sensitive to CMS, thus reflecting anhedonic-like behavior, and one group resilient to CMS [18], [25][27]. Another unique feature of the CMS model of depression is that only a fraction of the drug-treated rats respond to antidepressants, thus leaving the remaining rats in the anhedonic-like state (CMS-Esc-NR). This phenomenon is well known from the clinic as only a fraction of the depressive patients respond to antidepressant medication. The proportion of drug responding/non-responding rats constitutes approximately 50%, thus allowing us to investigate molecular mechanisms of drug-induced recovery and therapy refraction, respectively. This unique individual treatment response as seen in the present study is reproducible in all our CMS studies [19], [23][25], [27], and has not been reported in CMS studies elsewhere [17], [20], [23], [33]. The proportion of escitalopram responders in the CMS model of depression correlates to clinical response rate ranging from approximately 60% to 85% in patients suffering from moderate or severe depression disorders [34], [35]. The placebo effect observed in clinical studies is 25–35% and may explain the higher response rates observed in clinical studies when compared to animal studies [36]. As deduced from the present CMS experiment, there is a general effect of antidepressant treatment with escitalopram after four weeks of treatment with a mean increase in sucrose intake of 43% (Figure 1B). This increase is even more pronounced in the escitalopram responder group (with a mean increase in sucrose intake of 73%) as illustrated in Figure 1C. On the other hand, antidepressant treatment had no significant effect on rats not responding to treatment illustrated by unchanged drinking behavior in the sucrose consumption test. Escitalopram displayed ameliorating effects on CMS-induced anhedonia (in the responder group – Esc-R) after one week of treatment (not shown), while full recovery was obtained after four weeks of treatment. The variable treatment response seen is further validated in the hierarchical clustering analysis based on the most regulated genes (approximately 14.000 genes), where the escitalopram non-responder group, in specific, diverges from all other CMS-generated phenotypes, indicating major molecular aberrations in treatment-resistant rats. Thus, in addition to modeling recovery, the CMS model of depression also appears to imitate human therapy refraction.

In the present study, a genome wide microarray analysis in CMS- and/or treatment-induced phenotypes was performed to assess effects on global transcriptomic profiles. This approach was chosen to quantitatively assess stress-related behaviors and antidepressant treatment of genes in order to investigate the molecular role of the LHb. The aim of the study was to narrow down the set of genes that were differentially regulated by stress and antidepressant treatment. Although many genes were found to be regulated most of them were only modestly regulated (p<0.05 and maximal fold changes of 30–40%). Therefore, differentially regulated genes were investigated collectively according to biological pathways. Stress- and treatment-associated gene regulations at this level is a common phenomenon confirming the statement from Mirnics et al. [37], that gene expression changes in psychiatric traits are small and that psychiatric diseases may result from cumulative subtle differences. Additionally, we have previously identified small gene regulations in a similar transcriptomic study investigating CMS- and treatment-induced changes in a hippocampal subregion [25]. Considering the small gene regulations in the present study and the high risk for false positives of individual genes, our strategy was to search for biological pathways implicated in depression and antidepressant treatment. We believe that clusters of genes pointing towards specific biological processes is the most convincing way to identify targets related to recovery and drug response under the present circumstances.

Constraints were applied for analysis of differential pathway regulation profiles and thus for the identification of genes/pathways involved in depression etiology (anhedonia) and recovery (Figure 5). Hence, pathways were prioritized if they were significantly regulated after exposure to CMS, and concomitantly also significantly affected by chronic treatment with escitalopram in the responder group (Esc-R) (Table 3). Using the applied constraint reduces the risk of false positives and allows us to focus on treatment effects specifically associated with recovery, thus excluding the non-specific (non-anhedonia-related) drug effects of escitalopram. According to the applied constraints we found 11 distinct pathways to be involved in the etiology/recovery of depression, two of which are affected when comparing unchallenged controls (U-V) with healthy drug-recovered (Esc-R) rats (Table 3). Among pathways involved in depression pathology/recovery, we found several intracellular signal transduction pathways. One such pathway category is growth factor receptor signaling (vascular endothelial growth factor (VEGF) signaling, epidermal growth factor (EGF) signaling and angiogenesis). Among growth factors, more specifically, brain-derived neurotrophic factor (BDNF) has been suggested as a strong candidate modulating stress-associated pathology and antidepressant treatment. This have led to the development of a neurotropic hypothesis of depression suggesting that stress-induced decreases of BDNF expression, and possibly other growth factors, contributes to the progression of depression and that up-regulation of BDNF plays a role in the actions of antidepressant treatment [38][40]. VEGF is another growth factor, which has gained lots of attention during the last decade with respect to depression and antidepressant treatment. VEGF was first identified as an angiogenic factor essential for the formation of new blood vessels in physiological and pathological processes [41], [42]. In addition to stimulating endothelial cell proliferation, migration, and survival recent studies established a neurotrophic and neuroprotective role of VEGF [43], [44]. Just like BDNF, VEGF also stimulates adult neurogenesis and a decreased level of neurogenesis as a consequence of chronic stress associates with decreased VEGF levels [43], [45]. Furthermore, exposing rats to the CMS model of depression resulted in decreased VEGF expression in the hippocampus [23]. On the other hand multiple classes of antidepressants and ECS increase the expression of VEGF [46]-[48] and it has been show that blockage of VEGF expression blocks the beneficial effects of antidepressants in animals subjected to unpredictable chronic mild stress [49], [50]. In the present study we found that the VEGF signaling pathway was affected during depression and antidepressant treatment. Genes affecting this pathway includes vegfa and Flt1 (VEGF receptor 1), and downstream intracellular pathways of the growth factor receptor signaling pathway, like protein kinase C (Prkca, Prkd3 and Prkci) and signal transducers and activators of transcription 3, 4, 6 (Stat3,Stat4, Stat6)(Table 4). Another prominent intracellular pathway identified in the present study with respect to depression and antidepressant treatment is the Wnt signaling pathway. Wnt's are secreted glycoproteins that signal through the frizzled (Fz) receptors which couple to several signaling cascades. Wnt signaling is crucial for embryonic development and disruption of this cascade has been associated with various neurodevelopmental and neuropsychiatric disorders [51], [52]. More specific, a role for the Wnt signaling pathway in depression-related behavior has been reported in several studies. Stress is associated with imbalances of Wnt signaling molecules and several of these molecules are regulated by various classes of antidepressants [53]. Genes affecting the Wnt pathway in the present study includes the intracellular signaling molecules Axin1, secreted frizzeled-related protein 2 (sfrp2) and segment polarity protein dishelveled homolog DVL-1 (dvl1) (Table 4). In fact, Wnt's and growth factors/neurogenic factors share common intracellular downstream cascades, targeting both the mitogen-activated protein kinase/extracellular signal-regulated kinase (MAPK/ERK) and phosphatidylinositol-3-kinase (PI3K) cascades. Both of these intracellular pathways were shown to be disrupted in subjects with depression and in animal models of depression [53]. Integrins, transmembrane receptors that mediate the attachment between a cell and the tissue that surrounds it, also appeared as a pathway of interest under the present constraints. Research combining depression and integrin signal transduction is scarce; however, intracellular signals via integrins are related to neuroplastic mechanisms like cell growth and division, cellular differentiation and apoptosis [54]. The integrin-mediated signaling pathway converges with intracellular pathways mediated by growth factors, resulting in transcriptional stimulation via the MAPK/ERK cascade. In the present study, changes in integrin-mediated signaling was mediated via an integrin ligand (fibronectin, Fn1), integrin receptors (integrin alpha-7 light, Itga7 and integrin alpha-M, Itgam), and downstream signaling molecules like proto-oncogene tyrosine-protein kinase Fyn (Fyn) and son of sevenless homolog 1 (Sos1) (Table 4). All three pathways mentioned in the present study: growth factor receptor signaling, Wnt signaling and integrin-mediated signaling, converge on activating the transcription factor cAMP response element-binding protein (CREB), which is known to regulate the transcription of several genes implicated in depression pathology and treatment response [53], [54]. Furthermore, CREB has repeatedly been shown to be modulated in depressive-like states and after antidepressant treatment [55]. In fact, we also found that transcriptional regulation was significantly regulated as a pathway in the present study, including the transcription factor processed cyclic adenosine monophosphate (cAMP)-responsive element-binding protein 3-like protein 2 (creb3l2), indicating depression and treatment-associated changes even at a transcriptional level of signal transduction (Table 4).

thumbnail
Figure 5. Relationships between groups built up for PANTHER analysis.

Statistical enrichment of biological pathways was assessed by binomial statistics within the PANTHER system. (A and B) Relationships for classification of pathways implicated in depression and recovery; (C and D) Relationships for classification of pathways implicated in treatment response. Dashed lines indicate significant differences in pathways between groups (p<0.05), regular lines indicate no differences in pathways between groups. U-V, unchallenged control rats; CMS-V, CMS vehicle; Esc-R, CMS escitalopram responders; Esc-NR, CMS escitalopram non-responders.

https://doi.org/10.1371/journal.pone.0080666.g005

The segregation of drug-treated CMS rats into responder (Esc-R) and non-responder (Esc-NR) groups, allows us to search for molecular pathways implicated in the individual treatment response seen in the CMS model of depression. Pathways of interest will thus reflect mechanisms affected by antidepressant treatment (CMS-V versus Esc-R), but also mechanisms showing divergence between rats responding to escitalopram treatment and rats not responding to antidepressant treatment (Esc-R versus Esc-NR; Figure 5C-D). According to the applied constraints we found 24 distinct pathways to be involved in treatment response, seven of which are affected when comparing anhedonic-like rats (CMS-V) with escitalopram non-responders (Esc-NR) rats (Table 5). A major part of pathways implicated in treatment response belongs to the G-protein-coupled receptor-signaling pathway. Molecular pathways represented in this pathway category includes metabotropic glutamate receptor signaling, dopamine receptor-mediated signaling, 5HT1 type receptor-mediated signaling and GABA-B receptor signaling, among others (Table 5) – all mediating their physiological functions via G-proteins. Common for this category of pathways is that ligand binding to the G protein coupled receptors activate or inhibit the adenylyl cyclase (AC)-cAMP signal transduction pathway. AC is the critical enzyme to catalyze cAMP syntese. cAMP is the ubiquitous intracellular second messenger, which stimulates activities of cAMP-dependent protein kinases (e.g. protein kinase A – PKA) and CREB. In line with CREB alterations during depression, PKA levels are found to be reduced in depressed suicide subjects and in animal models of depression [56]. Evidence from a number of different studies has demonstrated that antidepressant treatments up-regulate the AC-cAMP signal transduction pathway, including increased levels of Gαs coupling to AC and increased levels of PKA leading to increased CREB function [57]. In fact, in the present study we found that Creb1 is affected by antidepressant treatment as well as cAMP-dependent protein-kinase type-I-alpha regulatory subunit (Prkar1a), adenylate cyclase 2, 3 and 5 (Adcy2, Adcy3, Adcy5) and several G proteins; Gnai3, Gnb1, Gng10 (Table 6). Additionally, we found other intracellular pathways implicated in treatment response, like hedgehog signaling, integrin signaling, growth factor receptor signaling, and chemokine and cytokine signaling (Table 5).

In general, the present microarray study led to molecular findings associated with depression and antidepressant treatment in the LHb. Due to low regulations at gene level, we applied a strategy searching for pathways implicated in the development and recovery from depression. Pathways involved in intracellular signal transduction are highly represented in the present study. As deduced from above, alterations in intracellular signaling and components involved in these processes have repeatedly been associated with the disorder of depression and its treatment. It is known that intracellular signaling is a mediator of neuroplasticity and molecular mechanisms of neuroplasticity indeed overlap molecular mechanisms of depression and antidepressant response. However, most depression-related studies focus on research on the frontal cortex or the hippocampus as dysfunctions in both of these brain regions are highly associated with depression and treatment. To our knowledge, the present study is the first to investigate molecular mechanisms of depression-like behavior and treatment in the LHb of rats exposed to an animal model of depression. Deduced from our study, it seems that stress and antidepressant treatment primarily modulate signal transduction pathways in LHb, leading us to speculate that increased intracellular signaling contributes to increases in neuronal activities during depressive-like states. The conclusions assumed from the present study are speculative and based on pathway analyses on significant gene alterations. In order to confirm the present results, it is necessary to validate that the individual intracellular signaling systems in fact are disturbed under depression and antidepressant treatment. Monoaminergic systems like serotonin, dopamine, and noradrenalin neurotransmitter systems are deregulated in depressive-like states. The LHb directly innervates these transmitter systems, indicating a central role for this brain region during depression and antidepressant treatment. A dysfunction of afferent and efferent LHb innervation might thus create imbalances in monoaminergic neurotransmission leading to depressive-like states or vice versa. Unraveling the role of the LHb in depression and treatment might lead us to new molecular targets involved in depression pathology and novel treatment strategies.

Supporting Information

Figure S1.

Coronal slice image with the selected lateral habenula area. DG, dentate gyrus; D3V, Dorsal third ventricle; mHb, medial habenula; LHb, lateral habenula.

https://doi.org/10.1371/journal.pone.0080666.s001

(TIF)

Figure S2.

Venn diagrams of NCBI identified genes. (A) In between groups comparisons for CMS-V versus U-V and CMS-V versus Esc-R; (B) In between groups comparisons for Esc-R versus CMS-V and Esc-R versus Esc-NR. In bold font are written symbols for genes up-regulated in the current comparison, in italic font are written symbols of genes down-regulated in the current comparison. * – significance of changes confirmed by qPCR for left-side comparison; + – significance of changes confirmed by qPCR for right-side comparison. CMS-V, vehicle treated anhedonic-like group; U-V, vehicle treated CMS unchallenged group; Esc-R, escitalopram treatment responders group; Esc-NR, escitalopram treatment non-responders group.

https://doi.org/10.1371/journal.pone.0080666.s002

(TIF)

Figure S3.

Stress and anti-depressant effects of genes measured by qPCR. (A) Stress-induced up-regulation of tested genes demonstrated as a ratio of gene expression between CMS-V and U-V groups; treatment-induced recovery demonstrated as a ratio of gene expression between Esc-R and U-V groups. (B) According to response to escitalopram treatment (ratio to Esc-R) all tested genes were up-regulated in both CMS-V and Esc-NR groups. *, **, *** – p<0.05, 0.01, 0.001, respectively, by Dunnet's ANOVA post test: A – in comparison with U-V; B – in comparison with Esc-R. CMS-V, vehicle treated anhedonic-like group; U-V, vehicle treated CMS unchallenged group; Esc-R, treatment responders group; Esc-NR, treatment non-responders group.

https://doi.org/10.1371/journal.pone.0080666.s003

(TIF)

Table S1.

Database related characteristics for 13 genes used to confirm gene expression results by qPCR.

https://doi.org/10.1371/journal.pone.0080666.s004

(DOCX)

Acknowledgments

We thank Stine Dhiin Heide Hansen for her skilled animal handling and care.

Author Contributions

Conceived and designed the experiments: TC LJ OW. Performed the experiments: TC LJ. Analyzed the data: TC EVB. Contributed reagents/materials/analysis tools: TC OW EVB. Wrote the paper: TC EVB.

References

  1. 1. Krishnan V, Nestler EJ (2008) The molecular neurobiology of depression. Nature 455: 894–902.
  2. 2. Hikosaka O, Sesack SR, Lecourtier L, Shepard PD (2008) Habenula: crossroad between the basal ganglia and the limbic system. J Neurosci 28: 11825–11829.
  3. 3. Hikosaka O (2010) The habenula: from stress evasion to value-based decision-making. Nat Rev Neurosci 11: 503–513.
  4. 4. Wirtshafter D, Asin KE, Pitzer MR (1994) Dopamine agonists and stress produce different patterns of Fos-like immunoreactivity in the lateral habenula. Brain Res 633: 21–26.
  5. 5. Caldecott-Hazard S, Mazziotta J, Phelps M (1988) Cerebral correlates of depressed behavior in rats, visualized using 14C-2-deoxyglucose autoradiography. J Neurosci 8: 1951–1961.
  6. 6. Shumake J, Edwards E, Gonzalez-Lima F (2003) Opposite metabolic changes in the habenula and ventral tegmental area of a genetic model of helpless behavior. Brain Res 963: 274–281.
  7. 7. Morris JS, Smith KA, Cowen PJ, Friston KJ, Dolan RJ (1999) Covariation of activity in habenula and dorsal raphe nuclei following tryptophan depletion. Neuroimage 10: 163–172.
  8. 8. Roiser JP, Levy J, Fromm SJ, Nugent AC, Talagala SL, et al. (2009) The effects of tryptophan depletion on neural responses to emotional words in remitted depression. Biol Psychiatry 66: 441–450.
  9. 9. Berton O, McClung CA, DiLeone RJ, Krishnan V, Renthal W, et al. (2006) Essential role of BDNF in the mesolimbic dopamine pathway in social defeat stress. Science 311: 864–868.
  10. 10. Yang LM, Hu B, Xia YH, Zhang BL, Zhao H (2008) Lateral habenula lesions improve the behavioral response in depressed rats via increasing the serotonin level in dorsal raphe nucleus. Behav Brain Res 188: 84–90.
  11. 11. Meng H, Wang Y, Huang M, Lin W, Wang S, et al. (2011) Chronic deep brain stimulation of the lateral habenula nucleus in a rat model of depression. Brain Res 1422: 32–38.
  12. 12. Winter C, Vollmayr B, Djodari-Irani A, Klein J, Sartorius A (2011) Pharmacological inhibition of the lateral habenula improves depressive-like behavior in an animal model of treatment resistant depression. Behav Brain Res 216: 463–465.
  13. 13. Sartorius A, Kiening KL, Kirsch P, von Gall CC, Haberkorn U, et al. (2010) Remission of major depression under deep brain stimulation of the lateral habenula in a therapy-refractory patient. Biol Psychiatry 67: e9–e11.
  14. 14. Willner P, Muscat R, Papp M (1992) Chronic mild stress-induced anhedonia: a realistic animal model of depression. NeurosciBiobehavRev 16: 525–534.
  15. 15. Willner P (2005) Chronic mild stress (CMS) revisited: consistency and behavioural-neurobiological concordance in the effects of CMS. Neuropsychobiology 52: 90–110.
  16. 16. Henningsen K, Andreasen JT, Bouzinova EV, Jayatissa MN, Jensen MS, et al. (2009) Cognitive deficits in the rat chronic mild stress model for depression: relation to anhedonic-like responses. BehavBrain Res 198: 136–141.
  17. 17. Grippo AJ, Francis J, Beltz TG, Felder RB, Johnson AK (2005) Neuroendocrine and cytokine profile of chronic mild stress-induced anhedonia. Physiol Behav 84: 697–706.
  18. 18. Christiansen S, Bouzinova EV, Palme R, Wiborg O (2012) Circadian activity of the hypothalamic-pituitary-adrenal axis is differentially affected in the rat chronic mild stress model of depression. Stress 15: 647–657.
  19. 19. Jayatissa MN, Bisgaard C, Tingstrîm A, Papp M, Wiborg O (2006) Hippocampal cytogenesis correlates to escitalopram-mediated recovery in a chronic mild stress rat model of depression. Neuropsychopharmacology 31: 2395–2404.
  20. 20. Banasr M, Chowdhury GM, Terwilliger R, Newton SS, Duman RS, et al. (2010) Glial pathology in an animal model of depression: reversal of stress-induced cellular, metabolic and behavioral deficits by the glutamate-modulating drug riluzole. Mol Psychiatry 15: 501–511.
  21. 21. Holm MM, Nieto-Gonzalez JL, Vardya I, Henningsen K, Jayatissa MN, et al. (2011) Hippocampal GABAergic dysfunction in a rat chronic mild stress model of depression. Hippocampus 21(4): 422–33.
  22. 22. Grønli J, Fiske E, Murison R, Bjorvatn B, Sõrensen E, et al. (2007) Extracellular levels of serotonin and GABA in the hippocampus after chronic mild stress in rats. A microdialysis study in an animal model of depression. BehavBrain Res 181: 42–51.
  23. 23. Bergström A, Jayatissa MN, Thykjær T, Wiborg O (2007) Molecular pathways associated with stress resilience and drug resistance in the chronic mild stress rat model of depression: a gene expression study. JMolNeurosci 33: 201–215.
  24. 24. Bisgaard CF, Jayatissa MN, Enghild JJ, Sanchez C, Artemychyn R, et al. (2007) Proteomic investigation of the ventral rat hippocampus links DRP-2 to escitalopram treatment resistance and SNAP to stress resilience in the chronic mild stress model of depression. JMolNeurosci 32: 132–144.
  25. 25. Christensen T, Bisgaard CF, Wiborg O (2011) Biomarkers of anhedonic-like behavior, antidepressant drug refraction, and stress resilience in a rat model of depression. Neuroscience 196: 66–79.
  26. 26. Henningsen K, Palmfeldt J, Christiansen S, Baiges I, Bak S, et al.. (2012) Candidate hippocampal biomarkers of susceptibility and resilience to stress in a rat model of depression. Mol Cell Proteomics 11: M111 016428.
  27. 27. Bisgaard CF, Bak S, Christensen T, Jensen ON, Enghild JJ, et al. (2012) Vesicular signalling and immune modulation as hedonic fingerprints: proteomic profiling in the chronic mild stress depression model. J Psychopharmacol 26: 1569–1583.
  28. 28. Berman RM, Narasimhan M, Charney DS (1997) Treatment-refractory depression: definitions and characteristics. Depress Anxiety 5: 154–164.
  29. 29. Fava M (2003) Diagnosis and definition of treatment-resistant depression. Biol Psychiatry 53: 649–659.
  30. 30. Paxinos G, Watson C (1998) The rat brain in stereotactic coordinates. Academic Press.
  31. 31. Christensen T, Bisgaard CF, Nielsen HB, Wiborg O (2010) Transcriptome differentiation along the dorso-ventral axis in laser-captured microdissected rat hippocampal granular cell layer. Neuroscience 170: 731–741.
  32. 32. Eisen MB, Spellman PT, Brown PO, Botstein D (1998) Cluster analysis and display of genome-wide expression patterns. ProcNatlAcadSciUSA 95: 14863–14868.
  33. 33. Li N, Liu RJ, Dwyer JM, Banasr M, Lee B, et al. (2011) Glutamate N-methyl-D-aspartate receptor antagonists rapidly reverse behavioral and synaptic deficits caused by chronic stress exposure. Biol Psychiatry 69: 754–761.
  34. 34. Ventura D, Armstrong EP, Skrepnek GH, Haim EM (2007) Escitalopram versus sertraline in the treatment of major depressive disorder: a randomized clinical trial. CurrMedResOpin 23: 245–250.
  35. 35. Ou JJ, Xun GL, Wu RR, Li LH, Fang MS, et al. (2011) Efficacy and safety of escitalopram versus citalopram in major depressive disorder: a 6-week, multicenter, randomized, double-blind, flexible-dose study. Psychopharmacology (Berl) 213: 639–646.
  36. 36. Quitkin FM, Rabkin JG, Gerald J, Davis JM, Klein DF (2000) Validity of clinical trials of antidepressants. AmJPsychiatry 157: 327–337.
  37. 37. Mirnics K, Middleton FA, Lewis DA, Levitt P (2001) Analysis of complex brain disorders with gene expression microarrays: schizophrenia as a disease of the synapse. Trends Neurosci 24: 479–486.
  38. 38. Duman RS, Heninger GR, Nestler EJ (1997) A molecular and cellular theory of depression. Arch Gen Psychiatry 54: 597–606.
  39. 39. Duman RS, Monteggia LM (2006) A neurotrophic model for stress-related mood disorders. Biol Psychiatry 59: 1116–1127.
  40. 40. Bath KG, Schilit A, Lee FS (2013) Stress effects on BDNF expression: Effects of age, sex, and form of stress. Neuroscience 239: 149–156.
  41. 41. Nowacka MM, Obuchowicz E (2012) Vascular endothelial growth factor (VEGF) and its role in the central nervous system: a new element in the neurotrophic hypothesis of antidepressant drug action. Neuropeptides 46: 1–10.
  42. 42. Ferrara N, Davis-Smyth T (1997) The biology of vascular endothelial growth factor. Endocr Rev 18: 4–25.
  43. 43. Ruiz de Almodovar C, Lambrechts D, Mazzone M, Carmeliet P (2009) Role and therapeutic potential of VEGF in the nervous system. Physiol Rev 89: 607–648.
  44. 44. Brockington A, Lewis C, Wharton S, Shaw PJ (2004) Vascular endothelial growth factor and the nervous system. Neuropathol Appl Neurobiol 30: 427–446.
  45. 45. Heine VM, Zareno J, Maslam S, Joels M, Lucassen PJ (2005) Chronic stress in the adult dentate gyrus reduces cell proliferation near the vasculature and VEGF and Flk-1 protein expression. Eur J Neurosci 21: 1304–1314.
  46. 46. Newton SS, Collier EF, Hunsberger J, Adams D, Terwilliger R, et al. (2003) Gene profile of electroconvulsive seizures: induction of neurotrophic and angiogenic factors. J Neurosci 23: 10841–10851.
  47. 47. Warner-Schmidt JL, Duman RS (2007) VEGF is an essential mediator of the neurogenic and behavioral actions of antidepressants. Proc Natl Acad Sci U S A 104: 4647–4652.
  48. 48. Warner-Schmidt JL, Duman RS (2008) VEGF as a potential target for therapeutic intervention in depression. Curr Opin Pharmacol 8: 14–19.
  49. 49. Greene J, Banasr M, Lee B, Warner-Schmidt J, Duman RS (2009) Vascular endothelial growth factor signaling is required for the behavioral actions of antidepressant treatment: pharmacological and cellular characterization. Neuropsychopharmacology 34: 2459–2468.
  50. 50. Fournier NM, Duman RS (2012) Role of vascular endothelial growth factor in adult hippocampal neurogenesis: implications for the pathophysiology and treatment of depression. Behav Brain Res 227: 440–449.
  51. 51. Logan CY, Nusse R (2004) The Wnt signaling pathway in development and disease. Annu Rev Cell Dev Biol 20: 781–810.
  52. 52. Inestrosa NC, Arenas E (2010) Emerging roles of Wnts in the adult nervous system. Nat Rev Neurosci 11: 77–86.
  53. 53. Voleti B, Duman RS (2012) The roles of neurotrophic factor and Wnt signaling in depression. Clin Pharmacol Ther 91: 333–338.
  54. 54. Wu X, Reddy DS (2012) Integrins as receptor targets for neurological disorders. Pharmacol Ther 134: 68–81.
  55. 55. Carlezon WA Jr, Duman RS, Nestler EJ (2005) The many faces of CREB. Trends Neurosci 28: 436–445.
  56. 56. Dwivedi Y, Pandey GN (2008) Adenylyl cyclase-cyclicAMP signaling in mood disorders: role of the crucial phosphorylating enzyme protein kinase A. Neuropsychiatr Dis Treat. 4: 161–176.
  57. 57. Pittenger C, Duman RS (2008) Stress, depression, and neuroplasticity: a convergence of mechanisms. Neuropsychopharmacology 33: 88–109.