Abstract:Abstract Background: Understanding whether and which microbes played a mediating role between an exposure and a disease outcome are essential for researchers to develop clinical interventions to treat the disease by modulating the microbes. Existing methods for mediation analysis of the microbiome are often limited to a global test of community-level mediation or selection of mediating microbes without control of the false discovery rate (FDR). Further, while the null hypothesis of no mediation at each microbe is a composite null that consists of three types of null (no exposure-microbe association, no microbe-outcome association given the exposure, or neither), most existing methods for the global test such as MedTest and MODIMA treat the microbes as if they are all under the same type of null. Results: We propose a new approach based on inverse regression that regresses the (possibly transformed) relative abundance of each taxon on the exposure and the exposure-adjusted outcome to assess the exposure-taxon and taxon-outcome associations simultaneously. Then the association p-values are used to test mediation at both the community and individual taxon levels. This approach fits nicely into our Linear Decomposition Model (LDM) framework, so our new method is implemented in the LDM and enjoys all the features of the LDM, i.e., allowing an arbitrary number of taxa to be tested, supporting continuous, discrete, or multivariate exposures and outcomes as well as adjustment of confounding covariates, accommodating clustered data, and offering analysis at the relative abundance or presence-absence scale. We refer to this new method as LDM-med. Using extensive simulations, we showed that LDM-med always controlled the type I error of the global test and had compelling power over existing methods; LDM-med always preserved the FDR of testing individual taxa and had much better sensitivity than alternative approaches. In contrast, MedTest and MODIMA had severely inflated type I error when different taxa were under different types of null. The flexibility of LDM-med for a variety of mediation analyses is illustrated by the application to a murine microbiome dataset, which identified a plausible mediator. Conclusions: Inverse regression coupled with the LDM is a strategy that performs well and is capable of handling mediation analysis in a wide variety of microbiome studies.

Mediation effect selection in high‐dimensional and compositional microbiome data

Estimating and testing the microbial causal mediation effect with high-dimensional and compositional microbiome data

Assessing mediating effects of high‐dimensional microbiome measurements in dietary intervention studies

HT-MMIOW: A Hypothesis Test approach for Microbiome Mediation using Inverse Odds Weighting

Testing Mediation Effects in High-Dimensional Microbiome Data with False Discovery Rate Control

Hypothesis-driven mediation analysis for compositional data: an application to gut microbiome

Inverse Probability Weighting-based Mediation Analysis for Microbiome Data

Estimation of total mediation effect for high-dimensional omics mediators

PhyloMed: a phylogeny-based test of mediation effect in microbiome

Causal Mediation Tree Model for Feature Identification on High-Dimensional Mediators

High-dimensional causal mediation analysis by partial sum statistic and sample splitting strategy in imaging genetics application

Mediation with External Summary Statistic Information (MESSI)

High-dimensional confounding in causal mediation: a comparison study of double machine learning and regularized partial correlation network

Mediation analysis in longitudinal study with high-dimensional methylation mediators

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

Regression Analysis for Microbiome Compositional Data

MASH: Mediation Analysis of Survival Outcome and High-dimensional Omics Mediators with Application to Complex Diseases

Multimodal data integration via mediation analysis with high‐dimensional exposures and mediators

Variable selection for sparse Dirichlet-multinomial regression with an application to microbiome data analysis

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

Comprehensive microbiome causal mediation analysis using MiMed on user-friendly web interfaces