Abstract:Expression quantitative trait loci (eQTL) studies are used to understand the regulatory function of non-coding genome-wide association study (GWAS) risk loci, but colocalization alone does not demonstrate a causal relationship of gene expression affecting a trait. Evidence for mediation, that perturbation of gene expression in a given tissue or developmental context will induce a change in the downstream GWAS trait, can be provided by two-sample Mendelian Randomization (MR). Here, we introduce a new statistical method, MRLocus, for Bayesian estimation of the gene-to-trait effect from eQTL and GWAS summary data for loci with evidence of allelic heterogeneity, that is, containing multiple causal variants. MRLocus makes use of a colocalization step applied to each nearly-LD-independent eQTL, followed by an MR analysis step across eQTLs. Additionally, our method involves estimation of the extent of allelic heterogeneity through a dispersion parameter, indicating variable mediation effects from each individual eQTL on the downstream trait. Our method is evaluated against other state-of-the-art methods for estimation of the gene-to-trait mediation effect, using an existing simulation framework. In simulation, MRLocus often has the highest accuracy among competing methods, and in each case provides more accurate estimation of uncertainty as assessed through interval coverage. MRLocus is then applied to five candidate causal genes for mediation of particular GWAS traits, where gene-to-trait effects are concordant with those previously reported. We find that MRLocus's estimation of the causal effect across eQTLs within a locus provides useful information for determining how perturbation of gene expression or individual regulatory elements will affect downstream traits. The MRLocus method is implemented as an R package available at . Author summary Genome-wide association studies (GWAS) have identified many loci associated with complex traits and diseases. Expression quantitative trait loci (eQTL) may help to explain mechanisms of GWAS associations, if the gene has a role as a mediator of the trait or disease. Loci that exhibit allelic heterogeneity, that is, loci containing multiple causal variants, offer the opportunity to investigate whether effects are concordant and proportional across eQTL and GWAS; if the gene is a partial mediator of the trait, the sign and size of the effects across distinct eQTL variants should be reflected in GWAS associations. Such a Mendelian Randomization (MR) analysis of individual loci is complicated by moderate sample sizes in eQTL studies and linkage disequilibrium (LD), resulting in complex patterns of estimated effect sizes for eQTL and GWAS. We develop a statistical model, MRLocus, with two steps: selection of eQTL SNPs to act as instruments in the MR analysis of a genetic locus, and estimation of the gene-to-trait mediation effect taking instrument uncertainty into account. In simulation, the method has higher accuracy and better uncertainty measures compared to other competing methods, and we compare its estimates on candidate causal gene-trait pairs from literature.

Mendelian randomization integrating GWAS and eQTL data reveals genetic determinants of complex and clinical traits

A robust two‐sample transcriptome‐wide Mendelian randomization method integrating GWAS with multi‐tissue eQTL summary statistics

A robust two-sample Mendelian Randomization method integrating GWAS with multi-tissue eQTL summary statistics

SMR-Portal: an online platform for integrative analysis of GWAS and xQTL data to identify complex trait genes

Multi-trait Transcriptome-Wide Association Studies with Probabilistic Mendelian Randomization

Unraveling the polygenic architecture of complex traits using blood eQTL metaanalysis

Powerful and robust inference of complex phenotypes' causal genes with dependent expression quantitative loci by a median-based Mendelian randomization.

Systematic identification of trans eQTLs as putative drivers of known disease associations

Epigenome-wide association study and integrative analysis with the transcriptome based on GWAS summary statistics

Bayesian Mendelian Randomization Analysis for Latent Exposures Leveraging GWAS Summary Statistics for Traits Co-Regulated by the Exposures

Testing and Controlling for Horizontal Pleiotropy with Probabilistic Mendelian Randomization in Transcriptome-Wide Association Studies

A Bayesian Framework for Multiple Trait Colocalization from Summary Association Statistics

Multi-omic QTL mapping in early developmental tissues reveals phenotypic and temporal complexity of regulatory variants underlying GWAS loci

Systematic comparative analysis of Mendelian randomization methods for inferring causal genes of complex phenotypes and the application to psychiatric diseases

An integrative multi-context Mendelian randomization method for identifying risk genes across human tissues

Large-scale cis- and trans-eQTL analyses identify thousands of genetic loci and polygenic scores that regulate blood gene expression

Adjusting for Genetic Confounders in Transcriptome-Wide Association Studies Leads to Reliable Detection of Causal Genes

Mendelian randomization for causal inference accounting for pleiotropy and sample structure using genome-wide summary statistics

Integrative analysis of omics summary data reveals putative mechanisms underlying complex traits

Adjusting for genetic confounders in transcriptome-wide association studies improves discovery of risk genes of complex traits

MRLocus: Identifying causal genes mediating a trait through Bayesian estimation of allelic heterogeneity