The evolutionary origin of the universal distribution of mutation fitness effect

Long-range linkage effects in adapting sexual populations

In sexual populations, closely-situated genes have linked evolutionary fates, while genes spaced far in genome are commonly thought to evolve independently due to recombination. In the case where evolution depends essentially on supply of new mutations, this assumption has been confirmed by mathematical modeling. Here I examine it in the case of pre-existing genetic variation, where mutation is not important. A haploid population with [Formula: see text] genomes, [Formula: see text] loci, a fixed selection coefficient, and a small initial frequency of beneficial alleles [Formula: see text] is simulated by a Monte-Carlo algorithm. When the number of loci, L, is larger than a critical value of [Formula: see text] simulation demonstrates a host of linkage effects that decrease neither with the distance between loci nor the number of recombination crossovers. Due to clonal interference, the beneficial alleles become extinct at a fraction of loci [Formula: see text]. Due to a genetic background effect, the substitution rate varies broadly between loci, with the fastest value exceeding the one-locus limit by the factor of [Formula: see text] Thus, the far-situated parts of a long genome in a sexual population do not evolve as independent blocks. A potential link between these findings and the emergence of new Variants of Concern of SARS-CoV-2 is discussed.

Evolutionary implications of SARS-CoV-2 vaccination for the future design of vaccination strategies

Once the first SARS-CoV-2 vaccine became available, mass vaccination was the main pillar of the public health response to the COVID-19 pandemic. It was very effective in reducing hospitalizations and deaths. Here, we discuss the possibility that mass vaccination might accelerate SARS-CoV-2 evolution in antibody-binding regions compared to natural infection at the population level. Using the evidence of strong genetic variation in antibody-binding regions and taking advantage of the similarity between the envelope proteins of SARS-CoV-2 and influenza, we assume that immune selection pressure acting on these regions of the two viruses is similar. We discuss the consequences of this assumption for SARS-CoV-2 evolution in light of mathematical models developed previously for influenza. We further outline the implications of this phenomenon, if our assumptions are confirmed, for the future design of SARS-CoV-2 vaccination strategies.

Measurement of selection coefficients from genomic samples of adapting populations by computer modeling

The existing protocols of measuring the selection coefficients of loci neglect linkage effects existing between loci. This protocol is free from this limitation. The protocol inputs a set of DNA sequences at three time points, removes conserved sites, and estimates selection coefficients. If the user wishes to test the accuracy, it can ask the protocol to generate mock data by computer simulation of evolution. The main limitation is the need for sequence samples isolated from 30-100 populations adapting in parallel. For complete details on the use and execution of this protocol, please refer to Barlukova and Rouzine (2021).

Novelty Search Promotes Antigenic Diversity in Microbial Pathogens

Driven by host–pathogen coevolution, cell surface antigens are often the fastest evolving parts of a microbial pathogen. The persistent evolutionary impetus for novel antigen variants suggests the utility of novelty-seeking algorithms in predicting antigen diversification in microbial pathogens. In contrast to traditional genetic algorithms maximizing variant fitness, novelty-seeking algorithms optimize variant novelty. Here, we designed and implemented three evolutionary algorithms (fitness-seeking, novelty-seeking, and hybrid) and evaluated their performances in 10 simulated and 2 empirically derived antigen fitness landscapes. The hybrid walks combining fitness- and novelty-seeking strategies overcame the limitations of each algorithm alone, and consistently reached global fitness peaks. Thus, hybrid walks provide a model for microbial pathogens escaping host immunity without compromising variant fitness. Biological processes facilitating novelty-seeking evolution in natural pathogen populations include hypermutability, recombination, wide dispersal, and immune-compromised hosts. The high efficiency of the hybrid algorithm improves the evolutionary predictability of novel antigen variants. We propose the design of escape-proof vaccines based on high-fitness variants covering a majority of the basins of attraction on the fitness landscape representing all potential variants of a microbial antigen.

Loss and fixation of strongly favoured new variants: Understanding and extending Haldane's result via the Wright-Fisher model

The implications of Haldane's analysis for the fixation of beneficial alleles lies at the heart of much of 'population genetic thinking' and underlies many approaches that have been tailored to the detection of positive selection. Within the framework of a branching process, Haldane gave an approximation for the probability that fixation ultimately occurs when the selective advantage of a beneficial allele is small (≪1). Here, we make no use of branching processes. Rather, we work solely within a finite-population Wright-Fisher framework. We use this framework to analyse where Haldane's result applies, and extend Haldane's analysis. In particular, we present results for: (i) the domain of applicability of Haldane's analysis; (ii) the probability that loss occurs up to a given time; (iii) the probabilities that loss and fixation ultimately occur; (iv) an analytic approximation associated with the probability of loss and fixation ultimately occurring; (v) quantification of the crossover from weak to strong selection; (vi) determination of the number of invasive alleles that have a significant probability (>0.95) of invading a novel population. We note that the results obtained for (ii), (iii) and (iv) hold for an arbitrary initial number of mutations, and for selection that can be arbitrarily strong. Our results have fundamental implications for population and conservation genetics, and open up new avenues to identify traces of historically beneficial alleles through comparative genomics.

An evolution-based high-fidelity method of epistasis measurement: Theory and application to influenza

Linkage effects in a multi-locus population strongly influence its evolution. The models based on the traveling wave approach enable us to predict the average speed of evolution and the statistics of phylogeny. However, predicting statistically the evolution of specific sites and pairs of sites in the multi-locus context remains a mathematical challenge. In particular, the effects of epistasis, the interaction of gene regions contributing to phenotype, is difficult to predict theoretically and detect experimentally in sequence data. A large number of false-positive interactions arises from stochastic linkage effects and indirect interactions, which mask true epistatic interactions. Here we develop a proof-of-principle method to filter out false-positive interactions. We start by demonstrating that the averaging of haplotype frequencies over multiple independent populations is necessary but not sufficient for epistatic detection, because it still leaves high numbers of false-positive interactions. To compensate for the residual stochastic noise, we develop a three-way haplotype method isolating true interactions. The fidelity of the method is confirmed analytically and on simulated genetic sequences evolved with a known epistatic network. The method is then applied to a large sequence database of neurominidase protein of influenza A H1N1 obtained from various geographic locations to infer the epistatic network responsible for the difference between the pre-pandemic virus and the pandemic strain of 2009. These results present a simple and reliable technique to measure epistatic interactions of any sign from sequence data.

Interactions between genomic sites create a fitness landscape. The knowledge of topology and strength of interactions is vital for predicting the escape of viruses from drugs and immune response and their passing through fitness valleys. Many efforts have been invested into measuring these interactions from DNA sequence sets. Unfortunately, reproducibility of the results remains low due partly to a very small fraction of interaction pairs and partly to stochastic linkage noise masking true interactions. Here we propose a method to separate stochastic linkage and indirect interactions from epistatic interactions and apply it to influenza virus sequence data.