A Human-Mouse Conserved Sex Bias in Amygdala Gene Expression Related to Circadian Clock and Energy Metabolism Li-Chun Lin, David a Lewis and Etienne Sibille*
Total Page:16
File Type:pdf, Size:1020Kb
Lin et al. Molecular Brain 2011, 4:18 http://www.molecularbrain.com/content/4/1/18 SHORTREPORT Open Access A human-mouse conserved sex bias in amygdala gene expression related to circadian clock and energy metabolism Li-Chun Lin, David A Lewis and Etienne Sibille* Abstract Background: 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. Findings: 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. Conclusions: 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. Introduction amygdala complex improves microarray sensitivity to Males and females differ in behavior and brain structure, detect differences in gene expression by reducing sample as well as in prevalence of neuropsychiatric disorders. heterogeneity, and could facilitate the subsequent identi- The greater prevalence of major depression and rapid fication of genetic sex differences. To characterize sex cycling bipolar disorder in women [1,2] highlights the differences in the intrinsic molecular properties of the importance of studying potential mechanisms for sex amygdala, we examined the presence of sexually differences. Sex chromosome- and hormone-linked dimorphic patterns of gene expression in the basolateral genes may directly affect sexually dimorphic neurobiol- nuclei of the amygdala using postmortem samples from ogy [3]. For example, sex hormones regulate the size of control human subjects and mice. the medial amygdala and bed nucleus of the stria termi- nalis during development [4]. The lateral (LA) and basal Methods (basolateral; throughout the article we refer to human Human Brain Samples BA and mouse BLA) nuclei of the amygdala have been Postmortem amygdala samples from 12 healthy male identified as critical sites for learned fear and emotion controls and 12 healthy female controls (age = 39-64 regulation in healthy subjects and in patients with mood years) were obtained from the University of Pittsburgh disorders [5,6], but little is known about their sexual Brain Donation Program (Table 1), as previously dimorphism. The isolation of subnuclei from the described [6]. Male and female groups were matched in pairs as closely as possible on the basis of RNA integrity, * Correspondence: [email protected] pH,post-morteminterval,age,andsex,andthegroups Translational Neuroscience Program, Department of Psychiatry, and Center did not differ on any of these parameters (p > 0.05; ’ for Neuroscience, 3811 O Hara Street, University of Pittsburgh, Pittsburgh, PA Table 1). All subjects died suddenly without prolonged 15213, USA © 2011 Lin et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Lin et al. Molecular Brain 2011, 4:18 Page 2 of 6 http://www.molecularbrain.com/content/4/1/18 Table 1 Characteristics of matched human subjects Female Group Male Group Age Interval Age Interval Pair Case (Years) Race (Hours) pH RIN Case (Years) Race (Hours) pH RIN 1 HU1282 39 W 24.5 6.8 7.5 HU604 39 W 19.3 7.1 8.6 2 HU1092 40 B 16.6 6.8 8.0 HU1047 43 W 12.0 6.6 9.0 3 HU627a,b 43 B 14.1 7.1 7.0 HU857a,b 48 W 16.6 6.7 8.9 4 HU567 46 W 15.0 6.8 8.9 HU1067b 49 W 6.0 6.6 8.2 5 HU1280a 50 W 23.5 6.7 7.7 HU1086a,b 51 W 24.2 6.8 8.1 6 HU1391a,b 51 W 7.8 6.6 7.1 HU852a,b 54 W 8.0 6.8 9.1 7 HU686 52 W 22.6 7.1 8.5 HU1031b 53 W 23.2 6.8 8.9 8 HU575 55 B 11.3 6.8 9.6 HU1122 55 W 15.4 6.7 7.9 9 HU1247 58 W 22.7 6.4 8.4 HU685 56 W 14.0 6.6 8.0 10 HU1318 58 W 18.8 6.7 7.4 HU852 54 W 8.0 6.8 9.1 11 HU568 60 W 9.5 6.9 8.7 HU551 61 W 16.4 6.6 8.3 12 HU1466b 64 B 20.0 6.7 8.8 HU615 62 W 7.2 6.4 7.8 Matched samples processed for microarray: AVG 51.3 17.2 6.8 8.1 52.1 14.2 6.7 8.5 STDEV 8.1 5.7 0.2 0.8 6.7 6.2 0.2 0.5 Matched LA for final array analysis: AVG 52.4 17.9 6.8 8.4 52.4 13.5 6.7 8.4 STDEV 8.9 5.2 0.2 0.7 7.7 5.8 0.2 0.5 Matched BA for final array analysis: AVG 50.9 18.3 6.8 8.3 52.9 13.2 6.7 8.4 STDEV 7.8 5.5 0.2 0.7 8.7 4.4 0.2 0.5 Samples from LA (marked a) and BA (marked b) processed but excluded due to lower microarray quality control; W = Caucasian; B = African American; RIN = RNA integrity number; AVG = average; STDVE = standard deviation. agonal periods. See [6] for additional details. In contrast LMD 6500 laser microdissection system. All mice were to the light acetylcholinesterase (AChE)-stained LA maintained on a 12-hour light-dark cycle (lights on across species, human BA and mouse BLA contained 06:00 and lights off on 18:00) with ad libitum food and the heavy stained AChE-positive neuropil in the amyg- water. All procedures were approved by the University dala [7,8], implying that human BA and LA anatomically of Pittsburgh Institutional Animal Care and Use Com- correspond to the mouse BLA and LA according to the mittee (protocol # 0911014A-2, Animal Assurance # well-defined AChE boundaries. Hence, based on the dis- A3187-01). tribution of AChE activity and thionin staining, human lateral(LA)andbasal(BA)amygdalanucleiwere RNA extraction for Microarray Samples defined and separately harvested from 20-μmcryostat Total RNA was extracted using the RNeasy Plus Micro sections using a clean RNase-free pipette tip. All proce- Kit (QIAGEN, Valencia, CA) and processed for microar- dures were approved by the University of Pittsburgh’s ray analysis according to manufacturer’sprotocol.RNA Committee for the Oversight of Research Involving the samples were amplified with Illumina TotalPrep RNA Dead and Institutional Review Board for Biomedical Amplification Kit (for human tissue; Illumina, San Research. Diego, CA) and Nugen WT-Ovation Pico RNA Amplifi- cation System (for mouse tissue; NuGen, San Carlos, Mouse Brain Samples CA). The fragmented labeled cRNA samples were pro- Brains from adult C57BL/6J mice (n = 5 per sex; pooled cessed and hybridized to microarrays (Human HT-12v3 by 2 from 10 mice/group, 3 months of age) were har- and MouseWG-6v2.0 BeadChips). vested, fresh frozen and cryo-sectioned at 20 μmon UV-treated laser-capture microscopy slides (Leica Array Data Analysis Microsystems, Wetzlar, Germany). Sections on the slides For statistical analysis, the log2-transformed signal inten- were dried and stored at -80°C until used. Lateral (LA) sities of probe-sets that passed the initial expression to basolateral (BLA) amygdala nuclei were visualized filters (< 30% missing values, detection p-value > 0.1), with rapid thionin staining and dissected with the Leica were extracted and renormalized to the same mean and Lin et al. Molecular Brain 2011, 4:18 Page 3 of 6 http://www.molecularbrain.com/content/4/1/18 variance across arrays in order to control for any effects mouse LA (p< 6.72E-02). Assembly of RNA Polymerase of batch. Discovery threshold for male/female differ- III Complex was the only identified canonical pathway ences were applied in both species (two-tailed t-test p- in human LA and was at trend level (p < 0.1). value < 0.05; > 20% effect size). Functional analyses Within the mitochondrial-related functions, genes were performed on the identified gene list using the encoding enzymes of the mitochondrial electron trans- Ingenuity Pathways Knowledge Base [9] and compared port chain (ETC) consistently displayed high male-biased across species. expression (Figure 2A). Also, two genes coding for mito- chondria-localized sirtuins (SIRT3, SIRT5), which serve Real-time Quantitative PCR (qPCR) Data Analysis as metabolic regulators of circadian rhythms by utilizing Separate RNA samples were reverse-transcribed using NAD+ [10], also showed similar high male-biased expres- QScript cDNA super mix (QuantaBioscience, Gaithers- sion in human BA and mouse BLA. Of related interest burg, MD) using a combination of oligo (dT) and ran- was the identification of sex-biased circadian rhythm sig- dom primers. Each PCR reaction was performed in naling in human BA (IPA, p< 5.39E-04)andmouseLA triplicate and compared to three internal controls (p<3.54E-02), which included conserved high female- (GAPDH, beta-actin and cyclophilin).