Common mitochondrial deletions in RNA-Seq: evaluation of bulk, single-cell, and spatial transcriptomic datasets

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.

The relationship between case-control differential gene expression from brain tissue and genetic associations in schizophrenia

Large numbers of genetic loci have been identified that are known to contain common risk alleles for schizophrenia, but linking associated alleles to specific risk genes remains challenging. Given that most alleles that influence liability to schizophrenia are thought to do so by altered gene expression, intuitively, case-control differential gene expression studies should highlight genes with a higher probability of being associated with schizophrenia and could help identify the most likely causal genes within associated loci. Here, we test this hypothesis by comparing transcriptome analysis of the dorsolateral prefrontal cortex from 563 schizophrenia cases and 802 controls with genome-wide association study (GWAS) data from the third wave study of the Psychiatric Genomics Consortium. Genes differentially expressed in schizophrenia were not enriched for common allelic association statistics compared with other brain-expressed genes, nor were they enriched for genes within associated loci previously reported to be prioritized by genetic fine-mapping. Genes prioritized by Summary-based Mendelian Randomization were underexpressed in cases compared to other genes in the same GWAS loci. However, the overall strength and direction of expression change predicted by SMR were not related to that observed in the differential expression data. Overall, this study does not support the hypothesis that genes identified as differentially expressed from RNA sequencing of bulk brain tissue are enriched for those that show evidence for genetic associations. Such data have limited utility for prioritizing genes in currently associated loci in schizophrenia.

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.

Antipsychotic drug use complicates assessment of gene expression changes associated with schizophrenia

Recent postmortem transcriptomic studies of schizophrenia (SCZ) have shown hundreds of differentially expressed genes. However, the extent to which these gene expression changes reflect antipsychotic drug (APD) exposure remains uncertain. We compared differential gene expression in the prefrontal cortex of SCZ patients who tested positive for APDs at the time of death with SCZ patients who did not. APD exposure was associated with numerous changes in the brain transcriptome, especially among SCZ patients on atypical APDs. Brain transcriptome data from macaques chronically treated with APDs showed that APDs affect the expression of many functionally relevant genes, some of which show expression changes in the same directions as those observed in SCZ. Co-expression modules enriched for synaptic function showed convergent patterns between SCZ and some of the APD effects, while those associated with inflammation and glucose metabolism exhibited predominantly divergent patterns between SCZ and APD effects. In contrast, major cell-type shifts inferred in SCZ were primarily unaffected by APD use. These results show that APDs may confound SCZ-associated gene expression changes in postmortem brain tissue. Disentangling these effects will help identify causal genes and improve our neurobiological understanding of SCZ.

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.

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.

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.

Mapping genomic loci implicates genes and synaptic biology in schizophrenia

Schizophrenia has a heritability of 60-80%1, much of which is attributable to common risk alleles. Here, in a two-stage genome-wide association study of up to 76,755 individuals with schizophrenia and 243,649 control individuals, we report common variant associations at 287 distinct genomic loci. Associations were concentrated in genes that are expressed in excitatory and inhibitory neurons of the central nervous system, but not in other tissues or cell types. Using fine-mapping and functional genomic data, we identify 120 genes (106 protein-coding) that are likely to underpin associations at some of these loci, including 16 genes with credible causal non-synonymous or untranslated region variation. We also implicate fundamental processes related to neuronal function, including synaptic organization, differentiation and transmission. Fine-mapped candidates were enriched for genes associated with rare disruptive coding variants in people with schizophrenia, including the glutamate receptor subunit GRIN2A and transcription factor SP4, and were also enriched for genes implicated by such variants in neurodevelopmental disorders. We identify biological processes relevant to schizophrenia pathophysiology; show convergence of common and rare variant associations in schizophrenia and neurodevelopmental disorders; and provide a resource of prioritized genes and variants to advance mechanistic studies.

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.

Genome-wide association study of more than 40,000 bipolar disorder cases provides new insights into the underlying biology

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.

GRK3 deficiency elicits brain immune activation and psychosis

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.

Genome-Wide Association Studies of Schizophrenia and Bipolar Disorder in a Diverse Cohort of US Veterans

BACKGROUND

Schizophrenia (SCZ) and bipolar disorder (BIP) are debilitating neuropsychiatric disorders, collectively affecting 2% of the world's population. Recognizing the major impact of these psychiatric disorders on the psychosocial function of more than 200 000 US Veterans, the Department of Veterans Affairs (VA) recently completed genotyping of more than 8000 veterans with SCZ and BIP in the Cooperative Studies Program (CSP) #572.

METHODS

We performed genome-wide association studies (GWAS) in CSP #572 and benchmarked the predictive value of polygenic risk scores (PRS) constructed from published findings. We combined our results with available summary statistics from several recent GWAS, realizing the largest and most diverse studies of these disorders to date.

RESULTS

Our primary GWAS uncovered new associations between CHD7 variants and SCZ, and novel BIP associations with variants in Sortilin Related VPS10 Domain Containing Receptor 3 (SORCS3) and downstream of PCDH11X. Combining our results with published summary statistics for SCZ yielded 39 novel susceptibility loci including CRHR1, and we identified 10 additional findings for BIP (28 326 cases and 90 570 controls). PRS trained on published GWAS were significantly associated with case-control status among European American (P < 10-30) and African American (P < .0005) participants in CSP #572.

CONCLUSIONS

We have demonstrated that published findings for SCZ and BIP are robustly generalizable to a diverse cohort of US veterans. Leveraging available summary statistics from GWAS of global populations, we report 52 new susceptibility loci and improved fine-mapping resolution for dozens of previously reported associations.

Contributions of common genetic variants to risk of schizophrenia among individuals of African and Latino ancestry

Schizophrenia is a common, chronic and debilitating neuropsychiatric syndrome affecting tens of millions of individuals worldwide. While rare genetic variants play a role in the etiology of schizophrenia, most of the currently explained liability is within common variation, suggesting that variation predating the human diaspora out of Africa harbors a large fraction of the common variant attributable heritability. However, common variant association studies in schizophrenia have concentrated mainly on cohorts of European descent. We describe genome-wide association studies of 6152 cases and 3918 controls of admixed African ancestry, and of 1234 cases and 3090 controls of Latino ancestry, representing the largest such study in these populations to date. Combining results from the samples with African ancestry with summary statistics from the Psychiatric Genomics Consortium (PGC) study of schizophrenia yielded seven newly genome-wide significant loci, and we identified an additional eight loci by incorporating the results from samples with Latino ancestry. Leveraging population differences in patterns of linkage disequilibrium, we achieve improved fine-mapping resolution at 22 previously reported and 4 newly significant loci. Polygenic risk score profiling revealed improved prediction based on trans-ancestry meta-analysis results for admixed African (Nagelkerke's R2 = 0.032; liability R2 = 0.017; P < 10-52), Latino (Nagelkerke's R2 = 0.089; liability R2 = 0.021; P < 10-58), and European individuals (Nagelkerke's R2 = 0.089; liability R2 = 0.037; P < 10-113), further highlighting the advantages of incorporating data from diverse human populations.

Preservation of global synaptic excitatory to inhibitory ratio during long postmortem intervals

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.

Psychiatric drugs impact mitochondrial function in brain and other tissues

Mitochondria have been linked to the etiology of schizophrenia (SZ). However, studies of mitochondria in SZ might be confounded by the effects of pharmacological treatment with antipsychotic drugs (APDs) and other common medications. This review summarizes findings on relevant mitochondria mechanisms underlying SZ, and the potential impact of psychoactive drugs including primarily APDs, but also antidepressants and anxiolytics. The summarized data suggest that APDs impair mitochondria function by decreasing Complex I activity and ATP production and dissipation of the mitochondria membrane potential. At the same time, in the brains of patients with SZ, antipsychotic drug treatment normalizes gene expression modules enriched in mitochondrial genes that are decreased in SZ. This indicates that APDs may have both positive and negative effects on mitochondria. The available evidence suggests three conclusions i) alterations in mitochondria functions in SZ exist prior to APD treatment, ii) mitochondria alterations in SZ can be reversed by APD treatment, and iii) APDs directly cause impairment of mitochondria function. Overall, the mechanisms of action of psychiatric drugs on mitochondria are both direct and indirect; we conclude the effects of APDs on mitochondria may contribute to both their therapeutic and metabolic side effects. These studies support the hypothesis that neuronal mitochondria are an etiological factor in SZ. Moreover, APDs and other drugs must be considered in the evaluation of this pathophysiological role of mitochondria in SZ. Considering these effects, pharmacological actions on mitochondria may be a worthwhile target for further APD development.

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.

Association of Myoinositol Transporters with Schizophrenia and Bipolar Disorder: Evidence from Human and Animal Studies

Evidence from animal and human studies has linked myo-inositol (MI) with the pathophysiology and/or treatment of psychiatric disorders such as schizophrenia and bipolar disorder. However, there is still controversy surrounding the definitive role of MI in these disorders. Given that brain MI is differentially regulated by three transporters - SMIT1, SMIT2 and/or HMIT (encoded by the genes: SLC5A3, SLC5A11, and SLC2A13, respectively) - we used available datasets to describe the distribution in mouse and human brain of the different MI transporters and to examine changes in mRNA expression of these transporters in patients with schizophrenia and bipolar disorder. We found a differential distribution of the mRNA of each of the three MI transporters in both human and mouse brain regions. Interestingly, while individual neurons express SMIT1 and HMIT, non-neuronal cells express SMIT2, thus partially accounting for different uptake levels of MI and concordance to downstream second messenger signaling pathways. We also found that the expression of MI transporters is significantly changed in schizophrenia and bipolar disorder in a diagnostic-, brain region- and subtype-specific manner. We then examined the effects of germline deletion in mice of Slc5a3 on behavioral phenotypes related to schizophrenia and bipolar disorder. This gene deletion produces behavioral deficits that mirror some specific symptoms of schizophrenia and bipolar disorder. Finally, chronic administration of MI was able to reverse particular, but not all, behavioral deficits in Slc5a3 knockout mice; MI itself induced some behavioral deficits. Our data support a strong correlation between the expression of MI transporters and schizophrenia and bipolar disorder, and suggest that brain region-specific aberration of one or more of these transporters determines the partial behavioral phenotypes and/or symptomatic pattern of these disorders.

Melanin Concentrating Hormone Signaling Deficits in Schizophrenia: Association With Memory and Social Impairments and Abnormal Sensorimotor Gating

BACKGROUND

Evidence from anatomical, pharmacological, and genetic studies supports a role for the neuropeptide melanin concentrating hormone system in modulating emotional and cognitive functions. Genome-wide association studies revealed a potential association between the melanin concentrating hormone receptor (MCHR1) gene locus and schizophrenia, and the largest genome-wide association study conducted to date shows a credible genome-wide association.

METHODS

We analyzed MCHR1 and pro-melanin concentrating hormone RNA-Seq expression in the prefrontal cortex in schizophrenia patients and healthy controls. Disruptions in the melanin concentrating hormone system were modeled in the mouse brain by germline deletion of MCHR1 and by conditional ablation of melanin concentrating hormone expressing neurons using a Cre-inducible diphtheria toxin system.

RESULTS

MCHR1 expression is decreased in the prefrontal cortex of schizophrenia samples (false discovery rate (FDR) P < .05, CommonMind and PsychEncode combined datasets, n = 901) while pro-melanin concentrating hormone is below the detection threshold. MCHR1 expression decreased with aging (P = 6.6E-57) in human dorsolateral prefrontal cortex. The deletion of MCHR1 was found to lead to behavioral abnormalities mimicking schizophrenia-like phenotypes: hyperactivity, increased stereotypic and repetitive behavior, social impairment, impaired sensorimotor gating, and disrupted cognitive functions. Conditional ablation of pro-melanin concentrating hormone neurons increased repetitive behavior and produced a deficit in sensorimotor gating.

CONCLUSIONS

Our study indicates that early disruption of the melanin concentrating hormone system interferes with neurodevelopmental processes, which may contribute to the pathogenesis of schizophrenia. Further neurobiological research on the developmental timing and circuits that are affected by melanin concentrating hormone may lead to a therapeutic target for early prevention of schizophrenia.

Novel Complex Interactions between Mitochondrial and Nuclear DNA in Schizophrenia and Bipolar Disorder

Mitochondrial dysfunction has been associated with schizophrenia (SZ) and bipolar disorder (BD). This review examines recent publications and novel associations between mitochondrial genes and SZ and BD. Associations of nuclear-encoded mitochondrial variants with SZ were found using gene- and pathway-based approaches. Two control region mitochondrial DNA (mtDNA) SNPs, T16519C and T195C, both showed an association with SZ and BD. A review of 4 studies of A15218G located in the cytochrome B oxidase gene (CYTB, SZ = 11,311, control = 35,735) shows a moderate association with SZ (p = 2.15E-03). Another mtDNA allele A12308G was nominally associated with psychosis in BD type I subjects and SZ. The first published study testing the epistatic interaction between nuclear-encoded and mitochondria-encoded genes demonstrated evidence for potential interactions between mtDNA and the nuclear genome for BD. A similar analysis for the risk of SZ revealed significant joint effects (34 nuclear-mitochondria SNP pairs with joint effect p ≤ 5E-07) and significant enrichment of projection neurons. The mitochondria-encoded gene CYTB was found in both the epistatic interactions for SZ and BD and the single SNP association of SZ. Future efforts considering population stratification and polygenic risk scores will test the role of mitochondrial variants in psychiatric disorders.

Super-Obese Patient-Derived iPSC Hypothalamic Neurons Exhibit Obesogenic Signatures and Hormone Responses

The hypothalamus contains neurons that integrate hunger and satiety endocrine signals from the periphery and are implicated in the pathophysiology of obesity. The limited availability of human hypothalamic neurons hampers our understanding of obesity disease mechanisms. To address this, we generated human induced pluripotent stem cells (hiPSCs) from multiple normal body mass index (BMI; BMI ≤ 25) subjects and super-obese (OBS) donors (BMI ≥ 50) with polygenic coding variants in obesity-associated genes. We developed a method to reliably differentiate hiPSCs into hypothalamic-like neurons (iHTNs) capable of secreting orexigenic and anorexigenic neuropeptides. Transcriptomic profiling revealed that, although iHTNs maintain a fetal identity, they respond appropriately to metabolic hormones ghrelin and leptin. Notably, OBS iHTNs retained disease signatures and phenotypes of high BMI, exhibiting dysregulated respiratory function, ghrelin-leptin signaling, axonal guidance, glutamate receptors, and endoplasmic reticulum (ER) stress pathways. Thus, human iHTNs provide a powerful platform to study obesity and gene-environment interactions.

A Comprehensive Analysis of Nuclear-Encoded Mitochondrial Genes in Schizophrenia

BACKGROUND

The genetic risk factors of schizophrenia (SCZ), a severe psychiatric disorder, are not yet fully understood. Multiple lines of evidence suggest that mitochondrial dysfunction may play a role in SCZ, but comprehensive association studies are lacking. We hypothesized that variants in nuclear-encoded mitochondrial genes influence susceptibility to SCZ.

METHODS

We conducted gene-based and gene-set analyses using summary association results from the Psychiatric Genomics Consortium Schizophrenia Phase 2 (PGC-SCZ2) genome-wide association study comprising 35,476 cases and 46,839 control subjects. We applied the MAGMA method to three sets of nuclear-encoded mitochondrial genes: oxidative phosphorylation genes, other nuclear-encoded mitochondrial genes, and genes involved in nucleus-mitochondria crosstalk. Furthermore, we conducted a replication study using the iPSYCH SCZ sample of 2290 cases and 21,621 control subjects.

RESULTS

In the PGC-SCZ2 sample, 1186 mitochondrial genes were analyzed, among which 159 had p values < .05 and 19 remained significant after multiple testing correction. A meta-analysis of 818 genes combining the PGC-SCZ2 and iPSYCH samples resulted in 104 nominally significant and nine significant genes, suggesting a polygenic model for the nuclear-encoded mitochondrial genes. Gene-set analysis, however, did not show significant results. In an in silico protein-protein interaction network analysis, 14 mitochondrial genes interacted directly with 158 SCZ risk genes identified in PGC-SCZ2 (permutation p = .02), and aldosterone signaling in epithelial cells and mitochondrial dysfunction pathways appeared to be overrepresented in this network of mitochondrial and SCZ risk genes.

CONCLUSIONS

This study provides evidence that specific aspects of mitochondrial function may play a role in SCZ, but we did not observe its broad involvement even using a large sample.

Examining the role of common and rare mitochondrial variants in schizophrenia

Oxidative phosphorylation within mitochondria is the main source of aerobic energy for neuronal functioning, and the key genes are located in mitochondrial DNA. Deficits in oxidative phosphorylation functioning have been reported for schizophrenia, but efforts in the identification of genetic markers within the mitochondrial DNA that predispose to schizophrenia have been limited. We genotyped a set of mitochondrial SNPs using Illumina HumanExome arrays and tested for association in the Swedish schizophrenia sample (N> 10,000). We developed a novel approach for mitochondrial DNA imputation in order to increase the number of common SNPs available for association analysis. The most significant findings were for the mitochondrial SNPs C15452A (GRCh38.p10; rs527236209; p = 0.007; gene MT-CYB; defining haplogroup JT); A11251G (rs869096886; p = 0.007; gene MT-ND4; defining haplogroup JT), and T4216C (rs1599988; p = 0.008, gene MT-ND1, defining haplogroup R2'JT). We also conducted rare variant burden analyses and obtained a p-value of 0.007. For multimarker haplotypes analysis, the most significant finding was for the J group (OR: 0.86, p = 0.02). We conducted the largest association study of mitochondrial DNA variants and schizophrenia but did not find an association that survived multiple testing correction. Analysis of a larger sample is required and will allow a better understanding of the role of mitochondria in schizophrenia.

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.

A Genome-Wide Search for Bipolar Disorder Risk Loci Modified by Mitochondrial Genome Variation

Mitochondrial DNA mutations have been reported to be associated with bipolar disorder (BD). In this study, we performed genome-wide analyses to assess mitochondrial single-nucleotide polymorphism (mtSNP) effects on BD risk and early-onset BD (EOBD) among BD patients, focusing on interaction effects between nuclear SNPs (nSNPs) and mtSNPs. Common nSNP and mtSNP data from European American BD cases (n = 1,001) and controls (n = 1,034) from the Genetic Association Information Network BD study were analyzed to assess the joint effect of nSNP and nSNP-mtSNP interaction on the risk of BD and EOBD. The effect of nSNP-mtSNP interactions was also assessed. For BD risk, the strongest evidence of an association was obtained for nSNP rs1880924 in MGAM and mtSNP rs3088309 in CytB (p_joint = 8.2 × 10^-8, p_int = 1.4 × 10^-4). Our results also suggest that the minor allele of the nSNP rs583990 in CTNNA2 increases the risk of EOBD among carriers of the mtSNP rs3088309 minor allele, while the nSNP has no effect among those carrying the mtSNP major allele (OR = 4.53 vs. 1.05, p_joint = 2.1 × 10^-7, p_int = 1.16 × 10^-6). While our results are not statistically significant after multiple testing correction and a large-sample replication is required, our exploratory study demonstrates the potential importance of considering the mitochondrial genome for identifying genetic factors associated with BD.

A Circadian Genomic Signature Common to Ketamine and Sleep Deprivation in the Anterior Cingulate Cortex

BACKGROUND

Conventional antidepressants usually require several weeks to achieve a full clinical response in patients with major depressive disorder, an illness associated with dysregulated circadian rhythms and a high incidence of suicidality. Two rapid-acting antidepressant strategies, low-dose ketamine (KT) and sleep deprivation (SD) therapies, dramatically reduce depressive symptoms within 24 hours in a subset of major depressive disorder patients. However, it is unknown whether they exert their actions through shared regulatory mechanisms. To address this question, we performed comparative transcriptomics analyses to identify candidate genes and relevant pathways common to KT and SD.

METHODS

We used the forced swim test, a standardized behavioral approach to measure antidepressant-like activity of KT and SD. We investigated gene expression changes using high-density microarrays and pathway analyses (Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Gene Set Enrichment Analysis) in KT- and SD-treated mice compared with saline-treated control male mice.

RESULTS

We show that KT and SD elicit common transcriptional responses implicating distinct elements of the circadian clock and processes involved in neuronal plasticity. There is an overlap of 64 genes whose expression is common in KT and SD. Specifically, there is downregulation of clock genes including Ciart, Per2, Npas4, Dbp, and Rorb in both KT- and SD-treated mice.

CONCLUSIONS

We demonstrate a potential involvement of the circadian clock in rapid antidepressant responses. These findings could open new research avenues to help design chronopharmacological strategies to treat major depressive disorder.

A genome-wide association study of kynurenic acid in cerebrospinal fluid: implications for psychosis and cognitive impairment in bipolar disorder

Elevated cerebrospinal fluid (CSF) levels of the glia-derived N-methyl-D-aspartic acid receptor antagonist kynurenic acid (KYNA) have consistently been implicated in schizophrenia and bipolar disorder. Here, we conducted a genome-wide association study based on CSF KYNA in bipolar disorder and found support for an association with a common variant within 1p21.3. After replication in an independent cohort, we linked this genetic variant-associated with reduced SNX7 expression-to positive psychotic symptoms and executive function deficits in bipolar disorder. A series of post-mortem brain tissue and in vitro experiments suggested SNX7 downregulation to result in a caspase-8-driven activation of interleukin-1β and a subsequent induction of the brain kynurenine pathway. The current study demonstrates the potential of using biomarkers in genetic studies of psychiatric disorders, and may help to identify novel drug targets in bipolar disorder.

Targets of polyamine dysregulation in major depression and suicide: Activity-dependent feedback, excitability, and neurotransmission

Major depressive disorder (MDD) is a leading cause of disability worldwide characterized by altered neuronal activity in brain regions involved in the control of stress and emotion. Although multiple lines of evidence suggest that altered stress-coping mechanisms underlie the etiology of MDD, the homeostatic control of neuronal excitability in MDD at the molecular level is not well established. In this review, we examine past and current evidence implicating dysregulation of the polyamine system as a central factor in the homeostatic response to stress and the etiology of MDD. We discuss the cellular effects of abnormal metabolism of polyamines in the context of their role in sensing and modulation of neuronal, electrical, and synaptic activity. Finally, we discuss evidence supporting an allostatic model of depression based on a chronic elevation in polyamine levels resulting in self-sustained stress response mechanisms maintained by maladaptive homeostatic mechanisms.

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.

Evidence of Mitochondrial Dysfunction within the Complex Genetic Etiology of Schizophrenia

Genetic evidence has supported the hypothesis that schizophrenia (SZ) is a polygenic disorder caused by the disruption in function of several or many genes. The most common and reproducible cellular phenotype associated with SZ is a reduction in dendritic spines within the neocortex, suggesting alterations in dendritic architecture may cause aberrant cortical circuitry and SZ symptoms. Here, we review evidence supporting a multifactorial model of mitochondrial dysfunction in SZ etiology and discuss how these multiple paths to mitochondrial dysfunction may contribute to dendritic spine loss and/or underdevelopment in some SZ subjects. The pathophysiological role of mitochondrial dysfunction in SZ is based upon genomic analyses of both the mitochondrial genome and nuclear genes involved in mitochondrial function. Previous studies and preliminary data suggest SZ is associated with specific alleles and haplogroups of the mitochondrial genome, and also correlates with a reduction in mitochondrial copy number and an increase in synonymous and nonsynonymous substitutions of mitochondrial DNA. Mitochondrial dysfunction has also been widely implicated in SZ by genome-wide association, exome sequencing, altered gene expression, proteomics, microscopy analyses, and induced pluripotent stem cell studies. Together, these data support the hypothesis that SZ is a polygenic disorder with an enrichment of mitochondrial targets.

Olanzapine Reversed Brain Gene Expression Changes Induced by Phencyclidine Treatment in Non-Human Primates

The NMDA receptor antagonist phencyclidine (PCP) creates schizophrenia-like symptoms in normal controls. The effect of PCP on non-human primate brain gene expression was examined and compared to changes induced by olanzapine treatment. Experimental studies of PCP and antipsychotic drugs have direct relevance to understanding the patho-physiology and treatment of schizophrenia. Genome-wide changes in prefrontal cortex gene expression revealed alterations of 146 transcripts in the PCP treatment group compared to vehicle controls. Dysregulated genes were enriched in identified classes implicated in neurological and genetic disorders, including schizophrenia genes from the Psychiatric Genomics Consortium 108 loci as well as cell death in PCP-treated primates. Canonical pathway analysis revealed a significant overrepresentation of several groups including synaptic long-term potentiation and calcium signaling. Olanzapine coadministered with PCP normalized 34% of the 146 PCP-induced probe set expression changes, and a network of 17 olanzapine-normalized genes was identified enriched in schizophrenia candidate genes containing RGS4, SYN1 and AKT as nodes. The results of this study support the use of PCP administration in non-human primates as a glutamatergic model of schizophrenia and suggest that a large number of PCP-induced expression differences can be reversed by olanzapine. The results of this study may be informative for identification of potential candidates for pharmacogenetics and biomarker research related to the treatment of schizophrenia.

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.

The somatic common deletion in mitochondrial DNA is decreased in schizophrenia

Large deletions in mitochondrial DNA (mtDNA) can occur during or result from oxidative stress leading to a vicious cycle that increases reactive oxygen species (ROS) damage and decreases mitochondrial function, thereby causing further oxidative stress. The objective of this study was to determine if disease specific brain differences of the somatic mtDNA common deletion (4977 bp) could be observed in major depressive disorder (MDD), bipolar disorder (BD), and schizophrenia (SZ) compared to a control group. The accumulation of the mtDNA common deletion was measured using a quantitative assay across 10 brain regions (anterior cingulate cortex, amygdala, caudate nucleus, dorsolateral prefrontal cortex, hippocampus, nucleus accumbens, orbitofrontal cortex, putamen, substantia nigra, and thalamus). The correlation with age of the mtDNA deletion was highly significant across brain regions as previously shown. A significant decrease in the global accumulation of common deletion in subjects with SZ compared to MDD, BD, and controls was observed after correcting for age, pH, PMI, and gender. The decreases in SZ were largest in dopaminergic regions. One potential side effect of antipsychotic drugs on mitochondria is the impairment of mitochondria function, which might explain these findings. The decreased global brain mtDNA common deletion levels suggests that mitochondrial function is impaired and might be part of an overall mitochondria dysfunction signature in subjects with schizophrenia.

Peripheral biomarkers revisited: integrative profiling of peripheral samples for psychiatric research

Peripheral samples, such as blood and skin, have been used for decades in psychiatric research as surrogates for central nervous system samples. Although the validity of the data obtained from peripheral samples has been questioned and other state-of-the-art techniques, such as human brain imaging, genomics, and induced pluripotent stem cells, seem to reduce the value of peripheral cells, accumulating evidence has suggested that revisiting peripheral samples is worthwhile. Here, we re-evaluate the utility of peripheral samples and argue that establishing an understanding of the common signaling and biological processes in the brain and peripheral samples is required for the validity of such models. First, we present an overview of the available types of peripheral cells and describe their advantages and disadvantages. We then briefly summarize the main achievements of omics studies, including epigenome, transcriptome, proteome, and metabolome analyses, as well as the main findings of functional cellular assays, the results of which imply that alterations in neurotransmission, metabolism, the cell cycle, and the immune system may be partially responsible for the pathophysiology of major psychiatric disorders such as schizophrenia. Finally, we discuss the future utility of peripheral samples for the development of biomarkers and tailor-made therapies, such as multimodal assays that are used as a battery of disease and trait pathways and that might be potent and complimentary tools for use in psychiatric research.

Identification of pathways for bipolar disorder: a meta-analysis

IMPORTANCE

Genome-wide investigations provide systematic information regarding the neurobiology of psychiatric disorders.

OBJECTIVE

To identify biological pathways that contribute to risk for bipolar disorder (BP) using genes with consistent evidence for association in multiple genome-wide association studies (GWAS).

DATA SOURCES

Four independent data sets with individual genome-wide data available in July 2011 along with all data sets contributed to the Psychiatric Genomics Consortium Bipolar Group by May 2012. A prior meta-analysis was used as a source for brain gene expression data.

STUDY SELECTION

The 4 published GWAS were included in the initial sample. All independent BP data sets providing genome-wide data in the Psychiatric Genomics Consortium were included as a replication sample.

DATA EXTRACTION AND SYNTHESIS

We identified 966 genes that contained 2 or more variants associated with BP at P < .05 in 3 of 4 GWAS data sets (n = 12,127 [5253 cases, 6874 controls]). Simulations using 10,000 replicates of these data sets corrected for gene size and allowed the calculation of an empirical P value for each gene; empirically significant genes were entered into a pathway analysis. Each of these pathways was then tested in the replication sample (n = 8396 [3507 cases, 4889 controls]) using gene set enrichment analysis for single-nucleotide polymorphisms. The 226 genes were also compared with results from a meta-analysis of gene expression in the dorsolateral prefrontal cortex.

MAIN OUTCOMES AND MEASURES

Empirically significant genes and biological pathways. RESULTS Among 966 genes, 226 were empirically significant (P < .05). Seventeen pathways were overrepresented in analyses of the initial data set. Six of the 17 pathways were associated with BP in both the initial and replication samples: corticotropin-releasing hormone signaling, cardiac β-adrenergic signaling, phospholipase C signaling, glutamate receptor signaling, endothelin 1 signaling, and cardiac hypertrophy signaling. Among the 226 genes, 9 differed in expression in the dorsolateral prefrontal cortex in patients with BP: CACNA1C, DTNA, FOXP1, GNG2, ITPR2, LSAMP, NPAS3, NCOA2, and NTRK3.

CONCLUSIONS AND RELEVANCE

Pathways involved in the genetic predisposition to BP include hormonal regulation, calcium channels, second messenger systems, and glutamate signaling. Gene expression studies implicate neuronal development pathways as well. These results tend to reinforce specific hypotheses regarding BP neurobiology and may provide clues for new approaches to treatment and prevention.

Next generation sequence analysis and computational genomics using graphical pipeline workflows

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.

Schizophrenia miR-137 locus risk genotype is associated with dorsolateral prefrontal cortex hyperactivation

BACKGROUND

miR-137 dysregulation has been implicated in the etiology of schizophrenia, but its functional role remains to be determined.

METHODS

Functional magnetic resonance imaging scans were acquired on 48 schizophrenia patients and 63 healthy volunteers (total sample size N = 111 subjects), with similar mean age and sex distribution, while subjects performed a Sternberg Item Response Paradigm with memory loads of one, three, and five numbers. Dorsolateral prefrontal cortex (DLPFC) retrieval activation for the working memory load of three numbers, for which hyperactivation had been shown in schizophrenia patients compared with control subjects, was extracted. The genome-wide association study confirmed schizophrenia risk single nucleotide polymorphism rs1625579 (miR-137 locus) was genotyped (schizophrenia: GG n = 0, GT n = 9, TT n = 39; healthy volunteers: GG = 2, GT n = 15, and TT n = 46). Fisher's exact test examined the effect of diagnosis on rs1625579 allele frequency distribution (p = nonsignificant). Mixed model regression analyses examined the effects of diagnosis and genotype on working memory performance measures and DLPFC activation.

RESULTS

Patients showed significantly higher left DLPFC retrieval activation on working memory load 3, lower working memory performance, and longer response times compared with controls. There was no effect of genotype on working memory performance or response times in either group. However, individuals with the rs1625579 TT genotype had significantly higher left DLPFC activation than those with the GG/GT genotypes.

CONCLUSIONS

Our study suggests that the rs1625579 TT (miR-137 locus) schizophrenia risk genotype is associated with the schizophrenia risk phenotype DLPFC hyperactivation commonly considered a measure of brain inefficiency.

The KMO allele encoding Arg452 is associated with psychotic features in bipolar disorder type 1, and with increased CSF KYNA level and reduced KMO expression

The kynurenine pathway metabolite kynurenic acid (KYNA), modulating glutamatergic and cholinergic neurotransmission, is increased in cerebrospinal fluid (CSF) of patients with schizophrenia or bipolar disorder type 1 with psychotic features. KYNA production is critically dependent on kynurenine 3-monooxygenase (KMO). KMO mRNA levels and activity in prefrontal cortex (PFC) are reduced in schizophrenia. We hypothesized that KMO expression in PFC would be reduced in bipolar disorder with psychotic features and that a functional genetic variant of KMO would associate with this disease, CSF KYNA level and KMO expression. KMO mRNA levels were reduced in PFC of bipolar disorder patients with lifetime psychotic features (P=0.005, n=19) or schizophrenia (P=0.02, n=36) compared with nonpsychotic patients and controls. KMO genetic association to psychotic features in bipolar disorder type 1 was studied in 493 patients and 1044 controls from Sweden. The KMO Arg(452) allele was associated with psychotic features during manic episodes (P=0.003). KMO Arg(452) was studied for association to CSF KYNA levels in an independent sample of 55 Swedish patients, and to KMO expression in 717 lymphoblastoid cell lines and 138 hippocampal biopsies. KMO Arg(452) associated with increased levels of CSF KYNA (P=0.03) and reduced lymphoblastoid and hippocampal KMO expression (P≤0.05). Thus, findings from five independent cohorts suggest that genetic variation in KMO influences the risk for psychotic features in mania of bipolar disorder patients. This provides a possible mechanism for the previous findings of elevated CSF KYNA levels in those bipolar patients with lifetime psychotic features and positive association between KYNA levels and number of manic episodes.

Molecular and bioenergetic differences between cells with African versus European inherited mitochondrial DNA haplogroups: implications for population susceptibility to diseases

The geographic origins of populations can be identified by their maternally inherited mitochondrial DNA (mtDNA) haplogroups. This study compared human cybrids (cytoplasmic hybrids), which are cell lines with identical nuclei but mitochondria from different individuals with mtDNA from either the H haplogroup or L haplogroup backgrounds. The most common European haplogroup is H while individuals of maternal African origin are of the L haplogroup. Despite lower mtDNA copy numbers, L cybrids had higher expression levels for nine mtDNA-encoded respiratory complex genes, decreased ATP (adenosine triphosphate) turnover rates and lower levels of reactive oxygen species production, parameters which are consistent with more efficient oxidative phosphorylation. Surprisingly, GeneChip arrays showed that the L and H cybrids had major differences in expression of genes of the canonical complement system (5 genes), dermatan/chondroitin sulfate biosynthesis (5 genes) and CCR3 (chemokine, CC motif, receptor 3) signaling (9 genes). Quantitative nuclear gene expression studies confirmed that L cybrids had (a) lower expression levels of complement pathway and innate immunity genes and (b) increased levels of inflammation-related signaling genes, which are critical in human diseases. Our data support the hypothesis that mtDNA haplogroups representing populations from different geographic origins may play a role in differential susceptibilities to diseases.

A splice donor mutation in NAA10 results in the dysregulation of the retinoic acid signalling pathway and causes Lenz microphthalmia syndrome

INTRODUCTION

Lenz microphthalmia syndrome (LMS) is a genetically heterogeneous X-linked disorder characterised by microphthalmia/anophthalmia, skeletal abnormalities, genitourinary malformations, and anomalies of the digits, ears, and teeth. Intellectual disability and seizure disorders are seen in about 60% of affected males. To date, no gene has been identified for LMS in the microphthalmia syndrome 1 locus (MCOPS1). In this study, we aim to find the disease-causing gene for this condition.

METHODS AND RESULTS

Using exome sequencing in a family with three affected brothers, we identified a mutation in the intron 7 splice donor site (c.471+2T→A) of the N-acetyltransferase NAA10 gene. NAA10 has been previously shown to be mutated in patients with Ogden syndrome, which is clinically distinct from LMS. Linkage studies for this family mapped the disease locus to Xq27-Xq28, which was consistent with the locus of NAA10. The mutation co-segregated with the phenotype and cDNA analysis showed aberrant transcripts. Patient fibroblasts lacked expression of full length NAA10 protein and displayed cell proliferation defects. Expression array studies showed significant dysregulation of genes associated with genetic forms of anophthalmia such as BMP4, STRA6, and downstream targets of BCOR and the canonical WNT pathway. In particular, STRA6 is a retinol binding protein receptor that mediates cellular uptake of retinol/vitamin A and plays a major role in regulating the retinoic acid signalling pathway. A retinol uptake assay showed that retinol uptake was decreased in patient cells.

CONCLUSIONS

We conclude that the NAA10 mutation is the cause of LMS in this family, likely through the dysregulation of the retinoic acid signalling pathway.

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.

Analysis of miR-137 expression and rs1625579 in dorsolateral prefrontal cortex

MicroRNAs (miRNAs) are small non-coding RNAs that act as potent regulators of gene expression. A recent GWAS reported the rs1625579 SNP, located downstream of miR-137, as the strongest new association with schizophrenia [Ripke S, Sanders AR, Kendler KS, Levinson DF, Sklar P, Holmans PA, et al. Genome-wide association study identifies five new schizophrenia loci. Nat Genet 2011;43:969-76.]. Prior to this GWAS finding, a schizophrenia imaging-genetic study found miR-137 target genes significantly enriched for association with activation in the dorsolateral prefrontal cortex (DLPFC) [Potkin SG, Macciardi F, Guffanti G, Fallon JH, Wang Q, Turner JA, et al. Identifying gene regulatory networks in schizophrenia. Neuroimage 2010;53:839-47.]. We investigated the expression levels of miR-137 and three candidate target genes (ZNF804A, CACNA1C, TCF4) in the DLPFC of postmortem brain tissue from 2 independent cohorts: (1) 26 subjects (10 control (CTR), 7 schizophrenia (SZ), 9 bipolar disorder (BD)) collected at the UCI brain bank; and (2) 99 subjects (33 CTR, 35 SZ, 31 BD) obtained from the Stanley Medical Research Institute (SMRI). MiR-137 expression in the DLPFC did not differ between diagnoses. We also explored the relationship between rs1625579 genotypes and miR-137 expression. Significantly lower miR-137 expression levels were observed in the homozygous TT subjects compared to TG and GG subjects in the control group (30% decrease, p-value = 0.03). Moreover, reduced miR-137 levels in TT subjects corresponded to increased levels of the miR-137 target gene TCF4. The miR-137 expression pattern in 9 brain regions was significant for regional effect (ANOVA p-value = 1.83E-12), with amygdala and hippocampus having the highest miR-137 expression level. In conclusion, decreased miR-137 expression is associated with the SZ risk allele of rs1625579, and potential regulation of TCF4, another SZ candidate gene. This study offers additional support for involvement of miR-137 and downstream targets as mechanisms of risk for psychiatric disorders.

Conserved chromosome 2q31 conformations are associated with transcriptional regulation of GAD1 GABA synthesis enzyme and altered in prefrontal cortex of subjects with schizophrenia

Little is known about chromosomal loopings involving proximal promoter and distal enhancer elements regulating GABAergic gene expression, including changes in schizophrenia and other psychiatric conditions linked to altered inhibition. Here, we map in human chromosome 2q31 the 3D configuration of 200 kb of linear sequence encompassing the GAD1 GABA synthesis enzyme gene locus, and we describe a loop formation involving the GAD1 transcription start site and intergenic noncoding DNA elements facilitating reporter gene expression. The GAD1-TSS(-50kbLoop) was enriched with nucleosomes epigenetically decorated with the transcriptional mark, histone H3 trimethylated at lysine 4, and was weak or absent in skin fibroblasts and pluripotent stem cells compared with neuronal cultures differentiated from them. In the prefrontal cortex of subjects with schizophrenia, GAD1-TSS(-50kbLoop) was decreased compared with controls, in conjunction with downregulated GAD1 expression. We generated transgenic mice expressing Gad2 promoter-driven green fluorescent protein-conjugated histone H2B and confirmed that Gad1-TSS(-55kbLoop), the murine homolog to GAD1-TSS(-50kbLoop), is a chromosomal conformation specific for GABAergic neurons. In primary neuronal culture, Gad1-TSS(-55kbLoop) and Gad1 expression became upregulated when neuronal activity was increased. We conclude that 3D genome architectures, including chromosomal loopings for promoter-enhancer interactions involved in the regulation of GABAergic gene expression, are conserved between the rodent and primate brain, and subject to developmental and activity-dependent regulation, and disordered in some cases with schizophrenia. More broadly, the findings presented here draw a connection between noncoding DNA, spatial genome architecture, and neuronal plasticity in development and disease.

An anatomically comprehensive atlas of the adult human brain transcriptome

Neuroanatomically precise, genome-wide maps of transcript distributions are critical resources to complement genomic sequence data and to correlate functional and genetic brain architecture. Here we describe the generation and analysis of a transcriptional atlas of the adult human brain, comprising extensive histological analysis and comprehensive microarray profiling of ∼900 neuroanatomically precise subdivisions in two individuals. Transcriptional regulation varies enormously by anatomical location, with different regions and their constituent cell types displaying robust molecular signatures that are highly conserved between individuals. Analysis of differential gene expression and gene co-expression relationships demonstrates that brain-wide variation strongly reflects the distributions of major cell classes such as neurons, oligodendrocytes, astrocytes and microglia. Local neighbourhood relationships between fine anatomical subdivisions are associated with discrete neuronal subtypes and genes involved with synaptic transmission. The neocortex displays a relatively homogeneous transcriptional pattern, but with distinct features associated selectively with primary sensorimotor cortices and with enriched frontal lobe expression. Notably, the spatial topography of the neocortex is strongly reflected in its molecular topography-the closer two cortical regions, the more similar their transcriptomes. This freely accessible online data resource forms a high-resolution transcriptional baseline for neurogenetic studies of normal and abnormal human brain function.

Mitochondrial mutations and polymorphisms in psychiatric disorders

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.