The Contribution of High-Order Metabolic Interactions to the Global Activity of a Four-Species Microbial Community

The activity of a biological community is the outcome of complex processes involving interactions between community members. It is often unclear how to accurately incorporate these interactions into predictive models. Previous work has shown a range of positive and negative metabolic pairwise interactions between species. Here we examine the ability of a modified general Lotka-Volterra model with cell-cell interaction coefficients to predict the overall metabolic rate of a well-mixed microbial community comprised of four heterotrophic natural isolates, experimentally quantifying the strengths of two, three, and four-species interactions. Within this community, interactions between any pair of microbial species were positive, while higher-order interactions, between 3 or more microbial species, slightly modulated community metabolism. For this simple community, the metabolic rate of can be well predicted only with taking into account pairwise interactions. Simulations using the experimentally determined interaction parameters revealed that spatial heterogeneity in the distribution of cells increased the importance of multispecies interactions in dictating function at both the local and global scales.