Integrating QTL mapping and genome scans towards the characterization of candidate loci under parallel selection in the lake whitefish (Coregonus clupeaformis)

As natural selection must act on underlying genetic variation, discovering the number and location of loci under the influence of selection is imperative towards understanding adaptive divergence in evolving populations. Studies employing genome scans have hypothesized that the action of divergent selection should reduce gene flow at the genomic locations implicated in adaptation and speciation among natural populations, yet once 'outlier' patterns of variation have been identified the function and role of such loci needs to be confirmed. We integrated adaptive QTL mapping and genomic scans among diverging sympatric pairs of the lake whitefish (Coregonus clupeaformis) species complex in order to test the hypothesis that differentiation between dwarf and normal ecotypes at growth-associated QTL was maintained by directional selection. We found evidence of significantly high levels of molecular divergence among eight growth QTL where two of the strongest candidate loci under the influence of directional selection exhibited parallel reductions of gene flow over multiple populations.     

Advances in Ecological Speciation: an integrative approach

The role of natural selection in promoting reproductive isolation has received substantial renewed interest within the last two decades. As a consequence, the study of ecological speciation has become an extremely productive research area in modern evolutionary biology. Recent innovations in sequencing technologies offer an unprecedented opportunity to study the mechanisms involved in ecological speciation. Genome scans provide significant insights but have some important limitations; efforts are needed to integrate them with other approaches to make full use of the sequencing data deluge. An international conference ‘Advances in Ecological Speciation’ organized by the University of Porto (Portugal) aimed to review current progress in ecological speciation. Using some of the examples presented at the conference, we highlight the benefits of integrating ecological and genomic data and discuss different mechanisms of parallel evolution. Finally, future avenues of research are suggested to advance our knowledge concerning the role of natural selection in the establishment of reproductive isolation during ecological speciation.     

Divergence hitchhiking and the spread of genomic isolation during ecological speciation-with-gene-flow

In allopatric populations, geographical separation simultaneously isolates the entire genome, allowing genetic divergence to accumulate virtually anywhere in the genome. In sympatric populations, however, the strong divergent selection required to overcome migration produces a genetic mosaic of divergent and non-divergent genomic regions. In some recent genome scans, each divergent genomic region has been interpreted as an independent incidence of migration/selection balance, such that the reduction of gene exchange is restricted to a few kilobases around each divergently selected gene. I propose an alternative mechanism, 'divergence hitchhiking' (DH), in which divergent selection can reduce gene exchange for several megabases around a gene under strong divergent selection. Not all genes/markers within a DH region are divergently selected, yet the entire region is protected to some degree from gene exchange, permitting genetic divergence from mechanisms other than divergent selection to accumulate secondarily. After contrasting DH and multilocus migration/selection balance (MM/SB), I outline a model in which genomic isolation at a given genomic location is jointly determined by DH and genome-wide effects of the progressive reduction in realized migration, then illustrate DH using data from several pairs of incipient species in the wild.     

Genomic signatures of divergent selection and speciation patterns in a ‘natural experiment’, the young parallel radiations of Nicaraguan crater lake cichlid fishes

Divergent selection is the main driving force in sympatric ecological speciation and may also play a strong role in divergence between allopatric populations. Characterizing the genome‐wide impact of divergent selection often constitutes a first step in unravelling the genetic bases underlying adaptation and ecological speciation. The Midas cichlid fish (Amphilophus citrinellus) species complex in Nicaragua is a powerful system for studying evolutionary processes. Independent colonizations of isolated young crater lakes by Midas cichlid populations from the older and great lakes of Nicaragua resulted in the repeated evolution of adaptive radiations by intralacustrine sympatric speciation. In this study we performed genome scans on two repeated radiations of crater lake species and their great lake source populations (1030 polymorphic AFLPs, n ～ 30 individuals per species). We detected regions under divergent selection (0.3% in the crater lake Xiloá flock and 1.7% in the older crater lake Apoyo radiation) that might be responsible for the sympatric diversifications. We find no evidence that the same genomic regions have been involved in the repeated evolution of parallel adaptations across crater lake flocks. However, there is some genetic parallelism apparent (seven out of 51 crater lake to great lake outlier loci are shared; 13.7%) that is associated with the allopatric divergence of both crater lake flocks. Interestingly, our results suggest that the number of outlier loci involved in sympatric and allopatric divergence increases over time. A phylogeny based on the AFLP data clearly supports the monophyly of both crater lake species flocks and indicates a parallel branching order with a primary split along the limnetic‐benthic axis in both radiations.     

Global analysis of genes involved in freshwater adaptation in threespine sticklebacks (Gasterosteus aculeatus)

Examples of parallel evolution of phenotypic traits have been repeatedly demonstrated in threespine sticklebacks (Gasterosteus aculeatus) across their global distribution. Using these as a model, we performed a targeted genome scan--focusing on physiologically important genes potentially related to freshwater adaptation--to identify genetic signatures of parallel physiological evolution on a global scale. To this end, 50 microsatellite loci, including 26 loci within or close to (<6 kb) physiologically important genes, were screened in paired marine and freshwater populations from six locations across the Northern Hemisphere. Signatures of directional selection were detected in 24 loci, including 17 physiologically important genes, in at least one location. Although no loci showed consistent signatures of selection in all divergent population pairs, several outliers were common in multiple locations. In particular, seven physiologically important genes, as well as reference ectodysplasin gene (EDA), showed signatures of selection in three or more locations. Hence, although these results give some evidence for consistent parallel molecular evolution in response to freshwater colonization, they suggest that different evolutionary pathways may underlie physiological adaptation to freshwater habitats within the global distribution of the threespine stickleback.     

Parallel and nonparallel aspects of ecological, phenotypic, and genetic divergence across replicate population pairs of lake and stream stickleback

Parallel (or convergent) evolution provides strong evidence for a deterministic role of natural selection: similar phenotypes evolve when independent populations colonize similar environments. In reality, however, independent populations in similar environments always show some differences: some nonparallel evolution is present. It is therefore important to explicitly quantify the parallel and nonparallel aspects of trait variation, and to investigate the ecological and genetic explanations for each. We performed such an analysis for threespine stickleback (Gasterosteus aculeatus) populations inhabiting lake and stream habitats in six independent watersheds. Morphological traits differed in the degree to which lake-stream divergence was parallel across watersheds. Some aspects of this variation were correlated with ecological variables related to diet, presumably reflecting the strength and specifics of divergent selection. Furthermore, a genetic scan revealed some markers that diverged between lakes and streams in many of the watersheds and some that diverged in only a few watersheds. Moreover, some of the lake-stream divergence in genetic markers was associated within some of the lake-stream divergence in morphological traits. Our results suggest that parallel evolution, and deviations from it, are primarily the result of natural selection, which corresponds in only some respects to the dichotomous habitat classifications frequently used in such studies.     

Genomic pattern of adaptive divergence in Arabidopsis halleri, a model species for tolerance to heavy metal

Pollution by heavy metals is one of the strongest environmental constraints in human-altered environments that only a handful of species can cope with. Identifying the genes conferring to those species the ability to grow in polluted areas is a first step towards a global understanding of the evolutionary processes involved and will eventually improve phytoremediation practices. We used a genome-scan approach to detect loci under divergent selection among four populations of Arabidopsis halleri growing on either polluted or nonpolluted habitats. Based on a high density of amplified fragment length polymorphism (AFLP) markers (820 AFLP markers, i.e. ~1 marker per 0.3 Mb), evidence for selection was found for some markers in every sampled population. Four loci departed from neutrality in both metallicolous populations and thus constitute high-quality candidates for general adaptation to pollution. Interestingly, some candidates differed between the two metallicolous populations, suggesting the possibility that different loci may be involved in adaptation in the different metallicolous populations.     

Scanning the genome for gene SNPs related to climate adaptation and estimating selection at the molecular level in boreal black spruce

Outlier detection methods were used to scan the genome of the boreal conifer black spruce (Picea mariana [Mill.] B.S.P.) for gene single-nucleotide polymorphisms (SNPs) potentially involved in adaptations to temperature and precipitation variations. The scan involved 583 SNPs from 313 genes potentially playing adaptive roles. Differentiation estimates among population groups defined following variation in temperature and precipitation were moderately high for adaptive quantitative characters such as the timing of budset or tree height (Q(ST) = 0.189-0.314). Average differentiation estimates for gene SNPs were null, with F(ST) values of 0.005 and 0.006, respectively, among temperature and precipitation population groups. Using two detection approaches, a total of 26 SNPs from 25 genes distributed among 11 of the 12 linkage groups of black spruce were detected as outliers with F(ST) as high as 0.078. Nearly half of the outlier SNPs were located in exons and half of those were nonsynonymous. The functional annotations of genes carrying outlier SNPs and regression analyses between the frequencies of these SNPs and climatic variables supported their implication in adaptive processes. Several genes carrying outlier SNPs belonged to gene families previously found to harbour outlier SNPs in a reproductively isolated but largely sympatric congeneric species, suggesting differential subfunctionalization of gene duplicates. Selection coefficient estimates (S) were moderate but well above the magnitude of drift (>1/N(e)), indicating that the signature of natural selection could be detected at the nucleotide level despite the recent establishment of these populations during the Holocene.     

An assessment of the spatial scale of local adaptation in brown trout (Salmo trutta L.): footprints of selection at microsatellite DNA loci

Local adaptation is considered a paradigm in studies of salmonid fish populations. Yet, little is known about the geographical scale of local adaptation. Is adaptive divergence primarily evident at the scale of regions or individual populations? Also, many salmonid populations are subject to spawning intrusion by farmed conspecifics that experience selection regimes fundamentally different from wild populations. This prompts the question if adaptive differences between wild populations and hatchery strains are more pronounced than between different wild populations? We addressed these issues by analyzing variation at 74 microsatellite loci (including anonymous and expressed sequence tag- and quantitative trait locus-linked markers) in 15 anadromous wild brown trout (Salmo trutta L.) populations, representing five geographical regions, along with two lake populations and two hatchery strains used for stocking some of the populations. F_ST-based outlier tests revealed more outlier loci between different geographical regions separated by 522±228 km (mean±s.d.) than between populations within regions separated by 117±79 km (mean±s.d.). A significant association between geographical distance and number of outliers between regions was evident. There was no evidence for more outliers in comparisons involving hatchery trout, but the loci under putative selection generally were not the same as those found to be outliers between wild populations. Our study supports the notion of local adaption being increasingly important at the scale of regions as compared with individual populations, and suggests that loci involved in adaptation to captive environments are not necessarily the same as those involved in adaptive divergence among wild populations.     

Admixture facilitates adaptation from standing variation in the European aspen (Populus tremula L.), a widespread forest tree

Adaptation to new environments can start from new mutations or from standing variation already present in natural populations. Whether admixture constrains or facilitates adaptation from standing variation is largely unknown, especially in ecological keystone or foundation species. We examined patterns of neutral and adaptive population divergence in Populus tremula L., a widespread forest tree, using mapped molecular genetic markers. We detected the genetic signature of postglacial admixture between a Western and an Eastern lineage of P. tremula in Scandinavia, an area suspected to represent a zone of postglacial contact for many species of animals and plants. Stringent divergence‐based neutrality tests provided clear indications for locally varying selection at the European scale. Six of 12 polymorphisms under selection were located less than 1 kb away from the nearest gene predicted by the Populus trichocarpa genome sequence. Few of these loci exhibited a signature of ‘selective sweeps’ in diversity‐based tests, which is to be expected if adaptation occurs primarily from standing variation. In Scandinavia, admixture explained genomic patterns of ancestry and the nature of clinal variation and strength of selection for bud set, a phenological trait of great adaptive significance in temperate trees, measured in a common garden trial. Our data provide a hitherto missing direct link between past range shifts because of climatic oscillations, and levels of standing variation currently available for selection and adaptation in a terrestrial foundation species.     

Population‐genomic approach reveals adaptive floral divergence in discrete populations of a hawk moth‐pollinated violet

Local adaptation to contrasting biotic or abiotic environments is an important evolutionary step that presumably precedes floral diversification at the species level, yet few studies have demonstrated the adaptive nature of intraspecific floral divergence in wild plant populations. We combine a population‐genomic approach with phenotypic information on floral traits to examine whether the differentiation in metric floral traits exhibited by 14 populations of the southern Spanish hawk moth‐pollinated violet Viola cazorlensis reflects adaptive divergence. Screening of many amplified fragment length polymorphism (AFLP) loci using a multiple‐marker‐based neutrality test identified nine outlier loci (2.6% of the total) that departed from neutral expectations and were potentially under selection. Generalized analysis of molecular variance revealed significant relationships between genetic distance and population divergence in three floral traits when genetic distance was based on outlier loci, but not when it was based on neutral ones. Population means of floral traits were closely correlated with population scores on the first principal coordinate axis of the genetic distance matrix using outlier loci, and with the allelic frequencies of four of the outlier loci. Results strongly support the adaptive nature of intraspecific floral divergence exhibited by V. cazorlensis and illustrate the potential of genome scans to identify instances of adaptive divergence when used in combination with phenotypic information.     

Reproductive isolation of sympatric morphs in a population of Darwin's finches

Recent research on speciation has identified a central role for ecological divergence, which can initiate speciation when (i) subsets of a species or population evolve to specialize on different ecological resources and (ii) the resulting phenotypic modes become reproductively isolated. Empirical evidence for these two processes working in conjunction, particularly during the early stages of divergence, has been limited. We recently described a population of the medium ground finch, Geospiza fortis, that features large and small beak morphs with relatively few intermediates. As in other Darwin's finches of the Galápagos Islands, these morphs presumably diverged in response to variation in local food availability and inter- or intraspecific competition. We here demonstrate that the two morphs show strong positive assortative pairing, a pattern that holds over three breeding seasons and during both dry and wet conditions. We also document restrictions on gene flow between the morphs, as revealed by genetic variation at 10 microsatellite loci. Our results provide strong support for the central role of ecology during the early stages of adaptive radiation.     

Evolutionary Conservation Genomics Reveals Recent Speciation and Local Adaptation in Threatened Takins

Incorrect species delimitation will lead to inappropriate conservation decisions, especially for threatened species. The takin (Budorcas taxicolor) is a large artiodactyl endemic to the Himalayan–Hengduan–Qinling Mountains and is well known for its threatened status and peculiar appearance. However, the speciation, intraspecies taxonomy, evolutionary history, and adaptive evolution of this species still remain unclear, which greatly hampers its scientific conservation. Here, we de novo assembled a high-quality chromosome-level genome of takin and resequenced the genomes of 75 wild takins. Phylogenomics revealed that takin was positioned at the root of Caprinae. Population genomics based on the autosome, X chromosome, and Y chromosome SNPs and mitochondrial genomes consistently revealed the existence of two phylogenetic species and recent speciation in takins: the Himalayan takin (B. taxicolor) and the Chinese takin (B. tibetana), with the support of morphological evidence. Two genetically divergent subspecies were identified in both takin species, rejecting three previously proposed taxonomical viewpoints. Furthermore, their distribution boundaries were determined, suggesting that large rivers play important roles in shaping the genetic partition. Compared with the other subspecies, the Qinling subspecies presented the lowest genomic diversity, higher linkage disequilibrium, inbreeding, and genetic load, thus is in urgent need of genetic management and protection. Moreover, coat color gene (PMEL) variation may be responsible for the adaptive coat color difference between the two species following Gloger’s rule. Our findings provide novel insights into the recent speciation, local adaptation, scientific conservation of takins, and biogeography of the Himalaya–Hengduan biodiversity hotspot.     

Molecular signatures of lineage‐specific adaptive evolution in a unique sea basin: the example of an anadromous goby Leucopsarion petersii

Climate changes on various time scales often shape genetic novelty and adaptive variation in many biotas. We explored molecular signatures of directional selection in populations of the ice goby Leucopsarion petersii inhabiting a unique sea basin, the Sea of Japan, where a wide variety of environments existed in the Pleistocene in relation to shifts in sea level by repeated glaciations. This species consisted of two historically allopatric lineages, the Japan Sea (JS) and Pacific Ocean (PO) lineages, and these have lived under contrasting marine environments that are expected to have imposed different selection regimes caused by past climatic and current oceanographic factors. We applied a limited genome‐scan approach using seven candidate genes for phenotypic differences between two lineages in combination with 100 anonymous microsatellite loci. Neuropeptide Y (NPY) gene, which is an important regulator of food intake and potent orexigenic agent, and three anonymous microsatellites were identified as robust outliers, that is, candidate loci potentially under directional selection, by multiple divergence‐ and diversity‐based outlier tests in comparisons focused on multiple populations of the JS vs. PO lineages. For these outlier loci, populations of the JS lineage had putative signals of selective sweeps. Additionally, real‐time quantitative PCR analysis using fish reared in a common environment showed a higher expression level for NPY gene in the JS lineage. Thus, this study succeeded in identifying candidate genomic regions under selection across populations of the JS lineage and provided evidence for lineage‐specific adaptive evolution in this unique sea basin.     

Comparative genomics and the study of evolution by natural selection

Genomics profoundly affects most areas of biology, including ecology and evolutionary biology. By examining genome sequences from multiple species, comparative genomics offers new insight into genome evolution and the way natural selection moulds DNA sequence evolution. Functional divergence, as manifested in the accumulation of nonsynonymous substitutions in protein-coding genes, differs among lineages in a manner seemingly related to population size. For example, the ratio of nonsynonymous to synonymous substitution (d_N/d_S) is higher in apes than in rodents, compatible with Ohta's nearly neutral theory of molecular evolution, which suggests that the fixation of slightly deleterious mutations contributes to protein evolution at an extent negatively correlated with effective population size. While this supports the idea that functional evolution is not necessarily adaptive, comparative genomics is uncovering a role for positive Darwinian selection in 10–40% of all genes in different lineages, estimates that are likely to increase when the addition of more genomes gives increased power. Again, population size seems to matter also in this context, with a higher proportion of fixed amino acid changes representing advantageous mutations in large populations. Genes that are particularly prone to be driven by positive selection include those involved with reproduction, immune response, sensory perception and apoptosis. Genetic innovations are also frequently obtained by the gain or loss of complete gene sequences. Moreover, it is increasingly realized, from comparative genomics, that purifying selection conserves much more than just the protein-coding part of the genome, and this points at an important role for regulatory elements in trait evolution. Finally, genome sequencing using outbred or multiple individuals has provided a wealth of polymorphism data that gives information on population history, demography and marker evolution.