The effect of bacterial recombination on adaptation on fitness landscapes with limited peak accessibility

Emergence and widespread circulation of a recombinant SARS-CoV-2 lineage in North America

Although recombination is a feature of coronavirus evolution, previously detected recombinant lineages of SARS-CoV-2 have shown limited circulation thus far. Here, we present a detailed phylogenetic analysis of four SARS-CoV-2 lineages to investigate the possibility of virus recombination among them. Our analyses reveal well-supported phylogenetic differences between the Orf1ab region encoding viral non-structural proteins and the rest of the genome, including Spike (S) protein and remaining reading frames. By accounting for several deletions in NSP6, Orf3a, and S, we conclude that the B.1.628 major cluster, now designated as lineage XB, originated from a recombination event between viruses of B.1.631 and B.1.634 lineages. This scenario is supported by the spatiotemporal distribution of these lineages across the USA and Mexico during 2021, suggesting that the recombination event originated in this geographical region. This event raises important questions regarding the role and potential effects of recombination on SARS-CoV-2 evolution.

Horizontal gene transfer and adaptive evolution in bacteria

Horizontal gene transfer (HGT) is arguably the most conspicuous feature of bacterial evolution. Evidence for HGT is found in most bacterial genomes. Although HGT can considerably alter bacterial genomes, not all transfer events may be biologically significant and may instead represent the outcome of an incessant evolutionary process that only occasionally has a beneficial purpose. When adaptive transfers occur, HGT and positive selection may result in specific, detectable signatures in genomes, such as gene-specific sweeps or increased transfer rates for genes that are ecologically relevant. In this Review, we first discuss the various mechanisms whereby HGT occurs, how the genetic signatures shape patterns of genomic variation and the distinct bioinformatic algorithms developed to detect these patterns. We then discuss the evolutionary theory behind HGT and positive selection in bacteria, and discuss the approaches developed over the past decade to detect transferred DNA that may be involved in adaptation to new environments.

Bacterial fitness landscapes stratify based on proteome allocation associated with discrete aero-types

The fitness landscape is a concept commonly used to describe evolution towards optimal phenotypes. It can be reduced to mechanistic detail using genome-scale models (GEMs) from systems biology. We use recently developed GEMs of Metabolism and protein Expression (ME-models) to study the distribution of Escherichia coli phenotypes on the rate-yield plane. We found that the measured phenotypes distribute non-uniformly to form a highly stratified fitness landscape. Systems analysis of the ME-model simulations suggest that this stratification results from discrete ATP generation strategies. Accordingly, we define "aero-types", a phenotypic trait that characterizes how a balanced proteome can achieve a given growth rate by modulating 1) the relative utilization of oxidative phosphorylation, glycolysis, and fermentation pathways; and 2) the differential employment of electron-transport-chain enzymes. This global, quantitative, and mechanistic systems biology interpretation of fitness landscape formed upon proteome allocation offers a fundamental understanding of bacterial physiology and evolution dynamics.

Recombination and mutational robustness in neutral fitness landscapes

Mutational robustness quantifies the effect of random mutations on fitness. When mutational robustness is high, most mutations do not change fitness or have only a minor effect on it. From the point of view of fitness landscapes, robust genotypes form neutral networks of almost equal fitness. Using deterministic population models it has been shown that selection favors genotypes inside such networks, which results in increased mutational robustness. Here we demonstrate that this effect is massively enhanced by recombination. Our results are based on a detailed analysis of mesa-shaped fitness landscapes, where we derive precise expressions for the dependence of the robustness on the landscape parameters for recombining and non-recombining populations. In addition, we carry out numerical simulations on different types of random holey landscapes as well as on an empirical fitness landscape. We show that the mutational robustness of a genotype generally correlates with its recombination weight, a new measure that quantifies the likelihood for the genotype to arise from recombination. We argue that the favorable effect of recombination on mutational robustness is a highly universal feature that may have played an important role in the emergence and maintenance of mechanisms of genetic exchange.

Ecology and evolution of facilitation among symbionts

Facilitation occurs when one species positively impacts the fitness of another, and has predominantly been studied in free-living species like plants. Facilitation can also occur among symbiont (mutualistic or parasitic) species or strains, but equivalent studies are scarce. To advance an integrated view of the effect of facilitation on symbiont ecology and evolution, we review empirical evidence and their underlying mechanisms, explore the factors favouring its emergence, and discuss its consequences for virulence and transmission. We argue that the facilitation concept can improve understanding of the evolutionary forces shaping symbiont communities and their effects on hosts.

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.

Recombination and peak jumping

We show that genetic recombination can be a powerful mechanism for escaping suboptimal peaks. Recent studies of empirical fitness landscapes reveal complex gene interactions and multiple peaks. However, classical work on recombination largely ignores the effect of complex gene interactions. Briefly, we restrict to fitness landscapes where the global peak is difficult to access. If the optimal genotype can be generated by shuffling genes present in the population, then recombination will produce the genotype. If, in addition, recombination is sufficiently rare, then the proportion of the genotype is expected to increase. Specifically, we consider landscapes where shuffling of suboptimal peak genotypes can produce the global peak genotype. The advantage of recombination we identify has no correspondence for 2-locus systems or for smooth landscapes. The effect of recombination indicated is sometimes extreme, also for rare recombination, in the sense that shutting off recombination could result in the organism failing to adapt. A standard question about recombination is whether the mechanism tends to accelerate or decelerate adaptation. However, we argue that extreme effects may be more important than how the majority falls. In a limited sense, our result can be considered a support for Sewall Wright’s view that adaptation sometimes works better in subdivided populations.

Weak Epistasis May Drive Adaptation in Recombining Bacteria

The impact of epistasis on the evolution of multi-locus traits depends on recombination. While sexually reproducing eukaryotes recombine so frequently that epistasis between polymorphisms is not considered to play a large role in short-term adaptation, many bacteria also recombine, some to the degree that their populations are described as "panmictic" or "freely recombining." However, whether this recombination is sufficient to limit the ability of selection to act on epistatic contributions to fitness is unknown. We quantify homologous recombination in five bacterial pathogens and use these parameter estimates in a multilocus model of bacterial evolution with additive and epistatic effects. We find that even for highly recombining species (e.g., Streptococcus pneumoniae or Helicobacter pylori), selection on weak interactions between distant mutations is nearly as efficient as for an asexual species, likely because homologous recombination typically transfers only short segments. However, for strong epistasis, bacterial recombination accelerates selection, with the dynamics dependent on the amount of recombination and the number of loci. Epistasis may thus play an important role in both the short- and long-term adaptive evolution of bacteria, and, unlike in eukaryotes, is not limited to strong effect sizes, closely linked loci, or other conditions that limit the impact of recombination.

Mutation and recombination in pathogen evolution: Relevance, methods and controversies

Mutation and recombination drive the evolution of most pathogens by generating the genetic variants upon which selection operates. Those variants can, for example, confer resistance to host immune systems and drug therapies or lead to epidemic outbreaks. Given their importance, diverse evolutionary studies have investigated the abundance and consequences of mutation and recombination in pathogen populations. However, some controversies persist regarding the contribution of each evolutionary force to the development of particular phenotypic observations (e.g., drug resistance). In this study, we revise the importance of mutation and recombination in the evolution of pathogens at both intra-host and inter-host levels. We also describe state-of-the-art analytical methodologies to detect and quantify these two evolutionary forces, including biases that are often ignored in evolutionary studies. Finally, we present some of our former studies involving pathogenic taxa where mutation and recombination played crucial roles in the recovery of pathogenic fitness, the generation of interspecific genetic diversity, or the design of centralized vaccines. This review also illustrates several common controversies and pitfalls in the analysis and in the evaluation and interpretation of mutation and recombination outcomes.

Evolution at 'Sutures' and 'Centers': Recombination Can Aid Adaptation of Spatially Structured Populations on Rugged Fitness Landscapes

Epistatic interactions among genes can give rise to rugged fitness landscapes, in which multiple “peaks” of high-fitness allele combinations are separated by “valleys” of low-fitness genotypes. How populations traverse rugged fitness landscapes is a long-standing question in evolutionary biology. Sexual reproduction may affect how a population moves within a rugged fitness landscape. Sex may generate new high-fitness genotypes by recombination, but it may also destroy high-fitness genotypes by shuffling the genes of a fit parent with a genetically distinct mate, creating low-fitness offspring. Either of these opposing aspects of sex require genotypic diversity in the population. Spatially structured populations may harbor more diversity than well-mixed populations, potentially amplifying both positive and negative effects of sex. On the other hand, spatial structure leads to clumping in which mating is more likely to occur between like types, diminishing the effects of recombination. In this study, we use computer simulations to investigate the combined effects of recombination and spatial structure on adaptation in rugged fitness landscapes. We find that spatially restricted mating and offspring dispersal may allow multiple genotypes inhabiting suboptimal peaks to coexist, and recombination at the “sutures” between the clusters of these genotypes can create genetically novel offspring. Sometimes such an offspring genotype inhabits a new peak on the fitness landscape. In such a case, spatially restricted mating allows this fledgling subpopulation to avoid recombination with distinct genotypes, as mates are more likely to be the same genotype. Such population “centers” can allow nascent peaks to establish despite recombination. Spatial structure may therefore allow an evolving population to enjoy the creative side of sexual recombination while avoiding its destructive side.

For a novel genotype to establish in a population, it must (1) be created, and (2) not be subsequently lost. Recombination is a double-edged sword in this process, potentially fostering creation, but also hastening loss as the novel genotype is being recombined with other genotypes, especially when rare. In this study, we find that spatial structure may affect both the creative and destructive aspects of recombination in rugged fitness landscapes. By slowing the spread of high-fitness genotypes, spatially restricted mating and dispersal may allow diverse subpopulations to arise. Reproduction across the borders of these subpopulations—at “sutures”—may create genetic novelty. Depending on the topography of the fitness landscape, such novelty may be in the domain of attraction of a new, higher peak; the population may “peak-jump” to an area of genotype space unlikely to be explored by mutation alone. Lineages founded by peak-jumping events are particularly prone to early extinction, as recombination with unlike genotypes may disrupt the rare allele combination and thereby produce low-fitness offspring. However, these fledgling peak lineages may be protected from early extinction by mating within small homotypic clusters—in “centers”. Thus, spatial structure may allow a population to create rare genotypes via recombination, and allow those rare genotypes to persist despite recombination.

Benefits of a Recombination-Proficient Escherichia coli System for Adaptive Laboratory Evolution

Adaptive laboratory evolution typically involves the propagation of organisms asexually to select for mutants with the desired phenotypes. However, asexual evolution is prone to competition among beneficial mutations (clonal interference) and the accumulation of hitchhiking and neutral mutations. The benefits of horizontal gene transfer toward overcoming these known disadvantages of asexual evolution were characterized in a strain of Escherichia coli engineered for superior sexual recombination (genderless). Specifically, we experimentally validated the capacity of the genderless strain to reduce the mutational load and recombine beneficial mutations. We also confirmed that inclusion of multiple origins of transfer influences both the frequency of genetic exchange throughout the chromosome and the linkage of donor DNA. We built a simple kinetic model to estimate recombination frequency as a function of transfer size and relative genotype enrichment in batch transfers; the model output correlated well with the experimental data. Our results provide strong support for the advantages of utilizing the genderless strain over its asexual counterpart during adaptive laboratory evolution for generating beneficial mutants with reduced mutational load.

IMPORTANCE

Over 80 years ago Fisher and Muller began a debate on the origins of sexual recombination. Although many aspects of sexual recombination have been examined at length, experimental evidence behind the behaviors of recombination in many systems and the means to harness it remain elusive. In this study, we sought to experimentally validate some advantages of recombination in typically asexual Escherichia coli and determine if a sexual strain of E. coli can become an effective tool for strain development.

Recombination and phenotype evolution dynamics of Helicobacter pylori in colonized hosts

The ample genetic diversity and variability of Helicobater pylori, and therefore its phenotypic evolution, relate not only to frequent mutation and selection but also to intra-specific recombination. Webb and Blaser applied a mathematical model to distinguish the role of selection and mutation for Lewis antigen phenotype evolution during long-term gastric colonization in infected animal hosts (mice and gerbils). To investigate the role of recombination in Lewis antigen phenotype evolution, we have developed a prior population dynamic by adding recombination term to the model. We simulate and interpret the new model simulation's results with a comparative analysis of biological aspects. The main conclusions are as follows: (i) the models and consequently the hosts with higher recombination rate require a longer time for stabilization; and (ii) recombination and mutation have opposite effects on the size of H. pylori populations with phenotypes in the range of the most-fit ones (i.e. those that have a selective advantage) due to natural selection, although both can increase phenotypic diversity.

Influence of mutation and recombination on HIV-1 in vitro fitness recovery

The understanding of the evolutionary processes underlying HIV-1 fitness recovery is fundamental for HIV-1 pathogenesis, antiretroviral treatment and vaccine design. It is known that HIV-1 can present very high mutation and recombination rates, however the specific contribution of these evolutionary forces in the "in vitro" viral fitness recovery has not been simultaneously quantified. To this aim, we analyzed substitution, recombination and molecular adaptation rates in a variety of HIV-1 biological clones derived from a viral isolate after severe population bottlenecks and a number of large population cell culture passages. These clones presented an overall but uneven fitness gain, mean of 3-fold, respect to the initial passage values. We found a significant relationship between the fitness increase and the appearance and fixation of mutations. In addition, these fixed mutations presented molecular signatures of positive selection through the accumulation of non-synonymous substitutions. Interestingly, viral recombination correlated with fitness recovery in most of studied viral quasispecies. The genetic diversity generated by these evolutionary processes was positively correlated with the viral fitness. We conclude that HIV-1 fitness recovery can be derived from the genetic heterogeneity generated through both mutation and recombination, and under diversifying molecular adaptation. The findings also suggest nonrandom evolutionary pathways for in vitro fitness recovery.

Extensive genomic variation within clonal bacterial groups resulted from homologous recombination

Due to divergence, genetic variation is generally believed to be high among distantly related strains, low among closely related ones and little or none within the same classified clonal groups. Several recent genome-wide studies, however, revealed that significant genetic variation resides in a considerable number of genes among strains with identical MLST (Multilocus sequence typing) types and much of the variation was introduced by homologous recombination. Recognizing and understanding genomic variation within clonal bacterial groups could shed new light on the evolutionary path of infectious agents and the emergence of particularly pathogenic or virulent variants. This commentary presents our recent contributions to this line of work.

Multidimensional epistasis and the transitory advantage of sex

Identifying and quantifying the benefits of sex and recombination is a long-standing problem in evolutionary theory. In particular, contradictory claims have been made about the existence of a benefit of recombination on high dimensional fitness landscapes in the presence of sign epistasis. Here we present a comparative numerical study of sexual and asexual evolutionary dynamics of haploids on tunably rugged model landscapes under strong selection, paying special attention to the temporal development of the evolutionary advantage of recombination and the link between population diversity and the rate of adaptation. We show that the adaptive advantage of recombination on static rugged landscapes is strictly transitory. At early times, an advantage of recombination arises through the possibility to combine individually occurring beneficial mutations, but this effect is reversed at longer times by the much more efficient trapping of recombining populations at local fitness peaks. These findings are explained by means of well-established results for a setup with only two loci. In accordance with the Red Queen hypothesis the transitory advantage can be prolonged indefinitely in fluctuating environments, and it is maximal when the environment fluctuates on the same time scale on which trapping at local optima typically occurs.

Recombination accelerates adaptation on a large-scale empirical fitness landscape in HIV-1

Recombination has the potential to facilitate adaptation. In spite of the substantial body of theory on the impact of recombination on the evolutionary dynamics of adapting populations, empirical evidence to test these theories is still scarce. We examined the effect of recombination on adaptation on a large-scale empirical fitness landscape in HIV-1 based on in vitro fitness measurements. Our results indicate that recombination substantially increases the rate of adaptation under a wide range of parameter values for population size, mutation rate and recombination rate. The accelerating effect of recombination is stronger for intermediate mutation rates but increases in a monotonic way with the recombination rates and population sizes that we examined. We also found that both fitness effects of individual mutations and epistatic fitness interactions cause recombination to accelerate adaptation. The estimated epistasis in the adapting populations is significantly negative. Our results highlight the importance of recombination in the evolution of HIV-I.

One of the most challenging issues in evolutionary biology concerns the question of why most organisms exchange genetic material with each other, e.g. during sexual reproduction. Gene shuffling can create genetic diversity that facilitates adaptation to new environments, but theory shows that this effect is highly dependent on how different genes interact in determining the fitness of an organism. Using a large data set of fitness values based on HIV-1, we provide evidence that shuffling of genetic material indeed raises the level of genetic diversity, and as a result accelerates adaptation. Our results also propose genetic shuffling as a mechanism utilized by HIV to accelerate the evolution of multi-drug-resistant strains.

The impact of natural transformation on adaptation in spatially structured bacterial populations

BACKGROUND

Recent studies have demonstrated that natural transformation and the formation of highly structured populations in bacteria are interconnected. In spite of growing evidence about this connection, little is known about the dynamics of natural transformation in spatially structured bacterial populations.

RESULTS

In this work, we model the interdependency between the dynamics of the bacterial gene pool and those of environmental DNA in space to dissect the effect of transformation on adaptation. Our model reveals that even with only a single locus under consideration, transformation with a free DNA fragment pool results in complex adaptation dynamics that do not emerge in previous models focusing only on the gene shuffling effect of transformation at multiple loci. We demonstrate how spatial restriction on population growth and DNA diffusion in the environment affect the impact of transformation on adaptation. We found that in structured bacterial populations intermediate DNA diffusion rates predominantly cause transformation to impede adaptation by spreading deleterious alleles in the population.

CONCLUSION

Overall, our model highlights distinctive evolutionary consequences of bacterial transformation in spatially restricted compared to planktonic bacterial populations.

Fitness landscapes emerging from pharmacodynamic functions in the evolution of multidrug resistance

Adaptive evolution often involves beneficial mutations at more than one locus. In this case, the trajectory and rate of adaptation is determined by the underlying fitness landscape, that is, the fitness values and mutational connectivity of all genotypes under consideration. Drug resistance, especially resistance to multiple drugs simultaneously, is also often conferred by mutations at several loci so that the concept of fitness landscapes becomes important. However, fitness landscapes underlying drug resistance are not static but dependent on drug concentrations, which means they are influenced by the pharmacodynamics of the drugs administered. Here, I present a mathematical framework for fitness landscapes of multidrug resistance based on Hill functions describing how drug concentrations affect fitness. I demonstrate that these 'pharmacodynamic fitness landscapes' are characterized by pervasive epistasis that arises through (i) fitness costs of resistance (even when these costs are additive), (ii) nonspecificity of resistance mutations to drugs, in particular cross-resistance, and (iii) drug interactions (both synergistic and antagonistic). In the latter case, reciprocal drug suppression may even lead to reciprocal sign epistasis, so that the doubly resistant genotype occupies a local fitness peak that may be difficult to access by evolution. Simulations exploring the evolutionary dynamics on some pharmacodynamic fitness landscapes with both constant and changing drug concentrations confirm the crucial role of epistasis in determining the rate of multidrug resistance evolution.

The evolution of natural competence: disentangling costs and benefits of sex in bacteria

One of the most challenging questions in evolutionary biology is how sex has evolved in the face of substantial fitness costs. In this study, we focus on the evolution of bacterial sex in the form of natural transformation, where cells take up exogenous DNA and integrate it into the genome. Besides the physiological cost of producing a DNA uptake system, transformation can potentially impose a genetic cost as a result of an overrepresentation of deleterious mutations in the extracellular DNA pool. On the other hand, the uptake of DNA can be beneficial not only because of genetic effects but also because of the immediate nutritional value of the DNA. To disentangle these fitness costs and benefits, we developed a mathematical model and competed three bacterial types during adaptation to a new environment: competent cells capable of DNA import and digestion; competent cells capable of DNA import, digestion, and recombination; and noncompetent cells. Our results indicate a complex interplay between several physiological and ecological factors, including the rate at which DNA is taken up, the rate of DNA decay in the medium, and the nutritional value of DNA. In finite populations, the recombining type is often favored through the Fisher-Muller effect.

Homologous recombination drives both sequence diversity and gene content variation in Neisseria meningitidis

The study of genetic and phenotypic variation is fundamental for understanding the dynamics of bacterial genome evolution and untangling the evolution and epidemiology of bacterial pathogens. Neisseria meningitidis (Nm) is among the most intriguing bacterial pathogens in genomic studies due to its dynamic population structure and complex forms of pathogenicity. Extensive genomic variation within identical clonal complexes (CCs) in Nm has been recently reported and suggested to be the result of homologous recombination, but the extent to which recombination contributes to genomic variation within identical CCs has remained unclear. In this study, we sequenced two Nm strains of identical serogroup (C) and multi-locus sequence type (ST60), and conducted a systematic analysis with an additional 34 Nm genomes. Our results revealed that all gene content variation between the two ST60 genomes was introduced by homologous recombination at the conserved flanking genes, and 94.25% or more of sequence divergence was caused by homologous recombination. Recombination was found in genes associated with virulence factors, antigenic outer membrane proteins, and vaccine targets, suggesting an important role of homologous recombination in rapidly altering the pathogenicity and antigenicity of Nm. Recombination was also evident in genes of the restriction and modification systems, which may undermine barriers to DNA exchange. In conclusion, homologous recombination can drive both gene content variation and sequence divergence in Nm. These findings shed new light on the understanding of the rapid pathoadaptive evolution of Nm and other recombinogenic bacterial pathogens.