Allele-specific transcription factor binding across human brain regions offers mechanistic insight into eQTLs.

Transcription factors (TFs) regulate gene expression by facilitating or disrupting the formation of transcription initiation machinery at particular genomic loci. Because TF occupancy is driven in part by recognition of DNA sequence, genetic variation can influence TF-DNA associations and gene regulation. To identify variants that impact TF binding in human brain tissues, we assessed allele-specific binding (ASB) at heterozygous variants for 94 TFs in nine brain regions from two donors. Leveraging graph genomes constructed from phased genomic sequence data, we compared ChIP-seq signals between alleles at heterozygous variants within each brain region and identified thousands of variants exhibiting ASB for at least one TF. ASB reproducibility was measured by comparisons between independent experiments both within and between donors. We found that rare alleles in the general population more frequently led to reduced TF binding, whereas common alleles had an equal likelihood of increasing or decreasing binding. Further, for ASB variants in predicted binding motifs, the favored allele tended to be the one with the stronger expected motif match, but this concordance was not observed within highly occupied sites. We also found that neuron-specific cis-regulatory elements (cCREs), in contrast with oligodendrocyte-specific cCREs, showed depletion of ASB variants. We identified 2670 ASB variants associated with evidence for allele-specific gene expression in the brain from GTEx data and observed increasing eQTL effect direction concordance as ASB significance increases. These results provide a valuable and unique resource for mechanistic analysis of cis-regulatory variation in human brain tissue.

The impact of COVID-19 on a college freshman sample reveals genetic and nongenetic forms of susceptibility and resilience to stress.

Using a longitudinal approach, we sought to define the interplay between genetic and nongenetic factors in shaping vulnerability or resilience to COVID-19 pandemic stress, as indexed by the emergence of symptoms of depression and/or anxiety. University of Michigan freshmen were characterized at baseline using multiple psychological instruments. Subjects were genotyped, and a polygenic risk score for depression (MDD-PRS) was calculated. Daily physical activity and sleep were captured. Subjects were sampled at multiple time points throughout the freshman year on clinical rating scales, including GAD-7 and PHQ-9 for anxiety and depression, respectively. Two cohorts (2019 to 2021) were compared to a pre-COVID-19 cohort to assess the impact of the pandemic. Across cohorts, 26 to 40% of freshmen developed symptoms of anxiety or depression (N = 331). Depression symptoms significantly increased in the pandemic years and became more chronic, especially in females. Physical activity was reduced, and sleep was increased by the pandemic, and this correlated with the emergence of mood symptoms. While low MDD-PRS predicted lower risk for depression during a typical freshman year, this genetic advantage vanished during the pandemic. Indeed, females with lower genetic risk accounted for the majority of the pandemic-induced rise in depression. We developed a model that explained approximately half of the variance in follow-up depression scores based on psychological trait and state characteristics at baseline and contributed to resilience in genetically vulnerable subjects. We discuss the concept of multiple types of resilience, and the interplay between genetic, sex, and psychological factors in shaping the affective response to different types of stressors.

Effect of Selective Lesions of Nucleus Accumbens µ-Opioid Receptor-Expressing Cells on Heroin Self-Administration in Male and Female Rats: A Study with Novel Oprm1-Cre Knock-in Rats.

The brain µ-opioid receptor (MOR) is critical for the analgesic, rewarding, and addictive effects of opioid drugs. However, in rat models of opioid-related behaviors, the circuit mechanisms of MOR-expressing cells are less known because of a lack of genetic tools to selectively manipulate them. We introduce a CRISPR-based Oprm1-Cre knock-in transgenic rat that provides cell type-specific genetic access to MOR-expressing cells. After performing anatomic and behavioral validation experiments, we used the Oprm1-Cre knock-in rats to study the involvement of NAc MOR-expressing cells in heroin self-administration in male and female rats. Using RNAscope, autoradiography, and FISH chain reaction (HCR-FISH), we found no differences in Oprm1 expression in NAc, dorsal striatum, and dorsal hippocampus, or MOR receptor density (except dorsal striatum) or function between Oprm1-Cre knock-in rats and wildtype littermates. HCR-FISH assay showed that iCre is highly coexpressed with Oprm1 (95%-98%). There were no genotype differences in pain responses, morphine analgesia and tolerance, heroin self-administration, and relapse-related behaviors. We used the Cre-dependent vector AAV1-EF1a-Flex-taCasp3-TEVP to lesion NAc MOR-expressing cells. We found that the lesions decreased acquisition of heroin self-administration in male Oprm1-Cre rats and had a stronger inhibitory effect on the effort to self-administer heroin in female Oprm1-Cre rats. The validation of an Oprm1-Cre knock-in rat enables new strategies for understanding the role of MOR-expressing cells in rat models of opioid addiction, pain-related behaviors, and other opioid-mediated functions. Our initial mechanistic study indicates that lesioning NAc MOR-expressing cells had different effects on heroin self-administration in male and female rats.SIGNIFICANCE STATEMENT The brain µ-opioid receptor (MOR) is critical for the analgesic, rewarding, and addictive effects of opioid drugs. However, in rat models of opioid-related behaviors, the circuit mechanisms of MOR-expressing cells are less known because of a lack of genetic tools to selectively manipulate them. We introduce a CRISPR-based Oprm1-Cre knock-in transgenic rat that provides cell type-specific genetic access to brain MOR-expressing cells. After performing anatomical and behavioral validation experiments, we used the Oprm1-Cre knock-in rats to show that lesioning NAc MOR-expressing cells had different effects on heroin self-administration in males and females. The new Oprm1-Cre rats can be used to study the role of brain MOR-expressing cells in animal models of opioid addiction, pain-related behaviors, and other opioid-mediated functions.

Exploratory locomotion, a predictor of addiction vulnerability, is oligogenic in rats selected for this phenotype.

Artificially selected model organisms can reveal hidden features of the genetic architecture of the complex disorders that they model. Addictions are disease phenotypes caused by different intermediate phenotypes and pathways and thereby are potentially highly polygenic. High responder (bHR) and low responder (bLR) rat lines have been selectively bred (b) for exploratory locomotion (EL), a behavioral phenotype correlated with novelty-seeking, impulsive response to reward, and vulnerability to addiction, and is inversely correlated with spontaneous anxiety and depression-like behaviors. The rapid response to selection indicates loci of large effect for EL. Using exome sequencing of HR and LR rats, we identified alleles in gene-coding regions that segregate between the two lines. Quantitative trait locus (QTL) analysis in F2 rats derived from a bHR × bLR intercross confirmed that these regions harbored genes affecting EL. The combined effects of the seven genome-wide significant QTLs accounted for approximately one-third of the total variance in EL, and two-thirds of the variance attributable to genetic factors, consistent with an oligogenic architecture of EL estimated both from the phenotypic distribution of F2 animals and rapid response to selection. Genetic association in humans linked APBA2, the ortholog of the gene at the center of the strongest QTL, with substance use disorders and related behavioral phenotypes. Our finding is also convergent with molecular and animal behavioral studies implicating Apba2 in locomotion. These results provide multilevel evidence for genes/loci influencing EL. They shed light on the genetic architecture of oligogenicity in animals artificially selected for a phenotype modeling a more complex disorder in humans.

Splice-Break: exploiting an RNA-seq splice junction algorithm to discover mitochondrial DNA deletion breakpoints and analyses of psychiatric disorders.

Deletions in the 16.6 kb mitochondrial genome have been implicated in numerous disorders that often display muscular and/or neurological symptoms due to the high-energy demands of these tissues. We describe a catalogue of 4489 putative mitochondrial DNA (mtDNA) deletions, including their frequency and relative read rate, using a combinatorial approach of mitochondria-targeted PCR, next-generation sequencing, bioinformatics, post-hoc filtering, annotation, and validation steps. Our bioinformatics pipeline uses MapSplice, an RNA-seq splice junction detection algorithm, to detect and quantify mtDNA deletion breakpoints rather than mRNA splices. Analyses of 93 samples from postmortem brain and blood found (i) the 4977 bp 'common deletion' was neither the most frequent deletion nor the most abundant; (ii) brain contained significantly more deletions than blood; (iii) many high frequency deletions were previously reported in MitoBreak, suggesting they are present at low levels in metabolically active tissues and are not exclusive to individuals with diagnosed mitochondrial pathologies; (iv) many individual deletions (and cumulative metrics) had significant and positive correlations with age and (v) the highest deletion burdens were observed in major depressive disorder brain, at levels greater than Kearns-Sayre Syndrome muscle. Collectively, these data suggest the Splice-Break pipeline can detect and quantify mtDNA deletions at a high level of resolution.

Uneven balance of power between hypothalamic peptidergic neurons in the control of feeding.

Two classes of peptide-producing neurons in the arcuate nucleus (Arc) of the hypothalamus are known to exert opposing actions on feeding: the anorexigenic neurons that express proopiomelanocortin (POMC) and the orexigenic neurons that express agouti-related protein (AgRP) and neuropeptide Y (NPY). These neurons are thought to arise from a common embryonic progenitor, but our anatomical and functional understanding of the interplay of these two peptidergic systems that contribute to the control of feeding remains incomplete. The present study uses a combination of optogenetic stimulation with viral and transgenic approaches, coupled with neural activity mapping and brain transparency visualization to demonstrate the following: (i) selective activation of Arc POMC neurons inhibits food consumption rapidly in unsated animals; (ii) activation of Arc neurons arising from POMC-expressing progenitors, including POMC and a subset of AgRP neurons, triggers robust feeding behavior, even in the face of satiety signals from POMC neurons; (iii) the opposing effects on food intake are associated with distinct neuronal projection and activation patterns of adult hypothalamic POMC neurons versus Arc neurons derived from POMC-expressing lineages; and (iv) the increased food intake following the activation of orexigenic neurons derived from POMC-expressing progenitors engages an extensive neural network that involves the endogenous opioid system. Together, these findings shed further light on the dynamic balance between two peptidergic systems in the moment-to-moment regulation of feeding behavior.

Adaptation to single housing is dynamic: Changes in hormone levels, gene expression, signaling in the brain, and anxiety-like behavior in adult male C57Bl/6J mice.

For basic research, rodents are often housed in individual cages prior to behavioral testing. However, aspects of the experimental design, such as duration of isolation and timing of animal manipulation, may unintentionally introduce variance into collected data. Thus, we examined temporal correlates of acclimation of C57Bl/6J mice to single housing in a novel environment following two commonly used experimental time periods (7 or 14 days, SH7 or SH14). We measured circulating stress hormones (adrenocorticotropic hormone and corticosterone), basally or after injection stress, hippocampal gene expression of transcripts implicated in stress and affect regulation: the glucocorticoid receptor (GR), the mineralocorticoid receptor (MR), including the MR/GR ratio, and fibroblast growth factor 2 (FGF2). We also measured signaling in the mammalian target of rapamycin (mTOR) pathway. The basal elevation of stress hormones in the SH14 group is accompanied by a blunting in the circadian rhythms of GR and FGF2 hippocampal gene expression, and the MR/GR ratio, that is observed in SH7 mice. Following mild stress, the endocrine response and hippocampal mTOR pathway signaling are decreased in the SH14 mice. These neural and endocrine changes at 14 days of single housing likely underlie increased anxiety-like behavior measured in an elevated plus maze test. We conclude that multiple measures of stress responsiveness change dynamically between one and two weeks of single housing. The ramifications of these alterations should be considered when designing animal experiments since such hidden sources of variance might cause lack of replicability and misinterpretation of data.

Genetic background and epigenetic modifications in the core of the nucleus accumbens predict addiction-like behavior in a rat model.

This study provides a demonstration in the rat of a clear genetic difference in the propensity for addiction-related behaviors following prolonged cocaine self-administration. It relies on the use of selectively bred high-responder (bHR) and low-responder (bLR) rat lines that differ in several characteristics associated with "temperament," including novelty-induced locomotion and impulsivity. We show that bHR rats exhibit behaviors reminiscent of human addiction, including persistent cocaine-seeking and increased reinstatement of cocaine seeking. To uncover potential underlying mechanisms of this differential vulnerability, we focused on the core of the nucleus accumbens and examined expression and epigenetic regulation of two transcripts previously implicated in bHR/bLR differences: fibroblast growth factor (FGF2) and the dopamine D2 receptor (D2). Relative to bHRs, bLRs had lower FGF2 mRNA levels and increased association of a repressive mark on histones (H3K9me3) at the FGF2 promoter. These differences were apparent under basal conditions and persisted even following prolonged cocaine self-administration. In contrast, bHRs had lower D2 mRNA under basal conditions, with greater association of H3K9me3 at the D2 promoter and these differences were no longer apparent following prolonged cocaine self-administration. Correlational analyses indicate that the association of H3K9me3 at D2 may be a critical substrate underlying the propensity to relapse. These findings suggest that low D2 mRNA levels in the nucleus accumbens core, likely mediated via epigenetic modifications, may render individuals more susceptible to cocaine addiction. In contrast, low FGF2 levels, which appear immutable even following prolonged cocaine exposure, may serve as a protective factor.

Dysregulated fibroblast growth factor (FGF) signaling in neurological and psychiatric disorders.

The role of the fibroblast growth factor (FGF) system in brain-related disorders has received considerable attention in recent years. To understand the role of this system in neurological and psychiatric disorders, it is important to identify the specific members of the FGF family that are implicated, their location and the various mechanisms they can be modulated. Each disorder appears to impact specific molecular players in unique anatomical locations, and all of these could conceivably become targets for treatment. In the last several years, the issue of how to target this system directly has become an area of increasing interest. To date, the most promising therapeutics are small molecule inhibitors and antibodies that modulate FGF receptor (FGFR) function. Beyond attempting to modify the primary players affected by a given brain disorder, it may prove useful to target molecules, such as membrane-bound or extracellular proteins that interact with FGF ligands or FGFRs to modulate signaling.

Fibroblast growth factor 9 is a novel modulator of negative affect.

Both gene expression profiling in postmortem human brain and studies using animal models have implicated the fibroblast growth factor (FGF) family in affect regulation and suggest a potential role in the pathophysiology of major depressive disorder (MDD). FGF2, the most widely characterized family member, is down-regulated in the depressed brain and plays a protective role in rodent models of affective disorders. By contrast, using three microarray analyses followed by quantitative RT-PCR confirmation, we show that FGF9 expression is up-regulated in the hippocampus of individuals with MDD, and that FGF9 expression is inversely related to the expression of FGF2. Because little is known about FGF9's function in emotion regulation, we used animal models to shed light on its potential role in affective function. We found that chronic social defeat stress, an animal model recapitulating some aspects of MDD, leads to a significant increase in hippocampal FGF9 expression, paralleling the elevations seen in postmortem human brain tissue. Chronic intracerebroventricular administration of FGF9 increased both anxiety- and depression-like behaviors. In contrast, knocking down FGF9 expression in the dentate gyrus of the hippocampus using a lentiviral vector produced a decrease in FGF9 expression and ameliorated anxiety-like behavior. Collectively, these results suggest that high levels of hippocampal FGF9 play an important role in the development or expression of mood and anxiety disorders. We propose that the relative levels of FGF9 in relation to other members of the FGF family may prove key to understanding vulnerability or resilience in affective disorders.

FGF2 is a target and a trigger of epigenetic mechanisms associated with differences in emotionality: partnership with H3K9me3.

Posttranslational modifications of histone tails in chromatin template can result from environmental experiences such as stress and substance abuse. However, the role of epigenetic modifications as potential predisposing factors in affective behavior is less well established. To address this question, we used our selectively bred lines of high responder (bHR) and low responder (bLR) rats that show profound and stable differences in affective responses, with bLRs being prone to anxiety- and depression-like behavior and bHRs prone to addictive behavior. We first asked whether these phenotypes are associated with basal differences in epigenetic profiles. Our results reveal broad between-group differences in basal levels of trimethylated histone protein H3 at lysine 9 (H3K9me3) in hippocampus (HC), amygdala, and nucleus accumbens. Moreover, levels of association of H3K9me3 at Glucocorticoid Receptor (GR) and Fibroblast growth Factor 2 (FGF2) promoters differ reciprocally between bHRs and bLRs in these regions, consistent with these genes' opposing levels of expression and roles in modulating anxiety behavior. Importantly, this basal epigenetic pattern is modifiable by FGF2, a factor that modulates anxiety behavior. Thus, early-life FGF2, which decreases anxiety, altered the levels of H3K9me3 and its binding at FGF2 and GR promoters of bLRs rendering them more similar to bHRs. Conversely, knockdown of HC FGF2 altered both anxiety behavior and levels of H3K9me3 in bHRs, rendering them more bLR-like. These findings implicate FGF2 as a modifier of epigenetic mechanisms associated with emotional responsiveness, and point to H3K9me3 as a key player in the regulation of affective vulnerability.

Protective effects of chronic mild stress during adolescence in the low-novelty responder rat.

Stress-elicited behavioral and physiologic responses vary widely across individuals and depend on a combination of environmental and genetic factors. Adolescence is an important developmental period when neural circuits that guide emotional behavior and stress reactivity are still maturing. A critical question is whether stress exposure elicits contrasting effects when it occurs during adolescence versus adulthood. We previously found that Sprague-Dawley rats selectively bred for low-behavioral response to novelty (bred Low Responders; bLRs) are particularly sensitive to chronic unpredictable mild stress (CMS) exposure in adulthood, which exacerbates their typically high levels of spontaneous depressive- and anxiety-like behavior. Given developmental processes known to occur during adolescence, we sought to determine whether the impact of CMS on bLR rats is equivalent when they are exposed to it during adolescence as compared with adulthood. Young bLR rats were either exposed to CMS or control condition from postnatal days 35-60. As adults, we found that CMS-exposed bLRs maintained high levels of sucrose preference and exhibited increased social exploration along with decreased immobility on the forced swim test compared with bLR controls. These data indicate a protective effect of CMS exposure during adolescence in bLR rats.

Fibroblast growth factor 2 alters the oxytocin receptor in a developmental model of anxiety-like behavior in male rat pups.

We aimed to determine the short-term effects of early-life stress in the form of maternal separation (MS) on anxiety-like behavior in male rat pups. In order to assess anxiety, we measured 40kHz separation-induced ultrasonic vocalizations (USV) on postnatal day (PND) 11. We further aimed to evaluate the potential involvement of two neurochemical systems known to regulate social and anxiety-like behaviors throughout life: oxytocin (OT) and fibroblast growth factor 2 (FGF2). For these purposes, we tested the effects of neonatal administration (on PND1) of an acute dose of FGF2 on USV and its potential interaction with MS. In addition, we validated the anxiolytic effects of OT and measured oxytocin receptor (OTR) gene expression, binding and epigenetic regulation via histone acetylation. Our results show that MS potentiated USV while acute administration of OT and FGF2 attenuated them. Further, we found that both FGF2 and MS increased OTR gene expression and the association of acH3K14 with the OTR promoter in the bed nucleus of the stria terminalis (BNST). Comparable changes, though not as pronounced, were also found for the central amygdala (CeA). Our findings suggest that FGF2 may exert its anxiolytic effects in male MS rats by a compensatory increase in the acetylation of the OTR promoter to overcome reduced OT levels in the BNST.

Circadian patterns of gene expression in the human brain and disruption in major depressive disorder.

A cardinal symptom of major depressive disorder (MDD) is the disruption of circadian patterns. However, to date, there is no direct evidence of circadian clock dysregulation in the brains of patients who have MDD. Circadian rhythmicity of gene expression has been observed in animals and peripheral human tissues, but its presence and variability in the human brain were difficult to characterize. Here, we applied time-of-death analysis to gene expression data from high-quality postmortem brains, examining 24-h cyclic patterns in six cortical and limbic regions of 55 subjects with no history of psychiatric or neurological illnesses ("controls") and 34 patients with MDD. Our dataset covered ~12,000 transcripts in the dorsolateral prefrontal cortex, anterior cingulate cortex, hippocampus, amygdala, nucleus accumbens, and cerebellum. Several hundred transcripts in each region showed 24-h cyclic patterns in controls, and >100 transcripts exhibited consistent rhythmicity and phase synchrony across regions. Among the top-ranked rhythmic genes were the canonical clock genes BMAL1(ARNTL), PER1-2-3, NR1D1(REV-ERBa), DBP, BHLHE40 (DEC1), and BHLHE41(DEC2). The phasing of known circadian genes was consistent with data derived from other diurnal mammals. Cyclic patterns were much weaker in the brains of patients with MDD due to shifted peak timing and potentially disrupted phase relationships between individual circadian genes. This transcriptome-wide analysis of the human brain demonstrates a rhythmic rise and fall of gene expression in regions outside of the suprachiasmatic nucleus in control subjects. The description of its breakdown in MDD suggests potentially important molecular targets for treatment of mood disorders.

Basal microRNA expression patterns in reward circuitry of selectively bred high-responder and low-responder rats vary by brain region and genotype.

Mental health disorders involving altered reward, emotionality, and anxiety are thought to result from the interaction of individual predisposition (genetic factors) and personal experience (environmental factors), although the mechanisms that contribute to an individual's vulnerability to these disorders remain poorly understood. We used an animal model of individual variation [inbred high-responder/low-responder (bHR-bLR) rodents] known to vary in reward, anxiety, and emotional processing to examine neuroanatomical expression patterns of microRNAs (miRNAs). Laser capture microdissection was used to dissect the prelimbic cortex and the nucleus accumbens core and shell prior to analysis of basal miRNA expression in bHR and bLR male rats. These studies identified 187 miRNAs differentially expressed by genotype in at least one brain region, 10 of which were validated by qPCR. Four of these 10 qPCR-validated miRNAs demonstrated differential expression across multiple brain regions, and all miRNAs with validated differential expression between genotypes had lower expression in bHR animals compared with bLR animals. microRNA (miR)-484 and miR-128a expression differences between the prelimbic cortex of bHR and bLR animals were validated by semiquantitative in situ hybridization. miRNA expression analysis independent of genotype identified 101 miRNAs differentially expressed by brain region, seven of which validated by qPCR. Dnmt3a mRNA, a validated target of miR-29b, varied in a direction opposite that of miR-29b's differential expression between bHR and bLR animals. These data provide evidence that basal central nervous system miRNA expression varies in the bHR-bLR model, implicating microRNAs as potential epigenetic regulators of key neural circuits and individual differences associated with mental health disorders.

High novelty-seeking rats are resilient to negative physiological effects of the early life stress.

Exposure to early life stress dramatically impacts adult behavior, physiology, and neuroendocrine function. Using rats bred for novelty-seeking differences and known to display divergent anxiety, depression, and stress vulnerability, we examined the interaction between early life adversity and genetic predisposition for high- versus low-emotional reactivity. Thus, bred Low Novelty Responder (bLR) rats, which naturally exhibit high anxiety- and depression-like behavior, and bred High Novelty Responder (bHR) rats, which show low anxiety/depression together with elevated aggression, impulsivity, and addictive behavior, were subjected to daily 3 h maternal separation (MS) stress postnatal days 1-14. We hypothesized that MS stress would differentially impact adult bHR/bLR behavior, physiology (stress-induced defecation), and neuroendocrine reactivity. While MS stress did not impact bHR and bLR anxiety-like behavior in the open field test and elevated plus maze, it exacerbated bLRs' already high physiological response to stress - stress-induced defecation. In both tests, MS bLR adult offspring showed exaggerated stress-induced defecation compared to bLR controls while bHR offspring were unaffected. MS also selectively impacted bLRs' (but not bHRs') neuroendocrine stress reactivity, producing an exaggerated corticosterone acute stress response in MS bLR versus control bLR rats. These findings highlight how genetic predisposition shapes individuals' response to early life stress. Future work will explore neural mechanisms underlying the distinct behavioral and neuroendocrine consequences of MS in bHR/bLR animals.

Fibroblast growth factor-2 (FGF2) augmentation early in life alters hippocampal development and rescues the anxiety phenotype in vulnerable animals.

Individuals with mood disorders exhibit alterations in the fibroblast growth factor system, including reduced hippocampal fibroblast growth factor-2 (FGF2). It is difficult, however, to pinpoint whether these alterations are a cause or consequence of the disorder. The present study asks whether FGF2 administered the day after birth has long-lasting effects on hippocampal development and emotionality. We show that early-life FGF2 shifts the pace of neurogenesis, with an early acceleration around weaning followed by a deceleration in adulthood. This, in turn, results in a denser dentate gyrus with more neurons. To assess the impact of early-life FGF2 on emotionality, we use rats selectively bred for differences in locomotor response to novelty. Selectively bred low-responder (bLR) rats show low levels of novelty-induced locomotion and exhibit high levels of anxiety- and depression-like behavior compared with their selectively bred high-responder counterparts. Early-life FGF2 decreased anxiety-like behavior in highly anxious bLRs without altering other behaviors and without affecting high-responder rats. Laser capture microscopy of the dentate gyrus followed by microarray analysis revealed genes that were differentially expressed in bLRs exposed to early-life FGF2 vs. vehicle-treated bLRs. Some of the differentially expressed genes that have been positively associated with anxiety were down-regulated, whereas genes that promote cell survival were up-regulated. Overall, these results show a key role for FGF2 in the developmental trajectory of the hippocampus as well as the modulation of anxiety-like behavior in adulthood, and they point to potential downstream targets for the treatment of anxiety disorders.

Resource: A curated database of brain-related functional gene sets (Brain.GMT).

Transcriptional profiling has become a common tool for investigating the nervous system. During analysis, differential expression results are often compared to functional ontology databases, which contain curated gene sets representing well-studied pathways. This dependence can cause neuroscience studies to be interpreted in terms of functional pathways documented in better studied tissues (e.g., liver) and topics (e.g., cancer), and systematically emphasizes well-studied genes, leaving other findings in the obscurity of the brain "ignorome". To address this issue, we compiled a curated database of 918 gene sets related to nervous system function, tissue, and cell types ("Brain.GMT") that can be used within common analysis pipelines (GSEA, limma, edgeR) to interpret results from three species (rat, mouse, human). Brain.GMT includes brain-related gene sets curated from the Molecular Signatures Database (MSigDB) and extracted from public databases (GeneWeaver, Gemma, DropViz, BrainInABlender, HippoSeq) and published studies containing differential expression results. Although Brain.GMT is still undergoing development and currently only represents a fraction of available brain gene sets, "brain ignorome" genes are already better represented than in traditional Gene Ontology databases. Moreover, Brain.GMT substantially improves the quantity and quality of gene sets identified as enriched with differential expression in neuroscience studies, enhancing interpretation. •We compiled a curated database of 918 gene sets related to nervous system function, tissue, and cell types ("Brain.GMT").•Brain.GMT can be used within common analysis pipelines (GSEA, limma, edgeR) to interpret neuroscience transcriptional profiling results from three species (rat, mouse, human).•Although Brain.GMT is still undergoing development, it substantially improved the interpretation of differential expression results within our initial use cases.

Resource: A Curated Database of Brain-Related Functional Gene Sets (Brain.GMT).

Transcriptional profiling has become a common tool for investigating the nervous system. During analysis, differential expression results are often compared to functional ontology databases, which contain curated gene sets representing well-studied pathways. This dependence can cause neuroscience studies to be interpreted in terms of functional pathways documented in better studied tissues (e.g., liver) and topics (e.g., cancer), and systematically emphasizes well-studied genes, leaving other findings in the obscurity of the brain "ignorome". To address this issue, we compiled a curated database of 918 gene sets related to nervous system function, tissue, and cell types ("Brain.GMT") that can be used within common analysis pipelines (GSEA, limma, edgeR) to interpret results from three species (rat, mouse, human). Brain.GMT includes brain-related gene sets curated from the Molecular Signatures Database (MSigDB) and extracted from public databases (GeneWeaver, Gemma, DropViz, BrainInABlender, HippoSeq) and published studies containing differential expression results. Although Brain.GMT is still undergoing development and currently only represents a fraction of available brain gene sets, "brain ignorome" genes are already better represented than in traditional Gene Ontology databases. Moreover, Brain.GMT substantially improves the quantity and quality of gene sets identified as enriched with differential expression in neuroscience studies, enhancing interpretation.

Adolescent environmental enrichment induces social resilience and alters neural gene expression in a selectively bred rodent model with anxious phenotype.

Stress is a major influence on mental health status; the ways that individuals respond to or copes with stressors determine whether they are negatively affected in the future. Stress responses are established by an interplay between genetics, environment, and life experiences. Psychosocial stress is particularly impactful during adolescence, a critical period for the development of mood disorders. In this study we compared two established, selectively-bred Sprague Dawley rat lines, the "internalizing" bred Low Responder (bLR) line versus the "externalizing" bred High Responder (bHR) line, to investigate how genetic temperament and adolescent environment impact future responses to social interactions and psychosocial stress, and how these determinants of stress response interact. Male bLR and bHR rats were exposed to social and environmental enrichment in adolescence prior to experiencing social defeat and were then assessed for social interaction and anxiety-like behavior. Adolescent enrichment caused rats to display more social interaction, as well as nominally less social avoidance, less submission during defeat, and resilience to the effects of social stress on corticosterone, in a manner that seemed more notable in bLRs. For bHRs, enrichment also caused greater aggression during a neutral social encounter and nominally during defeat, and decreased anxiety-like behavior. To explore the neurobiology underlying the development of social resilience in the anxious phenotype bLRs, RNA-seq was conducted on the hippocampus and nucleus accumbens, two brain regions that mediate stress regulation and social behavior. Gene sets previously associated with stress, social behavior, aggression and exploratory activity were enriched with differential expression in both regions, with a particularly large effect on gene sets that regulate social behaviors. Our findings provide further evidence that adolescent enrichment can serve as an inoculating experience against future stressors. The ability to induce social resilience in a usually anxious line of animals by manipulating their environment has translational implications, as it underscores the feasibility of intervention strategies targeted at genetically vulnerable adolescent populations.