Mitochondria provide most of the energy for brain cells by the process of oxidative phosphorylation. Mitochondrial abnormalities and deficiencies in oxidative phosphorylation have been reported in individuals with schizophrenia (SZ), bipolar disorder (BD), and major depressive disorder (MDD) in transcriptomic, proteomic, and metabolomic studies. Several mutations in mitochondrial DNA (mtDNA) sequence have been reported in SZ and BD patients.

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.

Whole-genome and exome sequencing have already proven to be essential and powerful methods to identify genes responsible for simple Mendelian inherited disorders. These methods can be applied to complex disorders as well, and have been adopted as one of the current mainstream approaches in population genetics. These achievements have been made possible by next generation sequencing (NGS) technologies, which require substantial bioinformatics resources to analyze the dense and complex sequence data. The huge analytical burden of data from genome sequencing might be seen as a bottleneck slowing the publication of NGS papers at this time, especially in psychiatric genetics. We review the existing methods for processing NGS data, to place into context the rationale for the design of a computational resource. We describe our method, the Graphical Pipeline for Computational Genomics (GPCG), to perform the computational steps required to analyze NGS data. The GPCG implements flexible workflows for basic sequence alignment, sequence data quality control, single nucleotide polymorphism analysis, copy number variant identification, annotation, and visualization of results. These workflows cover all the analytical steps required for NGS data, from processing the raw reads to variant calling and annotation. The current version of the pipeline is freely available at http://pipeline.loni.ucla.edu. These applications of NGS analysis may gain clinical utility in the near future (e.g., identifying miRNA signatures in diseases) when the bioinformatics approach is made feasible. Taken together, the annotation tools and strategies that have been developed to retrieve information and test hypotheses about the functional role of variants present in the human genome will help to pinpoint the genetic risk factors for psychiatric disorders.

We investigated the genome-wide distribution of CNVs in the Alzheimer's disease (AD) Neuroimaging Initiative (ADNI) sample (146 with AD, 313 with Mild Cognitive Impairment (MCI), and 181 controls). Comparison of single CNVs between cases (MCI and AD) and controls shows overrepresentation of large heterozygous deletions in cases (p-value<0.0001). The analysis of CNV-Regions identifies 44 copy number variable loci of heterozygous deletions, with more CNV-Regions among affected than controls (p=0.005). Seven of the 44 CNV-Regions are nominally significant for association with cognitive impairment. We validated and confirmed our main findings with genome re-sequencing of selected patients and controls. The functional pathway analysis of the genes putatively affected by deletions of CNV-Regions reveals enrichment of genes implicated in axonal guidance, cell-cell adhesion, neuronal morphogenesis and differentiation. Our findings support the role of CNVs in AD, and suggest an association between large deletions and the development of cognitive impairment.

Mitochondrial deficiencies with unknown causes have been observed in schizophrenia (SZ) and bipolar disorder (BD) in imaging and postmortem studies. Polymorphisms and somatic mutations in mitochondrial DNA (mtDNA) were investigated as potential causes with next generation sequencing of mtDNA (mtDNA-Seq) and genotyping arrays in subjects with SZ, BD, major depressive disorder (MDD), and controls. The common deletion of 4,977 bp in mtDNA was compared between SZ and controls in 11 different vulnerable brain regions and in blood samples, and in dorsolateral prefrontal cortex (DLPFC) of BD, SZ, and controls. In a separate analysis, association of mitochondria SNPs (mtSNPs) with SZ and BD in European ancestry individuals (n = 6,040) was tested using Genetic Association Information Network (GAIN) and Wellcome Trust Case Control Consortium 2 (WTCCC2) datasets. The common deletion levels were highly variable across brain regions, with a 40-fold increase in some regions (nucleus accumbens, caudate nucleus and amygdala), increased with age, and showed little change in blood samples from the same subjects. The common deletion levels were increased in the DLPFC for BD compared to controls, but not in SZ. Full mtDNA genome resequencing of 23 subjects, showed seven novel homoplasmic mutations, five were novel synonymous coding mutations. By logistic regression analysis there were no significant mtSNPs associated with BD or SZ after genome wide correction. However, nominal association of mtSNPs (p < 0.05) to SZ and BD were found in the hypervariable region of mtDNA to T195C and T16519C. The results confirm prior reports that certain brain regions accumulate somatic mutations at higher levels than blood. The study in mtDNA of common polymorphisms, somatic mutations, and rare mutations in larger populations may lead to a better understanding of the pathophysiology of psychiatric disorders.

Suicidal behaviors are frequent in mood disorders patients but only a subset of them ever complete suicide. Understanding predisposing factors for suicidal behaviors in high risk populations is of major importance for the prevention and treatment of suicidal behaviors. The objective of this project was to investigate gene expression changes associated with suicide in brains of mood disorder patients by microarrays (Affymetrix HG-U133 Plus2.0) in the dorsolateral prefrontal cortex (DLPFC: 6 Non-suicides, 15 suicides), the anterior cingulate cortex (ACC: 6NS, 9S) and the nucleus accumbens (NAcc: 8NS, 13S). ANCOVA was used to control for age, gender, pH and RNA degradation, with P ≤ 0.01 and fold change ± 1.25 as criteria for significance. Pathway analysis revealed serotonergic signaling alterations in the DLPFC and glucocorticoid signaling alterations in the ACC and NAcc. The gene with the lowest p-value in the DLPFC was the 5-HT2A gene, previously associated both with suicide and mood disorders. In the ACC 6 metallothionein genes were down-regulated in suicide (MT1E, MT1F, MT1G, MT1H, MT1X, MT2A) and three were down-regulated in the NAcc (MT1F, MT1G, MT1H). Differential expression of selected genes was confirmed by qPCR, we confirmed the 5-HT2A alterations and the global down-regulation of members of the metallothionein subfamilies MT 1 and 2 in suicide completers. MTs 1 and 2 are neuro-protective following stress and glucocorticoid stimulations, suggesting that in suicide victims neuroprotective response to stress and cortisol may be diminished. Our results thus suggest that suicide-specific expression changes in mood disorders involve both glucocorticoids regulated metallothioneins and serotonergic signaling in different regions of the brain.

Brain-derived neurotrophic factor (BDNF), a neurotrophin highly expressed in the hippocampus, plays crucial roles in cognition, neuroplasticity, synaptic function, and dendritic remodeling. The common human Val66Met polymorphism of BDNF has been implicated in the pathophysiology of neuropsychiatric and neurodegenerative disorders, and in the outcome of pro-adaptive and therapeutic treatments. Altered gene-expression profile has been previously shown in BDNF Val66Met knock-in mice, which recapitulate the phenotypic hallmarks of individuals carrying the BDNF Met allele. The aim of this study was to investigate the impact of the BDNF Val66Met polymorphism in the knock-in mouse model on two hippocampal epigenetic marks for transcriptional repression and activation, respectively: trimethylation of lysine 27 on histone H3 (H3K27me3) and acetylation of histone H3 (AcH3), using a genome-wide approach. Chromatin immunoprecipitation followed by deep sequencing of immunoprecipitated DNA (ChIP-Seq) was carried out with specific antibodies for H3K27me3 and AcH3. Our results revealed broad alteration of H3K27me3 and AcH3 marks association profiles in BDNFMet/Met , compared to BDNFVal/Val mice. Bioinformatics analysis showed changes in several biological functions and related pathways, affected by the presence of the polymorphism. In particular, a number of networks of functional interaction contained BDNF as central node. Quantitative PCR analysis confirmed epigenetically related significant changes in the expression of five genes: Dvl1, Nos3, Reln, Lypd6, and Sh3gl2. The first three are involved in dendrite and spine remodeling, morphological features altered in BDNFMet/Met mice. This work in homozygous knock-in mice shows that the human BDNF Val66Met polymorphism induces an array of histone H3 epigenetic modifications, in turn altering the expression of select genes crucial for structural and functional neuronal features.

A human variant in the BDNF gene (Val66Met; rs6265) is associated with impaired fear extinction. Using super-resolution imaging, we demonstrate that the BDNF Met prodomain disassembles dendritic spines and eliminates synapses in hippocampal neurons. In vivo, ventral CA1 (vCA1) hippocampal neurons undergo similar morphological changes dependent on their transient co-expression of a SorCS2/p75NTR receptor complex during peri-adolescence. BDNF Met prodomain infusion into the vCA1 during this developmental time frame reduces dendritic spine density and prelimbic (PL) projections, impairing cued fear extinction. Adolescent BdnfMet/Met mice display similar spine and PL innervation deficits. Using fiber photometry, we found that, in wild-type mice, vCA1 neurons projecting to the PL encode extinction by enhancing neural activity in threat anticipation and rapidly subsiding their response. This adaptation is absent in BDNFMet/Met mice. We conclude that the BDNF Met prodomain renders vCA1-PL projection neurons underdeveloped, preventing their capacity for subsequent circuit modulation necessary for fear extinction. VIDEO ABSTRACT.

Bipolar disorder is a heritable mental illness with complex etiology. We performed a genome-wide association study of 41,917 bipolar disorder cases and 371,549 controls of European ancestry, which identified 64 associated genomic loci. Bipolar disorder risk alleles were enriched in genes in synaptic signaling pathways and brain-expressed genes, particularly those with high specificity of expression in neurons of the prefrontal cortex and hippocampus. Significant signal enrichment was found in genes encoding targets of antipsychotics, calcium channel blockers, antiepileptics and anesthetics. Integrating expression quantitative trait locus data implicated 15 genes robustly linked to bipolar disorder via gene expression, encoding druggable targets such as HTR6, MCHR1, DCLK3 and FURIN. Analyses of bipolar disorder subtypes indicated high but imperfect genetic correlation between bipolar disorder type I and II and identified additional associated loci. Together, these results advance our understanding of the biological etiology of bipolar disorder, identify novel therapeutic leads and prioritize genes for functional follow-up studies.

Hedgehog signaling is essential for bone formation, including functioning as a means for the growth plate to drive skeletal mineralization. However, the mechanisms regulating hedgehog signaling specifically in bone-forming osteoblasts are largely unknown. Here, we identified SLIT and NTRK-like protein-5(Slitrk5), a transmembrane protein with few identified functions, as a negative regulator of hedgehog signaling in osteoblasts. Slitrk5 is selectively expressed in osteoblasts and loss of Slitrk5 enhanced osteoblast differentiation in vitro and in vivo. Loss of SLITRK5 in vitro leads to increased hedgehog signaling and overexpression of SLITRK5 in osteoblasts inhibits the induction of targets downstream of hedgehog signaling. Mechanistically, SLITRK5 binds to hedgehog ligands via its extracellular domain and interacts with PTCH1 via its intracellular domain. SLITRK5 is present in the primary cilium, and loss of SLITRK5 enhances SMO ciliary enrichment upon SHH stimulation. Thus, SLITRK5 is a negative regulator of hedgehog signaling in osteoblasts that may be attractive as a therapeutic target to enhance bone formation.

The G protein-coupled receptor kinase (GRK) family member protein GRK3 has been linked to the pathophysiology of schizophrenia and bipolar disorder. Expression, as well as protein levels, of GRK3 are reduced in post-mortem prefrontal cortex of schizophrenia subjects. Here, we investigate functional behavior and neurotransmission related to immune activation and psychosis using mice lacking functional Grk3 and utilizing a variety of methods, including behavioral, biochemical, electrophysiological, molecular, and imaging methods. Compared to wildtype controls, the Grk3-/- mice show a number of aberrations linked to psychosis, including elevated brain levels of IL-1β, increased turnover of kynurenic acid (KYNA), hyper-responsiveness to D-amphetamine, elevated spontaneous firing of midbrain dopamine neurons, and disruption in prepulse inhibition. Analyzing human genetic data, we observe a link between psychotic features in bipolar disorder, decreased GRK expression, and increased concentration of CSF KYNA. Taken together, our data suggest that Grk3-/- mice show face and construct validity relating to the psychosis phenotype with glial activation and would be suitable for translational studies of novel immunomodulatory agents in psychotic disorders.

The genomic effects of circulating glucocorticoids are particularly relevant in cortico-limbic structures, which express a high concentration of steroid hormone receptors. To date, no studies have investigated genomic differences in hippocampal subregions, namely the dorsal (dHPC) and ventral (vHPC) hippocampus, in preclinical models treated with exogenous glucocorticoids. Chronic oral corticosterone (CORT) in mouse is a pharmacological approach that disrupts the activity of the hypothalamic-pituitary-adrenal axis, increases affective behavior, and induces genomic changes after stress in the HPC of wildtype (WT) mice and mice heterozygous for the gene coding for brain-derived neurotrophic factor Val66Met (hMet), a variant associated with genetic susceptibility to stress. Using RNA-sequencing, we investigated the genomic signatures of oral CORT in the dHPC and vHPC of WT and hMet male and female mice, and examined sex and genotype differences in response to oral CORT. Males under CORT showed lower glycemia and increased anxiety- and depression-like behavior compared to females that showed instead opposite affective behavior in response to CORT. Rank-rank-hypergeometric overlap (RRHO) was used to identify genes from a continuous gradient of significancy that were concordant across groups. RRHO showed that CORT-induced differentially expressed genes (DEGs) in WT mice and hMet mice converged in the dHPC of males and females, while in the vHPC, DEGs converged in males and diverged in females. The vHPC showed a higher number of DEGs compared to the dHPC and exhibited sex differences related to glucocorticoid receptor (GR)-binding genes and epigenetic modifiers. Methyl-DNA-immunoprecipitation in the vHPC revealed differential methylation of the exons 1C and 1F of the GR gene (Nr3c1) in hMet females. Together, we report behavioral and endocrinological sex differences in response to CORT, as well as epigenetic signatures that i) differ in the dHPC and vHPC,ii) are distinct in males and females, and iii) implicate differential methylation of Nr3c1 selectively in hMet females.

The multifactorial etiology of stress-related disorders necessitates a constant interrogation of the molecular convergences in preclinical models of stress that use disparate paradigms as stressors spanning from environmental challenges to genetic predisposition to hormonal signaling. Using RNA-sequencing, we investigated the genomic signatures in the ventral hippocampus common to mouse models of stress. Chronic oral corticosterone (CORT) induced increased anxiety- and depression-like behavior in wild-type male mice and male mice heterozygous for the gene coding for brain-derived neurotrophic factor Val66Met, a variant associated with genetic susceptibility to stress. In a separate set of male mice, chronic social defeat stress (CSDS) led to a susceptible or a resilient population, whose proportion was dependent on housing conditions, namely standard housing or enriched environment. Rank-rank-hypergeometric overlap (RRHO), a threshold-free approach that ranks genes by their p value and effect size direction, was used to identify genes from a continuous gradient of significancy that were concordant across groups. In mice treated with CORT and in standard-housed susceptible mice, differentially expressed genes (DEGs) were concordant for gene networks involved in neurotransmission, cytoskeleton function, and vascularization. Weighted gene co-expression analysis generated 54 gene hub modules and revealed two modules in which both CORT and CSDS-induced enrichment in DEGs, whose function was concordant with the RRHO predictions, and correlated with behavioral resilience or susceptibility. These data showed transcriptional concordance across models in which the stress coping depends upon hormonal, environmental, or genetic factors revealing common genomic drivers that embody the multifaceted nature of stress-related disorders.

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.

Fatty acid amide hydrolase (FAAH) degrades the endocannabinoid anandamide. A polymorphism in FAAH (FAAH C385A) reduces FAAH expression, increases anandamide levels, and increases the risk of obesity. Nevertheless, some studies have found no association between FAAH C385A and obesity. We investigated whether the environmental context governs the impact of FAAH C385A on metabolic outcomes. Using a C385A knock-in mouse model, we found that FAAH A/A mice are more susceptible to glucocorticoid-induced hyperphagia, weight gain, and activation of hypothalamic AMP-activated protein kinase (AMPK). AMPK inhibition occluded the amplified hyperphagic response to glucocorticoids in FAAH A/A mice. FAAH knockdown exclusively in agouti-related protein (AgRP) neurons mimicked the exaggerated feeding response of FAAH A/A mice to glucocorticoids. FAAH A/A mice likewise presented exaggerated orexigenic responses to ghrelin, while FAAH knockdown in AgRP neurons blunted leptin anorectic responses. Together, the FAAH A/A genotype amplifies orexigenic responses and decreases anorexigenic responses, providing a putative mechanism explaining the diverging human findings.

Social hierarchies exert a powerful influence on behavior, but the neurobiological mechanisms that detect and regulate hierarchical interactions are not well understood, especially at the level of neural circuits. Here, we use fiber photometry and chemogenetic tools to record and manipulate the activity of nucleus accumbens-projecting cells in the ventromedial prefrontal cortex (vmPFC-NAcSh) during tube test social competitions. We show that vmPFC-NAcSh projections signal learned hierarchical relationships, and are selectively recruited by subordinate mice when they initiate effortful social dominance behavior during encounters with a dominant competitor from an established hierarchy. After repeated bouts of social defeat stress, this circuit is preferentially activated during social interactions initiated by stress resilient individuals, and plays a necessary role in supporting social approach behavior in subordinated mice. These results define a necessary role for vmPFC-NAcSh cells in the adaptive regulation of social interaction behavior based on prior hierarchical interactions.

Common mitochondrial DNA (mtDNA) deletions are large structural variants in the mitochondrial genome that accumulate in metabolically active tissues with age and have been investigated in various diseases. We applied the Splice-Break2 pipeline (designed for high-throughput quantification of mtDNA deletions) to human RNA-Seq datasets and describe the methodological considerations for evaluating common deletions in bulk, single-cell, and spatial transcriptomics datasets. A robust evaluation of 1570 samples from 14 RNA-Seq studies showed: (i) the abundance of some common deletions detected in PCR-amplified mtDNA correlates with levels observed in RNA-Seq data; (ii) RNA-Seq library preparation method has a strong effect on deletion detection; (iii) deletions had a significant, positive correlation with age in brain and muscle; (iv) deletions were enriched in cortical grey matter, specifically in layers 3 and 5; and (v) brain regions with dopaminergic neurons (i.e., substantia nigra, ventral tegmental area, and caudate nucleus) had remarkable enrichment of common mtDNA deletions.

Abnormalities in the circadian clockwork often characterize patients with major depressive and bipolar disorders. Circadian clock genes are targets of interest in these patients. CRY2 is a circadian gene that participates in regulation of the evening oscillator. This is of interest in mood disorders where a lack of switch from evening to morning oscillators has been postulated.

Multicellular organisms have co-evolved with complex consortia of viruses, bacteria, fungi and parasites, collectively referred to as the microbiota1. In mammals, changes in the composition of the microbiota can influence many physiologic processes (including development, metabolism and immune cell function) and are associated with susceptibility to multiple diseases2. Alterations in the microbiota can also modulate host behaviours-such as social activity, stress, and anxiety-related responses-that are linked to diverse neuropsychiatric disorders3. However, the mechanisms by which the microbiota influence neuronal activity and host behaviour remain poorly defined. Here we show that manipulation of the microbiota in antibiotic-treated or germ-free adult mice results in significant deficits in fear extinction learning. Single-nucleus RNA sequencing of the medial prefrontal cortex of the brain revealed significant alterations in gene expression in excitatory neurons, glia and other cell types. Transcranial two-photon imaging showed that deficits in extinction learning after manipulation of the microbiota in adult mice were associated with defective learning-related remodelling of postsynaptic dendritic spines and reduced activity in cue-encoding neurons in the medial prefrontal cortex. In addition, selective re-establishment of the microbiota revealed a limited neonatal developmental window in which microbiota-derived signals can restore normal extinction learning in adulthood. Finally, unbiased metabolomic analysis identified four metabolites that were significantly downregulated in germ-free mice and have been reported to be related to neuropsychiatric disorders in humans and mouse models, suggesting that microbiota-derived compounds may directly affect brain function and behaviour. Together, these data indicate that fear extinction learning requires microbiota-derived signals both during early postnatal neurodevelopment and in adult mice, with implications for our understanding of how diet, infection, and lifestyle influence brain health and subsequent susceptibility to neuropsychiatric disorders.

Premenstrual dysphoric disorder (PMDD) affects over 5% of women, with symptoms similar to anxiety and major depression, and is associated with differential sensitivity to circulating ovarian hormones. Little is known about the genetic and epigenetic factors that increase the risk to develop PMDD. We report that 17β-estradiol (E2) affects the behavior and the epigenome in a mouse model carrying a single-nucleotide polymorphism of the brain-derived neurotrophic factor gene (BDNF Val66Met), in a way that recapitulates the hallmarks of PMDD. Ovariectomized mice heterozygous for the BDNF Met allele (Het-Met) and their matched wild-type (WT) mice were administered estradiol or vehicle in drinking water for 6 weeks. Using the open field and the splash test, we show that E2 add-back induces anxiety-like and depression-like behavior in Het-Met mice, but not in WT mice. RNA-seq of the ventral hippocampus (vHpc) highlights that E2-dependent gene expression is markedly different between WT mice and Het-Met mice. Through a comparative whole-genome RNA-seq analysis between mouse vHpc and lymphoblastoid cell line cultures from control women and women with PMDD, we discovered common epigenetic biomarkers that transcend species and cell types. Those genes include epigenetic modifiers of the ESC/E(Z) complex, an effector of response to ovarian steroids. Although the BDNF Met genotype intersects the behavioral and transcriptional traits of women with PMDD, we suggest that these similarities speak to the epigenetic factors by which ovarian steroids produce negative behavioral effects.