Molecular evidence for BDNF- and GABA-related dysfunctions in the amygdala of female subjects with Major Depression (original) (raw)
. Author manuscript; available in PMC: 2013 May 1.
Published in final edited form as: Mol Psychiatry. 2011 Sep 13;17(11):1130–1142. doi: 10.1038/mp.2011.113
Abstract
Women are twice as likely as men to develop major depressive disorder (MDD) and are more prone to recurring episodes. Hence, we tested the hypothesis that the illness may associate with robust molecular changes in female subjects, and investigated large-scale gene expression in the postmortem brain of MDD subjects paired with matched controls (n=21 pairs). We focused on the lateral/basolateral/basomedian (LBNC) complex of the amygdala as a neural hub of mood regulation affected in MDD. Among the most robust findings were downregulated transcripts for genes coding for GABA interneuron-related peptides, including somatostatin (SST), tachykinin, neuropeptide Y (NPY) and cortistatin, in a pattern reminiscent to that previously reported in mice with low BDNF. Changes were confirmed by quantitative PCR and not explained by demographic, technical or known clinical parameters. BDNF itself was significantly downregulated at the RNA and protein levels in MDD subjects. Investigating putative mechanisms, we show that this core MDD-related gene profile (including SST, NPY, TAC1, RGS4, CORT) is recapitulated by complementary patterns in mice with constitutive (BDNF-heterozygous) or activity-dependent (Exon IV knockout) decreases in BDNF function, with a common effect on SST and NPY. Together, these results provide both direct (low RNA/protein) and indirect (low BDNF-dependent gene pattern) evidence for reduced BDNF function in the amygdala of female subjects with MDD. Supporting studies in mutant mice models suggest a complex mechanism of low constitutive and activity-dependent BDNF function in MDD, particularly affecting SST/NPY-related GABA neurons, thus linking the neurotrophic and GABA hypotheses of depression.
Keywords: Major depression, female, amygdala, BDNF, GABA, somatostatin, neurotrophic hypothesis
INTRODUCTION
Major depression is a multifactorial psychiatric disorder that is responsible for substantial disability worldwide. Despite a high morbidity and mortality, its etiology and pathophysiology are not well defined. Notably, women are more likely to become depressed than men, display higher severity and symptom number, have frequent co-morbid anxiety symptoms, and are more prone to recurring episodes. However, the underlying biological vulnerabilities are not fully characterized yet.
Low neurotrophic factor support has been proposed as a unifying hypothesis for reduced cell numbers in frontal cortex1,2 and amygdala,3,4 reduced hippocampal volume.5,6 This interpretation is supported by indirect evidence in rodents showing that antidepressant treatments increase BDNF (protein) and BDNF (mRNA) levels, hippocampal BDNF infusion is sufficient to produce antidepressant-like effects,7,8 and BDNF appears necessary for some behavioral responses to antidepressant treatments.9,10 Direct evidence is more sparse and includes low circulating peripheral BDNF levels which are normalized by antidepressant treatment,11,12 and one study showing reduced pro-BDNF13 and BDNF levels14 in post-mortem hippocampal tissue of depressed patients.15 Associations between BDNF signaling and suicide completion have been made, including reduced BDNF levels in post-mortem hippocampal or midbrain,16–18 decrease in TrkB.T1 expression19 hypermethylation of BDNF promoter/exon IV20 and an increased risk for violent suicide in subjects carrying a Met allele at the BDNF Val66Met polymorphism.21
Evidence in MDD patients also suggests an impaired excitation/inhibition balance that is potentially mediated by decreased GABA content, as observed by microarray,22 imaging studies in occipital and frontal cortices23,24 or by transcranial magnetic stimulation paradigms.25 We recently reported a downregulation of the neuropeptide somatostatin (SST) in the dorsolateral (DLPFC)26 and subgenual cingulate (sgACC) cortices27 of MDD subjects. Expression of SST is dependent on BDNF,28 and identifies a specific GABA interneuron subtype.29 Notably, the downregulation of SST was consistently more robust in female MDD subjects.
At the neural network level, changes in the structure, function and coordinated activity of cortical and subcortical brain regions are thought to underlie the mood regulation deficit in depression.30 Our study focused on the amygdala as a critical component of this corticolimbic circuit of mood regulation.31 Evidence for misregulated amygdala function includes decreased oligodendrocyte numbers and markers3,4 and altered/sustained amygdala function in depressed subjects exposed to negative emotionally salient stimuli.32 Here we investigated molecular changes in the amygdala of female subjects with MDD and report robust downregulation of BDNF and of several BDNF-dependent genes. Notably, this deregulated pattern was recreated by low constitutive or reduced activity-dependent BDNF function in mutant mice and robustly affected SST expression, hence providing evidence in support of low neurotrophic and altered GABA-related functions in MDD.
MATERIAL AND METHODS
Human postmortem subjects
For all subjects, consensus DSM-IV diagnoses of MDD were made by an independent committee of experienced clinical research scientists at a case conference utilizing information obtained from clinical records, toxicology exam and a standardized psychological autopsy.33 This latter incorporates a structured interview, conducted by a licensed clinical psychologist with family members of the index subject, to assess diagnosis, psychopathology, medical, social and family histories, as well as history of substance abuse. A symptom score reflecting disease severity was calculated based on the presence at time of death of nine MDD symptoms: depressed mood, anhedonia, appetite disturbance, sleep disturbance, psychomotor change, anergia, self-recrimination, diminished ability to concentrate or make decision, and suicidality.27,34 Each symptom is scored 0 (absence) or 1 (presence) and the symptom score is sum of the individual symptom scores.
21 pairs of subjects were analyzed (Table 1), consisting of female subjects with MDD and control subjects matched for sex, race, and as closely as possible for age, postmortem interval and brain pH. The pairing protocol has been previously validated as reducing signal variability and inducing higher overall correlations of gene transcript levels compared to non-matched pairs.34 All subjects died suddenly without prolonged agonal periods. All cases and controls except one pair were white Caucasian. Brains were analyzed for adequate brain pH (>6.4) and RNA integrity by optical density (OD≥ 1.6) and Agilent bioanalyzer analysis (Agilent Technologies, Palo Alto, CA; RIN expert scoring system ≥ 7) as described.34 Accordingly, subject groups did not differ in mean age, PMI, brain pH, RNA integrity number (RIN), or RNA ratio (p>0.05). Rates of death by suicide, disease recurrence, evidence for antidepressant treatment at time of death, and alcohol dependence in MDD subjects were recorded. Toxicological screens on peripheral fluids identified the presence of at least one antidepressant in 14 subjects. (Table 1), while no psychoactive substance was detected in controls.
Table 1. Major Depression and Control Cohorts.
Y: yes, N: no, U: unknown; H: History of treatment, T: toxicological screen at time of death; AD: antidepressant unspecified, Ami: amitriptyline, Bup: bupropion, Cit: citalopram, Dox: doxepin, Esc: escitalopram, Flv: fluvoxamine, Flx: fluoxetine, Nef: nefazodone, Li: lithium, Mirt: mirtazapine, Nor: nortriptyline, Ola: olanzapine, Prx: paroxetine, Quet: quetiapine, Risp: risperidone, Sert: sertraline, Traz: trazodone, Venl: venlafaxine, Zip: ziprasidone
| Major Depression Group | Control Group | ||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Pairs | Case | Age (Years) | Postmortem Interval(hours) | pH | RNA ratio | RNA Integrity number | Suicide | Recurrent Episode | Antidepressant | Alcohol Dependence atDeath | Case | Age (Years) | Postmortem Interval(hours) | pH | RNA ratio | RNA Integrity number | Alcohol Dependence atDeath |
| 1 | 564 | 56 | 16.8 | 7.02 | 1.94 | 9.2 | Y | N | FlxT, DoxT | N | 568 | 60 | 9.5 | 6.88 | 1.9 | 8.7 | N |
| 2 | 666 | 16 | 10 | 7.29 | 2 | 9.4 | N | N | PrxH | N | 10013 | 16 | 9.3 | 6.69 | 1.8 | 9 | N |
| 3 | 803 | 65 | 18 | 6.96 | 1.86 | 9 | N | Y | CitH, SertH, NorT | N | 1466 | 64 | 20 | 6.74 | 2.01 | 8.8 | N |
| 4 | 934 | 54 | 17.9 | 6.2 | 1.16 | 8.2 | N | Y | CitH,T, MirtH,T, AmiH,RispH | N | 1247 | 58 | 22.7 | 6.37 | 1.28 | 8.4 | N |
| 5 | 967 | 40 | 22.2 | 6.56 | 1.64 | 7.4 | N | Y | FlxH, DoxH, LiH, PrxH, SertH | Y | 1282 | 39 | 24.5 | 6.84 | 1.32 | 7.5 | N |
| 6 | 986 | 53 | 11.9 | 6.7 | 1.83 | 8.7 | N | Y | ADH, FlxT | N | 575 | 55 | 11.3 | 6.81 | 1.79 | 9.6 | N |
| 1041 | 52 | 10.3 | 6.51 | 1.48 | 8.4 | N | Y | FlvH, MirtH, NefH, OlaH, QuetH, RispH, TrazH, VenlH, ZipH | Y | 1391 | 51 | 7.8 | 6.57 | 1.59 | 7.1 | N | |
| 7 | 1157 | 26 | 13.4 | 6.4 | 1.48 | 7.8 | Y | Y | CitH,T, PrxH | N | 1034 | 23 | 8.5 | 6.11 | 1.96 | 7.8 | N |
| 8 | 1190 | 47 | 22.3 | 6.55 | 1.64 | 8 | Y | Y | FlxH, TrazH | Y | 567 | 46 | 15 | 6.77 | 2.26 | 8.9 | N |
| 9 | 1202 | 39 | 11.2 | 6.39 | 1.79 | 8 | N | Y | EscH, BupH | N | 840 | 41 | 15.4 | 6.8 | 1.98 | 9.1 | N |
| 10 | 1221 | 28 | 24.8 | 6.61 | 1.83 | 7.2 | N | Y | MirtH | N | 546 | 37 | 23.5 | 6.74 | 1.95 | 8.6 | N |
| 11 | 1249 | 40 | 11.2 | 6.52 | 2.04 | 9 | N | Y | FlxH,T, TrazH | N | 1092 | 40 | 16.6 | 6.83 | 1.68 | 8 | N |
| 12 | 1254 | 39 | 12.8 | 6.37 | 1.92 | 9 | Y | Y | MirtH, SertH | N | 1403 | 45 | 12.3 | 6.67 | 1.8 | 8.2 | N |
| 13 | 1271 | 50 | 18.6 | 6.35 | 1.82 | 8.6 | N | Y | AmiH,T | Y | 1318 | 58 | 18.8 | 6.69 | 1.95 | 7.4 | N |
| 14 | 1289 | 46 | 25 | 6.27 | 1.35 | 7.3 | N | N | U | N | 1280 | 50 | 23.5 | 6.65 | 1.33 | 7.7 | N |
| 15 | 1315 | 28 | 12.4 | 7 | 1.52 | 7.9 | Y | N | PrxH, BupH | Y | 1099 | 24 | 9.1 | 6.46 | 1.86 | 8.6 | N |
| 16 | 1332 | 46 | 17.5 | 6.49 | 1.6 | 8.9 | N | Y | FlxH | N | 627 | 43 | 14.1 | 7.09 | 1 | 7 | N |
| 17 | 1356 | 60 | 20.6 | 6.06 | 1.78 | 8.5 | N | Y | CitH,T, NorH, PrxH, NefH, TrazH, BupH, SertH | Y | 818 | 67 | 24 | 7.06 | 1.48 | 8.4 | N |
| 18 | 1360 | 59 | 18.1 | 6.43 | 1.41 | 7.6 | Y | N | AmiH, CitH, FlxH, TrazH | N | 1081 | 57 | 14.9 | 6.78 | 1.8 | 9 | N |
| 19 | 1408 | 37 | 15.5 | 6.56 | 1.58 | 7 | N | Y | VenlH, PrxH, SertH, BupT | Y | 1196 | 36 | 14.5 | 6.44 | 1.84 | 8.2 | Y |
| 20 | 10028 | 72 | 23.1 | 6.66 | 1.37 | 7 | Y | N | N | N | 1355 | 74 | 24.9 | 6.62 | 1.89 | 7 | N |
| Mean | 46.0 | 16.2 | 6.7 | 1.67 | 8.2 | 46.8 | 16.5 | 6.6 | 1.74 | 8.2 | |||||||
| SEM | 3.0 | 1.3 | 0.04 | 0.05 | 0.16 | 3.2 | 1.1 | 0.06 | 0.07 | 0.16 |
Microarray Samples
Rostral amygdala samples enriched in lateral, basolateral, and basomedian nuclei were delineated as described previously,34 dissected from frozen coronal blocks ~2cm caudal to the temporal pole. Total RNA was processed for Illumina HT12 microarray analysis according to the microarray manufacturer’s protocol (Illumina Inc, San Diego, CA).
Statistical analysis
Gene selection was made using a random-intercept statistical model (RIM) to adjust for clinical, demographic and technical variables (Details in supplements). Multiple testing was adjusted by the Benjamini-Hochberg method.35 This resulted in a stringent selection of 116 differentially expressed genes, termed "q-list" (Adjusted p-value=Q-value<0.05; mean change>20%). In addition, an exploratory list of genes was obtained using uncorrected parametric and non-parametric statistical tests (RIM, paired T-test and Wilcoxon test), using moderate statistical stringencies for gene selection (p<0.05 in at least one statistical test), as proposed previously.34 The resulting exploratory list (n=4131 genes; 307 genes with fold change >20%) may carry a higher rate of false positives at the single gene level, but allowed investigation of cumulative effects over larger sets of genes and pathways (See BDNF-related genes and Ingenuity Pathway analysis).
In situ hybridization
Antisense and sense riboprobes for human SST mRNA were transcribed in the presence of [35S]-CTP (Amersham Biosciences, Piscataway, NJ) as described.36 The sections and [14C]-standards were exposed on the same BioMax MR film (Kodak, Rochester, NY) for 3 days and analyzed with the MCID software.
Protein isolation, Prepro-SST and BDNF Immunoblotting
Acetone precipitation of proteins was carried out following RNA extraction. Western blot analysis was performed as described.34 Dual signals were detected using the LI-COR Odyssey Infrared imaging system (LI-COR Biosciences, Lincoln, NE, USA), and SST or BDNF signal ratios to β-actin were calculated. Samples were processed in matched pairs on the same gel and results were replicated for a total of three different western blots.
BDNF +/− and exon-IV KO mice
Female BDNF+/+ and +/− mice (3 to 4 months old) were bred on a mixed S129/Sv×C57BL/6 genetic background. Heterozygous adult female mice with one functional BDNF allele (BDNF_HZ_) exhibit 50% reduction of BDNF mRNA levels in the hippocampus.37 Female BDNF exon IV KO (BDNF_IV-KO_) mice were crossed on C57BL/6 as described.38 Brains were rapidly removed and flash-frozen on dry ice. Caudo-rostral coronal sections were obtained on a cryostat until the caudal amygdala was observed (Fig #48,39). Brains were then flipped and cut rostro-caudally until observing the rostral amygdala (Fig #40). Bi-lateral amygdalae were then micro-punctured using 0.5mm diameter punches and stored in Trizol at −80°C.
Real-time quantitative Polymerase Chain Reaction (qPCR)
Small PCR products were amplified in quadruplets on a Mastercycler real-time PCR machine (Eppendorf, Hamburg, Germany), using universal PCR conditions. Results were calculated as the geometric mean of relative intensities compared to three validated internal controls (actin, glyceraldehyde-3- phosphate dehydrogenase, and cyclophilin G).
See detailed methods in Supplements.
RESULTS
Altered gene expression in the amygdala of female subjects with MDD suggests a low BDNF-dependent transcriptional drive affecting SST GABA interneurons
Using a microarray approach in the amygdala of postmortem female MDD subjects compared to matched controls and taking clinical and technical parameters into account, we identified 116 gene transcripts as differentially expressed after controlling for multiple testing (q-value<0.05; effect size>20%; q-list, Table 2 and Supplementary Table ST1). Based on previous identification of relative glial/neuronal enrichments in gene transcript origin,40 23 gene transcript changes were of neuronal origin and mostly down-regulated (~80%), 56 of mixed neuronal/glial origin (22 down, 34 up) and 31 of enriched glial origin and mostly upregulated (~87%). Notably, among the most down-regulated genes in each subcategory were SST, NPY, RGS4, CORT, and TAC1 (Table 2), in a pattern strikingly reminiscent to that previously reported in the cortex of male mice with low BDNF,28 thus suggesting that a deficit in BDNF function may have contributed to the gene pattern observed in MDD. We further explored this putative BDNF link by investigating a broader list of genes identified as BDNF-related in the Ingenuity Pathway database (n=283 genes). 52 BDNF-related genes displayed significant uncorrected MDD-affected changes (Supplementary Table 2) and BDNF-related genes were significantly over-represented in the q-list of MDD-affected gene list (6% vs. 1.2% in Ingenuity database; χ2=17.4, df=1, p<0.0001). Further Ingenuity-based exploratory functional analyses resulted in buildings of gene networks and canonical pathways coherent with our findings (Supplementary Figure S1).
Table 2. Summary of the q-list Core Genes Significantly Affected in Female MDD subjects.
39 out of 116 genes differently expressed in MDD (MDD) versus control are listed (Full list in Supplementary Table ST2). "Neuronal", "Neuronal-Glial", and "Glial" refers to enrichments of transcript origin (Supplements).40
| IlluminaProbeset ID | GeneTitle | GeneSymbol | RIM p-value | RIM q-value | All MDDSubjects(alr) | WM/GMFoldchange |
|---|---|---|---|---|---|---|
| 13 out of the 24 Neuronal Genes of the full q-list | ||||||
| ILMN_1812824 | Somatostatin | SST | 4.09E-05 | 0.020 | −0.84 | −1.96 |
| ILMN_1765966 | Chromogranin B (secretogranin 1) | CHGB | 7.95E-04 | 0.034 | −0.54 | −3.76 |
| ILMN_1679984 | Zinc finger, CCHC domain containing 12 | ZCCHC12 | 1.64E-04 | 0.024 | −0.53 | −4.24 |
| ILMN_1697512 | Solute carrier family 32 (GABA vesicular transporter), member 1 | SLC32A1 | 4.91E-04 | 0.030 | −0.48 | −7.27 |
| ILMN_1685834 | Amphiphysin, transcript variant 1 | AMPH | 1.63E-03 | 0.044 | −0.39 | −2.12 |
| ILMN_2384409 | Tachykinin, precursor 1, transcript variant alpha, | TAC1 | 2.19E-03 | 0.048 | −0.38 | −1.50 |
| ILMN_1810604 | ELMO/CED-12 domain containing 1 | ELMOD1 | 2.17E-03 | 0.048 | −0.38 | −2.72 |
| ILMN_1806147 | Guanine nucleotide binding protein (G protein), gamma 3 | GNG3 | 1.35E-04 | 0.023 | −0.31 | −1.91 |
| ILMN_1654632 | Regulator of G-protein signaling 7 binding protein | RGS7BP | 2.30E-03 | 0.049 | −0.31 | −5.46 |
| ILMN_1765966 | Chromogranin B (secretogranin 1) | CHGB | 7.95E-04 | 0.034 | −0.54 | −3.76 |
| ILMN_2354547 | Tumor suppressor candidate 3, transcript variant 1 | TUSC3 | 9.76E-04 | 0.036 | −0.34 | −1.89 |
| ILMN_1736154 | ProSAPiP1 protein | ProSAPiP1 | 6.97E-04 | 0.033 | 0.31 | −2.60 |
| ILMN_1760798 | Ryanodine receptor 2 (cardiac) | RYR2 | 7.82E-04 | 0.034 | 0.40 | −2.65 |
| 13 out of the 57 Neuronal/Glial Genes of the full q-list | ||||||
| ILMN_2071186 | Cortistatin | CORT | 1.72E-04 | 0.025 | −0.58 | −1.06 |
| ILMN_1731062 | Neuropeptide Y | NPY | 3.70E-04 | 0.029 | −0.52 | −1.18 |
| ILMN_1729165 | Transcription elongation factor A (SII)-like 6 | TCEAL6 | 8.23E-04 | 0.035 | −0.40 | 1.41 |
| ILMN_1771286 | PREDICTED: similar to phosphodiesterase 4D interacting protein | LOC653513 | 3.22E-04 | 0.028 | −0.39 | 1.17 |
| ILMN_1690397 | Dynein, cytoplasmic 1, intermediate chain 1 | DYNC1I1 | 1.02E-03 | 0.037 | −0.35 | −1.20 |
| ILMN_2399304 | Neuron navigator 2, transcript variant 2, | NAV2 | 5.92E-04 | 0.032 | 0.34 | 1.27 |
| ILMN_1768962 | A kinase (PRKA) anchor protein 8-like | AKAP8L | 4.14E-04 | 0.030 | 0.35 | 1.02 |
| ILMN_1663042 | Syndecan 4 | SDC4 | 3.34E-04 | 0.029 | 0.36 | −1.00 |
| ILMN_1655611 | Teashirt zinc finger homeobox 2 | TSHZ2 | 2.02E-03 | 0.047 | 0.38 | −1.04 |
| ILMN_1696757 | Tetratricopeptide repeat domain 14, transcript variant 2, | TTC14 | 5.92E-04 | 0.032 | 0.38 | 1.02 |
| ILMN_1682775 | Endothelin 1 | EDN1 | 2.50E-04 | 0.027 | 0.42 | 1.26 |
| ILMN_2078547 | Hypothetical protein HSPC268 | HSPC268 | 7.69E-04 | 0.034 | 0.42 | 1.02 |
| ILMN_1744897 | Potassium intermediate/small conductance calcium-activated channel, subfamily N, member 3 | KCNN3 | 5.09E-04 | 0.030 | 0.45 | −1.39 |
| 13 out of the 32 Glial Genes of the full q-list | ||||||
| ILMN_1758067 | Regulator of G-protein signalling 4 | RGS4 | 1.04E-03 | 0.038 | −0.47 | 1.52 |
| ILMN_2320164 | Purinergic receptor P2Y, G-protein coupled 12, transcript variant 1, | P2RY12 | 8.38E-04 | 0.035 | −0.46 | 2.49 |
| ILMN_2402172 | Septin 4, transcript variant 3, | SEPT4 | 8.25E-04 | 0.035 | 0.35 | 2.45 |
| ILMN_1665686 | Family with sequence similarity 38, member B | FAM38B | 2.45E-04 | 0.027 | 0.36 | 5.72 |
| ILMN_1810420 | Dysferlin, limb girdle muscular dystrophy 2B | DYSF | 1.75E-03 | 0.046 | 0.36 | 2.95 |
| ILMN_1670881 | Carbohydrate (N-acetylglucosamine 6-O) sulfotransferase 6 | CHST6 | 9.30E-05 | 0.020 | 0.36 | 3.92 |
| ILMN_1697176 | Glial fibrillary acidic protein | GFAP | 9.40E-04 | 0.036 | 0.37 | 3.60 |
| ILMN_1659895 | Moesin | MSN | 2.53E-04 | 0.027 | 0.37 | 1.91 |
| ILMN_1737631 | Progestin and adipoQ receptor family member VI, transcript 1 | PAQR6 | 7.85E-04 | 0.034 | 0.37 | 2.47 |
| ILMN_2323508 | Chromosome 9 open reading frame 58, transcript variant 2 | C9orf58 | 2.20E-03 | 0.048 | 0.38 | 2.25 |
| ILMN_1752668 | Dishevelled associated activator of morphogenesis 2 | DAAM2 | 3.78E-04 | 0.029 | 0.42 | 2.25 |
| ILMN_1750271 | Myelin-associated oligodendrocyte basic protein | MOBP | 4.33E-04 | 0.030 | 0.51 | 2.46 |
| ILMN_2331544 | Myelin basic protein transcript variant 7 | MBP | 6.95E-04 | 0.033 | 0.64 | 2.55 |
Confirmation of altered gene expression by qPCR
Microarray results for 13 genes in the MDD q-list (AMPH, CDK5RAP2, CORT, GFAP, KCNG1, MBP, MOBP, NPY, RGS4, SLC32A1, SST, SSTR1, TAC1) were confirmed by qPCR (Supplementary Table ST3), with a high Array/qPCR concordance (Pearson correlation r=0.95, p<1e–6) and similar directionality of expression for all genes (Figure 1A).
Figure 1. qPCR validation of altered gene expression and reduced BDNF transcript/protein in the amygdala of female MDD subjects.
(A) Confirmation of microarray results by qPCR (See also Supplementary Table ST4). (B) BDNF-IXd coding exon transcript expression in relative expression (%) of the control group mean (*, p<0.05). (C) Mature BDNF (m-BDNF) and pro-BDNF protein relative immunoreactivity migrates at the expected ~15 and 25kD size respectively. Examples of 2 sample pairs of control (C) and MDD (D) on the same gel. (D) Relative level of m-BDNF in MDD in function of respective paired CTRL (◆ADD-treated MDD ◇ADD-free MDD
Average of population) (*, p<0.05).
Exploratory analysis (i.e. unadjusted p-values) identified additional downregulated GABA-related genes (GAD1, GAD2, GABRA1, GABRA5, calretinin (CALB2)) and genes coding for receptors (HTR3A) found on GABAergic neurons in the amygdala,41 but these findings were not reliably confirmed (Supplementary Table ST3), potentially reflecting array false positives or a diluting effect from other unaffected GABA neuron subsets in the qPCR cDNA samples. For instance, parvalbumin (PV), a marker of a different GABA neuron subset, was not affected in our cohort.
Low BDNF transcript and protein levels in the amygdala of female subjects with MDD
Due to low detection by array, we measured BDNF by qPCR and western blotting. Results indicate significant decreases in MDD subjects of mRNA levels corresponding to the BDNF coding sequence (BDNF exon IXd; −22%; p<0.05, Figure 1B) and protein level, for both the precursor (pro-) and mature (m-) forms (respectively, −27% and −30%; p<0.05, Figure 1C–E). Low BDNF was observed regardless of the presence of antidepressant treatment (Figure 1E).
BDNF-IXd expression was co-regulated across all samples with the four core genes involved in the BDNF-related signature (SST, NPY, RGS4 and TAC1; Average R=0.39, p<0.05), suggesting a functional link, however the coregulation with BDNF was lower than mutual coregulation within those four genes (Average R=0.57, p<0.01) (Supplementary Table ST4), suggesting the presence of additional regulatory events.
BDNF transcription is initiated from at least nine promoters that respond to differing stimuli and drive transcription of a short-noncoding exon spliced to a common coding exon.42 For example, promoter IV is highly responsive to neuronal activity and induces activity-dependent BDNF expression in vitro and in vivo.38 Using qPCR, we report no difference in expression levels for selected activity-dependent (I, IV and IXa) BDNF promoter/non-coding exons, although we note the limits of assessing short term and activity-dependent gene activities under variable postmortem intervals for brain collection. No changes were observed in mRNA coding for TrkB (full length or truncated forms; Supplementary Figure S2).
Could BDNF-related changes be explained by other clinical/demographic factors?
Although clinical and technical parameters were accounted for in the single-gene statistical models, we further investigated potential effects of selected parameters on the 52 BDNF-related genes as a whole (Supplementary Table ST2).
Age
We previously reported that disease-associated genes, including BDNF, are robustly modulated by age.43 Accordingly, we observed a correlation between age and the 52 BDNF related-gene set (averaged relative change and age; R=0.44, p=0.004; See Supplementary Figure S3); however, the effect of MDD on BDNF-related genes was still highly significant after ANCOVA analysis with age as cofactor (p=8.8e−4). The details of relation between age and gene expression in MDD will be discussed elsewhere.
Antidepressant treatment was detected in 2/3 of subjects at time of death, but data on long-term exposure and treatment efficacy was not available. BDNF levels in MDD subjects did not correspond to antidepressant treatment. Moreover the levels of expression changes for 52 BDNF-related genes with nominal MDD p-values<0.05 were not different between treated and untreated MDD subjects (Supplementary Table ST2; Pearson correlation, r=0.91, p<1e−7). We further investigated potential antidepressant effects on orthologous genes in a mouse model of antidepressant response. In that study, 4-week chronic fluoxetine exposure reversed the high emotionality phenotype induced by chronic mild stress in mice.44 With the exception of two genes (ACAT2 and RAPGEF6) antidepressant did not affect expression of orthologous genes from the MDD-related q-list (Table 2 and Supplementary Tables ST1–2).
Alcohol use was used as a covariate in all statistical analysis performed and did not interfere with any of the positive results provided herein.
Hormonal Status
The cohort encompasses a large age range and information on hormonal status and replacement therapy was largely not available. So, we hypothesized that the transcriptional activation of genes involved in regulating estrogen function (CYP19A1, ESR1, ESR2, PGR and SULT1A1) or for genes proposed as peripheral biomarkers of menopausal status (AMH),45 may represent a useful proxy assay for hormonal status. The expression levels of these genes did not correlate with BDNF-related transcript changes (Supplementary Figure S4).
Disease severity
For each depressed subject, a symptom score reflecting the severity of the disease was calculated based on the presence at time of death of nine MDD episode symptoms (Methods)34. Score values showed no correlation with changes in the expression of any of the BDNF-related gene, suggesting that intensity of changes were independent of disease' severity (Supplementary Table ST2).
Death by suicide
The disease effect sizes on BDNF-related genes were reduced in suicide versus non-suicide MDD subjects (Supplementary Figure S5), but BDNF mRNA and pro-BDNF protein levels (but not mBDNF) were further decreased in depressed suicide victims (mRNA: alrsuicide/ctrl=−0.41 vs alrnon-suicide/ctrl=−0.28, p=0.34; pro-BDNF protein: alrsuicide/ctrl =−0.62 vs alrnon-suicide/ctrl =−0.47, p=0.28), although not significantly, potentially due to reduced sample size.
Translational investigation of putative mechanisms in mouse models with genetically-altered changes in BDNF level or function
We next investigated if similar changes were observed downstream from low BDNF in the amygdala of mutant mice with altered BDNF function. To differentiate the putative contribution of constitutive and activity-dependent BDNF functions, we investigated gene transcript changes in mice heterozygous for a constitutive deletion of the BDNF gene (BDNF_HZ_), which display ~50% less BDNF;37 and mice with a targeted disruption of exon IV (BDNF_IV-KO_), which results in a near complete blockade of activity-dependent BDNF protein expression.38 Using qPCR, 6 core genes of interest were investigated, based on confirmed changes in MDD and known BDNF-dependency, including SNAP25, a BDNF-regulated gene member46 of the SNARE complex, also observed downregulated in MDD47 (Table 3). In addition, we evaluated the expression of 5 other GABA-related genes displaying less robust changes (array results not qPCR-confirmed), and parvalbumin (PV) as a control GABA marker not affected in our cohort. qPCR results confirmed the BDNF-dependency for most genes in amygdala, but neither line of mutant mice mimicked the full pattern of MDD changes. Instead, gene changes in BDNF_HZ_ and BDNF_IV-KO_ displayed complementary profiles (Table 3), which together recreated the MDD profile, suggesting a combined low activity-dependent and constitutive BDNF function in MDD. SST and NPY were identified at the intersection of the BDNF_HZ_ and BDNF_IV-KO_ profiles, hence identifying SST/NPY GABA interneurons as the most vulnerable GABA interneuron subtype to low BDNF function (Table 2).
Table 3. Expression profile of BDNF- and GABA-related genes in MDD subject compared to BDNF_HZ_ and BDNFK_IV-KO_ mice.
Alr, average log ratio (MDD/control, BDNFHZ/WT and BDNFIV-KO/WT). All unadjusted p-values in Array analyses were obtained after RIM analyses, except (#) was obtained by Wilcoxon Test
| Human | Mouse | |||||||
| BDNFHZ | BDNFKIV-KO | |||||||
| Array | qPCR | qPCR | qPCR | |||||
| Alr | p-value | Alr | p-value | Alr | p-value | Alr | p-value | |
| BDNF-dependent | ||||||||
| TAC1 | −0.38 | 2.19E-03 | −0.75 | 0.017 | −1.70 | 0.010 | 0.50 | 0.189 |
| RGS4 | −0.47 | 1.04E-03 | −0.68 | 0.002 | −0.59 | 0.050 | −0.13 | 0.344 |
| NPY | −0.52 | 3.70E-04 | −0.92 | 0.001 | −0.46 | 0.034 | −0.33 | 0.014 |
| SST | −0.84 | 4.09E-05 | −1.53 | 0.000 | −0.46 | 0.099 | −0.31 | 0.016 |
| CORT | −0.58 | 1.72E-04 | −0.73 | 0.006 | −0.08 | 0.320 | −0.92 | 0.001 |
| SNAP25 | −0.44 | 4.15E-03 | −0.24 | 0.091 | −0.04 | 0.459 | −0.35 | 0.034 |
| Other GABA markers | ||||||||
| CALB2 | −0.24 | 0,031# | −0.21 | 0.208 | −0.17 | 0.415 | −0.06 | 0.436 |
| GAD1 | −0.22 | 0.004 | −0.28 | 0.076 | −0.45 | 0.052 | 0.06 | 0.353 |
| GABRA1 | −0.34 | 0.019 | 0.14 | 0.252 | −0.10 | 0.406 | −0.18 | 0.182 |
| PVALB | −0.09 | 0.541 | 0.04 | 0.391 | 0.16 | 0.364 | −0.12 | 0.342 |
| SLC6A1 (GAT-1) | 0.06 | 0.471 | −0.08 | 0.362 | −0.36 | 0.038 | 0.08 | 0.443 |
Low SST transcript and protein levels in the amygdala of female subjects with MDD
Low SST was confirmed by qPCR (Supplementary Table ST3). Using in situ hybridization, SST mRNA expression was significantly decreased in the lateral (p<0.01) and basomedial (p<0.05) amygdaloid nuclei in MDD subjects compared to control subjects (Figure 2A–C). SST expression was overall lower in the basolateral nucleus and not significantly different between control and MDD subjects (p=0.146). Low SST was also observed at the precursor protein level by quantitative Western blot analysis (Figure 2D–E).48 Finally, the correlation between mRNA levels of BDNF with SST levels (r=0.56, p<0.0001) was also observed at the protein level (r=0.59, p<0.0001).
Figure 2. Reduced SST mRNA and protein levels in amygdala of female MDD subjects.
(A) SST mRNA expression assessed by autoradiographic optical density measures in LBNC sub-regions (*, p<0.05). (B) LBNC details78 (C) In situ hybridization film autoradiograms of a MDD subject (left) displaying robust SST mRNA downregulation compared to a matched control (right) (D) Immunoblot of prepro-SST in CTRL and MDD subjects. (E) Relative prepro-SST levels in MDD in function of respective paired CTRL (◆Antidepressant-treated MDD ◇Antidepressant-free MDD;
Average of population).
DISCUSSION
We report robust gene transcript changes in the postmortem amygdala of female subjects with MDD. The pattern of altered transcripts in the most significantly dysregulated genes was reminiscent of changes observed in mice with low BDNF. Accordingly, we detected in human MDD subjects a significant downregulation of BDNF at the RNA and protein levels, confirmed the presence of a broader profile of low BDNF-related changes, and independently verified results by qPCR for several genes, including SST and NPY, two markers for a specific subset of GABA interneurons. We show that the core profile of BDNF-dependent gene dysregulation is recreated by the union of complementary patterns of gene transcript changes downstream from either low constitutive or activity-dependent BDNF function in mice. This suggests a model where parsimonious BDNF deficits may result in distinct molecular profiles, cellular deficits including increased vulnerability of SST-bearing GABA neurons, leading to potentially combinatorial pathological phenotypes. The robust findings observed here in human and rodent cohorts, in concert with supporting SST-related female-biased findings in sgACC and DLPFC,26,27 suggest that female-specific factors (such as organizational or activational hormonal effects) contribute an added biological vulnerability to MDD. In summary, these results provide both direct (low BDNF) and indirect (low-BDNF-dependent gene profile) evidence in support of the low neurotrophin hypothesis of MDD, while also linking it to the low GABA hypothesis of depression.
Low BDNF as an intermediate molecular phenotype in complex brain illnesses
Multiple known etiological factors affect BDNF, including, genetic, naturalistic (aging,49 exercise,50 stress51) and biochemical (serotonin, IGF-1)52 factors. Hence, low BDNF activity is not specific to MDD, and also observed in schizophrenia,53 dementia and neurodegenerative disorders.54 So, low BDNF and associated downstream gene changes can be considered a true intermediate molecular phenotype, which is associated with specific upstream etiological factors and putative downstream molecular/cellular effects (Figure 3). Here, our studies in MDD suggest an origin/cause that includes a combination of activity-dependent and constitutive deficits in BDNF, while the downstream proximal effects include molecular (altered GABA markers) and potentially pathophysiological effects (altered GABA-related microcircuitry function).
Figure 3. A proposed model of a BDNF-SST-GABA molecular intermediate phenotype within the complex multi-module pathophysiology of MDD.
Our results suggest that the frequent observation of low SST in MDD is linked to altered BDNF function. The two components of this "molecular intermediate phenotype" (BDNF & SST) are known to be controlled by sets of unique and/or shared upstream etiological factors. Their downstream impact on integrated GABA function will be determined by the extent of SST-expressing cell dysregulation and by the degree of disease effects on other components of the GABA microcircuitry (e.g. PV, not affected here). So, altered GABA circuitry may be considered a one-scale-higher "cellular intermediate phenotype". The outcome of dysregulated BDNF-SST-GABA function in the amygdala and in an extended corticolimbic neural network is hypothesized to alter mood regulation and contribute to the "affect dysregulation" symptom domain of the illness. In this model, sets of additional molecular/cellular intermediate phenotype (or modules) are under the control of their respective etiological factors, display multiple levels of interactions, but potentially impinge on different neural networks, which in turn mediate distinct symptom domains.
This model can account for how the low function of a common and widely distributed neurotrophic factor such as BDNF can be implicated in multiple disorders, and yet may still lead to disease-specific outcomes. It also provides a perspective on why low BDNF may not be causative of a particular DSM-IV diagnosis, but instead contributes one set of changes to complex molecular pathologies of those illnesses, where “complexity” reflects an array of parallel and interacting biological pathologies (defined as pathogenic modules in Figure 3), each with their own etiological factors and mediating pathways. In MDD, such patterns may result from individual genetic liability, prior biological insults (i.e., infection, inflammation), altered developmental trajectories, or biochemical exposure (i.e., interferon or glucocorticoids) for instance.
Low SST as an obligate and/or critical downstream phenotype of mental illnesses?
Similar remarks could apply to SST downregulation, which has been reported in schizophrenia,36 bipolar disorder,26 Alzheimer’s55, during normal aging49, and across brain regions in MDD.26,27 Here, the rodent findings convincingly suggest that low SST represents a common phenotype downstream of altered BDNF function (Table 3), and could thus be considered an obligate molecular/cellular consequence across BDNF-linked illnesses. In that sense, low SST could be viewed as a consequence of disease-related events, but its contribution (i.e. cause), in the amygdala or other regions, to clinical phenotypes is not known. SST is expressed in a subgroup of ~20% of GABA interneurons that provide delayed and sustained inhibition onto principal pyramidal neurons through dense projections onto dendritic arbors, and to some extent may provide unspecific normalizing inhibition (at least in mouse cortex).56 Hence, low SST is unlikely to remain biologically silent at the GABA microcircuitry level. However, this will depend on (1) whether all SST neurons are affected or whether a portion of SST cells are missing, (2) whether GABA function itself is affected in SST neurons, or (3) if deficits are specific to altered neuromodulatory peptides within those cell population, as SST, NPY and CORT (all reported low here) share similarities in cellular origins and may be implicated in the disease process on their own.
Moreover, the impact of low SST is expected to vary based on the broader GABA microcircuitry pathological context. For instance, changes in PV-expressing GABA neurons which directly target the cell soma and axon (Figure 3) (or other subtypes), are poised to interact with SST neuron-mediated inhibition, and may accordingly translate into distinct pyramidal cell regulation, pathophysiological network activity and potentially symptom-related behavioral outcomes such as schizophrenia and bipolar disorder26,57. Upstream factors affecting SST (e.g. aging, stress, glutamate) also affect additional biological modules (Figure 3). Together, it is expected that the combinatorial recruitment of different modules through various (and potentially interacting) etiological pathways will together result in complex disease-related patterns. We propose that BDNF, BDNF-dependent SST and other related gene changes, represent one of those complex biological modules, with its own risk factors and modulators.
Finally, brain region-specificities in terms of cellular composition and intrinsic vulnerabilities may further determine pathophysiological outcomes. We speculate that the changes observed here in the amygdala, an area specialized in detecting and assessing emotional salience of incoming stimuli, will contribute to altered stimulus integration, and in turn to increased and sustained reactivity of the amygdala in MDD patients,58 although this hypothesis will need to be tested in genetic rodent systems.
Expected partial cues from BDNF- or SST-related genetic mouse models
Post-mortem analyses do not allow discriminating if observed alterations are causative and/or consequential to mental illness, but dissecting disease pathways is possible in mouse models. However, the prediction of the model described above (Figure 3) is that current genetic mutations in mice do not recapitulate the full context and complexity of molecular changes observed in illness and so that they will only partially test causal link with altered behavior. For instance, decreased BDNF has been observed in rodents after chronic stress,59 but not consistently60, and low BDNF by itself is not sufficient to generate an anxious/depressive-like phenotype,61,62 and in some regions can exert antidepressant-like effects.63 On the other hand, disruption of BDNF (global or forebrain-specific)64,65 leads to higher emotionality when combined with exposure to chronic stress in female mice,62,66 while male BDNF_HET_ were as vulnerable as WT to chronic stress.68 Dipping further in BDNF specificity, mice lacking the activity-dependent BDNF promoter-IV appear more prone to develop increased emotionality.67 Similarly, SST gene disruption in mice does not induce behavioral alterations other than motor impairment,68 suggesting that low SST by itself is not sufficient to induce anxious/depressive-like states. However, the behavioral studies in SSTKO mutant mice may not have been optimized to detect subtle baseline/trait changes, were performed in compromised mice with complete loss of SST function, and did not assess increased vulnerability to develop higher emotionality, such as after inducing stress protocols for instance.
Comments and Limitations
- SST and co-expressed genes (NPY, CORT, TAC1) were used as surrogates for GABA function, but it remains to be determined whether markers for GABA function within SST neurons or only neuropeptides are affected. Indeed, SST itself displays antidepressant activity that is GABA-independent. 69
- Sex and species differences. A direct comparison of female and male cohorts will be needed to fully assess the higher vulnerability of females in developing BDNF-related profiles. The goal of current study was to use the increased disease prevalence in female subjects as an enrichment strategy for a potential more homogeneous underlying molecular pathology. Our prior study in the amygdala of male subjects33 suggests that the phenotype may not be as robust in male subjects, although we cannot at this point determine whether this was due to cohort heterogeneity or to smaller sample size. Nevertheless, neuroimaging reports on the impact of altered BDNF function on amygdala response to stressful stimulation, and pathways to depression and anxiety70,71 suggest that the uncovered changes in BDNF pathway may be equally relevant to male and female subjects. Increased molecular profiles in female subjects are supported by our prior studies of SST in cingulate and DLPFC26,27, and BDNF-related profiles were not observed in the amygdala of male depressed subjects.34 Biological differences in disease processes or clinical specificities of that cohort (familial MDD) may underlie these results. Finally, it should be noted that numerous other genes were observed differentially expressed in this study (Table S2), including metabolic-, energetic-, vascular- and glial-related changes, which will need to be further investigated. We have focused on the impact of low BDNF on GABA markers, but our microarray study also confirmed genes previously reported affected in the context of stress regulation and/or mood disorders, such as NPY, TAC1 and RGS4. Indeed, individuals with low NPY display an increase in neural responsivity to negative stimuli and appear overrepresented in subjects with MDD.72,73 TAC1 has been associated with bipolar mood disorders,74 and the tachykininergic system may be a target for new therapeutic opportunities. Decreased RGS4 DLPFC levels have been reported in schizophrenia,75 and RGS4 polymorphisms appear associated with depression factors in this illness.76 Finally, while expression of SNAP25 seems higher in prefrontal cortex in schizophrenia,77 alterations in its levels have been reported in ventral hippocampus of depressed subjects,47 suggesting a disease-differential BDNF-activity regulation of this SNARE complex member.
Supplementary Material
1
2
3
4
Acknowledgements
This work was supported by National Institute of Mental Health (NIMH) MH084060 (ES), MH085111 (ES) and MH084053 (DAL). The funding agency had no role in the study design, data collection and analysis, decision to publish and preparation of the manuscript. The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIMH or the National Institutes of Health.
Footnotes
Conflict of Interest
David A. Lewis currently receives investigator-initiated research support from the BMS Foundation, Bristol-Myers Squibb, Curridium Ltd and Pfizer and in 2008–2010 served as a consultant in the areas of target identification and validation and new compound development to AstraZeneca, BioLine RX, Bristol-Myers Squibb, Merck, Neurogen and SK Life Science.
References
- 1.Cotter D, Mackay D, Landau S, Kerwin R, Everall I. Reduced glial cell density and neuronal size in the anterior cingulate cortex in major depressive disorder. Arch Gen Psychiatry. 2001;58(6):545–553. doi: 10.1001/archpsyc.58.6.545. [DOI] [PubMed] [Google Scholar]
- 2.Rajkowska G, Halaris A, Selemon LD. Reductions in neuronal and glial density characterize the dorsolateral prefrontal cortex in bipolar disorder. Biol Psychiatry. 2001;49(9):741–752. doi: 10.1016/s0006-3223(01)01080-0. [DOI] [PubMed] [Google Scholar]
- 3.Bowley MP, Drevets WC, Ongur D, Price JL. Low glial numbers in the amygdala in major depressive disorder. Biol Psychiatry. 2002;52(5):404–412. doi: 10.1016/s0006-3223(02)01404-x. [DOI] [PubMed] [Google Scholar]
- 4.Hamidi M, Drevets WC, Price JL. Glial reduction in amygdala in major depressive disorder is due to oligodendrocytes. Biol Psychiatry. 2004;55(6):563–569. doi: 10.1016/j.biopsych.2003.11.006. [DOI] [PubMed] [Google Scholar]
- 5.Videbech P, Ravnkilde B. Hippocampal volume and depression: a meta-analysis of MRI studies. Am J Psychiatry. 2004;161(11):1957–1966. doi: 10.1176/appi.ajp.161.11.1957. [DOI] [PubMed] [Google Scholar]
- 6.Campbell S, Marriott M, Nahmias C, MacQueen GM. Lower hippocampal volume in patients suffering from depression: a meta-analysis. Am J Psychiatry. 2004;161(4):598–607. doi: 10.1176/appi.ajp.161.4.598. [DOI] [PubMed] [Google Scholar]
- 7.Shirayama Y, Chen AC, Nakagawa S, Russell DS, Duman RS. Brain-derived neurotrophic factor produces antidepressant effects in behavioral models of depression. J Neurosci. 2002;22(8):3251–3261. doi: 10.1523/JNEUROSCI.22-08-03251.2002. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 8.Deltheil T, Tanaka K, Reperant C, Hen R, David DJ, Gardier AM. Synergistic neurochemical and behavioural effects of acute intrahippocampal injection of brain-derived neurotrophic factor and antidepressants in adult mice. Int J Neuropsychopharmacol. 2009:1–11. doi: 10.1017/S1461145709000017. [DOI] [PubMed] [Google Scholar]
- 9.Monteggia LM, Barrot M, Powell CM, Berton O, Galanis V, Gemelli T, et al. Essential role of brain-derived neurotrophic factor in adult hippocampal function. Proc Natl Acad Sci U S A. 2004;101(29):10827–10832. doi: 10.1073/pnas.0402141101. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 10.Saarelainen T, Hendolin P, Lucas G, Koponen E, Sairanen M, MacDonald E, et al. Activation of the TrkB neurotrophin receptor is induced by antidepressant drugs and is required for antidepressant-induced behavioral effects. J Neurosci. 2003;23(1):349–357. doi: 10.1523/JNEUROSCI.23-01-00349.2003. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 11.Molendijk ML, Bus BA, Spinhoven P, Penninx BW, Kenis G, Prickaerts J, et al. Serum levels of brain-derived neurotrophic factor in major depressive disorder: state-trait issues, clinical features and pharmacological treatment. Mol Psychiatry. 2010 doi: 10.1038/mp.2010.98. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 12.Sen S, Duman R, Sanacora G. Serum brain-derived neurotrophic factor, depression, and antidepressant medications: meta-analyses and implications. Biol Psychiatry. 2008;64(6):527–532. doi: 10.1016/j.biopsych.2008.05.005. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 13.Dunham JS, Deakin JF, Miyajima F, Payton A, Toro CT. Expression of hippocampal brain-derived neurotrophic factor and its receptors in Stanley consortium brains. J Psychiatr Res. 2009;43(14):1175–1184. doi: 10.1016/j.jpsychires.2009.03.008. [DOI] [PubMed] [Google Scholar]
- 14.Thompson Ray M, Weickert CS, Wyatt E, Webster MJ. Decreased BDNF, trkB-TK+ and GAD(67) mRNA expression in the hippocampus of individuals with schizophrenia and mood disorders. J Psychiatry Neurosci. 2011;36(1):100048. doi: 10.1503/jpn.100048. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 15.Chen B, Dowlatshahi D, MacQueen GM, Wang JF, Young LT. Increased hippocampal BDNF immunoreactivity in subjects treated with antidepressant medication. Biol Psychiatry. 2001;50(4):260–265. doi: 10.1016/s0006-3223(01)01083-6. [DOI] [PubMed] [Google Scholar]
- 16.Dwivedi Y, Rizavi HS, Conley RR, Roberts RC, Tamminga CA, Pandey GN. Altered gene expression of brain-derived neurotrophic factor and receptor tyrosine kinase B in postmortem brain of suicide subjects. Arch Gen Psychiatry. 2003;60(8):804–815. doi: 10.1001/archpsyc.60.8.804. [DOI] [PubMed] [Google Scholar]
- 17.Karege F, Vaudan G, Schwald M, Perroud N, La Harpe R. Neurotrophin levels in postmortem brains of suicide victims and the effects of antemortem diagnosis and psychotropic drugs. Brain Res Mol Brain Res. 2005;136(1–2):29–37. doi: 10.1016/j.molbrainres.2004.12.020. [DOI] [PubMed] [Google Scholar]
- 18.Kozicz T, Tilburg-Ouwens D, Faludi G, Palkovits M, Roubos E. Gender-related urocortin 1 and brain-derived neurotrophic factor expression in the adult human midbrain of suicide victims with major depression. Neuroscience. 2008;152(4):1015–1023. doi: 10.1016/j.neuroscience.2007.12.050. [DOI] [PubMed] [Google Scholar]
- 19.Ernst C, Deleva V, Deng X, Sequeira A, Pomarenski A, Klempan T, et al. Alternative splicing, methylation state, and expression profile of tropomyosin-related kinase B in the frontal cortex of suicide completers. Arch Gen Psychiatry. 2009;66(1):22–32. doi: 10.1001/archpsyc.66.1.22. [DOI] [PubMed] [Google Scholar]
- 20.Keller S, Sarchiapone M, Zarrilli F, Videtic A, Ferraro A, Carli V, et al. Increased BDNF promoter methylation in the Wernicke area of suicide subjects. Arch Gen Psychiatry. 2010;67(3):258–267. doi: 10.1001/archgenpsychiatry.2010.9. [DOI] [PubMed] [Google Scholar]
- 21.Pregelj P, Nedic G, Paska AV, Zupanc T, Nikolac M, Balazic J, et al. The association between brain-derived neurotrophic factor polymorphism (BDNF Val66Met) and suicide. J Affect Disord. 2011;128(3):287–290. doi: 10.1016/j.jad.2010.07.001. [DOI] [PubMed] [Google Scholar]
- 22.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(8):e6585. doi: 10.1371/journal.pone.0006585. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 23.Sanacora G, Gueorguieva R, Epperson CN, Wu YT, Appel M, Rothman DL, et al. Subtype-specific alterations of gamma-aminobutyric acid and glutamate in patients with major depression. Arch Gen Psychiatry. 2004;61(7):705–713. doi: 10.1001/archpsyc.61.7.705. [DOI] [PubMed] [Google Scholar]
- 24.Sanacora G, Mason GF, Rothman DL, Behar KL, Hyder F, Petroff OA, et al. Reduced cortical gamma-aminobutyric acid levels in depressed patients determined by proton magnetic resonance spectroscopy. Arch Gen Psychiatry. 1999;56(11):1043–1047. doi: 10.1001/archpsyc.56.11.1043. [DOI] [PubMed] [Google Scholar]
- 25.Levinson AJ, Fitzgerald PB, Favalli G, Blumberger DM, Daigle M, Daskalakis ZJ. Evidence of cortical inhibitory deficits in major depressive disorder. Biol Psychiatry. 2010;67(5):458–464. doi: 10.1016/j.biopsych.2009.09.025. [DOI] [PubMed] [Google Scholar]
- 26.Sibille E, Morris HM, Kota RS, Lewis DA. GABA-related transcripts in the dorsolateral prefrontal cortex in mood disorders. Int J Neuropsychopharmacol. 2011:1–14. doi: 10.1017/S1461145710001616. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 27.Tripp A, Kota RS, Lewis DA, Sibille E. Reduced somatostatin in subgenual anterior cingulate cortex in major depression. Neurobiol Dis. 2011;42(1):116–124. doi: 10.1016/j.nbd.2011.01.014. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 28.Glorioso C, Sabatini M, Unger T, Hashimoto T, Monteggia LM, Lewis DA, et al. Specificity and timing of neocortical transcriptome changes in response to BDNF gene ablation during embryogenesis or adulthood. Mol Psychiatry. 2006;11(7):633–648. doi: 10.1038/sj.mp.4001835. [DOI] [PubMed] [Google Scholar]
- 29.Muller JF, Mascagni F, McDonald AJ. Postsynaptic targets of somatostatin-containing interneurons in the rat basolateral amygdala. J Comp Neurol. 2007;500(3):513–529. doi: 10.1002/cne.21185. [DOI] [PubMed] [Google Scholar]
- 30.Seminowicz DA, Mayberg HS, McIntosh AR, Goldapple K, Kennedy S, Segal Z, et al. Limbic-frontal circuitry in major depression: a path modeling metanalysis. Neuroimage. 2004;22(1):409–418. doi: 10.1016/j.neuroimage.2004.01.015. [DOI] [PubMed] [Google Scholar]
- 31.Pezawas L, Meyer-Lindenberg A, Goldman AL, Verchinski BA, Chen G, Kolachana BS, et al. Evidence of biologic epistasis between BDNF and SLC6A4 and implications for depression. Mol Psychiatry. 2008;13(7):709–716. doi: 10.1038/mp.2008.32. [DOI] [PubMed] [Google Scholar]
- 32.Hamilton JP, Gotlib IH. Neural substrates of increased memory sensitivity for negative stimuli in major depression. Biol Psychiatry. 2008;63(12):1155–1162. doi: 10.1016/j.biopsych.2007.12.015. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 33.Glantz LA, Lewis DA. Reduction of synaptophysin immunoreactivity in the prefrontal cortex of subjects with schizophrenia. Regional and diagnostic specificity. Arch Gen Psychiatry. 1997;54(10):943–952. doi: 10.1001/archpsyc.1997.01830220065010. [DOI] [PubMed] [Google Scholar]
- 34.Sibille E, Wang Y, Joeyen-Waldorf J, Gaiteri C, Surget A, Oh S, et al. A molecular signature of depression in the amygdala. Am J Psychiatry. 2009;166(9):1011–1024. doi: 10.1176/appi.ajp.2009.08121760. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 35.Benjamini Y, Hochberg Y. Controlling the False Discovery Rate - a Practical and Powerful Approach to Multiple Testing. J Roy Stat Soc B Met. 1995;57(1):289–300. [Google Scholar]
- 36.Morris HM, Hashimoto T, Lewis DA. Alterations in somatostatin mRNA expression in the dorsolateral prefrontal cortex of subjects with schizophrenia or schizoaffective disorder. Cereb Cortex. 2008;18(7):1575–1587. doi: 10.1093/cercor/bhm186. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 37.Korte M, Carroll P, Wolf E, Brem G, Thoenen H, Bonhoeffer T. Hippocampal long-term potentiation is impaired in mice lacking brain-derived neurotrophic factor. Proc Natl Acad Sci U S A. 1995;92(19):8856–8860. doi: 10.1073/pnas.92.19.8856. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 38.Sakata K, Woo NH, Martinowich K, Greene JS, Schloesser RJ, Shen L, et al. Critical role of promoter IV-driven BDNF transcription in GABAergic transmission and synaptic plasticity in the prefrontal cortex. Proc Natl Acad Sci U S A. 2009;106(14):5942–5947. doi: 10.1073/pnas.0811431106. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 39.Franklin KBJ, Paxinos G. The mouse brain in stereotaxic coordinates. 3rd edn. New York: Academic Press; 2008. p. xxv. 1 v. (various pagings)pp. [Google Scholar]
- 40.Sibille E, Arango V, Joeyen-Waldorf J, Wang Y, Leman S, Surget A, et al. Large-scale estimates of cellular origins of mRNAs: enhancing the yield of transcriptome analyses. J Neurosci Methods. 2008;167(2):198–206. doi: 10.1016/j.jneumeth.2007.08.009. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 41.Bloom FE, Morales M. The central 5-HT3 receptor in CNS disorders. Neurochem Res. 1998;23(5):653–659. doi: 10.1023/a:1022486705184. [DOI] [PubMed] [Google Scholar]
- 42.Pruunsild P, Kazantseva A, Aid T, Palm K, Timmusk T. Dissecting the human BDNF locus: bidirectional transcription, complex splicing, and multiple promoters. Genomics. 2007;90(3):397–406. doi: 10.1016/j.ygeno.2007.05.004. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 43.Erraji-BenChekroun L, Underwood MD, Arango V, Galfalvy HC, Pavlidis P, Smyrniotopoulos P, et al. Molecular aging in human prefrontal cortex is selective and continuous throughout adult life. Biological Psychiatry. 2005;57:549–558. doi: 10.1016/j.biopsych.2004.10.034. [DOI] [PubMed] [Google Scholar]
- 44.Guilloux JP, Seney M, Edgar N, Sibille E. Integrated behavioral z-scoring increases the sensitivity and reliability of behavioral phenotyping in mice: Relevance to emotionality and sex. J Neurosci Methods. 2011 doi: 10.1016/j.jneumeth.2011.01.019. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 45.van Disseldorp J, Faddy MJ, Themmen AP, de Jong FH, Peeters PH, van der Schouw YT, et al. Relationship of serum antimullerian hormone concentration to age at menopause. J Clin Endocrinol Metab. 2008;93(6):2129–2134. doi: 10.1210/jc.2007-2093. [DOI] [PubMed] [Google Scholar]
- 46.Takei N, Sasaoka K, Inoue K, Takahashi M, Endo Y, Hatanaka H. Brain-derived neurotrophic factor increases the stimulation-evoked release of glutamate and the levels of exocytosis-associated proteins in cultured cortical neurons from embryonic rats. J Neurochem. 1997;68(1):370–375. doi: 10.1046/j.1471-4159.1997.68010370.x. [DOI] [PubMed] [Google Scholar]
- 47.Fatemi SH, Earle JA, Stary JM, Lee S, Sedgewick J. Altered levels of the synaptosomal associated protein SNAP-25 in hippocampus of subjects with mood disorders and schizophrenia. Neuroreport. 2001;12(15):3257–3262. doi: 10.1097/00001756-200110290-00023. [DOI] [PubMed] [Google Scholar]
- 48.Lewis DA, Campbell MJ, Morrison JH. An immunohistochemical characterization of somatostatin-28 and somatostatin-281-12 in monkey prefrontal cortex. J Comp Neurol. 1986;248(1):1–18. doi: 10.1002/cne.902480102. [DOI] [PubMed] [Google Scholar]
- 49.Glorioso C, Oh S, Douillard GG, Sibille E. Brain molecular aging, promotion of neurological disease and modulation by Sirtuin5 longevity gene polymorphism. Neurobiol Dis. 2011;41(2):279–290. doi: 10.1016/j.nbd.2010.09.016. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 50.Russo-Neustadt A, Beard RC, Cotman CW. Exercise, antidepressant medications, and enhanced brain derived neurotrophic factor expression. Neuropsychopharmacology. 1999;21(5):679–682. doi: 10.1016/S0893-133X(99)00059-7. [DOI] [PubMed] [Google Scholar]
- 51.Duman RS, Monteggia LM. A neurotrophic model for stress-related mood disorders. Biol Psychiatry. 2006;59(12):1116–1127. doi: 10.1016/j.biopsych.2006.02.013. [DOI] [PubMed] [Google Scholar]
- 52.Mattson MP, Maudsley S, Martin B. A neural signaling triumvirate that influences ageing and age-related disease: insulin/IGF-1, BDNF and serotonin. Ageing Res Rev. 2004;3(4):445–464. doi: 10.1016/j.arr.2004.08.001. [DOI] [PubMed] [Google Scholar]
- 53.Buckley PF, Mahadik S, Pillai A, Terry A., Jr Neurotrophins and schizophrenia. Schizophr Res. 2007;94(1–3):1–11. doi: 10.1016/j.schres.2007.01.025. [DOI] [PubMed] [Google Scholar]
- 54.Zuccato C, Cattaneo E. Brain-derived neurotrophic factor in neurodegenerative diseases. Nat Rev Neurol. 2009;5(6):311–322. doi: 10.1038/nrneurol.2009.54. [DOI] [PubMed] [Google Scholar]
- 55.Gahete MD, Rubio A, Duran-Prado M, Avila J, Luque RM, Castano JP. Expression of Somatostatin, cortistatin, and their receptors, as well as dopamine receptors, but not of neprilysin, are reduced in the temporal lobe of Alzheimer's disease patients. J Alzheimers Dis. 2010;20(2):465–475. doi: 10.3233/JAD-2010-1385. [DOI] [PubMed] [Google Scholar]
- 56.Fino E, Yuste R. Dense inhibitory connectivity in neocortex. Neuron. 2011;69(6):1188–1203. doi: 10.1016/j.neuron.2011.02.025. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 57.Hashimoto T, Volk DW, Eggan SM, Mirnics K, Pierri JN, Sun Z, et al. Gene expression deficits in a subclass of GABA neurons in the prefrontal cortex of subjects with schizophrenia. J Neurosci. 2003;23(15):6315–6326. doi: 10.1523/JNEUROSCI.23-15-06315.2003. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 58.Siegle GJ, Thompson W, Carter CS, Steinhauer SR, Thase ME. Increased amygdala and decreased dorsolateral prefrontal BOLD responses in unipolar depression: related and independent features. Biol Psychiatry. 2007;61(2):198–209. doi: 10.1016/j.biopsych.2006.05.048. [DOI] [PubMed] [Google Scholar]
- 59.Tsankova NM, Berton O, Renthal W, Kumar A, Neve RL, Nestler EJ. Sustained hippocampal chromatin regulation in a mouse model of depression and antidepressant action. Nat Neurosci. 2006;9(4):519–525. doi: 10.1038/nn1659. [DOI] [PubMed] [Google Scholar]
- 60.Surget A, Wang Y, Leman S, Ibarguen-Vargas Y, Edgar N, Griebel G, et al. Corticolimbic transcriptome changes are state-dependent and region-specific in a rodent model of depression and of antidepressant reversal. Neuropsychopharmacology. 2009;34(6):1363–1380. doi: 10.1038/npp.2008.76. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 61.Chourbaji S, Hellweg R, Brandis D, Zorner B, Zacher C, Lang UE, et al. Mice with reduced brain-derived neurotrophic factor expression show decreased choline acetyltransferase activity, but regular brain monoamine levels and unaltered emotional behavior. Brain Res Mol Brain Res. 2004;121(1–2):28–36. doi: 10.1016/j.molbrainres.2003.11.002. [DOI] [PubMed] [Google Scholar]
- 62.Ibarguen-Vargas Y, Surget A, Vourc'h P, Leman S, Andres CR, Gardier AM, et al. Deficit in BDNF does not increase vulnerability to stress but dampens antidepressant-like effects in the unpredictable chronic mild stress. Behav Brain Res. 2009;202(2):245–251. doi: 10.1016/j.bbr.2009.03.040. [DOI] [PubMed] [Google Scholar]
- 63.Berton O, McClung CA, DiLeone RJ, Krishnan V, Renthal W, Russo SJ, et al. Essential role of BDNF in the mesolimbic dopamine pathway in social defeat stress. Science. 2006;311(5762):864–868. doi: 10.1126/science.1120972. [DOI] [PubMed] [Google Scholar]
- 64.Autry AE, Adachi M, Cheng P, Monteggia LM. Gender-specific impact of brain-derived neurotrophic factor signaling on stress-induced depression-like behavior. Biol Psychiatry. 2009;66(1):84–90. doi: 10.1016/j.biopsych.2009.02.007. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 65.Monteggia LM, Luikart B, Barrot M, Theobold D, Malkovska I, Nef S, et al. Brain-derived neurotrophic factor conditional knockouts show gender differences in depression-related behaviors. Biol Psychiatry. 2007;61(2):187–197. doi: 10.1016/j.biopsych.2006.03.021. [DOI] [PubMed] [Google Scholar]
- 66.Advani T, Koek W, Hensler JG. Gender differences in the enhanced vulnerability of BDNF+/− mice to mild stress. Int J Neuropsychopharmacol. 2009;12(5):583–588. doi: 10.1017/S1461145709000248. [DOI] [PubMed] [Google Scholar]
- 67.Sakata K, Jin L, Jha S. Lack of promoter IV-driven BDNF transcription results in depression-like behavior. Genes Brain Behav. 2010;9(7):712–721. doi: 10.1111/j.1601-183X.2010.00605.x. [DOI] [PubMed] [Google Scholar]
- 68.Zeyda T, Diehl N, Paylor R, Brennan MB, Hochgeschwender U. Impairment in motor learning of somatostatin null mutant mice. Brain Res. 2001;906(1–2):107–114. doi: 10.1016/s0006-8993(01)02563-x. [DOI] [PubMed] [Google Scholar]
- 69.Engin E, Stellbrink J, Treit D, Dickson CT. Anxiolytic and antidepressant effects of intracerebroventricularly administered somatostatin: behavioral and neurophysiological evidence. Neuroscience. 2008;157(3):666–676. doi: 10.1016/j.neuroscience.2008.09.037. [DOI] [PubMed] [Google Scholar]
- 70.Gatt JM, Nemeroff CB, Dobson-Stone C, Paul RH, Bryant RA, Schofield PR, et al. Interactions between BDNF Val66Met polymorphism and early life stress predict brain and arousal pathways to syndromal depression and anxiety. Mol Psychiatry. 2009;14(7):681–695. doi: 10.1038/mp.2008.143. [DOI] [PubMed] [Google Scholar]
- 71.Lau JY, Goldman D, Buzas B, Hodgkinson C, Leibenluft E, Nelson E, et al. BDNF gene polymorphism (Val66Met) predicts amygdala and anterior hippocampus responses to emotional faces in anxious and depressed adolescents. Neuroimage. 2010;53(3):952–961. doi: 10.1016/j.neuroimage.2009.11.026. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 72.Mickey BJ, Zhou Z, Heitzeg MM, Heinz E, Hodgkinson CA, Hsu DT, et al. Emotion processing, major depression, and functional genetic variation of neuropeptide Y. Arch Gen Psychiatry. 2011;68(2):158–166. doi: 10.1001/archgenpsychiatry.2010.197. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 73.Zhou Z, Zhu G, Hariri AR, Enoch MA, Scott D, Sinha R, et al. Genetic variation in human NPY expression affects stress response and emotion. Nature. 2008;452(7190):997–1001. doi: 10.1038/nature06858. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 74.Ogden CA, Rich ME, Schork NJ, Paulus MP, Geyer MA, Lohr JB, et al. Candidate genes, pathways and mechanisms for bipolar (manic-depressive) and related disorders: an expanded convergent functional genomics approach. Mol Psychiatry. 2004;9(11):1007–1029. doi: 10.1038/sj.mp.4001547. [DOI] [PubMed] [Google Scholar]
- 75.Mirnics K, Middleton FA, Stanwood GD, Lewis DA, Levitt P. Disease-specific changes in regulator of G-protein signaling 4 (RGS4) expression in schizophrenia. Mol Psychiatry. 2001;6(3):293–301. doi: 10.1038/sj.mp.4000866. [DOI] [PubMed] [Google Scholar]
- 76.Rethelyi JM, Bakker SC, Polgar P, Czobor P, Strengman E, Pasztor PI, et al. Association study of NRG1, DTNBP1, RGS4, G72/G30, and PIP5K2A with schizophrenia and symptom severity in a Hungarian sample. Am J Med Genet B Neuropsychiatr Genet. 2010;153B(3):792–801. doi: 10.1002/ajmg.b.31049. [DOI] [PubMed] [Google Scholar]
- 77.Etain B, Dumaine A, Mathieu F, Chevalier F, Henry C, Kahn JP, et al. A SNAP25 promoter variant is associated with early-onset bipolar disorder and a high expression level in brain. Mol Psychiatry. 2010;15(7):748–755. doi: 10.1038/mp.2008.148. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 78.Mai JK, Assheuer J, Paxinos G. Atlas of the human brain. 2nd edn. Amsterdam ; Boston: Elsevier Academic Press; 2004. pp. viii–246. [Google Scholar]
Associated Data
This section collects any data citations, data availability statements, or supplementary materials included in this article.
Supplementary Materials
1
2
3
4


