Fitness Effects Of New Mutations Research Articles

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.

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.

Self-replicating artificial neural networks give rise to universal evolutionary dynamics.

In evolutionary models, mutations are exogenously introduced by the modeler, rather than endogenously introduced by the replicator itself. We present a new deep-learning based computational model, the self-replicating artificial neural network (SeRANN). We train it to (i) copy its own genotype, like a biological organism, which introduces endogenous spontaneous mutations; and (ii) simultaneously perform a classification task that determines its fertility. Evolving 1,000 SeRANNs for 6,000 generations, we observed various evolutionary phenomena such as adaptation, clonal interference, epistasis, and evolution of both the mutation rate and the distribution of fitness effects of new mutations. Our results demonstrate that universal evolutionary phenomena can naturally emerge in a self-replicator model when both selection and mutation are implicit and endogenous. We therefore suggest that SeRANN can be applied to explore and test various evolutionary dynamics and hypotheses.

Changing fitness effects of mutations through long-term bacterial evolution.

The distribution of fitness effects of new mutations shapes evolution, but it is challenging to observe how it changes as organisms adapt. Using Escherichia coli lineages spanning 50,000 generations of evolution, we quantify the fitness effects of insertion mutations in every gene. Macroscopically, the fraction of deleterious mutations changed little over time whereas the beneficial tail declined sharply, approaching an exponential distribution. Microscopically, changes in individual gene essentiality and deleterious effects often occurred in parallel; altered essentiality is only partly explained by structural variation. The identity and effect sizes of beneficial mutations changed rapidly over time, but many targets of selection remained predictable because of the importance of loss-of-function mutations. Taken together, these results reveal the dynamic-but statistically predictable-nature of mutational fitness effects.

Strong purifying selection in haploid tissue-specific genes of Scots pine supports the masking theory.

The masking theory states that genes expressed in haploid stage will be under more efficient selection. In contrast, selection will be less efficient in genes expressed in diploid stage, where the fitness effects of recessive deleterious or beneficial mutations can be hidden from selection in heterozygous form. This difference can influence several evolutionary processes such as the maintenance of genetic variation, adaptation rate, and genetic load. Masking theory expectations have been confirmed in single-cell haploid and diploid organisms. However, in multicellular organisms, such as plants, the effects of haploid selection are not clear-cut. In plants, the great majority of studies indicating haploid selection have been carried out using male haploid tissues in angiosperms. Hence, evidence in these systems is confounded with the effects of sexual selection and intra-specific competition. Evidence from other plant groups is scarce and results show no support for the masking theory. Here we have used a gymnosperm Scots pine megagametophyte, a maternally-derived seed haploid tissue, and four diploid tissues to test the strength of purifying selection on a set of genes with tissue-specific expression. By using targeted resequencing data of those genes, we obtained estimates of genetic diversity, the site frequency spectrum of 0-fold and 4-fold sites, and inferred the distribution of fitness effects (DFE) of new mutations in haploid and diploid tissue-specific genes. Our results show that purifying selection is stronger for tissue-specific genes expressed in the haploid megagametophyte tissue, and that this signal of strong selection is not an artifact driven by high expression levels.

Mutation saturation for fitness effects at human CpG sites.

Whole exome sequences have now been collected for millions of humans, with the related goals of identifying pathogenic mutations in patients and establishing reference repositories of data from unaffected individuals. As a result, we are approaching an important limit, in which datasets are large enough that, in the absence of natural selection, every highly mutable site will have experienced at least one mutation in the genealogical history of the sample. Here, we focus on CpG sites that are methylated in the germline and experience mutations to T at an elevated rate of ~10-7 per site per generation; considering synonymous mutations in a sample of 390,000 individuals, ~ 99 % of such CpG sites harbor a C/T polymorphism. Methylated CpG sites provide a natural mutation saturation experiment for fitness effects: as we show, at nt sample sizes, not seeing a non-synonymous polymorphism is indicative of strong selection against that mutation. We rely on this idea in order to directly identify a subset of CpG transitions that are likely to be highly deleterious, including ~27 % of possible loss-of-function mutations, and up to 20 % of possible missense mutations, depending on the type of functional site in which they occur. Unlike methylated CpGs, most mutation types, with rates on the order of 10-8 or 10-9, remain very far from saturation. We discuss what these findings imply for interpreting the potential clinical relevance of mutations from their presence or absence in reference databases and for inferences about the fitness effects of new mutations.

What does the distribution of fitness effects of new mutations reflect? Insights from plants.

The distribution of fitness effects (DFE) of new mutations plays a central role in molecular evolution. It is therefore crucial to be able to estimate it accurately from genomic data and to understand the factors that shape it. After a rapid overview of available methods to characterize the fitness effects of mutations, we review what is known on the factors affecting them in plants. Available data indicate that life history traits (e.g. mating system and longevity) have a major effect on the DFE. By contrast, the impact of demography within species appears to be more limited. These results remain to be confirmed, and methods to estimate the joint evolution of demography, life history traits, and the DFE need to be developed.

Selection on Accessible Chromatin Regions in Capsella grandiflora.

Accurate estimates of genome-wide rates and fitness effects of new mutations are essential for an improved understanding of molecular evolutionary processes. Although eukaryotic genomes generally contain a large noncoding fraction, functional noncoding regions and fitness effects of mutations in such regions are still incompletely characterized. A promising approach to characterize functional noncoding regions relies on identifying accessible chromatin regions (ACRs) tightly associated with regulatory DNA. Here, we applied this approach to identify and estimate selection on ACRs in Capsella grandiflora, a crucifer species ideal for population genomic quantification of selection due to its favorable population demography. We describe a population-wide ACR distribution based on ATAC-seq data for leaf samples of 16 individuals from a natural population. We use population genomic methods to estimate fitness effects and proportions of positively selected fixations (α) in ACRs and find that intergenic ACRs harbor a considerable fraction of weakly deleterious new mutations, as well as a significantly higher proportion of strongly deleterious mutations than comparable inaccessible intergenic regions. ACRs are enriched for expression quantitative trait loci (eQTL) and depleted of transposable element insertions, as expected if intergenic ACRs are under selection because they harbor regulatory regions. By integrating empirical identification of intergenic ACRs with analyses of eQTL and population genomic analyses of selection, we demonstrate that intergenic regulatory regions are an important source of nearly neutral mutations. These results improve our understanding of selection on noncoding regions and the role of nearly neutral mutations for evolutionary processes in outcrossing Brassicaceae species.

Hunting for Beneficial Mutations: Conditioning on SIFT Scores When Estimating the Distribution of Fitness Effect of New Mutations.

The distribution of fitness effects (DFE) of new mutations is a key parameter of molecular evolution. The DFE can in principle be estimated by comparing the site frequency spectra (SFS) of putatively neutral and functional polymorphisms. Unfortunately, the DFE is intrinsically hard to estimate, especially for beneficial mutations because these tend to be exceedingly rare. There is therefore a strong incentive to find out whether conditioning on properties of mutations that are independent of the SFS could provide additional information. In the present study, we developed a new measure based on SIFT scores. SIFT scores are assigned to nucleotide sites based on their level of conservation across a multispecies alignment: the more conserved a site, the more likely mutations occurring at this site are deleterious, and the lower the SIFT score. If one knows the ancestral state at a given site, one can assign a value to new mutations occurring at the site based on the change of SIFT score associated with the mutation. We called this new measure δ. We show that properties of the DFE as well as the flux of beneficial mutations across classes covary with δ and, hence, that SIFT scores are informative when estimating the fitness effect of new mutations. In particular, conditioning on SIFT scores can help to characterize beneficial mutations.

Inferring Parameters of the Distribution of Fitness Effects of New Mutations When Beneficial Mutations Are Strongly Advantageous and Rare.

Characterizing the distribution of fitness effects (DFE) for new mutations is central in evolutionary genetics. Analysis of molecular data under the McDonald-Kreitman test has suggested that adaptive substitutions make a substantial contribution to between-species divergence. Methods have been proposed to estimate the parameters of the distribution of fitness effects for positively selected mutations from the unfolded site frequency spectrum (uSFS). Such methods perform well when beneficial mutations are mildly selected and frequent. However, when beneficial mutations are strongly selected and rare, they may make little contribution to standing variation and will thus be difficult to detect from the uSFS. In this study, I analyze uSFS data from simulated populations subject to advantageous mutations with effects on fitness ranging from mildly to strongly beneficial. As expected, frequent, mildly beneficial mutations contribute substantially to standing genetic variation and parameters are accurately recovered from the uSFS. However, when advantageous mutations are strongly selected and rare, there are very few segregating in populations at any one time. Fitting the uSFS in such cases leads to underestimates of the strength of positive selection and may lead researchers to false conclusions regarding the relative contribution adaptive mutations make to molecular evolution. Fortunately, the parameters for the distribution of fitness effects for harmful mutations are estimated with high accuracy and precision. The results from this study suggest that the parameters of positively selected mutations obtained by analysis of the uSFS should be treated with caution and that variability at linked sites should be used in conjunction with standing variability to estimate parameters of the distribution of fitness effects in the future.

PolyDFE: Inferring the Distribution of Fitness Effects and Properties of Beneficial Mutations from Polymorphism Data.

The possible evolutionary trajectories a population can follow is determined by the fitness effects of new mutations. Their relative frequencies are best specified through a distribution of fitness effects (DFE) that spans deleterious, neutral, and beneficial mutations. As such, the DFE is key to several aspects of the evolution of a population, and particularly the rate of adaptive molecular evolution (α). Inference of DFE from patterns of polymorphism and divergence has been a longstanding goal of evolutionary genetics.polyDFE provides a flexible statistical framework to estimate the DFE and α from site frequency spectrum (SFS) data. Several probability distributions can be fitted to the data to model the DFE. The method also jointly estimates a series of nuisance parameters that model the effect of unknown demography as well data imperfections, in particular possible errors in polarizing SNPs. This chapter is organized as a tutorial for polyDFE. We start by briefly reviewing the concept of DFE, α, and the principles underlying the method, and then provide an example using central chimpanzees data (Tataru et al., Genetics 207(3):1103-1119, 2017; Bataillon et al., Genome Biol Evol 7(4):1122-1132, 2015) to guide the user through the different steps of an analysis: formatting the data as input to polyDFE, fitting different models, obtaining estimates of parameters uncertainty and performing statistical tests, as well as model averaging procedures to obtain robust estimates of model parameters.

Deleterious variation shapes the genomic landscape of introgression.

While it is appreciated that population size changes can impact patterns of deleterious variation in natural populations, less attention has been paid to how gene flow affects and is affected by the dynamics of deleterious variation. Here we use population genetic simulations to examine how gene flow impacts deleterious variation under a variety of demographic scenarios, mating systems, dominance coefficients, and recombination rates. Our results show that admixture between populations can temporarily reduce the genetic load of smaller populations and cause increases in the frequency of introgressed ancestry, especially if deleterious mutations are recessive. Additionally, when fitness effects of new mutations are recessive, between-population differences in the sites at which deleterious variants exist creates heterosis in hybrid individuals. Together, these factors lead to an increase in introgressed ancestry, particularly when recombination rates are low. Under certain scenarios, introgressed ancestry can increase from an initial frequency of 5% to 30–75% and fix at many loci, even in the absence of beneficial mutations. Further, deleterious variation and admixture can generate correlations between the frequency of introgressed ancestry and recombination rate or exon density, even in the absence of other types of selection. The direction of these correlations is determined by the specific demography and whether mutations are additive or recessive. Therefore, it is essential that null models of admixture include both demography and deleterious variation before invoking other mechanisms to explain unusual patterns of genetic variation.

Biased Inference of Selection Due to GC-Biased Gene Conversion and the Rate of Protein Evolution in Flycatchers When Accounting for It.

The rate of recombination impacts on rates of protein evolution for at least two reasons: it affects the efficacy of selection due to linkage and influences sequence evolution through the process of GC-biased gene conversion (gBGC). We studied how recombination, via gBGC, affects inferences of selection in gene sequences using comparative genomic and population genomic data from the collared flycatcher (Ficedula albicollis). We separately analyzed different mutation categories (“strong”-to-“weak,” “weak-to-strong,” and GC-conservative changes) and found that gBGC impacts on the distribution of fitness effects of new mutations, and leads to that the rate of adaptive evolution and the proportion of adaptive mutations among nonsynonymous substitutions are underestimated by 22–33%. It also biases inferences of demographic history based on the site frequency spectrum. In light of this impact, we suggest that inferences of selection (and demography) in lineages with pronounced gBGC should be based on GC-conservative changes only. Doing so, we estimate that 10% of nonsynonymous mutations are effectively neutral and that 27% of nonsynonymous substitutions have been fixed by positive selection in the flycatcher lineage. We also find that gene expression level, sex-bias in expression, and the number of protein–protein interactions, but not Hill–Robertson interference (HRI), are strong determinants of selective constraint and rate of adaptation of collared flycatcher genes. This study therefore illustrates the importance of disentangling the effects of different evolutionary forces and genetic factors in interpretation of sequence data, and from that infer the role of natural selection in DNA sequence evolution.

Inbreeding depression maintained by recessive lethal mutations interacting with stabilizing selection on quantitative characters in a partially self-fertilizing population.

The bimodal distribution of fitness effects of new mutations and standing genetic variation, due to early-acting strongly deleterious recessive mutations and late-acting mildly deleterious mutations, is analyzed using the Kondrashov model for lethals (K), with either the infinitesimal model for selfing (IMS) or the Gaussian allele model (GAM) for quantitative genetic variance under stabilizing selection. In the combined models (KIMS and KGAM) high genomic mutation rates to lethals and weak stabilizing selection on many characters create strong interactions between early and late inbreeding depression, by changing the distribution of lineages selfed consecutively for different numbers of generations. Alternative stable equilibria can exist at intermediate selfing rates for a given set of parameters. Evolution of quantitative genetic variance under multivariate stabilizing selection can strongly influence the purging of nearly recessive lethals, and sometimes vice versa. If the selfing rate at the purging threshold for quantitative genetic variance in IMS or GAM alone exceeds that for nearly recessive lethals in K alone, then in KIMS and KGAM stabilizing selection causes selective interference with purging of lethals, increasing the mean number of lethals compared to K; otherwise, stabilizing selection causes selective facilitation in purging of lethals, decreasing the mean number of lethals.

Evolution of haploid-diploid life cycles when haploid and diploid fitnesses are not equal.

Many organisms spend a significant portion of their life cycle as haploids and as diploids (a haploid-diploid life cycle). However, the evolutionary processes that could maintain this sort of life cycle are unclear. Most previous models of ploidy evolution have assumed that the fitness effects of new mutations are equal in haploids and homozygous diploids, however, this equivalency is not supported by empirical data. With different mutational effects, the overall (intrinsic) fitness of a haploid would not be equal to that of a diploid after a series of substitution events. Intrinsic fitness differences between haploids and diploids can also arise directly, for example because diploids tend to have larger cell sizes than haploids. Here, we incorporate intrinsic fitness differences into genetic models for the evolution of time spent in the haploid versus diploid phases, in which ploidy affects whether new mutations are masked. Life-cycle evolution can be affected by intrinsic fitness differences between phases, the masking of mutations, or a combination of both. We find parameter ranges where these two selective forces act and show that the balance between them can favor convergence on a haploid-diploid life cycle, which is not observed in the absence of intrinsic fitness differences.

The Role of Deleterious Substitutions in Crop Genomes.

Populations continually incur new mutations with fitness effects ranging from lethal to adaptive. While the distribution of fitness effects of new mutations is not directly observable, many mutations likely either have no effect on organismal fitness or are deleterious. Historically, it has been hypothesized that a population may carry many mildly deleterious variants as segregating variation, which reduces the mean absolute fitness of the population. Recent advances in sequencing technology and sequence conservation-based metrics for inferring the functional effect of a variant permit examination of the persistence of deleterious variants in populations. The issue of segregating deleterious variation is particularly important for crop improvement, because the demographic history of domestication and breeding allows deleterious variants to persist and reach moderate frequency, potentially reducing crop productivity. In this study, we use exome resequencing of 15 barley accessions and genome resequencing of 8 soybean accessions to investigate the prevalence of deleterious single nucleotide polymorphisms (SNPs) in the protein-coding regions of the genomes of two crops. We conclude that individual cultivars carry hundreds of deleterious SNPs on average, and that nonsense variants make up a minority of deleterious SNPs. Our approach identifies known phenotype-altering variants as deleterious more frequently than the genome-wide average, suggesting that putatively deleterious variants are likely to affect phenotypic variation. We also report the implementation of a SNP annotation tool BAD_Mutations that makes use of a likelihood ratio test based on alignment of all currently publicly available Angiosperm genomes.

Field measurements of genotype by environment interaction for fitness caused by spontaneous mutations in Arabidopsis thaliana.

As the ultimate source of genetic diversity, spontaneous mutation is critical to the evolutionary process. The fitness effects of spontaneous mutations are almost always studied under controlled laboratory conditions rather than under the evolutionarily relevant conditions of the field. Of particular interest is the conditionality of new mutations-that is, is a new mutation harmful regardless of the environment in which it is found? In other words, what is the extent of genotype-environment interaction for spontaneous mutations? We studied the fitness effects of 25 generations of accumulated spontaneous mutations in Arabidopsis thaliana in two geographically widely separated field environments, in Michigan and Virginia. At both sites, mean total fitness of mutation accumulation lines exceeded that of the ancestors, contrary to the expected decrease in the mean due to new mutations but in accord with prior work on these MA lines. We observed genotype-environment interactions in the fitness effects of new mutations, such that the effects of mutations in Michigan were a poor predictor of their effects in Virginia and vice versa. In particular, mutational variance for fitness was much larger in Virginia compared to Michigan. This strong genotype-environment interaction would increase the amount of genetic variation maintained by mutation-selection balance.

The implications of small stem cell niche sizes and the distribution of fitness effects of new mutations in aging and tumorigenesis.

Somatic tissue evolves over a vertebrate's lifetime due to the accumulation of mutations in stem cell populations. Mutations may alter cellular fitness and contribute to tumorigenesis or aging. The distribution of mutational effects within somatic cells is not known. Given the unique regulatory regime of somatic cell division, we hypothesize that mutational effects in somatic tissue fall into a different framework than whole organisms; one in which there are more mutations of large effect. Through simulation analysis, we investigate the fit of tumor incidence curves generated using exponential and power‐law distributions of fitness effects (DFE) to known tumorigenesis incidence. Modeling considerations include the architecture of stem cell populations, that is, a large number of very small populations, and mutations that do and do not fix neutrally in the stem cell niche. We find that the typically quantified DFE in whole organisms is sufficient to explain tumorigenesis incidence. Further, deleterious mutations are predicted to accumulate via genetic drift, resulting in reduced tissue maintenance. Thus, despite there being a large number of stem cells throughout the intestine, its compartmental architecture leads to the accumulation of deleterious mutations and significant aging, making the intestinal stem cell niche a prime example of Muller's Ratchet.

Fitness effects of new mutations in Chlamydomonas reinhardtii across two stress gradients.

Most spontaneous mutations affecting fitness are likely to be deleterious, but the strength of selection acting on them might be impacted by environmental stress. Such stress‐dependent selection could expose hidden genetic variation, which in turn might increase the adaptive potential of stressed populations. On the other hand, this variation might represent a genetic load and thus lead to population extinction under stress. Previous studies to determine the link between stress and mutational effects on fitness, however, have produced inconsistent results. Here, we determined the net change in fitness in 29 genotypes of the green algae Chlamydomonas reinhardtii that accumulated mutations in the near absence of selection for approximately 1000 generations across two stress gradients, increasing NaCl and decreasing phosphate. We found mutational effects to be magnified under extremely stressful conditions, but such effects were specific both to the type of stress and to the genetic background. The detection of stress‐dependent fitness effects of mutations depended on accurately scaling relative fitness measures by generation times, thus offering an explanation for the inconsistencies among previous studies.

Accelerating Mutational Load Is Not Due to Synergistic Epistasis or Mutator Alleles in Mutation Accumulation Lines of Yeast.

Much of our knowledge about the fitness effects of new mutations has been gained from mutation accumulation (MA) experiments. Yet the fitness effect of single mutations is rarely measured in MA experiments. This raises several issues, notably for inferring epistasis for fitness. The acceleration of fitness decline in MA lines has been taken as evidence for synergistic epistasis, but establishing the role of epistasis requires measuring the fitness of genotypes carrying known numbers of mutations. Otherwise, accelerating fitness loss could be explained by increased genetic mutation rates. Here we segregated mutations accumulated over 4800 generations in haploid and diploid MA lines of the yeast Saccharomyces cerevisiae. We found no correspondence between an accelerated fitness decline and synergistic epistasis among deleterious mutations in haploid lines. Pairs of mutations showed no overall epistasis. Furthermore, several lines of evidence indicate that genetic mutation rates did not increase in the MA lines. Crucially, segregant fitness analyses revealed that MA accelerated in both haploid and diploid lines, even though the fitness of diploid lines was nearly constant during the MA experiment. This suggests that the accelerated fitness decline in haploids was caused by cryptic environmental factors that increased mutation rates in all lines during the last third of the lines' transfers. In addition, we provide new estimates of deleterious mutation rates, including lethal mutations, and highlight that nearly all the mutational load we observed was due to one or two mutations having a large effect on fitness.