Microorganisms display a stunning metabolic diversity. Understanding the origin of this diversity requires understanding how macroevolutionary processes such as innovation and diversification play out in the microbial world. Metabolic networks, which govern microbial resource use, can evolve through different mechanisms, e.g., horizontal gene transfer or de novo evolution of enzymes and pathways. This process is governed by a combination of environmental factors, selective pressures, and the constraints imposed by the genetic architecture of metabolic networks. In addition, many independent results hint that the process of niche construction, by which organisms actively modify their own and each other's niches and selective pressures, could play a major role in microbial innovation and diversification. Yet, the general principles by which niche construction shapes microbial macroevolutionary patterns remain largely unexplored. Here, we discuss several new hypotheses and directions, and suggest metabolic modeling methods that could allow us to explore large-scale empirical genotype-phenotype-(G-P)-environment spaces in order to study the macroevolutionary effects of niche construction. We hope that this short piece will further stimulate a systematic and quantitative characterization of macroevolutionary patterns and processes in microbial metabolism.

For microbiome biology to become a more predictive science, we must identify which descriptive features of microbial communities are reproducible and predictable, which are not, and why. We address this question by experimentally studying parallelism and convergence in microbial community assembly in replicate glucose-limited habitats. Here, we show that the previously observed family-level convergence in these habitats reflects a reproducible metabolic organization, where the ratio of the dominant metabolic groups can be explained from a simple resource-partitioning model. In turn, taxonomic divergence among replicate communities arises from multistability in population dynamics. Multistability can also lead to alternative functional states in closed ecosystems but not in metacommunities. Our findings empirically illustrate how the evolutionary conservation of quantitative metabolic traits, multistability, and the inherent stochasticity of population dynamics, may all conspire to generate the patterns of reproducibility and variability at different levels of organization that are commonplace in microbial community assembly.

Genome-scale stoichiometric modeling of metabolism has become a standard systems biology tool for modeling cellular physiology and growth. Extensions of this approach are emerging as a valuable avenue for predicting, understanding and designing microbial communities. Computation of microbial ecosystems in time and space (COMETS) extends dynamic flux balance analysis to generate simulations of multiple microbial species in molecularly complex and spatially structured environments. Here we describe how to best use and apply the most recent version of COMETS, which incorporates a more accurate biophysical model of microbial biomass expansion upon growth, evolutionary dynamics and extracellular enzyme activity modules. In addition to a command-line option, COMETS includes user-friendly Python and MATLAB interfaces compatible with the well-established COBRA models and methods, as well as comprehensive documentation and tutorials. This protocol provides a detailed guideline for installing, testing and applying COMETS to different scenarios, generating simulations that take from a few minutes to several days to run, with broad applicability to microbial communities across biomes and scales.