- Short report
- Open Access
A human-mouse conserved sex bias in amygdala gene expression related to circadian clock and energy metabolism
Molecular Brain volume 4, Article number: 18 (2011)
Major depression affects twice as many women as men, but the underlying molecular mechanisms responsible for the heightened female vulnerability are not known. The amygdala, composed of heterogeneous subnuclei, participates in multiple functional circuits regulating emotional responses to stress. We hypothesized that sex differences in molecular structure may contribute to differential mood regulation and disease vulnerability.
Using gene arrays followed by quantitative PCR validation, we compared the transcriptome profiles between sexes in human and mouse amygdala. We now report sexually dimorphic features of transcriptomes in the basolateral nucleus of the amygdala, and these features are highly conserved across species. A functional analysis of differential gene expression showed that mitochondrial-related gene groups were identified as the top biological pathways associated with sexual dimorphism in both species.
These results suggest that the basolateral amygdala is a sexually dimorphic structure, featuring a regulatory cascade of mitochondrial function and circadian rhythm, potentially linked through sirtuins and hormone nuclear receptors. Hence, baseline differences in amygdalar circadian regulation of cellular metabolism may contribute to sex-related differences in mood regulation and vulnerability to major depression.
Males and females differ in behavior and brain structure, as well as in prevalence of neuropsychiatric disorders. The greater prevalence of major depression and rapid cycling bipolar disorder in women [1, 2] highlights the importance of studying potential mechanisms for sex differences. Sex chromosome- and hormone-linked genes may directly affect sexually dimorphic neurobiology . For example, sex hormones regulate the size of the medial amygdala and bed nucleus of the stria terminalis during development . The lateral (LA) and basal (basolateral; throughout the article we refer to human BA and mouse BLA) nuclei of the amygdala have been identified as critical sites for learned fear and emotion regulation in healthy subjects and in patients with mood disorders [5, 6], but little is known about their sexual dimorphism. The isolation of subnuclei from the amygdala complex improves microarray sensitivity to detect differences in gene expression by reducing sample heterogeneity, and could facilitate the subsequent identification of genetic sex differences. To characterize sex differences in the intrinsic molecular properties of the amygdala, we examined the presence of sexually dimorphic patterns of gene expression in the basolateral nuclei of the amygdala using postmortem samples from control human subjects and mice.
Human Brain Samples
Postmortem amygdala samples from 12 healthy male controls and 12 healthy female controls (age = 39-64 years) were obtained from the University of Pittsburgh Brain Donation Program (Table 1), as previously described . Male and female groups were matched in pairs as closely as possible on the basis of RNA integrity, pH, post-mortem interval, age, and sex, and the groups did not differ on any of these parameters (p > 0.05; Table 1). All subjects died suddenly without prolonged agonal periods. See  for additional details. In contrast to the light acetylcholinesterase (AChE)-stained LA across species, human BA and mouse BLA contained the heavy stained AChE-positive neuropil in the amygdala [7, 8], implying that human BA and LA anatomically correspond to the mouse BLA and LA according to the well-defined AChE boundaries. Hence, based on the distribution of AChE activity and thionin staining, human lateral (LA) and basal (BA) amygdala nuclei were defined and separately harvested from 20-μm cryostat sections using a clean RNase-free pipette tip. All procedures were approved by the University of Pittsburgh's Committee for the Oversight of Research Involving the Dead and Institutional Review Board for Biomedical Research.
Mouse Brain Samples
Brains from adult C57BL/6J mice (n = 5 per sex; pooled by 2 from 10 mice/group, 3 months of age) were harvested, fresh frozen and cryo-sectioned at 20 μm on UV-treated laser-capture microscopy slides (Leica Microsystems, Wetzlar, Germany). Sections on the slides were dried and stored at -80°C until used. Lateral (LA) to basolateral (BLA) amygdala nuclei were visualized with rapid thionin staining and dissected with the Leica LMD 6500 laser microdissection system. All mice were maintained on a 12-hour light-dark cycle (lights on 06:00 and lights off on 18:00) with ad libitum food and water. All procedures were approved by the University of Pittsburgh Institutional Animal Care and Use Committee (protocol # 0911014A-2, Animal Assurance # A3187-01).
RNA extraction for Microarray Samples
Total RNA was extracted using the RNeasy Plus Micro Kit (QIAGEN, Valencia, CA) and processed for microarray analysis according to manufacturer's protocol. RNA samples were amplified with Illumina TotalPrep RNA Amplification Kit (for human tissue; Illumina, San Diego, CA) and Nugen WT-Ovation Pico RNA Amplification System (for mouse tissue; NuGen, San Carlos, CA). The fragmented labeled cRNA samples were processed and hybridized to microarrays (Human HT-12v3 and MouseWG-6v2.0 BeadChips).
Array Data Analysis
For statistical analysis, the log2-transformed signal intensities of probe-sets that passed the initial expression filters (< 30% missing values, detection p-value > 0.1), were extracted and renormalized to the same mean and variance across arrays in order to control for any effects of batch. Discovery threshold for male/female differences were applied in both species (two-tailed t-test p-value < 0.05; > 20% effect size). Functional analyses were performed on the identified gene list using the Ingenuity Pathways Knowledge Base  and compared across species.
Real-time Quantitative PCR (qPCR) Data Analysis
Separate RNA samples were reverse-transcribed using QScript cDNA super mix (QuantaBioscience, Gaithersburg, MD) using a combination of oligo (dT) and random primers. Each PCR reaction was performed in triplicate and compared to three internal controls (GAPDH, beta-actin and cyclophilin). All designed primer sequences were designed with Primer3Plus (Additional File 1, table S1). Results were calculated as the geometric mean of relative intensities compared to the three internal controls. Validation of the array results was considered at p< 0.05, using one-tailed unpaired t-test.
Using exploratory criteria we generated large and informative parallel human/mouse transcriptome datasets from male and female BA (BLA in mice) and LA nuclei of the amygdala. As an internal control to assess the sensitivity of the array approach, all detected Y-chromosome genes (male-specific) displayed significantly male-biased expression, and expression of numerous X-linked genes showed a female bias (Figure 1). In humans, 1335 sex-biased autosomal genes were identified in BA (Additional File 2, table S2), and fewer (n = 165) in LA (Additional File 3, table S3). In mice, ~1% of autosomal genes displayed sex-biased expression (BLA, n = 515, Additional File 4, table S4; LA, n = 556, Additional File 5, table S5).
To reduce false discovery and focus on the most relevant changes, we characterized the functional and cross-species relatedness of autosomal sex-biased gene sets using the Ingenuity Pathways Analysis (IPA) database. Mitochondrial-related functions were consistently identified among the top four most-represented pathways in human BA (p< 1.57E-06), mouse BA (p< 8.37E-03) and mouse LA (p< 6.72E-02). Assembly of RNA Polymerase III Complex was the only identified canonical pathway in human LA and was at trend level (p < 0.1).
Within the mitochondrial-related functions, genes encoding enzymes of the mitochondrial electron transport chain (ETC) consistently displayed high male-biased expression (Figure 2A). Also, two genes coding for mitochondria-localized sirtuins (SIRT3, SIRT5), which serve as metabolic regulators of circadian rhythms by utilizing NAD+, also showed similar high male-biased expression in human BA and mouse BLA. Of related interest was the identification of sex-biased circadian rhythm signaling in human BA (IPA, p< 5.39E-04) and mouse LA (p< 3.54E-02), which included conserved high female-biased expression of selected circadian clock genes (e.g., Period 2 and trend-level for CLOCK; Figure 2A). Intriguingly, as nuclear receptors link the circadian clock to metabolism , retinoic acid orphan receptor β (RORB) and thyroid hormone receptor α (THRA) displayed high female-biased expression in human BA (RORB: average log ratio [alr] = 0.29, p < 0.026; THRA: alr = 0.31, p < 0.013).
The accuracy and reproducibility of the approach was assessed by qPCR on independent samples. 12 genes out of 18 tested were individually confirmed (p< 0.05; same direction; Figure 2B). Overall, the qPCR validation results were highly correlated with array data using Pearson tests for significance (r = 0.86, p < 2.5E-06; Figure 2D).
We next evaluated the anatomical extent of the observed sex-biased gene profiles in other brain areas using existing human array data in cortical regions . Numerous ETC and SIRT genes in dorsolateral prefrontal cortex (Brodman area 9 [BA9]; Figure 2C), and fewer in the ventrolateral prefrontal cortex (Brodman area 47 [BA47]; Figure 2C) and anterior cingulate cortex (ACC; Figure 2C), consistently displayed high male-biased expression, suggesting variable levels of regional conservation and amygdala enrichment.
Notably, within the limited cohort size and age range, we did not observe statistically significant age-related effects for genes confirmed by qPCR, except for male-biased ATP synthase genes (ATP5A1, r = -0.78, p < 0.005; ATP6VOD1, r = -0.58, p < 0.05) in human male amygdala. Since no information was available for menopause or estrogen replacement treatment for the human female subjects, we interpreted these results as indicative of sexual dimorphism independently of age or estrogen status, although this will need to be tested in larger cohorts.
We report sexually dimorphic gene expression profiles in the BA nucleus of the amygdala in human subjects. Changes at the level of single genes extend to homologous BLA nucleus in the mouse amygdala, albeit to a lesser extent. At the functional level, mitochondrial-related gene groups were consistently identified as the top biological pathways associated with male/female differences in both species. Together with prominent differences in other gene groups (Figure 2), our results suggest a sexual dimorphism in the expression level of an evolutionarily conserved pathway - circadian clock and mitochondrial function, potentially linked through sirtuin and nuclear hormone receptors.
In experimental models, circadian genes modulate energy metabolism by setting the periodic oscillations of expression and activity of enzymes to influence physiological and behavioral functions. Substantial evidence has shown that patients with major depression or bipolar disorder often display abnormalities in circadian rhythm, such as disturbances in the sleep/wake cycle and diurnal mood changes. Genetics studies have found associations with CLOCK, PER1-3, and CRY1, and Sirtuin genes in mood disorders [13–16]. In molecular/genetic studies, mitochondrial dysfunction has been associated with mood disorders [17, 18], and mood disorders are common in patients with mitochondrial disorders . Thus, our findings suggest that baseline sexual dimorphism in mitochondrial density or functions may influence energy availability in the amygdala. Notably, while mitochondrial function was systemically identified as a top sex-biased gene family across species, the species-specific orthologous genes were not necessarily and systematically the same, but very close functionally (Figure 2E). Because of identification of mitochondrial-related genes with closely related functions and locations across species, we speculate that the sex-biased mitochondria profile may be under conserved transcription control to display similar sexual dimorphism across species.
Limitations of this study are noted. The results will need to be confirmed at the protein and functional levels in independent and larger cohorts. Similarly, the potential contribution of race-specific amygdala gene expression will need to be assessed in larger cohorts. It is also unclear whether our findings in mid-life adults can be generalized to adolescent and elderly populations. Here we identified age effects in 2 of the 18 genes validated by qPCR, suggesting that aging may not contribute to the evolutionary conserved genetic sex differences in amygdala. In addition, this study sampled 10 intact female mice and 12 female human subjects from 39-64 years old, which should considerably lower the impact of variable estrogen levels on our findings of genetic sex difference, although we cannot exclude residual effects.
In conclusion, we present evidence for significant male/female differences in amygdala gene expression across species. Delineating the contribution of mitochondrial and circadian signaling in mood disorders will require additional work. For instance, sex-biased gene expression of circadian clock, nuclear hormone receptors, mitochondrial sirtuins and electron chain enzymes in an interlocking-feedback loop, may functionally translate into differential amygdala responsiveness (Figure 2E). This loop may lead to sex differences in metabolic capacity to meet energy demand, and might potentially contribute to the greater female vulnerability to mood disorders.
anterior cingulate cortex
Brodman area 9
Brodman area 47
electron transport chain
Ingenuity Pathways Analysis
Real-time Quantitative PCR
retinoic acid orphan receptor β
thyroid hormone receptor α.
Kendler KS: Gender differences in the genetic epidemiology of major depression. J Gend Specif Med. 1998, 1 (2): 28-31.
Tondo L, Baldessarini RJ: Rapid cycling in women and men with bipolar manic-depressive disorders. The American journal of psychiatry. 1998, 155 (10): 1434-1436.
Jazin E, Cahill L: Sex differences in molecular neuroscience: from fruit flies to humans. Nature reviews. 11 (1): 9-17.
Hines M, Allen LS, Gorski RA: Sex differences in subregions of the medial nucleus of the amygdala and the bed nucleus of the stria terminalis of the rat. Brain research. 1992, 579 (2): 321-326. 10.1016/0006-8993(92)90068-K.
Sheline YI, Barch DM, Donnelly JM, Ollinger JM, Snyder AZ, Mintun MA: Increased amygdala response to masked emotional faces in depressed subjects resolves with antidepressant treatment: an fMRI study. Biological psychiatry. 2001, 50 (9): 651-658. 10.1016/S0006-3223(01)01263-X.
Sibille E, Wang Y, Joeyen-Waldorf J, Gaiteri C, Surget A, Oh S, Belzung C, Tseng GC, Lewis DA: A molecular signature of depression in the amygdala. The American journal of psychiatry. 2009, 166 (9): 1011-1024. 10.1176/appi.ajp.2009.08121760.
Svendsen CN, Bird ED: Acetylcholinesterase staining of the human amygdala. Neuroscience letters. 1985, 54 (2-3): 313-318. 10.1016/S0304-3940(85)80097-5.
Kitt CA, Hohmann C, Coyle JT, Price DL: Cholinergic innervation of mouse forebrain structures. The Journal of comparative neurology. 1994, 341 (1): 117-129. 10.1002/cne.903410110.
The Ingenuity Pathways Knowledge Base. [http://www.ingenuity.com/]
Asher G, Gatfield D, Stratmann M, Reinke H, Dibner C, Kreppel F, Mostoslavsky R, Alt FW, Schibler U: SIRT1 regulates circadian clock gene expression through PER2 deacetylation. Cell. 2008, 134 (2): 317-328. 10.1016/j.cell.2008.06.050.
Yang X, Downes M, Yu RT, Bookout AL, He W, Straume M, Mangelsdorf DJ, Evans RM: Nuclear receptor expression links the circadian clock to metabolism. Cell. 2006, 126 (4): 801-810. 10.1016/j.cell.2006.06.050.
Galfalvy HC, Erraji-Benchekroun L, Smyrniotopoulos P, Pavlidis P, Ellis SP, Mann JJ, Sibille E, Arango V: Sex genes for genomic analysis in human brain: internal controls for comparison of probe level data extraction. BMC bioinformatics. 2003, 4: 37-10.1186/1471-2105-4-37.
Roybal K, Theobold D, Graham A, DiNieri JA, Russo SJ, Krishnan V, Chakravarty S, Peevey J, Oehrlein N, Birnbaum S, Vitaterna MH, Orsulak P, Takahashi JS, Nestler EJ, Carlezon WA, McClung CA: Mania-like behavior induced by disruption of CLOCK. Proceedings of the National Academy of Sciences of the United States of America. 2007, 104 (15): 6406-6411. 10.1073/pnas.0609625104.
Hampp G, Ripperger JA, Houben T, Schmutz I, Blex C, Perreau-Lenz S, Brunk I, Spanagel R, Ahnert-Hilger G, Meijer JH, Albrecht U: Regulation of monoamine oxidase A by circadian-clock components implies clock influence on mood. Curr Biol. 2008, 18 (9): 678-683. 10.1016/j.cub.2008.04.012.
Soria V, Martinez-Amoros E, Escaramis G, Valero J, Perez-Egea R, Garcia C, Gutierrez-Zotes A, Puigdemont D, Bayes M, Crespo JM, Martorell L, Vilella E, Labad A, Vallejo J, Pérez V, Menchón JM, Estivill X, Gratacòs M, Urretavizcaya M: Differential association of circadian genes with mood disorders: CRY1 and NPAS2 are associated with unipolar major depression and CLOCK and VIP with bipolar disorder. Neuropsychopharmacology. 35 (6): 1279-1289.
Abe N, Uchida S, Otsuki K, Hobara T, Yamagata H, Higuchi F, Shibata T, Watanabe Y: Altered sirtuin deacetylase gene expression in patients with a mood disorder. Journal of psychiatric research.
Konradi C, Eaton M, MacDonald ML, Walsh J, Benes FM, Heckers S: Molecular evidence for mitochondrial dysfunction in bipolar disorder. Archives of general psychiatry. 2004, 61 (3): 300-308. 10.1001/archpsyc.61.3.300.
Rollins B, Martin MV, Sequeira PA, Moon EA, Morgan LZ, Watson SJ, Schatzberg A, Akil H, Myers RM, Jones EG, Wallace DC, Bunney WE, Vawter MP: Mitochondrial variants in schizophrenia, bipolar disorder, and major depressive disorder. PloS one. 2009, 4 (3): e4913-10.1371/journal.pone.0004913.
Fattal O, Link J, Quinn K, Cohen BH, Franco K: Psychiatric comorbidity in 36 adults with mitochondrial cytopathies. CNS spectrums. 2007, 12 (6): 429-438.
We gratefully acknowledge support of this research from NIMH (MH077159, MH084060, MH085111), the University of Pittsburgh Brain Tissue Donation Program (MH084053), and Xingbin Wang for statistical support.
Dr. Lewis currently receives investigator-initiated research support from the BMS Foundation, Bristol-Myers Squibb, Curridium Ltd and Pfizer and in 2007-2010 served as a consultant in the areas of target identification and validation and new compound development to AstraZeneca, BioLine RX, Bristol-Myers Squibb, Hoffman-Roche, Lilly, Merck, Neurogen, and SK Life Science.
LCL and ES designed this study. LCL carried out the molecular experiments, performed the statistical analysis and drafted the manuscript. DAL provided technical and conceptual support for the human brain collection. The final manuscript was written by LCL and ES. All authors have read, edited and approved the final manuscript.