Bayesian Genome-wide TWAS Method to Leverage both cis- and trans-eQTL Information through Summary Statistics

Transcriptome-wide association studies (TWASs) have been widely used to integrate gene expression and genetic data for studying complex traits. Due to the computational burden, existing TWAS methods do not assess distant trans-expression quantitative trait loci (eQTL) that are known to explain important expression variation for most genes. We propose a Bayesian genome-wide TWAS (BGW-TWAS) method that leverages both cis- and trans-eQTL information for a TWAS. Our BGW-TWAS method is based on Bayesian variable selection regression, which not only accounts for cis- and trans-eQTL of the target gene but also enables efficient computation by using summary statistics from standard eQTL analyses. Our simulation studies illustrated that BGW-TWASs achieved higher power compared to existing TWAS methods that do not assess trans-eQTL information. We further applied BWG-TWAS to individual-level GWAS data (N = ∼3.3K), which identified significant associations between the genetically regulated gene expression (GReX) of ZC3H12B and Alzheimer dementia (AD) (p value = 5.42 × 10⁻¹³), neurofibrillary tangle density (p value = 1.89 × 10⁻⁶), and global measure of AD pathology (p value = 9.59 × 10⁻⁷). These associations for ZC3H12B were completely driven by trans-eQTL. Additionally, the GReX of KCTD12 was found to be significantly associated with β-amyloid (p value = 3.44 × 10⁻⁸) which was driven by both cis- and trans-eQTL. Four of the top driven trans-eQTL of ZC3H12B are located within APOC1, a known major risk gene of AD and blood lipids. Additionally, by applying BGW-TWAS with summary-level GWAS data of AD (N = ∼54K), we identified 13 significant genes including known GWAS risk genes HLA-DRB1 and APOC1, as well as ZC3H12B.