Abstract:Background Mapping expression quantitative trait loci (eQTLs) has provided insight into gene regulation. Compared to cis-eQTLs, the regulatory mechanisms of trans-eQTLs are less known. Previous studies suggest that trans-eQTLs may regulate expression of remote genes by altering the expression of nearby genes. Trans-association has been studied in the mediation analysis with a single mediator. However, prior applications with one mediator are prone to model misspecification due to correlations between genes. Motivated from the observation that trans-eQTLs are more likely to associate with more than one cis-gene than randomly selected SNPs in the GTEx dataset, we developed a computational method to identify trans-eQTLs that are mediated by multiple mediators. Results We proposed two hypothesis tests for testing the total mediation effect (TME) and the component-wise mediation effects (CME), respectively. We demonstrated in simulation studies that the type I error rates were controlled in both tests despite model misspecification. The TME test was more powerful than the CME test when the two mediation effects are in the same direction, while the CME test was more powerful than the TME test when the two mediation effects are in opposite direction. Multiple mediator analysis had increased power to detect mediated trans-eQTLs, especially in large samples. In the HapMap3 data, we identified 11 mediated trans-eQTLs that were not detected by the single mediator analysis in the combined samples of African populations. Moreover, the mediated trans-eQTLs in the HapMap3 samples are more likely to be trait-associated SNPs. In terms of computation, although there is no limit in the number of mediators in our model, analysis takes more time when adding additional mediators. In the analysis of the HapMap3 samples, we included at most 5 cis-gene mediators. Majority of the trios we considered have one or two mediators. Conclusions Trans-eQTLs are more likely to associate with multiple cis-genes than randomly selected SNPs. Mediation analysis with multiple mediators improves power of identification of mediated trans-eQTLs, especially in large samples.

A mediation analysis framework based on variance component to remove genetic confounding effect

Mediation Analysis with Mendelian Randomization and Efficient Multiple GWAS Integration

Generalized Multi-Snp Mediation Intersection-Union Test.

Dissecting the colocalized GWAS and eQTLs with mediation analysis for high-dimensional exposures and confounders

A novel statistical framework for meta-analysis of total mediation effect with high-dimensional omics mediators in large-scale genomic consortia

Estimation of total mediation effect for high-dimensional omics mediators

A Framework for Mediation Analysis with Multiple Exposures, Multivariate Mediators, and Non-Linear Response Models

Mediation analysis in longitudinal study with high-dimensional methylation mediators

Causal Mediation Analysis Leveraging Multiple Types of Summary Statistics Data

Speeding up interval estimation for R2-based mediation effect of high-dimensional mediators via cross-fitting

Joint analysis of SNP and gene expression data in genetic association studies of complex diseases

A new robust and accurate two-sample Mendelian randomization method with a large number of genetic variants

Causal Mediation Analysis for Integrating Exposure, Genomic, and Phenotype Data

Non-linear Mediation Analysis with High-dimensional Mediators whose Causal Structure is Unknown

Identification of trans-eQTLs using mediation analysis with multiple mediators

Causal Mediation Analysis in the Context of Clinical Research.

Mediation pathway selection with unmeasured mediator-outcome confounding

New proposal to address mediation analysis interrogations by using genetic variants as instrumental variables

Mediation Analysis Demonstrates That Trans-eQTLs Are Often Explained by Cis-Mediation: A Genome-Wide Analysis among 1,800 South Asians

Estimating and Testing High-Dimensional Mediation Effects in Epigenetic Studies

A Novel Method for Mendelian Randomization Analyses With Pleiotropy and Linkage Disequilibrium in Genetic Variants From Individual Data