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

Meta-analysis is used to aggregate the effects of interest across multiple studies, while its methodology is largely underexplored in mediation analysis, particularly in estimating the total mediation effect of high-dimensional omics mediators. Large-scale genomic consortia, such as the Trans-Omics for Precision Medicine (TOPMed) program, comprise multiple cohorts with diverse technologies to elucidate the genetic architecture and biological mechanisms underlying complex human traits and diseases. Leveraging the recent established asymptotic standard error of the R-squared (R2)-based mediation effect estimation for high-dimensional omics mediators, we have developed a novel meta-analysis framework requiring only summary statistics and allowing inter-study heterogeneity. Whereas the proposed meta-analysis can uniquely evaluate and account for potential effect heterogeneity across studies due to, for example, varying genomic profiling platforms, our extensive simulations showed that the developed method was more computationally efficient and yielded satisfactory operating characteristics comparable to analysis of the pooled individual-level data when there was no inter-study heterogeneity. We applied the developed method to 5 TOPMed studies with over 5800 participants to estimate the mediation effects of gene expression on age-related variation in systolic blood pressure and sex-related variation in high-density lipoprotein (HDL) cholesterol. The proposed method is available in R package MetaR2M on GitHub.