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.

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

Mediation effect selection in high‐dimensional and compositional microbiome data

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

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

Causal mediation analysis in presence of multiple mediators uncausally related

Inverse Probability Weighting-based Mediation Analysis for Microbiome Data

Bayesian Sparse Mediation Analysis with Targeted Penalization of Natural Indirect Effects

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

Latent Class Dynamic Mediation Model with Application to Smoking Cessation Data

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

Instrumental Variable Estimation for Compositional Treatments

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

Mediation with External Summary Statistic Information (MESSI)

Decomposition of the Total Effect for Two Mediators: A Natural Counterfactual Interaction Effect Framework

Mediation Analysis for Count and Zero-Inflated Count Data without Sequential Ignorability and Its Application in Dental Studies

Estimation of total mediation effect for high-dimensional omics mediators

Causal mediation analysis for stochastic interventions

Regression Analysis for Microbiome Compositional Data

Population‐average mediation analysis for zero‐inflated count outcomes

Mediation Analysis using Semi-parametric Shape-Restricted Regression with Applications