Genomic Analyses Capture the Human-Induced Demographic Collapse and Recovery in a Wide-Ranging Cervid.

The glacial cycles of the Quaternary heavily impacted species through successions of population contractions and expansions. Similarly, populations have been intensely shaped by human pressures such as unregulated hunting and land use changes. White-tailed and mule deer survived in different refugia through the Last Glacial Maximum, and their populations were severely reduced after the European colonization. Here, we analyzed 73 resequenced deer genomes from across their North American range to understand the consequences of climatic and anthropogenic pressures on deer demographic and adaptive history. We found strong signals of climate-induced vicariance and demographic decline; notably, multiple sequentially Markovian coalescent recovers a severe decline in mainland white-tailed deer effective population size (Ne) at the end of the Last Glacial Maximum. We found robust evidence for colonial overharvest in the form of a recent and dramatic drop in Ne in all analyzed populations. Historical census size and restocking data show a clear parallel to historical Ne estimates, and temporal Ne/Nc ratio shows patterns of conservation concern for mule deer. Signatures of selection highlight genes related to temperature, including a cold receptor previously highlighted in woolly mammoth. We also detected immune genes that we surmise reflect the changing land use patterns in North America. Our study provides a detailed picture of anthropogenic and climatic-induced decline in deer diversity and clues to understanding the conservation concerns of mule deer and the successful demographic recovery of white-tailed deer.

Island demographics and trait associations in white-tailed deer.

When a population is isolated and composed of few individuals, genetic drift is the paramount evolutionary force and results in the loss of genetic diversity. Inbreeding might also occur, resulting in genomic regions that are identical by descent, manifesting as runs of homozygosity (ROHs) and the expression of recessive traits. Likewise, the genes underlying traits of interest can be revealed by comparing fixed SNPs and divergent haplotypes between affected and unaffected individuals. Populations of white-tailed deer (Odocoileus virginianus) on islands of Saint Pierre and Miquelon (SPM, France) have high incidences of leucism and malocclusions, both considered genetic defects; on the Florida Keys islands (USA) deer exhibit smaller body sizes, a polygenic trait. Here we aimed to reconstruct island demography and identify the genes associated with these traits in a pseudo case-control design. The two island populations showed reduced levels of genomic diversity and a build-up of deleterious mutations compared to mainland deer; there was also significant genome-wide divergence in Key deer. Key deer showed higher inbreeding levels, but not longer ROHs, consistent with long-term isolation. We identified multiple trait-related genes in ROHs including LAMTOR2 which has links to pigmentation changes, and NPVF which is linked to craniofacial abnormalities. Our mixed approach of linking ROHs, fixed SNPs and haplotypes matched a high number (~50) of a-priori body size candidate genes in Key deer. This suite of biomarkers and candidate genes should prove useful for population monitoring, noting all three phenotypes show patterns consistent with a complex trait and non-Mendelian inheritance.

Speciation without gene-flow in hybridizing deer.

Under the ecological speciation model, divergent selection acts on ecological differences between populations, gradually creating barriers to gene flow and ultimately leading to reproductive isolation. Hybridisation is part of this continuum and can both promote and inhibit the speciation process. Here, we used white-tailed (Odocoileus virginianus) and mule deer (O. hemionus) to investigate patterns of speciation in hybridizing sister species. We quantified genome-wide historical introgression and performed genome scans to look for signatures of four different selection scenarios. Despite ample modern evidence of hybridisation, we found negligible patterns of ancestral introgression and no signatures of divergence with gene flow, rather localized patterns of allopatric and balancing selection were detected across the genome. Genes under balancing selection were related to immunity, MHC and sensory perception of smell, the latter of which is consistent with deer biology. The deficiency of historical gene-flow suggests that white-tailed and mule deer were spatially separated during the glaciation cycles of the Pleistocene and genome wide differentiation accrued via genetic drift. Dobzhansky-Muller incompatibilities and selection against hybrids are hypothesised to be acting, and diversity correlations to recombination rates suggests these sister species are far along the speciation continuum.

Characterizing nutritional phenotypes using experimental nutrigenomics: Is there nutrient-specificity to different types of dietary stress?

The ability to directly measure and monitor poor nutrition in individual animals and ecological communities is hampered by methodological limitations. In this study, we use nutrigenomics to identify nutritional biomarkers in a freshwater zooplankter, Daphnia pulex, a ubiquitous primary consumer in lakes and a sentinel of environmental change. We grew animals in six ecologically relevant nutritional treatments: nutrient replete, low carbon (food), low phosphorus, low nitrogen, low calcium and high Cyanobacteria. We extracted RNA for transcriptome sequencing to identify genes that were nutrient responsive and capable of predicting nutritional status with a high degree of accuracy. We selected a list of 125 candidate genes, which were subsequently pruned to 13 predictive potential biomarkers. Using a nearest-neighbour classification algorithm, we demonstrate that these potential biomarkers are capable of classifying our samples into the correct nutritional group with 100% accuracy. The functional annotation of the selected biomarkers revealed some specific nutritional pathways and supported our hypothesis that animal responses to poor nutrition are nutrient specific and not simply different presentations of slow growth or energy limitation. This is a key step in uncovering the causes and consequences of nutritional limitation in animal consumers and their responses to small- and large-scale changes in biogeochemical cycles.

Genomic health is dependent on long-term population demographic history.

Current genetic methods of population assessment in conservation biology have been challenged by genome-scale analyses due to their quantitatively novel insights. These analyses include assessments of runs-of-homozygosity (ROH), genomic evolutionary rate profiling (GERP), and mutational load. Here, we aim to elucidate the relationships between these measures using three divergent ungulates: white-tailed deer, caribou, and mountain goat. The white-tailed deer is currently expanding, while caribou are in the midst of a significant decline. Mountain goats remain stable, having suffered a large historical bottleneck. We assessed genome-wide signatures of inbreeding using the inbreeding coefficient F and %ROH (FROH ) and identified evolutionarily constrained regions with GERP. Mutational load was estimated by identifying mutations in highly constrained elements (CEs) and sorting intolerant from tolerant (SIFT) mutations. Our results showed that F and FROH are higher in mountain goats than in caribou and white-tailed deer. Given the extended bottleneck and low Ne of the mountain goat, this supports the idea that the genome-wide effects of demographic change take time to accrue. Similarly, we found that mountain goats possess more highly constrained CEs and the lowest dN/dS values, both of which are indicative of greater purifying selection; this is also reflected by fewer mutations in CEs and deleterious mutations identified by SIFT. In contrast, white-tailed deer presented the highest mutational load with both metrics, in addition to dN/dS, while caribou were intermediate. Our results demonstrate that extended bottlenecks may lead to reduced diversity and increased FROH in ungulates, but not necessarily an increase in mutational load, probably due to the purging of deleterious alleles in small populations.

Epigenetics and island-mainland divergence in an insectivorous small mammal.

Geographically isolated populations, specifically island-mainland counterparts, tend to exhibit phenotypic variation in many species. The so-called island syndrome occurs when different environmental pressures lead to insular divergence from mainland populations. This phenomenon can be seen in an island population of Nova Scotia masked shrews (Sorex cinereus), which have developed a specialized feeding habit and digestive enzyme compared to their mainland counterparts. Epigenetic modifications, such as DNA methylation (DNAm), can impact phenotypes by altering gene expression without changing the DNA sequence. Here, we used a de novo masked shrew genome assembly and a mammalian methylation array profiling 37 thousand conserved CpGs to investigate morphological and DNA methylation patterns between island and mainland populations. Island shrews were morphologically and epigenetically different than their mainland counterparts, exhibiting a smaller body size. A gene ontology enrichment analyses of differentially methylated CpGs implicated developmental and digestive system related pathways. Based on our shrew epigenetic clock, island shrews might also be aging faster than their mainland counterparts. This study provides novel insight on phenotypic and epigenetic divergence in island-mainland mammal populations and suggests an underlying role of methylation in island-mainland divergence.

Contrasting whole-genome and reduced representation sequencing for population demographic and adaptive inference: an alpine mammal case study.

Genomes capture the adaptive and demographic history of a species, but the choice of sequencing strategy and sample size can impact such inferences. We compared whole genome and reduced representation sequencing approaches to study the population demographic and adaptive signals of the North American mountain goat (Oreamnos americanus). We applied the restriction site-associated DNA sequencing (RADseq) approach to 254 individuals and whole genome resequencing (WGS) approach to 35 individuals across the species range at mid-level coverage (9X) and to 5 individuals at high coverage (30X). We used ANGSD to estimate the genotype likelihoods and estimated the effective population size (Ne), population structure, and explicitly modelled the demographic history with δaδi and MSMC2. The data sets were overall concordant in supporting a glacial induced vicariance and extremely low Ne in mountain goats. We evaluated a set of climatic variables and geographic location as predictors of genetic diversity using redundancy analysis. A moderate proportion of total variance (36% for WGS and 21% for RADseq data sets) was explained by geography and climate variables; both data sets support a large impact of drift and some degree of local adaptation. The empirical similarities of WGS and RADseq presented herein reassuringly suggest that both approaches will recover large demographic and adaptive signals in a population; however, WGS offers several advantages over RADseq, such as inferring adaptive processes and calculating runs-of-homozygosity estimates. Considering the predicted climate-induced changes in alpine environments and the genetically depauperate mountain goat, the long-term adaptive capabilities of this enigmatic species are questionable.

Genomic correlates for migratory direction in a free-ranging cervid.

Animal migrations are some of the most ubiquitous and one of the most threatened ecological processes globally. A wide range of migratory behaviours occur in nature, and this behaviour is not uniform among and within species, where even individuals in the same population can exhibit differences. While the environment largely drives migratory behaviour, it is necessary to understand the genetic mechanisms influencing migration to elucidate the potential of migratory species to cope with novel conditions and adapt to environmental change. In this study, we identified genes associated with a migratory trait by undertaking pooled genome-wide scans on a natural population of migrating mule deer. We identified genomic regions associated with variation in migratory direction, including FITM1, a gene linked to the formation of lipids, and DPPA3, a gene linked to epigenetic modifications of the maternal line. Such a genetic basis for a migratory trait contributes to the adaptive potential of the species and might affect the flexibility of individuals to change their behaviour in the face of changes in their environment.

Genomic diversity and demographic history of the Dromiciops genus (Marsupialia: Microbiotheriidae).

Three orders represent the South American fauna of marsupials. Of these, Microbiotheria was until recently known as a monotypic genus with the only surviving species Dromiciops gliroides (monito del monte). The recent proposal of a new Dromiciops species (Dromiciops bozinovici), together with new information on the origin and diversification of living microbioterians has changed the prevailing paradigm around the evolutionary history of these emblematic marsupials. Here, we used a RADseq approach to test for evidence of admixture and past or current gene flow among both species of Dromiciops and evaluate the genetic structure within D. gliroides. We analyzed 127 samples of Dromiciops distributed across the known distribution range of both species. We also inferred the joint demographic history of these lineages, thus corroborating the status of D. bozinovici as a distinct species. Demographic history reconstruction indicated that D. bozinovici diverged from D. gliroides around 4my ago and has remained isolated and demographically stable ever since. In contrast, D. gliroides is subdivided into three subclades that experienced recent expansions and moderate gene flow among them (mostly from north to south). Furthermore, genetic distances among populations within D. gliroides were significantly correlated with geographic distances. These results suggest that some of the D. gliroides populations would have survived in glacial refuges, with posterior expansions after ice retreat. Our results have important implications for the systematics of the genus and have profound conservation consequences for the new species, especially considering the fragmentation level of the temperate rainforest.

The biogeography of Dromiciops in southern South America: Middle Miocene transgressions, speciation and associations with Nothofagus.

The current distribution of the flora and fauna of southern South America is the result of drastic geological events that occurred during the last 20 million years, including marine transgressions, glaciations and active vulcanism. All these have been associated with fragmentation, isolation and subsequent expansion of the biota, south of 35°S, such as the temperate rainforest. This forest is mostly dominated by Nothofagus trees and is the habitat of the relict marsupial monito del monte, genus Dromiciops, sole survivor of the order Microbiotheria. Preliminary analyses using mtDNA proposed the existence of three main Dromiciops lineages, distributed latitudinally, whose divergence was initially attributed to recent Pleistocene glaciations. Using fossil-calibrated dating on nuclear and mitochondrial genes, here we reevaluate this hypothesis and report an older (Miocene) biogeographic history for the genus. We performed phylogenetic reconstructions using sequences from two mitochondrial DNA and four nuclear DNA genes in 159 specimens from 31 sites across Chile and Argentina. Our phylogenetic analysis resolved three main clades with discrete geographic distributions. The oldest and most differentiated clade corresponds to that of the northern distribution (35.2°S to 39.3°S), which should be considered a distinct species (D. bozinovici, sensu D'Elía et al. 2016). According to our estimations, this species shared a common ancestor with D. gliroides (southern clades) about ~13 million years ago. Divergence time estimates for the southern clades (39.6°S to 42.0°S) ranged from 9.57 to 6.5 Mya. A strong genetic structure was also detected within and between clades. Demographic analyses suggest population size stability for the northern clade (D. bozinovici), and recent demographic expansions for the central and southern clades. All together, our results suggest that the diversification of Dromiciops were initiated by the Middle Miocene transgression (MMT), the massive marine flooding that covered several lowlands of the western face of Los Andes between 37 and 48°S. The MMT resulted from an increase in global sea levels at the Miocene climatic optimum, which shaped the biogeographic origin of several species, including Nothofagus forests, the habitat of Dromiciops.

The complexity of co-limitation: nutrigenomics reveal non-additive interactions of calcium and phosphorus on gene expression in Daphnia pulex.

Many lakes across Canada and northern Europe have experienced declines in ambient phosphorus (P) and calcium (Ca) supply for over 20 years. While these declines might create or exacerbate nutrient limitation in aquatic food webs, our ability to detect and quantify different types of nutrient stress on zooplankton remains rudimentary. Here, we used growth bioassay experiments and whole transcriptome RNAseq, collectively nutrigenomics, to examine the nutritional phenotypes produced by low supplies of P and Ca separately and together in the freshwater zooplankter Daphnia pulex. We found that daphniids in all three nutrient-deficient categories grew slower and differed in their elemental composition. Our RNAseq results show distinct responses in singly limited treatments (Ca or P) and largely a mix of these responses in animals under low Ca and P conditions. Deeper investigation of effect magnitude and gene functional annotations reveals this patchwork of responses to cumulatively represent a co-limited nutritional phenotype. Linear discriminant analysis identified a significant separation between nutritional treatments based upon gene expression patterns with the expression patterns of just five genes needed to predict animal nutritional status with 99% accuracy. These data reveal how nutritional phenotypes are altered by individual and co-limitation of two highly important nutritional elements (Ca and P) and provide evidence that aquatic consumers can respond to limitation by more than one nutrient at a time by differentially altering their metabolism. This use of nutrigenomics demonstrates its potential to address many of the inherent complexities in studying interactions between multiple nutritional stressors in ecology and beyond.

Population demographic history and adaptability of the vulnerable Lolokou Sucker Frog.

Amphibians are experiencing worldwide declines due to increasing anthropogenetic disturbances. However, the genetic variability and hence adaptability are still unknown for most frogs. We integrated the mitochondrial (ND2 gene), nuclear (TYR gene) and major histocompatibility complex (MHC) loci, to clarify the demographic patterns and immune-gene diversity of the Lolokou Sucker Frog (Amolops loloensis). Demographic analysis of the ND2 and TYR genes suggested that the Lolokou Sucker Frog experienced a population expansion within the last 10,000 years. High MHC diversity was detected, which has likely resulted from positive selection, indicating the current diversity bodes well for the species' adaptive potential to pathogenic challenges. These findings broaden our knowledge on the population history and evolution adaptation of the reclusive torrent frog, and conservation implications are provided.

Heritability Estimates of Antler and Body Traits in White-Tailed Deer (Odocoileus virginianus) From Genomic-Relatedness Matrices.

Estimating heritability (h2) is required to predict the response to selection and is useful in species that are managed or farmed using trait information. Estimating h2 in free-ranging populations is challenging due to the need for pedigrees; genomic-relatedness matrices (GRMs) circumvent this need and can be implemented in nearly any system where phenotypic and genome-wide single-nucleotide polymorphism (SNP) data are available. We estimated the heritability of 5 body and 3 antler traits in a free-ranging population of white-tailed deer (Odocoileus virginianus) on Anticosti Island, Quebec, Canada. We generated classic and robust GRMs from >10,000 SNPs: hind foot length, dressed body mass, and peroneus muscle mass had high h2 values of 0.62, 0.44, and 0.55, respectively. Heritability in male-only antler features ranged from 0.07 to 0.33. We explored the influence of filtering by minor allele frequency and data completion on h2: GRMs derived from fewer SNPs had reduced h2 estimates and the relatedness coefficients significantly deviated from those generated with more SNPs. As a corollary, we discussed limitations to the application of GRMs in the wild, notably how skewed GRMs, specifically many unrelated individuals, can increase variance around h2 estimates. This is the first study to estimate h2 on a free-ranging population of white-tailed deer and should be informative for breeding designs and management as these traits could respond to selection.

Genome-scale sampling suggests cryptic epigenetic structuring and insular divergence in Canada lynx.

Determining the molecular signatures of adaptive differentiation is a fundamental component of evolutionary biology. A key challenge is to identify such signatures in wild organisms, particularly between populations of highly mobile species that undergo substantial gene flow. The Canada lynx (Lynx canadensis) is one species where mainland populations appear largely undifferentiated at traditional genetic markers, despite inhabiting diverse environments and displaying phenotypic variation. Here, we used high-throughput sequencing to investigate both neutral genetic structure and epigenetic differentiation across the distributional range of Canada lynx. Newfoundland lynx were identified as the most differentiated population at neutral genetic markers, with demographic modelling suggesting that divergence from the mainland occurred at the end of the last glaciation (20-33 KYA). In contrast, epigenetic structure revealed hidden levels of differentiation across the range coincident with environmental determinants including winter conditions, particularly in the peripheral Newfoundland and Alaskan populations. Several biological pathways related to morphology were differentially methylated between populations, suggesting that epigenetic modifications might explain morphological differences seen between geographically peripheral populations. Our results indicate that epigenetic modifications, specifically DNA methylation, are powerful markers to investigate population differentiation in wild and non-model systems.

Forecasting Ecological Genomics: High-Tech Animal Instrumentation Meets High-Throughput Sequencing.

Recent advancements in animal tracking technology and high-throughput sequencing are rapidly changing the questions and scope of research in the biological sciences. The integration of genomic data with high-tech animal instrumentation comes as a natural progression of traditional work in ecological genetics, and we provide a framework for linking the separate data streams from these technologies. Such a merger will elucidate the genetic basis of adaptive behaviors like migration and hibernation and advance our understanding of fundamental ecological and evolutionary processes such as pathogen transmission, population responses to environmental change, and communication in natural populations.

Disruptive selection without genome-wide evolution across a migratory divide.

Transcontinental migration is a fascinating example of how animals can respond to climatic oscillation. Yet, quantitative data on fitness components are scarce, and the resulting population genetic consequences are poorly understood. Migratory divides, hybrid zones with a transition in migratory behaviour, provide a natural setting to investigate the micro-evolutionary dynamics induced by migration under sympatric conditions. Here, we studied the effects of migratory programme on survival, trait evolution and genome-wide patterns of population differentiation in a migratory divide of European barn swallows. We sampled a total of 824 individuals from both allopatric European populations wintering in central and southern Africa, respectively, along with two mixed populations from within the migratory divide. While most morphological characters varied by latitude consistent with Bergmann's rule, wing length co-varied with distance to wintering grounds. Survival data collected during a 5-year period provided strong evidence that this covariance is repeatedly generated by disruptive selection against intermediate phenotypes. Yet, selection-induced divergence did not translate into genome-wide genetic differentiation as assessed by microsatellites, mtDNA and >20 000 genome-wide SNP markers; nor did we find evidence of local genomic selection between migratory types. Among breeding populations, a single outlier locus mapped to the BUB1 gene with a role in mitotic and meiotic organization. Overall, this study provides evidence for an adaptive response to variation in migration behaviour continuously eroded by gene flow under current conditions of nonassortative mating. It supports the theoretical prediction that population differentiation is difficult to achieve under conditions of gene flow despite measurable disruptive selection.

Cross-Species Application of SNP Chips is Not Suitable for Identifying Runs of Homozygosity.

Cross-species application of single-nucleotide polymorphism (SNP) chips is a valid, relatively cost-effective alternative to the high-throughput sequencing methods generally required to obtain a genome-wide sampling of polymorphisms. Kharzinova et al. (2015) examined the applicability of SNP chips developed in domestic bovids (cattle and sheep) to a semi-wild cervid (reindeer). The ancestors of bovids and cervids diverged between 20 and 30 million years ago (Hassanin and Douzery 2003; Bibi et al. 2013). Empirical work has shown that for a SNP chip developed in a bovid and applied to a cervid species, approximately 50% genotype success with 1% of the loci being polymorphic is expected (Miller et al. 2012). The genotyping of Kharzinova et al. (2015) follows this pattern; however, these data are not appropriate for identifying runs of homozygosity (ROH) and can be problematic for estimating linkage disequilibrium (LD) and we caution readers in this regard.

Demographic inferences using short-read genomic data in an approximate Bayesian computation framework: in silico evaluation of power, biases and proof of concept in Atlantic walrus.

Approximate Bayesian computation (ABC) is a powerful tool for model-based inference of demographic histories from large genetic data sets. For most organisms, its implementation has been hampered by the lack of sufficient genetic data. Genotyping-by-sequencing (GBS) provides cheap genome-scale data to fill this gap, but its potential has not fully been exploited. Here, we explored power, precision and biases of a coalescent-based ABC approach where GBS data were modelled with either a population mutation parameter (θ) or a fixed site (FS) approach, allowing single or several segregating sites per locus. With simulated data ranging from 500 to 50 000 loci, a variety of demographic models could be reliably inferred across a range of timescales and migration scenarios. Posterior estimates were informative with 1000 loci for migration and split time in simple population divergence models. In more complex models, posterior distributions were wide and almost reverted to the uninformative prior even with 50 000 loci. ABC parameter estimates, however, were generally more accurate than an alternative composite-likelihood method. Bottleneck scenarios proved particularly difficult, and only recent bottlenecks without recovery could be reliably detected and dated. Notably, minor-allele-frequency filters - usual practice for GBS data - negatively affected nearly all estimates. With this in mind, we used a combination of FS and θ approaches on empirical GBS data generated from the Atlantic walrus (Odobenus rosmarus rosmarus), collectively providing support for a population split before the last glacial maximum followed by asymmetrical migration and a high Arctic bottleneck. Overall, this study evaluates the potential and limitations of GBS data in an ABC-coalescence framework and proposes a best-practice approach.

Epigenetic clocks, sex markers and age-class diagnostics in three harvested large mammals.

The development of epigenetic clocks, or the DNA methylation-based inference of age, is an emerging tool for ageing in free ranging populations. In this study, we developed epigenetic clocks for three species of large mammals that are the focus of extensive management throughout their range in North America: white-tailed deer, black bear and mountain goat. We quantified differential DNA methylation patterns at over 30,000 cytosine-guanine sites (CpGs) from tissue samples of all three species (black bear n = 49; white-tailed deer n = 47; mountain goat n = 45). We used a penalized regression model (elastic net) to build explanatory (black bear r = .95; white-tailed deer r = .99; mountain goat r = .97) and robust (black bear Median Absolute Error or MAE = 1.33; white-tailed deer MAE = 0.29; mountain goat MAE = 0.61) models of age or clocks. We also characterized individual CpG sites within each species that demonstrated clear differences in methylation levels between age classes and sex, which can be used to develop a suite of accessible diagnostic markers. This tool has the potential to contribute to wildlife monitoring by providing easily obtainable representations of age structure in managed populations.

Widespread evidence for incipient ecological speciation: a meta-analysis of isolation-by-ecology.

Ecologically mediated selection has increasingly become recognised as an important driver of speciation. The correlation between neutral genetic differentiation and environmental or phenotypic divergence among populations, to which we collectively refer to as isolation-by-ecology (IBE), is an indicator of ecological speciation. In a meta-analysis framework, we determined the strength and commonality of IBE in nature. On the basis of 106 studies, we calculated a mean effect size of IBE with and without controlling for spatial autocorrelation among populations. Effect sizes were 0.34 (95% CI 0.24-0.42) and 0.26 (95% CI 0.13-0.37), respectively, indicating that an average of 5% of the neutral genetic differentiation among populations was explained purely by ecological contrast. Importantly, spatial autocorrelation reduced IBE correlations for environmental variables, but not for phenotypes. Through simulation, we showed how the influence of isolation-by-distance and spatial autocorrelation of ecological variables can result in false positives or underestimated correlations if not accounted for in the IBE model. Collectively, this meta-analysis showed that ecologically induced genetic divergence is pervasive across time-scales and taxa, and largely independent of the choice of molecular marker. We discuss the importance of these results in the context of adaptation and ecological speciation and suggest future research avenues.