Genomic and phenotypic evolution of Escherichia coli in a novel citrate-only resource environment

Genomic sequencing reveals convergent adaptation during experimental evolution in two budding yeast species

Convergent evolution is central in the origins of multicellularity. Identifying the basis for convergent multicellular evolution is challenging because of the diverse evolutionary origins and environments involved. Haploid Kluyveromyces lactis populations evolve multicellularity during selection for increased settling in liquid media. Strong genomic and phenotypic convergence is observed between K. lactis and previously selected S. cerevisiae populations under similar selection, despite their >100-million-year divergence. We find K. lactis multicellularity is conferred by mutations in genes ACE2 or AIM44, with ACE2 being predominant. They are a subset of the six genes involved in the S. cerevisiae multicellularity. Both ACE2 and AIM44 regulate cell division, indicating that the genetic convergence is likely due to conserved cellular replication mechanisms. Complex population dynamics involving multiple ACE2/AIM44 genotypes are found in most K. lactis lineages. The results show common ancestry and natural selection shape convergence while chance and contingency determine the degree of divergence.

'Evolutionary poker': an agent-based model of interactome emergence and epistasis tested against Lenski's long-term E. coli experiments

A simple agent-based model is presented that produces results matching the experimental data found by Lenski's group for ≤50,000 generations of Escherichia coli bacteria under continuous selective pressure. Although various mathematical models have been devised previously to model the Lenski data, the present model has advantages in terms of overall simplicity and conceptual accessibility. The model also clearly illustrates a number of features of the evolutionary process that are otherwise not obvious, such as the roles of epistasis and historical contingency in adaptation and why evolution is time irreversible ('Dollo's law'). The reason for this irreversibility is that genomes become increasingly integrated or organized, and this organization becomes a novel selective factor itself, against which future generations must compete. Selection for integrated or synergistic networks, systems or sets of mutations or traits, not for individual mutations, confers the main adaptive advantage. The result is a punctuated form of evolution that follows a logarithmic occurrence probability, in which evolution proceeds very quickly when interactomes begin to form but which slows as interactomes become more robust and the difficulty of integrating new mutations increases. Sufficient parameters exist in the game to suggest not only how equilibrium or stasis is reached but also the conditions in which it will be punctuated, the factors governing the rate at which genomic organization occurs and novel traits appear, and how population size, genome size and gene variability affect these.

Large-scale duplication events underpin population-level flexibility in tRNA gene copy number in Pseudomonas fluorescens SBW25

The complement of tRNA genes within a genome is typically considered to be a (relatively) stable characteristic of an organism. Here, we demonstrate that bacterial tRNA gene set composition can be more flexible than previously appreciated, particularly regarding tRNA gene copy number. We report the high-rate occurrence of spontaneous, large-scale, tandem duplication events in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The identified duplications are up to ∼1 Mb in size (∼15% of the wildtype genome) and are predicted to change the copy number of up to 917 genes, including several tRNA genes. The observed duplications are inherently unstable: they occur, and are subsequently lost, at extremely high rates. We propose that this unusually plastic type of mutation provides a mechanism by which tRNA gene set diversity can be rapidly generated, while simultaneously preserving the underlying tRNA gene set in the absence of continued selection. That is, if a tRNA set variant provides no fitness advantage, then high-rate segregation of the duplication ensures the maintenance of the original tRNA gene set. However, if a tRNA gene set variant is beneficial, the underlying duplication fragment(s) may persist for longer and provide raw material for further, more stable, evolutionary change.

Duplicated antibiotic resistance genes reveal ongoing selection and horizontal gene transfer in bacteria

Horizontal gene transfer (HGT) and gene duplication are often considered as separate mechanisms driving the evolution of new functions. However, the mobile genetic elements (MGEs) implicated in HGT can copy themselves, so positive selection on MGEs could drive gene duplications. Here, we use a combination of modeling and experimental evolution to examine this hypothesis and use long-read genome sequences of tens of thousands of bacterial isolates to examine its generality in nature. Modeling and experiments show that antibiotic selection can drive the evolution of duplicated antibiotic resistance genes (ARGs) through MGE transposition. A key implication is that duplicated ARGs should be enriched in environments associated with antibiotic use. To test this, we examined the distribution of duplicated ARGs in 18,938 complete bacterial genomes with ecological metadata. Duplicated ARGs are highly enriched in bacteria isolated from humans and livestock. Duplicated ARGs are further enriched in an independent set of 321 antibiotic-resistant clinical isolates. Our findings indicate that duplicated genes often encode functions undergoing positive selection and horizontal gene transfer in microbial communities.

Evolution of a cross-feeding interaction following a key innovation in a long-term evolution experiment with Escherichia coli

The evolution of a novel trait can profoundly change an organism's effects on its environment, which can in turn affect the further evolution of that organism and any coexisting organisms. We examine these effects and feedbacks following the evolution of a novel function in the Long-Term Evolution Experiment (LTEE) with Escherichia coli. A characteristic feature of E. coli is its inability to grow aerobically on citrate (Cit-). Nonetheless, a Cit+ variant with this capacity evolved in one LTEE population after 31 000 generations. The Cit+ clade then coexisted stably with another clade that retained the ancestral Cit- phenotype. This coexistence was shaped by the evolution of a cross-feeding relationship based on C4-dicarboxylic acids, particularly succinate, fumarate, and malate, that the Cit+ variants release into the medium. Both the Cit- and Cit+ cells evolved to grow on these excreted resources. The evolution of aerobic growth on citrate thus led to a transition from an ecosystem based on a single limiting resource, glucose, to one with at least five resources that were either shared or partitioned between the two coexisting clades. Our findings show that evolutionary novelties can change environmental conditions in ways that facilitate diversity by altering ecosystem structure and the evolutionary trajectories of coexisting lineages.

Best Practices in Microbial Experimental Evolution: Using Reporters and Long-Read Sequencing to Identify Copy Number Variation in Experimental Evolution

Copy number variants (CNVs), comprising gene amplifications and deletions, are a pervasive class of heritable variation. CNVs play a key role in rapid adaptation in both natural, and experimental, evolution. However, despite the advent of new DNA sequencing technologies, detection and quantification of CNVs in heterogeneous populations has remained challenging. Here, we summarize recent advances in the use of CNV reporters that provide a facile means of quantifying de novo CNVs at a specific locus in the genome, and nanopore sequencing, for resolving the often complex structures of CNVs. We provide guidance for the engineering and analysis of CNV reporters and practical guidelines for single-cell analysis of CNVs using flow cytometry. We summarize recent advances in nanopore sequencing, discuss the utility of this technology, and provide guidance for the bioinformatic analysis of these data to define the molecular structure of CNVs. The combination of reporter systems for tracking and isolating CNV lineages and long-read DNA sequencing for characterizing CNV structures enables unprecedented resolution of the mechanisms by which CNVs are generated and their evolutionary dynamics.

Revisiting the Design of the Long-Term Evolution Experiment with Escherichia coli

The long-term evolution experiment (LTEE) with Escherichia coli began in 1988 and it continues to this day, with its 12 populations having recently reached 75,000 generations of evolution in a simple, well-controlled environment. The LTEE was designed to explore open-ended questions about the dynamics and repeatability of phenotypic and genetic evolution. Here I discuss various aspects of the LTEE's experimental design that have enabled its stability and success, including the choices of the culture regime, growth medium, ancestral strain, and statistical replication. I also discuss some of the challenges associated with a long-running project, such as handling procedural errors (e.g., cross-contamination) and managing the expanding collection of frozen samples. The simplicity of the experimental design and procedures have supported the long-term stability of the LTEE. That stability-along with the inherent creativity of the evolutionary process and the emergence of new genomic technologies-provides a platform that has allowed talented students and collaborators to pose questions, collect data, and make discoveries that go far beyond anything I could have imagined at the start of the LTEE.

The origin of genetic and metabolic systems: Evolutionary structuralinsights

DNA is derived from reverse transcription and its origin is related to reverse transcriptase, DNA polymerase and integrase. The gene structure originated from the evolution of the first RNA polymerase. Thus, an explanation of the origin of the genetic system must also explain the evolution of these enzymes. This paper proposes a polymer structure model, termed the stable complex evolution model, which explains the evolution of enzymes and functional molecules. Enzymes evolved their functions by forming locally tightly packed complexes with specific substrates. A metabolic reaction can therefore be considered to be the result of adaptive evolution in this way when a certain essential molecule is lacking in a cell. The evolution of the primitive genetic and metabolic systems was thus coordinated and synchronized. According to the stable complex model, almost all functional molecules establish binding affinity and specific recognition through complementary interactions, and functional molecules therefore have the nature of being auto-reactive. This is thermodynamically favorable and leads to functional duplication and self-organization. Therefore, it can be speculated that biological systems have a certain tendency to maintain functional stability or are influenced by an inherent selective power. The evolution of dormant bacteria may support this hypothesis, and inherent selectivity can be unified with natural selection at the molecular level.

Virus Evolution on Fitness Landscapes

The landscape paradigm is revisited in the light of evolution in simple systems. A brief overview of different classes of fitness landscapes is followed by a more detailed discussion of the RNA model, which is currently the only evolutionary model that allows for a comprehensive molecular analysis of a fitness landscape. Neutral networks of genotypes are indispensable for the success of evolution. Important insights into the evolutionary mechanism are gained by considering the topology of sequence and shape spaces. The dynamic concept of molecular quasispecies is viewed in the light of the landscape paradigm. The distribution of fitness values in state space is mirrored by the population structures of mutant distributions. Two classes of thresholds for replication error or mutations are important: (i) the-conventional-genotypic error threshold, which separates ordered replication from random drift on neutral networks, and (ii) a phenotypic error threshold above which the molecular phenotype is lost. Empirical landscapes are reviewed and finally, the implications of the landscape concept for virus evolution are discussed.

Extracellular Vesicles from 50,000 Generation Clones of the Escherichia coli Long-Term Evolution Experiment

Extracellular vesicles (EVs) are critical elements of cell-cell communication. Here, we characterized the outer membrane vesicles (OMVs) released by specific clones of Escherichia coli isolated from the Long-Term Evolution Experiment after 50,000 generations (50K) of adaptation to glucose minimal medium. Compared with their ancestor, the evolved clones produce small OMVs but also larger ones which display variable amounts of both OmpA and LPS. Tracking ancestral, fluorescently labelled OMVs revealed that they fuse with both ancestral- and 50K-evolved cells, albeit in different proportions. We quantified that less than 2% of the cells from one 50K-evolved clone acquired the fluorescence delivered by OMVs from the ancestral strain but that one cell concomitantly fuses with several OMVs. Globally, our results showed that OMV production in E. coli is a phenotype that varies along bacterial evolution and question the contribution of OMVs-mediated interactions in bacterial adaptation.

Idiosyncratic Purifying Selection on Metabolic Enzymes in the Long-Term Evolution Experiment with Escherichia coli

Bacteria, Archaea, and Eukarya all share a common set of metabolic reactions. This implies that the function and topology of central metabolism has been evolving under purifying selection over deep time. Central metabolism may similarly evolve under purifying selection during long-term evolution experiments, although it is unclear how long such experiments would have to run (decades, centuries, millennia) before signs of purifying selection on metabolism appear. I hypothesized that central and superessential metabolic enzymes would show evidence of purifying selection in the long-term evolution experiment with Escherichia coli (LTEE). I also hypothesized that enzymes that specialize on single substrates would show stronger evidence of purifying selection in the LTEE than generalist enzymes that catalyze multiple reactions. I tested these hypotheses by analyzing metagenomic time series covering 62,750 generations of the LTEE. I find mixed support for these hypotheses, because the observed patterns of purifying selection are idiosyncratic and population-specific. To explain this finding, I propose the Jenga hypothesis, named after a children's game in which blocks are removed from a tower until it falls. The Jenga hypothesis postulates that loss-of-function mutations degrade costly, redundant, and non-essential metabolic functions. Replicate populations can therefore follow idiosyncratic trajectories of lost redundancies, despite purifying selection on overall function. I tested the Jenga hypothesis by simulating the evolution of 1,000 minimal genomes under strong purifying selection. As predicted, the minimal genomes converge to different metabolic networks. Strikingly, the core genes common to all 1,000 minimal genomes show consistent signatures of purifying selection in the LTEE.

Adaptation dynamics between copy-number and point mutations

Together, copy-number and point mutations form the basis for most evolutionary novelty, through the process of gene duplication and divergence. While a plethora of genomic data reveals the long-term fate of diverging coding sequences and their cis-regulatory elements, little is known about the early dynamics around the duplication event itself. In microorganisms, selection for increased gene expression often drives the expansion of gene copy-number mutations, which serves as a crude adaptation, prior to divergence through refining point mutations. Using a simple synthetic genetic reporter system that can distinguish between copy-number and point mutations, we study their early and transient adaptive dynamics in real time in Escherichia coli. We find two qualitatively different routes of adaptation, depending on the level of functional improvement needed. In conditions of high gene expression demand, the two mutation types occur as a combination. However, under low gene expression demand, copy-number and point mutations are mutually exclusive; here, owing to their higher frequency, adaptation is dominated by copy-number mutations, in a process we term amplification hindrance. Ultimately, due to high reversal rates and pleiotropic cost, copy-number mutations may not only serve as a crude and transient adaptation, but also constrain sequence divergence over evolutionary time scales.

The landscape of innovation in bacteria, battleships, and beyond

We draw lessons from microbial experimental evolution and naval warfare to improve the understanding of innovation in financial markets. Major financial innovations often arise without explicit societal planning because novel approaches can be favored by markets, in a manner strikingly parallel to natural selection. We utilize the concept of an adaptive landscape to characterize environments that increase the speed and magnitude of innovation. We apply this adaptive landscape framework to innovation in portfolio management. We create a general taxonomy for understanding and nurturing innovation.

Divergent Evolution of Mutation Rates and Biases in the Long-Term Evolution Experiment with Escherichia coli

All organisms encode enzymes that replicate, maintain, pack, recombine, and repair their genetic material. For this reason, mutation rates and biases also evolve by mutation, variation, and natural selection. By examining metagenomic time series of the Lenski long-term evolution experiment (LTEE) with Escherichia coli (Good BH, McDonald MJ, Barrick JE, Lenski RE, Desai MM. 2017. The dynamics of molecular evolution over 60,000 generations. Nature 551(7678):45–50.), we find that local mutation rate variation has evolved during the LTEE. Each LTEE population has evolved idiosyncratic differences in their rates of point mutations, indels, and mobile element insertions, due to the fixation of various hypermutator and antimutator alleles. One LTEE population, called Ara+3, shows a strong, symmetric wave pattern in its density of point mutations, radiating from the origin of replication. This pattern is largely missing from the other LTEE populations, most of which evolved missense, indel, or structural mutations in topA, fis, and dusB—loci that all affect DNA topology. The distribution of mutations in those genes over time suggests epistasis and historical contingency in the evolution of DNA topology, which may have in turn affected local mutation rates. Overall, the replicate populations of the LTEE have largely diverged in their mutation rates and biases, even though they have adapted to identical abiotic conditions.

Genomic evolution of antibiotic resistance is contingent on genetic background following a long-term experiment with Escherichia coli

Antibiotic resistance is a growing health concern. Efforts to control resistance would benefit from an improved ability to forecast when and how it will evolve. Epistatic interactions between mutations can promote divergent evolutionary trajectories, which complicates our ability to predict evolution. We recently showed that differences between genetic backgrounds can lead to idiosyncratic responses in the evolvability of phenotypic resistance, even among closely related Escherichia coli strains. In this study, we examined whether a strain's genetic background also influences the genotypic evolution of resistance. Do lineages founded by different genotypes take parallel or divergent mutational paths to achieve their evolved resistance states? We addressed this question by sequencing the complete genomes of antibiotic-resistant clones that evolved from several different genetic starting points during our earlier experiments. We first validated our statistical approach by quantifying the specificity of genomic evolution with respect to antibiotic treatment. As expected, mutations in particular genes were strongly associated with each drug. Then, we determined that replicate lines evolved from the same founding genotypes had more parallel mutations at the gene level than lines evolved from different founding genotypes, although these effects were more subtle than those showing antibiotic specificity. Taken together with our previous work, we conclude that historical contingency can alter both genotypic and phenotypic pathways to antibiotic resistance.

Changes in Cell Size and Shape during 50,000 Generations of Experimental Evolution with Escherichia coli

Bacteria adopt a wide variety of sizes and shapes, with many species exhibiting stereotypical morphologies. How morphology changes, and over what timescales, is less clear. Previous work examining cell morphology in an experiment with Escherichia coli showed that populations evolved larger cells and, in some cases, cells that were less rod-like. That experiment has now run for over two more decades. Meanwhile, genome sequence data are available for these populations, and new computational methods enable high-throughput microscopic analyses. In this study, we measured stationary-phase cell volumes for the ancestor and 12 populations at 2,000, 10,000, and 50,000 generations, including measurements during exponential growth at the last time point. We measured the distribution of cell volumes for each sample using a Coulter counter and microscopy, the latter of which also provided data on cell shape. Our data confirm the trend toward larger cells while also revealing substantial variation in size and shape across replicate populations. Most populations first evolved wider cells but later reverted to the ancestral length-to-width ratio. All but one population evolved mutations in rod shape maintenance genes. We also observed many ghost-like cells in the only population that evolved the novel ability to grow on citrate, supporting the hypothesis that this lineage struggles with maintaining balanced growth. Lastly, we show that cell size and fitness remain correlated across 50,000 generations. Our results suggest that larger cells are beneficial in the experimental environment, while the reversion toward ancestral length-to-width ratios suggests partial compensation for the less favorable surface area-to-volume ratios of the evolved cells.IMPORTANCE Bacteria exhibit great morphological diversity, yet we have only a limited understanding of how their cell sizes and shapes evolve and of how these features affect organismal fitness. This knowledge gap reflects, in part, the paucity of the fossil record for bacteria. In this study, we revived and analyzed samples extending over 50,000 generations from 12 populations of experimentally evolving Escherichia coli to investigate the relation between cell size, shape, and fitness. Using this "frozen fossil record," we show that all 12 populations evolved larger cells concomitant with increased fitness, with substantial heterogeneity in cell size and shape across the replicate lines. Our work demonstrates that cell morphology can readily evolve and diversify, even among populations living in identical environments.

Challenges Faced by Highly Polyploid Bacteria with Limits on DNA Inheritance

Most studies of bacterial reproduction have centered on organisms that undergo binary fission. In these models, complete chromosome copies are segregated with great fidelity into two equivalent offspring cells. All genetic material is passed on to offspring, including new mutations and horizontally acquired sequences. However, some bacterial lineages employ diverse reproductive patterns that require management and segregation of more than two chromosome copies. Epulopiscium spp., and their close relatives within the Firmicutes phylum, are intestinal symbionts of surgeonfish (family Acanthuridae). Each of these giant (up to 0.6 mm long), cigar-shaped bacteria contains tens of thousands of chromosome copies. Epulopiscium spp. do not use binary fission but instead produce multiple intracellular offspring. Only ∼1% of the genetic material in an Epulopiscium sp. type B mother cell is directly inherited by its offspring cells. And yet, even in late stages of offspring development, mother-cell chromosome copies continue to replicate. Consequently, chromosomes take on a somatic or germline role. Epulopiscium sp. type B is a strict anaerobe and while it is an obligate symbiont, its host has a facultative association with this intestinal microorganism. Therefore, Epulopiscium sp. type B populations face several bottlenecks that could endanger their diversity and resilience. Multilocus sequence analyses revealed that recombination is important to diversification in populations of Epulopiscium sp. type B. By employing mechanisms common to others in the Firmicutes, the coordinated timing of mother-cell lysis, offspring development and congression may facilitate the substantial recombination observed in Epulopiscium sp. type B populations.

Evolutionary dynamics and structural consequences of de novo beneficial mutations and mutant lineages arising in a constant environment

BACKGROUND

Microbial evolution experiments can be used to study the tempo and dynamics of evolutionary change in asexual populations, founded from single clones and growing into large populations with multiple clonal lineages. High-throughput sequencing can be used to catalog de novo mutations as potential targets of selection, determine in which lineages they arise, and track the fates of those lineages. Here, we describe a long-term experimental evolution study to identify targets of selection and to determine when, where, and how often those targets are hit.

RESULTS

We experimentally evolved replicate Escherichia coli populations that originated from a mutator/nonsense suppressor ancestor under glucose limitation for between 300 and 500 generations. Whole-genome, whole-population sequencing enabled us to catalog 3346 de novo mutations that reached > 1% frequency. We sequenced the genomes of 96 clones from each population when allelic diversity was greatest in order to establish whether mutations were in the same or different lineages and to depict lineage dynamics. Operon-specific mutations that enhance glucose uptake were the first to rise to high frequency, followed by global regulatory mutations. Mutations related to energy conservation, membrane biogenesis, and mitigating the impact of nonsense mutations, both ancestral and derived, arose later. New alleles were confined to relatively few loci, with many instances of identical mutations arising independently in multiple lineages, among and within replicate populations. However, most never exceeded 10% in frequency and were at a lower frequency at the end of the experiment than at their maxima, indicating clonal interference. Many alleles mapped to key structures within the proteins that they mutated, providing insight into their functional consequences.

CONCLUSIONS

Overall, we find that when mutational input is increased by an ancestral defect in DNA repair, the spectrum of high-frequency beneficial mutations in a simple, constant resource-limited environment is narrow, resulting in extreme parallelism where many adaptive mutations arise but few ever go to fixation.

Prokaryotic life finds a way: insights from evolutionary experimentation in bacteria

While evolution proceeds through the generation of random variant alleles, the application of selective pressures can select for subsets of mutations that confer fitness-improving physiological benefits. This, in essence, defines the process of adaptive evolution. The rapid replication rate of bacteria has allowed for the design of experiments to study these processes over a reasonable timeframe within a laboratory setting. This has been greatly assisted by advances in tractability of diverse microorganisms, next generation sequencing technologies and bioinformatic analysis pipelines. Examining the processes by which organisms adapt their genetic code to cope with sub-optimal growth conditions has yielded a wealth of molecular insight into diverse biological processes. Here we discuss how the study of adaptive evolutionary trajectories in bacteria has allowed for improved understanding of stress responses, revealed important insight into microbial physiology, allowed for the production of highly optimised strains for use in biotechnology and increased our knowledge of the role of genomic plasticity in chronic infections.

Genomic and phenotypic evolution of Escherichia coli in a novel citrate-only resource environment

Evolutionary innovations allow populations to colonize new ecological niches. We previously reported that aerobic growth on citrate (Cit⁺) evolved in an Escherichia coli population during adaptation to a minimal glucose medium containing citrate (DM25). Cit⁺ variants can also grow in citrate-only medium (DM0), a novel environment for E. coli. To study adaptation to this niche, we founded two sets of Cit⁺ populations and evolved them for 2500 generations in DM0 or DM25. The evolved lineages acquired numerous parallel mutations, many mediated by transposable elements. Several also evolved amplifications of regions containing the maeA gene. Unexpectedly, some evolved populations and clones show apparent declines in fitness. We also found evidence of substantial cell death in Cit⁺ clones. Our results thus demonstrate rapid trait refinement and adaptation to the new citrate niche, while also suggesting a recalcitrant mismatch between E. coli physiology and growth on citrate.