Schizophrenia (SZ) and bipolar disorder (BP) patients share overlapping neurocognitive deficits of varied magnitude. Neuroimaging in patients and postmortem gene expression analyses suggest that compromised cingulate gyrus GABA-ergic interneurons may contribute to cognitive impairments in SZ and BP. To address this, we used radioactive in situ hybridization to investigate potential gene expression signatures for SZ and BP using interneuron cell-type specific markers including glutamic acid decarboxylase (GAD67), parvalbumin (PV), somatostatin (SST), and vasoactive intestinal peptide (VIP) within specific Brodmann's areas (BA) of the cingulate gyrus. We report reduced GAD67 mRNA in anterior midcingulate cortex (aMCC) of BP subjects within BA24c', the most dysregulated subregion across disorders that also demonstrated reduced PV and VIP mRNA in the SZ group. In the retrosplenial (RSC) and ectosplenial (ESC) cortices, decreases in PV expression were shared by both SZ and BP subjects. Our results show unique and shared transcription signatures of two disorders in specific cingulate gyrus regions and cell types. SZ and BP displayed divergent aMCC gene expression reductions suggesting transcriptional changes are associated with disease-specific gene/subregion signatures, potentially underlying differential subregional dysregulation within areas associated with error detection/action monitoring and the salience network. In RSC/ESC, transcriptional changes are associated with more common expression patterns, possibly related to overlapping effects on visuospatial memory processing and allocation of attentional resources involving the default mode network.

Background: Chronic stress has a profound impact on circadian regulation of physiology. In turn, disruption of circadian rhythms increases the risk of developing both psychiatric and metabolic disorders. To explore the role of chronic stress in modulating the links between neural and metabolic rhythms, we characterized the circadian transcriptional regulation across different brain regions and the liver as well as serum metabolomics in mice exposed to chronic social defeat stress, a validated model for studying depressive-like behaviors. Methods: Male C57BL/6J mice underwent chronic social defeat stress, and subsequent social interaction screening identified distinct behavioral phenotypes associated with stress resilience and susceptibility. Stressed mice and their control littermates were sacrificed every 4 hours over the circadian cycle for comprehensive analyses of the circadian transcriptome in the hypothalamus, hippocampus, prefrontal cortex, and liver together with assessments of the circadian circulatory metabolome. Results: Our data demonstrate that stress adaptation was characterized by reprogramming of the brain as well as the hepatic circadian transcriptome. Stress resiliency was associated with an increase in cyclic transcription in the hypothalamus, hippocampus, and liver. Furthermore, cross-tissue analyses revealed that resilient mice had enhanced transcriptional coordination of circadian pathways between the brain and liver. Conversely, susceptibility to social stress resulted in a loss of cross-tissue coordination. Circadian serum metabolomic profiles corroborated the transcriptome data, highlighting that stress-resilient mice gained circadian rhythmicity of circulating metabolites, including bile acids and sphingomyelins. Conclusions: This study reveals that resilience to stress is characterized by enhanced metabolic rhythms and circadian brain-liver transcriptional coordination.

Transcription factors (TFs) orchestrate gene expression programs crucial for brain function, but we lack detailed information about TF binding in human brain tissue. We generated a multiomic resource (ChIP-seq, ATAC-seq, RNA-seq, DNA methylation) on bulk tissues and sorted nuclei from several postmortem brain regions, including binding maps for more than 100 TFs. We demonstrate improved measurements of TF activity, including motif recognition and gene expression modeling, upon identification and removal of high TF occupancy regions. Further, predictive TF binding models demonstrate a bias for these high-occupancy sites. Neuronal TFs SATB2 and TBR1 bind unique regions depleted for such sites and promote neuronal gene expression. Binding sites for TFs, including TBR1 and PKNOX1, are enriched for risk variants associated with neuropsychiatric disorders, predominantly in neurons. This work, titled BrainTF, is a powerful resource for future studies seeking to understand the roles of specific TFs in regulating gene expression in the human brain.

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.

Suicide is a condition resulting from complex environmental and genetic risks that affect millions of people globally. Both structural and functional studies identified the hippocampus as one of the vulnerable brain regions contributing to suicide risk. Here, we have identified the hippocampal transcriptomes, gene ontology, cell type proportions, dendritic spine morphology, and transcriptomic signature in iPSC-derived neuronal precursor cells (NPCs) and neurons in postmortem brain tissue from suicide deaths. The hippocampal tissue transcriptomic data revealed that NPAS4 gene expression was downregulated while ALDH1A2, NAAA, and MLXIPL gene expressions were upregulated in tissue from suicide deaths. The gene ontology identified 29 significant pathways including NPAS4-associated gene ontology terms "excitatory post-synaptic potential", "regulation of postsynaptic membrane potential" and "long-term memory" indicating alteration of glutamatergic synapses in the hippocampus of suicide deaths. The cell type deconvolution identified decreased excitatory neuron proportion and an increased inhibitory neuron proportion providing evidence of excitation/inhibition imbalance in the hippocampus of suicide deaths. In addition, suicide deaths had increased dendric spine density, due to an increase of thin (relatively unstable) dendritic spines, compared to controls. The transcriptomes of iPSC-derived hippocampal-like NPCs and neurons revealed 31 and 33 differentially expressed genes in NPC and neurons, respectively, of suicide deaths. The suicide-associated differentially expressed genes in NPCs were RELN, CRH, EMX2, OXTR, PARM1 and IFITM2 which overlapped with previously published results. The previously-known suicide-associated differentially expressed genes in differentiated neurons were COL1A1, THBS1, IFITM2, AQP1, and NLRP2. Together, these findings would help better understand the hippocampal neurobiology of suicide for identifying therapeutic targets to prevent suicide.

Objective: This study assessed the proteomic profiles of cytokines and chemokines in individuals with moderate to severe depression, with or without comorbid medical disorders, compared to healthy controls. Two proteomic multiplex platforms were employed for this purpose. Metods: An immunofluorescent multiplex platform and an aptamer-based method were used to evaluate 32 protein analytes from 153 individuals with moderate to severe major depressive disorder (MDD) and healthy controls (HCs). The study focused on determining the level of agreement between the two platforms and evaluating the ability of individual analytes and principal components (PCs) to differentiate between the MDD and HC groups. Additionally, the study investigated the relationship between PCs consisting of chemokines and cytokines and comorbid inflammatory and cardiometabolic diseases. Findings: Analysis revealed a small or moderate correlation between 47% of the analytes measured by the two platforms. Two proteomic profiles were identified that differentiated individuals with moderate to severe MDD from HCs. High eotaxin, age, BMI, IP-10, or IL-10 characterized profile 1. This profile was associated with several cardiometabolic risk factors, including hypertension, hyperlipidemia, and type 2 diabetes. Profile 2 is characterized by higher age, BMI, interleukins, and a strong negative loading for eotaxin. This profile was associated with inflammation but not cardiometabolic risk factors. Conclusion: This study provides further evidence that proteomic profiles can be used to identify potential biomarkers and pathways associated with MDD and comorbidities. Our findings suggest that MDD is associated with distinct profiles of proteins that are also associated with cardiometabolic risk factors, inflammation, and obesity. In particular, the chemokines eotaxin and IP-10 appear to play a role in the relationship between MDD and cardiometabolic risk factors. These findings suggest that a focus on the interplay between MDD and comorbidities may be useful in identifying potential targets for intervention and improving overall health outcomes.

Transcription factors (TFs) orchestrate gene expression programs crucial for cell physiology, but our knowledge of their function in the brain is limited. Using bulk tissues and sorted nuclei from multiple human post-mortem brain regions, we generated a multi-omic resource (1121 total experiments) that includes binding maps for more than 100 TFs. We demonstrate improved measurements of TF activity, including motif recognition and gene expression modeling, upon identification and removal of regions of high TF occupancy. Further, we find that predictive TF binding models demonstrate a bias for these high occupancy sites. Neuronal TFs SATB2 and TBR1 bind unique regions depleted for such sites and promote neuronal gene expression. Several TFs, including TBR1 and PKNOX1, are enriched for risk variants associated with neuropsychiatric disorders, predominantly in neurons. These data are a powerful resource for future studies seeking to understand the role of TFs in epigenetic regulation in the human brain.

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.

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.

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.

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.

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.

A major challenge for depression research is understanding the delicate balance of neuronal and glial systems regulating emotion. This presentation outlines several lines of evidence suggesting that glia play critical roles in affective vulnerability and temperament.

The effects of smoking in the human brain were explored at the molecular level in the nucleus accumbens (NAcc), a key brain region involved in tobacco addiction. Gene expression data from post-mortem NAcc were analyzed according to smoking habits: Never smokers, Former smokers and Current smokers at the time of death. The effect of smoking was determined using an ANCOVA model, controlling for potential confounders (psychiatric diagnosis, gender, age, post-mortem interval, and brain pH) followed by pair-wise post-hoc comparisons. Q-values (false discovery rate adjusted p-values) < 0.05 were used in combination with a fold change of > ±1.3 to identify the most relevant genes. The greatest number of differentially expressed genes (DEGs) were found in subjects with a recent history of smoking (Current smokers) compared to either Former or Never smokers. Only two genes were differentially expressed between Former and Never smokers, suggesting that the effects of smoking on gene expression in the brain may be transient. Ingenuity Pathway Analysis (IPA) of DEGs identified a significant over-representation of neurotransmitter system genes (glutamate, GABA) in Current smokers. IPA also revealed many genes associated with cancer in Current smokers compared to Former and Never smokers despite no known cancer in any subjects. Genes associated with neoplasms, glioblastoma, gliomas and tumor regulations are among the top 10 transcripts.Our findings show that active smokers have a significant increase in cancer-related genes and alterations in glutamate and GABA neurotransmitter systems in the NAcc. To our knowledge this is the first study to identify cancer-related genes in the NAcc in Current smokers who have no evidence of cancer. ### Competing Interest Statement The authors have declared no competing interest. ### Funding Statement The study was funded by the Pritzker Family Philanthropic Fund, and the American Foundation for Suicide Prevention Standard Research Grant SRG-0-121-18 (to A.S.) The authors declare no other financial or non-financial competing interests. ### Author Declarations I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained. Yes The details of the IRB/oversight body that provided approval or exemption for the research described are given below: The Institutional Review Board of the University of California, Irvine gave ethical approval for this work. I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals. Yes I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance). Yes I have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable. Yes All data produced in the present study are available upon reasonable request to the authors.

Schizophrenia (SZ) and bipolar disorder (BD) are psychiatric disorders with high genetic correlation and heritability. Mitochondrial dysfunctions in synaptic, energetic, and metabolic pathways have been associated with the neurobiology of both SZ and BD. A common variant of the mtDNA A15218G allele was reported in four studies (SZ = 10371, control = 34595) with suggestive association with SZ (p = 2.15E-03, OR = 1.2, 95% CI 1.07-1.35). Gene-based and pathway-based approaches show that genes in Complex I are of particular importance in informing the risk of SZ.

A growing number of epidemiological and experimental studies has established that circadian disruption is strongly associated with psychiatric disorders, including major depressive disorder (MDD). This association is becoming increasingly relevant considering that modern lifestyles, social zeitgebers (time cues) and genetic variants contribute to disrupting circadian rhythms that may lead to psychiatric disorders. Circadian abnormalities associated with MDD include dysregulated rhythms of sleep, temperature, hormonal secretions, and mood which are modulated by the molecular clock. Rapid-acting antidepressants such as subanesthetic ketamine and sleep deprivation therapy can improve symptoms within 24 h in a subset of depressed patients, in striking contrast to conventional treatments, which generally require weeks for a full clinical response. Importantly, animal data show that sleep deprivation and ketamine have overlapping effects on clock gene expression. Furthermore, emerging data implicate the circadian system as a critical component involved in rapid antidepressant responses via several intracellular signaling pathways such as GSK3β, mTOR, MAPK, and NOTCH to initiate synaptic plasticity. Future research on the relationship between depression and the circadian clock may contribute to the development of novel therapeutic strategies for depression-like symptoms. In this review we summarize recent evidence describing: (1) how the circadian clock is implicated in depression, (2) how clock genes may contribute to fast-acting antidepressants, and (3) the mechanistic links between the clock genes driving circadian rhythms and neuroplasticity.

The study of postsynaptic excitation to inhibition (E/I ratio) imbalances in human brain diseases, is a highly relevant functional measurement poorly investigated due to postmortem degradation of synaptic receptors. We show that near-simultaneous recording of microtransplanted synaptic receptors after simulated morgue conditions allows the determination of the postsynaptic E/I ratio for at least 120 h after death, expanding the availability and use of human diseased tissue stored in brain banks.

Clinical and preclinical studies suggest that some of the behavioral alterations observed in schizophrenia (SZ) may be mechanistically linked to synaptic dysfunction of glutamatergic signaling. Recent genetic and proteomic studies suggest alterations of cortical glutamate receptors of the AMPA-type (AMPARs), which are the predominant ligand-gated ionic channels of fast transmission at excitatory synapses. The impact of gene and protein alterations on the electrophysiological activity of AMPARs is not known in SZ. In this proof of principle work, using human postmortem brain synaptic membranes isolated from the dorsolateral prefrontal cortex (DLPFC), we combined electrophysiological analysis from microtransplanted synaptic membranes (MSM) with transcriptomic (RNA-Seq) and label-free proteomics data in 10 control and 10 subjects diagnosed with SZ. We observed in SZ a reduction in the amplitude of AMPARs currents elicited by kainate, an agonist of AMPARs that blocks the desensitization of the receptor. This reduction was not associated with protein abundance but with a reduction in kainate's potency to activate AMPARs. Electrophysiologically-anchored dataset analysis (EDA) was used to identify synaptosomal proteins that linearly correlate with the amplitude of the AMPARs responses, gene ontology functional annotations were then used to determine protein-protein interactions. Protein modules associated with positive AMPARs current increases were downregulated in SZ, while protein modules that were upregulated in SZ were associated with decreased AMPARs currents. Our results indicate that transcriptomic and proteomic alterations, frequently observed in the DLPFC in SZ, converge at the synaptic level producing a functional electrophysiological impairment of AMPARs.

Accumulating evidence points to a significant link between disrupted circadian rhythms and neuronal disfunctions, though the molecular mechanisms underlying this connection are virtually unexplored. The transcript Homer1a, an immediate early gene related to postsynaptic signaling, has been demonstrated to exhibit robust circadian oscillation in the brain, which supports the hypothesis that Homer1a mediates the communication between circadian inputs and neuronal activity. Here, we determined how the circadian clock is implicated in Homer1a gene regulation by using circadian clock Bmal1-mutant mice either in the presence or absence of stress stimulation. The Homer1 gene generates multiple transcripts, but only the short variant Homer1a responds to acute stress with sleep deprivation (SD) in mice. Chromatin immunoprecipitation assays revealed that both transcription factor CREB and the circadian clock component BMAL1 bind to the Homer1 promoter in mouse brain. Importantly, circadian Homer1a gene expression is unaltered in the absence of BMAL1, while its immediate early response to SD relies on BMAL1. Deletion of Bmal1 results in attenuated CREB activity in mouse brain, which appears to contribute to decreased expression of Homer1a in response to SD. In conclusion, Homer1a undergoes bimodal control by the circadian clock and CREB.