Genome-wide association studies have identified 108 schizophrenia risk loci, but biological mechanisms for individual loci are largely unknown. Using developmental, genetic and illness-based RNA ...sequencing expression analysis in human brain, we characterized the human brain transcriptome around these loci and found enrichment for developmentally regulated genes with novel examples of shifting isoform usage across pre- and postnatal life. We found widespread expression quantitative trait loci (eQTLs), including many with transcript specificity and previously unannotated sequence that were independently replicated. We leveraged this general eQTL database to show that 48.1% of risk variants for schizophrenia associate with nearby expression. We lastly found 237 genes significantly differentially expressed between patients and controls, which replicated in an independent dataset, implicated synaptic processes, and were strongly regulated in early development. These findings together offer genetics- and diagnosis-related targets for better modeling of schizophrenia risk. This resource is publicly available at http://eqtl.brainseq.org/phase1 .
The hippocampus formation, although prominently implicated in schizophrenia pathogenesis, has been overlooked in large-scale genomics efforts in the schizophrenic brain. We performed RNA-seq in ...hippocampi and dorsolateral prefrontal cortices (DLPFCs) from 551 individuals (286 with schizophrenia). We identified substantial regional differences in gene expression and found widespread developmental differences that were independent of cellular composition. We identified 48 and 245 differentially expressed genes (DEGs) associated with schizophrenia within the hippocampus and DLPFC, with little overlap between the brain regions. 124 of 163 (76.6%) of schizophrenia GWAS risk loci contained eQTLs in any region. Transcriptome-wide association studies in each region identified many novel schizophrenia risk features that were brain region-specific. Last, we identified potential molecular correlates of in vivo evidence of altered prefrontal-hippocampal functional coherence in schizophrenia. These results underscore the complexity and regional heterogeneity of the transcriptional correlates of schizophrenia and offer new insights into potentially causative biology.
•Dorsolateral prefrontal cortex and hippocampus gene expression across development•Novel region-specific schizophrenia genetic risk features•Decreased regional functional coherence in schizophrenia•Public brain gene expression and eQTL resource at http://eqtl.brainseq.org/phase2
Collado-Torres et al. describe the BrainSeq Phase II gene expression resource encompassing two brain regions from 551 genotyped individuals spanning the entire human lifespan (286 with schizophrenia). This resource can answer region-specific questions about development and schizophrenia and its genetic risk.
Genome-wide association studies (GWASs) have reported many single nucleotide polymorphisms (SNPs) associated with psychiatric disorders, but knowledge is lacking regarding molecular mechanisms. Here ...we show that risk alleles spanning multiple genes across the 10q24.32 schizophrenia-related locus are associated in the human brain selectively with an increase in the expression of both BLOC-1 related complex subunit 7 (BORCS7) and a previously uncharacterized, human-specific arsenite methyltransferase (AS3MT) isoform (AS3MT(d2d3)), which lacks arsenite methyltransferase activity and is more abundant in individuals with schizophrenia than in controls. Conditional-expression analysis suggests that BORCS7 and AS3MT(d2d3) signals are largely independent. GWAS risk SNPs across this region are linked with a variable number tandem repeat (VNTR) polymorphism in the first exon of AS3MT that is associated with the expression of AS3MT(d2d3) in samples from both Caucasians and African Americans. The VNTR genotype predicts promoter activity in luciferase assays, as well as DNA methylation within the AS3MT gene. Both AS3MT(d2d3) and BORCS7 are expressed in adult human neurons and astrocytes, and they are upregulated during human stem cell differentiation toward neuronal fates. Our results provide a molecular explanation for the prominent 10q24.32 locus association, including a novel and evolutionarily recent protein that is involved in early brain development and confers risk for psychiatric illness.
GABA signaling molecules are critical for both human brain development and the pathophysiology of schizophrenia. We examined the expression of transcripts derived from three genes related to GABA ...signaling GAD1 (GAD67 and GAD25), SLC12A2 (NKCC1), and SLC12A5 (KCC2) in the prefrontal cortex (PFC) and hippocampal formation of a large cohort of nonpsychiatric control human brains (n = 240) across the lifespan (from fetal week 14 to 80 years) and in patients with schizophrenia (n = 30-31), using quantitative RT-PCR. We also examined whether a schizophrenia risk-associated promoter SNP in GAD1 (rs3749034) is related to expression of these transcripts. Our studies revealed that development and maturation of both the PFC and hippocampal formation are characterized by progressive switches in expression from GAD25 to GAD67 and from NKCC1 to KCC2. Previous studies have demonstrated that the former leads to GABA synthesis, and the latter leads to switching from excitatory to inhibitory neurotransmission. In the hippocampal formation, GAD25/GAD67 and NKCC1/KCC2 ratios are increased in patients with schizophrenia, reflecting a potentially immature GABA physiology. Remarkably, GAD25/GAD67 and NKCC1/KCC2 expression ratios are associated with rs3749034 genotype, with risk alleles again predicting a relatively less mature pattern. These findings suggest that abnormalities in GABA signaling critical to brain development contribute to genetic risk for schizophrenia.
See Calafato and Bramon (doi:10.1093/brain/awy345) for a scientific commentary on this article.
Cognitive deficits accompany schizophrenia and are associated with genetic risk but the direction of ...causality is unclear. Using causal modelling, Toulopoulou et al. show that cognitive deficits partially mediate the effects of cumulative genetic risk for schizophrenia. Impaired cognition is an intermediate phenotype on the causal path to the disorder with implications for diagnosis and intervention.
Abstract
Cognitive deficit is thought to represent, at least in part, genetic mechanisms of risk for schizophrenia, with recent evidence from statistical modelling of twin data suggesting direct causality from the former to the latter. However, earlier evidence was based on inferences from twin not molecular genetic data and it is unclear how much genetic influence 'passes through' cognition on the way to diagnosis. Thus, we included direct measurements of genetic risk (e.g. schizophrenia polygenic risk scores) in causation models to assess the extent to which cognitive deficit mediates some of the effect of polygenic risk scores on the disorder. Causal models of family data tested relationships among key variables and allowed parsing of genetic variance components. Polygenic risk scores were calculated from summary statistics from the current largest genome-wide association study of schizophrenia and were represented as a latent trait. Cognition was also modelled as a latent trait. Participants were 1313 members of 1078 families: 416 patients with schizophrenia, 290 unaffected siblings, and 607 controls. Modelling supported earlier findings that cognitive deficit has a putatively causal role in schizophrenia. In total, polygenic risk score explained 8.07% confidence interval (CI) 5.45-10.74% of schizophrenia risk in our sample. Of this, more than a third (2.71%, CI 2.41-3.85%) of the polygenic risk score influence was mediated through cognition paths, exceeding the direct influence of polygenic risk score on schizophrenia risk (1.43%, CI 0.46-3.08%). The remainder of the polygenic risk score influence (3.93%, CI 2.37-4.48%) reflected reciprocal causation between schizophrenia liability and cognition (e.g. mutual influences in a cyclical manner). Analysis of genetic variance components of schizophrenia liability indicated that 26.87% (CI 21.45-32.57%) was associated with cognition-related pathways not captured by polygenic risk score. The remaining variance in schizophrenia was through pathways other than cognition-related and polygenic risk score. Although our results are based on inference through statistical modelling and do not provide an absolute proof of causality, we find that cognition pathways mediate a significant part of the influence of cumulative genetic risk on schizophrenia. We estimate from our model that 33.51% (CI 27.34-43.82%) of overall genetic risk is mediated through influences on cognition, but this requires further studies and analyses as the genetics of schizophrenia becomes better characterized.
Schizophrenia polygenic risk is plausibly manifested by complex transcriptional dysregulation in the brain, involving networks of co-expressed and functionally related genes. The main purpose of this ...study was to identify and prioritize co-expressed gene sets in a hierarchical manner, based on the strength of the relationships with clinical diagnosis and with polygenic risk for schizophrenia. Weighted Gene Co-expression Network Analysis (WGCNA) was applied to RNA-quality-adjusted DLPFC RNA-Seq data from the LIBD Postmortem Human Brain Repository (90 controls, 74 schizophrenia cases; all Caucasians) to construct co-expression networks and detect "modules" of co-expressed genes. After multiple internal and external validation procedures, modules of selected interest were tested for enrichment in biological ontologies, for association with schizophrenia polygenic risk scores (PRSs) and with diagnosis, and also for enrichment in genes within the significant GWAS loci reported by the Psychiatric Genomic Consortium (PGC2). The association between schizophrenia genetic signals and modules of co-expression converged on one module showing not only a significant association with both diagnosis and PRS but also significant overlap with 36 PGC2 loci genes, deemed the strongest candidates for drug targets. This module contained many genes involved in synaptic signaling and neuroplasticity. Fifty-three PGC2 genes were in modules associated only with diagnosis and 59 in modules unrelated to diagnosis or PRS. Our study highlights complex relationships between gene co-expression networks in the brain and clinical state and polygenic risk for SCZ and provides a strategy for using this information in selecting and prioritizing potentially targetable gene sets for therapeutic drug development.
Genome-wide association studies (GWAS) have identified many genomic loci associated with risk for schizophrenia, but unambiguous identification of the relationship between disease-associated variants ...and specific genes, and in particular their effect on risk conferring transcripts, has proven difficult. To better understand the specific molecular mechanism(s) at the schizophrenia locus in 11q25, we undertook cis expression quantitative trait loci (cis-eQTL) mapping for this 2 megabase genomic region using postmortem human brain samples. To comprehensively assess the effects of genetic risk upon local expression, we evaluated multiple transcript features: genes, exons, and exon-exon junctions in multiple brain regions-dorsolateral prefrontal cortex (DLPFC), hippocampus, and caudate. Genetic risk variants strongly associated with expression of SNX19 transcript features that tag multiple rare classes of SNX19 transcripts, whereas they only weakly affected expression of an exon-exon junction that tags the majority of abundant transcripts. The most prominent class of SNX19 risk-associated transcripts is predicted to be overexpressed, defined by an exon-exon splice junction between exons 8 and 10 (junc8.10) and that is predicted to encode proteins that lack the characteristic nexin C terminal domain. Risk alleles were also associated with either increased or decreased expression of multiple additional classes of transcripts. With RACE, molecular cloning, and long read sequencing, we found a number of novel SNX19 transcripts that further define the set of potential etiological transcripts. We explored epigenetic regulation of SNX19 expression and found that DNA methylation at CpG sites near the primary transcription start site and within exon 2 partially mediate the effects of risk variants on risk-associated expression. ATAC sequencing revealed that some of the most strongly risk-associated SNPs are located within a region of open chromatin, suggesting a nearby regulatory element is involved. These findings indicate a potentially complex molecular etiology, in which risk alleles for schizophrenia generate epigenetic alterations and dysregulation of multiple classes of SNX19 transcripts.
Genetic variation in neuregulin 1 (NRG1) is associated with schizophrenia. The disease-associated SNPs are noncoding, and their functional implications remain unknown. We hypothesized that ...differential expression of the NRG1 gene explains its association to the disease. We examined four of the disease-associated SNPs that make up the original risk haplotype in the 5' upstream region of the gene for their effects on mRNA abundance of NRG1 types I-IV in human postmortem hippocampus. Diagnostic comparisons revealed a 34% increase in type I mRNA in schizophrenia and an interaction of diagnosis and genotype (SNP8NRG221132) on this transcript. Of potentially greater interest, a single SNP within the risk haplotype (SNP8NRG243177) and a 22-kb block of this core haplotype are associated with mRNA expression for the novel type IV isoform in patients and controls. Bioinformatic promoter analyses indicate that both SNPs lead to a gain/loss of putative binding sites for three transcription factors, serum response factor, myelin transcription factor-1, and High Mobility Group Box Protein-1. These data implicate variation in isoform expression as a molecular mechanism for the genetic association of NRG1 with schizophrenia.
A variation in the BDNF gene (val66met) affects the function of BDNF in neurons, predicts variation in human memory, and is associated with several neurological and psychiatric disorders. Here, we ...show that, in magnetic resonance imaging scans of a large sample of normal individuals, this polymorphism affects the anatomy of the hippocampus and prefrontal cortex, identifying a genetic mechanism of variation in brain morphology related to learning and memory.
PDF
