Fitness Effects Of Mutations Research Articles

Two fitness inference schemes compared using allele frequencies from 1,068,391 sequences sampled in the UK during the COVID-19 pandemic.

Throughout the course of the SARS-CoV-2 pandemic, genetic variation has contributed to the spread and persistence of the virus. For example, various mutations have allowed SARS-CoV-2 to escape antibody neutralization or to bind more strongly to the receptors that it uses to enter human cells. Here, we compared two methods that estimate the fitness effects of viral mutations using the abundant sequence data gathered over the course of the pandemic. Both approaches are grounded in population genetics theory but with different assumptions. One approach, tQLE, features an epistatic fitness landscape and assumes that alleles are nearly in linkage equilibrium. Another approach, MPL, assumes a simple, additive fitness landscape, but allows for any level of correlation between alleles. We characterized differences in the distributions of fitness values inferred by each approach and in the ranks of fitness values that they assign to sequences across time. We find that in a large fraction of weeks the two methods are in good agreement as to their top-ranked sequences, \textit{i.e.} as to which sequences observed that week are most fit. We also find that agreement between the ranking of sequences varies with genetic unimodality in the population in a given week.&#xD.

Reconciling theories of dominance with the relative rates of adaptive substitution on sex chromosomes and autosomes.

The dominance of beneficial mutations is a key evolutionary parameter affecting the rate and genetic basis of adaptation, yet it is notoriously difficult to estimate. A leading method to infer it is to compare the relative rates of adaptive substitution for X-linked and autosomal genes, which-according to a classic model by Charlesworth et al. (1987)-is a simple function of the dominance of new beneficial mutations. Recent evidence that rates of adaptive substitution are faster for X-linked genes implies, accordingly, that beneficial mutations are usually recessive. However, this conclusion is incompatible with leading theories of dominance, which predict that beneficial mutations tend to be dominant or overdominant with respect to fitness. To address this incompatibility, we use Fisher's geometric model to predict the distribution of fitness effects of new mutations and the relative rates of positively selected substitution on the X and autosomes. Previous predictions of faster-X theory emerge as a special case of our model in which the phenotypic effects of mutations are small relative to the distance to the phenotypic optimum. But as mutational effects become large relative to the optimum, we observe an elevated tempo of positively selected substitutions on the X relative to the autosomes across a broader range of dominance conditions, including those predicted by theories of dominance. Our results imply that, contrary to previous models, dominant and overdominant beneficial mutations can plausibly generate patterns of faster-X adaptation. We discuss resulting implications for genomic studies of adaptation and inferences of dominance.

Further Evidence for Strong Nonneutrality of Yeast Synonymous Mutations.

Although synonymous mutations are commonly assumed neutral or nearly so, recent years have seen reports of fitness effects of synonymous mutations detected under laboratory conditions. In a previous study, we used genome editing to construct thousands of yeast mutants each carrying a synonymous or nonsynonymous mutation in one of 21 genes, and discovered that most synonymous and most nonsynonymous mutations are deleterious. A concern was raised that this observation could be caused by the fitness effects of potential CRISPR/Cas9 off-target edits and/or secondary mutations, and an experiment that would be refractory to such effects was proposed. Using genome sequencing, we here show that no CRISPR/Cas9 off-target editing occurred, although some mutants did carry secondary mutations. Analysis of mutants with negligible effects from secondary mutations and new data collected from the proposed experiment confirms the original conclusion. These findings, along with other reports of fitness effects of synonymous mutations from both case and systematic studies, necessitate a paradigm shift from assuming (near) neutrality of synonymous mutations.

HERMES: Holographic Equivariant neuRal network model for Mutational Effect and Stability prediction.

Predicting the stability and fitness effects of amino-acid mutations in proteins is a cornerstone of biological discovery and engineering. Various experimental techniques have been developed to measure mutational effects, providing us with extensive datasets across a diverse range of proteins. By training on these data, machine learning approaches have advanced significantly in predicting mutational effects. Here, we introduce HERMES, a 3D rotationally equivariant structure-based neural network model for mutation effect prediction. Pre-trained to predict amino-acid propensities from their surrounding 3D structure atomic environments, HERMES can be efficiently fine-tuned to predict mutational effects, thanks to its symmetry-aware parameterization of the output space. Benchmarking against other models demonstrates that HERMES often outperforms or matches their performance in predicting mutation effects on stability, binding, and fitness, using either computationally or experimentally resolved protein structures. HERMES offers a versatile suit of tools for evaluating mutation effects and can be easily fine-tuned for specific predictive objectives using our open-source code.

Predicting the first steps of evolution in randomly assembled communities

Microbial communities can self-assemble into highly diverse states with predictable statistical properties. However, these initial states can be disrupted by rapid evolution of the resident strains. When a new mutation arises, it competes for resources with its parent strain and with the other species in the community. This interplay between ecology and evolution is difficult to capture with existing community assembly theory. Here, we introduce a mathematical framework for predicting the first steps of evolution in large randomly assembled communities that compete for substitutable resources. We show how the fitness effects of new mutations and the probability that they coexist with their parent depends on the size of the community, the saturation of its niches, and the metabolic overlap between its members. We find that successful mutations are often able to coexist with their parent strains, even in saturated communities with low niche availability. At the same time, these invading mutants often cause extinctions of metabolically distant species. Our results suggest that even small amounts of evolution can produce distinct genetic signatures in natural microbial communities.

Emergent time scales of epistasis in protein evolution

We introduce a data-driven epistatic model of protein evolution, capable of generating evolutionary trajectories spanning very different time scales reaching from individual mutations to diverged homologs. Our in silico evolution encompasses random nucleotide mutations, insertions and deletions, and models selection using a fitness landscape, which is inferred via a generative probabilistic model for protein families. We show that the proposed framework accurately reproduces the sequence statistics of both short-time (experimental) and long-time (natural) protein evolution, suggesting applicability also to relatively data-poor intermediate evolutionary time scales, which are currently inaccessible to evolution experiments. Our model uncovers a highly collective nature of epistasis, gradually changing the fitness effect of mutations in a diverging sequence context, rather than acting via strong interactions between individual mutations. This collective nature triggers the emergence of a long evolutionary time scale, separating fast mutational processes inside a given sequence context, from the slow evolution of the context itself. The model quantitatively reproduces epistatic phenomena such as contingency and entrenchment, as well as the loss of predictability in protein evolution observed in deep mutational scanning experiments of distant homologs. It thereby deepens our understanding of the interplay between mutation and selection in shaping protein diversity and functions, allows one to statistically forecast evolution, and challenges the prevailing independent-site models of protein evolution, which are unable to capture the fundamental importance of epistasis.

Hill-Robertson interference may bias the inference of fitness effects of new mutations in highly selfing species.

The accurate estimation of the distribution of fitness effects (DFE) of new mutations is critical for population genetic inference but remains a challenging task. While various methods have been developed for DFE inference using the site frequency spectrum of putatively neutral and selected sites, their applicability in species with diverse life history traits and complex demographic scenarios is not well understood. Selfing is common among eukaryotic species and can lead to decreased effective recombination rates, increasing the effects of selection at linked sites, including interference between selected alleles. We employ forward simulations to investigate the limitations of current DFE estimation approaches in the presence of selfing and other model violations, such as linkage, departures from semidominance, population structure, and uneven sampling. We find that distortions of the site frequency spectrum due to Hill-Robertson interference in highly selfing populations lead to mis-inference of the deleterious DFE of new mutations. Specifically, when inferring the distribution of selection coefficients, there is an overestimation of nearly neutral and strongly deleterious mutations and an underestimation of mildly deleterious mutations when interference between selected alleles is pervasive. In addition, the presence of cryptic population structure with low rates of migration and uneven sampling across subpopulations leads to the false inference of a deleterious DFE skewed towards effectively neutral/mildly deleterious mutations. Finally, the proportion of adaptive substitutions estimated at high rates of selfing is substantially overestimated. Our observations apply broadly to species and genomic regions with little/no recombination and where interference might be pervasive.

Energy landscapes of peptide-MHC binding.

Molecules of the Major Histocompatibility Complex (MHC) present short protein fragments on the cell surface, an important step in T cell immune recognition. MHC-I molecules process peptides from intracellular proteins; MHC-II molecules act in antigen-presenting cells and present peptides derived from extracellular proteins. Here we show that the sequence-dependent energy landscapes of MHC-peptide binding encode class-specific nonlinearities (epistasis). MHC-I has a smooth landscape with global epistasis; the binding energy is a simple deformation of an underlying linear trait. This form of epistasis enhances the discrimination between strong-binding peptides. In contrast, MHC-II has a rugged landscape with idiosyncratic epistasis: binding depends on detailed amino acid combinations at multiple positions of the peptide sequence. The form of epistasis affects the learning of energy landscapes from training data. For MHC-I, a low-complexity problem, we derive a simple matrix model of binding energies that outperforms current models trained by machine learning. For MHC-II, higher complexity prevents learning by simple regression methods. Epistasis also affects the energy and fitness effects of mutations in antigen-derived peptides (epitopes). In MHC-I, large-effect mutations occur predominantly in anchor positions of strong-binding epitopes. In MHC-II, large effects depend on the background epitope sequence but are broadly distributed over the epitope, generating a bigger target for escape mutations due to loss of presentation. Together, our analysis shows how an energy landscape of protein-protein binding constrains the target of escape mutations from T cell immunity, linking the complexity of the molecular interactions to the dynamics of adaptive immune response.

Variable fitness effects of bacteriophage resistance mutations in Escherichia coli: implications for phage therapy.

This study isolated more than 50 phage-resistant mutants from both in vitro and in vivo conditions, exposing an extra-intestinal pathogenic Escherichia coli strain to a single virulent phage. The characterization of these clones revealed several key findings: (1) mutations occurring during phage treatment affect the same pathways as those identified in vitro; (2) the resistance mechanisms are associated with the modification of two cell-wall components, with one involving receptor deletion (phage-specific mechanism) and the other, less frequent, involving receptor masking (phage-nonspecific mechanism); (3) an in vivo virulence assay demonstrated that the absence of the receptor abolishes virulence while masking the receptor preserves it; and (4) clones with a resistance mechanism nonspecific to a particular phage can remain susceptible to other phages. This supports the idea of incorporating diverse phages into therapeutic cocktails designed to collectively target both wild-type and phage-resistant strains, including those with resistance mechanisms nonspecific to a phage.

Evolution of gene regulatory networks by means of selection and random genetic drift.

The evolution of a population by means of genetic drift and natural selection operating on a gene regulatory network (GRN) of an individual has not been scrutinized in depth. Thus, the relative importance of various evolutionary forces and processes on shaping genetic variability in GRNs is understudied. In this study, we implemented a simulation framework, called EvoNET, that simulates forward-in-time the evolution of GRNs in a population. The fitness effect of mutations is not constant, rather fitness of each individual is evaluated on the phenotypic level, by measuring its distance from an optimal phenotype. Each individual goes through a maturation period, where its GRN may reach an equilibrium, thus deciding its phenotype. Afterwards, individuals compete to produce the next generation. We examine properties of the GRN evolution, such as robustness against the deleterious effect of mutations and the role of genetic drift. We are able to confirm previous hypotheses regarding the effect of mutations and we provide new insights on the interplay between random genetic drift and natural selection.

Direct inference of the distribution of fitness effects of spontaneous mutations from recombinant inbred Caenorhabditis elegans mutation accumulation lines.

The distribution of fitness effects of new mutations plays a central role in evolutionary biology. Estimates of the distribution of fitness effect from experimental mutation accumulation lines are compromised by the complete linkage disequilibrium between mutations in different lines. To reduce the linkage disequilibrium, we constructed 2 sets of recombinant inbred lines from a cross of 2 Caenorhabditis elegans mutation accumulation lines. One set of lines ("RIAILs") was intercrossed for 10 generations prior to 10 generations of selfing; the second set of lines ("RILs") omitted the intercrossing. Residual linkage disequilibrium in the RIAILs is much less than in the RILs, which affects the inferred distribution of fitness effect when the sets of lines are analyzed separately. The best-fit model estimated from all lines (RIAILs + RILs) infers a large fraction of mutations with positive effects (∼40%); models that constrain mutations to have negative effects fit much worse. The conclusion is the same using only the RILs. For the RIAILs, however, models that constrain mutations to have negative effects fit nearly as well as models that allow positive effects. When mutations in high linkage disequilibrium are pooled into haplotypes, the inferred distribution of fitness effect becomes increasingly negative-skewed and leptokurtic. We conclude that the conventional wisdom-most mutations have effects near 0, a handful of mutations have effects that are substantially negative, and mutations with positive effects are very rare-is likely correct, and that unless it can be shown otherwise, estimates of the distribution of fitness effect that infer a substantial fraction of mutations with positive effects are likely confounded by linkage disequilibrium.

Quantifying the strength of viral fitness trade-offs between hosts: a meta-analysis of pleiotropic fitness effects

Abstract The range of hosts a given virus can infect is widely presumed to be limited by fitness trade-offs between alternative hosts. These fitness trade-offs may arise naturally due to antagonistic pleiotropy if mutations that increase fitness in one host tend to decrease fitness in alternate hosts. Yet there is also growing recognition that positive pleiotropy may be more common than previously appreciated. With positive pleiotropy, mutations have concordant fitness effects such that a beneficial mutation can simultaneously increase fitness in different hosts, providing a genetic mechanism by which selection can overcome fitness trade-offs. How readily evolution can overcome fitness trade-offs therefore depends on the overall distribution of mutational fitness effects between hosts, including the relative frequency of antagonistic versus positive pleiotropy. We therefore conducted a systematic meta-analysis of the pleiotropic fitness effects of viral mutations reported in different hosts. Our analysis indicates that while both antagonistic and positive pleiotropy are common, fitness effects are overall positively correlated between hosts and unconditionally beneficial mutations are not uncommon. Moreover, the relative frequency of antagonistic versus positive pleiotropy may simply reflect the underlying frequency of beneficial and deleterious mutations in individual hosts. Given a mutation is beneficial in one host, the probability that it is deleterious in another host is roughly equal to the probability that any mutation is deleterious, suggesting there is no natural tendency toward antagonistic pleiotropy. The widespread prevalence of positive pleiotropy suggests that many fitness trade-offs may be readily overcome by evolution given the right selection pressures.

Environmental memory alters the fitness effects of adaptive mutations in fluctuating environments.

Evolution in a static laboratory environment often proceeds via large-effect beneficial mutations that may become maladaptive in other environments. Conversely, natural settings require populations to endure environmental fluctuations. A sensible assumption is that the fitness of a lineage in a fluctuating environment is the time average of its fitness over the sequence of static conditions it encounters. However, transitions between conditions may pose entirely new challenges, which could cause deviations from this time average. To test this, we tracked hundreds of thousands of barcoded yeast lineages evolving in static and fluctuating conditions and subsequently isolated 900 mutants for pooled fitness assays in 15 environments. Here we find that fitness in fluctuating environments indeed often deviates from the time average, leading to fitness non-additivity. Moreover, closer examination reveals that fitness in one component of a fluctuating environment is often strongly influenced by the previous component. We show that this environmental memory is especially common for mutants with high variance in fitness across tested environments. We use a simple mathematical model and whole-genome sequencing to propose mechanisms underlying this effect, including lag time evolution and sensing mutations. Our results show that environmental fluctuations impact fitness and suggest that variance in static environments can explain these impacts.

Mapping mutational fitness effects across the coxsackievirus B3 proteome reveals distinct profiles of mutation tolerability.

RNA viruses have notoriously high mutation rates due to error-prone replication by their RNA polymerase. However, natural selection concentrates variability in a few key viral proteins. To test whether this stems from different mutation tolerance profiles among viral proteins, we measured the effect of >40,000 non-synonymous mutations across the full proteome of coxsackievirus B3 as well as >97% of all possible codon deletions in the nonstructural proteins. We find significant variation in mutational tolerance within and between individual viral proteins, which correlated with both general and protein-specific structural and functional attributes. Furthermore, mutational fitness effects remained stable across cell lines, suggesting selection pressures are mostly conserved across environments. In addition to providing a rich dataset for understanding virus biology and evolution, our results illustrate that incorporation of mutational tolerance data into druggable pocket discovery can aid in selecting targets with high barriers to drug resistance.

Parallel HIV-1 evolutionary dynamics in humans and rhesus macaques who develop broadly neutralizing antibodies.

Human immunodeficiency virus (HIV)-1 exhibits remarkable genetic diversity. For this reason, an effective HIV-1 vaccine must elicit antibodies that can neutralize many variants of the virus. While broadly neutralizing antibodies (bnAbs) have been isolated from HIV-1 infected individuals, a general understanding of the virus-antibody coevolutionary processes that lead to their development remains incomplete. We performed a quantitative study of HIV-1 evolution in two individuals who developed bnAbs. We observed strong selection early in infection for mutations affecting HIV-1 envelope glycosylation and escape from autologous strain-specific antibodies, followed by weaker selection for bnAb resistance later in infection. To confirm our findings, we analyzed data from rhesus macaques infected with viruses derived from the same two individuals. We inferred remarkably similar fitness effects of HIV-1 mutations in humans and macaques. Moreover, we observed a striking pattern of rapid HIV-1 evolution, consistent in both humans and macaques, that precedes the development of bnAbs. Our work highlights strong parallels between infection in rhesus macaques and humans, and it reveals a quantitative evolutionary signature of bnAb development.

Spontaneous single-nucleotide substitutions and microsatellite mutations have distinct distributions of fitness effects.

The fitness effects of new mutations determine key properties of evolutionary processes. Beneficial mutations drive evolution, yet selection is also shaped by the frequency of small-effect deleterious mutations, whose combined effect can burden otherwise adaptive lineages and alter evolutionary trajectories and outcomes in clonally evolving organisms such as viruses, microbes, and tumors. The small effect sizes of these important mutations have made accurate measurements of their rates difficult. In microbes, assessing the effect of mutations on growth can be especially instructive, as this complex phenotype is closely linked to fitness in clonally evolving organisms. Here, we perform high-throughput time-lapse microscopy on cells from mutation-accumulation strains to precisely infer the distribution of mutational effects on growth rate in the budding yeast, Saccharomyces cerevisiae. We show that mutational effects on growth rate are overwhelmingly negative, highly skewed towards very small effect sizes, and frequent enough to suggest that deleterious hitchhikers may impose a significant burden on evolving lineages. By using lines that accumulated mutations in either wild-type or slippage repair-defective backgrounds, we further disentangle the effects of 2 common types of mutations, single-nucleotide substitutions and simple sequence repeat indels, and show that they have distinct effects on yeast growth rate. Although the average effect of a simple sequence repeat mutation is very small (approximately 0.3%), many do alter growth rate, implying that this class of frequent mutations has an important evolutionary impact.

Experimental Evolution Studies in Φ6 Cystovirus.

Experimental evolution studies, in which biological populations are evolved in a specific environment over time, can address questions about the nature of spontaneous mutations, responses to selection, and the origins and maintenance of novel traits. Here, we review more than 30 years of experimental evolution studies using the bacteriophage (phage) Φ6 cystovirus. Similar to many lab-studied bacteriophages, Φ6 has a high mutation rate, large population size, fast generation time, and can be genetically engineered or cryogenically frozen, which facilitates its rapid evolution in the laboratory and the subsequent characterization of the effects of its mutations. Moreover, its segmented RNA genome, outer membrane, and capacity for multiple phages to coinfect a single host cell make Φ6 a good non-pathogenic model for investigating the evolution of RNA viruses that infect humans. We describe experiments that used Φ6 to address the fitness effects of spontaneous mutations, the consequences of evolution in the presence of coinfection, the evolution of host ranges, and mechanisms and consequences of the evolution of thermostability. We highlight open areas of inquiry where further experimentation on Φ6 could inform predictions for pathogenic viruses.

Predicting the First Steps of Evolution in Randomly Assembled Communities.

Microbial communities can self-assemble into highly diverse states with predictable statistical properties. However, these initial states can be disrupted by rapid evolution of the resident strains. When a new mutation arises, it competes for resources with its parent strain and with the other species in the community. This interplay between ecology and evolution is difficult to capture with existing community assembly theory. Here, we introduce a mathematical framework for predicting the first steps of evolution in large randomly assembled communities that compete for substitutable resources. We show how the fitness effects of new mutations and the probability that they coexist with their parent depends on the size of the community, the saturation of its niches, and the metabolic overlap between its members. We find that successful mutations are often able to coexist with their parent strains, even in saturated communities with low niche availability. At the same time, these invading mutants often cause extinctions of metabolically distant species. Our results suggest that even small amounts of evolution can produce distinct genetic signatures in natural microbial communities.

Mutations in yeast are deleterious on average regardless of the degree of adaptation to the testing environment.

The role of spontaneous mutations in evolution depends on the distribution of their effects on fitness. Despite a general consensus that new mutations are deleterious on average, a handful of mutation accumulation experiments in diverse organisms instead suggest that beneficial and deleterious mutations can have comparable fitness impacts, i.e. the product of their respective rates and effects can be roughly equal. We currently lack a general framework for predicting when such a pattern will occur. One idea is that beneficial mutations will be more evident in genotypes that are not well adapted to the testing environment. We tested this prediction experimentally in the laboratory yeast Saccharomyces cerevisiae by allowing nine replicate populations to adapt to novel environments with complex sets of stressors. After >1000 asexual generations interspersed with 41 rounds of sexual reproduction, we assessed the mean effect of induced mutations on yeast growth in both the environment to which they had been adapting and the alternative novel environment. The mutations were deleterious on average, with the severity depending on the testing environment. However, we found no evidence that the adaptive match between genotype and environment is predictive of mutational fitness effects.

The distribution of fitness effects during adaptive walks using a simple genetic network.

The tempo and mode of adaptation depends on the availability of beneficial alleles. Genetic interactions arising from gene networks can restrict this availability. However, the extent to which networks affect adaptation remains largely unknown. Current models of evolution consider additive genotype-phenotype relationships while often ignoring the contribution of gene interactions to phenotypic variance. In this study, we model a quantitative trait as the product of a simple gene regulatory network, the negative autoregulation motif. Using forward-time genetic simulations, we measure adaptive walks towards a phenotypic optimum in both additive and network models. A key expectation from adaptive walk theory is that the distribution of fitness effects of new beneficial mutations is exponential. We found that both models instead harbored distributions with fewer large-effect beneficial alleles than expected. The network model also had a complex and bimodal distribution of fitness effects among all mutations, with a considerable density at deleterious selection coefficients. This behavior is reminiscent of the cost of complexity, where correlations among traits constrain adaptation. Our results suggest that the interactions emerging from genetic networks can generate complex and multimodal distributions of fitness effects.