Search PubMedSearch

SEARCH · Search PubMed

Results for “mutational fitness effects”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

The distribution of fitness effects of nonsynonymous mutations varies phylogenetically across animals.

The distribution of fitness effects (DFE) describes the selection coefficients of newly arising mutations and fundamentally influences population genetic processes. However, the extent and mechanisms of differences in the DFE for non-synonymous mutations have not been systematically investigated across species with divergent phylogenetic histories and ecologies. Here, we inferred the DFE in natural populations of 11 animal (sub)species, including humans, mice, fin whales, vaquitas, wolves, collared flycatchers, pied flycatchers, halictid bees, Drosophila, and mosquitoes. We found that mammals have a higher proportion of strongly deleterious mutations (defined as s≤-0.01; 22% to 47% in mammals; 0.0% to 5.4% in insects and birds) and a lower proportion of weakly deleterious mutations than insects and birds. Further, the DFE co-varies with phylogeny, such that the mean mutation effects are more similar in closely related species (Pagel's λ = 0.84, P = 0.01). Next, we investigated whether various summary statistics of the DFE were related to variation in life-history traits across these organisms. We found some support for genome size, body mass, and long-term effective population size being correlated with the DFE. Overall, our findings are consistent with predictions derived independently from the Fisher's Geometric Model (FGM), which defines organismal complexity as the number of phenotypes under selection. FGM predicts that mutations are more deleterious in complex organisms, while strongly deleterious mutations occur more frequently in smaller populations. Our study demonstrates strong phylogenetic signal in the evolution of a fundamental population genetics parameter, and proposes that, through mechanisms of epistasis, long-term population size and organismal complexity could be underlying variation in the DFE across animals.

Journal Article

Synonymous mutations in essential genes infrequently produce fitness effects in human cell lines.

The assumption that synonymous mutations are fitness-neutral is central to many foundational results in the fields of genetics, genomics, evolutionary biology, and medicine. However, recent results suggest synonymous mutations have pervasive and strong fitness effects. These vigorously debated studies in non-human model systems have even suggested that the proportion of synonymous mutations and their fitness effect sizes are similar to non-synonymous mutations. To probe the fitness effect of synonymous mutations, we utilized recent advances in base editing to test 8558 potential synonymous mutations in 128 highly essential genes in human cell lines. Importantly, our library design excluded splice-proximal sites, ensuring a direct test of codon-level synonymous effects independent of splicing disruption. We find that synonymous mutations rarely have fitness effects on growth, occurring around 37.9-fold (95% CI: 22.16-81.48-fold) less frequently than missense mutations. In this experimental context, these findings demonstrate that synonymous mutations impact cellular fitness far less frequently than missense mutations. These results deviate from earlier reports of widespread synonymous fitness effects in yeast, yet they align with recent prime editing data observed in other human cell lines.

Humans

Epistasis and the changing fitness landscapes of SARS-CoV-2.

Since its emergence in late 2019, millions of SARS-CoV-2 genomes have been generated as part of global efforts to monitor the evolution and spread of the virus. This unprecedented volume of data provides a unique opportunity to study viral evolution at unparalleled resolution. In particular, individual genomic sites can be observed to have mutated independently thousands of times. These mutation counts have been used to estimate site-specific mutation rates and fitness effects for most mutations across the viral genome. Here, we use these data to investigate how the landscape of mutational fitness costs has changed over the course of the pandemic. SARS-CoV-2 evolution over the past 6 years has been characterized by the emergence of distinct variants separated by long branches corresponding to evolutionary saltations involving up to 50 mutations. We compare inferred fitness landscapes of the Spike protein across these variants and find that shifts in the estimated effects of non-synonymous mutations are linked to genetic differences between them. Sites with altered fitness costs are enriched near positions where the genetic backgrounds differ. To explain the observed changes, we introduce a model with pairwise epistatic interactions between mutations and residues that differ between variants. This model is able to explain about half of the variance in the shifts of fitness effects and suggests that each mismatch between variants substantially alters mutation effects at typically 1 to 3 additional positions.

SARS-CoV-2

Fitness effects of EMS-induced mutations on the X chromosome of Drosophila melanogaster. II. Hemizygous fitness effects.

X chromosomes mutagenized with EMS were tested for their effects on the fitness of hemizygous carriers. The tests were carried out in populations in which treated and untreated X chromosomes segregated from matings between males and attached-X females; the populations were maintained for several generations, during which time changes in the frequencies of the treated and untreated chromosomes were observed. From the rates at which the frequencies changed, the fitness effects of the treated chromosomes were determined. It was found that flies hemizygous for a mutagenized chromosome were 1.7% less fit per mM EMS treatment than those hemizygous for an untreated chromosome. Since the same flies were only 0.5% per mM less viable than their untreated counterparts, the total fitness effect of an X chromosome carrying EMS-induced mutants is three to four times greater than its viability effect. By comparing the heterozygous effect of a mutagenized X chromosome on fitness with the corresponding hemizygous effect, the dominance value for the chromosome is estimated to be about 0.25.

Animals

Fitness effects of EMS-induced mutations on the X chromosome of Drosophila melanogaster. I. Viability effects and heterozygous fitness effects.

Drosophila melanogaster X chromosomes were mutagenized by feeding males sucrose solutions containing ethyl methanesulfonate (EMS); the concentrations of EMS in the food were 2.5 mM, 5.0 mM, and 10.0 mM. Chromosomes were exposed to the mutagen up to three times by treating males in succeeding generations. After treatment, the effective exposures were 2.5, 5.0, 7.5, 10.0, 15.0, and 30.0 mM EMS. X chromosomes treated in this manner were tested for effects on fitness in both hemizygous and heterozygous conditions, and for effects on viability in hemizygous and homozygous conditions. In addition, untreated X chromosomes were available for study. The viability and heterozygous fitness effects are presented in this paper, and the hemizygous fitness effects are discussed in the accompanying one (MITCHELL and SIMMONS 1977). Hemizygous and homozygous viability effects were measured by segregation tests in vial cultures. For hemizygous males, viability was reduced 0.5 percent per mM EMS treatment; for homozygous females, it was reduced 0.7% per mM treatment. The decline in viability appeared to be a linear function of EMS dose. The viabilities of males and females were strongly correlated. Heterozygous fitness effects were measured by monitoring changes in the frequencies of treated and untreated X chromosomes in discrete generation populations which, through the use of an X-Y translocation, maintained them only in heterozygous condition. Flies that were heterozygous for a treated chromosome were found to be 0.4% less fit per mM EMS than flies heterozygous for an untreated one.

Animals

Sexual selection, genomic evolution and population fitness in Drosophila pseudoobscura.

Sexual selection shapes the genome in unique ways. It is also likely to have significant fitness consequences, such as purging deleterious mutations from the genome or conversely maintaining genetic load in a population via sexual conflict. Here, we examined what the influence of sexual selection has on genomic variation potentially underlying population fitness using experimentally evolved Drosophila pseudoobscura populations. Sexual selection was manipulated by keeping replicate lines in elevated polyandry or strict monogamy for approximately 200 generations followed by individual-based sequencing. Using pi (π), fixation index (Fst)and recombination rate measures, we confirmed signatures of selection were not dispersed but mainly localized to the third and X chromosome. Overall mutational load was similar between lines but our analysis of the distribution of fitness effects revealed considerable variation between lines and chromosomes. Furthermore, we found that the distribution of transposable elements differs between the lines, with a higher load in monogamous lines. Our results suggest that complex interactions between purifying selection and sexual conflict are shaping the genome, particularly on chromosome 3 and the sex chromosome; sexual selection influences divergence across chromosomes but in a more complex way than proposed by simple 'purging' of deleterious loci.

Animals

A gene-based model of fitness and its implications for genetic variation: Linkage disequilibrium.

A widely used model of the effects of mutations on fitness (the "sites" model) assumes that heterozygous recessive or partially recessive deleterious mutations at different sites in a gene complement each other, similarly to mutations in different genes. However, the general lack of complementation between major effect allelic mutations suggests an alternative possibility, which we term the "gene" model. This assumes that a pair of heterozygous deleterious mutations in trans behave effectively as homozygotes, so that the fitnesses of trans heterozygotes are lower than those of cis heterozygotes. We examine the properties of the two different models, using both analytical and simulation methods. We show that the gene model predicts positive linkage disequilibrium (LD) between deleterious variants within the coding sequence, under conditions when the sites model predicts zero or slightly negative LD. We also show that focussing on rare variants when examining patterns of LD, especially with Lewontin's´ measure, is likely to produce misleading results with respect to inferences concerning the causes of the sign of LD. Synergistic epistasis between pairs of mutations was also modeled; it is less likely to produce negative LD under the gene model than the sites model. The theoretical results are discussed in relation to patterns of LD in natural populations of several species.

complementation

The shape of fitness functions and the distribution of mutational effect sizes jointly limit adaptation by regulatory mutations.

Mutations in gene regulatory regions have been shown to play a role in rapid adaptation, but the factors determining their contribution are largely unknown. Here, using the metabolic enzyme cytosine deaminase of budding yeast, we examine whether adaptation to 5-fluorocytosine, which requires reduced cytosine deamination and can readily arise from amino acid substitutions, may be reached by single promoter mutations. We generated all single-nucleotide substitutions and indels in the FCY1 promoter and assayed the resulting mutants in presence of 5-fluorocytosine. This revealed that no promoter mutation is sufficient for adaptation to occur. We next investigated how this inaccessibility of adaptation arises by combining large-scale expression measurements with the experimental characterization of the corresponding expression-fitness function. These experiments showed that the shape of this function precludes single promoter mutations from being adaptive. Although 24% of mutations significantly affect expression, the fitness curve is flat around wild-type level. As such, adaptation can only emerge from a severe reduction of expression, which cannot occur from a single mutation in the promoter. Our results show that the contribution of regulatory mutations to rapid adaptation depends not only on the distribution of mutational effect sizes on expression level but also on the shape of the function linking fitness to expression levels.

Promoter Regions, Genetic

A model for background selection in non-equilibrium populations.

In many taxa, levels of genetic diversity are observed to vary along their genome. The framework of background selection models this variation in terms of linkage to constrained sites, and recent applications have been able to explain a large portion of the variation in human genomes. However, these studies have also yielded conflicting results, stemming from two key limitations. First, existing models are inaccurate in a critical region of parameter space (), where the local reduction in diversity is sharpest. Second, they assume a constant population size over time. Here, we develop predictions for diversity under background selection based on the Hill-Robertson system of two-locus statistics, which allows for population size changes. We treat the joint effect of multiple selected loci independently, but we show that interference among them is well captured through local rescaling of mutation, recombination and selection in an iterative procedure that converges quickly. We further accommodate existing background selection theory to non-equilibrium demography, bridging the gap between weak and strong selection. Simulations show that our predictions are accurate across the entire range of selection coefficients. We characterize the temporal dynamics of linked selection under population size changes and demonstrate that patterns of diversity can be misinterpreted by other models. Specifically, biases due to the incorrect assumption of equilibrium carry over to downstream inferences of the distribution of fitness effects and deleterious mutation rate. Jointly modeling demography and linked selection therefore improves our understanding of the genomic landscape of diversity, which will help refine inferences of linked selection in humans and other species.

Journal Article

Bayesian inference of fitness landscapes via tree-structured branching processes.

MOTIVATION: The complex dynamics of cancer evolution, driven by mutation and selection, underlies the molecular heterogeneity observed in tumors. The evolutionary histories of tumors of different patients can be encoded as mutation trees and reconstructed in high resolution from single-cell sequencing data, offering crucial insights for studying fitness effects of and epistasis among mutations. Existing models, however, either fail to separate mutation and selection or neglect the evolutionary histories encoded by the tumor phylogenetic trees. RESULTS: We introduce FiTree, a tree-structured multi-type branching process model with epistatic fitness parameterization and a Bayesian inference scheme to learn fitness landscapes from single-cell tumor mutation trees. Through simulations, we demonstrate that FiTree outperforms state-of-the-art methods in inferring the fitness landscape underlying tumor evolution. Applying FiTree to a single-cell acute myeloid leukemia dataset, we identify epistatic fitness effects consistent with known biological findings and quantify uncertainty in predicting future mutational events. The new model unifies probabilistic graphical models of cancer progression with population genetics, offering a principled framework for understanding tumor evolution and informing therapeutic strategies. AVAILABILITY AND IMPLEMENTATION: The Python package FiTree and the analysis workflows are available at https://github.com/cbg-ethz/FiTree.

Bayes Theorem

Evolutionary consequences of domestication on the selective effects of new amino acid changing mutations in canids.

The domestication of wild canids led to dogs no longer living in the wild but instead residing alongside humans. Extreme changes in behavior and diet associated with domestication may have led to the relaxation of the selective pressure on traits that may be less important in the domesticated context. Thus, here we hypothesize that strongly deleterious mutations may have become less deleterious in domesticated populations. We test this hypothesis by estimating the distribution of fitness effects (DFE) for new amino acid changing mutations using whole-genome sequence data from 24 gray wolves and 61 breed dogs. We find that the DFE is strikingly similar across canids, with 26-28% of new amino acid changing mutations being neutral/nearly neutral (|s| < 1e-5), and 41-48% under strong purifying selection (|s| > 1e-2). Our results are robust to different model assumptions suggesting that the DFE is stable across short evolutionary timescales, even in the face of putative drastic changes in the selective pressure caused by artificial selection during domestication and breed formation. On par with previous works describing DFE evolution, our data indicate that the DFE of amino acid changing mutations depends more strongly on genome structure and organismal characteristics, and less so on shifting selective pressures or environmental factors. Given the constant DFE and previous data showing that genetic variants that differentiate wolf and dog populations are enriched in regulatory elements, we speculate that domestication may have had a larger impact on regulatory variation than on amino acid changing mutations.

Journal Article

A novel dimerization site in non-structural protein 5A of hepatitis C virus regulates viral replication fitness.

We previously found that high genome replication fitness of the hepatitis C virus (HCV) was associated with severe disease in immunocompromised patients. Elevated replication fitness was mediated by accumulation of mutations in the replication enhancing domain (ReED) within domain (D) 2 of non-structural protein (NS) 5A. NS5A is a partially unstructured phosphoprotein lacking enzymatic activity but fulfilling a key role in HCV replication due to interacting with various cellular and viral proteins. It can exist in a variety of dimeric and oligomeric conformations mediated by NS5A D1 with clinically approved NS5A inhibitors proposed to exert their antiviral function by fixing these dimers in distinct conformations. In this study, we aimed at elucidating the ReED's mode of action. AlphaFold modelling indicated a so far unrecognized NS5A dimerization site in the ReED. Indeed, split nano luciferase assays revealed a significantly stronger NS5A dimerization of high replicator ReED variants, suggesting that high replication fitness is mediated by enforcement of NS5A self-interaction. This hypothesis was supported by the effect of low dose (1 pM) NS5A inhibitor treatment, increasing replication fitness and phenocopying the effects of ReED mutations. Furthermore, we found that HCV isolate JFH1, replicating with very high efficiency, is completely resistant to the regulatory function of the ReED. Chimeric replicons composed of ReED resistant JFH1 and the ReED sensitive isolate J6 identified NS3 helicase and NS5B polymerase as critical genetic elements mediating ReED sensitivity/resistance. Our data overall suggest that the ReED in NS5A is a negative regulator of HCV replication fitness with dimerization releasing the inhibitory interaction with helicase and/or polymerase, thereby likely facilitating initiation of RNA synthesis.

Viral Nonstructural Proteins

Stability of phage T4 lysozymes. II. Unfolding with guanidinium chloride.

The denaturation by guanidinium chloride of three phage lysozymes (wild type and two mutants) was investigated. The study of solvent denaturation permitted the investigation of the relative stabilities of the proteins at neutral pH, in contrast to thermal denaturation studies reported earlier which could only be performed in acid pH. The results were interpreted assuming that the free energy of solution of proteins is a linear function of denaturant concentration. Using standard thermodynamic formulas this permits the calculation of the stabilities of the three proteins in the absence of guanidinium chloride. The single point mutation Trp 138 leads to Tyr leads to relatively large changes in stability and the interaction of the protein with guanidinium chloride. The changes associated with the subsequent double mutation, Trp 126 leads to Tyr, Trp 158 leads to Tyr, are much smaller indicating a relatively smooth adjustment of the protein structure to the changed side chains. Models of the structural effects of point mutations are discussed. It is found that the mutation at position 138 does not fit a model in which the effect of a substitution is to introduce an energetic strain in the structure. It does fit a model in which there is a partial unravelling of the structure as a result of the mutation. However, there are no changes in the backbone circular dichroism spectra associated with the mutation. The two observations are not necessarily in conflict. Further physical studies are required for the resolution of the problem.

Guanidines

SARS-CoV-2 genomic diversity and within-host evolution in individuals with persistent infection in the UK: an observational, longitudinal, population-based surveillance study.

BACKGROUND: Persistent SARS-CoV-2 infections in hospitalised immunocompromised individuals are known to facilitate accelerated within-host viral evolution, potentially contributing to the emergence of highly divergent variants. However, little is known about the evolutionary dynamics and transmission risks of persistent infections in the general population. We aimed to characterise the within-host evolution of SARS-CoV-2 during persistent infections identified through a large community surveillance study. METHODS: We used data from the Office for National Statistics COVID-19 Infection Survey (ONS-CIS), a large-scale, longitudinal, population-based surveillance study conducted in the UK from April, 2020, to March, 2023. For this analysis, we focused on infections with high viral load (cycle threshold &#x2264;30) and available genome sequences, from seven major SARS-CoV-2 lineages (alpha, delta, BA.1, BA.2, BA.4, BA.5, and XBB). ONS-CIS participants were randomly selected from the general population and tested regularly by RT-PCR, regardless of symptoms. We defined persistent infections as those with sustained or rebounding high viral RNA titres for 26 days or longer. We examined associated host characteristics and used raw sequence data to identify de novo mutations and estimate within-host synonymous and non-synonymous evolutionary rates across the SARS-CoV-2 genome. FINDINGS: Between Nov 2, 2020, and March 21, 2023, we identified 576 persistent infections with at least two sequences, including 11 alpha, 106 delta, 102 BA.1, 204 BA.2, 16 BA.4, 133 BA.5, and 4 XBB. Persistent infections were more common in males than females (p<0&#xb7;0001) and individuals older than 60 years (p=0&#xb7;0027). The median within-host genome-wide evolutionary rate was 7&#xb7;9&#x2009;&#xd7;&#x2009;10-4 substitutions per site per year (IQR 7&#xb7;0-9&#xb7;0&#x2009;&#xd7;&#x2009;10-4), with high inter-individual variability driven largely by non-synonymous mutations, particularly in the N-terminal and receptor-binding domains of the spike protein. Longer infection duration was associated with higher evolutionary rates, while no associations were found with age, sex, vaccination status, previous infection, or virus lineage. We found no clear evidence of transmission beyond the first month of infection in any of the 84 persistent infections lasting 56 days or longer. In total, we identified 379 recurrent mutations, including many with known or predicted negative fitness effects and low prevalence at the population level, as well as de novo reversions to the Wuhan-Hu-1 reference sequence, which were likely under positive selection within those individuals. INTERPRETATION: This study highlights the heterogeneous nature of within-host SARS-CoV-2 evolution in individuals with persistent infection in the community. Notably, a small subset of persistent infections with high viral loads underwent accelerated viral evolution or recurrently acquired hallmark mutations found in novel variants. In addition, onward transmission from a persistent infection during the later stages of infection is likely to be rare. These insights have important implications for prioritising genomic surveillance and managing patients with persistent infections. FUNDING: Department of Health and Social Care.

Humans

Gene-level complexity explains genome-wide variation in the distribution of fitness effects.

The distribution of fitness effects (DFE)-describing how harmful, neutral, or beneficial new mutations are-is central to understanding how populations evolve. Although the DFE varies across genomes and species, it remains unclear which aspects of genomic organization drive this variation. Here, we inferred gene-level selective constraints across the genomes of Mus musculus castaneus, Drosophila melanogaster and Saccharomyces cerevisiae using a combination of population genetics and machine learning trained on diverse gene features. Many gene features were predictive of selective constraint, with conservation, gene structure, and expression being the most informative. These selective constraints delineated gene classes with distinct DFEs. Genes with higher connectivity and expression-features reflecting how many traits a gene influences-experienced stronger and less dispersed deleterious effects with increasing selective constraint. Between species, the rate of adaptation decreased with increasing organismal complexity, whereas across the genome it did not decrease monotonically with selective constraint, but tended to be higher at intermediate levels. While between-species comparisons of DFE parameters were less consistent with predictions of Fisher's geometric model (FGM) based on organismal complexity, variation in DFE parameters across the genome aligned more closely with FGM when complexity was considered at the gene level. Our results suggest that gene-level complexity, captured by genomic feature proxies, provides a more informative definition of complexity for DFE variation than organism-level labels, and highlight the value of using gene features collectively to link genomic architecture, fitness landscapes, and patterns of molecular evolution.

Animals

Type 1 interferon perturbates clonal competition by reshaping human blood development.

Inflammation accelerates evolutionary dynamics of hematopoietic stem cells (HSCs) in clonal hematopoiesis and myeloid neoplasms. We studied HSCs, progenitors and immune cells from patients with myeloproliferative neoplasms at baseline and following interferon-&#x3b1; (IFN&#x3b1;) treatment, the only therapy to deplete mutated stem cells. We deployed single-cell multiomics methods that distinguish the IFN&#x3b1; effects on mutated stem cells from the admixed wild-type HSCs, with respect to their differentiation, transcriptomes, immunophenotypes and chromatin accessibility. IFN&#x3b1; simultaneously activated HSCs into two polarized states: a lymphoid progenitor expansion associated with an anti-inflammatory state and an inflammatory myeloid progenitor state derived from HSCs. The augmented lymphoid differentiation balanced the typical myeloproliferative-neoplasm-induced myeloid bias, associated with normalized blood counts. Somatic mutations modified the effects of IFN&#x3b1; on HSC differentiation and cell cycle entry rates. Clonal fitness upon IFN&#x3b1; exposure was due to resistance of CALR- or JAK2-mutated stem cells to differentiate into inflammatory myeloid progenitors.

Journal Article

The distribution of fitness effects varies phylogenetically across animals.

The distribution of fitness effects (DFE) describes the selection coefficients () of newly arising mutations and fundamentally influences population genetic processes. However, the extent and mechanisms of DFE variation have not been systematically investigated across species with divergent phylogenetic histories and ecological functions. Here, we inferred the DFE in natural populations of eleven animal (sub)species, including humans, mice, fin whales, vaquitas, wolves, collared flycatchers, pied flycatchers, halictid bees, Drosophila, and mosquitoes. We find that the DFE co-varies with phylogeny, where the expected mutation effects are more similar in closely related species (). Additionally, mammals have a higher proportion of strongly deleterious mutations (22% to 47% in mammals; 0.0% to 5.4% in insects and birds) and a lower proportion of weakly deleterious mutations than insects and birds. Population size is significantly negatively correlated with the expected impact of new deleterious mutations (), and the proportion of new beneficial mutations (). These findings align with Fisher's Geometric Model (FGM), which defines organismal complexity as the number of phenotypes under selection. Consistent with the FGM's predictions, we observe that mutations are more deleterious in complex organisms, while beneficial mutations occur more frequently in smaller populations to compensate for the drift load. Our study demonstrates strong phylogenetic constraints in the evolution of a fundamental population genetics parameter, and proposes that, through mechanisms of global epistasis, long-term population size and organismal complexity drive variation in the DFE across animals.

Fisher&#x2019;s geometric model

Gene duplication in tetraploid fish: model for gene silencing at unlinked duplicated loci.

Several groups of fishes, including salmonids and catastomids, appear to have originated through genome duplication events. However, these two groups retain approximately 50% of the loci examined as functioning duplicates, despite the passage of 50 million years or more of mutation and selection. Although other effects are not excluded, this apparently slow rate of duplicate silencing can be explained in terms of the effects of selection against defective double homozygotes to unlinked duplicates. We have derived a computer simulation of genetic drift that affords direct evaluation of the effects of population size (N), mutation rate (micron), initial allele frequencies, back mutation, fitness, and time on the probability of fixation for null alleles at unlinked duplicate loci. The results show that this probability is approximately linearly related to population size for N greater than or equal to 10(3). Specifically, for naive populations, the time for 50% probability of gene silencing is approximately equal to 15N + micron-3/4 generations. The retention of 50% of the loci as functional duplicates may therefore result from the large effective size of salmonid and catastomid populations. The results also show that, under most conditions for populations of 2000--3000 or larger, unlinked duplicate loci will be sustained in the functional state longer than tandem (linked) duplicates and hence are available for evolution of new functions for a longer time.

Alleles