Abstract
MOTIVATION: 16S rRNA gene sequencing is the most frequent approach for the characterization of the human gut microbiota. Despite different efforts in the literature, the inference of functional and metabolic interpretations from 16S rRNA gene sequencing data is still a challenging task. High-quality metabolic reconstructions of the human gut microbiota, such as AGORA and AGREDA, constitute a curated resource to improve functional inference from 16S rRNA data, but they are not typically integrated into standard bioinformatics tools. RESULTS: Here, we present q2-metnet, a QIIME2 plugin that enables the contextualization of 16S rRNA gene sequencing data into AGORA and AGREDA. In particular, based on relative abundances of taxa, q2-metnet determines normalized activity scores for the reactions and subsystems involved in the selected metabolic reconstruction. Using these scores, q2-metnet allows the user to conduct differential activity analysis for reactions and subsystems, as well as exploratory analysis using PCA and hierarchical clustering. We apply q2-metnet to a dataset from our group that involves 16S rRNA data from stool samples from lean, allergic to cow's milk, obese and celiac children, and the Belgian Flemish Gut Flora Project cohort, which includes faecal 16S rRNA data from obese and normal-weight adult individuals. In the first case, q2-metnet outperforms existing algorithms in separating different clinical conditions based on predicted pathway abundances and subsystem scores. In the second case, q2-metnet complements competing approaches in predicting functional alterations in the gut microbiota of obese individuals. Overall, q2-metnet constitutes a powerful bioinformatics tool to provide metabolic context to 16S rRNA data from the human gut microbiota. AVAILABILITY AND IMPLEMENTATION: Python code of q2-metnet is available in https://github.com/PlanesLab/q2-metnet and https://figshare.com/articles/dataset/q2-metnet_package/26180446.