Global epistasis and the emergence of function in microbial consortia.

The many functions of microbial communities emerge from a complex web of interactions between organisms and their environment. This poses a significant obstacle to engineering microbial consortia, hindering our ability to harness the potential of microorganisms for biotechnological applications. In this study, we demonstrate that the collective effect of ecological interactions between microbes in a community can be captured by simple statistical models that predict how adding a new species to a community will affect its function. These predictive models mirror the patterns of global epistasis reported in genetics, and they can be quantitatively interpreted in terms of pairwise interactions between community members. Our results illuminate an unexplored path to quantitatively predicting the function of microbial consortia from their composition, paving the way to optimizing desirable community properties and bringing the tasks of predicting biological function at the genetic, organismal, and ecological scales under the same quantitative formalism.

Metabolic similarity and the predictability of microbial community assembly.

When microbial communities form, their composition is shaped by selective pressures imposed by the environment. Can we predict which communities will assemble under different environmental conditions? Here, we hypothesize that quantitative similarities in metabolic traits across metabolically similar environments lead to predictable similarities in community composition. To that end, we measured the growth rate and by-product profile of a library of proteobacterial strains in a large number of single nutrient environments. We found that growth rates and secretion profiles were positively correlated across environments when the supplied substrate was metabolically similar. By analyzing hundreds of in-vitro communities experimentally assembled in an array of different synthetic environments, we then show that metabolically similar substrates select for taxonomically similar communities. These findings lead us to propose and then validate a comparative approach for quantitatively predicting the effects of novel substrates on the composition of complex microbial consortia.

The Architecture of Metabolic Networks Constrains the Evolution of Microbial Resource Hierarchies.

Microbial strategies for resource use are an essential determinant of their fitness in complex habitats. When facing environments with multiple nutrients, microbes often use them sequentially according to a preference hierarchy, resulting in well-known patterns of diauxic growth. In theory, the evolutionary diversification of metabolic hierarchies could represent a mechanism supporting coexistence and biodiversity by enabling temporal segregation of niches. Despite this ecologically critical role, the extent to which substrate preference hierarchies can evolve and diversify remains largely unexplored. Here, we used genome-scale metabolic modeling to systematically explore the evolution of metabolic hierarchies across a vast space of metabolic network genotypes. We find that only a limited number of metabolic hierarchies can readily evolve, corresponding to the most commonly observed hierarchies in genome-derived models. We further show how the evolution of novel hierarchies is constrained by the architecture of central metabolism, which determines both the propensity to change ranks between pairs of substrates and the effect of specific reactions on hierarchy evolution. Our analysis sheds light on the genetic and mechanistic determinants of microbial metabolic hierarchies, opening new research avenues to understand their evolution, evolvability, and ecology.

A structured evaluation of genome-scale constraint-based modeling tools for microbial consortia.

Harnessing the power of microbial consortia is integral to a diverse range of sectors, from healthcare to biotechnology to environmental remediation. To fully realize this potential, it is critical to understand the mechanisms behind the interactions that structure microbial consortia and determine their functions. Constraint-based reconstruction and analysis (COBRA) approaches, employing genome-scale metabolic models (GEMs), have emerged as the state-of-the-art tool to simulate the behavior of microbial communities from their constituent genomes. In the last decade, many tools have been developed that use COBRA approaches to simulate multi-species consortia, under either steady-state, dynamic, or spatiotemporally varying scenarios. Yet, these tools have not been systematically evaluated regarding their software quality, most suitable application, and predictive power. Hence, it is uncertain which tools users should apply to their system and what are the most urgent directions that developers should take in the future to improve existing capacities. This study conducted a systematic evaluation of COBRA-based tools for microbial communities using datasets from two-member communities as test cases. First, we performed a qualitative assessment in which we evaluated 24 published tools based on a list of FAIR (Findability, Accessibility, Interoperability, and Reusability) features essential for software quality. Next, we quantitatively tested the predictions in a subset of 14 of these tools against experimental data from three different case studies: a) syngas fermentation by C. autoethanogenum and C. kluyveri for the static tools, b) glucose/xylose fermentation with engineered E. coli and S. cerevisiae for the dynamic tools, and c) a Petri dish of E. coli and S. enterica for tools incorporating spatiotemporal variation. Our results show varying performance levels of the best qualitatively assessed tools when examining the different categories of tools. The differences in the mathematical formulation of the approaches and their relation to the results were also discussed. Ultimately, we provide recommendations for refining future GEM microbial modeling tools.

Emergent coexistence in multispecies microbial communities.

Understanding the mechanisms that maintain microbial biodiversity is a critical aspiration in ecology. Past work on microbial coexistence has largely focused on species pairs, but it is unclear whether pairwise coexistence in isolation is required for coexistence in a multispecies community. To address this question, we conducted hundreds of pairwise competition experiments among the stably coexisting members of 12 different enrichment communities in vitro. To determine the outcomes of these experiments, we developed an automated image analysis pipeline to quantify species abundances. We found that competitive exclusion was the most common outcome, and it was strongly hierarchical and transitive. Because many species that coexist within a stable multispecies community fail to coexist in pairwise co-culture under identical conditions, we concluded that multispecies coexistence is an emergent phenomenon. This work highlights the importance of community context for understanding the origins of coexistence in complex ecosystems.

Global epistasis on fitness landscapes.

Epistatic interactions between mutations add substantial complexity to adaptive landscapes and are often thought of as detrimental to our ability to predict evolution. Yet, patterns of global epistasis, in which the fitness effect of a mutation is well-predicted by the fitness of its genetic background, may actually be of help in our efforts to reconstruct fitness landscapes and infer adaptive trajectories. Microscopic interactions between mutations, or inherent nonlinearities in the fitness landscape, may cause global epistasis patterns to emerge. In this brief review, we provide a succinct overview of recent work about global epistasis, with an emphasis on building intuition about why it is often observed. To this end, we reconcile simple geometric reasoning with recent mathematical analyses, using these to explain why different mutations in an empirical landscape may exhibit different global epistasis patterns-ranging from diminishing to increasing returns. Finally, we highlight open questions and research directions. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

The community-function landscape of microbial consortia.

Quantitatively linking the composition and function of microbial communities is a major aspiration of microbial ecology. Microbial community functions emerge from a complex web of molecular interactions between cells, which give rise to population-level interactions among strains and species. Incorporating this complexity into predictive models is highly challenging. Inspired by a similar problem in genetics of predicting quantitative phenotypes from genotypes, an ecological community-function (or structure-function) landscape could be defined that maps community composition and function. In this piece, we present an overview of our current understanding of these community landscapes, their uses, limitations, and open questions. We argue that exploiting the parallels between both landscapes could bring powerful predictive methodologies from evolution and genetics into ecology, providing a boost to our ability to engineer and optimize microbial consortia.

Functional attractors in microbial community assembly.

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.

A metabolic modeling platform for the computation of microbial ecosystems in time and space (COMETS).

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.

The Macroevolutionary Consequences of Niche Construction in Microbial Metabolism.

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.

The Impact of Gut Microbiota-Derived Metabolites in Autism Spectrum Disorders.

Autism Spectrum Disorder (ASD) is a set of neurodevelopmental disorders characterised by behavioural impairment and deficiencies in social interaction and communication. A recent study estimated that 1 in 89 children have developed some form of ASD in European countries. Moreover, there is no specific treatment and since ASD is not a single clinical entity, the identification of molecular biomarkers for diagnosis remains challenging. Besides behavioural deficiencies, individuals with ASD often develop comorbid medical conditions including intestinal problems, which may reflect aberrations in the bidirectional communication between the brain and the gut. The impact of faecal microbial composition in brain development and behavioural functions has been repeatedly linked to ASD, as well as changes in the metabolic profile of individuals affected by ASD. Since metabolism is one of the major drivers of microbiome-host interactions, this review aims to report emerging literature showing shifts in gut microbiota metabolic function in ASD. Additionally, we discuss how these changes may be involved in and/or perpetuate ASD pathology. These valuable insights can help us to better comprehend ASD pathogenesis and may provide relevant biomarkers for improving diagnosis and identifying new therapeutic targets.

Tolerance to NADH/NAD+ imbalance anticipates aging and anti-aging interventions.

Redox couples coordinate cellular function, but the consequences of their imbalances are unclear. This is somewhat associated with the limitations of their experimental quantification. Here we circumvent these difficulties by presenting an approach that characterizes fitness-based tolerance profiles to redox couple imbalances using an in silico representation of metabolism. Focusing on the NADH/NAD+ redox couple in yeast, we demonstrate that reductive disequilibria generate metabolic syndromes comparable to those observed in cancer cells. The tolerance of yeast mutants to redox disequilibrium can also explain 30% of the variability in their experimentally measured chronological lifespan. Moreover, by predicting the significance of some metabolites to help stand imbalances, we correctly identify nutrients underlying mechanisms of pathology, lifespan-protecting molecules, or caloric restriction mimetics. Tolerance to redox imbalances becomes, in this way, a sound framework to recognize properties of the aging phenotype while providing a consistent biological rationale to assess anti-aging interventions.

Complex yeast-bacteria interactions affect the yield of industrial ethanol fermentation.

Sugarcane ethanol fermentation represents a simple microbial community dominated by S. cerevisiae and co-occurring bacteria with a clearly defined functionality. In this study, we dissect the microbial interactions in sugarcane ethanol fermentation by combinatorically reconstituting every possible combination of species, comprising approximately 80% of the biodiversity in terms of relative abundance. Functional landscape analysis shows that higher-order interactions counterbalance the negative effect of pairwise interactions on ethanol yield. In addition, we find that Lactobacillus amylovorus improves the yeast growth rate and ethanol yield by cross-feeding acetaldehyde, as shown by flux balance analysis and laboratory experiments. Our results suggest that Lactobacillus amylovorus could be considered a beneficial bacterium with the potential to improve sugarcane ethanol fermentation yields by almost 3%. These data highlight the biotechnological importance of comprehensively studying microbial communities and could be extended to other microbial systems with relevance to human health and the environment.

Suboptimal Global Transcriptional Response Increases the Harmful Effects of Loss-of-Function Mutations.

The fitness impact of loss-of-function mutations is generally assumed to reflect the loss of specific molecular functions associated with the perturbed gene. Here, we propose that rewiring of the transcriptome upon deleterious gene inactivation is frequently nonspecific and mimics stereotypic responses to external environmental change. Consequently, transcriptional response to gene deletion could be suboptimal and incur an extra fitness cost. Analysis of the transcriptomes of ∼1,500 single-gene deletion Saccharomyces cerevisiae strains supported this scenario. First, most transcriptomic changes are not specific to the deleted gene but are rather triggered by perturbations in functionally diverse genes. Second, gene deletions that alter the expression of dosage-sensitive genes are especially harmful. Third, by elevating the expression level of downregulated genes, we could experimentally mitigate the fitness defect of gene deletions. Our work shows that rewiring of genomic expression upon gene inactivation shapes the harmful effects of mutations.

Artificially selecting bacterial communities using propagule strategies.

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.

The ecology and evolution of microbial metabolic strategies.

Free-living microbes are generally capable of growing on multiple different nutrients. Some of those nutrients are used simultaneously, while others are used sequentially. The pattern of nutrient preferences and co-utilization defines the metabolic strategy of a microorganism. Metabolic strategies can substantially affect ecological interactions between species, but their evolution and distribution across the tree of life remain poorly characterized. We discuss how the confluence of better computational models of genotype-phenotype maps and high-throughput experimental tools can help us fill gaps in our knowledge and incorporate metabolic strategies into quantitative predictive models of microbial consortia.

High-order interactions distort the functional landscape of microbial consortia.

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.

On the deformability of an empirical fitness landscape by microbial evolution.

A fitness landscape is a map between the genotype and its reproductive success in a given environment. The topography of fitness landscapes largely governs adaptive dynamics, constraining evolutionary trajectories and the predictability of evolution. Theory suggests that this topography can be deformed by mutations that produce substantial changes to the environment. Despite its importance, the deformability of fitness landscapes has not been systematically studied beyond abstract models, and little is known about its reach and consequences in empirical systems. Here we have systematically characterized the deformability of the genome-wide metabolic fitness landscape of the bacterium Escherichia coli Deformability is quantified by the noncommutativity of epistatic interactions, which we experimentally demonstrate in mutant strains on the path to an evolutionary innovation. Our analysis shows that the deformation of fitness landscapes by metabolic mutations rarely affects evolutionary trajectories in the short range. However, mutations with large environmental effects produce long-range landscape deformations in distant regions of the genotype space that affect the fitness of later descendants. Our results therefore suggest that, even in situations in which mutations have strong environmental effects, fitness landscapes may retain their power to forecast evolution over small mutational distances despite the potential attenuation of that power over longer evolutionary trajectories. Our methods and results provide an avenue for integrating adaptive and eco-evolutionary dynamics with complex genetics and genomics.

Emergent simplicity in microbial community assembly.

A major unresolved question in microbiome research is whether the complex taxonomic architectures observed in surveys of natural communities can be explained and predicted by fundamental, quantitative principles. Bridging theory and experiment is hampered by the multiplicity of ecological processes that simultaneously affect community assembly in natural ecosystems. We addressed this challenge by monitoring the assembly of hundreds of soil- and plant-derived microbiomes in well-controlled minimal synthetic media. Both the community-level function and the coarse-grained taxonomy of the resulting communities are highly predictable and governed by nutrient availability, despite substantial species variability. By generalizing classical ecological models to include widespread nonspecific cross-feeding, we show that these features are all emergent properties of the assembly of large microbial communities, explaining their ubiquity in natural microbiomes.

Deconstructing a multiple antibiotic resistance regulation through the quantification of its input function.

Many essential bacterial responses present complex transcriptional regulation of gene expression. To what extent can the study of these responses substantiate the logic of their regulation? Here, we show how the input function of the genes constituting the response, i.e., the information of how their transcription rates change as function of the signals acting on the regulators, can serve as a quantitative tool to deconstruct the corresponding regulatory logic. To demonstrate this approach, we consider the multiple antibiotic resistance (mar) response in Escherichia coli. By characterizing the input function of its representative genes in wild-type and mutant bacteria, we recognize a dual autoregulation motif as main determinant of the response, which is further adjusted by the interplay with other regulators. We show that basic attributes, like its reaction to a wide range of stress or its moderate expression change, are associated with a strong negative autoregulation, while others, like the buffering of metabolic signals or the lack of memory to previous stress, are related to a weak positive autoregulation. With a mathematical model of the input functions, we identify some constraints fixing the molecular attributes of the regulators, and also notice the relevance of the bicystronic architecture harboring the dual autoregulation that is unique in E. coli. The input function emerges then as a tool to disentangle the rationale behind most of the attributes defining the mar phenotype. Overall, the present study supports the value of characterizing input functions to deconstruct the complexity of regulatory architectures in prokaryotic and eukaryotic systems.