FISHER'S MODEL AND THE GENOMICS OF ADAPTATION: RESTRICTED PLEIOTROPY,HETEROGENOUS MUTATION,AND PARALLEL EVOLUTION

Genetic theories of adaptation generally overlook the genes in which beneficial substitutions occur, and the likely variation in their mutational effects. We investigate the consequences of heterogeneous mutational effects among loci on the genetics of adaptation. We use a generalization of Fisher's geometrical model, which assumes multivariate Gaussian stabilizing selection on multiple characters. In our model, mutation has a distinct variance–covariance matrix of phenotypic effects for each locus. Consequently, the distribution of selection coefficients s varies across loci. We assume each locus can only affect a limited number of independent linear combinations of phenotypic traits (restricted pleiotropy), which differ among loci, an effect we term “orientation heterogeneity.” Restricted pleiotropy can sharply reduce the overall proportion of beneficial mutations. Orientation heterogeneity has little impact on the shape of the genomic distribution, but can substantially increase the probability of parallel evolution (the repeated fixation of beneficial mutations at the same gene in independent populations), which is highest with low pleiotropy. We also consider variation in the degree of pleiotropy and in the mean s across loci. The latter impacts the genomic distribution of s, but has a much milder effect on parallel evolution. We discuss these results in the light of evolution experiments.     

Rapid genome‐wide evolution in Brassica rapa populations following drought revealed by sequencing of ancestral and descendant gene pools

There is increasing evidence that evolution can occur rapidly in response to selection. Recent advances in sequencing suggest the possibility of documenting genetic changes as they occur in populations, thus uncovering the genetic basis of evolution, particularly if samples are available from both before and after selection. Here, we had a unique opportunity to directly assess genetic changes in natural populations following an evolutionary response to a fluctuation in climate. We analysed genome‐wide differences between ancestors and descendants of natural populations of Brassica rapa plants from two locations that rapidly evolved changes in multiple phenotypic traits, including flowering time, following a multiyear late‐season drought in California. These ancestor‐descendant comparisons revealed evolutionary shifts in allele frequencies in many genes. Some genes showing evolutionary shifts have functions related to drought stress and flowering time, consistent with an adaptive response to selection. Loci differentiated between ancestors and descendants (F_ST outliers) were generally different from those showing signatures of selection based on site frequency spectrum analysis (Tajima's D), indicating that the loci that evolved in response to the recent drought and those under historical selection were generally distinct. Very few genes showed similar evolutionary responses between two geographically distinct populations, suggesting independent genetic trajectories of evolution yielding parallel phenotypic changes. The results show that selection can result in rapid genome‐wide evolutionary shifts in allele frequencies in natural populations, and highlight the usefulness of combining resurrection experiments in natural populations with genomics for studying the genetic basis of adaptive evolution.     

Recurrent adaptation in a low‐dispersal trait

The study of natural populations from contrasting environments has greatly enhanced our understanding of ecological‐dependent selection, adaptation and speciation. Cases of parallel evolution in particular have facilitated the study of the molecular and genetic basis of adaptive variation. This includes the type and number of genes underlying adaptive traits, as well as the extent to which these genes are exchanged among populations and contribute repeatedly to parallel evolution. Yet, surprisingly few studies provide a comprehensive view on the evolutionary history of adaptive traits from mutation to widespread adaptation. When did key mutations arise, how did they increase in frequency, and how did they spread? In this issue of Molecular Ecology, Van Belleghem et al. ( 2015 ) reconstruct the evolutionary history of a gene associated with wing size in the salt marsh beetle Pogonus chalceus. Screening the entire distribution range of this species, they found a single origin for the allele associated with the short‐winged ecotype. This allele seemingly evolved in an isolated population and rapidly introgressed into other populations. These findings suggest that the adaptive genetic variation found in sympatric short‐ and long‐winged populations has an allopatric origin, confirming that allopatric phases may be important at early stages of speciation.     

AN ADAPTIVE WALK BY HUMAN IMMUNODEFICIENCY VIRUS TYPE 1 THROUGH A FLUCTUATING FITNESS LANDSCAPE

The mutational landscape model of adaptive sequence evolution has been used to explain an unexpected strong positive linear relationship between marginal fitness and mean site‐specific amino acid frequency in the functionally important HIV‐1 gp120 V3 protein region. The model predicts a positive linear relationship between the probability that a particular beneficial allele, among several, is the next to spread to fixation during an adaptive walk, its transition probability, and the allele's selection coefficient. Here, stochastic simulation is used to confirm the intuition that the linear relationship between transition probabilities and selection coefficients, predicted by the model, should, under fluctuating selection, produce a linear relationship between allele frequency, averaged across populations, and fitness. In addition, these relationships hold for the effective population size and mutation rate of HIV‐1 and for the moderately strong selection observed for V3. A survey of the strength of mutation for diverse organisms suggests that these relationships may be widely applicable.     

Heterosis in age‐specific selected populations of a seed beetle: Sex differences in longevity and reproductive behavior

We tested mutation accumulation hypothesis for the evolution of senescence using short‐lived and long‐lived populations of the seed‐feeding beetle, Acanthoscelides obtectus (Say), obtained by selection on early‐ and late‐life for many generations. The expected consequence of the mutation accumulation hypothesis is that in short‐lived populations, where the force of natural selection is the strongest early in life, the late‐life fitness traits should decline due to genetic drift which increases the frequency of mutations with deleterious effects in later adult stages. Since it is unlikely that identical deleterious mutations will increase in several independent populations, hybrid vigor for late‐life fitness is expected in offspring obtained in crosses among populations selected for early‐life fitness traits. We tested longevity of both sexes, female fecundity and male reproductive behavior for hybrid vigor by comparing hybrid and nonhybrid short‐lived populations. Hybrid vigor was confirmed for male virility, mating speed and copulation duration, and longevity of both sexes at late ages. In contrast to males, the results on female fecundity in short‐lived populations did not support mutation accumulation as a genetic mechanism for the evolution of this trait. Contrary to the prediction of this hypothesis, male mating ability indices and female fecundity in long‐lived populations exhibited hybrid vigor at all assayed age classes. We demonstrate that nonhybrid long‐lived populations diverged randomly regarding female and male reproductive fitness, indicating that sexually antagonistic selection, when accompanied with genetic drift for female fecundity and male virility, might be responsible for overriding natural selection in the independently evolving long‐lived populations.     

FURTHER SIMULATION STUDIES ON EVOLUTION BY GENE DUPLICATION

In order to understand the origin of multigene families, Monte Carlo simulations were performed to see how a genetic system evolves under unequal crossing-over, mutation, random genetic drift and natural selection, starting from a single gene copy. Both haploid and diploid models were examined. Beneficial, neutral, and detrimental mutations were incorporated, and “positive” selection favors those chromosomes (haploid) or individuals (diploid) with more beneficial mutations than others. The same model for haploids was previously investigated with special reference to the evolution of gene organization, and the ratio of the numbers of beneficial genes to pseudogenes was found to be a rough indicator of the relative strengths of positive and negative (against deleterious alleles) natural selection (Ohta, 1987b). In the present paper, the evolution of gene organization and of sequence divergence among genes in the multigene family is examined. It is shown that positive selection accelerates the accumulation of arrays containing different beneficial mutations, but that total divergence including both neutral and beneficial mutations is not very sensitive to positive selection, under this model. The proportion of beneficial mutations in the total mutations accumulated is a better indicator of positive selection than is the total divergence. It is pointed out that various observed examples in which amino-acid substitutions are accelerated, as compared with synonymous substitutions in duplicated genes (Li, 1985), may reflect the effect of selection similar to the present scheme. The diploid model is shown to be more efficient for accumulating beneficial mutations in duplicated genes than the haploid one, and the relevance of this finding to the advantage of sexual reproduction is discussed.     

Enhanced Synonymous Site Divergence in Positively Selected VertebrateAntimicrobial Peptide Genes

Nonrandom patterns associated with adaptively evolving genes can shed light on how selection and mutation produce rapid changes in sequences. I examine such patterns in two independent families of antimicrobial peptide genes: those in frogs, which are known to have evolved under positive selection, and those in flatfishes, which I show have also evolved under positive selection. I address two recently proposed hypotheses about the molecular evolution of antimicrobial peptide genes. The first is that the mature peptide region is replicated by an error-prone polymerase that increases the mutation rate and the transversion/transition ratio compared to the signal sequence of the same genes. The second is that mature peptides evolve in a coordinated fashion with their propieces, such that a change in net charge in one molecular region prompts an opposite change in charge in the other region. I test these hypotheses using alternative methods that minimize alignment errors, correct for phylogenetic nonindependence, reduce sequence saturation, and account for differing selection pressures on different regions of the gene. In both gene families I show that divergence at both synonymous and nonsynonymous sites within the mature peptide region is enhanced. However, in neither gene family is there evidence of an increased mutational transversion/transition ratio or coordinated evolution. My observations are consistent with either an elevated mutation rate in an adaptively evolving gene region or widespread selection on “silent” sites. These hypotheses challenge the assumption that mutations are random and can be measured by the synonymous substitution rate. [Reviewing Editor: Dr. Willie J. Swanson]     

The diversification of mate preferences by natural and sexual selection

The evolution of sexual display traits or preferences for them in response to divergent natural selection will alter sexual selection within populations, yet the role of sexual selection in ecological speciation has received little empirical attention. We evolved 12 populations of Drosophila serrata in a two‐way factorial design to investigate the roles of natural and sexual selection in the evolution of female mate preferences for male cuticular hydrocarbons (CHCs). Mate preferences weakened in populations evolving under natural selection alone, implying a cost in the absence of their expression. Comparison of the vectors of linear sexual selection revealed that the populations diverged in the combination of male CHCs that females found most attractive, although this was not significant using a mixed modelling approach. Changes in preference direction tended to evolve when natural and sexual selection were unconstrained, suggesting that both processes may be the key to initial stages of ecological speciation. Determining the generality of this result will require data from various species across a range of novel environments.     

RAPID EVOLUTION OF CHEATING MITOCHONDRIAL GENOMES IN SMALL YEAST POPULATIONS

Outcrossed sex exposes genes to competition with their homologues, allowing alleles that transmit more often than their competitors to spread despite organismal fitness costs. Mitochondrial populations in species with biparental inheritance are thought to be especially susceptible to such cheaters because they lack strict transmission rules like meiosis or maternal inheritance. Yet the interaction between mutation and natural selection in the evolution of cheating mitochondrial genomes has not been tested experimentally. Using yeast experimental populations, we show that although cheaters were rare in a large sample of spontaneous respiratory‐deficient mitochondrial mutations (petites), cheaters evolve under experimentally enforced outcrossing even when mutation supply and selection are restricted by repeatedly bottlenecking populations.     

PLEIOTROPY CAUSES LONG-TERM GENETIC CONSTRAINTS ON LIFE-HISTORY EVOLUTION IN BRASSICA RAPA

Fundamental, long-term genetic trade-offs constrain life-history evolution in wild crucifer populations. I studied patterns of genetic constraint in Brassica rapa by estimating genetic correlations among life-history components by quantitative genetic analyses among ten wild populations, and within four populations. Genetic correlations between age and size at first reproduction were always greater than +0.8 within and among all populations studied. Although quantitative genetics might provide insight about genetic constraints if genetic parameters remain approximately constant, little evidence has been available to determine the constancy of genetic correlations. I found strong and consistent estimates of genetic correlations between life-history components, which were very similar within four natural populations. Population differentiation also showed these same trade-offs, resulting from long-term genetic constraint. For some traits, evolutionary changes among populations were incompatible with a model of genetic drift. Historical patterns of natural selection were inferred from population differentiation, suggesting that correlated response to selection has caused some traits to evolve opposite to the direct forces of natural selection. Comparison with Arabidopsis suggests that these life-history trade-offs are caused by genes that regulate patterns of resource allocation to different components of fitness. Ecological and energetic models may correctly predict these trade-offs because there is little additive genetic variation for rates of resource acquisition, but resource allocation is genetically variable.     

EVALUATING PATTERNS OF CONVERGENT EVOLUTION AND TRANS‐SPECIES POLYMORPHISM AT MHC IMMUNOGENES IN TWO SYMPATRIC STICKLEBACK SPECIES

The immunologically important major histocompatibility complex (MHC) harbors some of the most polymorphic genes in vertebrates. These genes presumably evolve under parasite‐mediated selection and frequently show inconsistent allelic genealogies, where some alleles are more similar between species than within species. This phenomenon is thought to arise either from convergent evolution under parallel selection or from the preservation of ancient allelic lineages beyond speciation events (trans‐species polymorphism, TSP). Here, we examine natural populations of two sympatric stickleback species (Gasterosteus aculeatus and Pungitius pungitius) to investigate the contribution of these two mechanisms to the evolution of inconsistent allelic genealogies at the MHC. Overlapping parasite taxa between the two host species in three different habitats suggest contemporary parallel selection on the MHC genes. Accordingly, we detected a lack of species‐specific phylogenetic clustering in the immunologically relevant antigen‐binding residues of the MHC IIB genes which contrasted with the rest of the coding and noncoding sequence. However, clustering was not habitat‐specific and a codon‐usage analysis revealed patterns of similarity by descent. In this light, common descent via TSP, in combination with intraspecies gene conversion, rather than convergent evolution is the more strongly supported scenario for the inconsistent genealogy at the MHC.     

Genomic evidence for parallel adaptation to cities

Urban evolutionary biology is the study of rapid evolutionary change in response to humans and our uses of land to support city dwellers. Because cities are relatively modern additions to the natural world, research on urban evolution tends to focus on microevolutionary change that has happened across a few to many hundreds of generations. These questions still fall under the broad purview of evolutionary ecology. However, the severity, rapidity and replication of environmental changes that drive evolution in this context make it worthy of specific attention. Urban evolution provides the opportunity to study the earliest stages of evolution in a context that is scientifically interesting and societally important. The newness of urban populations and their proximity to natural populations also creates challenges when trying to detect population genetic change. In a From the Cover article in this issue of Molecular Ecology, Mueller et al. use whole genome resequencing data to address some of these challenges while exploring genetic changes associated with urbanization in three replicate urban‐rural burrowing owl (Athene cunicularia) populations. Combining multiple approaches across these sample sites Mueller et al. find evidence for selection on genes whose function is related to synapses, neuron projections, brain connectivity and cognitive function in general. That selection was parallel suggests that phenotypes related to brain processes were probably particularly important for urban adaptation.     

Negative frequency‐dependent selection maintains coexisting genotypes during fluctuating selection

Natural environments are rarely static; rather selection can fluctuate on timescales ranging from hours to centuries. However, it is unclear how adaptation to fluctuating environments differs from adaptation to constant environments at the genetic level. For bacteria, one key axis of environmental variation is selection for planktonic or biofilm modes of growth. We conducted an evolution experiment with Burkholderia cenocepacia, comparing the evolutionary dynamics of populations evolving under constant selection for either biofilm formation or planktonic growth with populations in which selection fluctuated between the two environments on a weekly basis. Populations evolved in the fluctuating environment shared many of the same genetic targets of selection as those evolved in constant biofilm selection, but were genetically distinct from the constant planktonic populations. In the fluctuating environment, mutations in the biofilm‐regulating genes wspA and rpfR rose to high frequency in all replicate populations. A mutation in wspA first rose rapidly and nearly fixed during the initial biofilm phase but was subsequently displaced by a collection of rpfR mutants upon the shift to the planktonic phase. The wspA and rpfR genotypes coexisted via negative frequency‐dependent selection around an equilibrium frequency that shifted between the environments. The maintenance of coexisting genotypes in the fluctuating environment was unexpected. Under temporally fluctuating environments, coexistence of two genotypes is only predicted under a narrow range of conditions, but the frequency‐dependent interactions we observed provide a mechanism that can increase the likelihood of coexistence in fluctuating environments.     

Mutational neighbourhood and mutation supply rate constrain adaptation in Pseudomonas aeruginosa

Understanding adaptation by natural selection requires understanding the genetic factors that determine which beneficial mutations are available for selection. Here, using experimental evolution of rifampicin-resistant Pseudomonas aeruginosa, we show that different genotypes vary in their capacity for adaptation to the cost of antibiotic resistance. We then use sequence data to show that the beneficial mutations associated with fitness recovery were specific to particular genetic backgrounds, suggesting that genotypes had access to different sets of beneficial mutations. When we manipulated the supply rate of beneficial mutations, by altering effective population size during evolution, we found that it constrained adaptation in some selection lines by restricting access to rare beneficial mutations, but that the effect varied among the genotypes in our experiment. These results suggest that mutational neighbourhood varies even among genotypes that differ by a single amino acid change, and this determines their capacity for adaptation as well as the influence of population biology processes that alter mutation supply rate.     

The origins of weedy rice

Where do weeds come from? How do they evolve from nonweedy ancestors? In this issue of Molecular Ecology, Londo and Schaal examine the origin of weedy rice (Oryza sativa) populations in the USA. Analysing nuclear DNA sequence and microsatellite data, they show the importance of parallel evolution, hybridization, gene flow, and migration in the evolution of these weeds.     

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.     

Evidence of gene orthology and trans‐species polymorphism,but not of parallel evolution,despite high levels of concerted evolution in the major histocompatibility complex of flamingo species

The major histocompatibility complex (MHC) is a cornerstone in the study of adaptive genetic diversity. Intriguingly, highly polymorphic MHC sequences are often not more similar within species than between closely related species. Divergent selection of gene duplicates, balancing selection maintaining trans‐species polymorphism (TSP) that predate speciation and parallel evolution of species sharing similar selection pressures can all lead to higher sequence similarity between species. In contrast, high rates of concerted evolution increase sequence similarity of duplicated loci within species. Assessing these evolutionary models remains difficult as relatedness and ecological similarities are often confounded. As sympatric species of flamingos are more distantly related than allopatric species, flamingos represent an ideal model to disentangle these evolutionary models. We characterized MHC Class I exon 3, Class IIB exon 2 and exon 3 of the six extant flamingo species. We found up to six MHC Class I loci and two MHC Class IIB loci. As all six species shared the same number of MHC Class IIB loci, duplication appears to predate flamingo speciation. However, the high rate of concerted evolution has prevented the divergence of duplicated loci. We found high sequence similarity between all species regardless of codon position. The latter is consistent with balancing selection maintaining TSP, as under this mechanism amino acid sites under pathogen‐mediated selection should be characterized by fewer synonymous codons (due to their common ancestry) than under parallel evolution. Overall, balancing selection maintaining TSP appears to result in high MHC similarity between species regardless of species relatedness and geographical distribution.     

The Fixation Probability of Rare Mutators in Finite Asexual Populations

A mutator is an allele that increases the mutation rate throughout the genome by disrupting some aspect of DNA replication or repair. Mutators that increase the mutation rate by the order of 100-fold have been observed to spontaneously emerge and achieve high frequencies in natural populations and in long-term laboratory evolution experiments with Escherichia coli. In principle, the fixation of mutator alleles is limited by (i) competition with mutations in wild-type backgrounds, (ii) additional deleterious mutational load, and (iii) random genetic drift. Using a multiple-locus model and employing both simulation and analytic methods, we investigate the effects of these three factors on the fixation probability P_fix of an initially rare mutator as a function of population size N, beneficial and deleterious mutation rates, and the strength of mutations s. Our diffusion-based approximation for P_fix successfully captures effects ii and iii when selection is fast compared to mutation (). This enables us to predict the conditions under which mutators will be evolutionarily favored. Surprisingly, our simulations show that effect i is typically small for strong-effect mutators. Our results agree semiquantitatively with existing laboratory evolution experiments and suggest future experimental directions.     

The effect of population structure on the rate of evolution

Ecological factors exert a range of effects on the dynamics of the evolutionary process. A particularly marked effect comes from population structure, which can affect the probability that new mutations reach fixation. Our interest is in population structures, such as those depicted by ‘star graphs’, that amplify the effects of selection by further increasing the fixation probability of advantageous mutants and decreasing the fixation probability of disadvantageous mutants. The fact that star graphs increase the fixation probability of beneficial mutations has lead to the conclusion that evolution proceeds more rapidly in star-structured populations, compared with mixed (unstructured) populations. Here, we show that the effects of population structure on the rate of evolution are more complex and subtle than previously recognized and draw attention to the importance of fixation time. By comparing population structures that amplify selection with other population structures, both analytically and numerically, we show that evolution can slow down substantially even in populations where selection is amplified.