Gene expression patterns in the hippocampus and amygdala of endogenous depression and chronic stress models (original) (raw)

Introduction

Major depressive disorder (MDD) is the third most prevalent and costly disease worldwide, and it is projected to become number one by 2030 (WHO, GBD report, 2004). In addition to increasing the risk for suicide, MDD presents frequent comorbidity with other psychiatric disorders, and it has been identified as a risk factor for many illnesses, including obesity, cardiovascular, neurodegenerative and other diseases.1, 2, 3, 4 Understanding the pathophysiology of MDD, and providing better treatments, would improve mental health and quality of life of the largest psychiatric population in the world.5

There are two major hypotheses with regard to the etiology of MDD: the monoamine deficiency and stress hypotheses of depression. The monoamine deficiency hypothesis proposes that MDD can be etiologically explained by a deficiency in the monoamine neurotransmitters serotonin, norepinephrine and dopamine. Although this hypothesis is still prominent and is the basis of most antidepressant development to date, current opinion is that monoamine deficiency only partially explains MDD.6, 7, 8 On the basis of observations in MDD, the stress hypothesis relates pathological alterations in the stress-responsive hypothalamic-pituitary-adrenal (HPA) axis to causality of depression.9, 10, 11, 12, 13, 14, 15, 16, 17 Despite a large body of data supporting a role for stress in MDD, the nature and degree of its involvement remains uncertain.18 Knowing whether stress-induced depression shares any or all molecular mechanisms with that of endogenous depression may guide the development of different treatment alternatives.

Animal models that reproduce key symptoms of MDD offer a unique opportunity for experimental exploration unfeasible in human studies.13, 14, 17, 19, 20, 21, 22 Selective breeding has been used to generate animal models of depression with extreme traits.23, 24, 25, 26, 27, 28 We used this approach to create two substrains of the Wistar Kyoto (WKY) rat, an accepted model of MDD that replicates many behavioral and physiological endophenotypes present in major depression.29, 30, 31, 32 The two substrains of the WKY, identified as WKY More Immobile (WMI) and WKY Less Immobile (WLI), were generated by bidirectional selective breeding from the WKY based on depressive behavior in the forced swim test (FST) of the naïve, unstressed animals; therefore, WMIs represent an endogenous depression model.33 The chronic stress model employed four phylogenetically, physiologically and behaviorally different rat strains34, 35 to identify general, strain-independent molecular characteristics of the chronic stress response. We chose a chronic restraint stress (CRS) paradigm that has been shown to increase depression-like behaviors36, 37, 38, 39, 40, 41, 42 and reduce hippocampal volume.43, 44, 45

The hippocampus and amygdala were selected to investigate gene expression profiles, as these brain regions have been strongly implicated in the cause and consequences of both depression and chronic stress.46, 47, 48, 49, 50, 51 Our goal was to examine both the monoamine and stress hypotheses of depression, the former by identifying the contribution of genes implicated by the monoamine hypothesis to the two depression model expression profiles, and the latter by defining the overlap between the expression profiles of endogenous depression and chronic stress response. Our present results found no evidence for the involvement of monoamine neurotransmission-related genes in the expression profile of the endogenous depression model and found divergent expression profiles in endogenous depression and chronic stress models. However, novel target genes and pathways have emerged from this study, which have the potential to advance our knowledge with regard to the etiology of depression and its treatment.

Methods

Animal models

All procedures were approved by the Institutional Animal Care and Use Committee of Northwestern University. The WMI-WLI selective breeding commenced as described previously.33 Animals showing the most extreme FST behavior within each line were selected for breeding, specifically avoiding sibling mating until the G5 generation. Adult WMI and WLI male animals from the 13th generation of selective breeding were employed in this study. The experimental design included daily administration of desipramine (10 mg kg−1) or saline, subcutaneously for 14 days, to both WMI and WLI animals (_n_=9 per treatment per strain). Animals were killed by fast decapitation on the 15th day immediately after removal from the home cage.

For the CRS experiment, adult male Fischer 344 (F344), Brown-Norway (BN), Lewis (LEW) and WKY rats were obtained from Harlan Laboratories (Indianapolis, IN, USA) (_n_=9 per strain per treatment). Rats were either exposed to CRS in a breathable decapicone for 2 h per day for a 2-week period or remained in their home cage (no stress, NS). Body weight was monitored throughout the experiment. On the 15th day, both groups of rats were tested in the elevated plus maze test in parallel and killed by decapitation immediately following the 5-min test. Animals were ∼100 days old at the time of sample collections. Blood samples were collected for the determination of plasma corticosterone (CORT) levels and adrenal weight was also determined. The elevated plus maze test of anxiety was used as a confirmation of strain differences in anxiety measures basally (NS group) and after CRS: the results are not shown since the purpose of this study is to identify chronic stress-induced changes in gene expression independently of strain effects.

Behavioral testing

FST was performed as described previously.30, 52 Briefly, on day 1, rats were individually placed in the water tank (water temperature 22–24 °C) for 15 min. After 24 h, rats were once again placed in the tank for a 5-min test session. The test session was videotaped and scored by a trained observer, using the scoring system developed by Detke et al.53 The open-field test was performed as described previously,54 but the animals' movements were followed and analyzed by the TSE Videomot 2 version 5.75 software.

Radioimmunoassay for plasma CORT

Assays were carried out in duplicate, as described previously,29 using the CORT RIA kit (MP Biomedicals, Solon, OH, USA). The assay sensitivity was 2–4 pg per tube. The intra- and interassay coefficients of variation were 3.5 and 8%, respectively.

Microarray experiments

Brain regions were dissected immediately after decapitation, as described previously,55 and stored at −80 °C in RNAlater (Ambion, Austin, TX, USA). Individual hippocampi and amygdalae were homogenized in TRIzol (Invitrogen, Carlsbad, CA, USA), and RNA was isolated following the manufacturer's protocol. All RNA samples were treated with DNase1 (Qiagen, Valencia, CA, USA) according to the manufacturer's methods.

Total RNA isolated from WMI and WLI brain regions was reverse-transcribed, and-double stranded cDNA was synthesized with the GeneChip® Expression 3-Amplification One-cycle kit (Affymetrix, Santa Clara, CA, USA. In an in vitro reaction with T7 RNA polymerase, the cDNA was linearly amplified and labeled with biotinylated nucleotides (Affymetrix). Ten micrograms of biotin-labeled and fragmented cRNA was then hybridized onto Rat Genome 230 2.0 GeneChip arrays (Affymetrix).

Total RNA from the chronic stress experiment was reverse transcribed, followed by second-strand cDNA synthesis. For each sample, an in vitro transcription reaction was carried out incorporating biotinylated nucleotides according to the manufacturer's protocol for Illumina® Totalprep RNA amplification kit (Ambion). Biotin-labeled cRNA, 1.5 μg, was then hybridized onto RatRef-12 Expression BeadChips (Illumina, San Diego, CA, USA).

Statistical analysis for WMI-WLI experiment

Probe intensity data from Rat 230v2 Affymetrix GeneChip arrays were read into the R software environment (http://www.R-project.org) directly from CEL files using the R/affy package.56 Affymetrix data was normalized using the robust multiarray average method for probe set data.57 Data quality was assessed using histograms of signal intensities, scatter plots and hierarchical clustering of samples.

Analysis of variance (ANOVA) methods, performed with the R/maanova package,58, 59 were used to statistically resolve gene expression differences. The experimental design included chronic administration of desipramine or saline to both WMI and WLI animals, of which three animals per strain per treatment were randomly selected for the microarray analyses. As we were only interested in the strain effect in this experiment, treatment (desipramine or saline) was a covariate. For each probe set g, strain i, covariate j and replicate k, a linear model for the log-transformed expression measure, y ijkg, can be formulated as a sum of components that contribute to the overall intensity value:

where μ g is the mean intensity over all 12 samples, α ig is the effect of strain i (_i_=1, 2), β jg is the additive effect of covariate j (_j_=1, 2), (αβ)ijg is the interaction effect of strain i, covariate j and probe set g and ɛ ijkg is the residual error for strain i, covariate j, replicate k and probe set g, respectively. To identify transcripts with a high probability of being differentially expressed between strains, statistical tests for the null hypothesis H0: α _ig_=0 vs the alternative H0: α _ig_≠0 were performed. We excluded transcripts in which the strain effect was not consistent between different covariate levels, as identified by a statistical test for the null hypothesis H0: (αβ)_ijg_=0 vs the alternative H0: (αβ)_ijg_≠0. Both these tests were performed using Fs, a modified F statistic incorporating shrinkage estimates of residual variance.60 _P_-values were calculated by permuting model residuals 1000 times. Unless otherwise noted, transcripts with differences between strains were identified as those with _P_-values, or estimated false-positive rates, less than 0.01 for the hypothesis H0: α _ig_=0 and greater than 0.01 for the hypothesis H0: (αβ)_ijg_=0. As an adjustment for multiple testing, we used the _q_-value transformation of the _P_-values to estimate false discovery rate.61 To relate probe sets to genes, probe set IDs were mapped to symbols using NetAffx (http://www.affymetrix.com/analysis/netaffx/). When reporting fold changes, we consider WLI to be the reference strain. If the measured expression level is higher in the WMI strain than in the WLI strain, the fold change is positive, and if it is lower in WMI than in WLI, the fold change is negative.

For each experiment, we also applied models of combined data from both tissues to measure the differences in transcript abundance between tissue. We used the model

which is equivalent to model (1), except that there is an additional variable (γ kg) for tissue with indices _k_=1,2 for the amygdala and hippocampus. Because we applied the normalization separately by tissue, we mean centered the normalized samples before applying this model.

Statistical analysis for stress experiment

Probe intensity data from Illumina Chip arrays were read directly into the R software environment from bead summary files produced by BeadStudio using the R/beadarray package.62 Quantile normalization was applied to the Illumina bead summary data using the R/preprocessCore package.63 Data quality was assessed as described for the WMI-WLI experiment.

Analysis of variance methods were similarly applied for this experiment. Twelve rats were exposed to chronic stress conditions, and 12 rats were not exposed to chronic stress. The animals for each stress condition consisted of three rats each from four different strains (BN-SS, F344, LEW and WKY). The linear model of equation (1) was applied for each probe set g. In this case, α ig, is the effect of condition i (_i_=1, 2), β jg is the additive effect of strain j (_j_=1, 2, 3, 4), which is a covariate. To identify transcripts exhibiting consistent differential expression owing to condition across strains, we performed and applied statistical tests for the null hypotheses H0: α _ig_=0 vs H0: α _ig_≠0 and H0: (αβ)_ijg_=0 vs H0: (αβ)_ijg_≠0 that are analogous to those described above for the WLI-WMI comparison. Probe set IDs were mapped to symbols using Illumina's annotation resource (http://www.illumina.com/support/annotation_files.ilmn). When reporting fold changes, we consider the NS to be the reference condition. If the measured expression in the CRS condition is greater than in the NS condition, the fold change is positive. If it is less under CRS than under NS, the fold change is negative. To facilitate across-experiment comparison by gene, probes were mapped to gene symbols using current probe annotation files provided by the array manufacturers. The suffixes ‘_predicted’ and ‘_mapped’ were truncated from the Illumina symbol assignments. Where there were multiple probes for the same symbol, probes with the largest F statistic for the relevant hypothesis test were chosen to give a one-to-one mapping across experiments for 10 112 genes. The linear model of equation (2) was applied analogously for this experiment to obtain the distribution of estimated tissue effects.

Real-time RT-PCR

Real-time reverse transcription-polymerase chain reaction (RT-PCR) was used to confirm microarray results for a subset of randomly selected genes from those fulfilling the criteria of significant (P<0.01) expression differences. Reverse transcription for 2 μg of each sample was performed with Invitrogen's Superscript® III First-Strand kit (18080-051) according to the manufacturer's protocol. Primers were designed to amplify 80–150 bp regions and to contain a maximum amount of microarray probe sequence using default settings of ABI's Primer Express software (version 3.0, PE Applied Biosystems, Carlsbad, CA, USA). Primer sequences are listed in Supplementary Table S1. Forty nanograms of cDNA was amplified in 20 μl reactions (1 × SYBR green reaction mix (ABI, Carlsbad, CA, USA), 250 μM primers) in the ABI 7900HT PCR machine using the relative quantification (–ddCt) method, with 18s RNA as the internal control.

Results

Selective breeding and characterization of WMI and WLI strains

The WMI and WLI rat strains consistently maintained significant, dichotomous FST phenotypes throughout 21 generations of selective breeding (Figure 1a), with WMIs always exhibiting greater immobility scores than WLIs. In the open-field test of exploration/anxiety, WMI and WLI male animals showed similar level of exploration of the inner circle (Figure 1b), suggesting no differences in anxiety-related behaviors. In contrast, WMIs explored the arena significantly less than WLIs (Figure 1c), and the activity traces (Figure 1d) illustrate freezing-like behavior exhibited by the WMIs, a behavioral pattern very similarly to psychomotor retardation. There were no differences in the two strains' basal plasma thyroxine and CORT levels or their CORT responses to acute restraint stress (Supplementary Table S2). Furthermore, no additional behavioral differences were observed between the two strains, including anxiety-like behaviors as assessed by the elevated plus maze and defensive burying tests, confirming the results of the open-field test, and spatial learning and memory as measured in the Morris water maze (Supplementary Table S2). These data verify that the behavioral differences between the two substrains of WKYs are not fear or anxiety driven, but rather related to depressive state.

Figure 1

The endogenous depression model, the Wistar–Kyoto More Immobile (WMI) strain, shows depressive behavior not linked to fear/anxiety (a–d). Chronic restraint stress (CRS) increases adrenocortical function consistently in all four strains (e and f). (a) In the forced swim test (FST), immobility scores of the WMI and WKY Less Immobile (WLI) animals differ significantly across generations. (b) Time spent in the inner circle of the open-field test, (c) total distance traveled and (d) movement traces of representative WMI and WLI animals. (e) Plasma corticosterone levels were consistently elevated after CRS in all four strains and (f) adrenal weights were consistently greater in the CRS group of all strains. *P<0.01, **P<0.001.

PowerPoint slide

Full size image

CRS decreased body weight gain, increased plasma CORT levels and adrenal weight

Genetic polymorphisms between the strains range from 25.9% between F344 and LEW to 66% between BN and WKY, representing a substantial interstrain variation. However, CRS affected all strain of rats. Specifically, CRS resulted in a lower body weight gain compared with NS rats (F[1,65]=212.9, P<0.001), but there was no weight loss in response to the CRS paradigm (body weight on day 1: 239±5 g; on day 14 of CRS: 244±5 g). Plasma CORT levels were elevated in all rats exposed to the chronic stress procedure compared with those not subjected to stress (Figure 1e; strain: (F[1,65]=105.37, P<0.001) and those subjected to stress (F[3,65]=6.75, P<0.001; strain, stress interaction NS). As expected, chronic stress increased adrenal weight in all four rat strains (Figure 1f; strain: F[3,41]=1.34, NS; stress: F[1,41]=4.77, _P_=0.035). Taken together, these data established that all CRS rats were undergoing expected physiological changes associated with CRS exposure.

Differential gene expression profile in the amygdala and hippocampus of WMI and WLI male rats

Six hundred and thirty-eight genes in the amygdala and 463 in the hippocampus were differentially expressed between the WMI and WLI animals (P<0.01). Genes that additionally had fold changes above 1.4 (40% increase or decrease) are listed in Table 1. The complete data set is shown in Supplementary Tables S3 (amygdala) and S4 (hippocampus). Twenty-seven genes were differentially expressed in both brain regions between WMI and WLI in the same direction.

Table 1 Differentially expressed depression genes (top) and stress genes (bottom) in the amygdala and hippocampus

Full size table

Gene expression changes in response to chronic stress in the amygdala and hippocampus of four strains of male rats

To eliminate strain-specific effects in our selection attributable to differences in stress reactivity or potential cognitive differences between strains, genes were selected based on consistent, significant expression level differences in all four strains. We found 125 genes in the amygdala and 126 genes in the hippocampus that were differentially expressed between the CRS and NS conditions (P<0.01) (Supplementary Tables S5 and S6). From these, the 10 amygdalar and hippocampal genes with absolute fold changes above 1.4 are shown in Table 1. Twelve of the genes showed the same directional changes in both brain regions.

Minimal overlap between depression and stress response

Both the WMI-WLI and CRS-NS gene sets were examined for brain region-specific overlap using all known aliases to account for possible differences between microarray platform target nomenclatures. Only two genes in the amygdala (Kcnj14 and Mprl9) and four in the hippocampus (Chi3l1, collagen type I alpha 1 (Col1a1), matrix metallopeptidase 14 (Mmp14) and RGD1305680) showed differential expression in both the chronic stress response and endogenous depression comparisons. On the basis of the stress hypothesis of depression, transcripts with lowered expression levels in WMI animals should show similarly decreased expression in CRS animals, and increased expression in WMIs should correspond to increased expression in CRS. When we examined the directionality of fold changes for the six overlapping genes, only three genes (Col1a1, Mmp14 and RGD1305680), all in hippocampal tissue, showed these hypothesized directional changes. This provides suggestive evidence that the molecular signature of chronic stress and endogenous depressive behavior in these brain regions are not similar. However, since comparison of overlap between small sets of selected genes can be sensitive to thresholding, we also compared the differences between model effects for CRS and NS with the differences for WMI and WLI.

Expression differences in both models were relatively small (Figures 2a and b). If there is an underlying relationship between expression differences in response to CRS and endogenous depressive behavior for many genes, we should see significant positive correlation between the differences in gene expression across the two experiments. In particular, the differences between WMI and WLI model effects and differences between CRS and NS model effects ((α_2g−_α_1_g) in equation (1)) should show positive correlation. The Spearman's correlation (ρ) for these comparisons of differences were not significantly different than zero in either tissue (_ρ_=−0.12, amygdala; _ρ_=−0.02, hippocampus; Figures 2a and b). This lack of observed correlation could be explained by the different platforms or by the chronologies of the experiments. To address this concern, equation (2), which combined data of both brain regions and estimates the tissue effects, was employed. The tissue (interbrain region) effects showed significant positive correlation between experiments (_ρ_=0.48, P<0.001; Figure 2c). We also measured correlation for the subset of genes for which the magnitude of tissue effects were comparable to the CRS-NS and WMI-WLI effects, specifically those having a magnitude below 0.4 on both platforms. Within this subset, we still found significant positive correlation (_ρ_=0.31, P<0.001). We conclude that we should have observed correlation between the CRS-NS and WMI-WLI effects if it were present despite the platform differences. This suggests that, on a global scale, the effects of chronic stress and endogenous depression on the amygdalar and hippocampal transcriptome are dissimilar.

Figure 2

An across-experiment, gene-to-gene mapping was created for 10 112 gene symbols. Scatter plots are shown for observed effects obtained from the chronic restraint stress-no stress (CRS-NS) microarray experiment (vertical axis) and the Wistar–Kyoto More Immobile-WKY Less Immobile (WMI-WLI) microarray experiment (horizontal axis) by applying the linear models of Equation (1) and (2). The statistics are differences in condition effects, α_2_g_−_α_1_g, for the amygdala (a) and hippocampus (b), and differences in brain region-related effects, γ_2_g_−_γ_1_g (c).

PowerPoint slide

Full size image

Confirmation of the expression differences by real-time RT-PCR

We carried out quantitative RT-PCR confirmation of the microarray data using total RNA samples from both experiments and primer pairs described in Supplementary Table S1. Pearson's correlation of fold change from the Affymetrix microarray for 16 genes with relative quantification ratios from qPCR (RQ ratio) showed significant concurrence of data (_r_=0.720, _P_=0.002; Figure 3a). Pearson's correlation of fold change from the Illumina microarray for nine genes (three were confirmed in both brain regions) with relative quantification ratios from qPCR showed similarly significant correlation (_r_=0.725, _P_=0.008; Figure 3b). Positive and negative fold changes are defined in Methods.

Figure 3

(a) Validation of genes differentially expressed in Wistar–Kyoto More Immobile (WMI) and WKY Less Immobile (WLI) hippocampi or amygdala (_n_=6 per strain) by real-time reverse transcription-polymerase chain reaction (RT-PCR). The correlation between fold change in the Affymetrix microarray experiment and real-time RT-PCR determination of relative quantification ratios are shown (Pearson's correlation, _r_=0.720, _P_=0.002). (b) Validation of genes differentially expressed in chronic restraint stress (CRS) vs no stress (NS) (_n_=12 per treatment) hippocampi or amygdala by real-time RT-PCR. The correlation between fold change in the Illumina array experiment and real-time RT-PCR determination of relative quantification ratios are shown (Pearson's correlation, _r_=0.725, _P_=0.008). (c) Lack of correlation between absolute fold change and significant _P-_value (–log P) from the microarray analyses for genes with expression changes validated by real-time RT-PCR.

PowerPoint slide

Full size image

When we examined the validated transcript fold change in the microarrays in relation to their respective _P_-value, we found no significant correlation (_r_=0.113, _P_>0.05; Figure 3c). These data confirm that significance and fold change criteria cannot be used interchangeably, but that expression changes can be confirmed by real-time RT-PCR, even if the fold change is small.

The molecular pathways of depression and stress response genes

The lists of hippocampal and amygdalar genes with significant expression differences were combined into single gene lists for each of our two models, and the rat gene symbols were converted to their human orthologs to access a larger body of literature. These master ortholog lists were subjected to functional annotation using the PANTHER Classification System.64 In addition, we used the PANTHER binomial statistics tool to compare our annotated ortholog lists to the NCBI Homo sapiens reference gene list to determine significant over- or under-representation of molecular function, biological process or pathway classification terms.65

Table 2 shows significantly over- and under-represented terms from both gene sets. A group of genes identified to belong to the integrin signaling pathway (including Araf, Arf1, Arf4, Arhgap10, Col12a1, Col15a1, Col1a1, Col4a1, Col5a2, Col6a6, Crkl, Elmo1, Fnlb, Grb2, Itga11, Itga6, Itgae, Itgb6, Itgb8, Itgbl1, Mapk1, Pik3r1, Pitk2b, Rac1, Rhoa and Tln2) were most overly represented in the WMI-WLI comparison. From the CRS-NS comparison, the Huntington disease pathway showed the only significantly clustered group of genes (Capn9, Gapdh, Rhog, Tnfaip8, Tnfaip8l3 and tumor suppressor protein 53).

Table 2 Over- and under-represented gene ontology (GO) groups determined by the Panther software analysis

Full size table

The biological process terms ‘Cellular Process’, ‘Cell Communication’ and ‘Signal Transduction’ were significantly over-represented in both the CRS-NS and the WMI-WLI gene lists. It is of further interest that metabolic process terms are the most over-represented in the WMI-WLI comparison, while cellular process terms, mainly immune system-related, are over-represented in the CRS-NS gene list. There was no overlap in over- or under-represented molecular function terms between the stress response and endogenous depression gene lists.

Discussion

This study identified novel molecular signatures of chronic stress and of endogenous depressive behavior in animal models. However, absence of direct association between the gene expression profiles suggests that independent molecular mechanisms regulate these two states.

Among genes showing significantly altered transcription levels in the amygdala and hippocampus in the endogenous depression model, there was a complete omission of serotoninergic, adrenergic and dopaminergic neurotransmission-related genes. This finding, although surprising, is congruent with gene expression studies of human postmortem MDD brains.66, 67, 68 We did find that multiple newer antidepressant targets, including Grm5,69 and a number of phosphodiesterases,70 were differentially expressed between WMIs and WLIs. In addition to these previously described targets, we identified several novel depressive behavior-related genes and pathways. Physiologically intriguing findings include the decreased mRNA levels of prolactin receptor and angiotensin I-converting enzyme in the WMI hippocampus; these genes have known peripheral function, but they have not been associated with depressive behavior as yet. Other intriguing examples of novel implications for endogenous depression include the solute carrier family 4, sodium bicarbonate co-transporter and member 5 gene (Slc4a5), which has recently been associated with metabolic phenotypes.71 Interestingly, metabolic processes were the most over-represented biological process terms in the WMI-WLI gene ontology. The integrin signaling pathway, which was over-represented in the endogenous depression gene set, contains different members of the collagen family, some of which are also implicated in the chronic stress response. Finally, mRNA levels of dehydrodolichyl diphosphate synthase are increased in both the hippocampus and the amygdala of WMI. Dehydrodolichyl diphosphate synthase is responsible for synthesizing dolichol, which accumulates in the neuropathological human brain.72

The set of genes whose hippocampal or amygdalar expression patterns were altered by chronic stress in all four rat strains represent a generalizable molecular response to chronic stress. Many of these genes have shown matching directional changes by unpredictable chronic mild stress, and these include: Col1a1, tissue plasminogen activator, insulin-like growth factor binding protein 2, amyloid beta precursor protein and transthyretin in the amygdala and/or hippocampus.39, 73 In addition, stress has previously been found to increase the expression of amyloid beta precursor protein74, 75 and neuropeptide Y.76, 77, 78 Our findings, however, also implicate a number of novel ‘stress genes’. Among these, the decreased mRNA levels of tumor suppressor protein 53, insulin-like growth factor 2 and hemoglobin beta and alpha a1 are of particular interest. Tumor suppressor protein 53 has been linked to the pathology of Huntington disease,79 a finding that is connected to the over-representation of our chronic stress response genes involved in the Huntington pathway. Insulin-like growth factor 2 is an imprinted gene, described as having both neuroprotective and neurodegenerative properties,80 as well as involvement in metabolic disorders. Recently, hemoglobin beta and alpha transcripts and proteins have been found in cortical and hippocampal astrocytes, mature oligodendrocytes and a subpopulation of dopaminergic neurons.81 If the decreased expression of hemoglobin beta and alpha a1 is a manifestation of chronic stress-induced decreases in oxygen storage, these findings would implicate chronic stress in neurodegenerative processes. If, however, they are related to the neurotoxicity of hemoglobins,82 the decreased expression of hemoglobin beta and alpha a1 and of tumor suppressor protein 53 would propose the heretical idea of chronic stress being neuroprotective.

If depression in general is a result of chronic stress and stress-related factors, these genes should be shared with the endogenous ‘depression transcriptome’, as defined by WMI-WLI differences. However, a surprisingly small number of genes were shared between WMI-WLI and NS-CRS differentially expressed genes in the amygdala and hippocampus. Genes showing parallel directional expression changes in the hippocampus between the ‘endogenous depression’ and the ‘chronic stress’ gene sets include: Col1a1, Mmp14 and RGD1305680 (homologous to the KIAA0240 gene). Mmp14 (also known as MT1-MMP) degrades endogenous beta amyloids83 and is implicated in the pathogenesis of central nervous system inflammatory disorders. It is also a tethered membrane collagenase, which, together with the involvement of Col1a1, implies the significance of extracellular matrix remodeling in depressive behavior and chronic stress. Ontological analysis further and more broadly confirmed the dissimilarity between the endogenous depressive behavior and generic stress response expression profiles. Functional annotation of the genes involved in both states revealed few common biological processes and no common molecular functions between the two gene sets.

This study using animal models contends that endogenous depression and chronic stress response are regulated by independent molecular pathways. In humans, an inference of this finding would be that stress-induced and endogenous subtypes of depression are etiologically distinct. There are both human and animal studies supporting this assumption. For example, only a subset of depressive phenotypes can be causally attributed to stressful life events, and these phenotypes indicate only a moderately increased risk of developing depression.84, 85 In addition, although the chronic mild stress model is thought of as a model of depression, some animals are resistant to the anhedonic effect of stress.73, 86, 87, 88 Similarly, some animals are resilient to learned helplessness in the congenital learned helplessness model of depression,25, 89 suggesting that differences in vulnerability to stress may determine the development of depressive phenotype in these models. Thus, it is possible that stress-susceptible individuals could develop a subtype of depression distinct from endogenous depression; this subtype may show greater comorbidity with other stress-related disorders. The existence of different subtypes of depression is in further agreement with the findings of Krishnan and Nestler,84 which state that there is no ‘unified theory’ of depression. Should our findings translate to MDD, they suggest that different treatment strategies, dependent on the patient's depressive subtype, could be beneficial. Furthermore, as most animal models of depression are based on stress, the above findings predict that using various genetic animal models of depression could lead to novel drug targets that may achieve successful treatment of patients in which current treatment methods fail.

References

1. Rumsfeld JS, Ho PM . Depression and cardiovascular disease: a call for recognition. Circulation 2005; 111: 250–253.
   Article PubMed Google Scholar
2. Swaab DF, Bao AM, Lucassen PJ . The stress system in the human brain in depression and neurodegeneration. Ageing Res Rev 2005; 4: 141–194.
   Article CAS PubMed Google Scholar
3. Bornstein SR, Schuppenies A, Wong ML, Licinio J . Approaching the shared biology of obesity and depression: the stress axis as the locus of gene–environment interactions. Mol Psychiatry 2006; 11: 892–902.
   Article CAS PubMed Google Scholar
4. Kessler RC, Ormel J, Demler O, Stang PE . Comorbid mental disorders account for the role impairment of commonly occurring chronic physical disorders: results from the National Comorbidity Survey. J Occup Environ Med 2003; 45: 1257–1266.
   Article PubMed Google Scholar
5. Altar CA, Vawter MP, Ginsberg SD . Target identification for CNS diseases by transcriptional profiling. Neuropsychopharmacology 2009; 34: 18–54.
   Article CAS PubMed Google Scholar
6. Delgado PL . Depression: the case for a monoamine deficiency. J Clin Psychiatry 2000; 61 (Suppl 6): 7–11.
   CAS PubMed Google Scholar
7. Ruhe HG, Mason NS, Schene AH . Mood is indirectly related to serotonin, norepinephrine and dopamine levels in humans: a meta-analysis of monoamine depletion studies. Mol Psychiatry 2007; 12: 331–359.
   Article CAS PubMed Google Scholar
8. Fava M, Kendler KS . Major depressive disorder. Neuron 2000; 28: 335–341.
   Article CAS PubMed Google Scholar
9. Holsboer F . The stress hormone system is back on the map. Curr Psychiatry Rep 2000; 2: 454–456.
   Article CAS PubMed Google Scholar
10. Wong ML, Licinio J . Research and treatment approaches to depression. Nat Rev Neurosci 2001; 2: 343–351.
    Article CAS PubMed Google Scholar
11. Gold PW, Chrousos GP . Organization of the stress system and its dysregulation in melancholic and atypical depression: high vs low CRH/NE states. Mol Psychiatry 2002; 7: 254–275.
    Article CAS PubMed Google Scholar
12. Nestler EJ, Barrot M, DiLeone RJ, Eisch AJ, Gold SJ, Monteggia LM . Neurobiology of depression. Neuron 2002; 34: 13–25.
    Article CAS PubMed Google Scholar
13. Hasler G, Drevets WC, Manji HK, Charney DS . Discovering endophenotypes for major depression. Neuropsychopharmacology 2004; 29: 1765–1781.
    Article CAS PubMed Google Scholar
14. Cryan JF, Holmes A . The ascent of mouse: advances in modelling human depression and anxiety. Nat Rev Drug Discov 2005; 4: 775–790.
    Article CAS PubMed Google Scholar
15. de Kloet ER, Joels M, Holsboer F . Stress and the brain: from adaptation to disease. Nat Rev Neurosci 2005; 6: 463–475.
    Article CAS PubMed Google Scholar
16. Bale TL . Stress sensitivity and the development of affective disorders. Horm Behav 2006; 50: 529–533.
    Article CAS PubMed Google Scholar
17. Muller MB, Holsboer F . Mice with mutations in the HPA-system as models for symptoms of depression. Biol Psychiatry 2006; 59: 1104–1115.
    Article CAS PubMed Google Scholar
18. Monroe SM, Reid MW . Gene–environment interactions in depression research: genetic polymorphisms and life-stress polyprocedures. Psychol Sci 2008; 19: 947–956.
    Article PubMed Google Scholar
19. Jacobson LH, Cryan JF . Feeling strained? Influence of genetic background on depression-related behavior in mice: a review. Behav Genet 2007; 37: 171–213.
    Article CAS PubMed Google Scholar
20. Kas MJ, Fernandes C, Schalkwyk LC, Collier DA . Genetics of behavioural domains across the neuropsychiatric spectrum; of mice and men. Mol Psychiatry 2007; 12: 324–330.
    Article CAS PubMed Google Scholar
21. Urani A, Chourbaji S, Gass P . Mutant mouse models of depression: candidate genes and current mouse lines. Neurosci Biobehav Rev 2005; 29: 805–828.
    Article CAS PubMed Google Scholar
22. Tecott LH . The genes and brains of mice and men. Am J Psychiatry 2003; 160: 646–656.
    Article PubMed Google Scholar
23. Scott PA, Cierpal MA, Kilts CD, Weiss JM . Susceptibility and resistance of rats to stress-induced decreases in swim-test activity: a selective breeding study. Brain Res 1996; 725: 217–230.
    Article CAS PubMed Google Scholar
24. Steimer T, Driscoll P . Divergent stress responses and coping styles in psychogenetically selected Roman high-(RHA) and low-(RLA) avoidance rats: behavioural, neuroendocrine and developmental aspects. Stress 2003; 6: 87–100.
    Article CAS PubMed Google Scholar
25. Henn FA, Vollmayr B . Stress models of depression: forming genetically vulnerable strains. Neurosci Biobehav Rev 2005; 29: 799–804.
    Article PubMed Google Scholar
26. Weiss JM, Cierpial MA, West CHK . Selective breeding of rats for high and low motor activity in a swim test: toward a new animal model of depression. Pharmacol Biochem Behav 1998; 61: 49–66.
    Article CAS PubMed Google Scholar
27. Weiss JM, West CH, Emery MS, Bonsall RW, Moore JP, Boss-Williams KA . Rats selectively-bred for behavior related to affective disorders: proclivity for intake of alcohol and drugs of abuse, and measures of brain monoamines. Biochem Pharmacol 2008; 75: 134–159.
    Article CAS PubMed Google Scholar
28. El Yacoubi M, Bouali S, Popa D, Naudon L, Leroux-Nicollet I, Hamon M et al. Behavioral, neurochemical, and electrophysiological characterization of a genetic mouse model of depression. Proc Natl Acad Sci USA 2003; 100: 6227–6232.
    Article CAS PubMed PubMed Central Google Scholar
29. Redei E, Pare WP, Aird F, Kluczynski J . Strain differences in hypothalamic-pituitary-adrenal activity and stress ulcer. Am J Physiol 1994; 266 (Part 2): R353–R360.
    CAS PubMed Google Scholar
30. Redei EE, Ahmadiyeh N, Baum A, Sasso D, Slone J, Solberg LC et al. Novel animal models of affective disorders. In: Tucker GJ, McKinney W (eds). Seminars in Clinical Neuropsychiatry, vol. 6. WB Saunders Company: Philadelphia, 2001.
    Google Scholar
31. Solberg LC, Ahmadiyeh N, Baum AE, Vitaterna MH, Takahashi JS, Turek FW et al. Depressive-like behavior and stress reactivity are independent traits in a Wistar Kyoto × Fisher 344 cross. Mol Psychiatry 2003; 8: 423–433.
    Article CAS PubMed Google Scholar
32. Rittenhouse PA, Lopez-Rubalcava C, Stanwood GD, Lucki I . Amplified behavioral and endocrine responses to forced swim stress in the Wistar–Kyoto rat. Psychoneuroendocrinology 2002; 27: 303–318.
    Article PubMed Google Scholar
33. Will CC, Aird F, Redei EE . Selectively bred Wistar–Kyoto rats: an animal model of depression and hyper-responsiveness to antidepressants. Mol Psychiatry 2003; 8: 925–932.
    Article CAS PubMed Google Scholar
34. Thomas MA, Chen CF, Jensen-Seaman MI, Tonellato PJ, Twigger SN . Phylogenetics of rat inbred strains. Mamm Genome 2003; 14: 61–64.
    Article PubMed Google Scholar
35. Gomez F, Lahmame A, de Kloet ER, Armario A . Hypothalamic-pituitary-adrenal response to chronic stress in five inbred rat strains: differential responses are mainly located at the adrenocortical level. Neuroendocrinology 1996; 63: 327–337.
    Article CAS PubMed Google Scholar
36. Wood GE, Norris EH, Waters E, Stoldt JT, McEwen BS . Chronic immobilization stress alters aspects of emotionality and associative learning in the rat. Behav Neurosci 2008; 122: 282–292.
    Article PubMed Google Scholar
37. Bravo JA, Diaz-Veliz G, Mora S, Ulloa JL, Berthoud VM, Morales P et al. Desipramine prevents stress-induced changes in depressive-like behavior and hippocampal markers of neuroprotection. Behav Pharmacol 2009; 20: 273–285.
    Article PubMed Google Scholar
38. Gil M, Marti J, Armario A . Inhibition of catecholamine synthesis depresses behavior of rats in the holeboard and forced swim tests: influence of previous chronic stress. Pharmacol Biochem Behav 1992; 43: 597–601.
    Article CAS PubMed Google Scholar
39. Joo Y, Choi KM, Lee YH, Kim G, Lee DH, Roh GS et al. Chronic immobilization stress induces anxiety- and depression-like behaviors and decreases transthyretin in the mouse cortex. Neurosci Lett 2009; 461: 121–125.
    Article CAS PubMed Google Scholar
40. Kim KS, Han PL . Optimization of chronic stress paradigms using anxiety- and depression-like behavioral parameters. J Neurosci Res 2006; 83: 497–507.
    Article CAS PubMed Google Scholar
41. Haenisch B, Bilkei-Gorzo A, Caron MG, Bonisch H . Knockout of the norepinephrine transporter and pharmacologically diverse antidepressants prevent behavioral and brain neurotrophin alterations in two chronic stress models of depression. J Neurochem 2009; 111: 403–416.
    Article CAS PubMed PubMed Central Google Scholar
42. Veena J, Srikumar BN, Raju TR, Shankaranarayana Rao BS . Exposure to enriched environment restores the survival and differentiation of new born cells in the hippocampus and ameliorates depressive symptoms in chronically stressed rats. Neurosci Lett 2009; 455: 178–182.
    Article CAS PubMed Google Scholar
43. Lee T, Jarome T, Li SJ, Kim JJ, Helmstetter FJ . Chronic stress selectively reduces hippocampal volume in rats: a longitudinal magnetic resonance imaging study. Neuroreport 2009; 20: 1554–1558.
    Article PubMed PubMed Central Google Scholar
44. Vyas A, Pillai AG, Chattarji S . Recovery after chronic stress fails to reverse amygdaloid neuronal hypertrophy and enhanced anxiety-like behavior. Neuroscience 2004; 128: 667–673.
    Article CAS PubMed Google Scholar
45. Vyas A, Bernal S, Chattarji S . Effects of chronic stress on dendritic arborization in the central and extended amygdala. Brain Res 2003; 965: 290–294.
    Article CAS PubMed Google Scholar
46. Bremner JD, Narayan M, Anderson ER, Staib LH, Miller HL, Charney DS . Hippocampal volume reduction in major depression. Am J Psychiatry 2000; 157: 115–118.
    Article CAS PubMed Google Scholar
47. Sheline YI, Wang PW, Gado MH, Csernansky JG, Vannier MW . Hippocampal atrophy in recurrent major depression. Proc Natl Acad Sci USA 1996; 93: 3908–3913.
    Article CAS PubMed PubMed Central Google Scholar
48. Sheline YI, Gado MH, Price JL . Amygdala core nuclei volumes are decreased in recurrent major depression. NeuroReport 1998; 9: 2023–2028.
    Article CAS PubMed Google Scholar
49. Sheline YI, Sanghavi M, Mintun MA, Gado MH . Depression duration but not age predicts hippocampal volume loss in medically healthy women with recurrent major depression. J Neurosci 1999; 19: 5034–5043.
    Article CAS PubMed PubMed Central Google Scholar
50. Hastings RS, Parsey RV, Oquendo MA, Arango V, Mann JJ . Volumetric analysis of the prefrontal cortex, amygdala, and hippocampus in major depression. Neuropsychopharmacology 2004; 29: 952–959.
    Article PubMed Google Scholar
51. Kim JS, Shukla SD . Acute in vivo effect of ethanol (binge drinking) on histone H3 modifications in rat tissues. Alcohol Alcohol 2006; 41: 126–132.
    Article CAS PubMed Google Scholar
52. Solberg LC, Baum AE, Ahmadiyeh N, Shimomura K, Li R, Turek FW et al. Sex- and lineage-specific inheritance of depression-like behavior in the rat. Mamm Genome 2004; 15: 648–662.
    Article PubMed PubMed Central Google Scholar
53. Detke MJ, Rickels M, Lucki I . Active behaviors in the rat forced swimming test differentially produced by serotonergic and noradrenergic antidepressants. Psychopharmacology 1995; 121: 66–72.
    Article CAS PubMed Google Scholar
54. Nosek K, Dennis K, Andrus BM, Ahmadiyeh N, Baum AE, Woods LC et al. Context and strain-dependent behavioral response to stress. Behav Brain Funct 2008; 4: 23.
    Article PubMed PubMed Central Google Scholar
55. Wilcoxon JS, Kuo AG, Disterhoft JF, Redei EE . Behavioral deficits associated with fetal alcohol exposure are reversed by prenatal thyroid hormone treatment: a role for maternal thyroid hormone deficiency in FAE. Mol Psychiatry 2005; 10: 961–971.
    Article CAS PubMed Google Scholar
56. Gautier L, Cope L, Bolstad BM, Irizarry RA . Affy—analysis of Affymetrix GeneChip data at the probe level. Bioinformatics 2004; 20: 307–315.
    Article CAS PubMed Google Scholar
57. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U et al. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 2003; 4: 249–264.
    Article PubMed Google Scholar
58. Churchill GA . Using ANOVA to analyze microarray data. Biotechniques 2004; 37: 173–175, 177.
    Article CAS PubMed Google Scholar
59. Wu X, Pang ST, Sahlin L, Blanck A, Norstedt G, Flores-Morales A . Gene expression profiling of the effects of castration and estrogen treatment in the rat uterus. Biol Reprod 2003; 69: 1308–1317.
    Article CAS PubMed Google Scholar
60. Cui X, Hwang JT, Qiu J, Blades NJ, Churchill GA . Improved statistical tests for differential gene expression by shrinking variance components estimates. Biostatistics 2005; 6: 59–75.
    Article PubMed Google Scholar
61. Storey JD, Xiao W, Leek JT, Tompkins RG, Davis RW . Significance analysis of time course microarray experiments. Proc Natl Acad Sci USA 2005; 102: 12837–12842.
    Article CAS PubMed PubMed Central Google Scholar
62. Dunning MJ, Smith ML, Ritchie ME, Tavare S . beadarray: R classes and methods for Illumina bead-based data. Bioinformatics 2007; 23: 2183–2184.
    Article CAS PubMed Google Scholar
63. Bolstad BM, Irizarry RA, Astrand M, Speed TP . A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 2003; 19: 185–193.
    Article CAS PubMed Google Scholar
64. Thomas PD, Kejariwal A, Campbell MJ, Mi H, Diemer K, Guo N et al. PANTHER: a browseable database of gene products organized by biological function, using curated protein family and subfamily classification. Nucleic Acids Res 2003; 31: 334–341.
    Article CAS PubMed PubMed Central Google Scholar
65. Cho RJ, Campbell MJ . Transcription, genomes, function. Trends Genet 2000; 16: 409–415.
    Article CAS PubMed Google Scholar
66. Sibille E, Arango V, Galfalvy HC, Pavlidis P, Erraji-Benchekroun L, Ellis SP et al. Gene expression profiling of depression and suicide in human prefrontal cortex. Neuropsychopharmacology 2004; 29: 351–361.
    Article CAS PubMed Google Scholar
67. Sequeira A, Mamdani F, Ernst C, Vawter MP, Bunney WE, Lebel V et al. Global brain gene expression analysis links glutamatergic and GABAergic alterations to suicide and major depression. PLoS One 2009; 4: e6585.
    Article CAS PubMed PubMed Central Google Scholar
68. Sequeira A, Turecki G . Genome wide gene expression studies in mood disorders. OMICS 2006; 10: 444–454.
    Article CAS PubMed Google Scholar
69. Covington III HE, Vialou V, Nestler EJ . From synapse to nucleus: novel targets for treating depression. Neuropharmacology 2010; 58: 683–693.
    Article CAS PubMed Google Scholar
70. Esposito K, Reierson GW, Luo HR, Wu GS, Licinio J, Wong ML . Phosphodiesterase genes and antidepressant treatment response: a review. Ann Med 2009; 41: 177–185.
    Article CAS PubMed Google Scholar
71. Stutz AM, Teran-Garcia M, Rao DC, Rice T, Bouchard C, Rankinen T . Functional identification of the promoter of SLC4A5, a gene associated with cardiovascular and metabolic phenotypes in the HERITAGE Family Study. Eur J Hum Genet 2009; 17: 1481–1489.
    Article CAS PubMed PubMed Central Google Scholar
72. Endo S, Zhang YW, Takahashi S, Koyama T . Identification of human dehydrodolichyl diphosphate synthase gene. Biochim Biophys Acta 2003; 1625: 291–295.
    Article CAS PubMed Google Scholar
73. Bergstrom A, Jayatissa MN, Thykjaer T, Wiborg O . Molecular pathways associated with stress resilience and drug resistance in the chronic mild stress rat model of depression: a gene expression study. J Mol Neurosci 2007; 33: 201–215.
    Article CAS PubMed Google Scholar
74. Dong H, Csernansky JG . Effects of stress and stress hormones on amyloid-beta protein and plaque deposition. J Alzheimers Dis 2009; 18: 459–469.
    Article CAS PubMed PubMed Central Google Scholar
75. Lee KH, Yu DH, Lee YS . Gene expression profiling of rat cerebral cortex development using cDNA microarrays. Neurochem Res 2009; 34: 1030–1038.
    Article CAS PubMed Google Scholar
76. Zambello E, Fuchs E, Abumaria N, Rygula R, Domenici E, Caberlotto L . Chronic psychosocial stress alters NPY system: different effects in rat and tree shrew. Prog Neuropsychopharmacol Biol Psychiatry 2010; 34: 122–130.
    Article CAS PubMed Google Scholar
77. Thorsell A, Slawecki CJ, El Khoury A, Mathe AA, Ehlers CL . The effects of social isolation on neuropeptide Y levels, exploratory and anxiety-related behaviors in rats. Pharmacol Biochem Behav 2006; 83: 28–34.
    Article CAS PubMed Google Scholar
78. Heilig M, Widerlov E . Neurobiology and clinical aspects of neuropeptide Y. Crit Rev Neurobiol 1995; 9: 115–136.
    CAS PubMed Google Scholar
79. Feng Z, Jin S, Zupnick A, Hoh J, de Stanchina E, Lowe S et al. p53 tumor suppressor protein regulates the levels of huntingtin gene expression. Oncogene 2006; 25: 1–7.
    Article CAS PubMed Google Scholar
80. Dikkes P, D BJ, Guo WH, Chao C, Hemond P, Yoon K et al. IGF2 knockout mice are resistant to kainic acid-induced seizures and neurodegeneration. Brain Res 2007; 1175: 85–95.
    Article CAS PubMed Google Scholar
81. Biagioli M, Pinto M, Cesselli D, Zaninello M, Lazarevic D, Roncaglia P et al. Unexpected expression of alpha- and beta-globin in mesencephalic dopaminergic neurons and glial cells. Proc Natl Acad Sci USA 2009; 106: 15454–15459.
    Article CAS PubMed PubMed Central Google Scholar
82. Wang X, Mori T, Sumii T, Lo EH . Hemoglobin-induced cytotoxicity in rat cerebral cortical neurons: caspase activation and oxidative stress. Stroke 2002; 33: 1882–1888.
    Article CAS PubMed Google Scholar
83. Liao MC, Van Nostrand WE . Degradation of soluble and fibrillar amyloid beta-protein by matrix metalloproteinase (MT1-MMP) in vitro. Biochemistry 2010; 49: 1127–1136.
    Article CAS PubMed Google Scholar
84. Krishnan V, Nestler EJ . The molecular neurobiology of depression. Nature 2008; 455: 894–902.
    Article CAS PubMed PubMed Central Google Scholar
85. Kendler KS, Karkowski LM, Prescott CA . Causal relationship between stressful life events and the onset of major depression. Am J Psychiatry 1999; 156: 837–841.
    Article CAS PubMed Google Scholar
86. Strekalova T, Spanagel R, Bartsch D, Henn FA, Gass P . Stress-induced anhedonia in mice is associated with deficits in forced swimming and exploration. Neuropsychopharmacology 2004; 29: 2007–2017.
    Article PubMed Google Scholar
87. Jayatissa MN, Bisgaard C, Tingstrom A, Papp M, Wiborg O . Hippocampal cytogenesis correlates to escitalopram-mediated recovery in a chronic mild stress rat model of depression. Neuropsychopharmacology 2006; 31: 2395–2404.
    Article CAS PubMed Google Scholar
88. Bisgaard CF, Jayatissa MN, Enghild JJ, Sanchez C, Artemychyn R, Wiborg O . 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. J Mol Neurosci 2007; 32: 132–144.
    Article CAS PubMed Google Scholar
89. Schulz D, Mirrione MM, Henn FA . Cognitive aspects of congenital learned helplessness and its reversal by the monoamine oxidase (MAO)-B inhibitor deprenyl. Neurobiol Learn Mem 2010; 93: 291–301.
    Article CAS PubMed Google Scholar

Download references

Acknowledgements

This work was supported the Davee Foundation, the RD Foundation, and NIH grant MH077234. The authors also wish to thank Laura Sittig, Elif Tunc-Ozcan, and Timothy Ullmann for helpful suggestions on the manuscript and Drs. Frasier Aird and Claire Will for early work on the endogenous depression model.

Author information

Authors and Affiliations

1. Department of Psychiatry and Behavioral Sciences, Feinberg School of Medicine, Northwestern University, Chicago, IL, USA
   B M Andrus, K Blizinsky, K Dennis, P K Shukla, D J Schaffer, J Radulovic & E E Redei
2. The Jackson Laboratory, Bar Harbor, ME, USA
   P T Vedell & G A Churchill

Authors

1. B M Andrus
   You can also search for this author inPubMed Google Scholar
2. K Blizinsky
   You can also search for this author inPubMed Google Scholar
3. P T Vedell
   You can also search for this author inPubMed Google Scholar
4. K Dennis
   You can also search for this author inPubMed Google Scholar
5. P K Shukla
   You can also search for this author inPubMed Google Scholar
6. D J Schaffer
   You can also search for this author inPubMed Google Scholar
7. J Radulovic
   You can also search for this author inPubMed Google Scholar
8. G A Churchill
   You can also search for this author inPubMed Google Scholar
9. E E Redei
   You can also search for this author inPubMed Google Scholar

Corresponding author

Correspondence toE E Redei.

Ethics declarations

Competing interests

The authors declare no conflict of interest.

Additional information

Supplementary Information accompanies the paper on the Molecular Psychiatry website

Supplementary information

PowerPoint slides

Rights and permissions

This work is licensed under the Creative Commons Attribution-NonCommercial-No Derivative Works 3.0 Unported License. To view a copy of this license, visit http://creativecommons.org/licenses/by-nc-nd/3.0/

Reprints and permissions

About this article

Cite this article

Andrus, B., Blizinsky, K., Vedell, P. et al. Gene expression patterns in the hippocampus and amygdala of endogenous depression and chronic stress models.Mol Psychiatry 17, 49–61 (2012). https://doi.org/10.1038/mp.2010.119

Download citation

• Received: 06 May 2010
• Revised: 05 October 2010
• Accepted: 11 October 2010
• Published: 16 November 2010
• Issue Date: January 2012
• DOI: https://doi.org/10.1038/mp.2010.119

Keywords