Prefrontal Coexpression of Schizophrenia Risk Genes Is Associated With Treatment Response in Patients

BACKGROUND Gene coexpression networks are relevant to functional and clinical translation of schizophrenia risk genes. We hypothesized that schizophrenia risk genes converge into coexpression pathways that may be associated with gene regulation mechanisms and with response to treatment in patients with schizophrenia. METHODS We identified gene coexpression networks in two prefrontal cortex postmortem RNA sequencing datasets (n = 688) and replicated them in four more datasets (n = 1295). We identified and replicated (p values < .001) a single module enriched for schizophrenia risk loci (13 risk genes in 10 loci). In silico screening of potential regulators of the schizophrenia risk module via bioinformatic analyses identified two transcription factors and three microRNAs associated with the risk module. To translate postmortem information into clinical phenotypes, we identified polymorphisms predicting coexpression and combined them to obtain an index approximating module coexpression (Polygenic Coexpression Index [PCI]). RESULTS The PCI-coexpression association was successfully replicated in two independent brain transcriptome datasets (n = 131; p values < .05). Finally, we tested the association between the PCI and short-term treatment response in two independent samples of patients with schizophrenia treated with olanzapine (n = 167). The PCI was associated with treatment response in the positive symptom domain in both clinical cohorts (p values < .05). CONCLUSIONS In summary, our findings in 1983 samples of human postmortem prefrontal cortex show that coexpression of a set of genes enriched for schizophrenia risk genes is relevant to treatment response. This coexpression pathway may be coregulated by transcription factors and microRNA associated with it.

[1]  S. Kay,et al.  The positive and negative syndrome scale (PANSS) for schizophrenia. , 1987, Schizophrenia bulletin.

[2]  T. Insel.  The NIMH Research Domain Criteria (RDoC) Project: precision medicine for psychiatry. , 2014, The American journal of psychiatry.

[3]  J. Rappsilber,et al.  Expanded View Figures , 2017 .

[4]  G. Pergola,et al.  Combined effect of genetic variants in the GluN2B coding gene (GRIN2B) on prefrontal function during working memory performance , 2015, Psychological Medicine.

[5]  M. Chiang,et al.  Rosiglitazone promotes neurite outgrowth and mitochondrial function in N2A cells via PPARgamma pathway. , 2014, Mitochondrion.

[6]  Manolis Kellis,et al.  HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants , 2011, Nucleic Acids Res..

[7]  Thomas Ragg,et al.  The RIN: an RNA integrity number for assigning integrity values to RNA measurements , 2006, BMC Molecular Biology.

[8]  M. Owen,et al.  The implications of the shared genetics of psychiatric disorders , 2016, Nature Medicine.

[9]  P. Michalak.  Coexpression, coregulation, and cofunctionality of neighboring genes in eukaryotic genomes. , 2008, Genomics.

[10]  A. Malhotra,et al.  Schizophrenia Polygenic Risk Score as a Predictor of Antipsychotic Efficacy in First-Episode Psychosis. , 2019, The American journal of psychiatry.

[11]  H. Okano,et al.  Differential gene expression profiles in neurons generated from lymphoblastoid B-cell line-derived iPS cells from monozygotic twin cases with treatment-resistant schizophrenia and discordant responses to clozapine , 2017, Schizophrenia Research.

[12]  Giuseppe Blasi,et al.  Schizophrenia polygenic risk score predicts mnemonic hippocampal activity , 2018, Brain : a journal of neurology.

[13]  Jianxin Shi,et al.  Genome-wide association study of clinical dimensions of schizophrenia: polygenic effect on disorganized symptoms. , 2012, The American journal of psychiatry.

[14]  M. Boehnke,et al.  So many correlated tests, so little time! Rapid adjustment of P values for multiple correlated tests. , 2007, American journal of human genetics.

[15]  Ellen T. Gelfand,et al.  A Novel Approach to High-Quality Postmortem Tissue Procurement: The GTEx Project , 2015, Biopreservation and biobanking.

[16]  Paul J. Harrison,et al.  Schizophrenia genes, gene expression, and neuropathology: on the matter of their convergence , 2005, Molecular Psychiatry.

[17]  P. Sullivan,et al.  Schizophrenia as a complex trait: evidence from a meta-analysis of twin studies. , 2003, Archives of general psychiatry.

[18]  Alexander Gusev,et al.  Transcriptome-wide association study of schizophrenia and chromatin activity yields mechanistic disease insights , 2016, Nature Genetics.

[19]  Benjamin A. Logsdon,et al.  Gene Expression Elucidates Functional Impact of Polygenic Risk for Schizophrenia , 2016, Nature Neuroscience.

[20]  Jeffrey T Leek,et al.  Developmental and genetic regulation of the human cortex transcriptome illuminate schizophrenia pathogenesis , 2018, Nature Neuroscience.

[21]  Daniel R Weinberger,et al.  Identification and prioritization of gene sets associated with schizophrenia risk by co-expression network analysis in human brain , 2018, bioRxiv.

[22]  David R. O'Brien,et al.  Cell Type-Specific Expression Analysis to Identify Putative Cellular Mechanisms for Neurogenetic Disorders , 2014, The Journal of Neuroscience.

[23]  Gerome Breen,et al.  Genetic identification of brain cell types underlying schizophrenia , 2017, Nature Genetics.

[24]  Y Wang,et al.  Genetic evidence for role of integration of fast and slow neurotransmission in schizophrenia , 2017, Molecular Psychiatry.

[25]  Annie W Shieh,et al.  Transcriptome-wide isoform-level dysregulation in ASD, schizophrenia, and bipolar disorder , 2018, Science.

[26]  Rui Luo,et al.  Is My Network Module Preserved and Reproducible? , 2011, PLoS Comput. Biol..

[27]  D. Weinberger,et al.  DRD2 co-expression network and a related polygenic index predict imaging, behavioral and clinical phenotypes linked to schizophrenia , 2017, Translational Psychiatry.

[28]  S. Cichon,et al.  Identification of a Bipolar Disorder Vulnerable Gene CHDH at 3p21.1 , 2017, Molecular Neurobiology.

[29]  Paul J. Harrison,et al.  Genetic Neuropathology of Schizophrenia: New Approaches to an Old Question and New Uses for Postmortem Human Brains , 2011, Biological Psychiatry.

[30]  F. Sambataro,et al.  Thalamic connectivity measured with fMRI is associated with a polygenic index predicting thalamo-prefrontal gene co-expression , 2019, Brain Structure and Function.

[31]  John P. Rice,et al.  Genomic Dissection of Bipolar Disorder and Schizophrenia, Including 28 Subphenotypes , 2017, Cell.

[32]  Vincenzo De Luca,et al.  Polygenic risk score prediction of antipsychotic dosage in schizophrenia , 2016, Schizophrenia Research.

[33]  Prashant S. Emani,et al.  Comprehensive functional genomic resource and integrative model for the human brain , 2018, Science.

[34]  Giulio Genovese,et al.  Schizophrenia risk from complex variation of complement component 4 , 2016, Nature.

[35]  Peter B. Jones,et al.  Phenotypic Manifestation of Genetic Risk for Schizophrenia During Adolescence in the General Population. , 2016, JAMA psychiatry.

[36]  Yang I Li,et al.  An Expanded View of Complex Traits: From Polygenic to Omnigenic , 2017, Cell.

[37]  Peter Kirsch,et al.  Hippocampal and Frontolimbic Function as Intermediate Phenotype for Psychosis: Evidence from Healthy Relatives and a Common Risk Variant in CACNA1C , 2014, Biological Psychiatry.

[38]  Gene W. Yeo,et al.  MicroRNA-101 Regulates Multiple Developmental Programs to Constrain Excitation in Adult Neural Networks , 2016, Neuron.

[39]  S. Horvath,et al.  Functional organization of the transcriptome in human brain , 2008, Nature Neuroscience.

[40]  J. Rapoport,et al.  Dysregulated protocadherin-pathway activity as an intrinsic defect in iPSC-derived cortical interneurons from patients with schizophrenia , 2019, Nature Neuroscience.

[41]  A. Børglum,et al.  Schizophrenia risk variants affecting microRNA function and site-specific regulation of NT5C2 by miR-206 , 2016, European Neuropsychopharmacology.

[42]  C. Spencer,et al.  Biological Insights From 108 Schizophrenia-Associated Genetic Loci , 2014, Nature.

[43]  J. Lieberman,et al.  The National Institute of Mental Health Clinical Antipsychotic Trials of Intervention Effectiveness (CATIE) project: schizophrenia trial design and protocol development. , 2003, Schizophrenia bulletin.

[44]  S. Horvath,et al.  Shared molecular neuropathology across major psychiatric disorders parallels polygenic overlap , 2016, Science.

[45]  J. Kleinman,et al.  NURR1 and ERR1 Modulate the Expression of Genes of a DRD2 Coexpression Network Enriched for Schizophrenia Risk , 2019, The Journal of Neuroscience.

[46]  Joris M. Mooij,et al.  MAGMA: Generalized Gene-Set Analysis of GWAS Data , 2015, PLoS Comput. Biol..

[47]  Daniel H. Geschwind,et al.  Systems biology and gene networks in neurodevelopmental and neurodegenerative disorders , 2015, Nature Reviews Genetics.

[48]  Daniel R. Weinberger,et al.  Functional neuroimaging and schizophrenia: a view towards effective connectivity modeling and polygenic risk , 2013, Dialogues in clinical neuroscience.

[49]  R. Keefe,et al.  Polygenic signal for symptom dimensions and cognitive performance in patients with chronic schizophrenia , 2018, Schizophrenia Research: Cognition.

[50]  Jeffrey T Leek,et al.  qSVA framework for RNA quality correction in differential expression analysis , 2017, Proceedings of the National Academy of Sciences of the United States of America.

[51]  Graziano Pesole,et al.  Pscan: finding over-represented transcription factor binding site motifs in sequences from co-regulated or co-expressed genes , 2009, Nucleic Acids Res..

[52]  Terence P. Speed,et al.  Systematic noise degrades gene co-expression signals but can be corrected , 2015, BMC Bioinformatics.

[53]  Pietro Liò,et al.  The BioMart community portal: an innovative alternative to large, centralized data repositories , 2015, Nucleic Acids Res..

[54]  G. Tseng,et al.  Beyond modules and hubs: the potential of gene coexpression networks for investigating molecular mechanisms of complex brain disorders , 2014, Genes, brain, and behavior.

[55]  Takashi Uehara,et al.  Specific gene expression patterns of 108 schizophrenia-associated loci in cortex , 2016, Schizophrenia Research.

[56]  S. Horvath,et al.  Statistical Applications in Genetics and Molecular Biology , 2011 .

[57]  Monica Driscoll,et al.  A long journey to reproducible results , 2017, Nature.

[58]  Danny J. Smith,et al.  The use of polygenic risk scores to identify phenotypes associated with genetic risk of schizophrenia: Systematic review , 2017, Schizophrenia Research.

[59]  D. Geschwind,et al.  Histone Acetylome-wide Association Study of Autism Spectrum Disorder , 2016, Cell.

[60]  D. Weinberger,et al.  From neuropathology to neurodevelopment , 1995, The Lancet.

[61]  Qi Dong,et al.  The Effects of CACNA1C Gene Polymorphism on Spatial Working Memory in Both Healthy Controls and Patients with Schizophrenia or Bipolar Disorder , 2012, Neuropsychopharmacology.

[62]  Warren W. Kretzschmar,et al.  Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression , 2017, Nature Genetics.

[63]  N. Amoroso,et al.  Transcriptomic context of DRD1 is associated with prefrontal activity and behavior during working memory , 2018, Proceedings of the National Academy of Sciences.

[64]  Sisi Jiang,et al.  Evidence for Association of Cell Adhesion Molecules Pathway and NLGN1 Polymorphisms with Schizophrenia in Chinese Han Population , 2015, PloS one.

[65]  Patrick F. Sullivan,et al.  Comparative Genomic Evidence for the Involvement of Schizophrenia Risk Genes in Antipsychotic Effects , 2017, Molecular Psychiatry.

[66]  Jing Chen,et al.  ToppGene Suite for gene list enrichment analysis and candidate gene prioritization , 2009, Nucleic Acids Res..

[67]  A. Ramasamy,et al.  Quality control parameters on a large dataset of regionally dissected human control brains for whole genome expression studies , 2011, Journal of neurochemistry.

[68]  Jun Li,et al.  Five novel loci associated with antipsychotic treatment response in patients with schizophrenia: a genome-wide association study. , 2018, The lancet. Psychiatry.

[69]  Nishant Singh,et al.  Differentiation Induces Dramatic Changes in miRNA Profile, Where Loss of Dicer Diverts Differentiating SH-SY5Y Cells Toward Senescence , 2016, Molecular Neurobiology.

[70]  F. Sambataro,et al.  Genetic Variation of a DRD2 Co‐expression Network is Associated with Changes in Prefrontal Function After D2 Receptors Stimulation , 2019, Cerebral cortex.

[71]  Luis de la Torre Ubieta,et al.  Genome-wide changes in lncRNA, splicing, and regional gene expression patterns in autism , 2016, Nature.

[72]  J. Callicott,et al.  Interaction of COMT (Val(108/158)Met) genotype and olanzapine treatment on prefrontal cortical function in patients with schizophrenia. , 2004, The American journal of psychiatry.

[73]  C. Gasse,et al.  Polygenic Risk Score for Schizophrenia and Treatment-Resistant Schizophrenia , 2017, Schizophrenia bulletin.

[74]  Laura J. Scott,et al.  Psychiatric genome-wide association study analyses implicate neuronal, immune and histone pathways , 2015, Nature Neuroscience.

[75]  A. Meyer-Lindenberg,et al.  Intermediate phenotypes and genetic mechanisms of psychiatric disorders , 2006, Nature Reviews Neuroscience.

[76]  K. Kendler.  What psychiatric genetics has taught us about the nature of psychiatric illness and what is left to learn , 2013, Molecular Psychiatry.

[77]  Jun S. Liu,et al.  The Genotype-Tissue Expression (GTEx) pilot analysis: Multitissue gene regulation in humans , 2015, Science.

[78]  Tilo Kircher,et al.  The impact of a CACNA1C gene polymorphism on learning and hippocampal formation in healthy individuals: A diffusion tensor imaging study , 2014, NeuroImage.

[79]  P. McGuffin,et al.  Toward Behavioral Genomics , 2001, Science.

[80]  A. Ultsch,et al.  What do all the (human) micro-RNAs do? , 2014, BMC Genomics.

[81]  Yuan Wu,et al.  A Novel Schizophrenia Diagnostic Model Based on Statistically Significant Changes in Gene Methylation in Specific Brain Regions , 2020, BioMed research international.

[82]  Manolis Kellis,et al.  HaploReg v4: systematic mining of putative causal variants, cell types, regulators and target genes for human complex traits and disease , 2015, Nucleic Acids Res..

[83]  Doug Speed,et al.  Systems genetics identifies a convergent gene network for cognition and neurodevelopmental disease , 2015, Nature Neuroscience.