What paths do advantageous alleles take during short‐term evolutionary change?

Combining experimental evolution with whole‐genome resequencing is a promising new strategy for investigating the dynamics of evolutionary change. Published studies that have resequenced laboratory‐selected populations of sexual organisms have typically focused on populations sampled at the end of an evolution experiment. These studies have attempted to associate particular alleles with phenotypic change and attempted to distinguish between different theoretical models of adaptation. However, neither the population used to initiate the experiment nor multiple time points sampled during the evolutionary trajectory are generally available for examination. In this issue of Molecular Ecology, Orozco‐terWengel et al. (2012) take a significant step forward by estimating genome‐wide allele frequencies at the start, 15 generations into and at the end of a 37‐generation Drosophila experimental evolution study. The authors identify regions of the genome that have responded to laboratory selection and describe the temporal dynamics of allele frequency change. They identify two common trajectories for putatively adaptive alleles: alleles either gradually increase in frequency throughout the entire 37 generations or alleles plateau at a new frequency by generation 15. The identification of complex trajectories of alleles under selection contributes to a growing body of literature suggesting that simple models of adaptation, whereby beneficial alleles arise and increase in frequency unimpeded until they become fixed, may not adequately describe short‐term response to selection.     

EVIDENCE OF ADAPTATION FROM ANCESTRAL VARIATION IN YOUNG POPULATIONS OF BEACH MICE

To understand how organisms adapt to novel habitats, which involves both demographic and selective events, we require knowledge of the evolutionary history of populations and also selected alleles. There are still few cases in which the precise mutations (and hence, defined alleles) that contribute to adaptive change have been identified in nature; one exception is the genetic basis of camouflaging pigmentation of oldfield mice (Peromyscus polionotus) that have colonized the sandy dunes of Florida's Gulf Coast. To quantify the genomic impact of colonization as well as the signature of selection, we resequenced 5000 1.5‐kb noncoding loci as well as a 160‐kb genomic region surrounding the melanocortin‐1 receptor (Mc1r), a gene that contributes to pigmentation differences, in beach and mainland populations. Using a genome‐wide phylogenetic approach, we recovered a single monophyletic group comprised of beach mice, consistent with a single colonization event of the Gulf Coast. We also found evidence of a severe founder event, estimated to have occurred less than 3000 years ago. In this demographic context, we show that all beach subspecies share a single derived light Mc1r allele, which was likely selected from standing genetic variation that originated in the mainland. Surprisingly, we were unable to identify a clear signature of selection in the Mc1r region, despite independent evidence that this locus contributes to adaptive coloration. Nonetheless, these data allow us to reconstruct and compare the evolutionary history of populations and alleles to better understand how adaptive evolution, following the colonization of a novel habitat, proceeds in nature.     

Genome‐wide selection components analysis in a fish with male pregnancy

A major goal of evolutionary biology is to identify the genome‐level targets of natural and sexual selection. With the advent of next‐generation sequencing, whole‐genome selection components analysis provides a promising avenue in the search for loci affected by selection in nature. Here, we implement a genome‐wide selection components analysis in the sex role reversed Gulf pipefish, Syngnathus scovelli. Our approach involves a double‐digest restriction‐site associated DNA sequencing (ddRAD‐seq) technique, applied to adult females, nonpregnant males, pregnant males, and their offspring. An F_ST comparison of allele frequencies among these groups reveals 47 genomic regions putatively experiencing sexual selection, as well as 468 regions showing a signature of differential viability selection between males and females. A complementary likelihood ratio test identifies similar patterns in the data as the F_ST analysis. Sexual selection and viability selection both tend to favor the rare alleles in the population. Ultimately, we conclude that genome‐wide selection components analysis can be a useful tool to complement other approaches in the effort to pinpoint genome‐level targets of selection in the wild.     

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.     

Eco‐genomic analysis of the poleward range expansion of the wasp spider Argiope bruennichi shows rapid adaptation and genomic admixture

Poleward range expansions are commonly attributed to global change, but could alternatively be driven by rapid evolutionary adaptation. A well‐documented example of a range expansion during the past decades is provided by the European wasp spider Argiope bruennichi. Using ecological niche modeling, thermal tolerance experiments and a genome‐wide analysis of gene expression divergence, we show that invasive populations have adapted to novel climatic conditions in the course of their expansion. Their climatic niche shift is mirrored in an increased cold tolerance and a population‐specific and functionally differentiated gene expression response. We generated an Argiope reference genome sequence and used population genome resequencing to assess genomic changes associated with the new climatic adaptations. We find clear genetic differentiation and a significant admixture with alleles from East Asian populations in the invasive Northern European populations. Population genetic modeling suggests that at least some of these introgressing alleles have contributed to the new adaptations during the expansion. Our results thus confirm the notion that range expansions are not a simple consequence of climate change, but are accompanied by fast genetic changes and adaptations that may be fuelled through admixture between long separated lineages.     

Testing for a genetic response to sexual selection in a wild Drosophila population

In accordance with the consensus that sexual selection is responsible for the rapid evolution of display traits on macroevolutionary scales, microevolutionary studies suggest sexual selection is a widespread and often strong form of directional selection in nature. However, empirical evidence for the contemporary evolution of sexually selected traits via sexual rather than natural selection remains weak. In this study, we used a novel application of quantitative genetic breeding designs to test for a genetic response to sexual selection on eight chemical display traits from a field population of the fly, Drosophila serrata. Using our quantitative genetic approach, we were able to detect a genetically based difference in means between groups of males descended from fathers who had either successfully sired offspring or were randomly collected from the same wild population for one of these display traits, the diene (Z,Z)‐5,9‐C_27 : 2. Our experimental results, in combination with previous laboratory studies on this system, suggest that both natural and sexual selection may be influencing the evolutionary trajectories of these traits in nature, limiting the capacity for a contemporary evolutionary response.     

Sexual conflict and the maintenance of genetic variation in natural populations

Understanding the factors that maintain genetic variation in natural populations is a foundational goal of evolutionary biology. To this end, population geneticists have developed a variety of models that can produce stable polymorphisms. In one of the earliest models, Owen ( 1953 ) demonstrated that differences in selection pressures acting on males and females could maintain multiple alleles of a gene at a stable equilibrium. If the selection pressures act in opposite directions in males and females, we refer to this as (inter‐) sexual conflict or sexual antagonism (Arnqvist & Rowe, 2005 ). Testing if sexual conflict maintains genetic variation in natural populations is a tremendous challenge—it requires both identifying loci that harbor sexually antagonistic alleles and determining whether those alleles are maintained as stable polymorphisms (Mank, 2017 ). Doing so genome‐wide is even harder because it is not tractable to identify sexually antagonistic alleles and test for stable polymorphisms at all loci. Dutoit et al. ( 2018 ) confront this challenge in a paper published in this issue of Molecular Ecology. Using gene expression and population genomic data from the collared flycatcher, Dutoit et al. ( 2018 ) identify associations and correlations between genomic signatures of balanced polymorphisms and sexual conflict.     

A genomic selection component analysis characterizes migration‐selection balance

The genetic differentiation of populations in response to local selection pressures has long been studied by evolutionary biologists, but key details about the process remain obscure. How rapidly can local adaptation evolve, how extensive is the process across the genome, and how strong are the opposing forces of natural selection and gene flow? Here, we combine direct measurement of survival and reproduction with whole‐genome genotyping of a plant species (Mimulus guttatus) that has recently invaded a novel habitat (the Quarry population). We renovate the classic selection component method to accommodate genomic data and observe selection at SNPs throughout the genome. SNPs showing viability selection in Quarry exhibit elevated divergence from neighboring populations relative to neutral SNPs. We also find that nonsignificant SNPs exhibit a subtle, but still significant, change in allele frequency toward neighboring populations, a predicted effect of gene flow. Given that the Quarry population is most probably only 30–40 generations old, the alleles conferring local advantage are almost certainly older than the population itself. Thus, local adaptation owes to the recruitment of standing genetic variation.     

Evidence for widespread selection in shaping the genomic landscape during speciation of Populus

Increasing our understanding of how evolutionary processes drive the genomic landscape of variation is fundamental to a better understanding of the genomic consequences of speciation. However, genome‐wide patterns of within‐ and between‐ species variation have not been fully investigated in most forest tree species despite their global ecological and economic importance. Here, we use whole‐genome resequencing data from four Populus species spanning the speciation continuum to reconstruct their demographic histories and investigate patterns of diversity and divergence within and between species. Using Populus trichocarpa as an outgroup species, we further infer the genealogical relationships and estimate the extent of ancient introgression among the three aspen species (Populus tremula, Populus davidiana and Populus tremuloides) throughout the genome. Our results show substantial variation in these patterns along the genomes with this variation being strongly predicted by local recombination rates and the density of functional elements. This implies that the interaction between recurrent selection and intrinsic genomic features has dramatically sculpted the genomic landscape over long periods of time. In addition, our findings provide evidence that, apart from background selection, recent positive selection and long‐term balancing selection have also been crucial components in shaping patterns of genome‐wide variation during the speciation process.     

GENETIC HITCHHIKING AND THE DYNAMIC BUILDUP OF GENOMIC DIVERGENCE DURING SPECIATION WITH GENE FLOW

A major issue in evolutionary biology is explaining patterns of differentiation observed in population genomic data, as divergence can be due to both direct selection on a locus and genetic hitchhiking. “Divergence hitchhiking” (DH) theory postulates that divergent selection on a locus reduces gene flow at physically linked sites, facilitating the formation of localized clusters of tightly linked, diverged loci. “Genome hitchhiking” (GH) theory emphasizes genome‐wide effects of divergent selection. Past theoretical investigations of DH and GH focused on static snapshots of divergence. Here, we used simulations assessing a variety of strengths of selection, migration rates, population sizes, and mutation rates to investigate the relative importance of direct selection, GH, and DH in facilitating the dynamic buildup of genomic divergence as speciation proceeds through time. When divergently selected mutations were limiting, GH promoted divergence, but DH had little measurable effect. When populations were small and divergently selected mutations were common, DH enhanced the accumulation of weakly selected mutations, but this contributed little to reproductive isolation. In general, GH promoted reproductive isolation by reducing effective migration rates below that due to direct selection alone, and was important for genome‐wide “congealing” or “coupling” of differentiation (F_ST) across loci as speciation progressed.     

Screening of duplicated loci reveals hidden divergence patterns in a complex salmonid genome

A whole‐genome duplication (WGD) doubles the entire genomic content of a species and is thought to have catalysed adaptive radiation in some polyploid‐origin lineages. However, little is known about general consequences of a WGD because gene duplicates (i.e., paralogs) are commonly filtered in genomic studies; such filtering may remove substantial portions of the genome in data sets from polyploid‐origin species. We demonstrate a new method that enables genome‐wide scans for signatures of selection at both nonduplicated and duplicated loci by taking locus‐specific copy number into account. We apply this method to RAD sequence data from different ecotypes of a polyploid‐origin salmonid (Oncorhynchus nerka) and reveal signatures of divergent selection that would have been missed if duplicated loci were filtered. We also find conserved signatures of elevated divergence at pairs of homeologous chromosomes with residual tetrasomic inheritance, suggesting that joint evolution of some nondiverged gene duplicates may affect the adaptive potential of these genes. These findings illustrate that including duplicated loci in genomic analyses enables novel insights into the evolutionary consequences of WGDs and local segmental gene duplications.     

Genomics of invasion: diversity and selection in introduced populations of monkeyflowers (Mimulus guttatus)

Global trade and travel is irreversibly changing the distribution of species around the world. Because introduced species experience drastic demographic events during colonization and often face novel environmental challenges from their native range, introduced populations may undergo rapid evolutionary change. Genomic studies provide the opportunity to investigate the extent to which demographic, historical and selective processes shape the genomic structure of introduced populations by analysing the signature that these processes leave on genomic variation. Here, we use next‐generation sequencing to compare genome‐wide relationships and patterns of diversity in native and introduced populations of the yellow monkeyflower (Mimulus guttatus). Genome resequencing data from 10 introduced populations from the United Kingdom (UK) and 12 native M. guttatus populations in North America (NA) demonstrated reduced neutral genetic diversity in the introduced range and showed that UK populations are derived from a geographic region around the North Pacific. A selective‐sweep analysis revealed site frequency changes consistent with selection on five of 14 chromosomes, with genes in these regions showing reduced silent site diversity. While the target of selection is unknown, genes associated with flowering time and biotic and abiotic stresses were located within the swept regions. The future identification of the specific source of origin of introduced UK populations will help determining whether the observed selective sweeps can be traced to unsampled native populations or occurred since dispersal across the Atlantic. Our study demonstrates the general potential of genome‐wide analyses to uncover a range of evolutionary processes affecting invasive populations.     

Experimental evidence of genome‐wide impact of ecological selection during early stages of speciation‐with‐gene‐flow

Theory predicts that speciation‐with‐gene‐flow is more likely when the consequences of selection for population divergence transitions from mainly direct effects of selection acting on individual genes to a collective property of all selected genes in the genome. Thus, understanding the direct impacts of ecologically based selection, as well as the indirect effects due to correlations among loci, is critical to understanding speciation. Here, we measure the genome‐wide impacts of host‐associated selection between hawthorn and apple host races of Rhagoletis pomonella (Diptera: Tephritidae), a model for contemporary speciation‐with‐gene‐flow. Allele frequency shifts of 32 455 SNPs induced in a selection experiment based on host phenology were genome wide and highly concordant with genetic divergence between co‐occurring apple and hawthorn flies in nature. This striking genome‐wide similarity between experimental and natural populations of R. pomonella underscores the importance of ecological selection at early stages of divergence and calls for further integration of studies of eco‐evolutionary dynamics and genome divergence.     

Isolation and characterization of microsatellite loci from a commercial cultivar of Musa acuminata

A genomic library from the commercial diploid cultivar ‘Ouro’ (Musa acuminata), enriched for CT‐ and GT‐repeats, was used to isolate and characterize 23 microsatellite loci. These loci were tested in 10 Musa genotypes, representing various Musa genomic groups with distinct ploidy level. The number of alleles per locus ranged from one to seven, and 20 loci were highly informative. Four loci appeared to amplify B genome‐specific alleles, while three loci seemed to be absent in the B genome. The polymorphism revealed by these loci will be extremely useful for genetic mapping, marker‐assisted selection, germplasm characterization and evolutionary studies in Musa.     

Genome‐wide association analyses reveal polygenic genomic architecture underlying divergent shell morphology in Spanish Littorina saxatilis ecotypes

Gene flow between diverging populations experiencing dissimilar ecological conditions can theoretically constrain adaptive evolution. To minimize the effect of gene flow, alleles underlying traits essential for local adaptation are predicted to be located in linked genome regions with reduced recombination. Local reduction in gene flow caused by selection is expected to produce elevated divergence in these regions. The highly divergent crab‐adapted and wave‐adapted ecotypes of the marine snail Littorina saxatilis present a model system to test these predictions. We used genome‐wide association (GWA) analysis of geometric morphometric shell traits associated with microgeographic divergence between the two L. saxatilis ecotypes within three separate sampling sites. A total of 477 snails that had individual geometric morphometric data and individual genotypes at 4,066 single nucleotide polymorphisms (SNPs) were analyzed using GWA methods that corrected for population structure among the three sites. This approach allowed dissection of the genomic architecture of shell shape divergence between ecotypes across a wide geographic range, spanning two glacial lineages. GWA revealed 216 quantitative trait loci (QTL) with shell size or shape differences between ecotypes, with most loci explaining a small proportion of phenotypic variation. We found that QTL were evenly distributed across 17 linkage groups, and exhibited elevated interchromosomal linkage, suggesting a genome‐wide response to divergent selection on shell shape between the two ecotypes. Shell shape trait‐associated loci showed partial overlap with previously identified outlier loci under divergent selection between the two ecotypes, supporting the hypothesis of diversifying selection on these genomic regions. These results suggest that divergence in shell shape between the crab‐adapted and wave‐adapted ecotypes is produced predominantly by a polygenic genomic architecture with positive linkage disequilibrium among loci of small effect.     

THE GENOMIC TRAJECTORY OF HYBRID SWARMS: OUTCOMES OF REPEATED CROSSES BETWEEN POPULATIONS OF TIGRIOPUS CALIFORNICUS

Introgressive hybridization between genetically divergent populations is an important evolutionary process. The degree to which repeated hybridization events between the same parental taxa lead to similar genomic outcomes is unknown. This study addressed this question by following genomic trajectories of replicate hybrid swarms of the copepod Tigriopus californicus over many generations of free mating. Swarm composition was determined both by differential reproductive success of founder individuals and subsequent selection on hybrid genotypes. For one cross, between two populations showing differential fitness in the laboratory and no hybrid breakdown, the genetic trajectory was highly repeatable: replicates rapidly became dominated by alleles from the fitter parent. In a second cross, between two populations showing similar fitness and significant F2 hybrid breakdown, alleles from alternative populations dominated different replicates. Swarms exhibited a general temporal trend of decreasing cytonuclear mismatch. Some patterns of differential introgression across the genome were strikingly congruent amongst swarm replicates, both within and between cross types, and reflected patterns of segregation distortion previously observed within controlled crosses between the same parental populations. Differences in heterozygosity between the sexes, and evidence for a previously suspected sex‐distortion locus, suggest that complex interactions between sex and genotype influence hybrid swarm outcome.     

Model‐based demographic inference of introgression history in European whitefish species pairs'

Parallel phenotypic differentiation is generally attributed to parallel adaptive divergence as an evolutionary response to similar environmental contrasts. Such parallelism may actually originate from several evolutionary scenarios ranging from repeated parallel divergence caused by divergent selection to a unique divergence event followed by gene flow. Reconstructing the evolutionary history underlying parallel phenotypic differentiation is thus fundamental to understand the relative contribution of demography and selection on genomic divergence during speciation. In this study, we investigate the divergence history of replicate European whitefish (Coregonus lavaretus), limnetic and benthic species pairs from two lakes in Norway and two lakes in Switzerland. Demographic models accounting for semi‐permeability and linked selection were fitted to the unfolded joint allele frequency spectrum built from genome‐wide SNPs and compared to each other in each species pair. We found strong support for a model of asymmetrical post‐glacial secondary contact between glacial lineages in all four lakes. Moreover, our results suggest that heterogeneous genomic differentiation has been shaped by the joint action of linked selection accelerating lineage sorting during allopatry, and heterogeneous migration eroding divergence at different rates along the genome following secondary contact. Our analyses reveal how the interplay between demography, selection and historical contingency has influenced the levels of diversity observed in previous whitefish phylogeographic studies. This study thus provides new insights into the historical demographic and selective processes that shaped the divergence associated with ecological speciation in European whitefish.     

Genetic changes in a novel breeding population of Brassica napus synthesized from hundreds of crosses between B. rapa and B. carinata

Introgression of genomic variation between and within related crop species is a significant evolutionary approach for population differentiation, genome reorganization and trait improvement. Using the Illumina Infinium Brassica 60K SNP array, we investigated genomic changes in a panel of advanced generation new‐type Brassica napus breeding lines developed from hundreds of interspecific crosses between 122 Brassica rapa and 74 Brassica carinata accessions, and compared them with representative accessions of their three parental species. The new‐type B. napus population presented rich genetic diversity and abundant novel genomic alterations, consisting of introgressions from B. rapa and B. carinata, novel allelic combinations, reconstructed linkage disequilibrium patterns and haplotype blocks, and frequent deletions and duplications (nonrandomly distributed), particularly in the C subgenome. After a much shorter, but very intensive, selection history compared to traditional B. napus, a total of 15 genomic regions with strong selective sweeps and 112 genomic regions with putative signals of selective sweeps were identified. Some of these regions were associated with important agronomic traits that were selected for during the breeding process, while others were potentially associated with restoration of genome stability and fertility after interspecific hybridization. Our results demonstrate how a novel method for population‐based crop genetic improvement can lead to rapid adaptation, restoration of genome stability and positive responses to artificial selection.     

Genome‐wide association and genome partitioning reveal novel genomic regions underlying variation in gastrointestinal nematode burden in a wild bird

Identifying the genetic architecture underlying complex phenotypes is a notoriously difficult problem that often impedes progress in understanding adaptive eco‐evolutionary processes in natural populations. Host–parasite interactions are fundamentally important drivers of evolutionary processes, but a lack of understanding of the genes involved in the host's response to chronic parasite insult makes it particularly difficult to understand the mechanisms of host life history trade‐offs and the adaptive dynamics involved. Here, we examine the genetic basis of gastrointestinal nematode (Trichostrongylus tenuis) burden in 695 red grouse (Lagopus lagopus scotica) individuals genotyped at 384 genome‐wide SNPs. We first use genome‐wide association to identify individual SNPs associated with nematode burden. We then partition genome‐wide heritability to identify chromosomes with greater heritability than expected from gene content, due to harbouring a multitude of additive SNPs with individually undetectable effects. We identified five SNPs on five chromosomes that accounted for differences of up to 556 worms per bird, but together explained at best 4.9% of the phenotypic variance. These SNPs were closely linked to genes representing a range of physiological processes including the immune system, protein degradation and energy metabolism. Genome partitioning indicated genome‐wide heritability of up to 29% and three chromosomes with excess heritability of up to 4.3% (total 8.9%). These results implicate SNPs and novel genomic regions underlying nematode burden in this system and suggest that this phenotype is somewhere between being based on few large‐effect genes (oligogenic) and based on a large number of genes with small individual but large combined effects (polygenic).