Multiomic profiling of transcription factor binding and function in human brain.

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.

Proteomic profiles of cytokines and chemokines in moderate to severe depression: Implications for comorbidities and biomarker discovery.

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.

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.

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.

Altered transcriptomes, cell type proportions, and dendritic spine morphology in hippocampus of suicide deaths.

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.

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.

Functional impairment of cortical AMPA receptors in schizophrenia.

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.

Homer1a Undergoes Bimodal Transcriptional Regulation by CREB and the Circadian Clock.

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.

A study of psychological pain in substance use disorder and its relationship to treatment outcome.

Substance Use Disorder (SUD) is a major public health concern affecting an estimated 22.5 million individuals in the United States. The primary aim of this study was to characterize psychological pain in a cohort of patients participating in outpatient treatment for SUD. A secondary aim was to determine the relationships between pre-treatment assessments of psychological pain, depression, anxiety and hopelessness with treatment retention time and completion rates. Data was analyzed from 289 patients enrolled in an outpatient community drug treatment clinic in Southern California, U.S. A previously determined threshold score on the Mee-Bunney Psychological Pain Assessment Scale (MBP) was utilized to group patients into high and low-moderate scoring subgroups. The higher pain group scored higher on measures of anxiety, hopelessness and depression compared to those in the low-moderate pain group. Additionally, patients scoring in the higher psychological pain group exhibited reduced retention times in treatment and more than two-fold increased odds of dropout relative to patients with lower pre-treatment levels of psychological pain. Among all assessments, the correlation between psychological pain and treatment retention time was strongest. To our knowledge, this is the first study to demonstrate that psychological pain is an important construct which correlates with relevant clinical outcomes in SUD. Furthermore, pre-treatment screening for psychological pain may help target higher-risk patients for clinical interventions aimed at improving treatment retention and completion rates.

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.

Usage of Therapeutic Sleep Deprivation: A Survey in Psychiatric Hospitals in Austria, Germany, and Switzerland.

Objective: Therapeutic sleep deprivation (SD) is a nonpharmacological treatment that is used most often for depression. The aim of this study was to examine the pattern of use of SD in psychiatric hospitals in Austria, Germany, and Switzerland. Methods: A questionnaire about perceived usage of SD was sent by mail to all 511 psychiatric hospitals in the three countries. Nonresponders were asked to answer the questionnaire by phone. We achieved a response rate of 75.3%. Results: SD was recommended by 61.3% of all hospitals. Despite this degree of recommendation, nearly two thirds of the psychiatric hospitals had not treated a patient with SD during the last 12 months. Of the respondents, 59.5% considered SD to be indicated for major depressive disorder, 17.7% for bipolar depression, and 7.8% for other indications. SD was administered most frequently in inpatient settings and in combination with other therapies. Total SD (patients kept awake entire night) and partial late SD (patients kept awake in second half of night) were judged equally effective. Of the hospitals, 53.0% reported having seen hypomania and 13.2% manic episodes as side effects (rates do not represent actual incident rates). Conclusion: The lack of large controlled studies for SD with its different forms of treatment probably still hinders a broader use of the therapy. Therefore, further efforts should be undertaken to provide high-quality scientific evidence for the usage of SD.

Inference of cell type content from human brain transcriptomic datasets illuminates the effects of age, manner of death, dissection, and psychiatric diagnosis.

Psychiatric illness is unlikely to arise from pathology occurring uniformly across all cell types in affected brain regions. Despite this, transcriptomic analyses of the human brain have typically been conducted using macro-dissected tissue due to the difficulty of performing single-cell type analyses with donated post-mortem brains. To address this issue statistically, we compiled a database of several thousand transcripts that were specifically-enriched in one of 10 primary cortical cell types in previous publications. Using this database, we predicted the relative cell type content for 833 human cortical samples using microarray or RNA-Seq data from the Pritzker Consortium (GSE92538) or publicly-available databases (GSE53987, GSE21935, GSE21138, CommonMind Consortium). These predictions were generated by averaging normalized expression levels across transcripts specific to each cell type using our R-package BrainInABlender (validated and publicly-released on github). Using this method, we found that the principal components of variation in the datasets strongly correlated with the predicted neuronal/glial content of the samples. This variability was not simply due to dissection-the relative balance of brain cell types appeared to be influenced by a variety of demographic, pre- and post-mortem variables. Prolonged hypoxia around the time of death predicted increased astrocytic and endothelial gene expression, illustrating vascular upregulation. Aging was associated with decreased neuronal gene expression. Red blood cell gene expression was reduced in individuals who died following systemic blood loss. Subjects with Major Depressive Disorder had decreased astrocytic gene expression, mirroring previous morphometric observations. Subjects with Schizophrenia had reduced red blood cell gene expression, resembling the hypofrontality detected in fMRI experiments. Finally, in datasets containing samples with especially variable cell content, we found that controlling for predicted sample cell content while evaluating differential expression improved the detection of previously-identified psychiatric effects. We conclude that accounting for cell type can greatly improve the interpretability of transcriptomic data.

Connective Tissue Growth Factor Is a Novel Prodepressant.

BACKGROUND: While downregulation of several growth factors in major depressive disorder is well established, less attention has been paid to the upregulation of other growth factors. Yet, upregulated growth factors may offer better therapeutic targets. We show that connective tissue growth factor (CTGF) represents a target based on its upregulation in major depressive disorder and studies in animal models implicating it in negative affect. METHODS: CTGF gene expression was first evaluated in the postmortem human amygdala. The findings were followed up in outbred rats and in two rat lines that were selectively bred for differences in novelty-seeking and anxiety behavior (bred low responders and bred high responders). We studied the impact of social defeat and early-life treatment with fibroblast growth factor 2 on CTGF expression. Finally, we assessed the ability of an anti-CTGF antibody (FG-3019) to alter CTGF expression and emotionality. RESULTS: In the human amygdala, CTGF expression was significantly increased in major depressive disorder compared with control subjects. CTGF expression was also significantly increased in the dentate gyrus of adult bred low responders compared with bred high responders. Social defeat stress in bred low responders significantly increased CTGF expression in the dentate gyrus. Early-life fibroblast growth factor 2, a treatment that reduces anxiety-like behavior throughout life, decreased CTGF expression in the adult dentate gyrus. In outbred rats, CTGF administration increased depression-like behavior. Chronic treatment with FG-3019 decreased CTGF expression, and acute and chronic treatment was antidepressant. CONCLUSIONS: This study is the first to implicate CTGF as a prodepressant molecule that could serve as a target for the development of novel therapeutics.

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.

Post-mortem molecular profiling of three psychiatric disorders.

BACKGROUND: Psychiatric disorders are multigenic diseases with complex etiology that contribute significantly to human morbidity and mortality. Although clinically distinct, several disorders share many symptoms, suggesting common underlying molecular changes exist that may implicate important regulators of pathogenesis and provide new therapeutic targets. METHODS: We performed RNA sequencing on tissue from the anterior cingulate cortex, dorsolateral prefrontal cortex, and nucleus accumbens from three groups of 24 patients each diagnosed with schizophrenia, bipolar disorder, or major depressive disorder, and from 24 control subjects. We identified differentially expressed genes and validated the results in an independent cohort. Anterior cingulate cortex samples were also subjected to metabolomic analysis. ChIP-seq data were used to characterize binding of the transcription factor EGR1. RESULTS: We compared molecular signatures across the three brain regions and disorders in the transcriptomes of post-mortem human brain samples. The most significant disease-related differences were in the anterior cingulate cortex of schizophrenia samples compared to controls. Transcriptional changes were assessed in an independent cohort, revealing the transcription factor EGR1 as significantly down-regulated in both cohorts and as a potential regulator of broader transcription changes observed in schizophrenia patients. Additionally, broad down-regulation of genes specific to neurons and concordant up-regulation of genes specific to astrocytes was observed in schizophrenia and bipolar disorder patients relative to controls. Metabolomic profiling identified disruption of GABA levels in schizophrenia patients. CONCLUSIONS: We provide a comprehensive post-mortem transcriptome profile of three psychiatric disorders across three brain regions. We highlight a high-confidence set of independently validated genes differentially expressed between schizophrenia and control patients in the anterior cingulate cortex and integrate transcriptional changes with untargeted metabolite profiling.