Taxa-function robustness in microbial communities

Background: The species composition of a microbial community is rarely fixed and often experiences fluctuations of varying degrees and at varying frequencies. These perturbations to a community's taxonomic profile naturally also alter the community's functional profile-the aggregate set of genes encoded by community members-ultimately altering the community's overall functional capacities. The magnitude of such functional changes and the specific shift that will occur in each function, however, are strongly dependent on how genes are distributed across community members' genomes. This gene distribution, in turn, is determined by the taxonomic composition of the community and would markedly differ, for example, between communities composed of species with similar genomic content vs. communities composed of species whose genomes encode relatively distinct gene sets. Combined, these observations suggest that community functional robustness to taxonomic perturbations could vary widely across communities with different compositions, yet, to date, a systematic study of the inherent link between community composition and robustness is lacking. Results: In this study, we examined how a community's taxonomic composition influences the robustness of that community's functional profile to taxonomic perturbation (here termed taxa-function robustness) across a wide array of environments. Using a novel simulation-based computational model to quantify this taxa-function robustness in host-associated and non-host-associated communities, we find notable differences in robustness between communities inhabiting different body sites, including significantly higher robustness in gut communities compared to vaginal communities that cannot be attributed solely to differences in species richness. We additionally find between-site differences in the robustness of specific functions, some of which are potentially related to site-specific environmental conditions. These taxa-function robustness differences are most strongly associated with differences in overall functional redundancy, though other aspects of gene distribution also influence taxa-function robustness in certain body environments, and are sufficient to cluster communities by environment. Further analysis revealed a correspondence between our robustness estimates and taxonomic and functional shifts observed across human-associated communities. Conclusions: Our analysis approach revealed intriguing taxa-function robustness variation across environments and identified features of community and gene distribution that impact robustness. This approach could be further applied for estimating taxa-function robustness in novel communities and for informing the design of synthetic communities with specific robustness requirements.