An analysis of genetically regulated gene expression and the role of co-expression networks across 16 psychiatric and substance use phenotypes.

Genome-wide association studies (GWASs) have identified thousands of risk loci for psychiatric and substance use phenotypes, however the biological consequences of these loci remain largely unknown. We performed a transcriptome-wide association study of 10 psychiatric disorders and 6 substance use phenotypes (GWAS sample size range, N = 9725-807,553) using expression quantitative trait loci data from 532 prefrontal cortex samples. We estimated the correlation of genetically regulated expression between phenotype pairs, and compared the results with the genetic correlations. We identified 393 genes with at least one significant phenotype association, comprising 458 significant associations across 16 phenotypes. Overall, the transcriptomic correlations for phenotype pairs were significantly higher than the respective genetic correlations. For example, attention deficit hyperactivity disorder and autism spectrum disorder, both childhood developmental disorders, had significantly higher transcriptomic correlation (r = 0.84) than genetic correlation (r = 0.35). Finally, we tested the enrichment of phenotype-associated genes in gene co-expression networks built from human prefrontal cortex samples. Phenotype-associated genes were enriched in multiple gene co-expression modules and the implicated modules contained genes involved in mRNA splicing and glutamatergic receptors, among others. Together, our results highlight the utility of gene expression data in the understanding of functional gene mechanisms underlying psychiatric disorders and substance use phenotypes.