Differential gene regulatory pattern in the human brain from schizophrenia using transcriptomic-causal network

Background Common and complex traits are the consequence of the interaction and regulation of multiple genes simultaneously, therefore characterizing the interconnectivity of genes is essential to unravel the underlying biological networks. However, the focus of many studies is on the differential expression of individual genes or on co-expression analysis. Methods Going beyond analysis of one gene at a time, we systematically integrated transcriptomics, genotypes and Hi-C data to identify interconnectivities among individual genes as a causal network. We utilized different machine learning techniques to extract information from the network and identify differential regulatory pattern between cases and controls. We used data from the Allen Brain Atlas for replication. Results Employing the integrative systems approach on the data from CommonMind Consortium showed that gene transcription is controlled by genetic variants proximal to the gene (cis-regulatory factors), and transcribed distal genes (trans-regulatory factors). We identified differential gene regulatory patterns in SCZ-cases versus controls and novel SCZ-associated genes that may play roles in the disorder since some of them are primary expressed in human brain. In addition, we observed genes known associated with SCZ are not likely (OR = 0.59) to have high impacts (degree > 3) on the network. Conclusions Causal networks could reveal underlying patterns and the role of genes individually and as a group. Establishing principles that govern relationships between genes provides a mechanistic understanding of the dysregulated gene transcription patterns in SCZ and creates more efficient experimental designs for further studies. This information cannot be obtained by studying a single gene at the time.

[1]  Eric S. Lander,et al.  Hi-C: A Method to Study the Three-dimensional Architecture of Genomes. , 2010, Journal of visualized experiments : JoVE.

[2]  M. Tobin,et al.  Mendelian Randomisation and Causal Inference in Observational Epidemiology , 2008, PLoS medicine.

[3]  Pall I. Olason,et al.  Common variants conferring risk of schizophrenia , 2009, Nature.

[4]  Kathryn Roeder,et al.  Differential Activity of Transcribed Enhancers in the Prefrontal Cortex of 537 cases with Schizophrenia and Controls , 2018, Molecular Psychiatry.

[5]  H. Abdi,et al.  Multiple Correspondence Analysis , 2006 .

[6]  S. Rhee,et al.  Towards revealing the functions of all genes in plants. , 2014, Trends in plant science.

[7]  Juliana Costa-Silva,et al.  RNA-Seq differential expression analysis: An extended review and a software tool , 2017, PloS one.

[8]  Robert Freedman,et al.  The human CHRNA7 and CHRFAM7A genes: A review of the genetics, regulation, and function , 2015, Neuropharmacology.

[9]  B. L. Roux,et al.  Multiple Correspondence Analysis , 2009 .

[10]  Damian Szklarczyk,et al.  STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets , 2018, Nucleic Acids Res..

[11]  P. Visscher,et al.  Common polygenic variation contributes to risk of schizophrenia and bipolar disorder , 2009, Nature.

[12]  H. Nasrallah,et al.  Atypical antipsychotic-induced metabolic side effects: insights from receptor-binding profiles , 2008, Molecular Psychiatry.

[13]  A. Dawid FUNDAMENTALS OF STATISTICAL CAUSALITY , 2007 .

[14]  Ahmad Samiei,et al.  Generating a robust statistical causal structure over 13 cardiovascular disease risk factors using genomics data , 2016, J. Biomed. Informatics.

[15]  James Y. Zou Analysis of protein-coding genetic variation in 60,706 humans , 2015, Nature.

[16]  Evan Fletcher,et al.  White Matter Changes Compromise Prefrontal Cortex Function in Healthy Elderly Individuals , 2006 .

[17]  G. von Heijne,et al.  Tissue-based map of the human proteome , 2015, Science.

[18]  Hiroyuki Ogata,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 1999, Nucleic Acids Res..

[19]  Rachel B. Brem,et al.  Stitching together Multiple Data Dimensions Reveals Interacting Metabolomic and Transcriptomic Networks That Modulate Cell Regulation , 2012, PLoS biology.

[20]  Etienne Sibille,et al.  Biological substrates underpinning diagnosis of major depression. , 2013, The international journal of neuropsychopharmacology.

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

[22]  Michael Wainberg,et al.  Opportunities and challenges for transcriptome-wide association studies , 2019, Nature Genetics.

[23]  João Pedro de Magalhães,et al.  Gene co-expression analysis for functional classification and gene–disease predictions , 2017, Briefings Bioinform..

[24]  Constantin F. Aliferis,et al.  The max-min hill-climbing Bayesian network structure learning algorithm , 2006, Machine Learning.

[25]  Allan R. Jones,et al.  The Allen Human Brain Atlas Comprehensive gene expression mapping of the human brain , 2012, Trends in Neurosciences.

[26]  D. Mohr,et al.  Major depressive disorder , 2016, Nature Reviews Disease Primers.

[27]  Markus Perola,et al.  Metabonomic, transcriptomic, and genomic variation of a population cohort , 2010, Molecular systems biology.

[28]  Steve Horvath,et al.  Using genetic markers to orient the edges in quantitative trait networks: The NEO software , 2008, BMC Systems Biology.

[29]  Tatiana A. Tatusova,et al.  Entrez Gene: gene-centered information at NCBI , 2004, Nucleic Acids Res..

[30]  E M Wijsman,et al.  Meta-analysis of 32 genome-wide linkage studies of schizophrenia , 2009, Molecular Psychiatry.

[31]  R. Kikinis,et al.  A review of diffusion tensor imaging studies in schizophrenia. , 2007, Journal of psychiatric research.

[32]  Christos Pantelis,et al.  The ubiquitin proteasome system and schizophrenia. , 2020, The lancet. Psychiatry.

[33]  M. Ashburner,et al.  Gene Ontology: tool for the unification of biology , 2000, Nature Genetics.

[34]  Cathy H. Wu,et al.  UniProt: the Universal Protein knowledgebase , 2004, Nucleic Acids Res..

[35]  Susumu Goto,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 2000, Nucleic Acids Res..

[36]  Daning Lu,et al.  Chromosome conformation elucidates regulatory relationships in developing human brain , 2016, Nature.

[37]  Eytan Domany,et al.  Comprehensive post mortem brain samples analysis detects global reduction of multiple proteasome subunits expression in schizophrenia , 2019, bioRxiv.

[38]  David B. Dunson,et al.  A hybrid bayesian approach for genome-wide association studies on related individuals , 2015, Bioinform..

[39]  R. Yadav,et al.  FoxO transcription factors in cancer metabolism. , 2018, Seminars in cancer biology.

[40]  David Botstein,et al.  GO: : TermFinder--open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes , 2004, Bioinform..

[41]  Ahmad Samiei,et al.  Arachidonic acid as a target for treating hypertriglyceridemia reproduced by a causal network analysis and an intervention study , 2018, Metabolomics.

[42]  Nicola J. Rinaldi,et al.  Genetic effects on gene expression across human tissues , 2017, Nature.

[43]  Jianxin Shi,et al.  Common variants on chromosome 6p22.1 are associated with schizophrenia , 2009, Nature.

[44]  Peter Storz,et al.  Forkhead homeobox type O transcription factors in the responses to oxidative stress. , 2011, Antioxidants & redox signaling.

[45]  Eytan Domany,et al.  Comprehensive gene expression analysis detects global reduction of proteasome subunits in schizophrenia , 2020 .

[46]  George Davey Smith,et al.  Bayesian network analysis complements Mendelian randomization approaches for exploratory analysis of causal relationships in complex data , 2019, bioRxiv.