Species interactions constrain adaptation and preserve ecological stability in an experimental microbial community

Species loss within a microbial community can increase resource availability and spur adaptive evolution. Environmental shifts that cause species loss or fluctuations in community composition are expected to become more common, so it is important to understand the evolutionary forces that shape the stability and function of the emergent community. Here we study experimental cultures of a simple, ecologically stable community of Saccharomyces cerevisiae and Lactobacillus plantarum, in order to understand how the presence or absence of a species impacts coexistence over evolutionary timescales. We found that evolution in coculture led to drastically altered evolutionary outcomes for L. plantarum, but not S. cerevisiae. Both monoculture- and co-culture-evolved L. plantarum evolved dozens of mutations over 925 generations of evolution, but only L. plantarum that had evolved in isolation from S. cerevisiae lost the capacity to coexist with S. cerevisiae. We find that the evolutionary loss of ecological stability corresponds with fitness differences between monoculture-evolved L. plantarum and S. cerevisiae and genetic changes that repeatedly evolve across the replicate populations of L. plantarum. This work shows how coevolution within a community can prevent destabilising evolution in individual species, thereby preserving ecological diversity and stability, despite rapid adaptation.Subject terms: Evolution, Microbial ecology     

Evolutionary coexistence in a fluctuating environment by specialization on resource level

Microbial communities in fluctuating environments, such as oceans or the human gut, contain a wealth of diversity. This diversity contributes to the stability of communities and the functions they have in their hosts and ecosystems. To improve stability and increase production of beneficial compounds, we need to understand the underlying mechanisms causing this diversity. When nutrient levels fluctuate over time, one possibly relevant mechanism is coexistence between specialists on low and specialists on high nutrient levels. The relevance of this process is supported by the observations of coexistence in the laboratory, and by simple models, which show that negative frequency dependence of two such specialists can stabilize coexistence. However, as microbial populations are often large and fast growing, they evolve rapidly. Our aim is to determine what happens when species can evolve; whether evolutionary branching can create diversity or whether evolution will destabilize coexistence. We derive an analytical expression of the invasion fitness in fluctuating environments and use adaptive dynamics techniques to find that evolutionarily stable coexistence requires a special type of trade-off between growth at low and high nutrients. We do not find support for the necessary evolutionary trade-off in data available for the bacterium Escherichia coli and the yeast Saccharomyces cerevisiae on glucose. However, this type of data is scarce and might exist for other species or in different conditions. Moreover, we do find evidence for evolutionarily stable coexistence of the two species together. Since we find this coexistence in the scarce data that are available, we predict that specialization on resource level is a relevant mechanism for species diversity in microbial communities in fluctuating environments in natural settings.     

Engineering oxygen-independent biotin biosynthesis in Saccharomyces cerevisiae

An oxygen requirement for de novo biotin synthesis in Saccharomyces cerevisiae precludes the application of biotin-prototrophic strains in anoxic processes that use biotin-free media. To overcome this issue, this study explores introduction of the oxygen-independent Escherichia coli biotin-biosynthesis pathway in S. cerevisiae. Implementation of this pathway required expression of seven E. coli genes involved in fatty-acid synthesis and three E. coli genes essential for the formation of a pimelate thioester, key precursor of biotin synthesis. A yeast strain expressing these genes readily grew in biotin-free medium, irrespective of the presence of oxygen. However, the engineered strain exhibited specific growth rates 25% lower in biotin-free media than in biotin-supplemented media. Following adaptive laboratory evolution in anoxic cultures, evolved cell lines that no longer showed this growth difference in controlled bioreactors, were characterized by genome sequencing and proteome analyses. The evolved isolates exhibited a whole-genome duplication accompanied with an alteration in the relative gene dosages of biosynthetic pathway genes. These alterations resulted in a reduced abundance of the enzymes catalyzing the first three steps of the E. coli biotin pathway. The evolved pathway configuration was reverse engineered in the diploid industrial S. cerevisiae strain Ethanol Red. The resulting strain grew at nearly the same rate in biotin-supplemented and biotin-free media non-controlled batches performed in an anaerobic chamber. This study established an unique genetic engineering strategy to enable biotin-independent anoxic growth of S. cerevisiae and demonstrated its portability in industrial strain backgrounds.     

Rapid evolution destabilizes species interactions in a fluctuating environment

Positive species interactions underlie the functioning of ecosystems. Given their importance, it is crucial to understand the stability of positive interactions over evolutionary timescales, in both constant and fluctuating environments; e.g., environments interrupted with periods of competition. We addressed this question using a two-species microbial system in which we modulated interactions according to the nutrient provided. We evolved in parallel four experimental replicates of species growing in isolation or together in consortia for 200 generations in both a constant and fluctuating environment with daily changes between commensalism and competition. We sequenced full genomes of single clones isolated at different time points during the experiment. We found that the two species coexisted over 200 generations in the constant commensal environment. In contrast, in the fluctuating environment, coexistence broke down when one of the species went extinct in two out of four cases. We showed that extinction was highly deterministic: when we replayed the evolution experiment from an intermediate time point we repeatably reproduced species extinction. We further show that these dynamics were driven by adaptive mutations in a small set of genes. In conclusion, in a fluctuating environment, rapid evolution destabilizes the long-term stability of positive pairwise interactions.Subject terms: Microbial ecology, Bacterial genetics, Population genetics     

Synergistic epistasis enhances the co-operativity of mutualistic interspecies interactions

Early evolution of mutualism is characterized by big and predictable adaptive changes, including the specialization of interacting partners, such as through deleterious mutations in genes not required for metabolic cross-feeding. We sought to investigate whether these early mutations improve cooperativity by manifesting in synergistic epistasis between genomes of the mutually interacting species. Specifically, we have characterized evolutionary trajectories of syntrophic interactions of Desulfovibrio vulgaris (Dv) with Methanococcus maripaludis (Mm) by longitudinally monitoring mutations accumulated over 1000 generations of nine independently evolved communities with analysis of the genotypic structure of one community down to the single-cell level. We discovered extensive parallelism across communities despite considerable variance in their evolutionary trajectories and the perseverance within many evolution lines of a rare lineage of Dv that retained sulfate-respiration (SR+) capability, which is not required for metabolic cross-feeding. An in-depth investigation revealed that synergistic epistasis across pairings of Dv and Mm genotypes had enhanced cooperativity within SR− and SR+ assemblages, enabling their coexistence within the same community. Thus, our findings demonstrate that cooperativity of a mutualism can improve through synergistic epistasis between genomes of the interacting species, enabling the coexistence of mutualistic assemblages of generalists and their specialized variants.Subject terms: Microbial ecology, Population genetics, Symbiosis, Population dynamics, Molecular evolution     

Protection of Salmonella by ampicillin-resistant Escherichia coli in the presence of otherwise lethal drug concentrations

Microbial systems have become the preferred testing grounds for experimental work on the evolution of traits that benefit other group members. This work, based on conceptual and theoretical models of frequency-dependent selection within populations, has proven fruitful in terms of understanding the dynamics of group beneficial or ‘public goods’ traits within species. Here, we expand the scope of microbial work on the evolution of group-beneficial traits to the case of multi-species communities, particularly those that affect human health. We examined whether β-lactamase-producing Escherichia coli could protect ampicillin-sensitive cohorts of other species, particularly species that could cause human disease. Both β-lactamase-secreting E. coli and, surprisingly, those engineered to retain it, allowed for survival of a large number of ampicillin-sensitive cohorts of Salmonella enterica serovar Typhimurium, including both laboratory and clinical isolates. The Salmonella survivors, however, remained sensitive to ampicillin when re-plated onto solid medium and there was no evidence of gene transfer. Salmonella survival did not even require direct physical contact with the resistant E. coli. The observed phenomenon appears to involve increased release of β-lactamase from the E. coli when present with S. enterica. Significantly, these findings imply that resistant E. coli, that are not themselves pathogenic, may be exploited, even when they are normally selfish with respect to other E. coli. Thus, Salmonella can gain protection against antibiotics from E. coli without gene transfer, a phenomenon not previously known. As a consequence, antibiotic-resistant E. coli can play a decisive role in the survival of a species that causes disease and may thereby interfere with successful treatment.     

Genome-Scale Co-Expression Network Comparison across Escherichia coli and Salmonella enterica Serovar Typhimurium Reveals Significant Conservation at the Regulon Level of Local Regulators Despite Their Dissimilar Lifestyles

Availability of genome-wide gene expression datasets provides the opportunity to study gene expression across different organisms under a plethora of experimental conditions. In our previous work, we developed an algorithm called COMODO (COnserved MODules across Organisms) that identifies conserved expression modules between two species. In the present study, we expanded COMODO to detect the co-expression conservation across three organisms by adapting the statistics behind it. We applied COMODO to study expression conservation/divergence between Escherichia coli, Salmonella enterica, and Bacillus subtilis. We observed that some parts of the regulatory interaction networks were conserved between E. coli and S. enterica especially in the regulon of local regulators. However, such conservation was not observed between the regulatory interaction networks of B. subtilis and the two other species. We found co-expression conservation on a number of genes involved in quorum sensing, but almost no conservation for genes involved in pathogenicity across E. coli and S. enterica which could partially explain their different lifestyles. We concluded that despite their different lifestyles, no significant rewiring have occurred at the level of local regulons involved for instance, and notable conservation can be detected in signaling pathways and stress sensing in the phylogenetically close species S. enterica and E. coli. Moreover, conservation of local regulons seems to depend on the evolutionary time of divergence across species disappearing at larger distances as shown by the comparison with B. subtilis. Global regulons follow a different trend and show major rewiring even at the limited evolutionary distance that separates E. coli and S. enterica.     

Differential evolution of the Saccharomyces cerevisiae DUP240 paralogs and implication of recombination in phylogeny

Multigene families are observed in all genomes sequenced so far and are the reflection of key evolutionary mechanisms. The DUP240 family, identified in Saccharomyces cerevisiae strain S288C, is composed of 10 paralogs: seven are organized as two tandem repeats and three are solo ORFs. To investigate the evolution of the three solo paralogs, YAR023c, YCR007c and YHL044w, we performed a comparative analysis between 15 S.cerevisiae strains. These three ORFs are present in all strains and the conservation of synteny indicates that they are not frequently involved in chromosomal reshaping, in contrast to the DUP240 ORFs organized in tandem repeats. Our analysis of nucleotide and amino acid variations indicates that YAR023c and YHL044w fix mutations more easily than YCR007c, although they all belong to the same multigene family. This comparative analysis was also conducted with five arbitrarily chosen Ascomycetes-specific genes and five arbitrarily chosen common genes (genes that have a homolog in at least one non-Ascomycetes organism). Ascomycetes-specific genes appear to be diverging faster than common genes in the S.cerevisiae species, a situation that was previously described between different yeast species. Our results point to the strong contribution, during DNA sequence evolution, of allelic recombination besides nucleotide substitution.     

Toward an ecological synthesis: a case for habitat selection

Habitat selection, and its associated density and frequency-dependent evolution, has a profound influence on such vital phenomena as population regulation, species interactions, the assembly of ecological communities, and the origin and maintenance of biodiversity. Different strategies of habitat selection, and their importance in ecology and evolution, can often be revealed simply by plots of density in adjacent habitats. For individual species, the strategies are closely intertwined with mechanisms of population regulation, and with the persistence of populations through time. For interacting species, strategies of habitat selection are not only responsible for species coexistence, but provide one of the most convenient mechanisms for measuring competition, and the various community structures caused by competitive interactions. Other kinds of interactions, such as those between predators and prey, demonstrate that an understanding of the coevolution of habitat-selection strategies among strongly interacting species is essential to properly interpret their spatial and temporal dynamics. At the evolutionary scale, the frequency dependence associated with habitat selection may often allow populations to diverge and diversify into separate species. Habitat selection thereby demonstrates how we can map microevolutionary strategies in behavior onto their population and community consequences, and from there, onto macroevolutionary patterns of speciation and adaptive radiation. We can anticipate that future studies of habitat selection will not only help us complete those maps, but that they will also continue to enrich the panoply of ideas that shape evolutionary ecology.     

Antagonistic coevolution limits population persistence of a virus in a thermally deteriorating environment

Understanding the conditions under which rapid evolutionary adaptation can prevent population extinction in deteriorating environments (i.e. evolutionary rescue) is a crucial aim in the face of global climate change. Despite a rapidly growing body of work in this area, little attention has been paid to the importance of interspecific coevolutionary interactions. Antagonistic coevolution commonly observed between hosts and parasites is likely to retard evolutionary rescue because it often reduces population sizes, and results in the evolution of costly host defence and parasite counter-defence. We used experimental populations of a bacterium Pseudomonas fluorescens SBW25 and a bacteriophage virus (SBW25Φ2), to study how host-parasite coevolution impacts viral population persistence in the face of gradually increasing temperature, an environmental stress for the virus but not the bacterium. The virus persisted much longer when it evolved in the presence of an evolutionarily constant host genotype (i.e. in the absence of coevolution) than when the bacterium and virus coevolved. Further experiments suggest that both a reduction in population size and costly infectivity strategies contributed to viral extinction as a result of coevolution. The results highlight the importance of interspecific evolutionary interactions for the evolutionary responses of populations to global climate change.     

Probabilistic transition from unstable predator-prey interaction to stable coexistence of Dictyostelium discoideum and Escherichia coli

Predator-prey interactions have been found at all levels within ecosystems. Despite their ecological ubiquity and importance, the process of transition to a stable coexistent state has been poorly verified experimentally. To investigate the stabilization process of predator-prey interactions, we previously constructed a reproducible experimental predator-prey system between Dictyostelium discoideum and Escherichia coli, and showed that the phenotypically changed E. coli contributed to stabilization of the system. In the present study, we focused on the transition to stable coexistence of both species after the phenotypic change in E. coli. Analysis of E. coli cells isolated from co-culture plates as single colony enabled us to readily identify the appearance of phenotypically changed E. coli that differed in colony morphology and growth rate. It was also demonstrated that two types of viscous colony, i.e., the dense-type and sparse-type, differing in spatial distribution of both species emerged probabilistically and all of the viscous colonies maintained stably were of the sparse-type. These results suggest that the phenotypically changed E. coli may produce two types of viscous colonies probabilistically. The difference in spatial distribution would affect localized interactions between both species and then cause probabilistic stabilization of predator-prey interactions.     

Multiple Molecular Mechanisms Cause Reproductive Isolation between Three Yeast Species

Nuclear-mitochondrial conflict (cytonuclear incompatibility) is a specific form of Dobzhansky-Muller incompatibility previously shown to cause reproductive isolation in two yeast species. Here, we identified two new incompatible genes, MRS1 and AIM22, through a systematic study of F2 hybrid sterility caused by cytonuclear incompatibility in three closely related Saccharomyces species (S. cerevisiae, S. paradoxus, and S. bayanus). Mrs1 is a nuclear gene product required for splicing specific introns in the mitochondrial COX1, and Aim22 is a ligase encoded in the nucleus that is required for mitochondrial protein lipoylation. By comparing different species, our result suggests that the functional changes in MRS1 are a result of coevolution with changes in the COX1 introns. Further molecular analyses demonstrate that three nonsynonymous mutations are responsible for the functional differences of Mrs1 between these species. Functional complementation assays to determine when these incompatible genes altered their functions show a strong correlation between the sequence-based phylogeny and the evolution of cytonuclear incompatibility. Our results suggest that nuclear-mitochondrial incompatibility may represent a general mechanism of reproductive isolation during yeast evolution.     

Long-Term and Short-Term Evolutionary Impacts of Transposable Elements on Drosophila

Transposable elements (TEs) are considered to be genomic parasites and their interactions with their hosts have been likened to the coevolution between host and other nongenomic, horizontally transferred pathogens. TE families, however, are vertically inherited as integral segments of the nuclear genome. This transmission strategy has been suggested to weaken the selective benefits of host alleles repressing the transposition of specific TE variants. On the other hand, the elevated rates of TE transposition and high incidences of deleterious mutations observed during the rare cases of horizontal transfers of TE families between species could create at least a transient process analogous to the influence of horizontally transmitted pathogens. Here, we formally address this analogy, using empirical and theoretical analysis to specify the mechanism of how host–TE interactions may drive the evolution of host genes. We found that host TE-interacting genes actually have more pervasive evidence of adaptive evolution than immunity genes that interact with nongenomic pathogens in Drosophila. Yet, both our theoretical modeling and empirical observations comparing Drosophila melanogaster populations before and after the horizontal transfer of P elements, which invaded D. melanogaster early last century, demonstrated that horizontally transferred TEs have only a limited influence on host TE-interacting genes. We propose that the more prevalent and constant interaction with multiple vertically transmitted TE families may instead be the main force driving the fast evolution of TE-interacting genes, which is fundamentally different from the gene-for-gene interaction of host–pathogen coevolution.     

Additions, Losses, and Rearrangements on the Evolutionary Route from a Reconstructed Ancestor to the Modern Saccharomyces cerevisiae Genome

Comparative genomics can be used to infer the history of genomic rearrangements that occurred during the evolution of a species. We used the principle of parsimony, applied to aligned synteny blocks from 11 yeast species, to infer the gene content and gene order that existed in the genome of an extinct ancestral yeast about 100 Mya, immediately before it underwent whole-genome duplication (WGD). The reconstructed ancestral genome contains 4,703 ordered loci on eight chromosomes. The reconstruction is complete except for the subtelomeric regions. We then inferred the series of rearrangement steps that led from this ancestor to the current Saccharomyces cerevisiae genome; relative to the ancestral genome we observe 73 inversions, 66 reciprocal translocations, and five translocations involving telomeres. Some fragile chromosomal sites were reused as evolutionary breakpoints multiple times. We identified 124 genes that have been gained by S. cerevisiae in the time since the WGD, including one that is derived from a hAT family transposon, and 88 ancestral loci at which S. cerevisiae did not retain either of the gene copies that were formed by WGD. Sites of gene gain and evolutionary breakpoints both tend to be associated with tRNA genes and, to a lesser extent, with origins of replication. Many of the gained genes in S. cerevisiae have functions associated with ethanol production, growth in hypoxic environments, or the uptake of alternative nutrient sources.     

A genome-wide survey of changes in protein evolutionary rates across four closely related species of <Emphasis Type="Italic">Saccharomyces</Emphasis> sensu stricto group

Background  

Changes in protein evolutionary rates among lineages have been frequently observed during periods of notable phenotypic evolution. It is also known that, following gene duplication and loss, the protein evolutionary rates of genes involved in such events changed because of changes in functional constraints acting on the genes. However, in the evolution of closely related species, excluding the aforementioned situations, the frequency of changes in protein evolutionary rates is still not clear at the genome-wide level. Here we examine the constancy of protein evolutionary rates in the evolution of four closely related species of the Saccharomyces sensu stricto group (S. cerevisiae, S. paradoxus, S. mikatae and S. bayanus).