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.

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.

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.

High emotional reactivity is associated with activation of a molecularly distinct hippocampal-amygdala circuit modulated by the glucocorticoid receptor

Emotions are characterized not only by their valence but also by whether they are stable or labile. Yet, we do not understand the molecular or circuit mechanisms that control the dynamic nature of emotional responses. We have shown that glucocorticoid receptor overexpression in the forebrain (GRov) leads to a highly reactive mouse with increased anxiety behavior coupled with greater swings in emotional responses. This phenotype is established early in development and persists into adulthood. However, the neural circuitry mediating this lifelong emotional lability remains unknown. In the present study, optogenetic stimulation in ventral dentate gyrus (vDG) of GRov mice led to a greater range and a prolonged duration of anxiety behavior. cFos expression analysis showed that the amplified behavioral response to vDG activation in GRov mice is coupled to increased neuronal activity in specific brain regions. Relative to wild type mice, GRov mice displayed glutamatergic/GABAergic activation imbalance in ventral CA1 (vCA1) and selectively increased glutamatergic activation in the basal posterior amygdaloid complex. Moreover, forebrain GR overexpression led to increased activation of molecularly distinct subpopulations of neurons within the hippocampus and the posterior basolateral amygdala (pBLA) as evident from the increased cFos co-labeling in the calbindin1+ glutamatergic neurons in vCA1 and in the DARPP-32/Ppp1r1b+ glutamatergic neurons in pBLA. We propose that a molecularly distinct hippocampal-amygdala circuit is shaped by stress early in life and tunes the dynamics of emotional responses.

Allele biased transcription factor binding across human brain regions gives mechanistic insight into eQTLs

Transcription Factors (TFs) influence gene expression by facilitating or disrupting the formation of transcription initiation machinery at particular genomic loci. Because genomic localization of TFs is in part driven by TF recognition of DNA sequence, variation in TF binding sites can disrupt TF-DNA associations and affect gene regulation. To identify variants that impact TF binding in human brain tissues, we quantified allele bias for 93 TFs analyzed with ChIP-seq experiments of multiple structural brain regions from two donors. Using graph genomes constructed from phased genomic sequence data, we compared ChIP-seq signal between alleles at heterozygous variants within each tissue sample from each donor. Comparison of results from different brain regions within donors and the same regions between donors provided measures of allele bias reproducibility. We identified thousands of DNA variants that show reproducible bias in ChIP-seq for at least one TF. We found that alleles that are rarer in the general population were more likely than common alleles to exhibit large biases, and more frequently led to reduced TF binding. Combining ChIP-seq with RNA-seq, we identified TF-allele interaction biases with RNA bias in a phased allele linked to 6,709 eQTL variants identified in GTEx data, 3,309 of which were found in neural contexts. Our results provide insights into the effects of both common and rare variation on gene regulation in the brain. These findings can facilitate mechanistic understanding of cis-regulatory variation associated with biological traits, including disease.

Neurotransmission-related gene expression in the frontal pole is altered in subjects with bipolar disorder and schizophrenia

The frontal pole (Brodmann area 10, BA10) is the largest cytoarchitectonic region of the human cortex, performing complex integrative functions. BA10 undergoes intensive adolescent grey matter pruning prior to the age of onset for bipolar disorder (BP) and schizophrenia (SCHIZ), and its dysfunction is likely to underly aspects of their shared symptomology. In this study, we investigated the role of BA10 neurotransmission-related gene expression in BP and SCHIZ. We performed qPCR to measure the expression of 115 neurotransmission-related targets in control, BP, and SCHIZ postmortem samples (n = 72). We chose this method for its high sensitivity to detect low-level expression. We then strengthened our findings by performing a meta-analysis of publicly released BA10 microarray data (n = 101) and identified sources of convergence with our qPCR results. To improve interpretation, we leveraged the unusually large database of clinical metadata accompanying our samples to explore the relationship between BA10 gene expression, therapeutics, substances of abuse, and symptom profiles, and validated these findings with publicly available datasets. Using these convergent sources of evidence, we identified 20 neurotransmission-related genes that were differentially expressed in BP and SCHIZ in BA10. These results included a large diagnosis-related decrease in two important therapeutic targets with low levels of expression, HTR2B and DRD4, as well as other findings related to dopaminergic, GABAergic and astrocytic function. We also observed that therapeutics may produce a differential expression that opposes diagnosis effects. In contrast, substances of abuse showed similar effects on BA10 gene expression as BP and SCHIZ, potentially amplifying diagnosis-related dysregulation.

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.

Single nucleus multiomics identifies ZEB1 and MAFB as candidate regulators of Alzheimer's disease-specific cis-regulatory elements

Cell type-specific transcriptional differences between brain tissues from donors with Alzheimer's disease (AD) and unaffected controls have been well documented, but few studies have rigorously interrogated the regulatory mechanisms responsible for these alterations. We performed single nucleus multiomics (snRNA-seq plus snATAC-seq) on 105,332 nuclei isolated from cortical tissues from 7 AD and 8 unaffected donors to identify candidate cis-regulatory elements (CREs) involved in AD-associated transcriptional changes. We detected 319,861 significant correlations, or links, between gene expression and cell type-specific transposase accessible regions enriched for active CREs. Among these, 40,831 were unique to AD tissues. Validation experiments confirmed the activity of many regions, including several candidate regulators of APP expression. We identified ZEB1 and MAFB as candidate transcription factors playing important roles in AD-specific gene regulation in neurons and microglia, respectively. Microglia links were globally enriched for heritability of AD risk and previously identified active regulatory regions.

Genome-wide association study in a rat model of temperament identifies multiple loci for exploratory locomotion and anxiety-like traits

Common genetic factors likely contribute to multiple psychiatric diseases including mood and substance use disorders. Certain stable, heritable traits reflecting temperament, termed externalizing or internalizing, play a large role in modulating vulnerability to these disorders. To model these heritable tendencies, we selectively bred rats for high and low exploration in a novel environment [bred High Responders (bHR) vs. Low Responders (bLR)]. To identify genes underlying the response to selection, we phenotyped and genotyped 538 rats from an F₂ cross between bHR and bLR. Several behavioral traits show high heritability, including the selection trait: exploratory locomotion (EL) in a novel environment. There were significant phenotypic and genetic correlations between tests that capture facets of EL and anxiety. There were also correlations with Pavlovian conditioned approach (PavCA) behavior despite the lower heritability of that trait. Ten significant and conditionally independent loci for six behavioral traits were identified. Five of the six traits reflect different facets of EL that were captured by three behavioral tests. Distance traveled measures from the open field and the elevated plus maze map onto different loci, thus may represent different aspects of novelty-induced locomotor activity. The sixth behavioral trait, number of fecal boli, is the only anxiety-related trait mapping to a significant locus on chromosome 18 within which the Pik3c3 gene is located. There were no significant loci for PavCA. We identified a missense variant in the Plekhf1 gene on the chromosome 1:95 Mb QTL and Fancf and Gas2 as potential candidate genes that may drive the chromosome 1:107 Mb QTL for EL traits. The identification of a locomotor activity-related QTL on chromosome 7 encompassing the Pkhd1l1 and Trhr genes is consistent with our previous finding of these genes being differentially expressed in the hippocampus of bHR vs. bLR rats. The strong heritability coupled with identification of several loci associated with exploratory locomotion and emotionality provide compelling support for this selectively bred rat model in discovering relatively large effect causal variants tied to elements of internalizing and externalizing behaviors inherent to psychiatric and substance use disorders.

Large Common Mitochondrial DNA Deletions Are Associated with a Mitochondrial SNP T14798C Near the 3' Breakpoints

Introduction

Large somatic deletions of mitochondrial DNA (mtDNA) accumulate with aging in metabolically active tissues such as the brain. We have cataloged the breakpoints and frequencies of large mtDNA deletions in the human brain.

Methods

We quantified 112 high-frequency mtDNA somatic deletions across four human brain regions with the Splice-Break2 pipeline. In addition, we utilized PLINK/Seq to test the association of mitochondrial genotypes with the abundance of these high-frequency mtDNA deletions. A conservative p value threshold of 5E-08 was used to find the significant loci.

Results

One mtDNA SNP (T14798C) was significantly associated with mtDNA deletions in two brain regions, the dorsolateral prefrontal cortex (DLPFC) and the superior temporal gyrus. Since the DLPFC showed the most robust association between T14798C and two deletion breakpoints (7816-14807 and 5462-14807), this association was tested in the DLPFC of a replication sample and validated the first results. Incorporating the C allele at 14,798 bp increased the perfect/imperfect length of the repeat at the 3' breakpoint of the two associated deletions.

Conclusion

This is the first study to identify the association of mtDNA SNP with large mtDNA deletions in the human brain. The T14798C allele located in the MT-CYB gene is a common polymorphism that occurs in several mitochondrial haplogroups. We hypothesize that the T14798C association with two deletions occurs by extending the repeat length around the 3' deletion breakpoints. This simple mechanism suggests that mtDNA SNPs can affect the mitochondrial genome structure, especially in brain where high levels of reactive oxygen species lead to deletion accumulation with aging.

Mitochondria DNA copy number, mitochondria DNA total somatic deletions, Complex I activity, synapse number, and synaptic mitochondria number are altered in schizophrenia and bipolar disorder

Mitochondrial dysfunction is a neurobiological phenomenon implicated in the pathophysiology of schizophrenia and bipolar disorder that can synergistically affect synaptic neurotransmission. We hypothesized that schizophrenia and bipolar disorder share molecular alterations at the mitochondrial and synaptic levels. Mitochondria DNA (mtDNA) copy number (CN), mtDNA common deletion (CD), mtDNA total deletion, complex I activity, synapse number, and synaptic mitochondria number were studied in the postmortem human dorsolateral prefrontal cortex (DLPFC), superior temporal gyrus (STG), primary visual cortex (V1), and nucleus accumbens (NAc) of controls (CON), and subjects with schizophrenia (SZ), and bipolar disorder (BD). The results showed (i) the mtDNA CN is significantly higher in DLPFC of both SZ and BD, decreased in the STG of BD, and unaltered in V1 and NAc of both SZ and BD; (ii) the mtDNA CD is significantly higher in DLPFC of BD while unaltered in STG, V1, and NAc of both SZ and BD; (iii) The total deletion burden is significantly higher in DLPFC in both SZ and BD while unaltered in STG, V1, and NAc of SZ and BD; (iv) Complex I activity is significantly lower in DLPFC of both SZ and BD, which is driven by the presence of medications, with no alteration in STG, V1, and NAc. In addition, complex I protein concentration, by ELISA, was decreased across three cortical regions of SZ and BD subjects; (v) The number of synapses is decreased in DLPFC of both SZ and BD, while the synaptic mitochondria number was significantly lower in female SZ and female BD compared to female controls. Overall, these findings will pave the way to understand better the pathophysiology of schizophrenia and bipolar disorder for therapeutic interventions.

Identification of potential blood biomarkers associated with suicide in major depressive disorder

Suicides have increased to over 48,000 deaths yearly in the United States. Major depressive disorder (MDD) is the most common diagnosis among suicides, and identifying those at the highest risk for suicide is a pressing challenge. The objective of this study is to identify changes in gene expression associated with suicide in brain and blood for the development of biomarkers for suicide. Blood and brain were available for 45 subjects (53 blood samples and 69 dorsolateral prefrontal cortex (DLPFC) samples in total). Samples were collected from MDD patients who died by suicide (MDD-S), MDDs who died by other means (MDD-NS) and non-psychiatric controls. We analyzed gene expression using RNA and the NanoString platform. In blood, we identified 14 genes which significantly differentiated MDD-S versus MDD-NS. The top six genes differentially expressed in blood were: PER3, MTPAP, SLC25A26, CD19, SOX9, and GAR1. Additionally, four genes showed significant changes in brain and blood between MDD-S and MDD-NS; SOX9 was decreased and PER3 was increased in MDD-S in both tissues, while CD19 and TERF1 were increased in blood but decreased in DLPFC. To our knowledge, this is the first study to analyze matched blood and brain samples in a well-defined population of MDDs demonstrating significant differences in gene expression associated with completed suicide. Our results strongly suggest that blood gene expression is highly informative to understand molecular changes in suicide. Developing a suicide biomarker signature in blood could help health care professionals to identify subjects at high risk for suicide.

Rare coding variants in ten genes confer substantial risk for schizophrenia

Rare coding variation has historically provided the most direct connections between gene function and disease pathogenesis. By meta-analysing the whole exomes of 24,248 schizophrenia cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in 10 genes as conferring substantial risk for schizophrenia (odds ratios of 3-50, P < 2.14 × 10-6) and 32 genes at a false discovery rate of <5%. These genes have the greatest expression in central nervous system neurons and have diverse molecular functions that include the formation, structure and function of the synapse. The associations of the NMDA (N-methyl-D-aspartate) receptor subunit GRIN2A and AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptor subunit GRIA3 provide support for dysfunction of the glutamatergic system as a mechanistic hypothesis in the pathogenesis of schizophrenia. We observe an overlap of rare variant risk among schizophrenia, autism spectrum disorders1, epilepsy and severe neurodevelopmental disorders2, although different mutation types are implicated in some shared genes. Most genes described here, however, are not implicated in neurodevelopment. We demonstrate that genes prioritized from common variant analyses of schizophrenia are enriched in rare variant risk3, suggesting that common and rare genetic risk factors converge at least partially on the same underlying pathogenic biological processes. Even after excluding significantly associated genes, schizophrenia cases still carry a substantial excess of URVs, which indicates that more risk genes await discovery using this approach.

Differences in microglia morphological profiles reflect divergent emotional temperaments: insights from a selective breeding model

Microglia play critical roles in healthy brain development and function, as well as the neuropathology underlying a range of brain diseases. Despite evidence for a role of microglia in affective regulation and mood disorders, little is known regarding how variation in microglia status relates to individual differences in emotionality. Using a selective breeding model, we have generated rat lines with unique temperamental phenotypes that reflect broad emotional traits: bred low responder rats (bLRs) are novelty-averse and model a passive coping style, whereas bred high responder rats (bHRs) are highly exploratory and model an active coping style. To identify a functional role of microglia in these phenotypes, we administered minocycline, an antibiotic with potent microglia inhibiting properties and observed shifts in forced swim, sucrose preference, and social interaction behaviors in bLRs. Using detailed anatomical analyses, we compared hippocampal microglia profiles of bHRs and bLRs and found that although the lines had similar numbers of microglia, selective breeding was associated with a shift in the morphological features of these cells. Specifically, microglia from bLRs were characterized by a hyper-ramified morphology, with longer processes and more complicated branching patterns than microglia from bHRs. This morphology is thought to reflect an early stage of microglia activation and suggests that bLR microglia are in a reactive state even when animals are not overtly challenged. Taken together, our results provide novel evidence linking variation in inborn temperament with differences in the baseline status of microglia and implicate a role for microglia in shaping enduring emotional characteristics.

Characterizing the behavioral and neuroendocrine features of susceptibility and resilience to social stress

Evaluating and coping with stressful social events as they unfold is a critical strategy in overcoming them without long-lasting detrimental effects. Individuals display a wide range of responses to stress, which can manifest in a variety of outcomes for the brain as well as subsequent behavior. Chronic Social Defeat Stress (CSDS) in mice has been widely used to model individual variation following a social stressor. Following a course of repeated intermittent psychological and physical stress, mice diverge into separate populations of social reactivity: resilient (socially interactive) and susceptible (socially avoidant) animals. A rich body of work reveals distinct neurobiological and behavioral consequences of this experience that map onto the resilient and susceptible groups. However, the range of factors that emerge over the course of defeat have not been fully described. Therefore, in the current study, we focused on characterizing behavioral, physiological, and neuroendocrine profiles of mice in three separate phases: before, during, and following CSDS. We found that following CSDS, traditional read-outs of anxiety-like and depression-like behaviors do not map on to the resilient and susceptible groups. By contrast, behavioral coping strategies used during the initial social stress encounter better predict which mice will eventually become resilient or susceptible. In particular, mice that will emerge as susceptible display greater escape behavior on Day 1 of social defeat than those that will emerge as resilient, indicating early differences in coping mechanisms used between the two groups. We further show that the social avoidance phenotype in susceptible mice is specific to the aggressor strain and does not generalize to conspecifics or other strains, indicating that there may be features of threat discrimination that are specific to the susceptible mice. Our findings suggest that there are costs and benefits to both the resilient and susceptible outcomes, reflected in their ability to cope and adapt to the social stressor.

Adolescent cocaine differentially impacts psychomotor sensitization and epigenetic profiles in adult male rats with divergent affective phenotypes

Adolescent drug use reliably predicts increased addiction liability in adulthood, but not all individuals are equally impacted. To explore the biological bases of this differential reactivity to early life drug experience, we used a genetic rat model of temperament and evaluated the impact of adolescent cocaine exposure on adult psychomotor sensitization. Relative to adult bred low-responder (bLR) rats, bred high-responders (bHR) are more sensitive to the psychomotor-activating effects of cocaine and reinstate drug-seeking behavior more readily following prolonged cocaine exposure and/or abstinence. We found that a 7-day sensitizing cocaine regimen (15 mg/kg/day) during either adolescence or adulthood produced psychomotor sensitization in bHRs only, while a dual cocaine exposure prevented further sensitization, suggesting limits on neuroplasticity. By contrast, adolescent cocaine in bLRs shifted their resilient phenotype, rendering them more responsive to cocaine in adulthood following adolescent cocaine. To begin to explore the neural correlates of these behavioral phenotypes, we assessed two functionally opposite epigenetic chromatin modifications implicated in addiction liability, permissive acetylation (ac) and repressive tri-methylation (me3) on Histone 3 Lysine 9 (H3K9), in four striatal sub-regions. In bHRs, decreased H3K9me3 and increased acH3K9 in the nucleus accumbens (NAc) core associated with cocaine sensitization. In bLRs, the combination of cocaine exposure in adolescence and adulthood, which lead to an increased response to a cocaine challenge, also increased acH3K9 in the core. Thus, adolescent cocaine experience interacts with genetic background to elicit different behavioral profiles relevant to addiction in adulthood, with concurrent modifications in the epigenetic histone profiles in the NAc that associate with cocaine sensitization and with metaplasticity.

Genetic Liability for Internalizing Versus Externalizing Behavior Manifests in the Developing and Adult Hippocampus: Insight From a Meta-analysis of Transcriptional Profiling Studies in a Selectively Bred Rat Model

BACKGROUND

For more than 16 years, we have selectively bred rats for either high or low levels of exploratory activity within a novel environment. These bred high-responder (bHR) and bred low-responder (bLR) rats model temperamental extremes, exhibiting large differences in internalizing and externalizing behaviors relevant to mood and substance use disorders.

METHODS

We characterized persistent differences in gene expression related to bHR/bLR phenotype across development and adulthood in the hippocampus, a region critical for emotional regulation, by meta-analyzing 8 transcriptional profiling datasets (microarray and RNA sequencing) spanning 43 generations of selective breeding (postnatal day 7: n = 22; postnatal day 14: n = 49; postnatal day 21: n = 21; adult: n = 46; all male). We cross-referenced expression differences with exome sequencing within our colony to pinpoint candidates likely to mediate the effect of selective breeding on behavioral phenotype. The results were compared with hippocampal profiling from other bred rat models.

RESULTS

Genetic and transcriptional profiling results converged to implicate multiple candidate genes, including two previously associated with metabolism and mood: Trhr and Ucp2. Results also highlighted bHR/bLR functional differences in the hippocampus, including a network essential for neurodevelopmental programming, proliferation, and differentiation, centering on Bmp4 and Mki67. Finally, we observed differential expression related to microglial activation, which is important for synaptic pruning, including 2 genes within implicated chromosomal regions: C1qa and Mfge8.

CONCLUSIONS

These candidate genes and functional pathways may direct bHR/bLR rats along divergent developmental trajectories and promote a widely different reactivity to the environment.

Optimization and evaluation of fluorescence in situ hybridization chain reaction in cleared fresh-frozen brain tissues

Transcript labeling in intact tissues using in situ hybridization chain reaction has potential to provide vital spatiotemporal information for molecular characterization of heterogeneous neuronal populations. However, large tissue labeling in non-perfused or fresh-frozen rodent and postmortem human samples, which provide more flexible utilization than perfused tissues, is largely unexplored. In the present study, we optimized the combination of in situ hybridization chain reaction in fresh-frozen rodent brains and then evaluated the uniformity of neuronal labeling between two clearing methods, CLARITY and iDISCO⁺. We found that CLARITY yielded higher signal-to-noise ratios but more limited imaging depth and required longer clearing times, whereas, iDISCO⁺ resulted in better tissue clearing, greater imaging depth and a more uniform labeling of larger samples. Based on these results, we used iDISCO⁺-cleared fresh-frozen rodent brains to further validate this combination and map the expression of a few genes of interest pertaining to mood disorders. We then examined the potential of in situ hybridization chain reaction to label transcripts in cleared postmortem human brain tissues. The combination failed to produce adequate mRNA labeling in postmortem human cortical slices but produced visually adequate labeling in the cerebellum tissues. We next, investigated the multiplexing ability of in situ hybridization chain reaction in cleared tissues which revealed inconsistent fluorescence output depending upon the fluorophore conjugated to the hairpins. Finally, we applied our optimized protocol to assess the effect of glucocorticoid receptor overexpression on basal somatostatin expression in the mouse cortex. The constitutive glucocorticoid receptor overexpression resulted in lower number density of somatostatin-expressing neurons compared to wild type. Overall, the combination of in situ hybridization chain reaction with clearing methods, especially iDISCO⁺, may find broad application in the transcript analysis in rodent studies, but its limited use in postmortem human tissues can be improved by further optimizations.

Sepsis survivor mice exhibit a behavioral endocrine syndrome with ventral hippocampal dysfunction

Severe acute stressors are known to trigger mood disorders in humans. Sepsis represents one such stressor, and survivors often suffer long term from psychiatric morbidity. We hypothesized that sepsis leads to lasting changes in neural circuits involved in stress integration, altering affective behavior and the stress response. To investigate this hypothesis, sepsis was induced in male C57Bl/6 mice using cecal ligation and puncture (CLP), and control mice underwent sham surgery. Mice recovered from acute illness within 2 weeks, after which they exhibited increased avoidance behavior and behavioral despair compared with sham, with behavioral changes observed more than 5 weeks after recovery. Sepsis survivors also showed evidence of enhanced hypothalamic-pituitary-adrenal (HPA) axis activity, with increased corticosterone after a novel stressor and increased adrenal weight. In the brain, sepsis survivor mice showed decreased stress-induced cfos mRNA and increased glucocorticoid receptor immunoreactivity specifically in the ventral hippocampus, a brain region known to coordinate emotional behavior and HPA axis activity. We conclude that murine sepsis survivors exhibit a behavioral neuroendocrine syndrome of negative affective behavior and HPA axis hyperactivity, which could be explained by ventral hippocampal dysfunction. These findings could contribute to our understanding of the human post-intensive care syndrome.

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.

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.

Genome-wide association study identifies 30 loci associated with bipolar disorder

Bipolar disorder is a highly heritable psychiatric disorder. We performed a genome-wide association study (GWAS) including 20,352 cases and 31,358 controls of European descent, with follow-up analysis of 822 variants with P < 1 × 10-4 in an additional 9,412 cases and 137,760 controls. Eight of the 19 variants that were genome-wide significant (P < 5 × 10-8) in the discovery GWAS were not genome-wide significant in the combined analysis, consistent with small effect sizes and limited power but also with genetic heterogeneity. In the combined analysis, 30 loci were genome-wide significant, including 20 newly identified loci. The significant loci contain genes encoding ion channels, neurotransmitter transporters and synaptic components. Pathway analysis revealed nine significantly enriched gene sets, including regulation of insulin secretion and endocannabinoid signaling. Bipolar I disorder is strongly genetically correlated with schizophrenia, driven by psychosis, whereas bipolar II disorder is more strongly correlated with major depressive disorder. These findings address key clinical questions and provide potential biological mechanisms for bipolar disorder.

Cognitive Control as a 5-HT1A-Based Domain That Is Disrupted in Major Depressive Disorder

Heterogeneity within Major Depressive Disorder (MDD) has hampered identification of biological markers (e.g., intermediate phenotypes, IPs) that might increase risk for the disorder or reflect closer links to the genes underlying the disease process. The newer characterizations of dimensions of MDD within Research Domain Criteria (RDoC) domains may align well with the goal of defining IPs. We compare a sample of 25 individuals with MDD compared to 29 age and education matched controls in multimodal assessment. The multimodal RDoC assessment included the primary IP biomarker, positron emission tomography (PET) with a selective radiotracer for 5-HT_1A [(11C)WAY-100635], as well as event-related functional MRI with a Go/No-go task targeting the Cognitive Control network, neuropsychological assessment of affective perception, negative memory bias and Cognitive Control domains. There was also an exploratory genetic analysis with the serotonin transporter (5-HTTLPR) and monamine oxidase A (MAO-A) genes. In regression analyses, lower 5-HT_1A binding potential (BP) in the MDD group was related to diminished engagement of the Cognitive Control network, slowed resolution of interfering cognitive stimuli, one element of Cognitive Control. In contrast, higher/normative levels of 5-HT_1A BP in MDD (only) was related to a substantial memory bias toward negative information, but intact resolution of interfering cognitive stimuli and greater engagement of Cognitive Control circuitry. The serotonin transporter risk allele was associated with lower 1a BP and the corresponding imaging and cognitive IPs in MDD. Lowered 5HT 1a BP was present in half of the MDD group relative to the control group. Lowered 5HT 1a BP may represent a subtype including decreased engagement of Cognitive Control network and impaired resolution of interfering cognitive stimuli. Future investigations might link lowered 1a BP to neurobiological pathways and markers, as well as probing subtype-specific treatment targets.

Effects of early-life FGF2 on ultrasonic vocalizations (USVs) and the mu-opioid receptor in male Sprague-Dawley rats selectively-bred for differences in their response to novelty

In previous studies, early-life fibroblast growth factor-2 (FGF2) administration conferred resilience to developing anxiety-like behavior in vulnerable animals in adulthood. To follow up on this work, we administered FGF2 the day after birth to animals that differ in emotional behavior and further explored its long-term effects on affective behavior and circuitry. Selectively-bred "high responder" rats (bHRs) exhibit low levels of anxiety-like and depression-like behavior, whereas selectively-bred "low responders" (bLRs) display high levels of anxiety-like and depression-like behavior. We found that early-life administration of FGF2 decreased negative affect in bLRs during the early post-natal period, as indexed by 40 kHz ultrasonic vocalizations (USVs) in response to a brief maternal separation on PND11. FGF2 also increased positive affect during the juvenile period, as measured by 50 kHz USVs in response to heterospecific hand play ("tickling") after weaning. In general, we found that bHRs produced more 50 kHz USVs than bLRs. In adulthood, we measured opioid ligand and receptor expression in brain regions implicated in USV production and affect regulation by mRNA in situ hybridization. Within multiple affective brain regions, bHRs had greater expression of the mu opioid receptor than bLRs. FGF2 increased mu opioid expression in bLRs. The bLRs had more kappa and less delta receptor expression than bHRs, and FGF2 increased prodynorphin in bLRs. Our results provide support for further investigations into the role of growth factors and endogenous opioids in the treatment of disorders characterized by altered affect, such as anxiety and depression.

Neural cell adhesion molecule peptide mimetics modulate emotionality: pharmacokinetic and behavioral studies in rats and non-human primates

Recent evidence highlights the fibroblast growth factor (FGF) family in emotion modulation. Although ligands that activate FGF receptors have antidepressant and anxiolytic effects in animal models, FGF ligands have a broad range of actions both in the brain and the periphery. Therefore, identifying molecular partners that may function as allosteric modulators could offer new avenues for drug development. Since neural cell adhesion molecule (NCAM) activates FGF receptors, we asked whether peripherally administered NCAM peptide mimetics penetrate the brain and alter the behavior of standardized tests that have predictive validity for drug treatments of anxiety or depression. The NCAM peptide mimetic, plannexin, acutely increased and chronically decreased anxiety, but did not have antidepressant effects in rats. Another NCAM peptide mimetic, FGL_L, had acute anxiogenic effects and chronic antidepressant effects in rats. A related NCAM peptide mimetic, FGL_S, had antidepressant effects without modulating anxiety-like behavior, and these antidepressant effects were blocked by an AMPA receptor antagonist. Cisternal cerebrospinal fluid (CSF) levels of FGLs correlated with blood plasma levels in rats and non-human primates, and CSF-to-blood ratios of FGL_S were comparable in both species. Results indicate that NCAM peptide mimetics penetrate the brain and support the suggestion that FGL_S may be a candidate for further development as a novel treatment for major depressive disorder in humans.

Quantitative validation of immunofluorescence and lectin staining using reduced CLARITY acrylamide formulations

The CLARITY technique enables three-dimensional visualization of fluorescent-labeled biomolecules in clarified intact brain samples, affording a unique view of molecular neuroanatomy and neurocircuitry. It is therefore, essential to find the ideal combination for clearing tissue and detecting the fluorescent-labeled signal. This method requires the formation of a formaldehyde-acrylamide fixative-generated hydrogel mesh through which cellular lipid is removed with sodium dodecyl sulfate. Several laboratories have used differential acrylamide and detergent concentrations to achieve better tissue clearing and antibody penetration, but the potential effects upon fluorescent signal retention is largely unknown. In an effort to optimize CLARITY processing procedures we performed quantitative parvalbumin immunofluorescence and lectin-based vasculature staining using either 4 or 8% sodium dodecyl sulfate detergent in combination with different acrylamide formulas in mouse brain slices. Using both confocal and CLARITY-optimized lightsheet microscope-acquired images, we demonstrate that 2% acrylamide monomer combined with 0.0125% bis-acrylamide and cleared with 4% sodium dodecyl sulfate generally provides the most optimal signal visualization amongst various hydrogel monomer concentrations, lipid removal times, and detergent concentrations.

Mitochondrial Complex I Deficiency in Schizophrenia and Bipolar Disorder and Medication Influence

Subjects with schizophrenia (SZ) and bipolar disorder (BD) show decreased protein and transcript levels for mitochondrial complex I. In vitro results suggest antipsychotic and antidepressant drugs may be responsible. We measured complex I activity in BD, SZ, and controls and presence of antipsychotic and antidepressant medications, mitochondrial DNA (mtDNA) copy number, and the mtDNA "common deletion" in the brain. Complex I activity in the prefrontal cortex was decreased by 45% in SZ compared to controls (p = 0.02), while no significant difference was found in BD. Complex I activity was significantly decreased (p = 0.01) in pooled cases (SZ and BD) that had detectable psychotropic medications and drugs compared to pooled cases with no detectable levels. Subjects with age at onset in their teens and psychotropic medications showed decreased (p < 0.05) complex I activity compared to subjects with an adult age at onset. Both SZ and BD groups displayed significant increases (p < 0.05) in mtDNA copy number compared to controls; however, common deletion burden was not altered. Complex I deficiency is found in SZ brain tissue, and psychotropic medications may play a role in mitochondrial dysfunction. Studies of medication-free first-episode psychosis patients are needed to elucidate whether mitochondrial pathophysiology occurs independent of medication effects.

Stress amplifies sex differences in primate prefrontal profiles of gene expression

BACKGROUND

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.

METHODS

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.

RESULTS

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.

CONCLUSIONS

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.

Selectively Bred Rats Provide a Unique Model of Vulnerability to PTSD-Like Behavior and Respond Differentially to FGF2 Augmentation Early in Life

Individuals respond differently to traumatic experiences, including their propensity to develop posttraumatic stress disorder (PTSD). Understanding individual differences in PTSD vulnerability will allow the development of improved prevention and treatment options. Here we characterized fear conditioning and extinction in rats selectively bred for differences in their locomotor response to a novel environment. Selectively bred high-responder (bHR) and low-responder (bLR) male rats are known to differ in their emotional reactivity on a range of measures of spontaneous anxiety- and depressive-like behaviors. We demonstrate that bHRs have facilitated extinction learning and retention compared with outbred Sprague Dawley rats, whereas bLRs show reduced extinction learning and retention. This indicates that bLRs are more vulnerable to PTSD-like behavior. Fibroblast growth factor 2 (FGF2) has previously been implicated in the development of these behavioral phenotypes and facilitates extinction learning in outbred animals, therefore we examined the effects of early-life FGF2 on bHR and bLR behavior. FGF2 administered on the day after birth facilitated extinction learning and retention in bHRs, but not in bLRs or control rats, during adulthood. This indicates that, under the current fear conditioning paradigm, early-life FGF2 has protective effects only in resilient animals. This stands in contrast to FGF2's ability to protect vulnerable animals in milder tests of anxiety. These results provide a unique animal model of individual differences in PTSD-like behavior, allowing the study of genetic, developmental, and environmental factors in its expression.

Adolescent cocaine exposure enhances goal-tracking behavior and impairs hippocampal cell genesis selectively in adult bred low-responder rats

RATIONALE

Environmental challenges during adolescence, such as drug exposure, can cause enduring behavioral and molecular changes that contribute to life-long maladaptive behaviors, including addiction. Selectively bred high-responder (bHR) and low-responder (bLR) rats represent a unique model for assessing the long-term impact of adolescent environmental manipulations, as they inherently differ on a number of addiction-related traits. bHR rats are considered "addiction-prone," whereas bLR rats are "addiction-resilient," at least under baseline conditions. Moreover, relative to bLRs, bHR rats are more likely to attribute incentive motivational value to reward cues, or to "sign-track."

OBJECTIVES

We utilized bHR and bLR rats to determine whether adolescent cocaine exposure can alter their inborn behavioral and neurobiological profiles, with a specific focus on Pavlovian conditioned approach behavior (i.e., sign- vs. goal-tracking) and hippocampal neurogenesis.

METHODS

bHR and bLR rats were administered cocaine (15 mg/kg) or saline for 7 days during adolescence (postnatal day, PND 33-39) and subsequently tested for Pavlovian conditioned approach behavior in adulthood (PND 62-75), wherein an illuminated lever (conditioned stimulus) was followed by the response-independent delivery of a food pellet (unconditioned stimulus). Behaviors directed toward the lever and the food cup were recorded as sign- and goal-tracking, respectively. Hippocampal cell genesis was evaluated on PND 77 by immunohistochemistry.

RESULTS

Adolescent cocaine exposure impaired hippocampal cell genesis (proliferation and survival) and enhanced the inherent propensity to goal-track in adult bLR, but not bHR, rats.

CONCLUSIONS

Adolescent cocaine exposure elicits long-lasting changes in stimulus-reward learning and enduring deficits in hippocampal neurogenesis selectively in adult bLR rats.

Utilizing a unique animal model to better understand human temperament

Individual differences in temperament are associated with psychopathology in humans. Moreover, the relationship between temperament and anxiety-, depression-, PTSD- and addiction-related behaviors can be modeled in animals. This review will highlight these relationships with a focus on individual differences in the response to stressors, fear conditioning and drugs of abuse using animals that differ in their response to a novel environment. We will discuss behavioral and neurobiological commonalities amongst these behaviors with a focus on the hippocampus and, in particular, growth factors as promising novel targets for therapeutic intervention.

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.

The microRNA network is altered in anterior cingulate cortex of patients with unipolar and bipolar depression

MicroRNAs (miRNAs) are small, non-coding RNAs acting as post-transcriptional regulators of gene expression. Though implicated in multiple CNS disorders, miRNAs have not been examined in any psychiatric disease state in anterior cingulate cortex (AnCg), a brain region centrally involved in regulating mood. We performed qPCR analyses of 29 miRNAs previously implicated in psychiatric illness (major depressive disorder (MDD), bipolar disorder (BP) and/or schizophrenia (SZ)) in AnCg of patients with MDD and BP versus controls. miR-132, miR-133a and miR-212 were initially identified as differentially expressed in BP, miR-184 in MDD and miR-34a in both MDD and BP (although none survived multiple correction testing and must be considered preliminary). In silico target prediction algorithms identified putative targets of differentially expressed miRNAs. Nuclear Co-Activator 1 (NCOA1), Nuclear Co-Repressor 2 (NCOR2) and Phosphodiesterase 4B (PDE4B) were selected based upon predicted targeting by miR-34a (with NCOR2 and PDE4B both targeted by miR-184) and published relevance to psychiatric illness. Luciferase assays identified PDE4B as a target of miR-34a and miR-184, while NCOA1 and NCOR2 were targeted by miR-34a and 184, respectively. qPCR analyses were performed to determine whether changes in miRNA levels correlated with mRNA levels of validated targets. NCOA1 showed an inverse correlation with miR-34a in BP, while NCOR2 demonstrated a positive correlation. In sum, this is the first study to demonstrate miRNA changes in AnCg in psychiatric illness and validate miR-34a as differentially expressed in CNS in MDD. These findings support a mechanistic role for miRNAs in the regulation of stress-responsive genes disrupted in psychiatric illness.

Quantitative Trait Locus and Brain Expression of HLA-DPA1 Offers Evidence of Shared Immune Alterations in Psychiatric Disorders

Genome-wide association studies of schizophrenia encompassing the major histocompatibility locus (MHC) were highly significant following genome-wide correction. This broad region implicates many genes including the MHC complex class II. Within this interval we examined the expression of two MHC II genes (HLA-DPA1 and HLA-DRB1) in brain from individual subjects with schizophrenia (SZ), bipolar disorder (BD), major depressive disorder (MDD), and controls by differential gene expression methods. A third MHC II mRNA, CD74, was studied outside of the MHC II locus, as it interacts within the same immune complex. Exon microarrays were performed in anterior cingulate cortex (ACC) in BD compared to controls, and both HLA-DPA1 and CD74 were decreased in expression in BD. The expression of HLA-DPA1 and CD74 were both reduced in hippocampus, amygdala, and dorsolateral prefrontal cortex regions in SZ and BD compared to controls by specific qPCR assay. We found several novel HLA-DPA1 mRNA variants spanning HLA-DPA1 exons 2-3-4 as suggested by exon microarrays. The intronic rs9277341 SNP was a significant cis expression quantitative trait locus (eQTL) that was associated with the total expression of HLA-DPA1 in five brain regions. A biomarker study of MHC II mRNAs was conducted in SZ, BD, MDD, and control lymphoblastic cell lines (LCL) by qPCR assay of 87 subjects. There was significantly decreased expression of HLA-DPA1 and CD74 in BD, and trends for reductions in SZ in LCLs. The discovery of multiple splicing variants in brain for HLA-DPA1 is important as the HLA-DPA1 gene is highly conserved, there are no reported splicing variants, and the functions in brain are unknown. Future work on the function and localization of MHC Class II proteins in brain will help to understand the role of alterations in neuropsychiatric disorders. The HLA-DPA1 eQTL is located within a large linkage disequilibrium block that has an irrefutable association with schizophrenia. Future tests in a larger cohort are needed to determine the significance of this eQTL association with schizophrenia. Our findings support the long-held hypothesis that alterations in immune function are associated with the pathophysiology of psychiatric disorders.

Nucleus accumbens cocaine-amphetamine regulated transcript mediates food intake during novelty conflict

Obesity is a persistent and pervasive problem, particularly in industrialized nations. It has come to be appreciated that the metabolic health of an individual can influence brain function and subsequent behavioral patterns. To examine the relationship between metabolic phenotype and central systems that regulate behavior, we tested rats with divergent metabolic phenotypes (Low Capacity Runner: LCR vs. High Capacity Runner: HCR) for behavioral responses to the conflict between hunger and environmental novelty using the novelty suppressed feeding (NSF) paradigm. Additionally, we measured expression of mRNA, for peptides involved in energy management, in response to fasting. Following a 24-h fast, LCR rats showed lower latencies to begin eating in a novel environment compared to HCR rats. A 48-h fast equilibrated the latency to begin eating in the novel environment. A 24-h fast differentially affected expression of cocaine-amphetamine regulated transcript (CART) mRNA in the nucleus accumbens (NAc), where 24-h of fasting reduced CART mRNA in LCR rats. Bilateral microinjections of CART 55-102 peptide into the NAc increased the latency to begin eating in the NSF paradigm following a 24-h fast in LCR rats. These results indicate that metabolic phenotype influences how animals cope with the conflict between hunger and novelty, and that these differences are at least partially mediated by CART signaling in the NAc. For individuals with poor metabolic health who have to navigate food-rich and stressful environments, changes in central systems that mediate conflicting drives may feed into the rates of obesity and exacerbate the difficulty individuals have in maintaining weight loss.

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.

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.

Individual Differences in Cue-Induced Motivation and Striatal Systems in Rats Susceptible to Diet-Induced Obesity

Pavlovian cues associated with junk-foods (caloric, highly sweet, and/or fatty foods), like the smell of brownies, can elicit craving to eat and increase the amount of food consumed. People who are more susceptible to these motivational effects of food cues may have a higher risk for becoming obese. Further, overconsumption of junk-foods leading to the development of obesity may itself heighten attraction to food cues. Here, we used a model of individual susceptibility to junk-foods diet-induced obesity to determine whether there are pre-existing and/or diet-induced increases in attraction to and motivation for sucrose-paired cues (ie, incentive salience or 'wanting'). We also assessed diet- vs obesity-associated alterations in mesolimbic function and receptor expression. We found that rats susceptible to diet-induced obesity displayed heightened conditioned approach prior to the development of obesity. In addition, after junk-food diet exposure, those rats that developed obesity also showed increased willingness to gain access to a sucrose cue. Heightened 'wanting' was not due to individual differences in the hedonic impact ('liking') of sucrose. Neurobiologically, Mu opioid receptor mRNA expression was lower in striatal 'hot-spots' that generate eating or hedonic impact only in those rats that became obese. In contrast, prolonged exposure to junk-food resulted in cross-sensitization to amphetamine-induced locomotion and downregulation of striatal D2R mRNA regardless of the development of obesity. Together these data shed light on individual differences in behavioral and neurobiological consequences of exposure to junk-food diets and the potential contribution of incentive sensitization in susceptible individuals to greater food cue-triggered motivation.

Mitochondrial mutations in subjects with psychiatric disorders

A considerable body of evidence supports the role of mitochondrial dysfunction in psychiatric disorders and mitochondrial DNA (mtDNA) mutations are known to alter brain energy metabolism, neurotransmission, and cause neurodegenerative disorders. Genetic studies focusing on common nuclear genome variants associated with these disorders have produced genome wide significant results but those studies have not directly studied mtDNA variants. The purpose of this study is to investigate, using next generation sequencing, the involvement of mtDNA variation in bipolar disorder, schizophrenia, major depressive disorder, and methamphetamine use. MtDNA extracted from multiple brain regions and blood were sequenced (121 mtDNA samples with an average of 8,800x coverage) and compared to an electronic database containing 26,850 mtDNA genomes. We confirmed novel and rare variants, and confirmed next generation sequencing error hotspots by traditional sequencing and genotyping methods. We observed a significant increase of non-synonymous mutations found in individuals with schizophrenia. Novel and rare non-synonymous mutations were found in psychiatric cases in mtDNA genes: ND6, ATP6, CYTB, and ND2. We also observed mtDNA heteroplasmy in brain at a locus previously associated with schizophrenia (T16519C). Large differences in heteroplasmy levels across brain regions within subjects suggest that somatic mutations accumulate differentially in brain regions. Finally, multiplasmy, a heteroplasmic measure of repeat length, was observed in brain from selective cases at a higher frequency than controls. These results offer support for increased rates of mtDNA substitutions in schizophrenia shown in our prior results. The variable levels of heteroplasmic/multiplasmic somatic mutations that occur in brain may be indicators of genetic instability in mtDNA.

Circadian dysregulation of clock genes: clues to rapid treatments in major depressive disorder

Conventional antidepressants require 2-8 weeks for a full clinical response. In contrast, two rapidly acting antidepressant interventions, low-dose ketamine and sleep deprivation (SD) therapy, act within hours to robustly decrease depressive symptoms in a subgroup of major depressive disorder (MDD) patients. Evidence that MDD may be a circadian-related illness is based, in part, on a large set of clinical data showing that diurnal rhythmicity (sleep, temperature, mood and hormone secretion) is altered during depressive episodes. In a microarray study, we observed widespread changes in cyclic gene expression in six regions of postmortem brain tissue of depressed patients matched with controls for time-of-death (TOD). We screened 12 000 transcripts and observed that the core clock genes, essential for controlling virtually all rhythms in the body, showed robust 24-h sinusoidal expression patterns in six brain regions in control subjects. In MDD patients matched for TOD with controls, the expression patterns of the clock genes in brain were significantly dysregulated. Some of the most robust changes were seen in anterior cingulate (ACC). These findings suggest that in addition to structural abnormalities, lesion studies, and the large body of functional brain imaging studies reporting increased activation in the ACC of depressed patients who respond to a wide range of therapies, there may be a circadian dysregulation in clock gene expression in a subgroup of MDDs. Here, we review human, animal and neuronal cell culture data suggesting that both low-dose ketamine and SD can modulate circadian rhythms. We hypothesize that the rapid antidepressant actions of ketamine and SD may act, in part, to reset abnormal clock genes in MDD to restore and stabilize circadian rhythmicity. Conversely, clinical relapse may reflect a desynchronization of the clock, indicative of a reactivation of abnormal clock gene function. Future work could involve identifying specific small molecules capable of resetting and stabilizing clock genes to evaluate if they can rapidly relieve symptoms and sustain improvement.

Altered choroid plexus gene expression in major depressive disorder

Given the emergent interest in biomarkers for mood disorders, we assessed gene expression in the choroid plexus (CP), the region that produces cerebrospinal fluid (CSF), in individuals with major depressive disorder (MDD). Genes that are expressed in the CP can be secreted into the CSF and may be potential biomarker candidates. Given that we have previously shown that fibroblast growth factor family members are differentially expressed in post-mortem brain of subjects with MDD and the CP is a known source of growth factors in the brain, we posed the question whether growth factor dysregulation would be found in the CP of subjects with MDD. We performed laser capture microscopy of the CP at the level of the hippocampus in subjects with MDD and psychiatrically normal controls. We then extracted, amplified, labeled, and hybridized the cRNA to Illumina BeadChips to assess gene expression. In controls, the most highly abundant known transcript was transthyretin. Moreover, half of the 14 most highly expressed transcripts in controls encode ribosomal proteins. Using BeadStudio software, we identified 169 transcripts differentially expressed (p < 0.05) between control and MDD samples. Using pathway analysis we noted that the top network altered in subjects with MDD included multiple members of the transforming growth factor-beta (TGFβ) pathway. Quantitative real-time PCR (qRT-PCR) confirmed downregulation of several transcripts that interact with the extracellular matrix in subjects with MDD. These results suggest that there may be an altered cytoskeleton in the CP in MDD subjects that may lead to a disrupted blood-CSF-brain barrier.

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.

Antecedents and consequences of drug abuse in rats selectively bred for high and low response to novelty

Human genetic and epidemiological studies provide evidence that only a subset of individuals who experiment with potentially addictive drugs become addicts. What renders some individuals susceptible to addiction remains to be determined, but most would agree that there is no single trait underlying the disorder. However, there is evidence in humans that addiction liability has a genetic component, and that certain personality characteristics related to temperament (e.g. the sensation-seeking trait) are associated with individual differences in addiction liability. Consequently, we have used a selective breeding strategy based on locomotor response to a novel environment to generate two lines of rats with distinct behavioral characteristics. We have found that the resulting phenotypes differ on a number of neurobehavioral dimensions relevant to addiction. Relative to bred low-responder (bLR) rats, bred high-responder (bHR) rats exhibit increased exploratory behavior, are more impulsive, more aggressive, seek stimuli associated with rewards, and show a greater tendency to relapse. We therefore utilize this unique animal model to parse the genetic, neural and environmental factors that contribute to addiction liability. Our work shows that the glucocorticoid receptor (GR), dopaminergic molecules, and members of the fibroblast growth factor family are among the neurotransmitters and neuromodulators that play a role in both the initial susceptibility to addiction as well as the altered neural responses that follow chronic drug exposure. Moreover, our findings suggest that the hippocampus plays a major role in mediating vulnerability to addiction. It is hoped that this work will emphasize the importance of personalized treatment strategies and identify novel therapeutic targets for humans suffering from addictive disorders. This article is part of a Special Issue entitled 'NIDA 40th Anniversary Issue'.

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.

Evidence of allelic imbalance in the schizophrenia susceptibility gene ZNF804A in human dorsolateral prefrontal cortex

The rs1344706, an intronic SNP within the zinc-finger protein 804A gene (ZNF804A), was identified as one of the most compelling risk SNPs for schizophrenia (SZ) and bipolar disorder (BD). It is however not clear by which molecular mechanisms ZNF804A increases disease risk. We evaluated the role of ZNF804A in SZ and BD by genotyping the originally associated rs1344706 SNP and an exonic SNP (rs12476147) located in exon four of ZNF804A in a sample of 422 SZ, 382 BD, and 507 controls from the isolated population of the Costa Rica Central Valley. We also investigated the rs1344706 SNP for allelic specific expression (ASE) imbalance in the dorsolateral prefrontal cortex (DLPFC) of 46 heterozygous postmortem brains. While no significant association between rs1344706 and SZ or BD was observed in the Costa Rica sample, we observed an increased risk of SZ for the minor allele (A) of the exonic rs12476147 SNP (p=0.026). Our ASE assay detected a significant over-expression of the rs12476147 A allele in DLPFC of rs1344706 heterozygous subjects. Interestingly, cDNA allele ratios were significantly different according to the intronic rs1344706 genotypes (p-value=0.03), with the rs1344706 A allele associated with increased ZNF804A rs12476147 A allele expression (average 1.06, p-value=0.02, for heterozygous subjects vs. genomic DNA). In conclusion, we have demonstrated a significant association of rs12476147 with SZ, and using a powerful within-subject design, an allelic expression imbalance of ZNF804A exonic SNP rs12476147 in the DLPFC. Although this data does not preclude the possibility of other functional variants in ZNF804A, it provides evidence that the rs1344706 SZ risk allele is the cis-regulatory variant directly responsible for this allelic expression imbalance in adult cortex.

G protein-linked signaling pathways in bipolar and major depressive disorders

The G-protein linked signaling system (GPLS) comprises a large number of G-proteins, G protein-coupled receptors (GPCRs), GPCR ligands, and downstream effector molecules. G-proteins interact with both GPCRs and downstream effectors such as cyclic adenosine monophosphate (cAMP), phosphatidylinositols, and ion channels. The GPLS is implicated in the pathophysiology and pharmacology of both major depressive disorder (MDD) and bipolar disorder (BPD). This study evaluated whether GPLS is altered at the transcript level. The gene expression in the dorsolateral prefrontal (DLPFC) and anterior cingulate (ACC) were compared from MDD, BPD, and control subjects using Affymetrix Gene Chips and real time quantitative PCR. High quality brain tissue was used in the study to control for confounding effects of agonal events, tissue pH, RNA integrity, gender, and age. GPLS signaling transcripts were altered especially in the ACC of BPD and MDD subjects. Transcript levels of molecules which repress cAMP activity were increased in BPD and decreased in MDD. Two orphan GPCRs, GPRC5B and GPR37, showed significantly decreased expression levels in MDD, and significantly increased expression levels in BPD. Our results suggest opposite changes in BPD and MDD in the GPLS, “activated” cAMP signaling activity in BPD and “blunted” cAMP signaling activity in MDD. GPRC5B and GPR37 both appear to have behavioral effects, and are also candidate genes for neurodegenerative disorders. In the context of the opposite changes observed in BPD and MDD, these GPCRs warrant further study of their brain effects.

Interaction between cholecystokinin and the fibroblast growth factor system in the ventral tegmental area of selectively bred high- and low-responder rats

Individual differences in the locomotor response to novelty have been linked to basal differences in dopaminergic neurotransmission. Mesolimbic dopaminergic outputs are regulated by cholecystokinin (CCK), a neuropeptide implicated in anxiety. In turn, CCK expression is regulated by fibroblast growth factor-2 (FGF2), which has recently been identified as an endogenous regulator of anxiety. FGF2 binds to the high-affinity fibroblast growth factor receptor-1 (FGF-R1) to regulate the development and maintenance of dopamine neurons in the ventral tegmental area (VTA). However, the relationship between the FGF and CCK systems in the VTA is not well understood. Therefore, we utilized the selectively-bred low-responder (bLR; high-anxiety) and high-responder (bHR; low-anxiety) rats to examine the effects of repeated (21-day) FGF2 treatment on CCK and FGF-R1 mRNA in the rostral VTA (VTAr). In vehicle-treated controls, both CCK and FGF-R1 mRNA levels were increased in the VTAr of bLR rats relative to bHR rats. Following FGF2 treatment, however, bHR-bLR differences in CCK and FGF-R1 mRNA expression were eliminated, due to decreased CCK mRNA levels in the VTAr of bLR rats and increased FGF-R1 expression in bHR rats. Differences after FGF2 treatment may denote distinct interactions between the CCK and FGF systems in the VTAr of bHR vs. bLR rats. Indeed, significant correlations between CCK and FGF-R1 mRNA expression were found in bHR, but not bLR rats. Colocalization studies suggest that CCK and FGF-R1 are coexpressed in some VTAr neurons. Taken together, our findings suggest that the FGF system is poised to modulate both CCK and FGF-R1 expression in the VTAr, which may be associated with individual differences in mesolimbic pathways associated with anxiety-like behavior.

Long-term effects of cocaine experience on neuroplasticity in the nucleus accumbens core of addiction-prone rats

Repeated exposure to drugs of abuse is associated with structural plasticity in brain reward pathways. Rats selectively bred for locomotor response to novelty differ on a number of neurobehavioral dimensions relevant to addiction. This unique genetic animal model was used here to examine both pre-existing differences and long-term consequences of repeated cocaine treatment on structural plasticity. Selectively bred high-responder (bHR) and low-responder (bLR) rats received repeated saline or cocaine injections for 9 consecutive days. Escalating doses of cocaine (7.5, 15 and 30 mg/kg) were administered on the first (day 1) and last (day 9) days of treatment and a single injection of the intermediate dose (15 mg/kg) was given on days 2-8. Motor activity in response to escalating doses of cocaine was compared on the first and last days of treatment to assess the acute and sensitized response to the drug. Following prolonged cocaine abstinence (28 days), spine density was examined on terminal dendrites of medium spiny neurons in the nucleus accumbens core. Relative to bLRs, bHRs exhibited increased psychomotor activation in response to both the acute and repeated effects of cocaine. There were no differences in spine density between bHR and bLR rats under basal conditions or following repeated saline treatment. However, spine density differed markedly between these two lines following prolonged cocaine abstinence. All spine types were decreased in cocaine-treated bHRs, while only mushroom spines were decreased in bLRs that received cocaine. Changes in spine density occurred specifically near the branch point of terminal dendrites. These findings indicate that structural plasticity associated with prolonged cocaine abstinence varies markedly in two selected strains of rats that vary on numerous traits relevant to addiction. Thus, genetic factors that contribute to individual variation in the behavioral response to cocaine also influence cocaine-induced structural plasticity.

Glutamate transporters: a key piece in the glutamate puzzle of major depressive disorder

Glutamatergic therapies are emerging as the new path for the treatment of Major Depression Disorder. Recent reports reviewing the use of glutamate activity modulators in the treatment of resistant depression advocate the importance of understanding the alterations of the diverse components of this complex system in mood disorders. In this postmortem study we used in situ hybridization and microarray analysis to evaluate the gene expression of the membrane transporters SLC1A2 and SLCA3 and the vesicular transporter SLCA17A7 in the hippocampus of Major Depressive Disorder (MDD) and Bipolar Disorder (BPD) subjects. Samples from 8 controls, 11 MDD and 6 BPD subjects were processed for in situ hybridization using cRNA probes for SLC1A2, SLC1A3 and SLC17A7. Laser capture microdissection was used to collect tissue from adjacent sections for microarray analysis. The results showed that the expression of the membrane transporters SLC1A2 and SLC1A3 was diminished in the MDD group compared to controls. The expression of the vesicular glutamate transporter SLC17A7 on the other hand was increased in MDD subjects. As for the BPD group, all three transporters showed trends similar to those observed in MDD, but the changes observed did not reach significance. We hypothesize that the decreased expression of the membrane glutamate transporters and the increased expression of the vesicular transporter in the hippocampus would affect the balance of the glutamatergic circuitry of the hippocampus, and that this effect may be a major contributor to depressive symptoms.