Heterogeneous adaptive trajectories of small populations on complex fitness landscapes

Bacteria evolve macroscopic multicellularity by the genetic assimilation of phenotypically plastic cell clustering

The evolutionary transition from unicellularity to multicellularity was a key innovation in the history of life. Experimental evolution is an important tool to study the formation of undifferentiated cellular clusters, the likely first step of this transition. Although multicellularity first evolved in bacteria, previous experimental evolution research has primarily used eukaryotes. Moreover, it focuses on mutationally driven (and not environmentally induced) phenotypes. Here we show that both Gram-negative and Gram-positive bacteria exhibit phenotypically plastic (i.e., environmentally induced) cell clustering. Under high salinity, they form elongated clusters of ~ 2 cm. However, under habitual salinity, the clusters disintegrate and grow planktonically. We used experimental evolution with Escherichia coli to show that such clustering can be assimilated genetically: the evolved bacteria inherently grow as macroscopic multicellular clusters, even without environmental induction. Highly parallel mutations in genes linked to cell wall assembly formed the genomic basis of assimilated multicellularity. While the wildtype also showed cell shape plasticity across high versus low salinity, it was either assimilated or reversed after evolution. Interestingly, a single mutation could genetically assimilate multicellularity by modulating plasticity at multiple levels of organization. Taken together, we show that phenotypic plasticity can prime bacteria for evolving undifferentiated macroscopic multicellularity.

Impact of population size on early adaptation in rugged fitness landscapes

Owing to stochastic fluctuations arising from finite population size, known as genetic drift, the ability of a population to explore a rugged fitness landscape depends on its size. In the weak mutation regime, while the mean steady-state fitness increases with population size, we find that the height of the first fitness peak encountered when starting from a random genotype displays various behaviours versus population size, even among small and simple rugged landscapes. We show that the accessibility of the different fitness peaks is key to determining whether this height overall increases or decreases with population size. Furthermore, there is often a finite population size that maximizes the height of the first fitness peak encountered when starting from a random genotype. This holds across various classes of model rugged landscapes with sparse peaks, and in some experimental and experimentally inspired ones. Thus, early adaptation in rugged fitness landscapes can be more efficient and predictable for relatively small population sizes than in the large-size limit. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Experimental and In Silico Analysis of TEM β-Lactamase Adaptive Evolution

Multiple mutations often have non-additive (epistatic) phenotypic effects. Epistasis is of fundamental biological relevance but is not well understood mechanistically. Adaptive evolution, i.e., the evolution of new biochemical activities, is rich in epistatic interactions. To better understand the principles underlying epistasis during genetic adaptation, we studied the evolution of TEM-1 β-lactamase variants exhibiting cefotaxime resistance. We report the collection of a library of 487 observed evolutionary trajectories for TEM-1 and determine the epistasis status based on cefotaxime resistance phenotype for 206 combinations of 2-3 TEM-1 mutations involving 17 positions under adaptive selective pressure. Gain-of-function (GOF) mutations are gatekeepers for adaptation. To see if GOF phenotypes can be inferred based solely on sequence data, we calculated the enrichment of GOF mutations in the different categories of epistatic pairs. Our results suggest that this is possible because GOF mutations are particularly enriched in sign and reciprocal sign epistasis, which leave a major imprint on the sequence space accessible to evolution. We also used FoldX to explore the relationship between thermodynamic stability and epistasis. We found that mutations in observed evolutionary trajectories tend to destabilize the folded structure of the protein, albeit their cumulative effects are consistently below the protein's free energy of folding. The destabilizing effect is stronger for epistatic pairs, suggesting that modest or local alterations in folding stability can modulate catalysis. Finally, we report a significant relationship between epistasis and the degree to which two protein positions are structurally and dynamically coupled, even in the absence of ligand.

Biofilm antimicrobial susceptibility through an experimental evolutionary lens

Experimental evolution experiments in which bacterial populations are repeatedly exposed to an antimicrobial treatment, and examination of the genotype and phenotype of the resulting evolved bacteria, can help shed light on mechanisms behind reduced susceptibility. In this review we present an overview of why it is important to include biofilms in experimental evolution, which approaches are available to study experimental evolution in biofilms and what experimental evolution has taught us about tolerance and resistance in biofilms. Finally, we present an emerging consensus view on biofilm antimicrobial susceptibility supported by data obtained during experimental evolution studies.

Fitness of evolving bacterial populations is contingent on deep and shallow history but only shallow history creates predictable patterns

Long-term evolution experiments have tested the importance of genetic and environmental factors in influencing evolutionary outcomes. Differences in phylogenetic history, recent adaptation to distinct environments and chance events, all influence the fitness of a population. However, the interplay of these factors on a population's evolutionary potential remains relatively unexplored. We tracked the outcome of 2000 generations of evolution of four natural isolates of Escherichia coli bacteria that were engineered to also create differences in shallow history by adding previously identified mutations selected in a separate long-term experiment. Replicate populations started from each progenitor evolved in four environments. We found that deep and shallow phylogenetic histories both contributed significantly to differences in evolved fitness, though by different amounts in different selection environments. With one exception, chance effects were not significant. Whereas the effect of deep history did not follow any detectable pattern, effects of shallow history followed a pattern of diminishing returns whereby fitter ancestors had smaller fitness increases. These results are consistent with adaptive evolution being contingent on the interaction of several evolutionary forces but demonstrate that the nature of these interactions is not fixed and may not be predictable even when the role of chance is small.

Population bottleneck has only marginal effect on fitness evolution and its repeatability in dioecious Caenorhabditis elegans

The predictability of evolution is expected to depend on the relative contribution of deterministic and stochastic processes. This ratio is modulated by effective population size. Smaller effective populations harbor less genetic diversity and stochastic processes are generally expected to play a larger role, leading to less repeatable evolutionary trajectories. Empirical insight into the relationship between effective population size and repeatability is limited and focused mostly on asexual organisms. Here, we tested whether fitness evolution was less repeatable after a population bottleneck in obligately outcrossing populations of Caenorhabditis elegans. Replicated populations founded by 500, 50, or five individuals (no/moderate/strong bottleneck) were exposed to a novel environment with a different bacterial prey. As a proxy for fitness, population size was measured after one week of growth before and after 15 weeks of evolution. Surprisingly, we found no significant differences among treatments in their fitness evolution. Even though the strong bottleneck reduced the relative contribution of selection to fitness variation, this did not translate to a significant reduction in the repeatability of fitness evolution. Thus, although a bottleneck reduced the contribution of deterministic processes, we conclude that the predictability of evolution may not universally depend on effective population size, especially in sexual organisms.

Population size mediates the contribution of high-rate and large-benefit mutations to parallel evolution

Mutations with large fitness benefits and mutations occurring at high rates may both cause parallel evolution, but their contribution is predicted to depend on population size. Moreover, high-rate and large-benefit mutations may have different long-term adaptive consequences. We show that small and 100-fold larger bacterial populations evolve resistance to a β-lactam antibiotic by using similar numbers, but different types of mutations. Small populations frequently substitute similar high-rate structural variants and loss-of-function point mutations, including the deletion of a low-activity β-lactamase, and evolve modest resistance levels. Large populations more often use low-rate, large-benefit point mutations affecting the same targets, including mutations activating the β-lactamase and other gain-of-function mutations, leading to much higher resistance levels. Our results demonstrate the separation by clonal interference of mutation classes with divergent adaptive consequences, causing a shift from high-rate to large-benefit mutations with increases in population size.

Extreme value theory as a framework for understanding mutation frequency distribution in cancer genomes

Currently, the population dynamics of preclonal cancer cells before clonal expansion of tumors has not been sufficiently addressed thus far. By focusing on preclonal cancer cell population as a Darwinian evolutionary system, we formulated and analyzed the observed mutation frequency among tumors (MFaT) as a proxy for the hypothesized sequence read frequency and beneficial fitness effect of a cancer driver mutation. Analogous to intestinal crypts, we assumed that sample donor patients are separate culture tanks where proliferating cells follow certain population dynamics described by extreme value theory (EVT). To validate this, we analyzed three large-scale cancer genome datasets, each harboring > 10000 tumor samples and in total involving > 177898 observed mutation sites. We clarified the necessary premises for the application of EVT in the strong selection and weak mutation (SSWM) regime in relation to cancer genome sequences at scale. We also confirmed that the stochastic distribution of MFaT is likely of the Fréchet type, which challenges the well-known Gumbel hypothesis of beneficial fitness effects. Based on statistical data analysis, we demonstrated the potential of EVT as a population genetics framework to understand and explain the stochastic behavior of driver-mutation frequency in cancer genomes as well as its applicability in real cancer genome sequence data.

Compensatory adaptation and diversification subsequent to evolutionary rescue in a model adaptive radiation

Biological populations may survive lethal environmental stress through evolutionary rescue. The rescued populations typically suffer a reduction in growth performance and harbor very low genetic diversity compared with their parental populations. The present study addresses how population size and within-population diversity may recover through compensatory evolution, using the experimental adaptive radiation of bacterium Pseudomonas fluorescens. We exposed bacterial populations to an antibiotic treatment and then imposed a one-individual-size population bottleneck on those surviving the antibiotic stress. During the subsequent compensatory evolution, population size increased and leveled off very rapidly. The increase of diversity was of slower paces and persisted longer. In the very early stage of compensatory evolution, populations of large sizes had a greater chance to diversify; however, this productivity-diversification relationship was not observed in later stages. Population size and diversity from the end of the compensatory evolution was not contingent on initial population growth performance. We discussed the possibility that our results be explained by the emergence of a "holey" fitness landscape under the antibiotic stress.

Population Bottlenecks Strongly Affect the Evolutionary Dynamics of Antibiotic Persistence

Bacterial persistence is a potential cause of antibiotic therapy failure. Antibiotic-tolerant persisters originate from phenotypic differentiation within a susceptible population, occurring with a frequency that can be altered by mutations. Recent studies have proven that persistence is a highly evolvable trait and, consequently, an important evolutionary strategy of bacterial populations to adapt to high-dose antibiotic therapy. Yet, the factors that govern the evolutionary dynamics of persistence are currently poorly understood. Theoretical studies predict far-reaching effects of bottlenecking on the evolutionary adaption of bacterial populations, but these effects have never been investigated in the context of persistence. Bottlenecking events are frequently encountered by infecting pathogens during host-to-host transmission and antibiotic treatment. In this study, we used a combination of experimental evolution and barcoded knockout libraries to examine how population bottlenecking affects the evolutionary dynamics of persistence. In accordance with existing hypotheses, small bottlenecks were found to restrict the adaptive potential of populations and result in more heterogeneous evolutionary outcomes. Evolutionary trajectories followed in small-bottlenecking regimes additionally suggest that the fitness landscape associated with persistence has a rugged topography, with distinct trajectories toward increased persistence that are accessible to evolving populations. Furthermore, sequencing data of evolved populations and knockout libraries after selection reveal various genes that are potentially involved in persistence, including previously known as well as novel targets. Together, our results do not only provide experimental evidence for evolutionary theories, but also contribute to a better understanding of the environmental and genetic factors that guide bacterial adaptation to antibiotic treatment.

Robustness and predictability of evolution in bottlenecked populations

Deterministic and stochastic evolutionary processes drive adaptation in natural populations. The strength of each component process is determined by the population size: deterministic components prevail in very large populations, while stochastic components are the driving mechanisms in small ones. Many natural populations, however, experience intermittent periods of growth, moving through states in which either stochastic or deterministic processes prevail. This growth is often countered by population bottlenecks, which abound in both natural and laboratory populations. Here we investigate how population bottlenecks shape the process of adaptation. We demonstrate that adaptive trajectories in populations experiencing regular bottlenecks can be naturally scaled in time units of generations; with this scaling the time courses of adaptation, fitness variance, and genetic diversity all become relatively insensitive to the timing of population bottlenecks, provided the bottleneck size exceeds a few thousand individuals. We also include analyses at the genotype level to investigate the impact of population bottlenecks on the predictability and distribution of evolutionary pathways. Irrespective of the timing of population bottlenecks, we find that predictability increases with population size. We also find that predictability of the adaptive pathways increases in increasingly rugged fitness landscapes. Overall, our work reveals that both the adaptation rate and the predictability of evolutionary trajectories are relatively robust to population bottlenecks.

Uncovering patterns of the evolution of genomic sequence entropy and complexity

The lack of consensus concerning the biological meaning of entropy and complexity of genomes and the different ways to assess these data hamper conclusions concerning what are the causes of genomic entropy variation among species. This study aims to evaluate the entropy and complexity of genomic sequences of several species without using homologies to assess relationships among these variables and non-molecular data (e.g., the number of individuals) to seek a trigger of interspecific genomic entropy variation. The results indicate a relationship among genomic entropy, genome size, genomic complexity, and the number of individuals: species with a small number of individuals harbors large genome, and hence, low entropy but a higher complexity. We defined that the complexity of a genome relies on the entropy of each DNA segment within genome. Then, the entropy and complexity of a genome reflects its organization solely. Exons of vertebrates harbor smaller entropies than non-exon regions (likely by the repeats that accumulated from duplications), whereas other taxonomic groups do not present this pattern. Our findings suggest that small initial population might have defined current genomic entropy and complexity: actual genomes are less complex than ancestral ones. Besides, our data disagree with the relationship between phenotype and genomic entropies previously established. Finally, by establishing the relationship between genomic entropy/complexity with the number of individuals and genome size, under an evolutive perspective, ideas concerning the genomic variability may emerge.

Quantum aspects of evolution: a contribution towards evolutionary explorations of genotype networks via quantum walks

Quantum biology seeks to explain biological phenomena via quantum mechanisms, such as enzyme reaction rates via tunnelling and photosynthesis energy efficiency via coherent superposition of states. However, less effort has been devoted to study the role of quantum mechanisms in biological evolution. In this paper, we used transcription factor networks with two and four different phenotypes, and used classical random walks (CRW) and quantum walks (QW) to compare network search behaviour and efficiency at finding novel phenotypes between CRW and QW. In the network with two phenotypes, at temporal scales comparable to decoherence time TD, QW are as efficient as CRW at finding new phenotypes. In the case of the network with four phenotypes, the QW had a higher probability of mutating to a novel phenotype than the CRW, regardless of the number of mutational steps (i.e. 1, 2 or 3) away from the new phenotype. Before quantum decoherence, the QW probabilities become higher turning the QW effectively more efficient than CRW at finding novel phenotypes under different starting conditions. Thus, our results warrant further exploration of the QW under more realistic network scenarios (i.e. larger genotype networks) in both closed and open systems (e.g. by considering Lindblad terms).

The Ingenuity of Bacterial Genomes

The genomes of bacteria contain fewer genes and substantially less noncoding DNA than those of eukaryotes, and as a result, they have much less raw material to invent new traits. Yet, bacteria are vastly more taxonomically diverse, numerically abundant, and globally successful in colonizing new habitats compared to eukaryotes. Although bacterial genomes are generally considered to be optimized for efficient growth and rapid adaptation, nonadaptive processes have played a major role in shaping the size, contents, and compact organization of bacterial genomes and have allowed the establishment of deleterious traits that serve as the raw materials for genetic innovation.

What lies beneath? Population dynamics conceal pace-of-life and sex ratio variation, with implications for resilience to environmental change

Life-history and pace-of-life syndrome theory predict that populations are comprised of individuals exhibiting different reproductive schedules and associated behavioural and physiological traits, optimized to prevailing social and environmental factors. Changing weather and social conditions provide in situ cues altering this life-history optimality; nevertheless, few studies have considered how tactical, sex-specific plasticity over an individual's lifespan varies in wild populations and influences population resilience. We examined the drivers of individual life-history schedules using 31 years of trapping data and 28 years of pedigree for the European badger (Meles meles L.), a long-lived, iteroparous, polygynandrous mammal that exhibits heterochrony in the timing of endocrinological puberty in male cubs. Our top model for the effects of environmental (social and weather) conditions during a badger's first year on pace-of-life explained <10% of variance in the ratio of fertility to age at first reproduction (F/α) and lifetime reproductive success. Conversely, sex ratio (SR) and sex-specific density explained 52.8% (males) and 91.0% (females) of variance in adult F/α ratios relative to the long-term population median F/α. Weather primarily affected the sexes at different life-history stages, with energy constraints limiting the onset of male reproduction but playing a large role in female strategic energy allocation, particularly in relation to ongoing mean temperature increases. Furthermore, the effects of social factors on age of first reproduction and year-to-year reproductive success covaried differently with sex, likely due to sex-specific responses to potential mate availability. For females, low same-sex densities favoured early primiparity; for males, instead, up to 10% of yearlings successfully mated at high same-sex densities. We observed substantial SR dynamism relating to differential mortality of life-history strategists within the population, and propose that shifting ratios of 'fast' and 'slow' life-history strategists contribute substantially to population dynamics and resilience to changing conditions.

Adaptive walks on high-dimensional fitness landscapes and seascapes with distance-dependent statistics

The dynamics of evolution is intimately shaped by epistasis - interactions between genetic elements which cause the fitness-effect of combinations of mutations to be non-additive. Analyzing evolutionary dynamics that involves large numbers of epistatic mutations is intrinsically difficult. A crucial feature is that the fitness landscape in the vicinity of the current genome depends on the evolutionary history. A key step is thus developing models that enable study of the effects of past evolution on future evolution. In this work, we introduce a broad class of high-dimensional random fitness landscapes for which the correlations between fitnesses of genomes are a general function of genetic distance. Their Gaussian character allows for tractable computational as well as analytic understanding. We study the properties of these landscapes focusing on the simplest evolutionary process: random adaptive (uphill) walks. Conventional measures of "ruggedness" are shown to not much affect such adaptive walks. Instead, the long-distance statistics of epistasis cause all properties to be highly conditional on past evolution, determining the statistics of the local landscape (the distribution of fitness-effects of available mutations and combinations of these), as well as the global geometry of evolutionary trajectories. In order to further explore the effects of conditioning on past evolution, we model the effects of slowly changing environments. At long times, such fitness "seascapes" cause a statistical steady state with highly intermittent evolutionary dynamics: populations undergo bursts of rapid adaptation, interspersed with periods in which adaptive mutations are rare and the population waits for more new directions to be opened up by changes in the environment. Finally, we discuss prospects for studying more complex evolutionary dynamics and on broader classes of high-dimensional landscapes and seascapes.

The Relationship between Spatial Structure and the Maintenance of Diversity in Microbial Populations

Spatial structure is pervasive in the microbial world, yet we know little about how it influences the evolution of microbial populations. It is thought that spatial structure limits the scale of competitive interactions and protracts selective sweeps. This may allow microbial populations to simultaneously explore multiple evolutionary paths. But how structured a microbial population must be before this effect is realized is not known. We used empirical and simulation studies to explore the relationship between spatial structure and the maintenance of diversity. The degree of spatial structure experienced by Escherichia coli metapopulations was manipulated by varying the migration rate between its component subpopulations. Each subpopulation was inoculated with an equal number of two equally fit genotypes, and their frequencies in 12 subpopulations were determined during 150 generations of evolution. We observed that the frequency of the "loser" genotypes decreased exponentially as the migration rate between the subpopulations was increased and that higher frequencies of the loser genotypes were maintained in structured metapopulations. These results demonstrate that structured microbial populations can evolve along multiple evolutionary trajectories even when migration rates between the subpopulations are relatively high.

On the deformability of an empirical fitness landscape by microbial evolution

A fitness landscape is a map between the genotype and its reproductive success in a given environment. The topography of fitness landscapes largely governs adaptive dynamics, constraining evolutionary trajectories and the predictability of evolution. Theory suggests that this topography can be deformed by mutations that produce substantial changes to the environment. Despite its importance, the deformability of fitness landscapes has not been systematically studied beyond abstract models, and little is known about its reach and consequences in empirical systems. Here we have systematically characterized the deformability of the genome-wide metabolic fitness landscape of the bacterium Escherichia coli Deformability is quantified by the noncommutativity of epistatic interactions, which we experimentally demonstrate in mutant strains on the path to an evolutionary innovation. Our analysis shows that the deformation of fitness landscapes by metabolic mutations rarely affects evolutionary trajectories in the short range. However, mutations with large environmental effects produce long-range landscape deformations in distant regions of the genotype space that affect the fitness of later descendants. Our results therefore suggest that, even in situations in which mutations have strong environmental effects, fitness landscapes may retain their power to forecast evolution over small mutational distances despite the potential attenuation of that power over longer evolutionary trajectories. Our methods and results provide an avenue for integrating adaptive and eco-evolutionary dynamics with complex genetics and genomics.

Fitness Landscape of the Immune Compromised Favors the Emergence of Antibiotic Resistance

Antibiotic resistance can come at a high cost, both in terms of fitness for the pathogen and poorer outcomes for patients. The fitness landscape encountered by bacterial pathogens varies greatly throughout patient populations in terms of host immunity as well as the duration and spectrum of antibiotics encountered. Severely immunocompromised patients present a favorable environment for antibiotic resistance to emerge due to lack of immune-mediated competition and increased opportunities to evolve both on-target and compensatory mutations. Such patients may present unique pathways for antibiotic resistance to emerge.

Environmental pleiotropy and demographic history direct adaptation under antibiotic selection

Evolutionary rescue following environmental change requires mutations permitting population growth in the new environment. If change is severe enough to prevent most of the population reproducing, rescue becomes reliant on mutations already present. If change is sustained, the fitness effects in both environments, and how they are associated—termed ‘environmental pleiotropy’—may determine which alleles are ultimately favoured. A population’s demographic history—its size over time—influences the variation present. Although demographic history is known to affect the probability of evolutionary rescue, how it interacts with environmental pleiotropy during severe and sustained environmental change remains unexplored. Here, we demonstrate how these factors interact during antibiotic resistance evolution, a key example of evolutionary rescue fuelled by pre-existing mutations with pleiotropic fitness effects. We combine published data with novel simulations to characterise environmental pleiotropy and its effects on resistance evolution under different demographic histories. Comparisons among resistance alleles typically revealed no correlation for fitness—i.e., neutral pleiotropy—above and below the sensitive strain’s minimum inhibitory concentration. Resistance allele frequency following experimental evolution showed opposing correlations with their fitness effects in the presence and absence of antibiotic. Simulations demonstrated that effects of environmental pleiotropy on allele frequencies depended on demographic history. At the population level, the major influence of environmental pleiotropy was on mean fitness, rather than the probability of evolutionary rescue or diversity. Our work suggests that determining both environmental pleiotropy and demographic history is critical for predicting resistance evolution, and we discuss the practicalities of this during in vivo evolution.

Experimental Design, Population Dynamics, and Diversity in Microbial Experimental Evolution

In experimental evolution, laboratory-controlled conditions select for the adaptation of species, which can be monitored in real time. Despite the current popularity of such experiments, nature's most pervasive biological force was long believed to be observable only on time scales that transcend a researcher's life-span, and studying evolution by natural selection was therefore carried out solely by comparative means. Eventually, microorganisms' propensity for fast evolutionary changes proved us wrong, displaying strong evolutionary adaptations over a limited time, nowadays massively exploited in laboratory evolution experiments. Here, we formulate a guide to experimental evolution with microorganisms, explaining experimental design and discussing evolutionary dynamics and outcomes and how it is used to assess ecoevolutionary theories, improve industrially important traits, and untangle complex phenotypes. Specifically, we give a comprehensive overview of the setups used in experimental evolution. Additionally, we address population dynamics and genetic or phenotypic diversity during evolution experiments and expand upon contributing factors, such as epistasis and the consequences of (a)sexual reproduction. Dynamics and outcomes of evolution are most profoundly affected by the spatiotemporal nature of the selective environment, where changing environments might lead to generalists and structured environments could foster diversity, aided by, for example, clonal interference and negative frequency-dependent selection. We conclude with future perspectives, with an emphasis on possibilities offered by fast-paced technological progress. This work is meant to serve as an introduction to those new to the field of experimental evolution, as a guide to the budding experimentalist, and as a reference work to the seasoned expert.

Parallel Evolution of High-Level Aminoglycoside Resistance in Escherichia coli Under Low and High Mutation Supply Rates

Antibiotic resistance is a major concern in public health worldwide, thus there is much interest in characterizing the mutational pathways through which susceptible bacteria evolve resistance. Here we use experimental evolution to explore the mutational pathways toward aminoglycoside resistance, using gentamicin as a model, under low and high mutation supply rates. Our results show that both normo and hypermutable strains of Escherichia coli are able to develop resistance to drug dosages > 1,000-fold higher than the minimal inhibitory concentration for their ancestors. Interestingly, such level of resistance was often associated with changes in susceptibility to other antibiotics, most prominently with increased resistance to fosfomycin. Whole-genome sequencing revealed that all resistant derivatives presented diverse mutations in five common genetic elements: fhuA, fusA and the atpIBEFHAGDC, cyoABCDE, and potABCD operons. Despite the large number of mutations acquired, hypermutable strains did not pay, apparently, fitness cost. In contrast to recent studies, we found that the mutation supply rate mainly affected the speed (tempo) but not the pattern (mode) of evolution: both backgrounds acquired the mutations in the same order, although the hypermutator strain did it faster. This observation is compatible with the adaptive landscape for high-level gentamicin resistance being relatively smooth, with few local maxima; which might be a common feature among antibiotics for which resistance involves multiple loci.

Divergent evolution peaks under intermediate population bottlenecks during bacterial experimental evolution

There is growing evidence that parallel molecular evolution is common, but its causes remain poorly understood. Demographic parameters such as population bottlenecks are predicted to be major determinants of parallelism. Here, we test the hypothesis that bottleneck intensity shapes parallel evolution by elucidating the genomic basis of adaptation to antibiotic-supplemented media in hundreds of populations of the bacterium Pseudomonas fluorescens Pf0-1. As expected, bottlenecking decreased the rate of phenotypic and molecular adaptation. Surprisingly, bottlenecking had no impact on the likelihood of parallel adaptive molecular evolution at a genome-wide scale. However, bottlenecking had a profound impact on the genes involved in antibiotic resistance. Specifically, under either intense or weak bottlenecking, resistance predominantly evolved by strongly beneficial mutations which provide high levels of antibiotic resistance. In contrast with intermediate bottlenecking regimes, resistance evolved by a greater diversity of genetic mechanisms, significantly reducing the observed levels of parallel genetic evolution. Our results demonstrate that population bottlenecking can be a major predictor of parallel evolution, but precisely how may be more complex than many simple theoretical predictions.

Sex in a test tube: testing the benefits of in vitro recombination

The origin and evolution of sex, and the associated role of recombination, present a major problem in biology. Sex typically involves recombination of closely related DNA or RNA sequences, which is fundamentally a random process that creates but also breaks up beneficial allele combinations. Directed evolution experiments, which combine in vitro mutation and recombination protocols with in vitro or in vivo selection, have proved to be an effective approach for improving functionality of nucleic acids and enzymes. As this approach allows extreme control over evolutionary conditions and parameters, it also facilitates the detection of small or position-specific recombination benefits and benefits associated with recombination between highly divergent genotypes. Yet, in vitro approaches have been largely exploratory and motivated by obtaining improved end products rather than testing hypotheses of recombination benefits. Here, we review the various experimental systems and approaches used by in vitro studies of recombination, discuss what they say about the evolutionary role of recombination, and sketch their potential for addressing extant questions about the evolutionary role of sex and recombination, in particular on complex fitness landscapes. We also review recent insights into the role of 'extracellular recombination' during the origin of life.This article is part of the themed issue 'Weird sex: the underappreciated diversity of sexual reproduction'.

Adaptive benefits from small mutation supplies in an antibiotic resistance enzyme

Populations with large mutation supplies adapt via the "greedy" substitution of the fittest genotype available, leading to fast and repeatable short-term responses. At longer time scales, smaller mutation supplies may in theory lead to larger improvements when distant high-fitness genotypes more readily evolve from lower-fitness intermediates. Here we test for long-term adaptive benefits from small mutation supplies using in vitro evolution of an antibiotic-degrading enzyme in the presence of a novel antibiotic. Consistent with predictions, large mutant libraries cause rapid initial adaptation via the substitution of cohorts of mutations, but show later deceleration and convergence. Smaller libraries show on average smaller initial, but also more variable, improvements, with two lines yielding alleles with exceptionally high resistance levels. These two alleles share three mutations with the large-library alleles, which are known from previous work, but also have unique mutations. Replay evolution experiments and analyses of the adaptive landscape of the enzyme suggest that the benefit resulted from a combination of avoiding mutational cohorts leading to local peaks and chance. Our results demonstrate adaptive benefits from limited mutation supplies on a rugged fitness landscape, which has implications for artificial selection protocols in biotechnology and argues for a better understanding of mutation supplies in clinical settings.

Genetic bottlenecks in intraspecies virus transmission

Ultimately, viral evolution is a consequence of mutations that arise within and spread between infected hosts. The transmission bottleneck determines how much of the viral diversity generated in one host passes to another during transmission. It therefore plays a vital role in linking within-host processes to larger evolutionary trends. Although many studies suggest that transmission severely restricts the amount of genetic diversity that passes between individuals, there are important exceptions to this rule. In many cases, the factors that determine the size of the transmission bottleneck are only beginning to be understood. Here, we review how transmission bottlenecks are measured, how they arise, and their consequences for viral evolution.

Effect of Population Size and Mutation Rate on the Evolution of RNA Sequences on an Adaptive Landscape Determined by RNA Folding

The dynamics of populations evolving on an adaptive landscape depends on multiple factors, including the structure of the landscape, the rate of mutations, and effective population size. Existing theoretical work often makes ad hoc and simplifying assumptions about landscape structure, whereas experimental work can vary important parameters only to a limited extent. We here overcome some of these limitations by simulating the adaptive evolution of RNA molecules, whose fitness is determined by the thermodynamics of RNA secondary structure folding. We study the influence of mutation rates and population sizes on final mean population fitness, on the substitution rates of mutations, and on population diversity. We show that evolutionary dynamics cannot be understood as a function of mutation rate µ, population size N, or population mutation rate Nµ alone. For example, at a given mutation rate, clonal interference prevents the fixation of beneficial mutations as population size increases, but larger populations still arrive at a higher mean fitness. In addition, at the highest population mutation rates we study, mean final fitness increases with population size, because small populations are driven to low fitness by the relatively higher incidence of mutations they experience. Our observations show that mutation rate and population size can interact in complex ways to influence the adaptive dynamics of a population on a biophysically motivated fitness landscape.

Mutation supply and the repeatability of selection for antibiotic resistance

Whether evolution can be predicted is a key question in evolutionary biology. Here we set out to better understand the repeatability of evolution, which is a necessary condition for predictability. We explored experimentally the effect of mutation supply and the strength of selective pressure on the repeatability of selection from standing genetic variation. Different sizes of mutant libraries of antibiotic resistance gene TEM-1 β-lactamase in Escherichia coli, generated by error-prone PCR, were subjected to different antibiotic concentrations. We determined whether populations went extinct or survived, and sequenced the TEM gene of the surviving populations. The distribution of mutations per allele in our mutant libraries followed a Poisson distribution. Extinction patterns could be explained by a simple stochastic model that assumed the sampling of beneficial mutations was key for survival. In most surviving populations, alleles containing at least one known large-effect beneficial mutation were present. These genotype data also support a model which only invokes sampling effects to describe the occurrence of alleles containing large-effect driver mutations. Hence, evolution is largely predictable given cursory knowledge of mutational fitness effects, the mutation rate and population size. There were no clear trends in the repeatability of selected mutants when we considered all mutations present. However, when only known large-effect mutations were considered, the outcome of selection is less repeatable for large libraries, in contrast to expectations. We show experimentally that alleles carrying multiple mutations selected from large libraries confer higher resistance levels relative to alleles with only a known large-effect mutation, suggesting that the scarcity of high-resistance alleles carrying multiple mutations may contribute to the decrease in repeatability at large library sizes.

Stronger selection can slow down evolution driven by recombination on a smooth fitness landscape

Stronger selection implies faster evolution—that is, the greater the force, the faster the change. This apparently self-evident proposition, however, is derived under the assumption that genetic variation within a population is primarily supplied by mutation (i.e. mutation-driven evolution). Here, we show that this proposition does not actually hold for recombination-driven evolution, i.e. evolution in which genetic variation is primarily created by recombination rather than mutation. By numerically investigating population genetics models of recombination, migration and selection, we demonstrate that stronger selection can slow down evolution on a perfectly smooth fitness landscape. Through simple analytical calculation, this apparently counter-intuitive result is shown to stem from two opposing effects of natural selection on the rate of evolution. On the one hand, natural selection tends to increase the rate of evolution by increasing the fixation probability of fitter genotypes. On the other hand, natural selection tends to decrease the rate of evolution by decreasing the chance of recombination between immigrants and resident individuals. As a consequence of these opposing effects, there is a finite selection pressure maximizing the rate of evolution. Hence, stronger selection can imply slower evolution if genetic variation is primarily supplied by recombination.

Genotypic Complexity of Fisher's Geometric Model

Fisher's geometric model was originally introduced to argue that complex adaptations must occur in small steps because of pleiotropic constraints. When supplemented with the assumption of additivity of mutational effects on phenotypic traits, it provides a simple mechanism for the emergence of genotypic epistasis from the nonlinear mapping of phenotypes to fitness. Of particular interest is the occurrence of reciprocal sign epistasis, which is a necessary condition for multipeaked genotypic fitness landscapes. Here we compute the probability that a pair of randomly chosen mutations interacts sign epistatically, which is found to decrease with increasing phenotypic dimension n, and varies nonmonotonically with the distance from the phenotypic optimum. We then derive expressions for the mean number of fitness maxima in genotypic landscapes comprised of all combinations of L random mutations. This number increases exponentially with L, and the corresponding growth rate is used as a measure of the complexity of the landscape. The dependence of the complexity on the model parameters is found to be surprisingly rich, and three distinct phases characterized by different landscape structures are identified. Our analysis shows that the phenotypic dimension, which is often referred to as phenotypic complexity, does not generally correlate with the complexity of fitness landscapes and that even organisms with a single phenotypic trait can have complex landscapes. Our results further inform the interpretation of experiments where the parameters of Fisher's model have been inferred from data, and help to elucidate which features of empirical fitness landscapes can be described by this model.

Different Evolutionary Paths to Complexity for Small and Large Populations of Digital Organisms

A major aim of evolutionary biology is to explain the respective roles of adaptive versus non-adaptive changes in the evolution of complexity. While selection is certainly responsible for the spread and maintenance of complex phenotypes, this does not automatically imply that strong selection enhances the chance for the emergence of novel traits, that is, the origination of complexity. Population size is one parameter that alters the relative importance of adaptive and non-adaptive processes: as population size decreases, selection weakens and genetic drift grows in importance. Because of this relationship, many theories invoke a role for population size in the evolution of complexity. Such theories are difficult to test empirically because of the time required for the evolution of complexity in biological populations. Here, we used digital experimental evolution to test whether large or small asexual populations tend to evolve greater complexity. We find that both small and large-but not intermediate-sized-populations are favored to evolve larger genomes, which provides the opportunity for subsequent increases in phenotypic complexity. However, small and large populations followed different evolutionary paths towards these novel traits. Small populations evolved larger genomes by fixing slightly deleterious insertions, while large populations fixed rare beneficial insertions that increased genome size. These results demonstrate that genetic drift can lead to the evolution of complexity in small populations and that purifying selection is not powerful enough to prevent the evolution of complexity in large populations.

Fisher's geometrical model and the mutational patterns of antibiotic resistance across dose gradients

Fisher's geometrical model (FGM) has been widely used to depict the fitness effects of mutations. It is a general model with few underlying assumptions that gives a large and comprehensive view of adaptive processes. It is thus attractive in several situations, for example adaptation to antibiotics, but comes with limitations, so that more mechanistic approaches are often preferred to interpret experimental data. It might be possible however to extend FGM assumptions to better account for mutational data. This is theoretically challenging in the context of antibiotic resistance because resistance mutations are assumed to be rare. In this article, we show with Escherichia coli how the fitness effects of resistance mutations screened at different doses of nalidixic acid vary across a dose-gradient. We found experimental patterns qualitatively consistent with the basic FGM (rate of resistance across doses, gamma distributed costs) but also unexpected patterns such as a decreasing mean cost of resistance with increasing screen dose. We show how different extensions involving mutational modules and variations in trait covariance across environments, can be discriminated based on these data. Overall, simple extensions of the FGM accounted well for complex mutational effects of resistance mutations across antibiotic doses.

Competition along trajectories governs adaptation rates towards antimicrobial resistance

The increasing availability of genotype-phenotype maps for different combinations of mutations has empowered evolutionary biologists with the tools to interrogate the predictability of adaptive evolution, especially in the context of the evolution of antimicrobial resistance. Large microbial populations are known to generate competing beneficial mutations, but determining how these mutations contribute to the adaptive trajectories that are most likely to be followed remains a challenge. Despite a recognition that there may also be competition between successive alleles on the same trajectory, prior studies have not fully considered how this impacts adaptation rates along, or likelihood of following, individual trajectories. Here, we develop a metric that quantifies the competition between successive alleles along adaptive trajectories and show how this competition largely governs the rate of evolution in simulations on empirical fitness landscapes for proteins involved in drug resistance in two species of malaria (Plasmodium falciparum and P. vivax). Our findings reveal that a trajectory with a larger-than-average initial fitness increase may have smaller fitness increases in later steps, which slows adaptation. In some circumstances, these trajectories may be outcompeted by alleles on faster alternative trajectories that are being explored simultaneously. The ability to predict adaptation rates along accessible trajectories has implications for efforts to manage antimicrobial resistance in real-world settings and for the broader intellectual pursuit of predictive evolution in complex adaptive fitness landscapes for a variety of application domains.

Modeling Tumor Clonal Evolution for Drug Combinations Design

Cancer is a clonal evolutionary process. This presents challenges for effective therapeutic intervention, given the constant selective pressure towards drug resistance. Mathematical modeling from population genetics, evolutionary dynamics, and engineering perspectives are being increasingly employed to study tumor progression, intratumoral heterogeneity, drug resistance, and rational drug scheduling and combinations design. In this review, we discuss promising opportunities these inter-disciplinary approaches hold for advances in cancer biology and treatment. We propose that quantitative modeling perspectives can complement emerging experimental technologies to facilitate enhanced understanding of disease progression and improved capabilities for therapeutic drug regimen designs.

Experimental evolution in biofilm populations

Biofilms are a major form of microbial life in which cells form dense surface associated communities that can persist for many generations. The long-life of biofilm communities means that they can be strongly shaped by evolutionary processes. Here, we review the experimental study of evolution in biofilm communities. We first provide an overview of the different experimental models used to study biofilm evolution and their associated advantages and disadvantages. We then illustrate the vast amount of diversification observed during biofilm evolution, and we discuss (i) potential ecological and evolutionary processes behind the observed diversification, (ii) recent insights into the genetics of adaptive diversification, (iii) the striking degree of parallelism between evolution experiments and real-life biofilms and (iv) potential consequences of diversification. In the second part, we discuss the insights provided by evolution experiments in how biofilm growth and structure can promote cooperative phenotypes. Overall, our analysis points to an important role of biofilm diversification and cooperation in bacterial survival and productivity. Deeper understanding of both processes is of key importance to design improved antimicrobial strategies and diagnostic techniques.

This review paper provides an overview of (i) the different experimental models used to study biofilm evolution, (ii) the vast amount of diversification observed during biofilm evolution (including potential causes and consequences) and (iii) recent insights in how growth in biofilms can lead to the evolution of cooperative phenotypes.

Repeatability of adaptation in experimental populations of different sizes

The degree to which evolutionary trajectories and outcomes are repeatable across independent populations depends on the relative contribution of selection, chance and history. Population size has been shown theoretically and empirically to affect the amount of variation that arises among independent populations adapting to the same environment. Here, we measure the contribution of selection, chance and history in different-sized experimental populations of the unicellular alga Chlamydomonas reinhardtii adapting to a high salt environment to determine which component of evolution is affected by population size. We find that adaptation to salt is repeatable at the fitness level in medium (Ne = 5 × 10(4)) and large (Ne = 4 × 10(5)) populations because of the large contribution of selection. Adaptation is not repeatable in small (Ne = 5 × 10(3)) populations because of large constraints from history. The threshold between stochastic and deterministic evolution in this case is therefore between effective population sizes of 10(3) and 10(4). Our results indicate that diversity across populations is more likely to be maintained if they are small. Experimental outcomes in large populations are likely to be robust and can inform our predictions about outcomes in similar situations.

Negative Epistasis and Evolvability in TEM-1 β-Lactamase--The Thin Line between an Enzyme's Conformational Freedom and Disorder

Epistasis is a key factor in evolution since it determines which combinations of mutations provide adaptive solutions and which mutational pathways toward these solutions are accessible by natural selection. There is growing evidence for the pervasiveness of sign epistasis--a complete reversion of mutational effects, particularly in protein evolution--yet its molecular basis remains poorly understood. We describe the structural basis of sign epistasis between G238S and R164S, two adaptive mutations in TEM-1 β-lactamase--an enzyme that endows antibiotics resistance. Separated by 10 Å, these mutations initiate two separate trajectories toward increased hydrolysis rates and resistance toward second and third-generation cephalosporins antibiotics. Both mutations allow the enzyme's active site to adopt alternative conformations and accommodate the new antibiotics. By solving the corresponding set of crystal structures, we found that R164S causes local disorder whereas G238S induces discrete conformations. When combined, the mutations in 238 and 164 induce local disorder whereby nonproductive conformations that perturb the enzyme's catalytic preorganization dominate. Specifically, Asn170 that coordinates the deacylating water molecule is misaligned, in both the free form and the inhibitor-bound double mutant. This local disorder is not restored by stabilizing global suppressor mutations and thus leads to an evolutionary cul-de-sac. Conformational dynamism therefore underlines the reshaping potential of protein's structures and functions but also limits protein evolvability because of the fragility of the interactions networks that maintain protein structures.

Genetic code evolution reveals the neutral emergence of mutational robustness, and information as an evolutionary constraint

The standard genetic code (SGC) is central to molecular biology and its origin and evolution is a fundamental problem in evolutionary biology, the elucidation of which promises to reveal much about the origins of life. In addition, we propose that study of its origin can also reveal some fundamental and generalizable insights into mechanisms of molecular evolution, utilizing concepts from complexity theory. The first is that beneficial traits may arise by non-adaptive processes, via a process of "neutral emergence". The structure of the SGC is optimized for the property of error minimization, which reduces the deleterious impact of point mutations. Via simulation, it can be shown that genetic codes with error minimization superior to the SGC can emerge in a neutral fashion simply by a process of genetic code expansion via tRNA and aminoacyl-tRNA synthetase duplication, whereby similar amino acids are added to codons related to that of the parent amino acid. This process of neutral emergence has implications beyond that of the genetic code, as it suggests that not all beneficial traits have arisen by the direct action of natural selection; we term these "pseudaptations", and discuss a range of potential examples. Secondly, consideration of genetic code deviations (codon reassignments) reveals that these are mostly associated with a reduction in proteome size. This code malleability implies the existence of a proteomic constraint on the genetic code, proportional to the size of the proteome (P), and that its reduction in size leads to an "unfreezing" of the codon - amino acid mapping that defines the genetic code, consistent with Crick's Frozen Accident theory. The concept of a proteomic constraint may be extended to propose a general informational constraint on genetic fidelity, which may be used to explain variously, differences in mutation rates in genomes with differing proteome sizes, differences in DNA repair capacity and genome GC content between organisms, a selective pressure in the evolution of sexual reproduction, and differences in translational fidelity. Lastly, the utility of the concept of an informational constraint to other diverse fields of research is explored.

The competition between simple and complex evolutionary trajectories in asexual populations

BACKGROUND

On rugged fitness landscapes where sign epistasis is common, adaptation can often involve either individually beneficial "uphill" mutations or more complex mutational trajectories involving fitness valleys or plateaus. The dynamics of the evolutionary process determine the probability that evolution will take any specific path among a variety of competing possible trajectories. Understanding this evolutionary choice is essential if we are to understand the outcomes and predictability of adaptation on rugged landscapes.

RESULTS

We present a simple model to analyze the probability that evolution will eschew immediately uphill paths in favor of crossing fitness valleys or plateaus that lead to higher fitness but less accessible genotypes. We calculate how this probability depends on the population size, mutation rates, and relevant selection pressures, and compare our analytical results to Wright-Fisher simulations.

CONCLUSION

We find that the probability of valley crossing depends nonmonotonically on population size: intermediate size populations are most likely to follow a "greedy" strategy of acquiring immediately beneficial mutations even if they lead to evolutionary dead ends, while larger and smaller populations are more likely to cross fitness valleys to reach distant advantageous genotypes. We explicitly identify the boundaries between these different regimes in terms of the relevant evolutionary parameters. Above a certain threshold population size, we show that the probability that the population finds the more distant peak depends only on a single simple combination of the relevant parameters.

Useful insights from evolutionary biology for developing perennial grain crops

Annual grain crops dominate agricultural landscapes and provide the majority of calories consumed by humanity. Perennial grain crops could potentially ameliorate the land degradation and off-site impacts associated with annual grain cropping. However, herbaceous perennial plants with constitutively high allocation to harvestable seeds are rare to absent in nature. Recent trade-off theory models suggest that rugged fitness landscapes may explain the absence of this form better than sink competition models. Artificial selection for both grain production and multiyear lifespan can lead to more rapid progress in the face of fitness and genetic trade-offs than natural selection but is likely to result in plant types that differ substantially from all current domestic crops. Perennial grain domestication is also likely to require the development of selection strategies that differ from published crop breeding methods, despite their success in improving long-domesticated crops; for this purpose, we have reviewed literature in the areas of population and evolutionary genetics, domestication, and molecular biology. Rapid domestication will likely require genes with large effect that are expected to exhibit strong pleiotropy and epistasis. Cryptic genetic variation will need to be deliberately exposed both to purge mildly deleterious alleles and to generate novel agronomic phenotypes. We predict that perennial grain domestication programs will benefit from population subdivision followed by selection for simple traits in each subpopulation, the evaluation of very large populations, high selection intensity, rapid cycling through generations, and heterosis. The latter may be particularly beneficial in the development of varieties with stable yield and tolerance to crowding.

Length of adaptive walk on uncorrelated and correlated fitness landscapes

We consider the adaptation dynamics of an asexual population that walks uphill on a rugged fitness landscape which is endowed with a large number of local fitness peaks. We work in a parameter regime where only those mutants that are a single mutation away are accessible, as a result of which the population eventually gets trapped at a local fitness maximum and the adaptive walk terminates. We study how the number of adaptive steps taken by the population before reaching a local fitness peak depends on the initial fitness of the population, the extreme value distribution of the beneficial mutations, and correlations among the fitnesses. Assuming that the relative fitness difference between successive steps is small, we analytically calculate the average walk length for both uncorrelated and correlated fitnesses in all extreme value domains for a given initial fitness. We present numerical results for the model where the fitness differences can be large and find that the walk length behavior differs from that in the former model in the Fréchet domain of extreme value theory. We also discuss the relevance of our results to microbial experiments.

Quantifying the role of population subdivision in evolution on rugged fitness landscapes

Natural selection drives populations towards higher fitness, but crossing fitness valleys or plateaus may facilitate progress up a rugged fitness landscape involving epistasis. We investigate quantitatively the effect of subdividing an asexual population on the time it takes to cross a fitness valley or plateau. We focus on a generic and minimal model that includes only population subdivision into equivalent demes connected by global migration, and does not require significant size changes of the demes, environmental heterogeneity or specific geographic structure. We determine the optimal speedup of valley or plateau crossing that can be gained by subdivision, if the process is driven by the deme that crosses fastest. We show that isolated demes have to be in the sequential fixation regime for subdivision to significantly accelerate crossing. Using Markov chain theory, we obtain analytical expressions for the conditions under which optimal speedup is achieved: valley or plateau crossing by the subdivided population is then as fast as that of its fastest deme. We verify our analytical predictions through stochastic simulations. We demonstrate that subdivision can substantially accelerate the crossing of fitness valleys and plateaus in a wide range of parameters extending beyond the optimal window. We study the effect of varying the degree of subdivision of a population, and investigate the trade-off between the magnitude of the optimal speedup and the width of the parameter range over which it occurs. Our results, obtained for fitness valleys and plateaus, also hold for weakly beneficial intermediate mutations. Finally, we extend our work to the case of a population connected by migration to one or several smaller islands. Our results demonstrate that subdivision with migration alone can significantly accelerate the crossing of fitness valleys and plateaus, and shed light onto the quantitative conditions necessary for this to occur.

Experimental evidence has recently been accumulating to suggest that fitness landscape ruggedness is common in a variety of organisms. Rugged landscapes arise from interactions between genetic variants, called epistasis, which can lead to fitness valleys or plateaus. The time needed to cross such fitness valleys or plateaus exhibits a rich dependence on population size, since stochastic effects have higher importance in small populations, increasing the probability of fixation of neutral or deleterious mutants. This may lead to an advantage of population subdivision, a possibility which has been strongly debated for nearly one hundred years. In this work, we quantitatively determine when, and to what extent, population subdivision accelerates valley and plateau crossing. Using the simple model of an asexual population subdivided into identical demes connected by gobal migration, we derive the conditions under which crossing by a subdivided population is driven by its fastest deme, thus giving rise to the maximal speedup. Our analytical predictions are verified using stochastic simulations. We investigate the effect of varying the degree of subdivision of a population. We generalize our results to weakly beneficial intermediates and to different population structures. We discuss the magnitude and robustness of the effect for realistic parameter values.

Variation in the fitness effects of mutations with population density and size in Escherichia coli

The fitness effects of mutations are context specific and depend on both external (e.g., environment) and internal (e.g., cellular stress, genetic background) factors. The influence of population size and density on fitness effects are unknown, despite the central role population size plays in the supply and fixation of mutations. We addressed this issue by comparing the fitness of 92 Keio strains (Escherichia coli K12 single gene knockouts) at comparatively high (1.2×10(7) CFUs/mL) and low (2.5×10(2) CFUs/mL) densities, which also differed in population size (high: 1.2×10(8); low: 1.25×10(3)). Twenty-eight gene deletions (30%) exhibited a fitness difference, ranging from 5 to 174% (median: 35%), between the high and low densities. Our analyses suggest this variation among gene deletions in fitness responses reflected in part both gene orientation and function, of the gene properties we examined (genomic position, length, orientation, and function). Although we could not determine the relative effects of population density and size, our results suggest fitness effects of mutations vary with these two factors, and this variation is gene-specific. Besides being a mechanism for density-dependent selection (r-K selection), the dependence of fitness effects on population density and size has implications for any population that varies in size over time, including populations undergoing evolutionary rescue, species invasions into novel habitats, and cancer progression and metastasis. Further, combined with recent advances in understanding the roles of other context-specific factors in the fitness effects of mutations, our results will help address theoretical and applied biological questions more realistically.

Empirical fitness landscapes and the predictability of evolution

The genotype-fitness map (that is, the fitness landscape) is a key determinant of evolution, yet it has mostly been used as a superficial metaphor because we know little about its structure. This is now changing, as real fitness landscapes are being analysed by constructing genotypes with all possible combinations of small sets of mutations observed in phylogenies or in evolution experiments. In turn, these first glimpses of empirical fitness landscapes inspire theoretical analyses of the predictability of evolution. Here, we review these recent empirical and theoretical developments, identify methodological issues and organizing principles, and discuss possibilities to develop more realistic fitness landscape models.

Non-additive costs and interactions alter the competitive dynamics of co-occurring ecologically distinct plasmids

Plasmids play an important role in shaping bacterial evolution and adaptation to heterogeneous environments. As modular genetic elements that are often conjugative, the selective pressures that act on plasmid-borne genes are distinct from those that act on the chromosome. Many bacteria are co-infected by multiple plasmids that impart niche-specific phenotypes. Thus, in addition to host-plasmid dynamics, interactions between co-infecting plasmids are likely to be important drivers of plasmid population dynamics, evolution and ecology. Agrobacterium tumefaciens is a facultative plant pathogen that commonly harbours two distinct megaplasmids. Virulence depends on the presence of the tumour-inducing (Ti) plasmid, with benefits that are primarily restricted to the disease environment. Here, we demonstrate that a second megaplasmid, the At plasmid, confers a competitive advantage in the rhizosphere. To assess the individual and interactive costs of these plasmids, we generated four isogenic derivatives: plasmidless, pAt only, pTi only and pAtpTi, and performed pairwise competitions under carbon-limiting conditions. These studies reveal a low cost to the virulence plasmid when outside of the disease environment, and a strikingly high cost to the At plasmid. In addition, the costs of pAt and pTi in the same host were significantly lower than predicted based on single plasmid costs, signifying the first demonstration of non-additivity between naturally occurring co-resident plasmids. Based on these empirically demonstrated costs and benefits, we developed a resource-consumer model to generate predictions about the frequencies of these genotypes in relevant environments, showing that non-additivity between co-residing plasmids allows for their stable coexistence across environments.

The dynamics of genetic draft in rapidly adapting populations

The accumulation of beneficial mutations on competing genetic backgrounds in rapidly adapting populations has a striking impact on evolutionary dynamics. This effect, known as clonal interference, causes erratic fluctuations in the frequencies of observed mutations, randomizes the fixation times of successful mutations, and leaves distinct signatures on patterns of genetic variation. Here, we show how this form of "genetic draft" affects the forward-time dynamics of site frequencies in rapidly adapting asexual populations. We calculate the probability that mutations at individual sites shift in frequency over a characteristic timescale, extending Gillespie's original model of draft to the case where many strongly selected beneficial mutations segregate simultaneously. We then derive the sojourn time of mutant alleles, the expected fixation time of successful mutants, and the site frequency spectrum of beneficial and neutral mutations. Finally, we show how this form of draft affects inferences in the McDonald-Kreitman test and how it relates to recent observations that some aspects of genetic diversity are described by the Bolthausen-Sznitman coalescent in the limit of very rapid adaptation.