Competitive resource allocation to metabolic pathways contributes to overflow metabolisms and emergent properties in cross-feeding microbial consortia

Cell Geometry and Membrane Protein Crowding Constrain Growth Rate, Overflow Metabolism, Respiration, and Maintenance Energy

A metabolic theory is presented for predicting maximum growth rate, overflow metabolism, respiration efficiency, and maintenance energy flux based on the intersection of cell geometry, membrane protein crowding, and metabolism. The importance of cytosolic macromolecular crowding on phenotype has been established in the literature but the importance of surface area has been largely overlooked due to incomplete knowledge of membrane properties. We demonstrate that the capacity of the membrane to host proteins increases with growth rate offsetting decreases in surface area-to-volume ratios (SA:V). This increase in membrane protein is hypothesized to be essential to competitive Escherichia coli phenotypes. The presented membrane-centric theory uses biophysical properties and metabolic systems analysis to successfully predict the phenotypes of E. coli K-12 strains, MG1655 and NCM3722, which are genetically similar but have SA:V ratios that differ up to 30%, maximum growth rates on glucose media that differ by 40%, and overflow phenotypes that start at growth rates that differ by 80%. These analyses did not consider cytosolic macromolecular crowding, highlighting the distinct properties of the presented theory. Cell geometry and membrane protein crowding are significant biophysical constraints on phenotype and provide a theoretical framework for improved understanding and control of cell biology.

Metabolic cross-feeding between the competent degrader Rhodococcus sp. strain p52 and an incompetent partner during catabolism of dibenzofuran: Understanding the leading and supporting roles

Microbial interactions, particularly metabolic cross-feeding, play important roles in removing recalcitrant environmental pollutants; however, the underlying mechanisms involved in this process remain unclear. Thus, this study aimed to elucidate the mechanism by which metabolic cross-feeding occurs during synergistic dibenzofuran degradation between a highly efficient degrader, Rhodococcus sp. strain p52, and a partner incapable of utilizing dibenzofuran. A bottom-up approach combined with pairwise coculturing was used to examine metabolic cross-feeding between strain p52 and Arthrobacter sp. W06 or Achromobacter sp. D10. Pairwise coculture not only promoted bacterial pair growth but also facilitated dibenzofuran degradation. Specifically, strain p52, acting as a donor, released dibenzofuran metabolic intermediates, including salicylic acid and gentisic acid, for utilization and growth, respectively, by the partner strains W06 and D10. Both salicylic acid and gentisic acid exhibited biotoxicity, and their accumulation inhibited dibenzofuran degradation. The transcriptional activity of the genes responsible for the catabolism of dibenzofuran and its metabolic intermediates was coordinately regulated in strain p52 and its cocultivated partners, thus achieving synergistic dibenzofuran degradation. This study provides insights into microbial metabolic cross-feeding during recalcitrant environmental pollutant removal.

Using metabolic networks to predict cross-feeding and competition interactions between microorganisms

Understanding the interactions between microorganisms and their impact on bacterial behavior at the community level is a key research topic in microbiology. Different methods, relying on experimental or mathematical approaches based on the diverse properties of bacteria, are currently employed to study these interactions. Recently, the use of metabolic networks to understand the interactions between bacterial pairs has increased, highlighting the relevance of this approach in characterizing bacteria. In this study, we leverage the representation of bacteria through their metabolic networks to build a predictive model aimed at reducing the number of experimental assays required for designing bacterial consortia with specific behaviors. Our novel method for predicting cross-feeding or competition interactions between pairs of microorganisms utilizes metabolic network features. Machine learning classifiers are employed to determine the type of interaction from automatically reconstructed metabolic networks. Several algorithms were assessed and selected based on comprehensive testing and careful separation of manually compiled data sets obtained from literature sources. We used different classification algorithms, including K Nearest Neighbors, XGBoost, Support Vector Machine, and Random Forest, tested different parameter values, and implemented several data curation approaches to reduce the biological bias associated with our data set, ultimately achieving an accuracy of over 0.9. Our method holds substantial potential to advance the understanding of community behavior and contribute to the development of more effective approaches for consortia design.IMPORTANCEUnderstanding bacterial interactions at the community level is critical for microbiology, and leveraging metabolic networks presents an efficient and effective approach. The introduction of this novel method for predicting interactions through machine learning classifiers has the potential to advance the field by reducing the number of experimental assays required and contributing to the development of more effective bacterial consortia.

Deciphering anaerobic ethanol metabolic pathways shaped by operational modes

Efficient anaerobic digestion requires the syntrophic cooperation among diverse microorganisms with various metabolic pathways. In this study, two operational modes, i.e., the sequencing batch reactor (SBR) and the continuous-flow reactor (CFR), were adopted in ethanol-fed systems with or without the supplement of powdered activated carbon (PAC) to examine their effects on ethanol metabolic pathways. Notably, the operational mode of SBR and the presence of CO2 facilitated ethanol metabolism towards propionate production. This was further evidenced by the dominance of Desulfobulbus, and the increased relative abundances of enzymes (EC: 1.2.7.1 and 1.2.7.11) involved in CO2 metabolism in SBRs. Moreover, SBRs exhibited superior biomass-based rates of ethanol degradation and methanogenesis, surpassing those in CFRs by 53.1% and 22.3%, respectively. Remarkably, CFRs with the extended solids retention time enriched high relative abundances of Geobacter of 71.7% and 70.4% under conditions with and without the addition of PAC, respectively. Although both long-term and short-term PAC additions led to the increased sludge conductivity and a reduced methanogenic lag phase, only the long-term PAC addition resulted in enhanced rates of ethanol degradation and propionate production/degradation. The strategies by adjusting operational mode and PAC addition could be adopted for modulating the anaerobic ethanol metabolic pathway and enriching Geobacter.

Synthesizing microbial biodiversity

The diversity of microbial ecosystems is linked to crucial ecological processes and functions. Despite its significance, the ecological mechanisms responsible for the initiation and maintenance of microbiome diversity are still not fully understood. The primary challenge lies in the difficulty of isolating, monitoring, and manipulating the complex and interrelated ecological processes that modulate the diversity of microbial communities in their natural habitats. Synthetic ecology experiments provide a suitable alternative for investigating the mechanisms behind microbial biodiversity in controlled laboratory settings, as the environment can be systematically and modularly manipulated by adding and removing components. This enables the testing of hypotheses and the advancement of predictive theories. In this review, we present an overview of recent progress toward achieving this goal.

Environment Constrains Fitness Advantages of Division of Labor in Microbial Consortia Engineered for Metabolite Push or Pull Interactions

Fitness benefits from division of labor are well documented in microbial consortia, but the dependency of the benefits on environmental context is poorly understood. Two synthetic Escherichia coli consortia were built to test the relationships between exchanged organic acid, local environment, and opportunity costs of different metabolic strategies. Opportunity costs quantify benefits not realized due to selecting one phenotype over another. The consortia catabolized glucose and exchanged either acetic or lactic acid to create producer-consumer food webs. The organic acids had different inhibitory properties and different opportunity costs associated with their positions in central metabolism. The exchanged metabolites modulated different consortial dynamics. The acetic acid-exchanging (AAE) consortium had a “push” interaction motif where acetic acid was secreted faster by the producer than the consumer imported it, while the lactic acid-exchanging (LAE) consortium had a “pull” interaction motif where the consumer imported lactic acid at a comparable rate to its production. The LAE consortium outperformed wild-type (WT) batch cultures under the environmental context of weakly buffered conditions, achieving a 55% increase in biomass titer, a 51% increase in biomass per proton yield, an 86% increase in substrate conversion, and the complete elimination of by-product accumulation all relative to the WT. However, the LAE consortium had the trade-off of a 42% lower specific growth rate. The AAE consortium did not outperform the WT in any considered performance metric. Performance advantages of the LAE consortium were sensitive to environment; increasing the medium buffering capacity negated the performance advantages compared to WT.

IMPORTANCE Most naturally occurring microorganisms persist in consortia where metabolic interactions are common and often essential to ecosystem function. This study uses synthetic ecology to test how different cellular interaction motifs influence performance properties of consortia. Environmental context ultimately controlled the division of labor performance as shifts from weakly buffered to highly buffered conditions negated the benefits of the strategy. Understanding the limits of division of labor advances our understanding of natural community functioning, which is central to nutrient cycling and provides design rules for assembling consortia used in applied bioprocessing.

Comparative Population Genetics in the Human Gut Microbiome

Genetic variation in the human gut microbiome is responsible for conferring a number of crucial phenotypes like the ability to digest food and metabolize drugs. Yet, our understanding of how this variation arises and is maintained remains relatively poor. Thus, the microbiome remains a largely untapped resource, as the large number of coexisting species in the microbiome presents a unique opportunity to compare and contrast evolutionary processes across species to identify universal trends and deviations. Here we outline features of the human gut microbiome that, while not unique in isolation, as an assemblage make it a system with unparalleled potential for comparative population genomics studies. We consciously take a broad view of comparative population genetics, emphasizing how sampling a large number of species allows researchers to identify universal evolutionary dynamics in addition to new genes, which can then be leveraged to identify exceptional species that deviate from general patterns. To highlight the potential power of comparative population genetics in the microbiome, we reanalyze patterns of purifying selection across ∼40 prevalent species in the human gut microbiome to identify intriguing trends which highlight functional categories in the microbiome that may be under more or less constraint.

Resource-allocation constraint governs structure and function of microbial communities in metabolic modeling

Predictive modeling tools for assessing microbial communities are important for realizing transformative capabilities of microbiomes in agriculture, ecology, and medicine. Constraint-based community-scale metabolic modeling is unique in its potential for making mechanistic predictions regarding both the structure and function of microbial communities. However, accessing this potential requires an understanding of key physicochemical constraints, which are typically considered on a per-species basis. What is needed is a means of incorporating global constraints relevant to microbial ecology into community models. Resource-allocation constraint, which describes how limited resources should be distributed to different cellular processes, sets limits on the efficiency of metabolic and ecological processes. In this study, we investigate the implications of resource-allocation constraints in community-scale metabolic modeling through a simple mechanism-agnostic implementation of resource-allocation constraints directly at the flux level. By systematically performing single-, two-, and multi-species growth simulations, we show that resource-allocation constraints are indispensable for predicting the structure and function of microbial communities. Our findings call for a scalable workflow for implementing a mechanistic version of resource-allocation constraints to ultimately harness the full potential of community-scale metabolic modeling tools.

Microbial Systems Ecology to Understand Cross-Feeding in Microbiomes

Understanding how microorganism-microorganism interactions shape microbial assemblages is a key to deciphering the evolution of dependencies and co-existence in complex microbiomes. Metabolic dependencies in cross-feeding exist in microbial communities and can at least partially determine microbial community composition. To parry the complexity and experimental limitations caused by the large number of possible interactions, new concepts from systems biology aim to decipher how the components of a system interact with each other. The idea that cross-feeding does impact microbiome assemblages has developed both theoretically and empirically, following a systems biology framework applied to microbial communities, formalized as microbial systems ecology (MSE) and relying on integrated-omics data. This framework merges cellular and community scales and offers new avenues to untangle microbial coexistence primarily by metabolic modeling, one of the main approaches used for mechanistic studies. In this mini-review, we first give a concise explanation of microbial cross-feeding. We then discuss how MSE can enable progress in microbial research. Finally, we provide an overview of a MSE framework mostly based on genome-scale metabolic-network reconstruction that combines top-down and bottom-up approaches to assess the molecular mechanisms of deterministic processes of microbial community assembly that is particularly suitable for use in synthetic biology and microbiome engineering.

Prospects for multi-omics in the microbial ecology of water engineering

Advances in high-throughput sequencing technologies and bioinformatics approaches over almost the last three decades have substantially increased our ability to explore microorganisms and their functions - including those that have yet to be cultivated in pure isolation. Genome-resolved metagenomic approaches have enabled linking powerful functional predictions to specific taxonomical groups with increasing fidelity. Additionally, related developments in both whole community gene expression surveys and metabolite profiling have permitted for direct surveys of community-scale functions in specific environmental settings. These advances have allowed for a shift in microbiome science away from descriptive studies and towards mechanistic and predictive frameworks for designing and harnessing microbial communities for desired beneficial outcomes. Water engineers, microbiologists, and microbial ecologists studying activated sludge, anaerobic digestion, and drinking water distribution systems have applied various (meta)omics techniques for connecting microbial community dynamics and physiologies to overall process parameters and system performance. However, the rapid pace at which new omics-based approaches are developed can appear daunting to those looking to apply these state-of-the-art practices for the first time. Here, we review how modern genome-resolved metagenomic approaches have been applied to a variety of water engineering applications from lab-scale bioreactors to full-scale systems. We describe integrated omics analysis across engineered water systems and the foundations for pairing these insights with modeling approaches. Lastly, we summarize emerging omics-based technologies that we believe will be powerful tools for water engineering applications. Overall, we provide a framework for microbial ecologists specializing in water engineering to apply cutting-edge omics approaches to their research questions to achieve novel functional insights. Successful adoption of predictive frameworks in engineered water systems could enable more economically and environmentally sustainable bioprocesses as demand for water and energy resources increases.

Stoichiometric Modeling of Artificial String Chemistries Reveals Constraints on Metabolic Network Structure

Uncovering the general principles that govern the structure of metabolic networks is key to understanding the emergence and evolution of living systems. Artificial chemistries can help illuminate this problem by enabling the exploration of chemical reaction universes that are constrained by general mathematical rules. Here, we focus on artificial chemistries in which strings of characters represent simplified molecules, and string concatenation and splitting represent possible chemical reactions. We developed a novel Python package, ARtificial CHemistry NEtwork Toolbox (ARCHNET), to study string chemistries using tools from the field of stoichiometric constraint-based modeling. In addition to exploring the topological characteristics of different string chemistry networks, we developed a network-pruning algorithm that can generate minimal metabolic networks capable of producing a specified set of biomass precursors from a given assortment of environmental nutrients. We found that the composition of these minimal metabolic networks was influenced more strongly by the metabolites in the biomass reaction than the identities of the environmental nutrients. This finding has important implications for the reconstruction of organismal metabolic networks and could help us better understand the rise and evolution of biochemical organization. More generally, our work provides a bridge between artificial chemistries and stoichiometric modeling, which can help address a broad range of open questions, from the spontaneous emergence of an organized metabolism to the structure of microbial communities.

Pseudomonas aeruginosa reverse diauxie is a multidimensional, optimized, resource utilization strategy

Pseudomonas aeruginosa is a globally-distributed bacterium often found in medical infections. The opportunistic pathogen uses a different, carbon catabolite repression (CCR) strategy than many, model microorganisms. It does not utilize a classic diauxie phenotype, nor does it follow common systems biology assumptions including preferential consumption of glucose with an 'overflow' metabolism. Despite these contradictions, P. aeruginosa is competitive in many, disparate environments underscoring knowledge gaps in microbial ecology and systems biology. Physiological, omics, and in silico analyses were used to quantify the P. aeruginosa CCR strategy known as 'reverse diauxie'. An ecological basis of reverse diauxie was identified using a genome-scale, metabolic model interrogated with in vitro omics data. Reverse diauxie preference for lower energy, nonfermentable carbon sources, such as acetate or succinate over glucose, was predicted using a multidimensional strategy which minimized resource investment into central metabolism while completely oxidizing substrates. Application of a common, in silico optimization criterion, which maximizes growth rate, did not predict the reverse diauxie phenotypes. This study quantifies P. aeruginosa metabolic strategies foundational to its wide distribution and virulence including its potentially, mutualistic interactions with microorganisms found commonly in the environment and in medical infections.

Enabling rational gut microbiome manipulations by understanding gut ecology through experimentally-evidenced in silico models

The gut microbiome has emerged as a contributing factor in non-communicable disease, rendering it a target of health-promoting interventions. Yet current understanding of the host-microbiome dynamic is insufficient to predict the variation in intervention outcomes across individuals. We explore the mechanisms that underpin the gut bacterial ecosystem and highlight how a more complete understanding of this ecology will enable improved intervention outcomes. This ecology varies within the gut over space and time. Interventions disrupt these processes, with cascading consequences throughout the ecosystem. In vivo studies cannot isolate and probe these processes at the required spatiotemporal resolutions, and in vitro studies lack the representative complexity required. However, we highlight that, together, both approaches can inform in silico models that integrate cellular-level dynamics, can extrapolate to explain bacterial community outcomes, permit experimentation and observation over ecological processes at high spatiotemporal resolution, and can serve as predictive platforms on which to prototype interventions. Thus, it is a concerted integration of these techniques that will enable rational targeted manipulations of the gut ecosystem.

Intraspecies Signaling between Common Variants of Pseudomonas aeruginosa Increases Production of Quorum-Sensing-Controlled Virulence Factors

The opportunistic pathogen Pseudomonas aeruginosa damages hosts through the production of diverse secreted products, many of which are regulated by quorum sensing (QS). The lasR gene, which encodes a central QS regulator, is frequently mutated in clinical isolates from chronic infections, and loss of LasR function (LasR-) generally impairs the activity of downstream QS regulators RhlR and PqsR. We found that in cocultures containing LasR+ and LasR- strains, LasR- strains hyperproduce the RhlR/RhlI-regulated antagonistic factors pyocyanin and rhamnolipids in diverse models and media and in different strain backgrounds. Diffusible QS autoinducers produced by the wild type were not required for this effect. Using transcriptomics, genetics, and biochemical approaches, we uncovered a reciprocal interaction between wild-type and lasR mutant pairs wherein the iron-scavenging siderophore pyochelin produced by the lasR mutant induced citrate release and cross-feeding from the wild type. Citrate, a metabolite often secreted in low iron environments, stimulated RhlR signaling and RhlI levels in LasR-but not in LasR+ strains. These studies reveal the potential for complex interactions between recently diverged, genetically distinct isolates within populations from single chronic infections.IMPORTANCE Coculture interactions between lasR loss-of-function and LasR+ Pseudomonas aeruginosa strains may explain the worse outcomes associated with the presence of LasR- strains. More broadly, this report illustrates how interactions within a genotypically diverse population, similar to those that frequently develop in natural settings, can promote unpredictably high virulence factor production.

Pseudomonad reverse carbon catabolite repression, interspecies metabolite exchange, and consortial division of labor

Microorganisms acquire energy and nutrients from dynamic environments, where substrates vary in both type and abundance. The regulatory system responsible for prioritizing preferred substrates is known as carbon catabolite repression (CCR). Two broad classes of CCR have been documented in the literature. The best described CCR strategy, referred to here as classic CCR (cCCR), has been experimentally and theoretically studied using model organisms such as Escherichia coli. cCCR phenotypes are often used to generalize universal strategies for fitness, sometimes incorrectly. For instance, extremely competitive microorganisms, such as Pseudomonads, which arguably have broader global distributions than E. coli, have achieved their success using metabolic strategies that are nearly opposite of cCCR. These organisms utilize a CCR strategy termed 'reverse CCR' (rCCR), because the order of preferred substrates is nearly reverse that of cCCR. rCCR phenotypes prefer organic acids over glucose, may or may not select preferred substrates to optimize growth rates, and do not allocate intracellular resources in a manner that produces an overflow metabolism. cCCR and rCCR have traditionally been interpreted from the perspective of monocultures, even though most microorganisms live in consortia. Here, we review the basic tenets of the two CCR strategies and consider these phenotypes from the perspective of resource acquisition in consortia, a scenario that surely influenced the evolution of cCCR and rCCR. For instance, cCCR and rCCR metabolism are near mirror images of each other; when considered from a consortium basis, the complementary properties of the two strategies can mitigate direct competition for energy and nutrients and instead establish cooperative division of labor.

Lytic bacteriophage have diverse indirect effects in a synthetic cross-feeding community

Bacteriophage shape the composition and function of microbial communities. Yet it remains difficult to predict the effect of phage on microbial interactions. Specifically, little is known about how phage influence mutualisms in networks of cross-feeding bacteria. We mathematically modeled the impacts of phage in a synthetic microbial community in which Escherichia coli and Salmonella enterica exchange essential metabolites. In this model, independent phage attack of either species was sufficient to temporarily inhibit both members of the mutualism; however, the evolution of phage resistance facilitated yields similar to those observed in the absence of phage. In laboratory experiments, attack of S. enterica with P22vir phage followed these modeling expectations of delayed community growth with little change in the final yield of bacteria. In contrast, when E. coli was attacked with T7 phage, S. enterica, the nonhost species, reached higher yields compared with no-phage controls. T7 infection increased nonhost yield by releasing consumable cell debris, and by driving evolution of partially resistant E. coli that secreted more carbon. Our results demonstrate that phage can have extensive indirect effects in microbial communities, that the nature of these indirect effects depends on metabolic and evolutionary mechanisms, and that knowing the degree of evolved resistance leads to qualitatively different predictions of bacterial community dynamics in response to phage attack.

Reverse diauxie phenotype in Pseudomonas aeruginosa biofilm revealed by exometabolomics and label-free proteomics

Microorganisms enhance fitness by prioritizing catabolism of available carbon sources using a process known as carbon catabolite repression (CCR). Planktonically grown Pseudomonas aeruginosa is known to prioritize the consumption of organic acids including lactic acid over catabolism of glucose using a CCR strategy termed "reverse diauxie." P. aeruginosa is an opportunistic pathogen with well-documented biofilm phenotypes that are distinct from its planktonic phenotypes. Reverse diauxie has been described in planktonic cultures, but it has not been documented explicitly in P. aeruginosa biofilms. Here a combination of exometabolomics and label-free proteomics was used to analyze planktonic and biofilm phenotypes for reverse diauxie. P. aeruginosa biofilm cultures preferentially consumed lactic acid over glucose, and in addition, the cultures catabolized the substrates completely and did not exhibit the acetate secreting "overflow" metabolism that is typical of many model microorganisms. The biofilm phenotype was enabled by changes in protein abundances, including lactate dehydrogenase, fumarate hydratase, GTP cyclohydrolase, L-ornithine N(5)-monooxygenase, and superoxide dismutase. These results are noteworthy because reverse diauxie-mediated catabolism of organic acids necessitates a terminal electron acceptor like O2, which is typically in low supply in biofilms due to diffusion limitation. Label-free proteomics identified dozens of proteins associated with biofilm formation including 16 that have not been previously reported, highlighting both the advantages of the methodology utilized here and the complexity of the proteomic adaptation for P. aeruginosa biofilms. Documenting the reverse diauxic phenotype in P. aeruginosa biofilms is foundational for understanding cellular nutrient and energy fluxes, which ultimately control growth and virulence.

Evolutionary causes and consequences of metabolic division of labour: why anaerobes do and aerobes don't

Metabolic division of the labour of organic matter decomposition into several steps carried out by different types of microbes is typical for many anoxic - but not oxic environments. An explanation of this well-known pattern is proposed based on the combination of three key insights: (i) well-studied anoxic environments are high flux environments: they are only anoxic because their high organic matter influx leads to oxygen depletion; (ii) shorter, incomplete catabolic pathways provide the capacity for higher flux, but this capacity is only advantageous in high flux environments; (iii) longer, complete catabolic pathways have energetic happy ends but only with high redox potential electron acceptors. Thus, aerobic environments favour longer pathways. Bioreactors, in contrast, are high flux environments and therefore favour division of catabolic labour even if aeration keeps them aerobic; therefore, host strains and feeding strategies must be carefully engineered to resist this pull.

In Silico Metabolic Design of Two-Strain Biofilm Systems Predicts Enhanced Biomass Production and Biochemical Synthesis

Engineered biofilm consortia have the potential to solve important biotechnological problems that have proved difficult for monoculture biofilms and planktonic consortia, such as conversion of lignocellulosic material to useful biochemicals. While considerable experimental progress has been reported for engineering and characterizing biofilm consortia, the field still lacks in silico tools for simulation, design, and optimization of stable, robust, and productive designed consortia. We developed biofilm consortia metabolic models for two coculture systems centered around the ecological design motif of a primary cell type that utilizes a supplied electron donor and secretes acetate as a byproduct and a secondary cell type that consumes the acetate, relieving byproduct inhibition on the primary cell type and enhancing overall system biomass. The models presented in this paper predict that distinct metabolic niches for the two cell types could be established by supplying electron donors and acceptors at opposite ends of the biofilm and that acetate consumption by the secondary cell type could increase total biomass accumulation and the synthesis of valuable biochemicals, such as isobutanol, by the primary cell type. System tunability is enhanced when each cell type is supplied with a unique terminal electron acceptor at opposite ends of the biofilm rather than competing for a common electron acceptor. Our model provides good qualitative agreement with data for a synthetic Escherichia coli coculture system, suggesting that the proposed design rules may have wide applicability to engineered biofilm consortia.

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.

Evolution of bidirectional costly mutualism from byproduct consumption

Mutualisms are essential for life, yet it is unclear how they arise. A two-stage process has been proposed for the evolution of mutualisms that involve exchanges of two costly resources. First, costly provisioning by one species may be selected for if that species gains a benefit from costless byproducts generated by a second species, and cooperators get disproportionate access to byproducts. Selection could then drive the second species to provide costly resources in return. Previously, a synthetic consortium evolved the first stage of this scenario: Salmonella enterica evolved costly production of methionine in exchange for costless carbon byproducts generated by an auxotrophic Escherichia coli Growth on agar plates localized the benefits of cooperation around methionine-secreting S. enterica Here, we report that further evolution of these partners on plates led to hypercooperative E. coli that secrete the sugar galactose. Sugar secretion arose repeatedly across replicate communities and is costly to E. coli producers, but enhances the growth of S. enterica The tradeoff between individual costs and group benefits led to maintenance of both cooperative and efficient E. coli genotypes in this spatially structured environment. This study provides an experimental example of de novo, bidirectional costly mutualism evolving from byproduct consumption. The results validate the plausibility of costly cooperation emerging from initially costless exchange, a scenario widely used to explain the origin of the mutualistic species interactions that are central to life on Earth.