Sex Differences in the Human Brain Transcriptome of Cases With Schizophrenia

While schizophrenia differs between males and females in age of onset, symptomatology and the course of the disease, the molecular mechanisms underlying these differences remain uncharacterized. In order to address questions about the sex-specific effects of schizophrenia, we performed a large-scale transcriptome analysis of RNA-seq data from 437 controls and 341 cases from two distinct cohorts from the CommonMind Consortium. Analysis across the cohorts identifies a reproducible gene expression signature of schizophrenia that is highly concordant with previous work. Differential expression across sex is reproducible across cohorts and identifies X- and Y-linked genes, as well as those involved in dosage compensation. Intriguingly, the sex expression signature is also enriched for genes involved in neurexin family protein binding and synaptic organization. Differential expression analysis testing a sex-by-diagnosis interaction effect did not identify any genome-wide signature after multiple testing corrections. Gene coexpression network analysis was performed to reduce dimensionality and elucidate interactions among genes. We found enrichment of co-expression modules for sex-by-diagnosis differential expression signatures, which were highly reproducible across the two cohorts and involve a number of diverse pathways, including neural nucleus development, neuron projection morphogenesis, and regulation of neural precursor cell proliferation. Overall, our results indicate that the effect size of sex differences in schizophrenia gene expression signatures is small and underscore the challenge of identifying robust sex-by-diagnosis signatures, which will require future analyses in larger cohorts.

[1]  N. Wray,et al.  Examining Sex-Differentiated Genetic Effects Across Neuropsychiatric and Behavioral Traits , 2020, Biological Psychiatry.

[2]  Gabriel E. Hoffman,et al.  dream: Powerful differential expression analysis for repeated measures designs , 2018, bioRxiv.

[3]  S. Ripke,et al.  Mapping genomic loci prioritises genes and implicates synaptic biology in schizophrenia , 2020, medRxiv.

[4]  Ayellet V. Segrè,et al.  The impact of sex on gene expression across human tissues , 2020, Science.

[5]  Janice M. Fullerton,et al.  Sex-Dependent Shared and Non-Shared Genetic Architecture Across Mood and Psychotic Disorders , 2020, bioRxiv.

[6]  P. Roussos,et al.  dream: Powerful differential expression analysis for repeated measures designs. , 2020, Bioinformatics.

[7]  K. Roche,et al.  A Cluster of Autism-Associated Variants on X-Linked NLGN4X Functionally Resemble NLGN4Y , 2020, Neuron.

[8]  K. Schlett,et al.  Dendritic spine morphology and memory formation depend on postsynaptic Caskin proteins , 2019, Scientific Reports.

[9]  B. Stranger,et al.  SHARED SEX-DEPENDENT GENETIC EFFECTS ACROSS NEUROPSYCHIATRIC AND BEHAVIORAL TRAITS , 2019, European Neuropsychopharmacology.

[10]  Kelsey S. Montgomery,et al.  CommonMind Consortium provides transcriptomic and epigenomic data for Schizophrenia and Bipolar Disorder , 2019, Scientific Data.

[11]  Emily E. Burke,et al.  Regional Heterogeneity in Gene Expression, Regulation, and Coherence in the Frontal Cortex and Hippocampus across Development and Schizophrenia , 2019, Neuron.

[12]  James C. Hu,et al.  The Gene Ontology Resource: 20 years and still GOing strong , 2019 .

[13]  Gregory J. Hunt,et al.  Dtangle: Accurate and Robust Cell Type Deconvolution , 2018, Bioinform..

[14]  The Gene Ontology Consortium,et al.  The Gene Ontology Resource: 20 years and still GOing strong , 2018, Nucleic Acids Res..

[15]  B. Stranger,et al.  The role of sex in the genomics of human complex traits , 2018, Nature Reviews Genetics.

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

[17]  M. Stephens,et al.  Flexible statistical methods for estimating and testing effects in genomic studies with multiple conditions , 2016, Nature Genetics.

[18]  Masahiko Watanabe,et al.  Distribution of Caskin1 protein and phenotypic characterization of its knockout mice using a comprehensive behavioral test battery , 2018, Molecular Brain.

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

[20]  C. Kellendonk,et al.  Can we use mice to study schizophrenia? , 2018, Philosophical Transactions of the Royal Society B: Biological Sciences.

[21]  Paolo Di Tommaso,et al.  Nextflow enables reproducible computational workflows , 2017, Nature Biotechnology.

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

[23]  E. Schadt,et al.  variancePartition: interpreting drivers of variation in complex gene expression studies , 2016, bioRxiv.

[24]  Matthew Stephens,et al.  False discovery rates: a new deal , 2016, bioRxiv.

[25]  Ralf Zimmer,et al.  Bioconductor’s EnrichmentBrowser: seamless navigation through combined results of set- & network-based enrichment analysis , 2016, BMC Bioinformatics.

[26]  L. Feuk,et al.  Spatial sexual dimorphism of X and Y homolog gene expression in the human central nervous system during early male development , 2016, Biology of Sex Differences.

[27]  E. Koonin,et al.  Substantial DNA methylation differences between two major neuronal subtypes in human brain , 2015, Nucleic acids research.

[28]  Claudia Biermann,et al.  Mathematical Methods Of Statistics , 2016 .

[29]  Cong Liu,et al.  Meta-analysis of sex differences in gene expression in schizophrenia , 2016, BMC Systems Biology.

[30]  Bin Zhang,et al.  Multiscale Embedded Gene Co-expression Network Analysis , 2015, PLoS Comput. Biol..

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

[32]  Matthew E. Ritchie,et al.  limma powers differential expression analyses for RNA-sequencing and microarray studies , 2015, Nucleic acids research.

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

[34]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

[35]  Samuel Bernard,et al.  Neurogenesis in the Striatum of the Adult Human Brain , 2014, Cell.

[36]  Wei Shi,et al.  featureCounts: an efficient general purpose program for assigning sequence reads to genomic features , 2013, Bioinform..

[37]  Charity W. Law,et al.  voom: precision weights unlock linear model analysis tools for RNA-seq read counts , 2014, Genome Biology.

[38]  K. Gerald van den Boogaart,et al.  Analyzing Compositional Data with R , 2013 .

[39]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[40]  A. Balazs,et al.  Complex formation of EphB1/Nck/Caskin1 leads to tyrosine phosphorylation and structural changes of the Caskin1 SH3 domain , 2012, Cell Communication and Signaling.

[41]  C. Disteche Dosage compensation of the sex chromosomes. , 2012, Annual review of genetics.

[42]  G. Smyth,et al.  Camera: a competitive gene set test accounting for inter-gene correlation , 2012, Nucleic acids research.

[43]  C A Jones,et al.  Animal models of schizophrenia , 2011, British journal of pharmacology.

[44]  H. Broihier,et al.  The Cytoplasmic Adaptor Protein Caskin Mediates Lar Signal Transduction during Drosophila Motor Axon Guidance , 2011, The Journal of Neuroscience.

[45]  J. Goldstein,et al.  Sex differences in schizophrenia , 2010, International review of psychiatry.

[46]  Wolfgang Viechtbauer,et al.  Conducting Meta-Analyses in R with the metafor Package , 2010 .

[47]  Mark D. Robinson,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[48]  M. Robinson,et al.  A scaling normalization method for differential expression analysis of RNA-seq data , 2010, Genome Biology.

[49]  Robert Kiss,et al.  High levels of structural disorder in scaffold proteins as exemplified by a novel neuronal protein, CASK‐interactive protein1 , 2009, The FEBS journal.

[50]  G. Strauss,et al.  Sex differences in schizophrenia and other psychotic disorders: a 20-year longitudinal study of psychosis and recovery. , 2008, Comprehensive psychiatry.

[51]  T. Südhof Neuroligins and neurexins link synaptic function to cognitive disease , 2008, Nature.

[52]  J. DeFelipe,et al.  Gender differences in human cortical synaptic density , 2008, Proceedings of the National Academy of Sciences.

[53]  T. Miyakawa,et al.  Schizophrenia-Relevant Behavioral Testing in Rodent Models: A Uniquely Human Disorder? , 2006, Biological Psychiatry.

[54]  Hidetoshi Shimodaira,et al.  Pvclust: an R package for assessing the uncertainty in hierarchical clustering , 2006, Bioinform..

[55]  C. Woolley,et al.  Sexually Dimorphic Synaptic Organization of the Medial Amygdala , 2005, The Journal of Neuroscience.

[56]  Pablo Tamayo,et al.  Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[57]  John D. Storey,et al.  Statistical significance for genomewide studies , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[58]  T. Südhof,et al.  CASK Participates in Alternative Tripartite Complexes in which Mint 1 Competes for Binding with Caskin 1, a Novel CASK-Binding Protein , 2002, The Journal of Neuroscience.

[59]  W. Huttner,et al.  Chlorate--a potent inhibitor of protein sulfation in intact cells. , 1986, Biochemical and biophysical research communications.

[60]  G. Schwarz Estimating the Dimension of a Model , 1978 .