Integrating phylogenetic and functional data in microbiome studies

Abstract

Microbiome functional data are frequently analyzed to identify associations between microbial gene families and sample groups of interest. This is most often performed with approaches focused on the metagenome-wide relative abundance of microbial functions. Although such approaches can provide valuable insights, it is challenging to distinguish between different possible explanations for variation in community-wide functional profiles by considering functions alone. To help address this problem, we have developed a novel, phylogeny-aware framework to expand taxonomic balance tree-based approaches to identify enriched functions more robustly. The key focus of our approach, termed POMS, is on identifying functions that are consistently enriched in sample groups across independent taxonomic lineages. Based on simulated data we demonstrate that POMS can more accurately identify gene families that confer a selective advantage compared with commonly used differential abundance approaches. We also show that POMS can identify enriched functions in real-world metagenomics datasets that are potential targets of strong selection on multiple members of the microbiome. While this framework may not be able to identify all potential functional enrichments, the enrichments it does identify are more interpretable and conservative compared with those identified by existing differential abundance approaches. More generally, POMS is a novel approach for exploring microbiome functional data, which could be used to complement standard analyses. POMS is freely available as an R package at: https://github.com/gavinmdouglas/POMS.