Genetic Constraints on Protein Evolution

ABSTRACT

Evolution requires the generation and optimization of new traits (“adaptation”) and involves the selection of mutations that improve cellular function. These mutations were assumed to arise by selection of neutral mutations present at all times in the population. Here we review recent evidence that indicates that deleterious mutations are more frequent in the population than previously recognized and that these mutations play a significant role in protein evolution through continuous positive selection. Positively selected mutations include adaptive mutations, i.e. mutations that directly affect enzymatic function, and compensatory mutations, which suppress the pleiotropic effects of adaptive mutations. Compensatory mutations are by far the most frequent of the two and would allow potentially adaptive but deleterious mutations to persist long enough in the population to be positively selected during episodes of adaptation. Compensatory mutations are, by definition, context-dependent and thus constrain the paths available for evolution. This provides a mechanistic basis for the examples of highly constrained evolutionary landscapes and parallel evolution reported in natural and experimental populations. The present review article describes these recent advances in the field of protein evolution and discusses their implications for understanding the genetic basis of disease and for protein engineering in vitro.     

Rate of adaptive peak shifts with partial genetic robustness

How adaptive evolution occurs with individually deleterious but jointly beneficial mutations has been one of the major problems in population genetics theory. Adaptation in this case is commonly described as a population's escape from a local peak to a higher peak on Sewall Wright's fitness landscape. Recent molecular genetic and computational studies have suggested that genetic robustness can facilitate such peak shifts. If phenotypic expressions of new mutations are suppressed under genetic robustness, mutations that are otherwise deleterious can accumulate in the population as neutral variants. When the robustness is perturbed by an environmental change or a major mutation, these variants become exposed to natural selection. It is argued that this process promotes adaptation because allelic combinations enriched under genetic robustness can then be positively selected. Here, I propose simple two- and three-locus models of adaptation with partial genetic robustness as suggested by recent studies. The waiting time until the fixation of an adaptive haplotype was observed in stochastic simulations and compared to the expectation without robustness. It is shown that peak shifts can be delayed or accelerated depending on the conditions of genetic robustness. The evolutionary significance of these processes is discussed.     

Evidence for hitchhiking of deleterious mutations within the human genome

Deleterious mutations present a significant obstacle to adaptive evolution. Deleterious mutations can inhibit the spread of linked adaptive mutations through a population; conversely, adaptive substitutions can increase the frequency of linked deleterious mutations and even result in their fixation. To assess the impact of adaptive mutations on linked deleterious mutations, we examined the distribution of deleterious and neutral amino acid polymorphism in the human genome. Within genomic regions that show evidence of recent hitchhiking, we find fewer neutral but a similar number of deleterious SNPs compared to other genomic regions. The higher ratio of deleterious to neutral SNPs is consistent with simulated hitchhiking events and implies that positive selection eliminates some deleterious alleles and increases the frequency of others. The distribution of disease-associated alleles is also altered in hitchhiking regions. Disease alleles within hitchhiking regions have been associated with auto-immune disorders, metabolic diseases, cancers, and mental disorders. Our results suggest that positive selection has had a significant impact on deleterious polymorphism and may be partly responsible for the high frequency of certain human disease alleles.     

Haldane's sieve and adaptation from the standing genetic variation

We consider populations that adapt to a sudden environmental change by fixing alleles found at mutation-selection balance. In particular, we calculate probabilities of fixation for previously deleterious alleles, ignoring the input of new mutations. We find that "Haldane's sieve"--the bias against the establishment of recessive beneficial mutations--does not hold under these conditions. Instead probabilities of fixation are generally independent of dominance. We show that this result is robust to patterns of sex expression for both X-linked and autosomal loci. We further show that adaptive evolution is invariably slower at X-linked than autosomal loci when evolution begins from mutation-selection balance. This result differs from that obtained when adaptation uses new mutations, a finding that may have some bearing on recent attempts to distinguish between hitchhiking and background selection by contrasting the molecular population genetics of X-linked vs. autosomal loci. Last, we suggest a test to determine whether adaptation used new mutations or previously deleterious alleles from the standing genetic variation.     

The rate of adaptation in asexuals

I study the population genetics of adaptation in asexuals. I show that the rate of adaptive substitution in an asexual species or nonrecombining chromosome region is a bell-shaped function of the mutation rate: at some point, increasing the mutation rate decreases the rate of substitution. Curiously, the mutation rate that maximizes the rate of adaptation depends solely on the strength of selection against deleterious mutations. In particular, adaptation is fastest when the genomic rate of mutation, U, equals the harmonic mean of selection coefficients against deleterious mutations, where we assume that selection for favorable alleles is milder than that against deleterious ones. This simple result is independent of the shape of the distribution of effects among favorable and deleterious mutations, population size, and the action of clonal interference. In the course of this work, I derive an approximation to the probability of fixation of a favorable mutation in an asexual genome or nonrecombining chromosome region in which both favorable and deleterious mutations occur.     

The cost of evolved constitutive lac gene expression is usually,but not always,maintained during evolution of generalist populations

Beneficial mutations can become costly following an environmental change. Compensatory mutations can relieve these costs, while not affecting the selected function, so that the benefits are retained if the environment shifts back to be similar to the one in which the beneficial mutation was originally selected. Compensatory mutations have been extensively studied in the context of antibiotic resistance, responses to specific genetic perturbations, and in the determination of interacting gene network components. Few studies have focused on the role of compensatory mutations during more general adaptation, especially as the result of selection in fluctuating environments where adaptations to different environment components may often involve trade‐offs. We examine whether costs of a mutation in lacI, which deregulated the expression of the lac operon in evolving populations of Escherichia coli bacteria, were compensated. This mutation occurred in multiple replicate populations selected in environments that fluctuated between growth on lactose, where the mutation was beneficial, and on glucose, where it was deleterious. We found that compensation for the cost of the lacI mutation was rare, but, when it did occur, it did not negatively affect the selected benefit. Compensation was not more likely to occur in a particular evolution environment. Compensation has the potential to remove pleiotropic costs of adaptation, but its rarity indicates that the circumstances to bring about the phenomenon may be peculiar to each individual or impeded by other selected mutations.     

Frequent beneficial mutations during single-colony serial transfer of Streptococcus pneumoniae

The appearance of new mutations within a population provides the raw material for evolution. The consistent decline in fitness observed in classical mutation accumulation studies has provided support for the long-held view that deleterious mutations are more common than beneficial mutations. Here we present results of a study using a mutation accumulation design with the bacterium Streptococcus pneumoniae in which the fitness of the derived populations increased. This rise in fitness was associated specifically with adaptation to survival during brief stationary phase periods between single-colony population bottlenecks. To understand better the population dynamics behind this unanticipated adaptation, we developed a maximum likelihood model describing the processes of mutation and stationary-phase selection in the context of frequent population bottlenecks. Using this model, we estimate that the rate of beneficial mutations may be as high as 4.8×10(-4) events per genome for each time interval corresponding to the pneumococcal generation time. This rate is several orders of magnitude higher than earlier estimates of beneficial mutation rates in bacteria but supports recent results obtained through the propagation of small populations of Escherichia coli. Our findings indicate that beneficial mutations may be relatively frequent in bacteria and suggest that in S. pneumoniae, which develops natural competence for transformation, a steady supply of such mutations may be available for sampling by recombination.     

Compensatory evolution in diploid populations

Compensatory mutations are individually deleterious but harmless in appropriate combinations either at more than two sites within a gene or on separate genes. Considering that dominance effects of selection and heterodimer formation of gene products may affect the rate of compensatory evolution, we investigate compensatory neutral mutation models for diploid populations. Our theoretical analysis on the average time until fixation of compensatory mutations shows that these factors play an important role in reducing the fixation time of compensatory mutations if mutation rates are not low. Compensatory evolution of heterodimers is shown to occur more easily if the deleterious effects of single mutants are recessive.     

Evolution by Gene Duplication and Compensatory Advantageous Mutations

Relaxation of selective constraint is thought to play an important role for evolution by gene duplication, in connection with compensatory advantageous mutant substitutions. Models were investigated by incorporating gene duplication by unequal crossing over, selection, mutation and random genetic drift into Monte Carlo simulations. Compensatory advantageous mutations were introduced, and simulations were carried out with and without relaxation, when genes are redundant on chromosomes. Relaxation was introduced by assuming that deleterious mutants have no effect on fitness, so long as one or more genes free of such mutations remain in the array. Compensatory mutations are characterized by the intermediate deleterious step of their substitutions, and therefore relaxation by gene redundancy is important. Through extensive Monte Carlo simulations, it was found that compensatory mutant substitutions require relaxation in addition to gene duplication, when mutant effects are large. However when mutant effects are small, such that the product of selection coefficient and population size is around unity, evolution by compensatory mutation is enhanced by gene duplication even without relaxation.     

Mistranslation Reduces Mutation Load in Evolving Proteins through Negative Epistasis with DNA Mutations

Translational errors during protein synthesis cause phenotypic mutations that are several orders of magnitude more frequent than DNA mutations. Such phenotypic mutations may affect adaptive evolution through their interactions with DNA mutations. To study how mistranslation may affect the adaptive evolution of evolving proteins, we evolved populations of green fluorescent protein (GFP) in either high-mistranslation or low-mistranslation Escherichia coli hosts. In both hosts, we first evolved GFP under purifying selection for the ancestral phenotype green fluorescence, and then under directional selection toward the new phenotype yellow fluorescence. High-mistranslation populations evolved modestly higher yellow fluorescence during each generation of evolution than low-mistranslation populations. We demonstrate by high-throughput sequencing that elevated mistranslation reduced the accumulation of deleterious DNA mutations under both purifying and directional selection. It did so by amplifying the fitness effects of deleterious DNA mutations through negative epistasis with phenotypic mutations. In contrast, mistranslation did not affect the incidence of beneficial mutations. Our findings show that phenotypic mutations interact epistatically with DNA mutations. By reducing a population’s mutation load, mistranslation can affect an important determinant of evolvability.     

How does adaptation sweep through the genome? Insights from long-term selection experiments

A major goal in evolutionary biology is to understand the origins and fates of adaptive mutations. Natural selection may act to increase the frequency of de novo beneficial mutations, or those already present in the population as standing genetic variation. These beneficial mutations may ultimately reach fixation in a population, or they may stop increasing in frequency once a particular phenotypic state has been achieved. It is not yet well understood how different features of population biology, and/or different environmental circumstances affect these adaptive processes. Experimental evolution is a promising technique for studying the dynamics of beneficial alleles, as populations evolving in the laboratory experience natural selection in a replicated, controlled manner. Whole-genome sequencing, regularly obtained over the course of sustained laboratory selection, could potentially reveal insights into the mutational dynamics that most likely occur in natural populations under similar circumstances. To date, only a few evolution experiments for which whole-genome data are available exist. This review describes results from these resequenced laboratory-selected populations, in systems with and without sexual recombination. In asexual systems, adaptation from new mutations can be studied, and results to date suggest that the complete, unimpeded fixation of these mutations is not always observed. In sexual systems, adaptation from standing genetic variation can be studied, and in the admittedly few examples we have, the complete fixation of standing variants is not always observed. To date, the relative frequency of adaptation from new mutations versus standing variation has not been tested using a single experimental system, but recent studies using Caenorhabditis elegans and Saccharomyces cerevisiae suggest that this a realistic future goal.     

Compensatory mutations are repeatable and clustered within proteins

Compensatory mutations improve fitness in genotypes that contain deleterious mutations but have no beneficial effects otherwise. As such, compensatory mutations represent a very specific form of epistasis. We show that intragenic compensatory mutations occur non-randomly over gene sequence. Compensatory mutations are more likely to appear at some sites than others. Moreover, the sites of compensatory mutations are more likely than expected by chance to be near the site of the original deleterious mutation. Furthermore, compensatory mutations tend to occur more commonly in certain regions of the protein even when controlling for clustering around the site of the deleterious mutation. These results suggest that compensatory evolution at the protein level is partially predictable and may be convergent.     

Emergent neutrality in adaptive asexual evolution

In nonrecombining genomes, genetic linkage can be an important evolutionary force. Linkage generates interference interactions, by which simultaneously occurring mutations affect each other's chance of fixation. Here, we develop a comprehensive model of adaptive evolution in linked genomes, which integrates interference interactions between multiple beneficial and deleterious mutations into a unified framework. By an approximate analytical solution, we predict the fixation rates of these mutations, as well as the probabilities of beneficial and deleterious alleles at fixed genomic sites. We find that interference interactions generate a regime of emergent neutrality: all genomic sites with selection coefficients smaller in magnitude than a characteristic threshold have nearly random fixed alleles, and both beneficial and deleterious mutations at these sites have nearly neutral fixation rates. We show that this dynamic limits not only the speed of adaptation, but also a population's degree of adaptation in its current environment. We apply the model to different scenarios: stationary adaptation in a time-dependent environment and approach to equilibrium in a fixed environment. In both cases, the analytical predictions are in good agreement with numerical simulations. Our results suggest that interference can severely compromise biological functions in an adapting population, which sets viability limits on adaptive evolution under linkage.     

Deleterious Passengers in Adapting Populations

Most new mutations are deleterious and are eventually eliminated by natural selection. But in an adapting population, the rapid amplification of beneficial mutations can hinder the removal of deleterious variants in nearby regions of the genome, altering the patterns of sequence evolution. Here, we analyze the interactions between beneficial “driver” mutations and linked deleterious “passengers” during the course of adaptation. We derive analytical expressions for the substitution rate of a deleterious mutation as a function of its fitness cost, as well as the reduction in the beneficial substitution rate due to the genetic load of the passengers. We find that the fate of each deleterious mutation varies dramatically with the rate and spectrum of beneficial mutations and the deleterious substitution rate depends nonmonotonically on the population size and the rate of adaptation. By quantifying this dependence, our results allow us to estimate which deleterious mutations will be likely to fix and how many of these mutations must arise before the progress of adaptation is significantly reduced.     

Adaptive co‐evolution of mitochondria and the Y‐chromosome: A resolution to conflict between evolutionary opponents

In most species with motile sperm, male fertility depends upon genes located on the Y‐chromosome and in the mitochondrial genome. Coordinated adaptive evolution for the function of male fertility between genes on the Y and the mitochondrion is hampered by their uniparental inheritance in opposing sexes: The Y‐chromosome is inherited uniparentally, father to son, and the mitochondrion is inherited maternally, mother to offspring. Preserving male fertility is problematic, because maternal inheritance permits mitochondrial mutations advantageous to females, but deleterious to male fertility, to accumulate in a population. Although uniparental inheritance with sex‐restricted adaptation also affects genes on the Y‐chromosome, females lack a Y‐chromosome and escape the potential maladaptive consequences of male‐limited selection. Evolutionary models have shown that mitochondrial mutations deleterious to male fertility can be countered by compensatory evolution of Y‐linked mutations that restore it. However, direct adaptive coevolution of Y‐ and mitochondrial gene combinations has not yet been mathematically characterized. We use population genetic models to show that adaptive coevolution of Y and mitochondrial genes are possible when Y‐mt gene combinations have positive effects on male fertility and populations are inbred.     

Intermolecular Interactions Drive Protein Adaptive and Coadaptive Evolution at Both Species and Population Levels

Proteins are the building blocks for almost all the functions in cells. Understanding the molecular evolution of proteins and the forces that shape protein evolution is essential in understanding the basis of function and evolution. Previous studies have shown that adaptation frequently occurs at the protein surface, such as in genes involved in host–pathogen interactions. However, it remains unclear whether adaptive sites are distributed randomly or at regions associated with particular structural or functional characteristics across the genome, since many proteins lack structural or functional annotations. Here, we seek to tackle this question by combining large-scale bioinformatic prediction, structural analysis, phylogenetic inference, and population genomic analysis of Drosophila protein-coding genes. We found that protein sequence adaptation is more relevant to function-related rather than structure-related properties. Interestingly, intermolecular interactions contribute significantly to protein adaptation. We further showed that intermolecular interactions, such as physical interactions, may play a role in the coadaptation of fast-adaptive proteins. We found that strongly differentiated amino acids across geographic regions in protein-coding genes are mostly adaptive, which may contribute to the long-term adaptive evolution. This strongly indicates that a number of adaptive sites tend to be repeatedly mutated and selected throughout evolution in the past, present, and maybe future. Our results highlight the important roles of intermolecular interactions and coadaptation in the adaptive evolution of proteins both at the species and population levels.     

Mutators, population size, adaptive landscape and the adaptation of asexual populations of bacteria.

Selection of mutator alleles, increasing the mutation rate up to 10, 000-fold, has been observed during in vitro experimental evolution. This spread is ascribed to the hitchhiking of mutator alleles with favorable mutations, as demonstrated by a theoretical model using selective parameters corresponding to such experiments. Observations of unexpectedly high frequencies of mutators in natural isolates suggest that the same phenomenon could occur in the wild. But it remains questionable whether realistic in natura parameter values could also result in selection of mutators. In particular, the main parameters of adaptation, the size of the adapting population and the height and steepness of the adaptive peak characterizing adaptation, are very variable in nature. By simulation approach, we studied the effect of these parameters on the selection of mutators in asexual populations, assuming additive fitness. We show that the larger the population size, the more likely the fixation of mutator alleles. At a large population size, at least four adaptive mutations are needed for mutator fixation; moreover, under stronger selection stronger mutators are selected. We propose a model based on multiple mutations to illustrate how second-order selection can optimize population fitness when few favorable mutations are required for adaptation.     

Ohno's model revisited: measuring the frequency of potentially adaptive mutations under various mutational drifts

The divergence of new gene functions is described by various scenarios that involve gene duplication, albeit, at fundamentally different stages. We performed experimental measurements and developed a subsequent model, aimed at predicting the rate of divergence under different scenarios. We used gene libraries of TEM-1 beta-lactamase that were drifted under purifying selection toward the original penicillinase activity or under no selection at all. The frequency of genes conferring a new function (degradation of a cephalosporin antibiotic) was measured at various stages of the drift, and a model that accounts for the differences in the observed adaptation dynamics of the drifting TEM-1 populations was derived. The results indicate that rapid nonfunctionalization in the population relieved from selection (Ohno's model) afforded only a narrow window of adaptation to cefotaxime (neofunctionalization). The trade-off between TEM-1's original function and the new evolving function also disfavored the "gene sharing" model. The rate of adaptation was maximal when selection for the original function was partially relieved to enable the accumulation of potentially adaptive mutations, while still purging a large fraction of otherwise deleterious mutations. Altogether, scenarios of subfunctionalization seem more feasible, whereby sustaining the original function by two copies facilitates the accumulation potentially adaptive mutations while purging nonfunctionalization mutations.     

Evolution of mutation rates in bacteria

Evolutionary success of bacteria relies on the constant fine-tuning of their mutation rates, which optimizes their adaptability to constantly changing environmental conditions. When adaptation is limited by the mutation supply rate, under some conditions, natural selection favours increased mutation rates by acting on allelic variation of the genetic systems that control fidelity of DNA replication and repair. Mutator alleles are carried to high frequency through hitchhiking with the adaptive mutations they generate. However, when fitness gain no longer counterbalances the fitness loss due to continuous generation of deleterious mutations, natural selection favours reduction of mutation rates. Selection and counter-selection of high mutation rates depends on many factors: the number of mutations required for adaptation, the strength of mutator alleles, bacterial population size, competition with other strains, migration, and spatial and temporal environmental heterogeneity. Such modulations of mutation rates may also play a role in the evolution of antibiotic resistance.     

Mitochondrial-nuclear interactions and accelerated compensatory evolution: evidence from the primate cytochrome C oxidase complex

Accelerated rates of mitochondrial protein evolution have been proposed to reflect Darwinian coadaptation for efficient energy production for mammalian flight and brain activity. However, several features of mammalian mtDNA (absence of recombination, small effective population size, and high mutation rate) promote genome degradation through the accumulation of weakly deleterious mutations. Here, we present evidence for "compensatory" adaptive substitutions in nuclear DNA- (nDNA) encoded mitochondrial proteins to prevent fitness decline in primate mitochondrial protein complexes. We show that high mutation rate and small effective population size, key features of primate mitochondrial genomes, can accelerate compensatory adaptive evolution in nDNA-encoded genes. We combine phylogenetic information and the 3D structure of the cytochrome c oxidase (COX) complex to test for accelerated compensatory changes among interacting sites. Physical interactions among mtDNA- and nDNA-encoded components are critical in COX evolution; amino acids in close physical proximity in the 3D structure show a strong tendency for correlated evolution among lineages. Only nuclear-encoded components of COX show evidence for positive selection and adaptive nDNA-encoded changes tend to follow mtDNA-encoded amino acid changes at nearby sites in the 3D structure. This bias in the temporal order of substitutions supports compensatory weak selection as a major factor in accelerated primate COX evolution.