Abstract:: The impact of genetic variants on gene expression has been intensely studied at the transcription level, yielding invaluable insights into the association between genes and the risk of complex disorders, such as schizophrenia (SCZ). However, the downstream impact of these variants and the molecular mechanisms connecting transcription variation to disease risk are not well understood. : We quantitated ribosome occupancy in prefrontal cortex samples of the BrainGVEX cohort. Together with transcriptomics and proteomics data from the same cohort, we performed cis-Quantitative Trait Locus (QTL) mapping and identified 3,253 expression QTLs (eQTLs), 1,344 ribosome occupancy QTLs (rQTLs), and 657 protein QTLs (pQTLs) out of 7,458 genes from 185 samples. Of the eQTLs identified, only 34% have their effects propagated to the protein level. Further analysis on the effect size of prefrontal cortex eQTLs identified from an independent dataset clearly replicated the post-transcriptional attenuation of eQTL effects. We identified omics-specific QTLs and investigated their potential in driving disease risks. Using a variant based approach, we found expression-specific QTLs (esQTLs) for 1,553 genes, ribosome-occupancy-specific QTLs (rsQTLs) for 155 genes, and protein-specific QTLs (psQTLs) for 161 genes. Among these omics-specific QTL, 38 showed strong colocalization with brain associated disorder GWAS signals, 29 of them are esQTLs. Because a gene could contain multiple QTL signals, each could either be shared across omics or omics-specific, we aggregated QTL signals from each omics for each gene and found 11 brain associated disorder risk genes that are driven predominantly by omics-specific QTL, all of them are driven by variants impacting transcriptional regulation. This gene-based approach also enabled us to categorize risk genes containing both omics-specific and shared QTL signals. The limited number of GWAS colocalization discoveries from gene-based omics-specific mapping, however, prompted us to take a complementary approach to investigate the functional relevance of genes driven predominantly by attenuated eQTL signals. Using S-PrediXcan we identified 74 SCZ risk genes across the three omics, 30% of which were novel, and 67% of these risk genes were confirmed to be causal in a MR-Egger test using data from the corresponding omics. Notably, 52 out of the 74 risk genes were identified using eQTL data and 68% of these SCZ-risk-gene-driving eQTLs show little to no evidence of driving corresponding variations at the protein level. : The effect of eQTLs on gene expression in the prefrontal cortex is commonly attenuated post-transcriptionally. Many of the attenuated eQTLs still correlate with GWAS signals of brain associated complex disorders, indicating the possibility that these eQTL variants drive disease risk through mechanisms other than regulating protein expression level. Further investigation is needed to elucidate the mechanistic link between attenuated eQTLs and brain associated complex disorders.

Brain transcriptome-wide association study implicates novel risk genes underlying schizophrenia risk

Joint-Tissue Integrative Analysis Identified Hundreds of Schizophrenia Risk Genes

The impact of common variants on gene expression in the human brain: from RNA to protein to schizophrenia risk

Common Variants on 17q25 and Gene–gene Interactions Conferring Risk of Schizophrenia in Han Chinese Population and Regulating Gene Expressions in Human Brain

Genome-wide association study followed by trans-ancestry meta-analysis identify 17 new risk loci for schizophrenia

Developmental and genetic regulation of the human cortex transcriptome illuminate schizophrenia pathogenesis

Brain eQTLs of European, African American, and Asian ancestry improve interpretation of schizophrenia GWAS

Transcriptome-wide Association Study Identifies New Susceptibility Genes and Pathways for Depression

Mendelian Randomization Study Using Dopaminergic Neuron-Specific eQTL Identifies Novel Risk Genes for Schizophrenia

Integration analysis of methylation quantitative trait loci and GWAS identify three schizophrenia risk variants

A genome-wide association study of attention and executive function in patients with schizophrenia

Novel genetic susceptibility loci identified by family based whole exome sequencing in Han Chinese schizophrenia patients

Integrated Transcriptome Analysis Reveals Novel Molecular Signatures for Schizophrenia Characterization.

Gene Expression Elucidates Functional Impact of Polygenic Risk for Schizophrenia

Comprehensive and integrative analyses identify TYW5 as a schizophrenia risk gene

Gene expression imputation across multiple brain regions provides insights into schizophrenia risk

Network-wide risk convergence in gene co-expression identifies reproducible genetic hubs of schizophrenia risk

Genetic architecture of RNA editing, splicing and gene expression in schizophrenia

Integrating Genome-Wide Association Study and Expression Quantitative Trait Locus Study Identifies Multiple Genes and Gene Sets Associated with Schizophrenia

Whole transcriptome analysis reveals dysregulation of molecular networks in schizophrenia

Regulatory variants at 2q33.1 confer schizophrenia risk by modulating distal gene TYW5 expression