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.The authors thank reviewers for comments that substantially improved this manuscript. BG and DS were partially supported by grants from the US Department of Energy (DE-SC0004962) and NIH (R01GM089978 and R01GM103502). (DE-SC0004962 - US Department of Energy; R01GM089978 - NIH; R01GM103502 - NIH)Published versio

Multiple long-term, experimentally-evolved populations of Escherichia coli acquire dependence upon citrate as an iron chelator for optimal growth on glucose

Background: Specialization for ecological niches is a balance of evolutionary adaptation and its accompanying tradeoffs. Here we focus on the Lenski Long-Term Evolution Experiment, which has maintained cultures of Escherichia coli in the same defined seasonal environment for 50,000 generations. Over this time, much adaptation and specialization to the environment has occurred. The presence of citrate in the growth media selected one lineage to gain the novel ability to utilize citrate as a carbon source after 31,000 generations. Here we test whether other strains have specialized to rely on citrate after 50,000 generations. Results: We show that in addition to the citrate-catabolizing strain, three other lineages evolving in parallel have acquired a dependence on citrate for optimal growth on glucose. None of these strains were stimulated indirectly by the sodium present in disodium citrate, nor exhibited even partial utilization of citrate as a carbon source. Instead, all three of these citrate-stimulated populations appear to rely on it as a chelator of iron. Conclusions: The strains we examine here have evolved specialization to their environment through apparent loss of function. Our results are most consistent with the accumulation of mutations in iron transport genes that were obviated by abundant citrate. The results present another example where a subtle decision in the design of an evolution experiment led to unexpected evolutionary outcomes.Organismic and Evolutionary Biolog

The Ability of Flux Balance Analysis to Predict Evolution of Central Metabolism Scales with the Initial Distance to the Optimum

The most powerful genome-scale framework to model metabolism, flux balance analysis (FBA), is an evolutionary optimality model. It hypothesizes selection upon a proposed optimality criterion in order to predict the set of internal fluxes that would maximize fitness. Here we present a direct test of the optimality assumption underlying FBA by comparing the central metabolic fluxes predicted by multiple criteria to changes measurable by a 13C-labeling method for experimentally-evolved strains. We considered datasets for three Escherichia coli evolution experiments that varied in their length, consistency of environment, and initial optimality. For ten populations that were evolved for 50,000 generations in glucose minimal medium, we observed modest changes in relative fluxes that led to small, but significant decreases in optimality and increased the distance to the predicted optimal flux distribution. In contrast, seven populations evolved on the poor substrate lactate for 900 generations collectively became more optimal and had flux distributions that moved toward predictions. For three pairs of central metabolic knockouts evolved on glucose for 600–800 generations, there was a balance between cases where optimality and flux patterns moved toward or away from FBA predictions. Despite this variation in predictability of changes in central metabolism, two generalities emerged. First, improved growth largely derived from evolved increases in the rate of substrate use. Second, FBA predictions bore out well for the two experiments initiated with ancestors with relatively sub-optimal yield, whereas those begun already quite optimal tended to move somewhat away from predictions. These findings suggest that the tradeoff between rate and yield is surprisingly modest. The observed positive correlation between rate and yield when adaptation initiated further from the optimum resulted in the ability of FBA to use stoichiometric constraints to predict the evolution of metabolism despite selection for rate

From Parasite to Mutualist: Rapid Evolution of Wolbachia in Natural Populations of Drosophila

Wolbachia are maternally inherited bacteria that commonly spread through host populations by causing cytoplasmic incompatibility, often expressed as reduced egg hatch when uninfected females mate with infected males. Infected females are frequently less fecund as a consequence of Wolbachia infection. However, theory predicts that because of maternal transmission, these “parasites” will tend to evolve towards a more mutualistic association with their hosts. Drosophila simulans in California provided the classic case of a Wolbachia infection spreading in nature. Cytoplasmic incompatibility allowed the infection to spread through individual populations within a few years and from southern to northern California (more than 700 km) within a decade, despite reducing the fecundity of infected females by 15%–20% under laboratory conditions. Here we show that the Wolbachia in California D. simulans have changed over the last 20 y so that infected females now exhibit an average 10% fecundity advantage over uninfected females in the laboratory. Our data suggest smaller but qualitatively similar changes in relative fecundity in nature and demonstrate that fecundity-increasing Wolbachia variants are currently polymorphic in natural populations

Identification of the potentiating mutations and synergistic epistasis that enabled the evolution of inter-species cooperation

Microbes often engage in cooperation through releasing biosynthetic compounds required by other species to grow. Given that production of costly biosynthetic metabolites is generally subjected to multiple layers of negative feedback, single mutations may frequently be insufficient to generate cooperative phenotypes. Synergistic epistatic interactions between multiple coordinated changes may thus often underlie the evolution of cooperation through overproduction of metabolites. To test the importance of synergistic mutations in cooperation we used an engineered bacterial consortium of an Escherichia coli methionine auxotroph and Salmonella enterica. S. enterica relies on carbon by-products from E. coli if lactose is the only carbon source. Directly selecting wild-type S. enterica in an environment that favored cooperation through secretion of methionine only once led to a methionine producer, and this producer both took a long time to emerge and was not very effective at cooperating. On the other hand, when an initial selection for resistance of S. enterica to a toxic methionine analog, ethionine, was used, subsequent selection for cooperation with E. coli was rapid, and the resulting double mutants were much more effective at cooperation. We found that potentiating mutations in metJ increase expression of metA, which encodes the first step of methionine biosynthesis. This increase in expression is required for the previously identified actualizing mutations in metA to generate cooperation. This work highlights that where biosynthesis of metabolites involves multiple layers of regulation, significant secretion of those metabolites may require multiple mutations, thereby constraining the evolution of cooperation

Parallel Mutations Result in a Wide Range of Cooperation and Community Consequences in a Two-Species Bacterial Consortium

Multi-species microbial communities play a critical role in human health, industry, and waste remediation. Recently, the evolution of synthetic consortia in the laboratory has enabled adaptation to be addressed in the context of interacting species. Using an engineered bacterial consortium, we repeatedly evolved cooperative genotypes and examined both the predictability of evolution and the phenotypes that determine community dynamics. Eight Salmonella enterica serovar Typhimurium strains evolved methionine excretion sufficient to support growth of an Escherichia coli methionine auxotroph, from whom they required excreted growth substrates. Non-synonymous mutations in metA, encoding homoserine trans-succinylase (HTS), were detected in each evolved S. enterica methionine cooperator and were shown to be necessary for cooperative consortia growth. Molecular modeling was used to predict that most of the non-synonymous mutations slightly increase the binding affinity for HTS homodimer formation. Despite this genetic parallelism and trend of increasing protein binding stability, these metA alleles gave rise to a wide range of phenotypic diversity in terms of individual versus group benefit. The cooperators with the highest methionine excretion permitted nearly two-fold faster consortia growth and supported the highest fraction of E. coli, yet also had the slowest individual growth rates compared to less cooperative strains. Thus, although the genetic basis of adaptation was quite similar across independent origins of cooperative phenotypes, quantitative measurements of metabolite production were required to predict either the individual-level growth consequences or how these propagate to community-level behavior

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

Population Dynamics Constrain the Cooperative Evolution of Cross-Feeding

Cross-feeding is the exchange of nutrients among species of microbes. It has two potential evolutionary origins, one as an exchange of metabolic wastes or byproducts among species, the other as a form of cooperation known as reciprocal altruism. This paper explores the conditions favoring the origin of cooperative cross-feeding between two species. There is an extensive literature on the evolution of cooperation, and some of the requirements for the evolution of cooperative cross-feeding follow from this prior work–specifically the requirement that interactions be limited to small groups of individuals, such as colonies in a spatially structured environment. Evolution of cooperative cross-feeding by a species also requires that cross-feeding from the partner species already exists, so that the cooperating mutant will automatically be reciprocated for its actions. Beyond these considerations, some unintuitive dynamical constraints apply. In particular, the benefit of cooperative cross-feeding applies only in the range of intermediate cell densities. At low density, resource concentrations are too low to offset the cost of cooperation. At high density, resources shared by both species become limiting, and the two species become competitors. These considerations suggest that the evolution of cooperative cross-feeding in nature may be more challenging than for other types of cooperation. However, the principles identified here may enable the experimental evolution of cross-feeding, as born out by a recent study