A universal niche geometry governs the response of ecosystems to environmental perturbations.

How ecosystems respond to environmental perturbations is a fundamental question in ecology, made especially challenging due to the strong coupling between species and their environment. Here, we introduce a theoretical framework for calculating the linear response of ecosystems to environmental perturbations in generalized consumer-resource models. Our construction is applicable to a wide class of systems, including models with non-reciprocal interactions, cross-feeding, and non-linear growth/consumption rates. Within our framework, all ecological variables are embedded into four distinct vector spaces and ecological interactions are represented by geometric transformations between these spaces. We show that near a steady state, such geometric transformations directly map environmental perturbations - in resource availability and mortality rates - to shifts in niche structure. We illustrate these ideas in a variety of settings including a minimal model for pH-induced toxicity in bacterial denitrification.

The proteome is a terminal electron acceptor.

Microbial metabolism is impressively flexible, enabling growth even when available nutrients differ greatly from biomass in redox state. E. coli, for example, rearranges its physiology to grow on reduced and oxidized carbon sources through several forms of fermentation and respiration. To understand the limits on and evolutionary consequences of metabolic flexibility, we developed a mathematical model coupling redox chemistry with principles of cellular resource allocation. Our integrated model clarifies key phenomena, including demonstrating that autotrophs grow slower than heterotrophs because of constraints imposed by intracellular production of reduced carbon. Our model further indicates that growth is improved by adapting the redox state of biomass to nutrients, revealing an unexpected mode of evolution where proteins accumulate mutations benefiting organismal redox balance.

Closed ecosystems extract energy through self-organized nutrient cycles.

Our planet is a self-sustaining ecosystem powered by light energy from the sun, but roughly closed to matter. Many ecosystems on Earth are also approximately closed to matter and recycle nutrients by self-organizing stable nutrient cycles, e.g., microbial mats, lakes, open ocean gyres. However, existing ecological models do not exhibit the self-organization and dynamical stability widely observed in such planetary-scale ecosystems. Here, we advance a conceptual model that explains the self-organization, stability, and emergent features of closed microbial ecosystems. Our model incorporates the bioenergetics of metabolism into an ecological framework. By studying this model, we uncover a crucial thermodynamic feedback loop that enables metabolically diverse communities to almost always stabilize nutrient cycles. Surprisingly, highly diverse communities self-organize to extract [Formula: see text]10[Formula: see text] of the maximum extractable energy, or [Formula: see text]100 fold more than randomized communities. Further, with increasing diversity, distinct ecosystems show strongly correlated fluxes through nutrient cycles. However, as the driving force from light increases, the fluxes of nutrient cycles become more variable and species-dependent. Our results highlight that self-organization promotes the efficiency and stability of complex ecosystems at extracting energy from the environment, even in the absence of any centralized coordination.

Matching metabolic supply to demand optimizes microbial growth.

Recent research has strengthened the notion that microbes allocate their biosynthetic capacity to maximize the growth rate, λ. Yet many microbes can grow substantially faster after laboratory evolution. Chure and Cremer advance a resource-allocation model, which they derive from first principles, that offers resolution to this conundrum.

Annotation-free discovery of functional groups in microbial communities.

Recent studies have shown that microbial communities are composed of groups of functionally cohesive taxa whose abundance is more stable and better-associated with metabolic fluxes than that of any individual taxon. However, identifying these functional groups in a manner that is independent of error-prone functional gene annotations remains a major open problem. Here we tackle this structure-function problem by developing a novel unsupervised approach that coarse-grains taxa into functional groups, solely on the basis of the patterns of statistical variation in species abundances and functional read-outs. We demonstrate the power of this approach on three distinct datasets. On data of replicate microcosms with heterotrophic soil bacteria, our unsupervised algorithm recovered experimentally validated functional groups that divide metabolic labour and remain stable despite large variation in species composition. When leveraged against the ocean microbiome data, our approach discovered a functional group that combines aerobic and anaerobic ammonia oxidizers whose summed abundance tracks closely with nitrate concentrations in the water column. Finally, we show that our framework can enable the detection of species groups that are probably responsible for the production or consumption of metabolites abundant in animal gut microbiomes, serving as a hypothesis-generating tool for mechanistic studies. Overall, this work advances our understanding of structure-function relationships in complex microbiomes and provides a powerful approach to discover functional groups in an objective and systematic manner.

Community interactions drive the evolution of antibiotic tolerance in bacteria.

The emergence of antibiotic tolerance (prolonged survival against exposure) in natural bacterial populations is a major concern. Since it has been studied primarily in isogenic populations, we do not yet understand how ecological interactions in a diverse community impact the evolution of tolerance. To address this, we studied the evolutionary dynamics of a synthetic bacterial community composed of two interacting strains. In this community, an antibiotic-resistant strain protected the other, susceptible strain by degrading the antibiotic ampicillin in the medium. Surprisingly, we found that in the presence of antibiotics, the susceptible strain evolved tolerance. Tolerance was typified by an increase in survival as well as an accompanying decrease in the growth rate, highlighting a trade-off between the two. A simple mathematical model explained that the observed decrease in the death rate, even when coupled with a decreased growth rate, is beneficial in a community with weak protective interactions. In the presence of strong interactions, the model predicted that the trade-off would instead be detrimental, and tolerance would not emerge, which we experimentally verified. By whole genome sequencing the evolved tolerant isolates, we identified two genetic hot spots which accumulated mutations in parallel lines, suggesting their association with tolerance. Our work highlights that ecological interactions can promote antibiotic tolerance in bacterial communities, which has remained understudied.

Fine-scale diversity of microbial communities due to satellite niches in boom and bust environments.

Recent observations have revealed that closely related strains of the same microbial species can stably coexist in natural and laboratory settings subject to boom and bust dynamics and serial dilutions, respectively. However, the possible mechanisms enabling the coexistence of only a handful of strains, but not more, have thus far remained unknown. Here, using a consumer-resource model of microbial ecosystems, we propose that by differentiating along Monod parameters characterizing microbial growth rates in high and low nutrient conditions, strains can coexist in patterns similar to those observed. In our model, boom and bust environments create satellite niches due to resource concentrations varying in time. These satellite niches can be occupied by closely related strains, thereby enabling their coexistence. We demonstrate that this result is valid even in complex environments consisting of multiple resources and species. In these complex communities, each species partitions resources differently and creates separate sets of satellite niches for their own strains. While there is no theoretical limit to the number of coexisting strains, in our simulations, we always find between 1 and 3 strains coexisting, consistent with known experiments and observations.

Horizontal gene transfer drives the evolution of dependencies in bacteria.

Many naturally occurring bacteria lead a lifestyle of metabolic dependency for crucial resources. We do not understand what factors drive bacteria toward this lifestyle and how. Here, we systematically show the crucial role of horizontal gene transfer (HGT) in dependency evolution in bacteria. Across 835 bacterial species, we map gene gain-loss dynamics on a deep evolutionary tree and assess the impact of HGT and gene loss on metabolic networks. Our analyses suggest that HGT-enabled gene gains can affect which genes are later lost. HGT typically adds new catabolic routes to bacterial metabolic networks, leading to new metabolic interactions between bacteria. We also find that gaining new routes can promote the loss of ancestral routes ("coupled gains and losses", CGLs). Phylogenetic patterns indicate that both dependencies-mediated by CGLs and those purely by gene loss-are equally likely. Our results highlight HGT as an important driver of metabolic dependency evolution in bacteria.

Particle foraging strategies promote microbial diversity in marine environments.

Microbial foraging in patchy environments, where resources are fragmented into particles or pockets embedded in a large matrix, plays a key role in natural environments. In the oceans and freshwater systems, particle-associated bacteria can interact with particle surfaces in different ways: some colonize only during short transients, while others form long-lived, stable colonies. We do not yet understand the ecological mechanisms by which both short- and long-term colonizers can coexist. Here, we address this problem with a mathematical model that explains how marine populations with different detachment rates from particles can stably coexist. In our model, populations grow only while on particles, but also face the increased risk of mortality by predation and sinking. Key to coexistence is the idea that detachment from particles modulates both net growth and mortality, but in opposite directions, creating a trade-off between them. While slow-detaching populations show the highest growth return (i.e., produce more net offspring), they are more susceptible to suffer higher rates of mortality than fast-detaching populations. Surprisingly, fluctuating environments, manifesting as blooms of particles (favoring growth) and predators (favoring mortality) significantly expand the likelihood that populations with different detachment rates can coexist. Our study shows how the spatial ecology of microbes in the ocean can lead to a predictable diversification of foraging strategies and the coexistence of multiple taxa on a single growth-limiting resource.

Interactions between strains govern the eco-evolutionary dynamics of microbial communities.

Genomic data has revealed that genotypic variants of the same species, that is, strains, coexist and are abundant in natural microbial communities. However, it is not clear if strains are ecologically equivalent, and at what characteristic genetic distance they might exhibit distinct interactions and dynamics. Here, we address this problem by tracking 10 taxonomically diverse microbial communities from the pitcher plant Sarracenia purpurea in the laboratory for more than 300 generations. Using metagenomic sequencing, we reconstruct their dynamics over time and across scales, from distant phyla to closely related genotypes. We find that most strains are not ecologically equivalent and exhibit distinct dynamical patterns, often being significantly more correlated with strains from another species than their own. Although even a single mutation can affect laboratory strains, on average, natural strains typically decouple in their dynamics beyond a genetic distance of 100 base pairs. Using mathematical consumer-resource models, we show that these taxonomic patterns emerge naturally from ecological interactions between community members, but only if the interactions are coarse-grained at the level of strains, not species. Finally, by analyzing genomic differences between strains, we identify major functional hubs such as transporters, regulators, and carbohydrate-catabolizing enzymes, which might be the basis for strain-specific interactions. Our work suggests that fine-scale genetic differences in natural communities could be created and stabilized via the rapid diversification of ecological interactions between strains.

Complementary resource preferences spontaneously emerge in diauxic microbial communities.

Many microbes grow diauxically, utilizing the available resources one at a time rather than simultaneously. The properties of communities of microbes growing diauxically remain poorly understood, largely due to a lack of theory and models of such communities. Here, we develop and study a minimal model of diauxic microbial communities assembling in a serially diluted culture. We find that unlike co-utilizing communities, diauxic community assembly repeatably and spontaneously leads to communities with complementary resource preferences, namely communities where species prefer different resources as their top choice. Simulations and theory explain that the emergence of complementarity is driven by the disproportionate contribution of the top choice resource to the growth of a diauxic species. Additionally, we develop a geometric approach for analyzing serially diluted communities, with or without diauxie, which intuitively explains several additional emergent community properties, such as the apparent lack of species which grow fastest on a resource other than their most preferred resource. Overall, our work provides testable predictions for the assembly of natural as well as synthetic communities of diauxically shifting microbes.

Resource-diversity relationships in bacterial communities reflect the network structure of microbial metabolism.

The relationship between the number of available nutrients and community diversity is a central question in ecological research that remains unanswered. Here we studied the assembly of hundreds of soil-derived microbial communities on a wide range of well-defined resource environments, from single carbon sources to combinations of up to 16. We found that, while single resources supported multispecies communities varying from 8 to 40 taxa, mean community richness increased only one-by-one with additional resources. Cross-feeding could reconcile these seemingly contrasting observations, with the metabolic network seeded by the supplied resources explaining the changes in richness due to both the identity and the number of resources, as well as the distribution of taxa across different communities. By using a consumer-resource model incorporating the inferred cross-feeding network, we provide further theoretical support to our observations and a framework to link the type and number of environmental resources to microbial community diversity.

Ecology-guided prediction of cross-feeding interactions in the human gut microbiome.

Understanding a complex microbial ecosystem such as the human gut microbiome requires information about both microbial species and the metabolites they produce and secrete. These metabolites are exchanged via a large network of cross-feeding interactions, and are crucial for predicting the functional state of the microbiome. However, till date, we only have information for a part of this network, limited by experimental throughput. Here, we propose an ecology-based computational method, GutCP, using which we predict hundreds of new experimentally untested cross-feeding interactions in the human gut microbiome. GutCP utilizes a mechanistic model of the gut microbiome with the explicit exchange of metabolites and their effects on the growth of microbial species. To build GutCP, we combine metagenomic and metabolomic measurements from the gut microbiome with optimization techniques from machine learning. Close to 65% of the cross-feeding interactions predicted by GutCP are supported by evidence from genome annotations, which we provide for experimental testing. Our method has the potential to greatly improve existing models of the human gut microbiome, as well as our ability to predict the metabolic profile of the gut.

Evidence for a multi-level trophic organization of the human gut microbiome.

The human gut microbiome is a complex ecosystem, in which hundreds of microbial species and metabolites coexist, in part due to an extensive network of cross-feeding interactions. However, both the large-scale trophic organization of this ecosystem, and its effects on the underlying metabolic flow, remain unexplored. Here, using a simplified model, we provide quantitative support for a multi-level trophic organization of the human gut microbiome, where microbes consume and secrete metabolites in multiple iterative steps. Using a manually-curated set of metabolic interactions between microbes, our model suggests about four trophic levels, each characterized by a high level-to-level metabolic transfer of byproducts. It also quantitatively predicts the typical metabolic environment of the gut (fecal metabolome) in approximate agreement with the real data. To understand the consequences of this trophic organization, we quantify the metabolic flow and biomass distribution, and explore patterns of microbial and metabolic diversity in different levels. The hierarchical trophic organization suggested by our model can help mechanistically establish causal links between the abundances of microbes and metabolites in the human gut.

Metabolic adaptations underlying genome flexibility in prokaryotes.

Even across genomes of the same species, prokaryotes exhibit remarkable flexibility in gene content. We do not know whether this flexible or "accessory" content is mostly neutral or adaptive, largely due to the lack of explicit analyses of accessory gene function. Here, across 96 diverse prokaryotic species, I show that a considerable fraction (~40%) of accessory genomes harbours beneficial metabolic functions. These functions take two forms: (1) they significantly expand the biosynthetic potential of individual strains, and (2) they help reduce strain-specific metabolic auxotrophies via intra-species metabolic exchanges. I find that the potential of both these functions increases with increasing genome flexibility. Together, these results are consistent with a significant adaptive role for prokaryotic pangenomes.

Multiple stable states in microbial communities explained by the stable marriage problem.

Experimental studies of microbial communities routinely reveal that they have multiple stable states. While each of these states is generally resilient, certain perturbations such as antibiotics, probiotics, and diet shifts, result in transitions to other states. Can we reliably both predict such stable states as well as direct and control transitions between them? Here we present a new conceptual model-inspired by the stable marriage problem in game theory and economics-in which microbial communities naturally exhibit multiple stable states, each state with a different species' abundance profile. Our model's core ingredient is that microbes utilize nutrients one at a time while competing with each other. Using only two ranked tables, one with microbes' nutrient preferences and one with their competitive abilities, we can determine all possible stable states as well as predict inter-state transitions, triggered by the removal or addition of a specific nutrient or microbe. Further, using an example of seven Bacteroides species common to the human gut utilizing nine polysaccharides, we predict that mutual complementarity in nutrient preferences enables these species to coexist at high abundances.

Diversity, Stability, and Reproducibility in Stochastically Assembled Microbial Ecosystems.

Microbial ecosystems are remarkably diverse, stable, and usually consist of a mixture of core and peripheral species. Here we propose a conceptual model exhibiting all these emergent properties in quantitative agreement with real ecosystem data, specifically species abundance and prevalence distributions. Resource competition and metabolic commensalism drive the stochastic ecosystem assembly in our model. We demonstrate that even when supplied with just one resource, ecosystems can exhibit high diversity, increasing stability, and partial reproducibility between samples.

In silico Evolution of Lysis-Lysogeny Strategies Reproduces Observed Lysogeny Propensities in Temperate Bacteriophages.

Bacteriophages are the most abundant organisms on the planet and both lytic and temperate phages play key roles as shapers of ecosystems and drivers of bacterial evolution. Temperate phages can choose between (i) lysis: exploiting their bacterial hosts by producing multiple phage particles and releasing them by lysing the host cell, and (ii) lysogeny: establishing a potentially mutually beneficial relationship with the host by integrating their chromosome into the host cell's genome. Temperate phages exhibit lysogeny propensities in the curiously narrow range of 5-15%. For some temperate phages, the propensity is further regulated by the multiplicity of infection, such that single infections go predominantly lytic while multiple infections go predominantly lysogenic. We ask whether these observations can be explained by selection pressures in environments where multiple phage variants compete for the same host. Our models of pairwise competition, between phage variants that differ only in their propensity to lysogenize, predict the optimal lysogeny propensity to fall within the experimentally observed range. This prediction is robust to large variation in parameters such as the phage infection rate, burst size, decision rate, as well as bacterial growth rate, and initial phage to bacteria ratio. When we compete phage variants whose lysogeny strategies are allowed to depend upon multiplicity of infection, we find that the optimal strategy is one which switches from full lysis for single infections to full lysogeny for multiple infections. Previous attempts to explain lysogeny propensity have argued for bet-hedging that optimizes the response to fluctuating environmental conditions. Our results suggest that there is an additional selection pressure for lysogeny propensity within phage populations infecting a bacterial host, independent of environmental conditions.