Stress amplifies sex differences in primate prefrontal profiles of gene expression
© The Author(s). 2017
Received: 25 July 2017
Accepted: 23 October 2017
Published: 2 November 2017
Stress is a recognized risk factor for mood and anxiety disorders that occur more often in women than men. Prefrontal brain regions mediate stress coping, cognitive control, and emotion. Here, we investigate sex differences and stress effects on prefrontal cortical profiles of gene expression in squirrel monkey adults.
Dorsolateral, ventrolateral, and ventromedial prefrontal cortical regions from 18 females and 12 males were collected after stress or no-stress treatment conditions. Gene expression profiles were acquired using HumanHT-12v4.0 Expression BeadChip arrays adapted for squirrel monkeys.
Extensive variation between prefrontal cortical regions was discerned in the expression of numerous autosomal and sex chromosome genes. Robust sex differences were also identified across prefrontal cortical regions in the expression of mostly autosomal genes. Genes with increased expression in females compared to males were overrepresented in mitogen-activated protein kinase and neurotrophin signaling pathways. Many fewer genes with increased expression in males compared to females were discerned, and no molecular pathways were identified. Effect sizes for sex differences were greater in stress compared to no-stress conditions for ventromedial and ventrolateral prefrontal cortical regions but not dorsolateral prefrontal cortex.
Stress amplifies sex differences in gene expression profiles for prefrontal cortical regions involved in stress coping and emotion regulation. Results suggest molecular targets for new treatments of stress disorders in human mental health.
Mood and anxiety disorders occur more often in women than men [1–4]. These mental health disorders are associated with stress-induced modifications in prefrontal cortex [5–7]. Neuroscience research has until recently focused on males , and sex differences in stress effects on primate prefrontal profiles of gene expression are not known. Differential gene expression profiles have been reported for diverse brain regions in humans [9, 10] and various monkey species [11, 12], but comparisons across prefrontal cortical regions in primates are uncommon. Prefrontal regions are of interest as elements in brain circuits that mediate stress coping, cognitive control, and emotion [13–16]. Human prefrontal neurobiology is modeled more often in nonhuman primates compared to other animals because of species differences in this region of brain .
Sex differences in gene expression profiles were initially thought to be limited in human brain . Recent evidence suggests that many genes may be differentially expressed between the sexes for various brain regions in mice [18–20], monkeys , and humans [9, 22, 23]. Certain sex differences may reflect species-specific adaptations , but sex differences in stress neurobiology commonly occur in many mammals  including squirrel monkeys . Sex differences and stress effects in squirrel monkey prefrontal cortical profiles of gene expression may provide insights on the biology of sex differences in human mental health.
Here, we examine prefrontal cortical gene expression profiles in female and male squirrel monkey adults randomly selected from a larger sample randomized to social separation stress or no-stress treatment conditions. Social separations consistently increase plasma levels of the stress hormone cortisol in female  and male  squirrel monkeys. Both sexes respond to separations from same-sex conspecifics with prolonged increases in cortisol, but females and males have not been directly compared. In separate studies, females show diminished neuroendocrine stress responses  and males show increased hippocampal neurogenesis  measured 2–4 weeks after exposure to same-sex social separations. Long-term effects are likewise assessed here for prefrontal cortical profiles of gene expression in both sexes. Gene expression profiling represents discovery research  that is not guided by hypotheses about specific genes. Analyses of gene expression profiles demand statistical sophistication , but we aim to balance statistics with biological implications in this report.
Twenty two male and 48 female squirrel monkey (Saimiri sciureus) adults served as subjects. These monkeys were approximately 6–17 years of age determined from birth records or estimated from dentition. All monkeys were housed in climate controlled rooms at ~ 26 °C on 12:12-h light/dark cycles with lights on at 07:00. Cages were cleaned daily and provisioned with fresh drinking water, commercial monkey chow, and a variety of fruit and vegetable supplements. Various toys, swinging perches, and simulated foraging activities were provided for environmental enrichment.
In naturalistic and semi-naturalistic conditions, squirrel monkeys live in sexually segregated groups comprised of multiple females and many fewer males . To simulate this social organization in our control condition, we randomly assigned 10 males to live with a familiar same-sex conspecific in undisturbed pairs and 24 females to live in undisturbed groups each comprised of four familiar same-sex conspecifics. Another age-matched sample of 12 males and 24 females were randomized to the stress treatment condition in which monkeys housed as described above were intermittently removed from same-sex conspecifics at 12–15-week intervals for 6–7 repeated social separations. During each 3-week separation session, monkeys were housed individually in cages that allowed visual, auditory, olfactory, and limited tactile contact between adjacent animals.
Prefrontal brain tissue collection
Brains from randomly selected monkeys in the stress (n = 6 males and 6 females) and no-stress (n = 6 males and 12 females) conditions were collected 4 weeks after completion of the treatment conditions when all monkeys were socially housed with same-sex conspecifics. Collection time was selected to assess long-term effects of stress. After euthanasia, brains were collected using established procedures  between 08:00 and 10:00 to control for circadian effects. All brains were collected during nonbreeding seasons when circulating sex steroid hormones remain stable at low levels in these seasonally breeding primates .
BeadChip array hybridization
Tissue samples were homogenized with a motorized pellet pestle, and total RNA was extracted using AllPrep DNA/RNA Micro kits (Qiagen). Samples were amplified (RiboAmp Plus 1.5-round RNA Amplification, Applied Biosystems) for the production of biotinylated cRNA (Bioarray High Yield RNA Labeling, ENZO) and subsequently hybridized to HumanHT-12v4.0 Expression BeadChip arrays (Illumina) that were scanned on a BeadStation system following manufacturer’s instructions. A total of 90 BeadChip arrays were used to generate gene expression profiles for dorsolateral, ventrolateral, and ventromedial prefrontal cortex from 18 females and 12 males. No pooling of samples across monkeys or regions was necessary. Prefrontal cortical regions, sex, and stress treatment conditions were randomly counterbalanced across separate BeadChip arrays.
To verify array results, we focused on HSD11B1 and CACNG2 expression in each of the three prefrontal cortical regions assessed by quantitative real-time polymerase chain reaction (qPCR) with ACTB as the reference gene. Primers designed with Primer-BLAST (NCBI) for Sybr Green assays with parameters set as described elsewhere  are as follows: HSD11B1 forward 5′ ATGTGGTGGTGACAGCGAG 3′; HSD11B1 reverse 5′ TATAGTGCGCTGAGGCTGCT 3′; CACNG2 forward 5′ GGCCCTGTCCTTCATCATC 3′; CACNG2 reverse 5′ GCAGAGGCCTGGAGGTAGT 3′; ACTB forward 5′ CAAGGCCAATCGTGAGAAGA 3′; and ACTB reverse 5′ AGAGGCGTACAAGGAAAGCA 3′. The same samples analyzed by array were used for qPCR reactions conducted following manufacturer’s instructions (Biorad, SsoAdvanced Universal SYBR Green) under these conditions: 250 nM primers, 50 mM Na+, 3 mM Mg++, and 1.2 mM dNTP, using the MxPro3000 (Stratagene). Results from qPCR were analyzed as described elsewhere .
Array data processing
Each BeadChip array has more than 47,000 probes designed for the human genome. To identify probes suited for squirrel monkeys, we compared probe sequences against the squirrel monkey genome (Broad Institute, GCA_000235385.1) using BLAT . Selected probes were required to match a single continuous segment of the squirrel monkey genome at 95% homology or greater. For the 11,209 selected probes, expression values were imported into R , background corrected using maximum likelihood estimation , log2 transformed, and quantile normalized. Normalized data were analyzed for outliers using box plots, sample-sample correlation matrices, and principle components analysis. Two males from the stress condition were excluded as outliers with a median sample-sample correlation of 0.946 compared to 0.980 for all remaining samples. Probes not associated with known genes were discarded. Using these procedures, we analyzed 8853 probes interrogating 7299 unique genes in 18 females and 10 males for 3 prefrontal cortical regions.
A multi-level regression model also known as a hierarchical linear model or a mixed-effects model  was used to assess stress, sex, and prefrontal cortical region main effects in expression data for each probe. These analyses were performed in R (R v.3.3.0: https://cran.r-project.org, RStudio v.0.99.896: http://www.rstudio.com/) using the lme() function in the nlme package (https://cran.r-project.org/package=nlme) with a default correlation structure (no within-group correlation) and method (REML: maximizing the restricted log-likelihood). Our model included the fixed effects of stress (reference level = control condition), sex (reference level = female), and prefrontal region (reference level = dorsolateral prefrontal cortex) with repeated measures accounted for by including individual subjects as a random effects variable (random intercept). Age was unevenly distributed across stress and sex groupings, and in humans, age alters frontal cortical gene expression profiles . We therefore included age class as a fixed effect covariate in our model for squirrel monkeys (younger adult = − 1, middle-age adult = 0, older adult = 1). Age class categories were used to accommodate age estimates based on dentition. All test statistics were evaluated with two-tailed probabilities, and Benjamini-Hochberg false discovery rate corrected q values were computed for each nominal P value to address the multiple comparisons issue .
To further explore regional differences in gene expression profiles, a single summary statistic for regional effects on expression of each probe was extracted using a likelihood ratio approach. Specifically, we used anova() function in R to compare the full multilevel model described above with a model that lacked the fixed effect term for region . For probes with significant regional effects (P < 0.05) that survived multiple testing correction (q < 0.05), hierarchical cluster analysis with average linkage was used to assess relations between regions in SYSTAT13 (http://systatsoftware.com/).
We also examined how regional differences in gene expression profiles of prefrontal cortex in squirrel monkeys compare to humans. To identify regions of human prefrontal cortex with gene expression profiles similar to regions of squirrel monkey prefrontal cortex, we compared monkey results to Allen Brain Atlas human brain microarray data (http://human.brain-map.org, downloaded 12/2015). Allen Brain Atlas microarray data span 160 cortical and subcortical brain regions in high-quality post-mortem tissue from healthy adult men and women . Preprocessed microarray data from the Allen Brain Atlas website were downloaded as z-scores with log2 expression data for each probe transformed so that the mean for the entire brain was equal to 0 and the standard deviation equal to 1. For our purposes, we then extracted human expression data from 20 specific regions of frontal or anterior cingulate cortex listed in Additional file 1: Table S1 and re-centered and rescaled these data for each probe. We also centered and scaled our squirrel monkey data for each probe to make similar formats for monkeys and humans, and averaged data for each probe by region. Average regional gene expression signatures for both species were then compared using a region-region correlation matrix and hierarchical clustering.
To identify molecular pathways in monkey gene expression data, probes with significant sex differences (P < 0.05) that survived multiple testing correction (q < 0.05) were uploaded to DAVID6.7  for annotation with Kyoto Encyclopedia of Genes and Genomes (KEGG) terms. The complete set of 8853 probes that targeted 7299 unique genes identified as suitable for squirrel monkeys served as background for all pathway analyses. Chromosome locations are unknown for many squirrel monkey genes, and therefore we used a comparative approach with Ensembl BioMart to identify genes on sex chromosomes in humans and marmoset monkeys (Callithrix jacchus, C_jacchus3.2.1). Like squirrel monkeys, marmosets are also New World primates. Effect sizes for sex differences in gene expression were calculated using Eta squared η 2  and assessed by analysis of variance in SYSTAT13 with stress as a between-subjects factor and region considered a within-subjects factor.
Prefrontal cortical region main effects were discerned for 5884 probes that targeted 5155 unique genes (2.6 × 10−33 ≤ P ≤ 0.05). Nominal P values for 5557 of these probes that targeted 4905 genes survived multiple testing correction (2.3 × 10−29 ≤ q ≤ 0.05; Additional file 2: Table S2). More than 240 probes with region main effects (P < 0.05) that survived multiple testing correction (q < 0.05) targeted X chromosome genes in humans and marmoset monkeys (Additional file 2: Table S2). None of the probes in our gene expression profiling analyses of squirrel monkeys targeted Y chromosome genes in humans or marmosets.
Cross-species comparisons indicated that squirrel monkey dorsolateral prefrontal cortex resembled middle frontal gyrus (R = 0.197, P = 4.98 × 10−87) and precentral gyrus (R = 0.183, P = 1.06 × 10−74) in humans portrayed in Allen Brain Atlas microarray data. Squirrel monkey ventromedial prefrontal cortex resembled parolfactory (or subcallosal) gyri (R = 0.111, P = 1.71 × 10−28) and gyrus rectus (R = 0.107, P = 2.00 × 10−26) in humans. Squirrel monkey ventrolateral prefrontal cortex resembled posterior orbital gyrus (R = 0.047, P = 2.49 × 10−6) and frontal pole (R = 0.063, P = 2.49 × 10−6) in humans.
Sex difference main effects in squirrel monkeys were discerned for 1270 probes that targeted 1219 unique genes (8.5 × 10−8 < P < 0.05). Nominal P values for 82 of these probes that targeted 80 unique genes survived multiple testing correction (3.8 × 10−4 ≤ q ≤ 0.05; Additional file 3: Table S3). Of these 82 probes, 19 had nominally significant age class effects (P < 0.05) that did not survive multiple testing correction (Additional file 3: Table S3). For these 19 probes, unbalanced age distributions could potentially contribute to differences between the sexes but this possibility does not likely explain sex differences for the remaining 63 probes.
Most genes expressed differently between the sexes across squirrel monkey prefrontal cortical regions did not reside on sex chromosomes. None of the probes targeted Y chromosome genes and only 2 of 82 probes with sex difference main effects (P < 0.05) that survived multiple testing correction (q < 0.05) targeted X chromosome genes in humans and marmoset monkeys. Both of these genes (SLC25A5 and TSC22D3) were expressed at higher levels in squirrel monkey females compared to males across all three prefrontal cortical regions.
Many probes expressed at higher levels in squirrel monkey females compared to males targeted autosomal genes. For example, HSD11B1 is located on chromosome 1 in humans and this gene was expressed at higher levels in squirrel monkey females compared to males across prefrontal cortical regions discerned by two probes on the array (Additional file 3: Table S3). Array results for both probes correlated with HSD11B1 expression determined by qPCR (Additional file 4: Table S4) and qPCR determinations were greater in females compared to males (F(1,24) = 17.2, P < 0.001 for the sex main effect). A region main effect was also discerned by qPCR (F(2,48) = 21.3, P < 0.001) and by both array probes (Additional file 2: Table S2) for HSD11B1 expression. Regional differences do not likely explain correlations because qPCR correlated with array measures for HSD11B1 regardless of whether prefrontal cortical regions were analyzed separately or combined (Additional file 4: Table S4). Comparable correlations were also identified for another autosomal gene (CACNG2 in Additional file 4: Table S4) that differed significantly between prefrontal cortical regions determined by array (Additional file 2: Table S2) and by qPCR (F(2,48) = 6.8, P = 0.002 for the region main effect).
Of the 82 probes with sex difference main effects (P < 0.05) that survived multiple testing correction (q < 0.05), 64 (78%) were increased across prefrontal cortical regions in squirrel monkey females compared to males. Conversely, 18 of 82 probes (22%) were increased across prefrontal cortical regions in males compared to females. Molecular pathways of genes targeted by probes with sex difference main effects were assessed using DAVID6.7 .
Genes in two molecular pathways. Specific genes overrepresented in neurotrophin and MAPK signaling pathways with greater expression in squirrel monkey females compared to males (mean ± SEM) across prefrontal cortical regions
7.63 ± 0.08
7.16 ± .012
8.77 ± 0.08
8.24 ± .011
7.94 ± 0.07
7.38 ± 0.10
8.31 ± 0.08
7.87 ± 0.12
Neurotrophin and MAPK signaling
6.54 ± 0.06
6.19 ± 0.09
11.88 ± 0.05
11.58 ± 0.06
11.85 ± 0.05
11.55 ± 0.06
5.93 ± 0.06
5.72 ± 0.07
12.09 ± 0.04
11.82 ± 0.05
To test for greater variation in gene expression associated with estrous cycle-related hormonal fluctuations, we followed a statistical approach described elsewhere by other investigators . Coefficients of variation (CV) were calculated separately for no-stress control females and control males for each of the 82 probes with sex difference main effects (P < 0.05) that survived multiple testing correction (q < 0.05). Twelve females from the control condition were examined to increase chances of assessing different stages of the estrous cycle. Sex-specific CVs were not significantly greater in females compared to males for any prefrontal cortical region.
Top 25 genes with the largest disparity in effect sizes for sex differences in stress versus no-stress conditions. Effect sizes (η 2) for sex differences in stress and no-stress conditions are presented for ventrolateral (vlPFC) and ventromedial (vmPFC) prefrontal cortical regions. F > M signifies genes with greater expression in females compared to males, and M > F signifies the reverse
F < M
F < M
F < M
F < M
F < M
F < M
F < M
F < M
F < M
F < M
F < M
F < M
F < M
F < M
F < M
F < M
F < M
F < M
F < M
M > F
F < M
F < M
F < M
F < M
We examined gene expression profiles for dorsolateral, ventrolateral, and ventromedial prefrontal cortical regions from female and male squirrel monkey adults randomly selected from a larger sample randomized to stress or no-stress conditions. Human-specific probes were aligned to the monkey genome and those with ≥ 95% sequence homology were selected for analysis. Results indicate extensive variation between prefrontal cortical regions in the expression of numerous autosomal and sex chromosome genes. Robust sex differences were also discerned across prefrontal cortical regions in the expression of mostly autosomal genes. Genes with increased expression in females compared to males were overrepresented in MAPK and neurotrophin signaling pathways. Many fewer genes with increased expression in males compared to females were discerned, and no molecular pathways were identified. Effect sizes for sex differences were greater in stress compared to no-stress conditions for ventromedial and ventrolateral prefrontal cortical regions but not dorsolateral prefrontal cortex.
Differential gene expression profiles have been reported for diverse brain regions in humans [9, 10] and various monkey species [11, 12] but comparisons between primate prefrontal cortical regions are uncommon. We found that gene expression profiles for prefrontal cortical regions in squirrel monkeys resembled profiles of anatomically comparable prefrontal regions in humans. Squirrel monkey dorsolateral and ventromedial gene expression profiles more closely resembled one another compared to ventrolateral prefrontal cortex. Dorsolateral prefrontal cortex is involved in working memory and cognitive control whereas ventral regions play a role in emotion regulation [50–52]. How prefrontal cortical gene expression profiles mediate behavioral functions remain to be explored but primate prefrontal profiles are spatially distinct.
Robust sex differences were identified across prefrontal cortical regions despite initial suggestions that sex differences in gene expression profiles are limited in human brain . Recent results instead suggest that many genes are differentially expressed between the sexes for various brain regions in mice [18, 19], monkeys , and humans [9, 22, 23]. Our results agree and confirm that 82 probes that targeted 80 unique genes are differentially expressed in male and female squirrel monkeys across prefrontal cortical regions.
Similar sex differences in different species raise the possibility of a conserved sexual signature in gene expression profiles of brain. Reinius et al.  found 85 genes differentially expressed between the sexes for occipital cortex in humans and cynomologus macaque monkeys. Only two of these genes were differentially expressed in occipital cortex of female and male marmoset monkeys . These two genes were not examined in our study of squirrel monkeys but another gene in squirrel monkeys overlapped with the Reinius et al.  list for humans and macaques. In all three species, SLC6A1 is upregulated in cerebral cortex of females compared to males. Molecular genetic phylogenies suggest that macaques more closely resemble humans than do either squirrel monkeys or marmosets  but macaques and squirrel monkeys are polygamous and sexually dimorphic in size whereas marmoset males and females are monogamous and similar in size .
SLC6A1 encodes gamma-aminobutyric acid (GABA) transporter 1 (GAT-1) which terminates GABA neurotransmission via synaptic reuptake specifically in cerebral cortex . Selective GAT-1 inhibitors such as tiagabine have anxiolytic and antidepressant effects in mice , and novel GAT-1 inhibitors are now being tested as potential treatments for anxiety and depression . Studies of GAT-1 inhibitors have thus far focused primarily on males [56, 57], but our results suggest that sex differences warrant attention.
Another autosomal gene expressed at higher levels in squirrel monkey females compared to males across prefrontal cortical regions is HSD11B1. This gene encodes an enzyme that regulates glucocorticoid exposure by catalyzing conversion of inactive to active glucocorticoids . A recent study of 134 healthy humans between the ages of 20 and 81 years reported increased skeletal muscle HSD11B1 expression that correlated with age in women but not men . We did not find significant correlations between age and HSD11B1 expression in any prefrontal cortical region for our smaller sample of female or male squirrel monkeys.
In humans with major depression, we recently identified FGF9 as a gene of interest . Administration of FGF9 to male rats increases anxiety and depression-like behavior , and FGF9 stimulates expression of SPRY4 in vitro . In squirrel monkeys, FGF9 did not differ between the sexes but SPRY4 was increased across prefrontal cortical regions in females compared to males. SPRY4 encodes a negative regulator that suppresses neurotrophic functions [62, 63]. Inhibition of negative regulators that suppress neurotrophic functions may provide a novel approach for the development of new antidepressant medications .
Differential expression of single genes is noteworthy, but collective differences in multiple genes along functional pathways provide further insights . Both MAPK and neurotrophin signaling pathways were significantly enriched with genes expressed at higher levels in squirrel monkey females compared to males across prefrontal cortical regions. More than half of the genes enriched in these pathways were associated with both MAPK and neurotrophin signaling. Neurotrophins activate MAPK signaling via extracellular signal-regulated kinases, and MAPK signaling is also activated by cJun N-terminal kinases elicited by cellular stress and pro-inflammatory conditions . These findings agree with evidence that women may be more vulnerable to the depressogenic effects of inflammation compared to men .
In human dorsal prefrontal cortex, the majority of genes differentially expressed between the sexes are autosomal and not sex chromosome genes . Most sex differences in squirrel monkey prefrontal cortical regions are likewise autosomal, but two X chromosome genes were identified. Both SLC25A5 and TSC22D3 may have escaped X chromosome inactivation [67, 68] insofar as higher expression levels were discerned in squirrel monkey females compared to males across prefrontal cortical regions. TSC22D3 is induced by stress levels of glucocorticoids  and encodes a leucine zipper protein linked to MAPK signaling  and neuroplasticity in depression . SLC25A5 encodes a mitochondrial protein  associated with anxiety in mice .
None of the stress main effects nor interactions in our study survived multiple testing correction possibly because of exposure to a mild intermittent social stressor. Nevertheless, effect sizes for sex differences were greater in stress compared to no-stress conditions for ventromedial and ventrolateral prefrontal cortical regions but not dorsolateral prefrontal cortex. Functional differences between regions suggest that stress modulates effect sizes for sex differences in prefrontal cortical regions involved in emotion regulation and not working memory nor cognitive control [50–52]. Two genes provide examples. Effect sizes for sex differences in SLC6A1 and SPRY4 expression were two to sevenfold greater in stress compared to no-stress conditions for ventromedial and ventrolateral prefrontal cortical regions. Increased SLC6A1 and SPRY4 expression in females compared to males during stressful conditions may respectively accelerate termination of GABA neurotransmission and inhibit neurotrophic functions. As discussed above, GABA neurotransmission and neurotrophic functions have been both implicated in mood and anxiety disorders [56, 60].
Greater effect sizes for sex differences in stress compared to no-stress conditions have additional implications for calls to include both sexes in preclinical animal research . Studies not explicitly focused on stress often include manipulations or conditions that are inherently stressful. Our findings suggest that certain sex differences are amplified in stressful situations. If our findings from monkeys hold true for humans and other animals, then studies that include unintentional stressors with males alone may overlook significant sources of biological variability.
Our results should be interpreted in the context of potential limitations. Objectively similar stressors may be viewed differently by the sexes and small samples restricted our capacity to detect stress, sex, and region interactions. Although we analyzed 8853 probes that targeted 7299 unique genes, many probes designed for humans do not match the squirrel monkey genome and were not considered. Our array results require further validation, and bulk tissue samples do not necessarily reflect specific cell types or cellular stress states. Tissue samples were collected during non-breeding seasons when circulating sex steroid hormones remain stable at low levels, but organizational effects of sex steroids on brain development were not considered.
Prefrontal cortical regions in squirrel monkey adults differed in the expression of numerous autosomal and sex chromosome genes. Robust sex differences were also identified across prefrontal cortical regions in the expression of mostly autosomal genes. Stress amplified sex difference effect sizes for gene expression profiles in prefrontal cortical regions involved in stress coping and emotion regulation. Results suggest new molecular targets for sex-specific treatments of stress disorders in human mental health.
We are grateful to I. Birt, A. Pankonin, and D. Romero Vargas for help with array data processing.
This research was supported by National Institutes of Health MH47573 and the Pritzker Neuropsychiatric Disorders Research Consortium Fund LLC.
Availability of data and materials
Data are available for public use at https://pritzkerneuropsych.org/www/data/.
DA collected the data. MH, KEM, and TLB analyzed the data. RMM, SJW, HA, and AFS designed the study and interpreted data. AGL and DML worked on all aspects of this research. All authors read and approved the final manuscript.
All procedures with animals were conducted in accordance with state and federal laws, standards of the Department of Health and Human Services, and were approved by the Institutional Animal Care and Use Committee at Stanford University.
Drs. Morrison and Bale report no biomedical financial interests or potential conflicts of interest related to this work. All other authors are members of the Pritzker Neuropsychiatric Disorders Research Consortium, which is supported by Pritzker Neuropsychiatric Disorders Research Fund, LLC. A shared intellectual property agreement exists between the academic and philanthropic entities of the Consortium. The Pritzker Neuropsychiatric Disorders Research Fund had no role in study design, data collection or analysis, decision to publish, or preparation of the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Altemus M, Sarvaiya N, Neill EC. Sex differences in anxiety and depression clinical perspectives. Front Neuroendocrinol. 2014;35:320–30.View ArticlePubMedPubMed CentralGoogle Scholar
- McHenry J, Carrier N, Hull E, Kabbaj M. Sex differences in anxiety and depression: role of testosterone. Front Neuroendocrinol. 2014;35:42–57.View ArticlePubMedGoogle Scholar
- Seedat S, Scott KM, Angermeyer MC, Berglund P, Bromet EJ, Brugha TS, Demyttenaere K, de Girolamo G, Haro JM, Jin R, et al. Cross-national associations between gender and mental disorders in the World Health Organization World Mental Health Surveys. Arch Gen Psychiatry. 2009;66:785–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Bale TL, Epperson CN. Sex differences and stress across the lifespan. Nat Neurosci. 2015;18:1413–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Oikonomidis L, Santangelo AM, Shiba Y, Clarke FH, Robbins TW, Roberts AC. A dimensional approach to modeling symptoms of neuropsychiatric disorders in the marmoset monkey. Dev Neurobiol. 2017;77:328–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Lueken U, Hahn T. Functional neuroimaging of psychotherapeutic processes in anxiety and depression: from mechanisms to predictions. Curr Opin Psychiatry. 2016;29:25–31.View ArticlePubMedGoogle Scholar
- McEwen BS, Morrison JH. The brain on stress: vulnerability and plasticity of the prefrontal cortex over the life course. Neuron. 2013;79:16–29.View ArticlePubMedPubMed CentralGoogle Scholar
- Beery AK, Zucker I. Sex bias in neuroscience and biomedical research. Neurosci Biobehav Rev. 2011;35:565–72.View ArticlePubMedGoogle Scholar
- Trabzuni D, Ramasamy A, Imran S, Walker R, Smith C, Weale ME, Hardy J, Ryten M. North American brain expression C. Widespread sex differences in gene expression and splicing in the adult human brain. Nat Commun. 2013;4:2771.View ArticlePubMedPubMed CentralGoogle Scholar
- Hawrylycz M, Miller JA, Menon V, Feng D, Dolbeare T, Guillozet-Bongaarts AL, Jegga AG, Aronow BJ, Lee CK, Bernard A, et al. Canonical genetic signatures of the adult human brain. Nat Neurosci. 2015;18:1832–44.View ArticlePubMedPubMed CentralGoogle Scholar
- Bernard A, Lubbers LS, Tanis KQ, Luo R, Podtelezhnikov AA, Finney EM, McWhorter MM, Serikawa K, Lemon T, Morgan R, et al. Transcriptional architecture of the primate neocortex. Neuron. 2012;73:1083–99.View ArticlePubMedPubMed CentralGoogle Scholar
- Fukuoka T, Sumida K, Yamada T, Higuchi C, Nakagaki K, Nakamura K, Kohsaka S, Saito K, Oeda K. Gene expression profiles in the common marmoset brain determined using a newly developed common marmoset-specific DNA microarray. Neurosci Res. 2010;66:62–85.View ArticlePubMedGoogle Scholar
- Ochsner KN, Silvers JA, Buhle JT. Functional imaging studies of emotion regulation: a synthetic review and evolving model of the cognitive control of emotion. Ann N Y Acad Sci. 2012;1251:E1–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Ray RD, Zald DH. Anatomical insights into the interaction of emotion and cognition in the prefrontal cortex. Neurosci Biobehav Rev. 2012;36:479–501.View ArticlePubMedGoogle Scholar
- Collins KA, Mendelsohn A, Cain CK, Schiller D. Taking action in the face of threat: neural synchronization predicts adaptive coping. J Neurosci. 2014;34:14733–8.View ArticlePubMedPubMed CentralGoogle Scholar
- McKlveen JM, Myers B, Herman JP. The medial prefrontal cortex: coordinator of autonomic, neuroendocrine and behavioural responses to stress. J Neuroendocrinol. 2015;27:446–56.View ArticlePubMedPubMed CentralGoogle Scholar
- Rinn JL, Snyder M. Sexual dimorphism in mammalian gene expression. Trends Genet. 2005;21:298–305.View ArticlePubMedGoogle Scholar
- Yang X, Schadt EE, Wang S, Wang H, Arnold AP, Ingram-Drake L, Drake TA, Lusis AJ. Tissue-specific expression and regulation of sexually dimorphic genes in mice. Genome Res. 2006;16:995–1004.View ArticlePubMedPubMed CentralGoogle Scholar
- Vied C, Ray S, Badger CD, Bundy JL, Arbeitman MN, Nowakowski RS. Transcriptomic analysis of the hippocampus from six inbred strains of mice suggests basis for sex-specific susceptibility and severity of neurological disorders. J Comp Neurol. 2016;524:2696–710.View ArticlePubMedGoogle Scholar
- Morgan CP, Bale TL. Sex differences in microRNA regulation of gene expression: no smoke, just miRs. Biol Sex Differ. 2012;3:22.View ArticlePubMedPubMed CentralGoogle Scholar
- Reinius B, Saetre P, Leonard JA, Blekhman R, Merino-Martinez R, Gilad Y, Jazin E. An evolutionarily conserved sexual signature in the primate brain. PLoS Genet. 2008;4:e1000100.View ArticlePubMedPubMed CentralGoogle Scholar
- Xu H, Wang F, Liu Y, Yu Y, Gelernter J, Zhang H. Sex-biased methylome and transcriptome in human prefrontal cortex. Hum Mol Genet. 2014;23:1260–70.View ArticlePubMedGoogle Scholar
- Mayne BT, Bianco-Miotto T, Buckberry S, Breen J, Clifton V, Shoubridge C, Roberts CT. Large scale gene expression meta-analysis reveals tissue-specific, sex-biased gene expression in humans. Front Genet. 2016;7:183.View ArticlePubMedPubMed CentralGoogle Scholar
- Crews D. Diversity of hormone-behavior relations in reproductive behavior. In: Becker JB, Breedlove SM, Crews D, McCarthy MM, editors. Behavioral endocrinology. 2nd ed. Cambridge: MIT Press; 2002. p. 223–87.Google Scholar
- Palanza P, Parmigiani S. How does sex matter? Behavior, stress and animal models of neurobehavioral disorders. Neurosci Biobehav Rev. 2017;76(Pt A):134–43.View ArticlePubMedGoogle Scholar
- Mendoza SP, Lyons DM, Saltzman W. Sociophysiology of squirrel monkeys. Am J Primatol. 1991;23:37–54.View ArticleGoogle Scholar
- Mendoza SP, Hennessy MB, Lyons DM. Distinct immediate and prolonged effects of separation on plasma cortisol in female squirrel monkeys. Psychobiology. 1992;20:300–6.Google Scholar
- Karssen AM, Her S, Li JZ, Patel PD, Meng F, Bunney WE Jr, Jones EG, Watson SJ, Akil H, Myers RM, et al. Stress-induced changes in primate prefrontal profiles of gene expression. Mol Psychiatry. 2007;12:1089–102.View ArticlePubMedGoogle Scholar
- Lee AG, Buckmaster CL, Yi E, Schatzberg AF, Lyons DM. Coping and glucocorticoid receptor regulation by stress inoculation. Psychoneuroendocrinology. 2014;49:272–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Lyons DM, Buckmaster PS, Lee AG, Wu C, Mitra R, Duffey LM, Buckmaster CL, Her S, Patel PD, Schatzberg AF. Stress coping stimulates hippocampal neurogenesis in adult monkeys. Proc Natl Acad Sci U S A. 2010;107:14823–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Akil H. Scientific strategy in neuroscience: discovery science versus hypothesis-driven research. Neurosci Quarterly. 2003; Spring:4–5, https://www.sfn.org/news-and-calendar/neuroscience-quarterly/past-issues. Accessed 1 Oct 2017.
- Quackenbush J. Computational approaches to analysis of DNA microarray data. Yearb Med Inform. 2006;45:91–103.Google Scholar
- Lyons DM, Mendoza SP, Mason WA. Sexual segregation in squirrel monkeys (Saimiri sciureus): a transactional analysis of adult social dynamics. J Comp Psychol. 1992;106:323–30.View ArticlePubMedGoogle Scholar
- Patel PD, Katz M, Karssen AM, Lyons DM. Stress-induced changes in corticosteroid receptor expression in primate hippocampus and prefrontal cortex. Psychoneuroendocrinology. 2008;33:360–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Schiml PA, Mendoza SP, Saltzman W, Lyons DM, Mason WA. Annual physiological changes in individually housed squirrel monkeys (Saimiri sciureus). Am J Primatol. 1999;47:93–103.View ArticlePubMedGoogle Scholar
- Gergen JA, MacLean PD. A stereotaxic atlas of the squirrel monkey's brain (Saimiri sciureus) vol. 933. Bethesda: Public Health Service; 1962.View ArticleGoogle Scholar
- Thornton B, Basu C. Real-time PCR (qPCR) primer design using free online software. Biochem Mol Biol Educ. 2011;39:145–54.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25:402–8.View ArticlePubMedGoogle Scholar
- Kent WJ. BLAT—the BLAST-like alignment tool. Genome Res. 2002;12:656–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Core Team R. R: a language and environment for statistical computing. Vienna. Austria: R Foundation for Statistical Computing; 2013.Google Scholar
- Xie Y, Wang X, Story M. Statistical methods of background correction for Illumina BeadArray data. Bioinformatics. 2009;25:751–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Maxwell SE, Delaney HD. Designing experiments and analyzing data: a model comparison perspective. 2nd ed. Mahwah: Lawerence Erlbaum Associates; 2004. p. 763–820.Google Scholar
- Lu T, Pan Y, Kao SY, Li C, Kohane I, Chan J, Yankner BA. Gene regulation and DNA damage in the ageing human brain. Nature. 2004;429:883–91.View ArticlePubMedGoogle Scholar
- Storey JD, Tibshirani R. Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003;100:9440–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Bates D, Machler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.View ArticleGoogle Scholar
- Hawrylycz MJ, Lein ES, Guillozet-Bongaarts AL, Shen EH, Ng L, Miller JA, van de Lagemaat LN, Smith KA, Ebbert A, Riley ZL, et al. An anatomically comprehensive atlas of the adult human brain transcriptome. Nature. 2012;489:391–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.View ArticlePubMedGoogle Scholar
- Lakens D. Calculating and reporting effect sizes to facilitate cumulative science: a practical primer for t-tests and ANOVAs. Front Psychol. 2013;4:863.View ArticlePubMedPubMed CentralGoogle Scholar
- Becker JB, Prendergast BJ, Liang JW. Female rats are not more variable than male rats: a meta-analysis of neuroscience studies. Biol Sex Differ. 2016;7:34.View ArticlePubMedPubMed CentralGoogle Scholar
- Mitchell DG. The nexus between decision making and emotion regulation: a review of convergent neurocognitive substrates. Behav Brain Res. 2011;217:215–31.View ArticlePubMedGoogle Scholar
- Smolker HR, Depue BE, Reineberg AE, Orr JM, Banich MT. Individual differences in regional prefrontal gray matter morphometry and fractional anisotropy are associated with different constructs of executive function. Brain Struct Funct. 2015;220:1291–306.View ArticlePubMedGoogle Scholar
- Mueller SC. The influence of emotion on cognitive control: relevance for development and adolescent psychopathology. Front Psychol. 2011;2:327.View ArticlePubMedPubMed CentralGoogle Scholar
- Shoshani J, Groves CP, Simons EL, Gunnell GF. Primate phylogeny: morphological vs. molecular results. Mol Phylogenet Evol. 1996;5:102–54.View ArticlePubMedGoogle Scholar
- Smuts BR, Cheney DL, Seyfarth RM, Wrangham RW, Struhsaker TT. Primate societies. Chicago: University of Chicago Press; 1986.Google Scholar
- Conti F, Minelli A, Melone M. GABA transporters in the mammalian cerebral cortex: localization, development and pathological implications. Brain Res Brain Res Rev. 2004;45:196–212.View ArticlePubMedGoogle Scholar
- Thoeringer CK, Erhardt A, Sillaber I, Mueller MB, Ohl F, Holsboer F, Keck ME. Long-term anxiolytic and antidepressant-like behavioural effects of tiagabine, a selective GABA transporter-1 (GAT-1) inhibitor, coincide with a decrease in HPA system activity in C57BL/6 mice. J Psychopharmacol. 2010;24:733–43.View ArticlePubMedGoogle Scholar
- Salat K, Podkowa A, Malikowska N, Kern F, Pabel J, Wojcieszak E, Kulig K, Wanner KT, Strach B, Wyska E. Novel, highly potent and in vivo active inhibitor of GABA transporter subtype 1 with anticonvulsant, anxiolytic, antidepressant and antinociceptive properties. Neuropharmacology. 2017;113(Pt A):331–42.View ArticlePubMedGoogle Scholar
- Wyrwoll CS, Holmes MC, Seckl JR. 11beta-hydroxysteroid dehydrogenases and the brain: from zero to hero, a decade of progress. Front Neuroendocrinol. 2011;32:265–86.View ArticlePubMedPubMed CentralGoogle Scholar
- Hassan-Smith ZK, Morgan SA, Sherlock M, Hughes B, Taylor AE, Lavery GG, Tomlinson JW, Stewart PM. Gender-specific differences in skeletal muscle 11beta-HSD1 expression across healthy aging. J Clin Endocrinol Metab. 2015;100:2673–81.View ArticlePubMedGoogle Scholar
- Aurbach EL, Inui EG, Turner CA, Hagenauer MH, Prater KE, Li JZ, Absher D, Shah N, Blandino P Jr, Bunney WE, et al. Fibroblast growth factor 9 is a novel modulator of negative affect. Proc Natl Acad Sci U S A. 2015;112:11953–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Sylvestersen KB, Herrera PL, Serup P, Rescan C. FGF9 signalling stimulates Spred and Sprouty expression in embryonic mouse pancreas mesenchyme. Gene Expr Patterns. 2011;11:105–11.View ArticlePubMedGoogle Scholar
- Katoh Y, Katoh M. FGF signaling inhibitor, SPRY4, is evolutionarily conserved target of WNT signaling pathway in progenitor cells. Int J Mol Med. 2006;17:529–32.PubMedGoogle Scholar
- Kim HJ, Bar-Sagi D. Modulation of signalling by Sprouty: a developing story. Nat Rev Mol Cell Biol. 2004;5:441–50.View ArticlePubMedGoogle Scholar
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Cargnello M, Roux PP. Activation and function of the MAPKs and their substrates, the MAPK-activated protein kinases. Microbiol Mol Biol Rev. 2011;75:50–83.View ArticlePubMedPubMed CentralGoogle Scholar
- Bekhbat M, Neigh GN. Sex differences in the neuro-immune consequences of stress: focus on depression and anxiety. Brain Behav Immun. 2017; doi.org/10.1016/j.bbi.2017.02.006.
- Carrel L, Willard HF. X-inactivation profile reveals extensive variability in X-linked gene expression in females. Nature. 2005;434:400–4.View ArticlePubMedGoogle Scholar
- Peeters SB, Cotton AM, Brown CJ. Variable escape from X-chromosome inactivation: identifying factors that tip the scales towards expression. BioEssays. 2014;36:746–56.View ArticlePubMedPubMed CentralGoogle Scholar
- Ayroldi E, Riccardi C. Glucocorticoid-induced leucine zipper (GILZ): a new important mediator of glucocorticoid action. FASEB J. 2009;23:3649–58.View ArticlePubMedGoogle Scholar
- Clark AR, Lasa M. Crosstalk between glucocorticoids and mitogen-activated protein kinase signalling pathways. Curr Opin Pharmacol. 2003;3:404–11.View ArticlePubMedGoogle Scholar
- Frodl T, Carballedo A, Hughes MM, Saleh K, Fagan A, Skokauskas N, McLoughlin DM, Meaney J, O'Keane V, Connor TJ. Reduced expression of glucocorticoid-inducible genes GILZ and SGK-1: high IL-6 levels are associated with reduced hippocampal volumes in major depressive disorder. Transl Psychiatry. 2012;2:e88.View ArticlePubMedPubMed CentralGoogle Scholar
- Clemencon B, Babot M, Trezeguet V. The mitochondrial ADP/ATP carrier (SLC25 family): pathological implications of its dysfunction. Mol Asp Med. 2013;34:485–93.View ArticleGoogle Scholar
- Filiou MD, Zhang Y, Teplytska L, Reckow S, Gormanns P, Maccarrone G, Frank E, Kessler MS, Hambsch B, Nussbaumer M, et al. Proteomics and metabolomics analysis of a trait anxiety mouse model reveals divergent mitochondrial pathways. Biol Psychiatry. 2011;70:1074–82.View ArticlePubMedGoogle Scholar
- McCullough LD, de Vries GJ, Miller VM, Becker JB, Sandberg K, McCarthy MM. NIH initiative to balance sex of animals in preclinical studies: generative questions to guide policy, implementation, and metrics. Biol Sex Differ. 2014;5:15.View ArticlePubMedPubMed CentralGoogle Scholar