Key issues review: evolution on rugged adaptive landscapes

Have you tried turning it off and on again? Oscillating selection to enhance fitness-landscape traversal in adaptive laboratory evolution experiments

Adaptive Laboratory Evolution (ALE) is a powerful tool for engineering and understanding microbial physiology. ALE relies on the selection and enrichment of mutations that enable survival or faster growth under a selective condition imposed by the experimental setup. Phenotypic fitness landscapes are often underpinned by complex genotypes involving multiple genes, with combinatorial positive and negative effects on fitness. Such genotype relationships result in mutational fitness landscapes with multiple local fitness maxima and valleys. Traversing local maxima to find a global maximum often requires an individual or sub-population of cells to traverse fitness valleys. Traversing involves gaining mutations that are not adaptive for a given local maximum but are necessary to 'peak shift' to another local maximum, or eventually a global maximum. Despite these relatively well understood evolutionary principles, and the combinatorial genotypes that underlie most metabolic phenotypes, the majority of applied ALE experiments are conducted using constant selection pressures. The use of constant pressure can result in populations becoming trapped within local maxima, and often precludes the attainment of optimum phenotypes associated with global maxima. Here, we argue that oscillating selection pressures is an easily accessible mechanism for traversing fitness landscapes in ALE experiments, and provide theoretical and practical frameworks for implementation.

The metabolic, virulence and antimicrobial resistance profiles of colonising Streptococcus pneumoniae shift after PCV13 introduction in urban Malawi

Streptococcus pneumoniae causes substantial mortality among children under 5-years-old worldwide. Polysaccharide conjugate vaccines (PCVs) are highly effective at reducing vaccine serotype disease, but emergence of non-vaccine serotypes and persistent nasopharyngeal carriage threaten this success. We investigated the hypothesis that following vaccine, adapted pneumococcal genotypes emerge with the potential for vaccine escape. We genome sequenced 2804 penumococcal isolates, collected 4-8 years after introduction of PCV13 in Blantyre, Malawi. We developed a pipeline to cluster the pneumococcal population based on metabolic core genes into "Metabolic genotypes" (MTs). We show that S. pneumoniae population genetics are characterised by emergence of MTs with distinct virulence and antimicrobial resistance (AMR) profiles. Preliminary in vitro and murine experiments revealed that representative isolates from emerging MTs differed in growth, haemolytic, epithelial infection, and murine colonisation characteristics. Our results suggest that in the context of PCV13 introduction, pneumococcal population dynamics had shifted, a phenomenon that could further undermine vaccine control and promote spread of AMR.

Extinction scenarios in evolutionary processes: a multinomial Wright-Fisher approach

We study a discrete-time multi-type Wright-Fisher population process. The mean-field dynamics of the stochastic process is induced by a general replicator difference equation. We prove several results regarding the asymptotic behavior of the model, focusing on the impact of the mean-field dynamics on it. One of the results is a limit theorem that describes sufficient conditions for an almost certain path to extinction, first eliminating the type which is the least fit at the mean-field equilibrium. The effect is explained by the metastability of the stochastic system, which under the conditions of the theorem spends almost all time before the extinction event in a neighborhood of the equilibrium. In addition to the limit theorems, we propose a maximization principle for a general deterministic replicator dynamics and study its implications for the stochastic model.

Phenotype switching of the mutation rate facilitates adaptive evolution

The mutation rate plays an important role in adaptive evolution. It can be modified by mutator and anti-mutator alleles. Recent empirical evidence hints that the mutation rate may vary among genetically identical individuals: evidence from bacteria suggests that the mutation rate can be affected by expression noise of a DNA repair protein and potentially also by translation errors in various proteins. Importantly, this non-genetic variation may be heritable via a transgenerational epigenetic mode of inheritance, giving rise to a mutator phenotype that is independent from mutator alleles. Here, we investigate mathematically how the rate of adaptive evolution is affected by the rate of mutation rate phenotype switching. We model an asexual population with two mutation rate phenotypes, non-mutator and mutator. An offspring may switch from its parental phenotype to the other phenotype. We find that switching rates that correspond to so-far empirically described non-genetic systems of inheritance of the mutation rate lead to higher rates of adaptation on both artificial and natural fitness landscapes. These switching rates can maintain within the same individuals both a mutator phenotype and intermediary mutations, a combination that facilitates adaptation. Moreover, non-genetic inheritance increases the proportion of mutators in the population, which in turn increases the probability of hitchhiking of the mutator phenotype with adaptive mutations. This in turns facilitates the acquisition of additional adaptive mutations. Our results rationalize recently observed noise in the expression of proteins that affect the mutation rate and suggest that non-genetic inheritance of this phenotype may facilitate evolutionary adaptive processes.

Evolutionary modeling suggests that addictions may be driven by competition-induced microbiome dysbiosis

Recent studies revealed mechanisms by which the microbiome affects its host's brain, behavior and wellbeing, and that dysbiosis - persistent microbiome-imbalance - is associated with the onset and progress of various chronic diseases, including addictive behaviors. Yet, understanding of the ecological and evolutionary processes that shape the host-microbiome ecosystem and affect the host state, is still limited. Here we propose that competition dynamics within the microbiome, associated with host-microbiome mutual regulation, may promote dysbiosis and aggravate addictive behaviors. We construct a mathematical framework, modeling the dynamics of the host-microbiome ecosystem in response to alterations. We find that when this ecosystem is exposed to substantial perturbations, the microbiome may shift towards a composition that reinforces the new host state. Such a positive feedback loop augments post-perturbation imbalances, hindering attempts to return to the initial equilibrium, promoting relapse episodes and prolonging addictions. We show that the initial microbiome composition is a key factor: a diverse microbiome enhances the ecosystem's resilience, whereas lower microbiome diversity is more prone to lead to dysbiosis, exacerbating addictions. This framework provides evolutionary and ecological perspectives on host-microbiome interactions and their implications for host behavior and health, while offering verifiable predictions with potential relevance to clinical treatments.

Antifungal alternation can be beneficial for durability but at the cost of generalist resistance

The evolution of resistance to pesticides is a major burden in agriculture. Resistance management involves maximizing selection pressure heterogeneity, particularly by combining active ingredients with different modes of action. We tested the hypothesis that alternation may delay the build-up of resistance not only by spreading selection pressure over longer periods, but also by decreasing the rate of evolution of resistance to alternated fungicides, by applying an experimental evolution approach to the economically important crop pathogen Zymoseptoria tritici. Our results show that alternation is either neutral or slows the overall resistance evolution rate, relative to continuous fungicide use, but results in higher levels of generalism in evolved lines. We demonstrate that the nature of the fungicides, and therefore their relative intrinsic risk of resistance may underly this trade-off, more so than the number of fungicides and the rhythm of alternation. This trade-off is also dynamic over the course of resistance evolution. These findings open up new possibilities for tailoring resistance management effectively while optimizing interplay between alternation components.

Intragenic compensation through the lens of deep mutational scanning

A significant fraction of mutations in proteins are deleterious and result in adverse consequences for protein function, stability, or interaction with other molecules. Intragenic compensation is a specific case of positive epistasis when a neutral missense mutation cancels effect of a deleterious mutation in the same protein. Permissive compensatory mutations facilitate protein evolution, since without them all sequences would be extremely conserved. Understanding compensatory mechanisms is an important scientific challenge at the intersection of protein biophysics and evolution. In human genetics, intragenic compensatory interactions are important since they may result in variable penetrance of pathogenic mutations or fixation of pathogenic human alleles in orthologous proteins from related species. The latter phenomenon complicates computational and clinical inference of an allele's pathogenicity. Deep mutational scanning is a relatively new technique that enables experimental studies of functional effects of thousands of mutations in proteins. We review the important aspects of the field and discuss existing limitations of current datasets. We reviewed ten published DMS datasets with quantified functional effects of single and double mutations and described rates and patterns of intragenic compensation in eight of them.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-022-01005-w.

Fitness landscape analysis of a tRNA gene reveals that the wild type allele is sub-optimal, yet mutationally robust

Fitness landscape mapping and the prediction of evolutionary trajectories on these landscapes are major tasks in evolutionary biology research. Evolutionary dynamics is tightly linked to the landscape topography, but this relation is not straightforward. Here, we analyze a fitness landscape of a yeast tRNA gene, previously measured under four different conditions. We find that the wild type allele is sub-optimal, and 8–10% of its variants are fitter. We rule out the possibilities that the wild type is fittest on average on these four conditions or located on a local fitness maximum. Notwithstanding, we cannot exclude the possibility that the wild type might be fittest in some of the many conditions in the complex ecology that yeast lives at. Instead, we find that the wild type is mutationally robust (“flat”), while more fit variants are typically mutationally fragile. Similar observations of mutational robustness or flatness have been so far made in very few cases, predominantly in viral genomes.

Systematic investigation of the link between enzyme catalysis and cold adaptation

Cold temperature is prevalent across the biosphere and slows the rates of chemical reactions. Increased catalysis has been predicted to be a dominant adaptive trait of enzymes to reduced temperature, and this expectation has informed physical models for enzyme catalysis and influenced bioprospecting strategies. To systematically test rate enhancement as an adaptive trait to cold, we paired kinetic constants of 2223 enzyme reactions with their organism's optimal growth temperature (T_Growth) and analyzed trends of rate constants as a function of T_Growth. These data do not support a general increase in rate enhancement in cold adaptation. In the model enzyme ketosteroid isomerase (KSI), there is prior evidence for temperature adaptation from a change in an active site residue that results in a tradeoff between activity and stability. Nevertheless, we found that little of the rate constant variation for 20 KSI variants was accounted for by T_Growth. In contrast, and consistent with prior expectations, we observed a correlation between stability and T_Growth across 433 proteins. These results suggest that temperature exerts a weaker selection pressure on enzyme rate constants than stability and that evolutionary forces other than temperature are responsible for the majority of enzymatic rate constant variation.

An Antigenic Thrift-Based Approach to Influenza Vaccine Design

The antigenic drift theory states that influenza evolves via the gradual accumulation of mutations, decreasing a host’s immune protection against previous strains. Influenza vaccines are designed accordingly, under the premise of antigenic drift. However, a paradox exists at the centre of influenza research. If influenza evolved primarily through mutation in multiple epitopes, multiple influenza strains should co-circulate. Such a multitude of strains would render influenza vaccines quickly inefficacious. Instead, a single or limited number of strains dominate circulation each influenza season. Unless additional constraints are placed on the evolution of influenza, antigenic drift does not adequately explain these observations. Here, we explore the constraints placed on antigenic drift and a competing theory of influenza evolution – antigenic thrift. In contrast to antigenic drift, antigenic thrift states that immune selection targets epitopes of limited variability, which constrain the variability of the virus. We explain the implications of antigenic drift and antigenic thrift and explore their current and potential uses in the context of influenza vaccine design.

Gamblers: An Antibiotic-Induced Evolvable Cell Subpopulation Differentiated by Reactive-Oxygen-Induced General Stress Response

Antibiotics can induce mutations that cause antibiotic resistance. Yet, despite their importance, mechanisms of antibiotic-promoted mutagenesis remain elusive. We report that the fluoroquinolone antibiotic ciprofloxacin (cipro) induces mutations by triggering transient differentiation of a mutant-generating cell subpopulation, using reactive oxygen species (ROS). Cipro-induced DNA breaks activate the Escherichia coli SOS DNA-damage response and error-prone DNA polymerases in all cells. However, mutagenesis is limited to a cell subpopulation in which electron transfer together with SOS induce ROS, which activate the sigma-S (σS) general-stress response, which allows mutagenic DNA-break repair. When sorted, this small σS-response-"on" subpopulation produces most antibiotic cross-resistant mutants. A U.S. Food and Drug Administration (FDA)-approved drug prevents σS induction, specifically inhibiting antibiotic-promoted mutagenesis. Further, SOS-inhibited cell division, which causes multi-chromosome cells, promotes mutagenesis. The data support a model in which within-cell chromosome cooperation together with development of a "gambler" cell subpopulation promote resistance evolution without risking most cells.

Computational Complexity as an Ultimate Constraint on Evolution

Experiments show that evolutionary fitness landscapes can have a rich combinatorial structure due to epistasis. For some landscapes, this structure can produce a computational constraint that prevents evolution from finding local fitness optima-thus overturning the traditional assumption that local fitness peaks can always be reached quickly if no other evolutionary forces challenge natural selection. Here, I introduce a distinction between easy landscapes of traditional theory where local fitness peaks can be found in a moderate number of steps, and hard landscapes where finding local optima requires an infeasible amount of time. Hard examples exist even among landscapes with no reciprocal sign epistasis; on these semismooth fitness landscapes, strong selection weak mutation dynamics cannot find the unique peak in polynomial time. More generally, on hard rugged fitness landscapes that include reciprocal sign epistasis, no evolutionary dynamics-even ones that do not follow adaptive paths-can find a local fitness optimum quickly. Moreover, on hard landscapes, the fitness advantage of nearby mutants cannot drop off exponentially fast but has to follow a power-law that long-term evolution experiments have associated with unbounded growth in fitness. Thus, the constraint of computational complexity enables open-ended evolution on finite landscapes. Knowing this constraint allows us to use the tools of theoretical computer science and combinatorial optimization to characterize the fitness landscapes that we expect to see in nature. I present candidates for hard landscapes at scales from single genes, to microbes, to complex organisms with costly learning (Baldwin effect) or maintained cooperation (Hankshaw effect). Just how ubiquitous hard landscapes (and the corresponding ultimate constraint on evolution) are in nature becomes an open empirical question.

Molecular Fitness Landscapes from High-Coverage Sequence Profiling

The function of fitness (or molecular activity) in the space of all possible sequences is known as the fitness landscape. Evolution is a random walk on the fitness landscape, with a bias toward climbing hills. Mapping the topography of real fitness landscapes is fundamental to understanding evolution, but previous efforts were hampered by the difficulty of obtaining large, quantitative data sets. The accessibility of high-throughput sequencing (HTS) has transformed this study, enabling large-scale enumeration of fitness for many mutants and even complete sequence spaces in some cases. We review the progress of high-throughput studies in mapping molecular fitness landscapes, both in vitro and in vivo, as well as opportunities for future research. Such studies are rapidly growing in number. HTS is expected to have a profound effect on the understanding of real molecular fitness landscapes.

With a little help from my friends: cooperation can accelerate the rate of adaptive valley crossing

Background

Natural selection favors changes that lead to genotypes possessing high fitness. A conflict arises when several mutations are required for adaptation, but each mutation is separately deleterious. The process of a population evolving from a genotype encoding for a local fitness maximum to a higher fitness genotype is termed an adaptive peak shift.

Results

Here we suggest cooperative behavior as a factor that can facilitate adaptive peak shifts. We model cooperation in a public goods scenario, wherein each individual contributes resources that are later equally redistributed among all cooperating individuals. We use mathematical modeling and stochastic simulations to study the effect of cooperation on peak shifts in both panmictic and structured populations. Our results show that cooperation can substantially affect the rate of complex adaptation. Furthermore, we show that cooperation increases the population diversity throughout the peak shift process, thus increasing the robustness of the population to sudden environmental changes.

Conclusions

We provide a new explanation to adaptive valley crossing in natural populations and suggest that the long term evolution of a species depends on its social behavior.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0983-2) contains supplementary material, which is available to authorized users.