Genetic association of FMRP targets with psychiatric disorders.

Genes encoding the mRNA targets of fragile X mental retardation protein (FMRP) are enriched for genetic association with psychiatric disorders. However, many FMRP targets possess functions that are themselves genetically associated with psychiatric disorders, including synaptic transmission and plasticity, making it unclear whether the genetic risk is truly related to binding by FMRP or is alternatively mediated by the sampling of genes better characterised by another trait or functional annotation. Using published common variant, rare coding variant and copy number variant data, we examined the relationship between FMRP binding and genetic association with schizophrenia, major depressive disorder and bipolar disorder. High-confidence targets of FMRP, derived from studies of multiple tissue types, were enriched for common schizophrenia risk alleles, as well as rare loss-of-function and de novo nonsynonymous variants in schizophrenia cases. Similarly, through common variation, FMRP targets were associated with major depressive disorder, and we present novel evidence of association with bipolar disorder. These relationships could not be explained by other functional annotations known to be associated with psychiatric disorders, including those related to synaptic structure and function. This study reinforces the evidence that targeting by FMRP captures a subpopulation of genes enriched for genetic association with a range of psychiatric disorders.

De novo mutations in schizophrenia implicate synaptic networks.

Inherited alleles account for most of the genetic risk for schizophrenia. However, new (de novo) mutations, in the form of large chromosomal copy number changes, occur in a small fraction of cases and disproportionally disrupt genes encoding postsynaptic proteins. Here we show that small de novo mutations, affecting one or a few nucleotides, are overrepresented among glutamatergic postsynaptic proteins comprising activity-regulated cytoskeleton-associated protein (ARC) and N-methyl-d-aspartate receptor (NMDAR) complexes. Mutations are additionally enriched in proteins that interact with these complexes to modulate synaptic strength, namely proteins regulating actin filament dynamics and those whose messenger RNAs are targets of fragile X mental retardation protein (FMRP). Genes affected by mutations in schizophrenia overlap those mutated in autism and intellectual disability, as do mutation-enriched synaptic pathways. Aligning our findings with a parallel case-control study, we demonstrate reproducible insights into aetiological mechanisms for schizophrenia and reveal pathophysiology shared with other neurodevelopmental disorders.

The synapse in schizophrenia.

It has been several decades since synaptic dysfunction was first suggested to play a role in schizophrenia, but only in the last few years has convincing evidence been obtained as progress has been made in elucidating the genetic underpinnings of the disorder. In the intervening years much has been learned concerning the complex macromolecular structure of the synapse itself, and genetic studies are now beginning to draw upon these advances. Here we outline our current understanding of the genetic architecture of schizophrenia and examine the evidence for synaptic involvement. A strong case can now be made that disruption of glutamatergic signalling pathways regulating synaptic plasticity contributes to the aetiology of schizophrenia.

Evolutionary expansion and anatomical specialization of synapse proteome complexity.

Understanding the origins and evolution of synapses may provide insight into species diversity and the organization of the brain. Using comparative proteomics and genomics, we examined the evolution of the postsynaptic density (PSD) and membrane-associated guanylate kinase (MAGUK)-associated signaling complexes (MASCs) that underlie learning and memory. PSD and MASC orthologs found in yeast carry out basic cellular functions to regulate protein synthesis and structural plasticity. We observed marked changes in signaling complexity at the yeast-metazoan and invertebrate-vertebrate boundaries, with an expansion of key synaptic components, notably receptors, adhesion/cytoskeletal proteins and scaffold proteins. A proteomic comparison of Drosophila and mouse MASCs revealed species-specific adaptation with greater signaling complexity in mouse. Although synaptic components were conserved amongst diverse vertebrate species, mapping mRNA and protein expression in the mouse brain showed that vertebrate-specific components preferentially contributed to differences between brain regions. We propose that the evolution of synapse complexity around a core proto-synapse has contributed to invertebrate-vertebrate differences and to brain specialization.

Transcriptional programs regulating neuronal differentiation are disrupted in DLG2 knockout human embryonic stem cells and enriched for schizophrenia and related disorders risk variants.

Coordinated programs of gene expression drive brain development. It is unclear which transcriptional programs, in which cell-types, are affected in neuropsychiatric disorders such as schizophrenia. Here we integrate human genetics with transcriptomic data from differentiation of human embryonic stem cells into cortical excitatory neurons. We identify transcriptional programs expressed during early neurogenesis in vitro and in human foetal cortex that are down-regulated in DLG2-/- lines. Down-regulation impacted neuronal differentiation and maturation, impairing migration, morphology and action potential generation. Genetic variation in these programs is associated with neuropsychiatric disorders and cognitive function, with associated variants predominantly concentrated in loss-of-function intolerant genes. Neurogenic programs also overlap schizophrenia GWAS enrichment previously identified in mature excitatory neurons, suggesting that pathways active during prenatal cortical development may also be associated with mature neuronal dysfunction. Our data from human embryonic stem cells, when combined with analysis of available foetal cortical gene expression data, de novo rare variants and GWAS statistics for neuropsychiatric disorders and cognition, reveal a convergence on transcriptional programs regulating excitatory cortical neurogenesis.

Polygenic risk for schizophrenia and subcortical brain anatomy in the UK Biobank cohort.

Research has shown differences in subcortical brain volumes between participants with schizophrenia and healthy controls. However, none of these differences have been found to associate with schizophrenia polygenic risk. Here, in a large sample (n = 14,701) of unaffected participants from the UK Biobank, we test whether schizophrenia polygenic risk scores (PRS) limited to specific gene-sets predict subcortical brain volumes. We compare associations with schizophrenia PRS at the whole genome level ('genomic', including all SNPs associated with the disorder at a p-value threshold < 0.05) with 'genic' PRS (based on SNPs in the vicinity of known genes), 'intergenic' PRS (based on the remaining SNPs), and genic PRS limited to SNPs within 7 gene-sets previously found to be enriched for genetic association with schizophrenia ('abnormal behaviour,' 'abnormal long-term potentiation,' 'abnormal nervous system electrophysiology,' 'FMRP targets,' '5HT2C channels,' 'CaV2 channels' and 'loss-of-function intolerant genes'). We observe a negative association between the 'abnormal behaviour' gene-set PRS and volume of the right thalamus that survived correction for multiple testing (ß = -0.031, pFDR = 0.005) and was robust to different schizophrenia PRS p-value thresholds. In contrast, the only association with genomic PRS surviving correction for multiple testing was for right pallidum, which was observed using a schizophrenia PRS p-value threshold < 0.01 (ß = -0.032, p = 0.0003, pFDR = 0.02), but not when using other PRS P-value thresholds. We conclude that schizophrenia PRS limited to functional gene sets may provide a better means of capturing differences in subcortical brain volume than whole genome PRS approaches.

De novo mutations identified by exome sequencing implicate rare missense variants in SLC6A1 in schizophrenia.

Schizophrenia is a highly polygenic disorder with important contributions from both common and rare risk alleles. We analyzed exome sequencing data for de novo variants (DNVs) in a new sample of 613 schizophrenia trios and combined this with published data to give a total of 3,444 trios. In this new data, loss-of-function (LoF) DNVs were significantly enriched among 3,471 LoF-intolerant genes, which supports previous findings. In the full dataset, genes associated with neurodevelopmental disorders (n = 159) were significantly enriched for LoF DNVs. Within these neurodevelopmental disorder genes, SLC6A1, which encodes a γ-aminobutyric acid transporter, was associated with missense-damaging DNVs. In 1,122 trios for which genome-wide common variant data were available, schizophrenia and bipolar disorder polygenic risk were significantly overtransmitted to probands. Probands carrying LoF or deletion DNVs in LoF-intolerant or neurodevelopmental disorder genes had significantly less overtransmission of schizophrenia polygenic risk than did non-carriers, which provides a second robust line of evidence that these DNVs increase liability to schizophrenia.

Novel Insight Into the Etiology of Autism Spectrum Disorder Gained by Integrating Expression Data With Genome-wide Association Statistics.

BACKGROUND: A recent genome-wide association study (GWAS) of autism spectrum disorder (ASD) (ncases = 18,381, ncontrols = 27,969) has provided novel opportunities for investigating the etiology of ASD. Here, we integrate the ASD GWAS summary statistics with summary-level gene expression data to infer differential gene expression in ASD, an approach called transcriptome-wide association study (TWAS). METHODS: Using FUSION software, ASD GWAS summary statistics were integrated with predictors of gene expression from 16 human datasets, including adult and fetal brains. A novel adaptation of established statistical methods was then used to test for enrichment within candidate pathways and specific tissues and at different stages of brain development. The proportion of ASD heritability explained by predicted expression of genes in the TWAS was estimated using stratified linkage disequilibrium score regression. RESULTS: This study identified 14 genes as significantly differentially expressed in ASD, 13 of which were outside of known genome-wide significant loci (±500 kb). XRN2, a gene proximal to an ASD GWAS locus, was inferred to be significantly upregulated in ASD, providing insight into the functional consequence of this associated locus. One novel transcriptome-wide significant association from this study is the downregulation of PDIA6, which showed minimal evidence of association in the GWAS, and in gene-based analysis using MAGMA. Predicted gene expression in this study accounted for 13.0% of the total ASD single nucleotide polymorphism heritability. CONCLUSIONS: This study has implicated several genes as significantly up/downregulated in ASD, providing novel and useful information for subsequent functional studies. This study also explores the utility of TWAS-based enrichment analysis and compares TWAS results with a functionally agnostic approach.

Pharmacogenomic Variants and Drug Interactions Identified Through the Genetic Analysis of Clozapine Metabolism.

OBJECTIVE: Clozapine is the only effective medication for treatment-resistant schizophrenia, but its worldwide use is still limited because of its complex titration protocols. While the discovery of pharmacogenomic variants of clozapine metabolism may improve clinical management, no robust findings have yet been reported. This study is the first to adopt the framework of genome-wide association studies (GWASs) to discover genetic markers of clozapine plasma concentrations in a large sample of patients with treatment-resistant schizophrenia. METHODS: The authors used mixed-model regression to combine data from multiple assays of clozapine metabolite plasma concentrations from a clozapine monitoring service and carried out a genome-wide analysis of clozapine, norclozapine, and their ratio on 10,353 assays from 2,989 individuals. These analyses were adjusted for demographic factors known to influence clozapine metabolism, although it was not possible to adjust for all potential mediators given the available data. GWAS results were used to pinpoint specific enzymes and metabolic pathways and compounds that might interact with clozapine pharmacokinetics. RESULTS: The authors identified four distinct genome-wide significant loci that harbor common variants affecting the metabolism of clozapine or its metabolites. Detailed examination pointed to coding and regulatory variants at several CYP* and UGT* genes as well as corroborative evidence for interactions between the metabolism of clozapine, coffee, and tobacco. Individual effects of single single-nucleotide polymorphisms (SNPs) fine-mapped from these loci were large, such as the minor allele of rs2472297, which was associated with a reduction in clozapine concentrations roughly equivalent to a decrease of 50 mg/day in clozapine dosage. On their own, these single SNPs explained from 1.15% to 9.48% of the variance in the plasma concentration data. CONCLUSIONS: Common genetic variants with large effects on clozapine metabolism exist and can be found via genome-wide approaches. Their identification opens the way for clinical studies assessing the use of pharmacogenomics in the clinical management of patients with treatment-resistant schizophrenia.

Dynamic expression of genes associated with schizophrenia and bipolar disorder across development.

Common genetic variation contributes a substantial proportion of risk for both schizophrenia and bipolar disorder. Furthermore, there is evidence of significant, but not complete, overlap in genetic risk between the two disorders. It has been hypothesised that genetic variants conferring risk for these disorders do so by influencing brain development, leading to the later emergence of symptoms. The comparative profile of risk gene expression for schizophrenia and bipolar disorder across development over different brain regions however remains unclear. Using genotypes derived from genome-wide associations studies of the largest available cohorts of patients and control subjects, we investigated whether genes enriched for schizophrenia and bipolar disorder association show a bias for expression across any of 13 developmental stages in prefrontal cortical and subcortical brain regions. We show that genetic association with schizophrenia is positively correlated with expression in the prefrontal cortex during early midfetal development and early infancy, and negatively correlated with expression during late childhood, which stabilises in adolescence. In contrast, risk-associated genes for bipolar disorder did not exhibit a bias towards expression at any prenatal stage, although the pattern of postnatal expression was similar to that of schizophrenia. These results highlight the dynamic expression of genes harbouring risk for schizophrenia and bipolar disorder across prefrontal cortex development and support the hypothesis that prenatal neurodevelopmental events are more strongly associated with schizophrenia than bipolar disorder.

Targeted Sequencing of 10,198 Samples Confirms Abnormalities in Neuronal Activity and Implicates Voltage-Gated Sodium Channels in Schizophrenia Pathogenesis.

BACKGROUND: Sequencing studies have pointed to the involvement in schizophrenia of rare coding variants in neuronally expressed genes, including activity-regulated cytoskeleton-associated protein (ARC) and N-methyl-D-aspartate receptor (NMDAR) complexes; however, larger samples are required to reveal novel genes and specific biological mechanisms. METHODS: We sequenced 187 genes, selected for prior evidence of association with schizophrenia, in a new dataset of 5207 cases and 4991 controls. Included among these genes were members of ARC and NMDAR postsynaptic protein complexes, as well as voltage-gated sodium and calcium channels. We performed a rare variant meta-analysis with published sequencing data for a total of 11,319 cases, 15,854 controls, and 1136 trios. RESULTS: While no individual gene was significantly associated with schizophrenia after genome-wide correction for multiple testing, we strengthen the evidence that rare exonic variants in the ARC (p = 4.0 × 10-4) and NMDAR (p = 1.7 × 10-5) synaptic complexes are risk factors for schizophrenia. In addition, we found that loss-of-function variants and missense variants at paralog-conserved sites were enriched in voltage-gated sodium channels, particularly the alpha subunits (p = 8.6 × 10-4). CONCLUSIONS: In one of the largest sequencing studies of schizophrenia to date, we provide novel evidence that multiple voltage-gated sodium channels are involved in schizophrenia pathogenesis and confirm the involvement of ARC and NMDAR postsynaptic complexes.

Publisher Correction: Common schizophrenia alleles are enriched in mutation-intolerant genes and in regions under strong background selection.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Common schizophrenia alleles are enriched in mutation-intolerant genes and in regions under strong background selection.

Schizophrenia is a debilitating psychiatric condition often associated with poor quality of life and decreased life expectancy. Lack of progress in improving treatment outcomes has been attributed to limited knowledge of the underlying biology, although large-scale genomic studies have begun to provide insights. We report a new genome-wide association study of schizophrenia (11,260 cases and 24,542 controls), and through meta-analysis with existing data we identify 50 novel associated loci and 145 loci in total. Through integrating genomic fine-mapping with brain expression and chromosome conformation data, we identify candidate causal genes within 33 loci. We also show for the first time that the common variant association signal is highly enriched among genes that are under strong selective pressures. These findings provide new insights into the biology and genetic architecture of schizophrenia, highlight the importance of mutation-intolerant genes and suggest a mechanism by which common risk variants persist in the population.

The proteomes of neurotransmitter receptor complexes form modular networks with distributed functionality underlying plasticity and behaviour.

Neuronal synapses play fundamental roles in information processing, behaviour and disease. Neurotransmitter receptor complexes, such as the mammalian N-methyl-D-aspartate receptor complex (NRC/MASC) comprising 186 proteins, are major components of the synapse proteome. Here we investigate the organisation and function of NRC/MASC using a systems biology approach. Systematic annotation showed that the complex contained proteins implicated in a wide range of cognitive processes, synaptic plasticity and psychiatric diseases. Protein domains were evolutionarily conserved from yeast, but enriched with signalling domains associated with the emergence of multicellularity. Mapping of protein-protein interactions to create a network representation of the complex revealed that simple principles underlie the functional organisation of both proteins and their clusters, with modularity reflecting functional specialisation. The known functional roles of NRC/MASC proteins suggest the complex co-ordinates signalling to diverse effector pathways underlying neuronal plasticity. Importantly, using quantitative data from synaptic plasticity experiments, our model correctly predicts robustness to mutations and drug interference. These studies of synapse proteome organisation suggest that molecular networks with simple design principles underpin synaptic signalling properties with important roles in physiology, behaviour and disease.

Arc Requires PSD95 for Assembly into Postsynaptic Complexes Involved with Neural Dysfunction and Intelligence.

Arc is an activity-regulated neuronal protein, but little is known about its interactions, assembly into multiprotein complexes, and role in human disease and cognition. We applied an integrated proteomic and genetic strategy by targeting a tandem affinity purification (TAP) tag and Venus fluorescent protein into the endogenous Arc gene in mice. This allowed biochemical and proteomic characterization of native complexes in wild-type and knockout mice. We identified many Arc-interacting proteins, of which PSD95 was the most abundant. PSD95 was essential for Arc assembly into 1.5-MDa complexes and activity-dependent recruitment to excitatory synapses. Integrating human genetic data with proteomic data showed that Arc-PSD95 complexes are enriched in schizophrenia, intellectual disability, autism, and epilepsy mutations and normal variants in intelligence. We propose that Arc-PSD95 postsynaptic complexes potentially affect human cognitive function.

Novel Findings from CNVs Implicate Inhibitory and Excitatory Signaling Complexes in Schizophrenia.

We sought to obtain novel insights into schizophrenia pathogenesis by exploiting the association between the disorder and chromosomal copy number (CNV) burden. We combined data from 5,745 cases and 10,675 controls with other published datasets containing genome-wide CNV data. In this much-enlarged sample of 11,355 cases and 16,416 controls, we show for the first time that case CNVs are enriched for genes involved in GABAergic neurotransmission. Consistent with non-genetic reports of GABAergic deficits in schizophrenia, our findings now show disrupted GABAergic signaling is of direct causal relevance, rather than a secondary effect or due to confounding. Additionally, we independently replicate and greatly extend previous findings of CNV enrichment among genes involved in glutamatergic signaling. Given the strong functional links between the major inhibitory GABAergic and excitatory glutamatergic systems, our findings converge on a broad, coherent set of pathogenic processes, providing firm foundations for studies aimed at dissecting disease mechanisms.

Clustered coding variants in the glutamate receptor complexes of individuals with schizophrenia and bipolar disorder.

Current models of schizophrenia and bipolar disorder implicate multiple genes, however their biological relationships remain elusive. To test the genetic role of glutamate receptors and their interacting scaffold proteins, the exons of ten glutamatergic 'hub' genes in 1304 individuals were re-sequenced in case and control samples. No significant difference in the overall number of non-synonymous single nucleotide polymorphisms (nsSNPs) was observed between cases and controls. However, cluster analysis of nsSNPs identified two exons encoding the cysteine-rich domain and first transmembrane helix of GRM1 as a risk locus with five mutations highly enriched within these domains. A new splice variant lacking the transmembrane GPCR domain of GRM1 was discovered in the human brain and the GRM1 mutation cluster could perturb the regulation of this variant. The predicted effect on individuals harbouring multiple mutations distributed in their ten hub genes was also examined. Diseased individuals possessed an increased load of deleteriousness from multiple concurrent rare and common coding variants. Together, these data suggest a disease model in which the interplay of compound genetic coding variants, distributed among glutamate receptors and their interacting proteins, contribute to the pathogenesis of schizophrenia and bipolar disorders.

Neurotransmitters drive combinatorial multistate postsynaptic density networks.

The mammalian postsynaptic density (PSD) comprises a complex collection of approximately 1100 proteins. Despite extensive knowledge of individual proteins, the overall organization of the PSD is poorly understood. Here, we define maps of molecular circuitry within the PSD based on phosphorylation of postsynaptic proteins. Activation of a single neurotransmitter receptor, the N-methyl-D-aspartate receptor (NMDAR), changed the phosphorylation status of 127 proteins. Stimulation of ionotropic and metabotropic glutamate receptors and dopamine receptors activated overlapping networks with distinct combinatorial phosphorylation signatures. Using peptide array technology, we identified specific phosphorylation motifs and switching mechanisms responsible for the integration of neurotransmitter receptor pathways and their coordination of multiple substrates in these networks. These combinatorial networks confer high information-processing capacity and functional diversity on synapses, and their elucidation may provide new insights into disease mechanisms and new opportunities for drug discovery.

Reconstructing protein complexes: from proteomics to systems biology.

Modern high throughput technologies in biological science often create lists of interesting molecules. The challenge is to reconstruct a descriptive model from these lists that reflects the underlying biological processes as accurately as possible. Once we have such a model or network, what can we learn from it? Specifically, given that we are interested in some biological process associated with the model, what new properties can we predict and subsequently test? Here, we describe, at an introductory level, a range of bioinformatics techniques that can be systematically applied to proteomic datasets. When combined, these methods give us a global overview of the network and the properties of the proteins and their interactions. These properties can then be used to predict functional pathways within the network and to examine substructure. To illustrate the application of these methods, we draw upon our own work concerning a complex of 186 proteins found in neuronal synapses in mammals. The techniques discussed are generally applicable and could be used to examine lists of proteins involved with the biological response to electric or magnetic fields.