The interspecies exchange of metabolites plays a key role in the spatiotemporal dynamics of microbial communities. This raises the question of whether ecosystem-level behavior of structured communities can be predicted using genome-scale metabolic models for multiple organisms. We developed a modeling framework that integrates dynamic flux balance analysis with diffusion on a lattice and applied it to engineered communities. First, we predicted and experimentally confirmed the species ratio to which a two-species mutualistic consortium converges and the equilibrium composition of a newly engineered three-member community. We next identified a specific spatial arrangement of colonies, which gives rise to what we term the "eclipse dilemma": does a competitor placed between a colony and its cross-feeding partner benefit or hurt growth of the original colony? Our experimentally validated finding that the net outcome is beneficial highlights the complex nature of metabolic interactions in microbial communities while at the same time demonstrating their predictability.

The oral microbiome, the complex ecosystem of microbes inhabiting the human mouth, harbors several thousands of bacterial types. The proliferation of pathogenic bacteria within the mouth gives rise to periodontitis, an inflammatory disease known to also constitute a risk factor for cardiovascular disease. While much is known about individual species associated with pathogenesis, the system-level mechanisms underlying the transition from health to disease are still poorly understood. Through the sequencing of the 16S rRNA gene and of whole community DNA we provide a glimpse at the global genetic, metabolic, and ecological changes associated with periodontitis in 15 subgingival plaque samples, four from each of two periodontitis patients, and the remaining samples from three healthy individuals. We also demonstrate the power of whole-metagenome sequencing approaches in characterizing the genomes of key players in the oral microbiome, including an unculturable TM7 organism. We reveal the disease microbiome to be enriched in virulence factors, and adapted to a parasitic lifestyle that takes advantage of the disrupted host homeostasis. Furthermore, diseased samples share a common structure that was not found in completely healthy samples, suggesting that the disease state may occupy a narrow region within the space of possible configurations of the oral microbiome. Our pilot study demonstrates the power of high-throughput sequencing as a tool for understanding the role of the oral microbiome in periodontal disease. Despite a modest level of sequencing (∼2 lanes Illumina 76 bp PE) and high human DNA contamination (up to ∼90%) we were able to partially reconstruct several oral microbes and to preliminarily characterize some systems-level differences between the healthy and diseased oral microbiomes.

Metabolic exchange mediates interactions among microbes, helping explain diversity in microbial communities. As these interactions often involve a fitness cost, it is unclear how stable cooperation can emerge. Here we use genome-scale metabolic models to investigate whether the release of “costless” metabolites (i.e. those that cause no fitness cost to the producer), can be a prominent driver of intermicrobial interactions. By performing over 2 million pairwise growth simulations of 24 species in a combinatorial assortment of environments, we identify a large space of metabolites that can be secreted without cost, thus generating ample cross-feeding opportunities. In addition to providing an atlas of putative interactions, we show that anoxic conditions can promote mutualisms by providing more opportunities for exchange of costless metabolites, resulting in an overrepresentation of stable ecological network motifs. These results may help identify interaction patterns in natural communities and inform the design of synthetic microbial consortia.

Interactions between microbial species are sometimes mediated by the exchange of small molecules, secreted by one species and metabolized by another. Both one-way (commensal) and two-way (mutualistic) interactions may contribute to complex networks of interdependencies. Understanding these interactions constitutes an open challenge in microbial ecology, with applications ranging from the human microbiome to environmental sustainability. In parallel to natural communities, it is possible to explore interactions in artificial microbial ecosystems, e.g. pairs of genetically engineered mutualistic strains. Here we computationally generate artificial microbial ecosystems without re-engineering the microbes themselves, but rather by predicting their growth on appropriately designed media. We use genome-scale stoichiometric models of metabolism to identify media that can sustain growth for a pair of species, but fail to do so for one or both individual species, thereby inducing putative symbiotic interactions. We first tested our approach on two previously studied mutualistic pairs, and on a pair of highly curated model organisms, showing that our algorithms successfully recapitulate known interactions, robustly predict new ones, and provide novel insight on exchanged molecules. We then applied our method to all possible pairs of seven microbial species, and found that it is always possible to identify putative media that induce commensalism or mutualism. Our analysis also suggests that symbiotic interactions may arise more readily through environmental fluctuations than genetic modifications. We envision that our approach will help generate microbe-microbe interaction maps useful for understanding microbial consortia dynamics and evolution, and for exploring the full potential of natural metabolic pathways for metabolic engineering applications. Microbial metabolism affects biogeochemical cycles and human health. In most natural environments, multiple microbial species interact with each other, forming complex ecosystems whose properties are poorly understood. In an effort to understand inter-microbial interactions, and to explore new metabolic engineering avenues, researchers have started building artificial microbial ecosystems, e.g. pairs of genetically engineered strains that require each other for survival. Here we computationally explore the possibility of creating artificial microbial ecosystems without re-engineering the microbes themselves, but rather by manipulating the environment in which they grow. Specifically, using the framework of flux balance analysis, we predict environments in which either one or both microbes in a pair would not be able to grow without the other, inducing commensal (one-way) or mutualistic (two-way) interactions, respectively. Our algorithms can successfully recapitulate known inter-microbial interactions, and predict millions of new ones across any pair amongst different microbial species. Surprisingly, we find that it is always possible to identify conditions that induce mutualistic or commensal interactions between any two species. Hence, our method should help in mapping naturally occurring microbe-microbe interactions, and in engineering new ones through a novel, environment-driven branch of synthetic ecology.