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.

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.

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.

•

We simulate how imbalances in NADH/NAD+ ratio modify cellular metabolic behavior

•

This reveals a mechanism to understand metabolic alterations at low growth rates

•

Tolerance to imbalance explains experimentally measured lifespan in yeast

•

We predict lifespan-protecting metabolites in yeast, animal, and human models

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.

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.

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.

Many cellular responses result from the integration of numerous regulatory signals. To deconstruct the regulation of one of these responses, which enables resistance to multiple antibiotics in Escherichia coli, a team led by Juan F. Poyatos at the National Center for Biotechnology in Madrid studied the response input function by combining theoretical models and experiments. This function quantifies the rate of transcription of the genes constituting the response with respect to the signals acting on its cognate regulators. By examining how the shape of the function changes in different situations, e.g., when a given regulator is mutated, the team identified the implications for the specificity and dynamics of the response of a dual autoregulation at the core of the control architecture. The use of input functions as quantitative tools allows us to reverse engineer the complex regulations that dictate essential physiological functions in both prokaryotic and eukaryotic cells.

Heart rate dynamics in doxorubicin-induced cardiomyopathy

The clinical use of doxorubicin, an effective chemotherapeutic is hampered by the development of irreversible cardiotoxicity. Here we test time-frequency analysis of heart rate (HR) variability (HRV) for early detection of doxorubicin-induced cardiotoxicity. Experiments were conducted in adult male Wistar rats treated for 15 days with doxorubicin (DOXO, total dose 15 mg kg(-1), i.p.) or saline (CONT). DOXO rats exhibited cardiotoxicity confirmed by histological examination without developing heart failure as estimated by echocardiography. However, HR variability increase reflected subtle microscopic changes of cardiac toxicity in DOXO rats. The results recommend time-frequency analysis of HRV for early detection of doxorubicin-induced cardiomyopathy.

Rewiring of genetic networks in response to modification of genetic background

Genome-scale genetic interaction networks are progressively contributing to map the molecular circuitry that determines cellular behavior. To what extent this mapping changes in response to different environmental or genetic conditions is, however, largely unknown. Here, we assembled a genetic network using an in silico model of metabolism in yeast to explicitly ask how separate genetic backgrounds alter network structure. Backgrounds defined by single deletions of metabolically active enzymes induce strong rewiring when the deletion corresponds to a catabolic gene, evidencing a broad redistribution of fluxes to alternative pathways. We also show how change is more pronounced in interactions linking genes in distinct functional modules and in those connections that present weak epistasis. These patterns reflect overall the distributed robustness of catabolism. In a second class of genetic backgrounds, in which a number of neutral mutations accumulate, we dominantly observe modifications in the negative interactions that together with an increase in the number of essential genes indicate a global reduction in buffering. Notably, neutral trajectories that originate considerable changes in the wild-type network comprise mutations that diminished the environmental plasticity of the corresponding metabolism, what emphasizes a mechanistic integration of genetic and environmental buffering. More generally, our work demonstrates how the specific mechanistic causes of robustness influence the architecture of multiconditional genetic interaction maps.

Balancing noise and plasticity in eukaryotic gene expression

BACKGROUND

Coupling the control of expression stochasticity (noise) to the ability of expression change (plasticity) can alter gene function and influence adaptation. A number of factors, such as transcription re-initiation, strong chromatin regulation or genome neighboring organization, underlie this coupling. However, these factors do not necessarily combine in equivalent ways and strengths in all genes. Can we identify then alternative architectures that modulate in distinct ways the linkage of noise and plasticity?

RESULTS

Here we first show that strong chromatin regulation, commonly viewed as a source of coupling, can lead to plasticity without noise. The nature of this regulation is relevant too, with plastic but noiseless genes being subjected to general activators whereas plastic and noisy genes experience more specific repression. Contrarily, in genes exhibiting poor transcriptional control, it is translational efficiency what separates noise from plasticity, a pattern related to transcript length. This additionally implies that genome neighboring organization -as modifier- appears only effective in highly plastic genes. In this class, we confirm bidirectional promoters (bipromoters) as a configuration capable to reduce coupling by abating noise but also reveal an important trade-off, since bipromoters also decrease plasticity. This presents ultimately a paradox between intergenic distances and modulation, with short intergenic distances both associated and disassociated to noise at different plasticity levels.

CONCLUSIONS

Balancing the coupling among different types of expression variability appears as a potential shaping force of genome regulation and organization. This is reflected in the use of different control strategies at genes with different sets of functional constraints.