Ecosystems biology of microbial metabolism

Review: Towards the next-generation models of the rumen microbiome for enhancing predictive power and guiding sustainable production strategies

The rumen ecosystem harbours a galaxy of microbes working in syntrophy to carry out a metabolic cascade of hydrolytic and fermentative reactions. This fermentation process allows ruminants to harvest nutrients from a wide range of feedstuff otherwise inaccessible to the host. The interconnection between the ruminant and its rumen microbiota shapes key animal phenotypes such as feed efficiency and methane emissions and suggests the potential of reducing methane emissions and enhancing feed conversion into animal products by manipulating the rumen microbiota. Whilst significant technological progress in omics techniques has increased our knowledge of the rumen microbiota and its genome (microbiome), translating omics knowledge into effective microbial manipulation strategies remains a great challenge. This challenge can be addressed by modelling approaches integrating causality principles and thus going beyond current correlation-based approaches applied to analyse rumen microbial genomic data. However, existing rumen models are not yet adapted to capitalise on microbial genomic information. This gap between the rumen microbiota available omics data and the way microbial metabolism is represented in the existing rumen models needs to be filled to enhance rumen understanding and produce better predictive models with capabilities for guiding nutritional strategies. To fill this gap, the integration of computational biology tools and mathematical modelling frameworks is needed to translate the information of the metabolic potential of the rumen microbes (inferred from their genomes) into a mathematical object. In this paper, we aim to discuss the potential use of two modelling approaches for the integration of microbial genomic information into dynamic models. The first modelling approach explores the theory of state observers to integrate microbial time series data into rumen fermentation models. The second approach is based on the genome-scale network reconstructions of rumen microbes. For a given microorganism, the network reconstruction produces a stoichiometry matrix of the metabolism. This matrix is the core of the so-called genome-scale metabolic models which can be exploited by a plethora of methods comprised within the constraint-based reconstruction and analysis approaches. We will discuss how these methods can be used to produce the next-generation models of the rumen microbiome.

Construction of Environmental Synthetic Microbial Consortia: Based on Engineering and Ecological Principles

In synthetic biology, engineering principles are applied to system design. The development of synthetic microbial consortia represents the intersection of synthetic biology and microbiology. Synthetic community systems are constructed by co-cultivating two or more microorganisms under certain environmental conditions, with broad applications in many fields including ecological restoration and ecological theory. Synthetic microbial consortia tend to have high biological processing efficiencies, because the division of labor reduces the metabolic burden of individual members. In this review, we focus on the environmental applications of synthetic microbial consortia. Although there are many strategies for the construction of synthetic microbial consortia, we mainly introduce the most widely used construction principles based on cross-feeding. Additionally, we propose methods for constructing synthetic microbial consortia based on traits and spatial structure from the perspective of ecology to provide a basis for future work.

Bottom-up synthetic ecology study of microbial consortia to enhance lignocellulose bioconversion

Lignocellulose is the most abundant organic carbon polymer on the earth. Its decomposition and conversion greatly impact the global carbon cycle. Furthermore, it provides feedstock for sustainable fuel and other value-added products. However, it continues to be underutilized, due to its highly recalcitrant and heterogeneric structure. Microorganisms, which have evolved versatile pathways to convert lignocellulose, undoubtedly are at the heart of lignocellulose conversion. Numerous studies that have reported successful metabolic engineering of individual strains to improve biological lignin valorization. Meanwhile, the bottleneck of single strain modification is becoming increasingly urgent in the conversion of complex substrates. Alternatively, increased attention has been paid to microbial consortia, as they show advantages over pure cultures, e.g., high efficiency and robustness. Here, we first review recent developments in microbial communities for lignocellulose bioconversion. Furthermore, the emerging area of synthetic ecology, which is an integration of synthetic biology, ecology, and computational biology, provides an opportunity for the bottom-up construction of microbial consortia. Then, we review different modes of microbial interaction and their molecular mechanisms, and discuss considerations of how to employ these interactions to construct synthetic consortia via synthetic ecology, as well as highlight emerging trends in engineering microbial communities for lignocellulose bioconversion.

Trophic Interactions and the Drivers of Microbial Community Assembly

Despite numerous surveys of gene and species content in heterotrophic microbial communities, such as those found in animal guts, oceans, or soils, it is still unclear whether there are generalizable biological or ecological processes that control their dynamics and function. Here, we review experimental and theoretical advances to argue that networks of trophic interactions, in which the metabolic excretions of one species are the primary resource for another, constitute the central drivers of microbial community assembly. Trophic interactions emerge from the deconstruction of complex forms of organic matter into a wealth of smaller metabolic intermediates, some of which are released to the environment and serve as a nutritional buffet for the community. The structure of the emergent trophic network and the rate at which primary resources are supplied control many features of microbial community assembly, including the relative contributions of competition and cooperation and the emergence of alternative community states. Viewing microbial community assembly through the lens of trophic interactions also has important implications for the spatial dynamics of communities as well as the functional redundancy of taxonomic groups. Given the ubiquity of trophic interactions across environments, they impart a common logic that can enable the development of a more quantitative and predictive microbial community ecology.

Recurrent Urinary Tract Infections: Unraveling the Complicated Environment of Uncomplicated rUTIs

Urinary tract infections (UTIs) are frequent in humans, affecting the upper and lower urinary tract. Present diagnosis relies on the positive culture of uropathogenic bacteria from urine and clinical markers of inflammation of the urinary tract. The bladder is constantly challenged by adverse environmental stimuli which influence urinary tract physiology, contributing to a dysbiotic environment. Simultaneously, pathogens are primed by environmental stressors such as antibiotics, favoring recurrent UTIs (rUTIs), resulting in chronic illness. Due to different confounders for UTI onset, a greater understanding of the fundamental environmental mechanisms and microbial ecology of the human urinary tract is required. Such advancements could promote the tandem translation of bench and computational studies for precision treatments and clinical management of UTIs. Therefore, there is an urgent need to understand the ecological interactions of the human urogenital microbial communities which precede rUTIs. This review aims to outline the mechanistic aspects of rUTI ecology underlying dysbiosis between both the human microbiome and host physiology which predisposes humans to rUTIs. By assessing the applications of next generation and systems level methods, we also recommend novel approaches to elucidate the systemic consequences of rUTIs which requires an integrated approach for successful treatment. To this end, we will provide an outlook towards the so-called 'uncomplicated environment of UTIs', a holistic and systems view that applies ecological principles to define patient-specific UTIs. This perspective illustrates the need to withdraw from traditional reductionist perspectives in infection biology and instead, a move towards a systems-view revolving around patient-specific pathophysiology during UTIs.

Non-additive microbial community responses to environmental complexity

Environmental composition is a major, though poorly understood, determinant of microbiome dynamics. Here we ask whether general principles govern how microbial community growth yield and diversity scale with an increasing number of environmental molecules. By assembling hundreds of synthetic consortia in vitro, we find that growth yield can remain constant or increase in a non-additive manner with environmental complexity. Conversely, taxonomic diversity is often much lower than expected. To better understand these deviations, we formulate metrics for epistatic interactions between environments and use them to compare our results to communities simulated with experimentally-parametrized consumer resource models. We find that key metabolic and ecological factors, including species similarity, degree of specialization, and metabolic interactions, modulate the observed non-additivity and govern the response of communities to combinations of resource pools. Our results demonstrate that environmental complexity alone is not sufficient for maintaining community diversity, and provide practical guidance for designing and controlling microbial ecosystems.

A multidimensional perspective on microbial interactions

Beyond being simply positive or negative, beneficial or inhibitory, microbial interactions can involve a diverse set of mechanisms, dependencies and dynamical properties. These more nuanced features have been described in great detail for some specific types of interactions, (e.g. pairwise metabolic cross-feeding, quorum sensing or antibiotic killing), often with the use of quantitative measurements and insight derived from modeling. With a growing understanding of the composition and dynamics of complex microbial communities for human health and other applications, we face the challenge of integrating information about these different interactions into comprehensive quantitative frameworks. Here, we review the literature on a wide set of microbial interactions, and explore the potential value of a formal categorization based on multidimensional vectors of attributes. We propose that such an encoding can facilitate systematic, direct comparisons of interaction mechanisms and dependencies, and we discuss the relevance of an atlas of interactions for future modeling and rational design efforts.

NetMet: A Network-Based Tool for Predicting Metabolic Capacities of Microbial Species and their Interactions

Metabolic conversions allow organisms to produce a set of essential metabolites from the available nutrients in an environment, frequently requiring metabolic exchanges among co-inhabiting organisms. Genomic-based metabolic simulations are being increasingly applied for exploring metabolic capacities, considering different environments and different combinations of microorganisms. NetMet is a web-based tool and a software package for predicting the metabolic performances of microorganisms and their corresponding combinations in user-defined environments. The algorithm takes, as input, lists of (i) species-specific enzymatic reactions (EC numbers), and (ii) relevant metabolic environments. The algorithm generates, as output, lists of (i) compounds that individual species can produce in each given environment, and (ii) compounds that are predicted to be produced through complementary interactions. The tool is demonstrated in two case studies. First, we compared the metabolic capacities of different haplotypes of the obligatory fruit and vegetable pathogen Candidatus Liberibacter solanacearum to those of their culturable taxonomic relative Liberibacter crescens. Second, we demonstrated the potential production of complementary metabolites by pairwise combinations of co-occurring endosymbionts of the plant phloem-feeding whitefly Bemisia tabaci.

Disruption of cross-feeding interactions by invading taxa can cause invasional meltdown in microbial communities

The strength of biotic interactions within an ecological community affects the susceptibility of the community to invasion by introduced taxa. In microbial communities, cross-feeding is a widespread type of biotic interaction that has the potential to affect community assembly and stability. Yet, there is little understanding of how the presence of cross-feeding within a community affects invasion risk. Here, I develop a metabolite-explicit model where native microbial taxa interact through both cross-feeding and competition for metabolites. I use this model to study how the strength of biotic interactions, especially cross-feeding, influence whether an introduced taxon can join the community. I found that stronger cross-feeding and competition led to much lower invasion risk, as both types of biotic interactions lead to greater metabolite scarcity for the invader. I also evaluated the impact of a successful invader on community composition and structure. The effect of invaders on the native community was greatest at intermediate levels of cross-feeding; at this ‘critical’ level of cross-feeding, successful invaders generally cause decreased diversity, decreased productivity, greater metabolite availability, and decreased quantities of metabolites exchanged among taxa. Furthermore, these changes resulting from a successful primary invader made communities further susceptible to future invaders. The increase in invasion risk was greatest when the network of metabolite exchange between taxa was minimally redundant. Thus, this model demonstrates a case of invasional meltdown that is mediated by initial invaders disrupting the metabolite exchange networks of the native community.

Microcins in Enterobacteriaceae: Peptide Antimicrobials in the Eco-Active Intestinal Chemosphere

Microcins are low-molecular-weight, ribosomally produced, highly stable, bacterial-inhibitory molecules involved in competitive, and amensalistic interactions between Enterobacteriaceae in the intestine. These interactions take place in a highly complex chemical landscape, the intestinal eco-active chemosphere, composed of chemical substances that positively or negatively influence bacterial growth, including those originated from nutrient uptake, and those produced by the action of the human or animal host and the intestinal microbiome. The contribution of bacteria results from their effect on the host generated molecules, on food and digested food, and organic substances from microbial origin, including from bacterial degradation. Here, we comprehensively review the main chemical substances present in the human intestinal chemosphere, particularly of those having inhibitory effects on microorganisms. With this background, and focusing on Enterobacteriaceae, the most relevant human pathogens from the intestinal microbiota, the microcin’s history and classification, mechanisms of action, and mechanisms involved in microcin’s immunity (in microcin producers) and resistance (non-producers) are reviewed. Products from the chemosphere likely modulate the ecological effects of microcin activity. Several cross-resistance mechanisms are shared by microcins, colicins, bacteriophages, and some conventional antibiotics, which are expected to produce cross-effects. Double-microcin-producing strains (such as microcins MccM and MccH47) have been successfully used for decades in the control of pathogenic gut organisms. Microcins are associated with successful gut colonization, facilitating translocation and invasion, leading to bacteremia, and urinary tract infections. In fact, Escherichia coli strains from the more invasive phylogroups (e.g., B2) are frequently microcinogenic. A publicly accessible APD3 database http://aps.unmc.edu/AP/ shows particular genes encoding microcins in 34.1% of E. coli strains (mostly MccV, MccM, MccH47, and MccI47), and much less in Shigella and Salmonella (<2%). Some 4.65% of Klebsiella pneumoniae are microcinogenic (mostly with MccE492), and even less in Enterobacter or Citrobacter (mostly MccS). The high frequency and variety of microcins in some Enterobacteriaceae indicate key ecological functions, a notion supported by their dominance in the intestinal microbiota of biosynthetic gene clusters involved in the synthesis of post-translationally modified peptide microcins.

Synthetic microbial consortia for biosynthesis and biodegradation: promises and challenges

Functional differentiation and metabolite exchange enable microbial consortia to perform complex metabolic tasks and efficiently cycle the nutrients. Inspired by the cooperative relationships in environmental microbial consortia, synthetic microbial consortia have great promise for studying the microbial interactions in nature and more importantly for various engineering applications. However, challenges coexist with promises, and the potential of consortium-based technologies is far from being fully harnessed. Thorough understanding of the underlying molecular mechanisms of microbial interactions is greatly needed for the rational design and optimization of defined consortia. These knowledge gaps could be potentially filled with the assistance of the ongoing revolution in systems biology and synthetic biology tools. As current fundamental and technical obstacles down the road being removed, we would expect new avenues with synthetic microbial consortia playing important roles in biological and environmental engineering processes such as bioproduction of desired chemicals and fuels, as well as biodegradation of persistent contaminants.

System-level analysis of metabolic trade-offs during anaerobic photoheterotrophic growth in Rhodopseudomonas palustris

BACKGROUND

Living organisms need to allocate their limited resources in a manner that optimizes their overall fitness by simultaneously achieving several different biological objectives. Examination of these biological trade-offs can provide invaluable information regarding the biophysical and biochemical bases behind observed cellular phenotypes. A quantitative knowledge of a cell system's critical objectives is also needed for engineering of cellular metabolism, where there is interest in mitigating the fitness costs that may result from human manipulation.

RESULTS

To study metabolism in photoheterotrophs, we developed and validated a genome-scale model of metabolism in Rhodopseudomonas palustris, a metabolically versatile gram-negative purple non-sulfur bacterium capable of growing phototrophically on various carbon sources, including inorganic carbon and aromatic compounds. To quantitatively assess trade-offs among a set of important biological objectives during different metabolic growth modes, we used our new model to conduct an 8-dimensional multi-objective flux analysis of metabolism in R. palustris. Our results revealed that phototrophic metabolism in R. palustris is light-limited under anaerobic conditions, regardless of the available carbon source. Under photoheterotrophic conditions, R. palustris prioritizes the optimization of carbon efficiency, followed by ATP production and biomass production rate, in a Pareto-optimal manner. To achieve maximum carbon fixation, cells appear to divert limited energy resources away from growth and toward CO2 fixation, even in the presence of excess reduced carbon. We also found that to achieve the theoretical maximum rate of biomass production, anaerobic metabolism requires import of additional compounds (such as protons) to serve as electron acceptors. Finally, we found that production of hydrogen gas, of potential interest as a candidate biofuel, lowers the cellular growth rates under all circumstances.

CONCLUSIONS

Photoheterotrophic metabolism of R. palustris is primarily regulated by the amount of light it can absorb and not the availability of carbon. However, despite carbon's secondary role as a regulating factor, R. palustris' metabolism strives for maximum carbon efficiency, even when this increased efficiency leads to slightly lower growth rates.

Designing Metabolic Division of Labor in Microbial Communities

Microbes face a trade-off between being metabolically independent and relying on neighboring organisms for the supply of some essential metabolites. This balance of conflicting strategies affects microbial community structure and dynamics, with important implications for microbiome research and synthetic ecology. A "gedanken" (thought) experiment to investigate this trade-off would involve monitoring the rise of mutual dependence as the number of metabolic reactions allowed in an organism is increasingly constrained. The expectation is that below a certain number of reactions, no individual organism would be able to grow in isolation and cross-feeding partnerships and division of labor would emerge. We implemented this idealized experiment using in silico genome-scale models. In particular, we used mixed-integer linear programming to identify trade-off solutions in communities of Escherichia coli strains. The strategies that we found revealed a large space of opportunities in nuanced and nonintuitive metabolic division of labor, including, for example, splitting the tricarboxylic acid (TCA) cycle into two separate halves. The systematic computation of possible solutions in division of labor for 1-, 2-, and 3-strain consortia resulted in a rich and complex landscape. This landscape displayed a nonlinear boundary, indicating that the loss of an intracellular reaction was not necessarily compensated for by a single imported metabolite. Different regions in this landscape were associated with specific solutions and patterns of exchanged metabolites. Our approach also predicts the existence of regions in this landscape where independent bacteria are viable but are outcompeted by cross-feeding pairs, providing a possible incentive for the rise of division of labor. IMPORTANCE Understanding how microbes assemble into communities is a fundamental open issue in biology, relevant to human health, metabolic engineering, and environmental sustainability. A possible mechanism for interactions of microbes is through cross-feeding, i.e., the exchange of small molecules. These metabolic exchanges may allow different microbes to specialize in distinct tasks and evolve division of labor. To systematically explore the space of possible strategies for division of labor, we applied advanced optimization algorithms to computational models of cellular metabolism. Specifically, we searched for communities able to survive under constraints (such as a limited number of reactions) that would not be sustainable by individual species. We found that predicted consortia partition metabolic pathways in ways that would be difficult to identify manually, possibly providing a competitive advantage over individual organisms. In addition to helping understand diversity in natural microbial communities, our approach could assist in the design of synthetic consortia.

Probabilistic Modeling of Microbial Metabolic Networks for Integrating Partial Quantitative Knowledge Within the Nitrogen Cycle

Understanding the interactions between microbial communities and their environment sufficiently to predict diversity on the basis of physicochemical parameters is a fundamental pursuit of microbial ecology that still eludes us. However, modeling microbial communities is problematic, because (i) communities are complex, (ii) most descriptions are qualitative, and (iii) quantitative understanding of the way communities interact with their surroundings remains incomplete. One approach to overcoming such complications is the integration of partial qualitative and quantitative descriptions into more complex networks. Here we outline the development of a probabilistic framework, based on Event Transition Graph (ETG) theory, to predict microbial community structure across observed chemical data. Using reverse engineering, we derive probabilities from the ETG that accurately represent observations from experiments and predict putative constraints on communities within dynamic environments. These predictions can feedback into the future development of field experiments by emphasizing the most important functional reactions, and associated microbial strains, required to characterize microbial ecosystems.

Integrated proteomic and metabolomic analysis of a reconstructed three-species microbial consortium for one-step fermentation of 2-keto-l-gulonic acid, the precursor of vitamin C

Microbial consortia, with the merits of strong stability, robustness, and multi-function, played critical roles in human health, bioenergy, and food manufacture, etc. On the basis of ‘build a consortium to understand it’, a novel microbial consortium consisted of Gluconobacter oxydans, Ketogulonicigenium vulgare and Bacillus endophyticus was reconstructed to produce 2-keto-l-gulonic acid (2-KGA), the precursor of vitamin C. With this synthetic consortium, 73.7 g/L 2-KGA was obtained within 30 h, which is comparable to the conventional industrial method. A combined time-series proteomic and metabolomic analysis of the fermentation process was conducted to further investigate the cell–cell interaction. The results suggested that the existence of B. endophyticus and G. oxydans together promoted the growth of K. vulgare by supplying additional nutrients, and promoted the 2-KGA production by supplying more substrate. Meanwhile, the growth of B. endophyticus and G. oxydans was compromised from the competition of the nutrients by K. vulgare, enabling the efficient production of 2-KGA. This study provides valuable guidance for further study of synthetic microbial consortia.

Universal scaling across biochemical networks on Earth

The application of network science to biology has advanced our understanding of the metabolism of individual organisms and the organization of ecosystems but has scarcely been applied to life at a planetary scale. To characterize planetary-scale biochemistry, we constructed biochemical networks using a global database of 28,146 annotated genomes and metagenomes and 8658 cataloged biochemical reactions. We uncover scaling laws governing biochemical diversity and network structure shared across levels of organization from individuals to ecosystems, to the biosphere as a whole. Comparing real biochemical reaction networks to random reaction networks reveals that the observed biological scaling is not a product of chemistry alone but instead emerges due to the particular structure of selected reactions commonly participating in living processes. We show that the topology of biochemical networks for the three domains of life is quantitatively distinguishable, with >80% accuracy in predicting evolutionary domain based on biochemical network size and average topology. Together, our results point to a deeper level of organization in biochemical networks than what has been understood so far.

Model Microbial Consortia as Tools for Understanding Complex Microbial Communities

A major biological challenge in the postgenomic era has been untangling the composition and functions of microbes that inhabit complex communities or microbiomes. Multi-omics and modern bioinformatics have provided the tools to assay molecules across different cellular and community scales; however, mechanistic knowledge over microbial interactions often remains elusive. This is due to the immense diversity and the essentially undiminished volume of not-yet-cultured microbes. Simplified model communities hold some promise in enabling researchers to manage complexity so that they can mechanistically understand the emergent properties of microbial community interactions. In this review, we surveyed several approaches that have effectively used tractable model consortia to elucidate the complex behavior of microbial communities. We go further to provide some perspectives on the limitations and new opportunities with these approaches and highlight where these efforts are likely to lead as advances are made in molecular ecology and systems biology.

Cell-to-cell bacterial interactions promoted by drier conditions on soil surfaces

Bacterial cell-to-cell interactions are in the core of evolutionary and ecological processes in soil and other environments. Under most conditions, natural soils are unsaturated where the fragmented aqueous habitats and thin liquid films confine bacterial cells within small volumes and close proximity for prolonged periods. We report effects of a range of hydration conditions on bacterial cell-level interactions that are marked by plasmid transfer between donor and recipient cells within populations of the soil bacterium Pseudomonas putida Using hydration-controlled sand microcosms, we demonstrate that the frequency of cell-to-cell contacts under prescribed hydration increases with lowering water potential values (i.e., under drier conditions where the aqueous phase shrinks and fragments). These observations were supported using a mechanistic individual-based model for linking macroscopic soil water potential to microscopic distribution of liquid phase and explicit bacterial cell interactions in a simplified porous medium. Model results are in good agreement with observations and inspire confidence in the underlying mechanisms. The study highlights important physical factors that control short-range bacterial cell interactions in soil and on surfaces, specifically, the central role of the aqueous phase in mediating bacterial interactions and conditions that promote genetic information transfer in support of soil microbial diversity.

An enormous potential for niche construction through bacterial cross-feeding in a homogeneous environment

Microorganisms modify their environment by excreting by-products of metabolism, which can create new ecological niches that can help microbial populations diversify. A striking example comes from experimental evolution of genetically identical Escherichia coli populations that are grown in a homogeneous environment with the single carbon source glucose. In such experiments, stable communities of genetically diverse cross-feeding E. coli cells readily emerge. Some cells that consume the primary carbon source glucose excrete a secondary carbon source, such as acetate, that sustains other community members. Few such cross-feeding polymorphisms are known experimentally, because they are difficult to screen for. We studied the potential of bacterial metabolism to create new ecological niches based on cross-feeding. To do so, we used genome scale models of the metabolism of E. coli and metabolisms of similar complexity, to identify unique pairs of primary and secondary carbon sources in these metabolisms. We then combined dynamic flux balance analysis with analytical calculations to identify which pair of carbon sources can sustain a polymorphic cross-feeding community. We identified almost 10,000 such pairs of carbon sources, each of them corresponding to a unique ecological niche. Bacterial metabolism shows an immense potential for the construction of new ecological niches through cross feeding.

Biodiversity can emerge in a completely homogeneous environment from populations with initially genetically identical individuals. This striking observation comes from experimental evolution of bacteria, which create new ecological niches when they excrete nutrient-rich waste products that can sustain the life of other bacteria. It is difficult to estimate the potential of any one organism for such metabolic niche construction experimentally, because it is challenging to screen for novel metabolic abilities on a large scale. We therefore used experimentally validated models of bacterial metabolism to predict how many novel niches organisms like Escherichia coli can construct, if a novel niche must be able to sustain a stable community of microbes that differ in the nutrients they consume. We identify thousands of such niches. They differ in their primary carbon source and a secondary carbon source that is excreted by some microbes and used by others. Because we restricted ourselves to chemically simple environments, we may even have underestimated the enormous potential of microbes for niche construction.

Host-microbe interaction in the gastrointestinal tract

The gastrointestinal tract is a highly complex organ in which multiple dynamic physiological processes are tightly coordinated while interacting with a dense and extremely diverse microbial population. From establishment in early life, through to host-microbe symbiosis in adulthood, the gut microbiota plays a vital role in our development and health. The effect of the microbiota on gut development and physiology is highlighted by anatomical and functional changes in germ-free mice, affecting the gut epithelium, immune system and enteric nervous system. Microbial colonisation promotes competent innate and acquired mucosal immune systems, epithelial renewal, barrier integrity, and mucosal vascularisation and innervation. Interacting or shared signalling pathways across different physiological systems of the gut could explain how all these changes are coordinated during postnatal colonisation, or after the introduction of microbiota into germ-free models. The application of cell-based in-vitro experimental systems and mathematical modelling can shed light on the molecular and signalling pathways which regulate the development and maintenance of homeostasis in the gut and beyond.

Modeling trophic dependencies and exchanges among insects' bacterial symbionts in a host-simulated environment

BACKGROUND

Individual organisms are linked to their communities and ecosystems via metabolic activities. Metabolic exchanges and co-dependencies have long been suggested to have a pivotal role in determining community structure. In phloem-feeding insects such metabolic interactions with bacteria enable complementation of their deprived nutrition. The phloem-feeding whitefly Bemisia tabaci (Hemiptera: Aleyrodidae) harbors an obligatory symbiotic bacterium, as well as varying combinations of facultative symbionts. This well-defined bacterial community in B. tabaci serves here as a case study for a comprehensive and systematic survey of metabolic interactions within the bacterial community and their associations with documented occurrences of bacterial combinations. We first reconstructed the metabolic networks of five common B. tabaci symbionts genera (Portiera, Rickettsia, Hamiltonella, Cardinium and Wolbachia), and then used network analysis approaches to predict: (1) species-specific metabolic capacities in a simulated bacteriocyte-like environment; (2) metabolic capacities of the corresponding species' combinations, and (3) dependencies of each species on different media components.

RESULTS

The predictions for metabolic capacities of the symbionts in the host environment were in general agreement with previously reported genome analyses, each focused on the single-species level. The analysis suggests several previously un-reported routes for complementary interactions and estimated the dependency of each symbiont in specific host metabolites. No clear association was detected between metabolic co-dependencies and co-occurrence patterns.

CONCLUSIONS

The analysis generated predictions for testable hypotheses of metabolic exchanges and co-dependencies in bacterial communities and by crossing them with co-occurrence profiles, contextualized interaction patterns into a wider ecological perspective.

From Genes to Ecosystems in Microbiology: Modeling Approaches and the Importance of Individuality

Models are important tools in microbial ecology. They can be used to advance understanding by helping to interpret observations and test hypotheses, and to predict the effects of ecosystem management actions or a different climate. Over the past decades, biological knowledge and ecosystem observations have advanced to the molecular and in particular gene level. However, microbial ecology models have changed less and a current challenge is to make them utilize the knowledge and observations at the genetic level. We review published models that explicitly consider genes and make predictions at the population or ecosystem level. The models can be grouped into three general approaches, i.e., metabolic flux, gene-centric and agent-based. We describe and contrast these approaches by applying them to a hypothetical ecosystem and discuss their strengths and weaknesses. An important distinguishing feature is how variation between individual cells (individuality) is handled. In microbial ecosystems, individual heterogeneity is generated by a number of mechanisms including stochastic interactions of molecules (e.g., gene expression), stochastic and deterministic cell division asymmetry, small-scale environmental heterogeneity, and differential transport in a heterogeneous environment. This heterogeneity can then be amplified and transferred to other cell properties by several mechanisms, including nutrient uptake, metabolism and growth, cell cycle asynchronicity and the effects of age and damage. For example, stochastic gene expression may lead to heterogeneity in nutrient uptake enzyme levels, which in turn results in heterogeneity in intracellular nutrient levels. Individuality can have important ecological consequences, including division of labor, bet hedging, aging and sub-optimality. Understanding the importance of individuality and the mechanism(s) underlying it for the specific microbial system and question investigated is essential for selecting the optimal modeling strategy.

Stable Isotope-Assisted Metabolic Profiling Reveals Growth Mode Dependent Differential Metabolism and Multiple Catabolic Pathways of l-Phenylalanine in Rubrivivax benzoatilyticus JA2

Anoxygenic phototrophic bacteria are metabolically versatile and survive under different growth modes using diverse organic compounds, yet their metabolic diversity is largely unexplored. In the present study, we employed stable-isotope-assisted metabolic profiling to unravel the l-phenylalanine catabolism in Rubrivivax benzoatilyticus JA2 under varying growth modes. Strain JA2 grows under anaerobic and aerobic conditions by utilizing l-phenylalanine as a nitrogen source. Furthermore, ring-labeled ¹³C₆-phenylalanine feeding followed by liquid chromatography-mass spectrometry exometabolite profiling revealed 60 labeled metabolic features (M + 6, M + 12, and M + 18) derived solely from l-phenylalanine, of which 11 were identified, 7 putatively identified, and 42 unidentified under anaerobic and aerobic conditions. However, labeled metabolites were significantly higher in aerobic compared to anaerobic conditions. Furthermore, detected metabolites and enzyme activities indicated multiple l-phenylalanine catabolic routes mainly Ehrlich, homogentisate-dependent melanin, benzenoid, and unidentified pathways operating under anaerobic and aerobic conditions in strain JA2. Interestingly, the study indicated l-phenylalanine-dependent and independent benzenoid biosynthesis in strain JA2 and a differential flux of l-phenylalanine to Ehrlich and benzenoid pathways under anaerobic and aerobic conditions. Additionally, unidentified labeled metabolites strongly suggest the presence of unknown phenylalanine catabolic routes in strain JA2. Overall, the study uncovered the l-phenylalanine catabolic diversity in strain JA2 and demonstrated the potential of stable isotope-assisted metabolomics in unraveling the hidden metabolic repertoire.

Analysis of Microbial Functions in the Rhizosphere Using a Metabolic-Network Based Framework for Metagenomics Interpretation

Advances in metagenomics enable high resolution description of complex bacterial communities in their natural environments. Consequently, conceptual approaches for community level functional analysis are in high need. Here, we introduce a framework for a metagenomics-based analysis of community functions. Environment-specific gene catalogs, derived from metagenomes, are processed into metabolic-network representation. By applying established ecological conventions, network-edges (metabolic functions) are assigned with taxonomic annotations according to the dominance level of specific groups. Once a function-taxonomy link is established, prediction of the impact of dominant taxa on the overall community performances is assessed by simulating removal or addition of edges (taxa associated functions). This approach is demonstrated on metagenomic data describing the microbial communities from the root environment of two crop plants – wheat and cucumber. Predictions for environment-dependent effects revealed differences between treatments (root vs. soil), corresponding to documented observations. Metabolism of specific plant exudates (e.g., organic acids, flavonoids) was linked with distinct taxonomic groups in simulated root, but not soil, environments. These dependencies point to the impact of these metabolite families as determinants of community structure. Simulations of the activity of pairwise combinations of taxonomic groups (order level) predicted the possible production of complementary metabolites. Complementation profiles allow formulating a possible metabolic role for observed co-occurrence patterns. For example, production of tryptophan-associated metabolites through complementary interactions is unique to the tryptophan-deficient cucumber root environment. Our approach enables formulation of testable predictions for species contribution to community activity and exploration of the functional outcome of structural shifts in complex bacterial communities. Understanding community-level metabolism is an essential step toward the manipulation and optimization of microbial function. Here, we introduce an analysis framework addressing three key challenges of such data: producing quantified links between taxonomy and function; contextualizing discrete functions into communal networks; and simulating environmental impact on community performances. New technologies will soon provide a high-coverage description of biotic and a-biotic aspects of complex microbial communities such as these found in gut and soil. This framework was designed to allow the integration of high-throughput metabolomic and metagenomic data toward tackling the intricate associations between community structure, community function, and metabolic inputs.

Compartmentalized metabolic network reconstruction of microbial communities to determine the effect of agricultural intervention on soils

Soil microbial communities are responsible for a wide range of ecological processes and have an important economic impact in agriculture. Determining the metabolic processes performed by microbial communities is crucial for understanding and managing ecosystem properties. Metagenomic approaches allow the elucidation of the main metabolic processes that determine the performance of microbial communities under different environmental conditions and perturbations. Here we present the first compartmentalized metabolic reconstruction at a metagenomics scale of a microbial ecosystem. This systematic approach conceives a meta-organism without boundaries between individual organisms and allows the in silico evaluation of the effect of agricultural intervention on soils at a metagenomics level. To characterize the microbial ecosystems, topological properties, taxonomic and metabolic profiles, as well as a Flux Balance Analysis (FBA) were considered. Furthermore, topological and optimization algorithms were implemented to carry out the curation of the models, to ensure the continuity of the fluxes between the metabolic pathways, and to confirm the metabolite exchange between subcellular compartments. The proposed models provide specific information about ecosystems that are generally overlooked in non-compartmentalized or non-curated networks, like the influence of transport reactions in the metabolic processes, especially the important effect on mitochondrial processes, as well as provide more accurate results of the fluxes used to optimize the metabolic processes within the microbial community.

Limitations of a metabolic network-based reverse ecology method for inferring host-pathogen interactions

Background

Host–pathogen interactions are important in a wide range of research fields. Given the importance of metabolic crosstalk between hosts and pathogens, a metabolic network-based reverse ecology method was proposed to infer these interactions. However, the validity of this method remains unclear because of the various explanations presented and the influence of potentially confounding factors that have thus far been neglected.

Results

We re-evaluated the importance of the reverse ecology method for evaluating host–pathogen interactions while statistically controlling for confounding effects using oxygen requirement, genome, metabolic network, and phylogeny data. Our data analyses showed that host–pathogen interactions were more strongly influenced by genome size, primary network parameters (e.g., number of edges), oxygen requirement, and phylogeny than the reserve ecology-based measures.

Conclusion

These results indicate the limitations of the reverse ecology method; however, they do not discount the importance of adopting reverse ecology approaches altogether. Rather, we highlight the need for developing more suitable methods for inferring host–pathogen interactions and conducting more careful examinations of the relationships between metabolic networks and host–pathogen interactions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1696-7) contains supplementary material, which is available to authorized users.

Community metabolic modeling approaches to understanding the gut microbiome: Bridging biochemistry and ecology

Interest in the human microbiome is at an all time high. The number of human microbiome studies is growing exponentially, as are reported associations between microbial communities and disease. However, we have not been able to translate the ever-growing amount of microbiome sequence data into better health. To do this, we need a practical means of transforming a disease-associated microbiome into a health-associated microbiome. This will require a framework that can be used to generate predictions about community dynamics within the microbiome under different conditions, predictions that can be tested and validated. In this review, using the gut microbiome to illustrate, we describe two classes of model that are currently being used to generate predictions about microbial community dynamics: ecological models and metabolic models. We outline the strengths and weaknesses of each approach and discuss the insights into the gut microbiome that have emerged from modeling thus far. We then argue that the two approaches can be combined to yield a community metabolic model, which will supply the framework needed to move from high-throughput omics data to testable predictions about how prebiotic, probiotic, and nutritional interventions affect the microbiome. We are confident that with a suitable model, researchers and clinicians will be able to harness the stream of sequence data and begin designing strategies to make targeted alterations to the microbiome and improve health.

Model-based quantification of metabolic interactions from dynamic microbial-community data

An important challenge in microbial ecology is to infer metabolic-exchange fluxes between growing microbial species from community-level data, concerning species abundances and metabolite concentrations. Here we apply a model-based approach to integrate such experimental data and thereby infer metabolic-exchange fluxes. We designed a synthetic anaerobic co-culture of Clostridium acetobutylicum and Wolinella succinogenes that interact via interspecies hydrogen transfer and applied different environmental conditions for which we expected the metabolic-exchange rates to change. We used stoichiometric models of the metabolism of the two microorganisms that represents our current physiological understanding and found that this understanding - the model - is sufficient to infer the identity and magnitude of the metabolic-exchange fluxes and it suggested unexpected interactions. Where the model could not fit all experimental data, it indicates specific requirement for further physiological studies. We show that the nitrogen source influences the rate of interspecies hydrogen transfer in the co-culture. Additionally, the model can predict the intracellular fluxes and optimal metabolic exchange rates, which can point to engineering strategies. This study therefore offers a realistic illustration of the strengths and weaknesses of model-based integration of heterogenous data that makes inference of metabolic-exchange fluxes possible from community-level experimental data.

A multi-objective constraint-based approach for modeling genome-scale microbial ecosystems

Interplay within microbial communities impacts ecosystems on several scales, and elucidation of the consequent effects is a difficult task in ecology. In particular, the integration of genome-scale data within quantitative models of microbial ecosystems remains elusive. This study advocates the use of constraint-based modeling to build predictive models from recent high-resolution -omics datasets. Following recent studies that have demonstrated the accuracy of constraint-based models (CBMs) for simulating single-strain metabolic networks, we sought to study microbial ecosystems as a combination of single-strain metabolic networks that exchange nutrients. This study presents two multi-objective extensions of CBMs for modeling communities: multi-objective flux balance analysis (MO-FBA) and multi-objective flux variability analysis (MO-FVA). Both methods were applied to a hot spring mat model ecosystem. As a result, multiple trade-offs between nutrients and growth rates, as well as thermodynamically favorable relative abundances at community level, were emphasized. We expect this approach to be used for integrating genomic information in microbial ecosystems. Following models will provide insights about behaviors (including diversity) that take place at the ecosystem scale.

Metabolic modelling in a dynamic evolutionary framework predicts adaptive diversification of bacteria in a long-term evolution experiment

BACKGROUND

Predicting adaptive trajectories is a major goal of evolutionary biology and useful for practical applications. Systems biology has enabled the development of genome-scale metabolic models. However, analysing these models via flux balance analysis (FBA) cannot predict many evolutionary outcomes including adaptive diversification, whereby an ancestral lineage diverges to fill multiple niches. Here we combine in silico evolution with FBA and apply this modelling framework, evoFBA, to a long-term evolution experiment with Escherichia coli.

RESULTS

Simulations predicted the adaptive diversification that occurred in one experimental population and generated hypotheses about the mechanisms that promoted coexistence of the diverged lineages. We experimentally tested and, on balance, verified these mechanisms, showing that diversification involved niche construction and character displacement through differential nutrient uptake and altered metabolic regulation.

CONCLUSION

The evoFBA framework represents a promising new way to model biochemical evolution, one that can generate testable predictions about evolutionary and ecosystem-level outcomes.

Metabolic Network Modeling of Microbial Interactions in Natural and Engineered Environmental Systems

We review approaches to characterize metabolic interactions within microbial communities using Stoichiometric Metabolic Network (SMN) models for applications in environmental and industrial biotechnology. SMN models are computational tools used to evaluate the metabolic engineering potential of various organisms. They have successfully been applied to design and optimize the microbial production of antibiotics, alcohols and amino acids by single strains. To date however, such models have been rarely applied to analyze and control the metabolism of more complex microbial communities. This is largely attributed to the diversity of microbial community functions, metabolisms, and interactions. Here, we firstly review different types of microbial interaction and describe their relevance for natural and engineered environmental processes. Next, we provide a general description of the essential methods of the SMN modeling workflow including the steps of network reconstruction, simulation through Flux Balance Analysis (FBA), experimental data gathering, and model calibration. Then we broadly describe and compare four approaches to model microbial interactions using metabolic networks, i.e., (i) lumped networks, (ii) compartment per guild networks, (iii) bi-level optimization simulations, and (iv) dynamic-SMN methods. These approaches can be used to integrate and analyze diverse microbial physiology, ecology and molecular community data. All of them (except the lumped approach) are suitable for incorporating species abundance data but so far they have been used only to model simple communities of two to eight different species. Interactions based on substrate exchange and competition can be directly modeled using the above approaches. However, interactions based on metabolic feedbacks, such as product inhibition and synthropy require extensions to current models, incorporating gene regulation and compounding accumulation mechanisms. SMN models of microbial interactions can be used to analyze complex “omics” data and to infer and optimize metabolic processes. Thereby, SMN models are suitable to capitalize on advances in high-throughput molecular and metabolic data generation. SMN models are starting to be applied to describe microbial interactions during wastewater treatment, in-situ bioremediation, microalgae blooms methanogenic fermentation, and bioplastic production. Despite their current challenges, we envisage that SMN models have future potential for the design and development of novel growth media, biochemical pathways and synthetic microbial associations.

Principles for designing synthetic microbial communities

Advances in synthetic biology to build microbes with defined and controllable properties are enabling new approaches to design and program multispecies communities. This emerging field of synthetic ecology will be important for many areas of biotechnology, bioenergy and bioremediation. This endeavor draws upon knowledge from synthetic biology, systems biology, microbial ecology and evolution. Fully realizing the potential of this discipline requires the development of new strategies to control the intercellular interactions, spatiotemporal coordination, robustness, stability and biocontainment of synthetic microbial communities. Here, we review recent experimental, analytical and computational advances to study and build multi-species microbial communities with defined functions and behavior for various applications. We also highlight outstanding challenges and future directions to advance this field.

Visualization of Metabolic Interaction Networks in Microbial Communities Using VisANT 5.0

The complexity of metabolic networks in microbial communities poses an unresolved visualization and interpretation challenge. We address this challenge in the newly expanded version of a software tool for the analysis of biological networks, VisANT 5.0. We focus in particular on facilitating the visual exploration of metabolic interaction between microbes in a community, e.g. as predicted by COMETS (Computation of Microbial Ecosystems in Time and Space), a dynamic stoichiometric modeling framework. Using VisANT's unique metagraph implementation, we show how one can use VisANT 5.0 to explore different time-dependent ecosystem-level metabolic networks. In particular, we analyze the metabolic interaction network between two bacteria previously shown to display an obligate cross-feeding interdependency. In addition, we illustrate how a putative minimal gut microbiome community could be represented in our framework, making it possible to highlight interactions across multiple coexisting species. We envisage that the "symbiotic layout" of VisANT can be employed as a general tool for the analysis of metabolism in complex microbial communities as well as heterogeneous human tissues. VisANT is freely available at: http://visant.bu.edu and COMETS at http://comets.bu.edu.

Increasing Metagenomic Resolution of Microbiome Interactions Through Functional Phylogenomics and Bacterial Sub-Communities

The genomic composition of the microbiome and its relationship with the environment is an exciting open question in biology. Metagenomics is a useful tool in the discovery of previously unknown taxa, but its use to understand the functional and ecological capacities of the microbiome is limited until taxonomy and function are understood in the context of the community. We suggest that this can be achieved using a combined functional phylogenomics and co-culture-based experimental strategy that can increase our capacity to measure sub-community interactions. Functional phylogenomics can identify and partition the genome such that hidden gene functions and gene clusters with unique evolutionary signals are revealed. We can test these phylogenomic predictions using an experimental model based on sub-community populations that represent a subset of the diversity directly obtained from environmental samples. These populations increase the detection of mechanisms that drive functional forces in the assembly of the microbiome, in particular the role of metabolites from key taxa in community interactions. Our combined approach leverages the potential of metagenomics to address biological questions from ecological systems.

Putative bacterial interactions from metagenomic knowledge with an integrative systems ecology approach

Following the trend of studies that investigate microbial ecosystems using different metagenomic techniques, we propose a new integrative systems ecology approach that aims to decipher functional roles within a consortium through the integration of genomic and metabolic knowledge at genome scale. For the sake of application, using public genomes of five bacterial strains involved in copper bioleaching: Acidiphilium cryptum, Acidithiobacillus ferrooxidans, Acidithiobacillus thiooxidans, Leptospirillum ferriphilum, and Sulfobacillus thermosulfidooxidans, we first reconstructed a global metabolic network. Next, using a parsimony assumption, we deciphered sets of genes, called Sets from Genome Segments (SGS), that (1) are close on their respective genomes, (2) take an active part in metabolic pathways and (3) whose associated metabolic reactions are also closely connected within metabolic networks. Overall, this SGS paradigm depicts genomic functional units that emphasize respective roles of bacterial strains to catalyze metabolic pathways and environmental processes. Our analysis suggested that only few functional metabolic genes are horizontally transferred within the consortium and that no single bacterial strain can accomplish by itself the whole copper bioleaching. The use of SGS pinpoints a functional compartmentalization among the investigated species and exhibits putative bacterial interactions necessary for promoting these pathways.

Calibration and analysis of genome-based models for microbial ecology

Microbial ecosystem modeling is complicated by the large number of unknown parameters and the lack of appropriate calibration tools. Here we present a novel computational framework for modeling microbial ecosystems, which combines genome-based model construction with statistical analysis and calibration to experimental data. Using this framework, we examined the dynamics of a community of Escherichia coli strains that emerged in laboratory evolution experiments, during which an ancestral strain diversified into two coexisting ecotypes. We constructed a microbial community model comprising the ancestral and the evolved strains, which we calibrated using separate monoculture experiments. Simulations reproduced the successional dynamics in the evolution experiments, and pathway activation patterns observed in microarray transcript profiles. Our approach yielded detailed insights into the metabolic processes that drove bacterial diversification, involving acetate cross-feeding and competition for organic carbon and oxygen. Our framework provides a missing link towards a data-driven mechanistic microbial ecology.

Metabolic network modeling of microbial communities

Genome-scale metabolic network reconstructions and constraint-based analyses are powerful methods that have the potential to make functional predictions about microbial communities. Genome-scale metabolic networks are used to characterize the metabolic functions of microbial communities via several techniques including species compartmentalization, separating species-level and community-level objectives, dynamic analysis, the 'enzyme-soup' approach, multiscale modeling, and others. There are many challenges in the field, including a need for tools that accurately assign high-level omics signals to individual community members, the need for improved automated network reconstruction methods, and novel algorithms for integrating omics data and engineering communities. As technologies and modeling frameworks improve, we expect that there will be corresponding advances in the fields of ecology, health science, and microbial community engineering.

The proportion of genes in a functional category is linked to mass-specific metabolic rate and lifespan

Metabolic rate and lifespan are important biological parameters that are studied in a wide range of research fields. They are known to correlate with body mass, but their association with gene (protein) functions is poorly understood. In this study, we collected data on the metabolic rate and lifespan of various organisms and investigated the relationship of these parameters with their genomes. We showed that the proportion of genes in a functional category, but not genome size, was correlated with mass-specific metabolic rate and maximal lifespan. In particular, the proportion of genes in oxic reactions (which occur in the presence of oxygen) was significantly associated with these two biological parameters. Additionally, we found that temperature, taxonomy, and mode-of-life traits had little effect on the observed associations. Our findings emphasize the importance of considering the biological functions of genes when investigating the relationships between genome, metabolic rate, and lifespan. Moreover, this provides further insights into these relationships, and may be useful for estimating metabolic rate and lifespan in individuals and the ecosystem using a combination of body mass measurements and genomic data.

Approaches and developments in studying the human microbiome network

The human microbiome is dynamic and unique to each individual, and its role is being increasingly recognized in healthy physiology and in disease, including gastrointestinal and neuropsychiatric disorders. Therefore, characterizing the human microbiome and the factors that shape its bacterial population, how they are related to host-specific attributes, and understanding the ways in which it can be manipulated and the phenotypic consequences of such manipulations are of great importance. Characterization of the microbiome so far has been mostly based on compositional studies alone, where relative abundances of different species are compared in different conditions, such as health and disease. However, inter-relationships among the bacterial species, such as competition and cooperation over metabolic resources, may be an important factor that affects the structure and function of the microbiome. Here we review the network-based approaches in answering such questions and explore the first attempts that focus on the interactions facet, complementing compositional studies, towards understanding the microbiome structure and its complex relationship with the human host.

Species interactions differ in their genetic robustness

Conflict and cooperation between bacterial species drive the composition and function of microbial communities. Stability of these emergent properties will be influenced by the degree to which species' interactions are robust to genetic perturbations. We use genome-scale metabolic modeling to computationally analyze the impact of genetic changes when Escherichia coli and Salmonella enterica compete, or cooperate. We systematically knocked out in silico each reaction in the metabolic network of E. coli to construct all 2583 mutant stoichiometric models. Then, using a recently developed multi-scale computational framework, we simulated the growth of each mutant E. coli in the presence of S. enterica. The type of interaction between species was set by modulating the initial metabolites present in the environment. We found that the community was most robust to genetic perturbations when the organisms were cooperating. Species ratios were more stable in the cooperative community, and community biomass had equal variance in the two contexts. Additionally, the number of mutations that have a substantial effect is lower when the species cooperate than when they are competing. In contrast, when mutations were added to the S. enterica network the system was more robust when the bacteria were competing. These results highlight the utility of connecting metabolic mechanisms and studies of ecological stability. Cooperation and conflict alter the connection between genetic changes and properties that emerge at higher levels of biological organization.

Evolutionary limits to cooperation in microbial communities

Microbes produce many compounds that are costly to a focal cell but promote the survival and reproduction of neighboring cells. This observation has led to the suggestion that microbial strains and species will commonly cooperate by exchanging compounds. Here, we examine this idea with an ecoevolutionary model where microbes make multiple secretions, which can be exchanged among genotypes. We show that cooperation between genotypes only evolves under specific demographic regimes characterized by intermediate genetic mixing. The key constraint on cooperative exchanges is a loss of autonomy: strains become reliant on complementary genotypes that may not be reliably encountered. Moreover, the form of cooperation that we observe arises through mutual exploitation that is related to cheating and "Black Queen" evolution for a single secretion. A major corollary is that the evolution of cooperative exchanges reduces community productivity relative to an autonomous strain that makes everything it needs. This prediction finds support in recent work from synthetic communities. Overall, our work suggests that natural selection will often limit cooperative exchanges in microbial communities and that, when exchanges do occur, they can be an inefficient solution to group living.

Dynamic flux balance analysis for synthetic microbial communities

Dynamic flux balance analysis (DFBA) is an extension of classical flux balance analysis that allows the dynamic effects of the extracellular environment on microbial metabolism to be predicted and optimised. Recently this computational framework has been extended to microbial communities for which the individual species are known and genome-scale metabolic reconstructions are available. In this review, the authors provide an overview of the emerging DFBA approach with a focus on two case studies involving the conversion of mixed hexose/pentose sugar mixtures by synthetic microbial co-culture systems. These case studies illustrate the key requirements of the DFBA approach, including the incorporation of individual species metabolic reconstructions, formulation of extracellular mass balances, identification of substrate uptake kinetics, numerical solution of the coupled linear program/differential equations and model adaptation for common, suboptimal growth conditions and identified species interactions. The review concludes with a summary of progress to date and possible directions for future research.

Metabolic resource allocation in individual microbes determines ecosystem interactions and spatial dynamics

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.

Syntrophic exchange in synthetic microbial communities

Metabolic crossfeeding is an important process that can broadly shape microbial communities. However, little is known about specific crossfeeding principles that drive the formation and maintenance of individuals within a mixed population. Here, we devised a series of synthetic syntrophic communities to probe the complex interactions underlying metabolic exchange of amino acids. We experimentally analyzed multimember, multidimensional communities of Escherichia coli of increasing sophistication to assess the outcomes of synergistic crossfeeding. We find that biosynthetically costly amino acids including methionine, lysine, isoleucine, arginine, and aromatics, tend to promote stronger cooperative interactions than amino acids that are cheaper to produce. Furthermore, cells that share common intermediates along branching pathways yielded more synergistic growth, but exhibited many instances of both positive and negative epistasis when these interactions scaled to higher dimensions. In more complex communities, we find certain members exhibiting keystone species-like behavior that drastically impact the community dynamics. Based on comparative genomic analysis of >6,000 sequenced bacteria from diverse environments, we present evidence suggesting that amino acid biosynthesis has been broadly optimized to reduce individual metabolic burden in favor of enhanced crossfeeding to support synergistic growth across the biosphere. These results improve our basic understanding of microbial syntrophy while also highlighting the utility and limitations of current modeling approaches to describe the dynamic complexities underlying microbial ecosystems. This work sets the foundation for future endeavors to resolve key questions in microbial ecology and evolution, and presents a platform to develop better and more robust engineered synthetic communities for industrial biotechnology.

Synthetic microbial ecosystems for biotechnology

Most highly controlled and specific applications of microorganisms in biotechnology involve pure cultures. Maintaining single strain cultures is important for industry as contaminants can reduce productivity and lead to longer "down-times" during sterilisation. However, microbes working together provide distinct advantages over pure cultures. They can undertake more metabolically complex tasks, improve efficiency and even expand applications to open systems. By combining rapidly advancing technologies with ecological theory, the use of microbial ecosystems in biotechnology will inevitably increase. This review provides insight into the use of synthetic microbial communities in biotechnology by applying the engineering paradigm of measure, model, manipulate and manufacture, and illustrate the emerging wider potential of the synthetic ecology field. Systems to improve biofuel production using microalgae are also discussed.

Age-correlated stress resistance improves fitness of yeast: support from agent-based simulations

BACKGROUND

Resistance to stress is often heterogeneous among individuals within a population, which helps protect against intermittent stress (bet hedging). This is also the case for heat shock resistance in the budding yeast Saccharomyces cerevisiae. Interestingly, the resistance appears to be continuously distributed (vs. binary, switch-like) and correlated with replicative age (vs. random). Older, slower-growing cells are more resistant than younger, faster-growing ones. Is there a fitness benefit to age-correlated stress resistance?

RESULTS

Here this hypothesis is explored using a simple agent-based model, which simulates a population of individual cells that grow and replicate. Cells age by accumulating damage, which lowers their growth rate. They synthesize trehalose at a metabolic cost, which helps protect against heat shock. Proteins Tsl1 and Tps3 (trehalose synthase complex regulatory subunit TSL1 and TPS3) represent the trehalose synthesis complex and they are expressed using constant, age-dependent and stochastic terms. The model was constrained by calibration and comparison to data from the literature, including individual-based observations obtained using high-throughput microscopy and flow cytometry. A heterogeneity network was developed, which highlights the predominant sources and pathways of resistance heterogeneity. To determine the best trehalose synthesis strategy, model strains with different Tsl1/Tps3 expression parameters were placed in competition in an environment with intermittent heat shocks.

CONCLUSIONS

For high severities and low frequencies of heat shock, the winning strain used an age-dependent bet hedging strategy, which shows that there can be a benefit to age-correlated stress resistance. The study also illustrates the utility of combining individual-based observations and modeling to understand mechanisms underlying population heterogeneity, and the effect on fitness.

Computational approaches for understanding energy metabolism

There has been a surge of interest in understanding the regulation of metabolic networks involved in disease in recent years. Quantitative models are increasingly being used to interrogate the metabolic pathways that are contained within this complex disease biology. At the core of this effort is the mathematical modeling of central carbon metabolism involving glycolysis and the citric acid cycle (referred to as energy metabolism). Here, we discuss several approaches used to quantitatively model metabolic pathways relating to energy metabolism and discuss their formalisms, successes, and limitations.