Understanding the link between community composition and function is a major challenge in microbial population biology, with implications for the management of natural microbiomes and the design of synthetic consortia. Specifically, it is poorly understood whether community functions can be quantitatively predicted from traits of species in monoculture. Inspired by the study of complex genetic interactions, we have examined how the amylolytic rate of combinatorial assemblages of six starch-degrading soil bacteria depend on the separate functional contributions from each species and their interactions. Filtering our results through the theory of biochemical kinetics, we show that this simple function is additive in the absence of interactions among community members. For about half of the combinatorially assembled consortia, the amylolytic function is dominated by pairwise and higher-order interactions. For the other half, the function is additive despite the presence of strong competitive interactions. We explain the mechanistic basis of these findings and propose a quantitative framework that allows us to separate the effect of behavioral and population dynamics interactions. Our results suggest that the functional robustness of a consortium to pairwise and higher-order interactions critically affects our ability to predict and bottom-up engineer ecosystem function in complex communities.

Artificial selection is a promising approach to manipulate microbial communities. Here, we report the outcome of two artificial selection experiments at the microbial community level. Both used "propagule" selection strategies, whereby the best-performing communities are used as the inocula to form a new generation of communities. Both experiments were contrasted to a random selection control. The first experiment used a defined set of strains as the starting inoculum, and the function under selection was the amylolytic activity of the consortia. The second experiment used multiple soil communities as the starting inocula, and the function under selection was the communities' cross-feeding potential. In both experiments, the selected communities reached a higher mean function than the control. In the first experiment, this was caused by a decline in function of the control, rather than an improvement of the selected line. In the second experiment, this response was fueled by the large initial variance in function across communities, and stopped when the top-performing community "fixed" in the metacommunity. Our results are in agreement with basic expectations from breeding theory, pointing to some of the limitations of community-level selection experiments that can inform the design of future studies.

Francisella tularensis is a Gram-negative bacterium whose ability to replicate within macrophages and cause disease is strictly dependent upon the coordinate activities of three transcription regulators called MglA, SspA, and PigR. MglA and SspA form a complex that associates with RNA polymerase (RNAP), whereas PigR is a putative DNA-binding protein that functions by contacting the MglA-SspA complex. Most transcription activators that bind the DNA are thought to occupy only those promoters whose activities they regulate. Here we show using chromatin immunoprecipitation coupled with high-throughput DNA sequencing (ChIP-Seq) that PigR, MglA, and SspA are found at virtually all promoters in F. tularensis and not just those of regulated genes. Furthermore, we find that the ability of PigR to associate with promoters is dependent upon the presence of MglA, suggesting that interaction with the RNAP-associated MglA-SspA complex is what directs PigR to promoters in F. tularensis. Finally, we present evidence that the ability of PigR (and thus MglA and SspA) to positively control the expression of genes is dictated by a specific 7 base pair sequence element that is present in the promoters of regulated genes. The three principal regulators of virulence gene expression in F. tularensis therefore function in a non-classical manner with PigR interacting with the RNAP-associated MglA-SspA complex at the majority of promoters but only activating transcription from those that contain a specific sequence element. Our findings reveal how transcription factors can exert regulatory effects at a restricted set of promoters despite being associated with most or all. This distinction between occupancy and regulatory effect uncovered by our data may be relevant to the study of RNAP-associated transcription regulators in other pathogenic bacteria.