Robust Microbial Signature Discovery via Post-Selection Inference for Microbiome Compositions
提出PoDA方法,通过后选择推断解决微生物组数据的高维、组成性和偏差问题,在帕金森病研究中识别出可复现的微生物特征。
Identifying taxa associated with host phenotypes is crucial for understanding host-microbe interactions and their underlying molecular mechanisms. However, analyzing microbiome data presents unique challenges, as the observed abundances of taxa are high-dimensional, compositional, and subject to both sample-specific and taxon-specific biases. Many existing methods for differential abundance testing struggle to balance false discovery rate control with statistical power. In this paper, we propose PoDA, a post-selection inference method for differential abundance analysis to address the limitations of the existing methods. PoDA begins by selecting a subset of taxa likely associated with the phenotype using penalized regression under a mean-shift model. It then leverages the unselected taxa to correct for sample-specific bias and assign p-values to the selected taxa. To ensure valid inference after the selection process, PoDA employs an information-splitting procedure, which is repeated to enhance stability and power. Comprehensive simulation studies demonstrate the superiority of PoDA over existing methods. We further applied PoDA to two microbiome case-control studies of Parkinson’s disease (PD). The method identified a set of candidate microbial signatures associated with PD and showed improved replicability across the two independent datasets. PoDA is implemented in R and available at https://github.com/weihaowang01/PoDA.