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.     

FstSNP-HapMap3: a database of SNPs with high population differentiation for HapMap3

The International HapMap Project has recently made available genotypes and frequency data for phase 3 (NCBI build 36, dbSNPb129) of the HapMap providing an enriched genotype dataset for approximately 1.6 million single nucleotide polymorphisms (SNPs) from 1,115 individuals with ancestry from parts of Africa, Asia, Europe, North America and Mexico. In the present study, we aim to facilitate pharmacogenetics studies by providing a database of SNPs with high population differentiation through a genomewide test on allele frequency variation among 11 HapMap3 samples. Common SNPs with minor allele frequency greater than 5¢ from each of 11 HapMap3 samples were included in the present analysis. The population differentiation is measured in terms of fixation index (Fst), and the SNPs with Fst values over 0.5 were defined as highly differentiated SNPs. Our tests were carried out between all pairs of the 11 HapMap3 samples or among subgroups with the same continental ancestries. Altogether we carried out 64 genomewide Fst tests and identified 28,215 highly differentiated SNPs for 49 different combinations of HapMap3 samples in the current database.     

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.     

Natural selection in the water: freshwater invasion and adaptation by water colour in the Amazonian pufferfish

Natural selection and ecological adaptation are ultimately responsible for much of the origin of biodiversity. Yet, the identification of divergent natural selection has been hindered by the spatial complexity of natural systems, the difficulty in identifying genes under selection and their relationship to environment, and the confounding genomic effects of time. Here, we employed genome scans, population genetics and sequence-based phylogeographic methods to identify divergent natural selection on population boundaries in a freshwater invader, the Amazonian pufferfish, Colomesus asellus. We sampled extensively across markedly different hydrochemical settings in the Amazon Basin and use 'water colour' to test for ecological isolation. We distinguish the relative contribution of natural selection across hydrochemical gradients from biogeographic history in the origin and maintenance of population boundaries within a single species and across a complex ecosystem. We show that spatially distinct population structure generated by multiple forces (i.e. water colour and vicariant biogeographic history) can be identified if the confounding effects of genetic drift have not accumulated between selective populations. Our findings have repercussions for studies aimed at identifying engines of biodiversity and assessing their temporal progression in understudied and ecologically complex tropical ecosystems.     

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.     

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.     

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.     

Genetic differentiation,hybridization and adaptive divergence in two subspecies of the acorn barnacle Tetraclita japonica in the northwestern Pacific

Two acorn barnacles, Tetraclita japonica japonica and Tetraclita japonica formosana, have been recently reclassified as two subspecies, because they are morphologically similar and genetically indistinguishable in mitochondrial DNA sequences. The two barnacles are distinguishable by parietes colour and exhibit parapatric distributions, coexisting in Japan, where T. j. formosana is very low in abundance. Here we investigated the genetic differentiation between the subspecies using 209 polymorphic amplified fragment length polymorphism markers and 341 individuals from 12 locations. The subspecies are genetically highly differentiated (Φ_CT = 0.267). Bayesian analysis and principal component analysis indicate the presence of hybrids in T. j. formosana samples from Japan. Strong differentiation between the northern and southern populations of T. j. japonica was revealed, and a break between Taiwan and Okinawa was also found in T. j. formosana. The differentiation between the two taxa at individual loci does not deviate from neutral expectation, suggesting that the oceanographic pattern which restricts larval dispersal is a more important factor than divergent selection in maintaining genetic and phenotypic differentiation. The T. j. formosana in Japan are probably recent migrants from Okinawa, and their presence in Japan may represent a poleward range shift driven by global warming. This promotes hybridization and might lead to a breakdown of the boundary between the subspecies. However, both local adaptation and larval dispersal are crucial in determining the population structure within each subspecies. Our study provides new insights into the interplay of local adaptation and dispersal in determining the distribution and genetic structure of intertidal biota and the biogeography of the northwestern Pacific.     

The n = 1 constraint in population genomics

A key objective of population genomics is to identify portions of the genome that have been shaped by natural selection rather than by neutral divergence. A previously recognized but underappreciated challenge to this objective is that observations of allele frequencies across genomes in natural populations often correspond to a single, unreplicated instance of the outcome of evolution. This is because the composition of each individual genomic region and population is expected to be the outcome of a unique array of evolutionary processes. Given a single observation, inference of the evolutionary processes that led to the observed state of a locus is associated with considerable uncertainty. This constraint on inference can be ameliorated by utilizing multi-allelic (e.g. DNA haplotypes) rather than bi-allelic markers, by analysing two or more populations with certain models and by utilizing studies of replicated experimental evolution. Future progress in population genomics will follow from research that recognizes the 'n = 1 constraint' and that utilizes appropriate and explicit evolutionary models for analysis.