Search

ABSTRACT Migration in animals and associated adaptations to contrasting environments are underpinned by complex genetic architecture. Here, we explore the genomic basis of facultative anadromy in brown trout (Salmo trutta), wherein some individuals migrate to sea while others remain resident in natal rivers, to better understand how alternative migratory tactics (AMTs) are maintained evolutionarily. To identify genomic variants associated with AMTs, we sequenced whole genomes for 194 individual trout from five anadromous–resident population pairs, situated above and below waterfalls, in five different Irish rivers. These waterfalls act as natural barriers to upstream migration and hence we predicted that loci underpinning AMTs should be under similar divergent selection across these replicate pairs. A sliding windows based analysis revealed a highly polygenic adaptive divergence between anadromous and resident populations, encompassing 329 differentiated genomic regions. These regions were associated with 292 genes involved in various processes crucial for AMTs, including energy homeostasis, reproduction, osmoregulation, immunity, circadian rhythm and neural function. Furthermore, examining patterns of diversity we were able to link specific genes and biological processes to putative AMT trait classes: migratory‐propensity, migratory‐lifestyle and residency. Importantly, AMT outlier regions possessed higher genetic diversity than the background genome, particularly in the anadromous group, suggesting balancing selection may play a role in maintaining genetic variation. Overall, the results from this study provide important insights into the genetic architecture of migration and the evolutionary mechanisms shaping genomic diversity within and across populations. [ABSTRACT FROM AUTHOR]

Abstract Local adaptations to the environment are an important aspect of the diversity of a species and their discovery, description and quantification has important implications for the fields of taxonomy, evolutionary and conservation biology. In this study, we scan genomes from several populations across the distributional range of the Eurasian lynx, with the objective of finding genomic windows under positive selection which may underlie local adaptations to different environments. A total of 394 genomic windows are found to be associated to local environmental conditions, and they are enriched for genes involved in metabolism, behaviour, synaptic organization and neural development. Adaptive genetic structure, reconstructed from SNPs in candidate windows, is considerably different than the neutral genetic structure of the species. A widespread adaptively homogeneous group is recovered occupying areas of harsher snow and temperature climatic conditions in the north‐western, central and eastern parts of the distribution. Adaptively divergent populations are recovered in the westernmost part of the range, especially within the Baltic population, but also predicted for different patches in the western and southern part of the range, associated with different snow and temperature regimes. Adaptive differentiation driven by climate does not correlate much with the subspecies taxonomic delimitations, suggesting that subspecific divergences are mostly driven by neutral processes of genetic drift and gene flow. Our results will aid the selection of source populations for assisted gene flow or genetic rescue programs by identifying what climatic patterns to look for as predictors of pre‐adaptation of individuals. Particularly, the Carpathian population is confirmed as the best source of individuals for the genetic rescue of the endangered, isolated and genetically eroded Balkan population. Additionally, reintroductions in central and western Europe, currently based mostly on Carpathian lynxes, could consider the Baltic population as an additional source to increase adaptive variation and likely improve adaptation to their milder climate.

Gene function conservation is crucial in molecular ecology, especially for key traits like growth and maturation in teleost fish. The vgll3 and six6 genes are known to influence age‐at‐maturity in Atlantic salmon, but their impact on other fish species is poorly understood. Here, we investigated the association of vgll3 and six6 in the domestication of gilthead seabream and European seabass, both undergoing selective breeding for growth‐related traits in the Mediterranean. We analysed two different sets of samples using two different genotyping approaches. The first dataset comprised farmed and wild populations from Greece, genotyped for SNPs within the two genes ('gene‐level genotyping'). The second dataset examined 300–600 k SNPs located in the chromosomes of the two genes, derived from a meta‐analysis of a Pool‐Seq experiment involving farmed and wild populations distributed widely across the Mediterranean ('chromosome‐level genotyping'). The gene‐level analysis revealed a statistically significant allele frequency differences between farmed and wild populations on both genes in each species. This finding was partially supported by the chromosome‐level analysis, identifying highly differentiated regions may be involved in the domestication process at varying distances from the candidate genes. Noteworthy genomic features were found, such as a CpG island in gilthead seabream and novel candidate genes in European seabass, warranting further investigation. These findings support a putative role of vgll3 and six6 in the maturation and growth of gilthead seabream and European seabass, emphasizing the need for further research on their conserved function. [ABSTRACT FROM AUTHOR]

Local adaptations to the environment are an important aspect of the diversity of a species and their discovery, description and quantification has important implications for the fields of taxonomy, evolutionary and conservation biology. In this study, we scan genomes from several populations across the distributional range of the Eurasian lynx, with the objective of finding genomic windows under positive selection which may underlie local adaptations to different environments. A total of 394 genomic windows are found to be associated to local environmental conditions, and they are enriched for genes involved in metabolism, behaviour, synaptic organization and neural development. Adaptive genetic structure, reconstructed from SNPs in candidate windows, is considerably different than the neutral genetic structure of the species. A widespread adaptively homogeneous group is recovered occupying areas of harsher snow and temperature climatic conditions in the north‐western, central and eastern parts of the distribution. Adaptively divergent populations are recovered in the westernmost part of the range, especially within the Baltic population, but also predicted for different patches in the western and southern part of the range, associated with different snow and temperature regimes. Adaptive differentiation driven by climate does not correlate much with the subspecies taxonomic delimitations, suggesting that subspecific divergences are mostly driven by neutral processes of genetic drift and gene flow. Our results will aid the selection of source populations for assisted gene flow or genetic rescue programs by identifying what climatic patterns to look for as predictors of pre‐adaptation of individuals. Particularly, the Carpathian population is confirmed as the best source of individuals for the genetic rescue of the endangered, isolated and genetically eroded Balkan population. Additionally, reintroductions in central and western Europe, currently based mostly on Carpathian lynxes, could consider the Baltic population as an additional source to increase adaptive variation and likely improve adaptation to their milder climate. [ABSTRACT FROM AUTHOR]

The detection of selective sweeps from population genomic data often relies on the premise that the beneficial mutations in question have fixed very near the sampling time. As it has been previously shown that the power to detect a selective sweep is strongly dependent on the time since fixation as well as the strength of selection, it is naturally the case that strong, recent sweeps leave the strongest signatures. However, the biological reality is that beneficial mutations enter populations at a rate, one that partially determines the mean wait time between sweep events and hence their age distribution. An important question thus remains about the power to detect recurrent selective sweeps when they are modeled by a realistic mutation rate and as part of a realistic distribution of fitness effects, as opposed to a single, recent, isolated event on a purely neutral background as is more commonly modeled. Here we use forward-in-time simulations to study the performance of commonly used sweep statistics, within the context of more realistic evolutionary baseline models incorporating purifying and background selection, population size change, and mutation and recombination rate heterogeneity. Results demonstrate the important interplay of these processes, necessitating caution when interpreting selection scans; specifically, false-positive rates are in excess of true-positive across much of the evaluated parameter space, and selective sweeps are often undetectable unless the strength of selection is exceptionally strong. [ABSTRACT FROM AUTHOR]

Evolutionary genetics has long struggled with understanding how functional genes under selection remain polymorphic in natural populations. Taking as a starting point that natural selection is ultimately a manifestation of ecological processes, we spotlight an underemphasized and potentially ubiquitous ecological effect that may have fundamental effects on the maintenance of genetic variation. Negative frequency dependency is a well‐established emergent property of density dependence in ecology, because the relative profitability of different modes of exploiting or utilizing limiting resources tends to be inversely proportional to their frequency in a population. We suggest that this may often generate negative frequency‐dependent selection (NFDS) on major effect loci that affect rate‐dependent physiological processes, such as metabolic rate, that are phenotypically manifested as polymorphism in pace‐of‐life syndromes. When such a locus under NFDS shows stable intermediate frequency polymorphism, this should generate epistatic selection potentially involving large numbers of loci with more minor effects on life‐history (LH) traits. When alternative alleles at such loci show sign epistasis with a major effect locus, this associative NFDS will promote the maintenance of polygenic variation in LH genes. We provide examples of the kind of major effect loci that could be involved and suggest empirical avenues that may better inform us on the importance and reach of this process. [ABSTRACT FROM AUTHOR]

Abstract Gene function conservation is crucial in molecular ecology, especially for key traits like growth and maturation in teleost fish. The vgll3 and six6 genes are known to influence age‐at‐maturity in Atlantic salmon, but their impact on other fish species is poorly understood. Here, we investigated the association of vgll3 and six6 in the domestication of gilthead seabream and European seabass, both undergoing selective breeding for growth‐related traits in the Mediterranean. We analysed two different sets of samples using two different genotyping approaches. The first dataset comprised farmed and wild populations from Greece, genotyped for SNPs within the two genes (‘gene‐level genotyping’). The second dataset examined 300–600 k SNPs located in the chromosomes of the two genes, derived from a meta‐analysis of a Pool‐Seq experiment involving farmed and wild populations distributed widely across the Mediterranean (‘chromosome‐level genotyping’). The gene‐level analysis revealed a statistically significant allele frequency differences between farmed and wild populations on both genes in each species. This finding was partially supported by the chromosome‐level analysis, identifying highly differentiated regions may be involved in the domestication process at varying distances from the candidate genes. Noteworthy genomic features were found, such as a CpG island in gilthead seabream and novel candidate genes in European seabass, warranting further investigation. These findings support a putative role of vgll3 and six6 in the maturation and growth of gilthead seabream and European seabass, emphasizing the need for further research on their conserved function.

Wheat was one of the crops domesticated in the Fertile Crescent region approximately 10,000 years ago. Despite undergoing recent polyploidization, hull-to-free-thresh transition events, and domestication bottlenecks, wheat is now grown in over 130 countries and accounts for a quarter of the world's cereal production. The main reason for its widespread success is its broad genetic diversity that allows it to thrive in different environments. To trace historical selection and hybridization signatures, genome scans were performed on two datasets: approximately 113K SNPs from 921 predominantly bread wheat accessions and approximately 110K SNPs from about 400 wheat accessions representing all ploidy levels. To identify environmental factors associated with the loci, a genome–environment association (GEA) was also performed. The genome scans on both datasets identified a highly differentiated region on chromosome 4A where accessions in the first dataset were dichotomized into a group (n = 691), comprising nearly all cultivars, wild emmer, and most landraces, and a second group (n = 230), dominated by landraces and spelt accessions. The grouping of cultivars is likely linked to their potential ancestor, bread wheat cv. Norin-10. The 4A region harbored important genes involved in adaptations to environmental conditions. The GEA detected loci associated with latitude and temperature. The genetic signatures detected in this study provide insight into the historical selection and hybridization events in the wheat genome that shaped its current genetic structure and facilitated its success in a wide spectrum of environmental conditions. The genome scans and GEA approaches applied in this study can help in screening the germplasm housed in gene banks for breeding, and for conservation purposes. [ABSTRACT FROM AUTHOR]

Adaptive genetic divergence occurs when selection imposed by the environment causes the genomic component of the phenotype to differentiate. However, genomic signatures of natural selection are usually identified without information on which trait is responding to selection by which selective agent(s). Here, we integrate whole‐genome sequencing with phenomics and measures of putative selective agents to assess the extent of adaptive divergence in threespine stickleback occupying the highly heterogeneous lake Mývatn, NE Iceland. We find negligible genome wide divergence, yet multiple traits (body size, gill raker structure and defence traits) were divergent along known ecological gradients (temperature, predatory bird densities and water depth). SNP based heritability of all measured traits was high (h2 = 0.42–0.65), indicating adaptive potential for all traits. Environment‐association analyses further identified thousands of loci putatively involved in selection, related to genes linked to, for instance, neuron development and protein phosphorylation. Finally, we found that loci linked to water depth were concurrently associated with pelvic spine length variation ‐ supporting the conclusion that divergence in pelvic spine length occurred in the face of gene flow. Our results suggest that whilst there is substantial genetic variation in the traits measured, phenotypic divergence of Mývatn stickleback is mostly weakly associated with environmental gradients, potentially as a result of substantial gene flow. Our study illustrates the value of integrative studies that combine genomic assays of multivariate trait variation with landscape genomics. [ABSTRACT FROM AUTHOR]

The blue shark Prionace glauca is a top predator with one of the widest geographical distributions of any shark species. It is classified as Critically Endangered in the Mediterranean Sea, and Near Threatened globally. Previous genetic studies did not reject the null hypothesis of a single global population. The blue shark was proposed as a possible archetype of the "grey zone of population differentiation," coined to designate cases where population structure may be too recent or too faint to be detected using a limited set of markers. Here, blue shark samples collected throughout its global range were sequenced using a specific RAD method (DArTseq), which recovered 37,655 genome‐wide single nucleotide polymorphisms (SNPs). Two main groups emerged, with Mediterranean Sea and northern Atlantic samples (Northern population) differentiated significantly from the Indo‐west Pacific samples (Southern population). Significant pairwise FST values indicated further genetic differentiation within the Atlantic Ocean, and between the Atlantic Ocean and the Mediterranean Sea. Reconstruction of recent demographic history suggested divergence between Northern and Southern populations occurred about 500 generations ago and revealed a drastic reduction in effective population size from a large ancestral population. Our results illustrate the power of genome scans to detect population structure and reconstruct demographic history in highly migratory marine species. Given that the management plans of the blue shark (targeted or bycatch) fisheries currently assume panmictic regional stocks, we strongly recommend that the results presented here be considered in future stock assessments and conservation strategies. [ABSTRACT FROM AUTHOR]

Understanding the contribution of neutral and adaptive evolutionary processes to population differentiation is often necessary for better informed management and conservation of rare species. In this study, we focused on Pinus torreyana Parry (Torrey pine), one of the world's rarest pines, endemic to one island and one mainland population in California. Small population size, low genetic diversity, and susceptibility to abiotic and biotic stresses suggest Torrey pine may benefit from interpopulation genetic rescue to preserve the species' evolutionary potential. We leveraged reduced representation sequencing to tease apart the respective contributions of stochastic and deterministic evolutionary processes to population differentiation. We applied these data to model spatial and temporal demographic changes in effective population sizes and genetic connectivity, to identify loci possibly under selection, and evaluate genetic rescue as a potential conservation strategy. Overall, we observed exceedingly low standing variation within both Torrey pine populations, reflecting consistently low effective population sizes across time, and limited genetic differentiation, suggesting maintenance of gene flow between populations following divergence. However, genome scans identified more than 2000 candidate SNPs potentially under divergent selection. Combined with previous observations indicating population phenotypic differentiation, this indicates natural selection has probably contributed to the evolution of population genetic differences. Thus, while reduced genetic diversity, small effective population size, and genetic connectivity between populations suggest genetic rescue could mitigate the adverse effects of rarity, evidence for adaptive differentiation suggests genetic mixing could disrupt adaptation. Further work evaluating the fitness consequences of inter‐population admixture is necessary to empirically evaluate the trade‐offs associated with genetic rescue in Torrey pine. [ABSTRACT FROM AUTHOR]

Abstract Parasitism is a pervasive phenomenon in nature with the relationship between species driving evolution in both parasite and host. Due to their host‐dependent lifestyle, parasites may adapt to the abiotic environment in ways that differ from their hosts or from free‐living relatives; yet rarely has this been assessed. Here, we test two competing hypotheses related to whether putatively adaptive genetic variation in a specialist mistletoe associates with the same, or different, climatic variables as its host species. We sampled 11 populations of the specialist mistletoe Amyema gibberula var. tatei (n = 154) and 10 populations of its associated host Hakea recurva subsp. recurva (n = 160). Reduced‐representation sequencing was used to obtain genome‐wide markers and putatively adaptive variation detected using genome scan methods. Climate associations were identified using generalized dissimilarity modelling, and these were mapped geographically to visualize the spatial patterns of genetic composition. Our results supported the hypothesis of parasites and host species responding differently to climatic variables. Temperature was relatively more important in predicting allelic turnover in the specialist mistletoe while precipitation was more important for the host. This suggests that parasitic plants and host species may respond differently to selective pressures, potentially as a result of differing nutrient acquisition strategies. Specifically, mistletoes acquire water from hosts (rather than the abiotic environment), which may provide a buffer to precipitation as a selective pressure. This work deepens and complements the physiological and other ecological studies of adaptation and provides a window into the evolutionary processes that underlie previously observed phenomena. Applying these methods to a comparative study in a host–parasite system has also highlighted factors that affect the study of selection pressure on nonmodel organisms, such as differing adaptation rates and lack of reference genomes.

Understanding the evolution of local adaptations is a central aim of evolutionary biology and key for the identification of unique populations and lineages of conservation relevance. By combining RAD sequencing and whole‐genome sequencing, we identify genetic signatures of local adaptation in mountain hares (Lepus timidus) from isolated and distinctive habitats of its wide distribution: Ireland, the Alps and Fennoscandia. Demographic modelling suggested that the split of these mountain hares occurred around 20 thousand years ago, providing the opportunity to study adaptive evolution over a short timescale. Using genome‐wide scans, we identified signatures of extreme differentiation among hares from distinct geographic areas that overlap with area‐specific selective sweeps, suggesting targets for local adaptation. Several identified candidate genes are associated with traits related to the uniqueness of the different environments inhabited by the three groups of mountain hares, including coat colour, ability to live at high altitudes and variation in body size. In Irish mountain hares, a variant of ASIP, a gene previously implicated in introgression‐driven winter coat colour variation in mountain and snowshoe hares (L. americanus), may underlie brown winter coats, reinforcing the repeated nature of evolution at ASIP moulding adaptive seasonal colouration. Comparative genomic analyses across several hare species suggested that mountain hares' adaptive variants appear predominantly species‐specific. However, using coalescent simulations, we also show instances where the candidate adaptive variants have been introduced via introgressive hybridization. Our study shows that standing adaptive variation, including that introgressed from other species, was a crucial component of the post‐glacial dynamics of species. [ABSTRACT FROM AUTHOR]

Abstract Spatially varying selection promotes variance in allele frequencies, increasing genetic differentiation between the demes of a metapopulation. For that reason, outliers in the genome‐wide distribution of summary statistics measuring genetic differentiation, such as FST, are often interpreted as evidence for alleles that contribute to local adaptation. However, theoretical studies have shown that in spatially structured populations the spread of beneficial mutations with spatially uniform fitness effects can also induce transient genetic differentiation. In recent years, numerous empirical studies have suggested that such species‐wide, or global, adaptation makes a substantial contribution to molecular evolution. In this perspective, we discuss how commonly such global adaptation may influence the genome‐wide distribution of FST and generate genetic differentiation patterns, which could be mistaken for local adaptation. To illustrate this, we use forward‐in‐time population genetic simulations assuming parameters for the rate and strength of beneficial mutations consistent with estimates from natural populations. We demonstrate that the spread of globally beneficial mutations in parapatric populations may frequently generate FST outliers, which could be misinterpreted as evidence for local adaptation. The spread of beneficial mutations causes selective sweeps at flanking sites, so in some cases, the effects of global versus local adaptation may be distinguished by examining patterns of nucleotide diversity within and between populations in addition to FST. However, when local adaptation has been only recently established, it may be much more difficult to distinguish from global adaptation, due to less accumulation of linkage disequilibrium at flanking sites. Through our discussion, we conclude that a large fraction of FST outliers that are presumed to arise from local adaptation may instead be due to global adaptation.

Uncovering the genetic and evolutionary basis of local adaptation is a major focus of evolutionary biology. The recent development of cost-effective methods for obtaining high-quality genome-scale data makes it possible to identify some of the loci responsible for adaptive differences among populations. Two basic approaches for identifying putatively locally adaptive loci have been developed and are broadly used: one that identifies loci with unusually high genetic differentiation among populations (differentiation outlier methods) and one that searches for correlations between local population allele frequencies and local environments (genetic-environment association methods). Here, we review the promises and challenges of these genome scan methods, including correcting for the confounding influence of a species' demographic history, biases caused by missing aspects of the genome, matching scales of environmental data with population structure, and other statistical considerations. In each case, we make suggestions for best practices for maximizing the accuracy and efficiency of genome scans to detect the underlying genetic basis of local adaptation. With attention to their current limitations, genome scan methods can be an important tool in finding the genetic basis of adaptive evolutionary change.

Parasitism is a pervasive phenomenon in nature with the relationship between species driving evolution in both parasite and host. Due to their host‐dependent lifestyle, parasites may adapt to the abiotic environment in ways that differ from their hosts or from free‐living relatives; yet rarely has this been assessed. Here, we test two competing hypotheses related to whether putatively adaptive genetic variation in a specialist mistletoe associates with the same, or different, climatic variables as its host species. We sampled 11 populations of the specialist mistletoe Amyema gibberula var. tatei (n = 154) and 10 populations of its associated host Hakea recurva subsp. recurva (n = 160). Reduced‐representation sequencing was used to obtain genome‐wide markers and putatively adaptive variation detected using genome scan methods. Climate associations were identified using generalized dissimilarity modelling, and these were mapped geographically to visualize the spatial patterns of genetic composition. Our results supported the hypothesis of parasites and host species responding differently to climatic variables. Temperature was relatively more important in predicting allelic turnover in the specialist mistletoe while precipitation was more important for the host. This suggests that parasitic plants and host species may respond differently to selective pressures, potentially as a result of differing nutrient acquisition strategies. Specifically, mistletoes acquire water from hosts (rather than the abiotic environment), which may provide a buffer to precipitation as a selective pressure. This work deepens and complements the physiological and other ecological studies of adaptation and provides a window into the evolutionary processes that underlie previously observed phenomena. Applying these methods to a comparative study in a host–parasite system has also highlighted factors that affect the study of selection pressure on nonmodel organisms, such as differing adaptation rates and lack of reference genomes. [ABSTRACT FROM AUTHOR]

Abstract Long‐term field studies coupled with quantitative genomics offer a powerful means to understand the genetic bases underlying quantitative traits and their evolutionary changes. However, analyzing and interpreting the time scales at which adaptive evolution occurs is challenging. First, while evolution is predictable in the short term, with strikingly rapid phenotypic changes in data series, it remains unpredictable in the long term. Second, while the temporal dynamics of some loci with large effects on phenotypic variation and fitness have been characterized, this task can be complicated in cases of highly polygenic trait architecture implicating numerous small effect size loci, or when statistical tests are sensitive to the heterogeneity of some key characteristics of the genome, like variation in recombination rate along the chromosomes. After introducing these aforementioned challenges, we discuss a recent investigation of the genomic architecture and spatio‐temporal variation in great tit bill length, which was related to the recent use of bird feeders. We discuss how this case study illustrates the importance of considering different temporal scales and evolutionary mechanisms both while analyzing trait temporal trends and when searching for and interpreting the signals of putative genomic footprints of selection. More generally this commentary discusses interesting challenges for unraveling the time scale at which adaptive traits evolve and their genomic bases.

Spatially varying selection promotes variance in allele frequencies, increasing genetic differentiation between the demes of a metapopulation. For that reason, outliers in the genome‐wide distribution of summary statistics measuring genetic differentiation, such as FST, are often interpreted as evidence for alleles that contribute to local adaptation. However, theoretical studies have shown that in spatially structured populations the spread of beneficial mutations with spatially uniform fitness effects can also induce transient genetic differentiation. In recent years, numerous empirical studies have suggested that such species‐wide, or global, adaptation makes a substantial contribution to molecular evolution. In this perspective, we discuss how commonly such global adaptation may influence the genome‐wide distribution of FST and generate genetic differentiation patterns, which could be mistaken for local adaptation. To illustrate this, we use forward‐in‐time population genetic simulations assuming parameters for the rate and strength of beneficial mutations consistent with estimates from natural populations. We demonstrate that the spread of globally beneficial mutations in parapatric populations may frequently generate FST outliers, which could be misinterpreted as evidence for local adaptation. The spread of beneficial mutations causes selective sweeps at flanking sites, so in some cases, the effects of global versus local adaptation may be distinguished by examining patterns of nucleotide diversity within and between populations in addition to FST. However, when local adaptation has been only recently established, it may be much more difficult to distinguish from global adaptation, due to less accumulation of linkage disequilibrium at flanking sites. Through our discussion, we conclude that a large fraction of FST outliers that are presumed to arise from local adaptation may instead be due to global adaptation. [ABSTRACT FROM AUTHOR]

Abstract Background The investigation of the genetic basis of local adaptation in non-model species is an interesting focus of evolutionary biologists and molecular ecologists. Identifying these adaptive genetic variabilities on the genome responsible can provide insight into the genetic mechanism of local adaptation. Results We investigated the spatial distribution of genetic variation in 22 natural populations of Pterocarya stenoptera across its distribution area in China to provide insights into the complex interplay between multiple environmental variables and adaptive genetic differentiation. The Bayesian analysis of population structure showed that the 22 populations of P. stenoptera were subdivided into two groups. Redundancy analysis demonstrated that this genetic differentiation was caused by the divergent selection of environmental difference. A total of 44 outlier loci were mutually identified by Arlequin and BayeScan, 43 of which were environment-associated loci (EAL). The results of latent factor mixed model analysis showed that solar radiation in June (Sr6), minimum temperature of the coldest month (Bio6), temperature seasonality (Bio4), and water vapor pressure in January (Wvp1) were associated with the highest numbers of EAL. Sr6 was associated with the ecological habitat of “prefered light”, and Bio6 and Wvp1 were associated with the ecological habitat of “warm and humid environment”. Conclusions Our results provided empirical evidence that environmental variables related to the ecological habitats of species play key roles in driving adaptive differentiation of species genome.

It has long been discussed to what extent related species develop similar genetic mechanisms to adapt to similar environments. Most studies documenting such convergence have either used different lineages within species or surveyed only a limited portion of the genome. Here, we investigated whether similar or different sets of orthologous genes were involved in genetic adaptation of natural populations of three related plant species to similar environmental gradients in the Alps. We used whole‐genome pooled population sequencing to study genome‐wide SNP variation in 18 natural populations of three Brassicaceae (Arabis alpina, Arabidopsis halleri, and Cardamine resedifolia) from the Swiss Alps. We first de novo assembled draft reference genomes for all three species. We then ran population and landscape genomic analyses with ~3 million SNPs per species to look for shared genomic signatures of selection and adaptation in response to similar environmental gradients acting on these species. Genes with a signature of convergent adaptation were found at significantly higher numbers than expected by chance. The most closely related species pair showed the highest relative over‐representation of shared adaptation signatures. Moreover, the identified genes of convergent adaptation were enriched for nonsynonymous mutations, suggesting functional relevance of these genes, even though many of the identified candidate genes have hitherto unknown or poorly described functions based on comparison with Arabidopsis thaliana. We conclude that adaptation to heterogeneous Alpine environments in related species is partly driven by convergent evolution, but that most of the genomic signatures of adaptation remain species‐specific. [ABSTRACT FROM AUTHOR]

In recent years, genome‐scan methods have been extensively used to detect local signatures of selection and introgression. Most of these methods are either designed for one or the other case, which may impair the study of combined cases. Here, we introduce a series of versatile genome‐scan methods applicable for both cases, the detection of selection and introgression. The proposed approaches are based on nonparametric k‐nearest neighbour (kNN) techniques, while incorporating pairwise Fixation Index (FST) and pairwise nucleotide differences (dxy) as features. We benchmark our methods using a wide range of simulation scenarios, with varying parameters, such as recombination rates, population background histories, selection strengths, the proportion of introgression and the time of gene flow. We find that kNN‐based methods perform remarkably well compared with the state‐of‐the‐art. Finally, we demonstrate how to perform kNN‐based genome scans on real‐world genomic data using the population genomics R‐package popgenome. [ABSTRACT FROM AUTHOR]

Convergent evolution, where independent lineages evolve similar phenotypes in response to similar challenges, can provide valuable insight into how selection operates and the limitations it encounters. However, it has only recently become possible to explore how convergent evolution is reflected at the genomic level. The overlapping outlier approach (OOA), where genome scans of multiple independent lineages are used to find outliers that overlap and therefore identify convergently evolving loci, is becoming popular. Here, we present a quantitative analysis of 34 studies that used this approach across many sampling designs, taxa, and sampling intensities. We found that OOA studies with increased biological sampling power within replicates have increased likelihood of finding overlapping, "convergent" signals of adaptation between them. When identifying convergent loci as overlapping outliers, it is tempting to assume that any false‐positive outliers derived from individual scans will fail to overlap across replicates, but this cannot be guaranteed. We highlight how population demographics and genomic context can contribute toward both true convergence and false positives in OOA studies. We finish with an exploration of emerging methods that couple genome scans with phenotype and environmental measures, leveraging added information from genome data to more directly test hypotheses of the likelihood of convergent evolution. [ABSTRACT FROM AUTHOR]

Sympatric tree species are subject to similar climatic drivers, posing a question as to whether they display comparable adaptive responses. However, no study has explicitly examined local adaptation of co‐occurring parasitic and autotrophic plant species to the abiotic environment. Here we test the hypotheses that a generalist parasitic tree would display a weaker signal of selection and that genomic variation would associate with fewer climatic variables (particularly precipitation) but have similar spatial patterns to a sympatric autotrophic tree species. To test these hypotheses, we collected samples from 17 sites across the range of two tree species, the hemiparasite Nuytsia floribunda (n = 264) and sympatric autotroph Melaleuca rhaphiophylla (n = 272). We obtained 5,531 high‐quality genome‐wide single nucleotide polymorphisms (SNPs) for M. rhaphiophylla and 6,727 SNPs for N. floribunda using DArTseq genome scan technology. Population differentiation and environmental association approaches were used to identify signals of selection. Generalized dissimilarly modelling was used to detect climatic and spatial patterns of local adaptation across climatic gradients. Overall, 322 SNPs were identified as putatively adaptive for the autotroph, while only 57 SNPs were identified for the parasitic species. We found genomic variation to associate with different sets of bioclimatic variables for each species, with precipitation relatively less important for the parasite. Spatial patterns of predicted adaptive variability were different and indicate that co‐occurring species with disparate life history traits may not respond equally to selective pressures (i.e., temperature and precipitation). Together, these findings provide insight into local adaptation of sympatric parasitic and autotrophic tree species to abiotic environments. [ABSTRACT FROM AUTHOR]

Summary: Latitudinal and elevational gradients provide valuable experimental settings for studies of the potential impact of global warming on forest tree species. The availability of long‐term phenological surveys in common garden experiments for traits associated with climate, such as bud flushing for sessile oaks (Quercus petraea), provide an ideal opportunity to investigate this impact.We sequenced 18 sessile oak populations and used available sequencing data for three other closely related European white oak species (Quercus pyrenaica, Quercus pubescens, and Quercus robur) to explore the evolutionary processes responsible for shaping the genetic variation across latitudinal and elevational gradients in extant sessile oaks. We used phenotypic surveys in common garden experiments and climatic data for the population of origin to perform genome‐wide scans for population differentiation and genotype–environment and genotype–phenotype associations.The inferred historical relationships between Q. petraea populations suggest that interspecific gene flow occurred between Q. robur and Q. petraea populations from cooler or wetter areas. A genome‐wide scan of differentiation between Q. petraea populations identified single nucleotide polymorphisms (SNPs) displaying strong interspecific relative divergence between these two species. These SNPs followed genetic clines along climatic or phenotypic gradients, providing further support for the likely contribution of introgression to the adaptive divergence of Q. petraea populations.Overall, the results indicate that outliers and associated SNPs are Q. robur ancestry‐informative. We discuss the results of this study in the framework of the postglacial colonization scenario, in which introgression and diversifying selection have been proposed as essential drivers of Q. petraea microevolution. See also the Editorial on this article by Plomion & Martin, 226: 943–946. [ABSTRACT FROM AUTHOR]

The identification of the genetic basis of domestication in fish species is of timely importance for the aquaculture industry in order to increase productivity, quality, and the welfare of farmed fish. The goal of this study is to investigate the largely unknown aquaculture-induced evolution in gilthead seabream, which is one of the most important farmed fish in the Mediterranean region. We used a panel of 1159 genome-wide SNPs, and genotyped 956 fish from 23 wild populations of Mediterranean-wide distribution and 362 farmed fish from five Greek hatcheries. We assessed the genetic diversity of the sampled populations and contrasted the results of four different approaches of outlier detection methods. We recognized one very strong candidate and two good candidate SNPs with evidence for aquaculture-induced evolution in gilthead seabream. The annotation of these SNPs revealed neighboring genes with biological roles from stress tolerance and disease resistance to sexual maturation that may explain our observations. In conclusion, we demonstrate that the genome of gilthead seabream, despite the fact that the species is often suggested to be in the early stages of the domestication process, shows evidence of aquaculture-induced evolution. We report on a list of genes that may explain our observations and that may be investigated further. We anticipate that our findings will stimulate additional research with the use of SNP panels of higher density that can elucidate the genomic architecture of domestication in this species of high aquacultural interest.

Long‐term field studies coupled with quantitative genomics offer a powerful means to understand the genetic bases underlying quantitative traits and their evolutionary changes. However, analyzing and interpreting the time scales at which adaptive evolution occurs is challenging. First, while evolution is predictable in the short term, with strikingly rapid phenotypic changes in data series, it remains unpredictable in the long term. Second, while the temporal dynamics of some loci with large effects on phenotypic variation and fitness have been characterized, this task can be complicated in cases of highly polygenic trait architecture implicating numerous small effect size loci, or when statistical tests are sensitive to the heterogeneity of some key characteristics of the genome, like variation in recombination rate along the chromosomes. After introducing these aforementioned challenges, we discuss a recent investigation of the genomic architecture and spatio‐temporal variation in great tit bill length, which was related to the recent use of bird feeders. We discuss how this case study illustrates the importance of considering different temporal scales and evolutionary mechanisms both while analyzing trait temporal trends and when searching for and interpreting the signals of putative genomic footprints of selection. More generally this commentary discusses interesting challenges for unraveling the time scale at which adaptive traits evolve and their genomic bases. [ABSTRACT FROM AUTHOR]

Identifying the evolutionary and developmental bases of adaptive phenotypes is of central interest in evolutionary biology. Cichlid fishes have been a useful research model due to their extraordinary phenotypic diversity reflecting adaptations to often very narrow niches. Among them, the scale-eating Perissodus microlepis is considered to be a textbook example for balanced polymorphism: its asymmetric head and handed behavior is thought to be maintained by negative frequency-dependent selection via prey-predator interactions. However, several contradictory findings and open questions have emerged in recent years, challenging our understanding of this model. Here, we review existing evidence for both genetic and non-genetic effects influencing head asymmetry, the association between morphological asymmetry and behavioral laterality, and the identification of signatures of balancing selection. Recent technological and theoretical developments have opened new exciting research avenues that can help identifying the drivers of adaptive traits in P. microlepis and other nonmodel organisms, and we discuss promising directions worth exploring. We highlight the importance of using integrative approaches that analyze genetic, environmental, and epigenetic variation in natural populations to aid a comprehensive understanding of why cichlids are so diverse and how evolution has produced and continues to generate such a vibrant and often complex phenotypic diversity. [ABSTRACT FROM AUTHOR]

High Andean wetlands are crucial for biodiversity conservation, climate change mitigation and adaptation, and the well-being of local communities in the Andes. Although community genetics provides an evolutionary viewpoint on ecological processes and allows for potential new understandings of the mechanisms that influence biodiversity and ecosystem function and their response to environmental changes, no research has yet explored the ecological effects of genetic variation in high Andean species. Here, we searched for single-nucleotide-polymorphism (SNP) loci of Carex gayana , an abundant plant species that fulfills crucial ecological functions in high Andean wetlands, that are associated with plant and aquatic macroinvertebrate diversity, and then investigated the spatial distribution of their alleles across a latitudinal gradient of ca. 600 km in Chile's Norte Chico. To identify SNP loci of C. gayana associated with community diversity, we analyzed a genomic dataset of SNP loci obtained through genotyping-by-sequencing of 158 individuals from 17 Andean wetlands using a genome scan approach. We found 30 structure and spatial outlier loci in C. gayana that were also linked with community parameters. We performed spatially-based causal modelling to infer the most probable origin of the observed relationships between SNP allele frequencies and community diversity parameters, resulting in 11 SNP candidates with possible extended effects on communities. The prevalence of alleles linked to greater taxonomic diversity varied greatly across wetlands. Some wetlands exhibited high frequencies of these alleles, while others had very low frequencies. These findings suggest that certain sites may be less resilient to global changes, while others could potentially serve as a source for restoration initiatives. Functional analysis of candidate genes using reference genomes of closely related Carex species allowed the identification of genes involved mainly in plant metabolism, responses to biotic and abiotic stress, defense mechanisms, growth and development, with functions commonly found in high-altitude plant species, and suggests that extended genetic effects of C. gayana genes to high Andean communities could be related to flower/reproductive development, defense mechanisms against herbivore insects and stress response. [Display omitted] • High Andean wetlands are important ecosystems impacted by global change. • Keystone species likely impose genetic-based effects on interacting communities. • Eleven SNP loci of Carex gayana with possible community genetic effects were found. • Candidate genes suggest biological processes vital for high altitude plants. [ABSTRACT FROM AUTHOR]

Ordination is a common tool in ecology that aims at representing complex biological information in a reduced space. In landscape genetics, ordination methods such as principal component analysis (PCA) have been used to detect adaptive variation based on genomic data. Taking advantage of environmental data in addition to genotype data, redundancy analysis (RDA) is another ordination approach that is useful to detect adaptive variation. This study aims at proposing a test statistic based on RDA to search for loci under selection. We compare redundancy analysis to pcadapt, which is a nonconstrained ordination method, and to a latent factor mixed model (LFMM), which is a univariate genotype–environment association method. Individual‐based simulations identify evolutionary scenarios where RDA genome scans have a greater statistical power than genome scans based on PCA. By constraining the analysis with environmental variables, RDA performs better than PCA in identifying adaptive variation when selection gradients are weakly correlated with population structure. In addition, we show that if RDA and LFMM have a similar power to identify genetic markers associated with environmental variables, the RDA‐based procedure has the advantage to identify the main selective gradients as a combination of environmental variables. To give a concrete illustration of RDA in population genomics, we apply this method to the detection of outliers and selective gradients on an SNP data set of Populus trichocarpa (Geraldes et al.,). The RDA‐based approach identifies the main selective gradient contrasting southern and coastal populations to northern and continental populations in the north‐western American coast. [ABSTRACT FROM AUTHOR]

As is true of virtually every realm of the biological sciences, our understanding of speciation is increasingly informed by the genomic revolution of the past decade. Investigators can ask detailed questions relating to both the extrinsic (e.g. inter- and intra-population and ecological interactions) and intrinsic (e.g. genome content and architecture) forces that drive speciation. Technologies ranging from restriction-site associated DNA sequencing (RADseq), to whole genome sequencing and assembly, to transcriptomics, to CRISPR are revolutionizing the means by which investigators can both frame and test hypotheses of lineage diversification. Our review aims to examine both extrinsic and intrinsic aspects of speciation. Genome-scale data have already served to fundamentally clarify the role of gene flow during (and after) speciation, although we predict that the differential propensity for speciation among phylogenetic lineages will be one of the most exciting frontiers for future genomic investigation. We propose that a unified theory of speciation will take into account the idiosyncratic features of genomic architecture examined in the light of each organism’s biology and ecology drawn from across the full breadth of the Tree of Life. [ABSTRACT FROM AUTHOR]

Genomic approaches are opening avenues for understanding all aspects of biological life, especially as they begin to be applied to multiple individuals and populations. However, these approaches typically depend on the availability of a sequenced genome for the species of interest. While the number of genomes being sequenced is exploding, one group that has lagged behind are weeds. Although the power of genomic approaches for weed science has been recognized, what is needed to implement these approaches is unfamiliar to many weed scientists. In this review we attempt to address this problem by providing a primer on genome sequencing and provide examples of how genomics can help answer key questions in weed science such as: (1) Where do agricultural weeds come from; (2) what genes underlie herbicide resistance; and, more speculatively, (3) can we alter weed populations to make them easier to control? This review is intended as an introduction to orient weed scientists who are thinking about initiating genome sequencing projects to better understand weed populations, to highlight recent publications that illustrate the potential for these methods, and to provide direction to key tools and literature that will facilitate the development and execution of weed genomic projects.

Understanding the genetic underpinnings of schistosome host preferences is critical. Luo et al. recently identified genes associated with intermediate and definitive host-switching based on a new chromosome-level genome for Schistosoma japonicum , population genetic comparisons, and follow-up experiments. This represents a guide to fully map-selected schistosome genes using population genetics. [ABSTRACT FROM AUTHOR]

Urbanization significantly alters natural ecosystems and has accelerated globally. Urban wildlife populations are often highly fragmented by human infrastructure, and isolated populations may adapt in response to local urban pressures. However, relatively few studies have identified genomic signatures of adaptation in urban animals. We used a landscape genomic approach to examine signatures of selection in urban populations of white-footed mice ( Peromyscus leucopus) in New York City. We analysed 154,770 SNPs identified from transcriptome data from 48 P. leucopus individuals from three urban and three rural populations and used outlier tests to identify evidence of urban adaptation. We accounted for demography by simulating a neutral SNP data set under an inferred demographic history as a null model for outlier analysis. We also tested whether candidate genes were associated with environmental variables related to urbanization. In total, we detected 381 outlier loci and after stringent filtering, identified and annotated 19 candidate loci. Many of the candidate genes were involved in metabolic processes and have well-established roles in metabolizing lipids and carbohydrates. Our results indicate that white-footed mice in New York City are adapting at the biomolecular level to local selective pressures in urban habitats. Annotation of outlier loci suggests selection is acting on metabolic pathways in urban populations, likely related to novel diets in cities that differ from diets in less disturbed areas. [ABSTRACT FROM AUTHOR]

With the highly dense genomic data available nowadays, ignoring linkage between genes would result in a huge loss of information. One way to prevent such a loss is to focus on the blocks of chromosomes shared identical by descent (IBD) in populations. The development of the theoretical framework modelling IBD processes is essential to support the advent of new tools such as haplotype phasing, imputation, inferring population structure and demographic history, mapping loci or detecting signatures of selection. This article aims to present the relevant models used in this context, and specify the underlying definitions of identity by descent that are yet to be gathered at one place. In light of this, we derived a general expression for the expected IBD block length, for any population model at any generation after founding. [ABSTRACT FROM AUTHOR]

Ecological processes, non-ecological processes or a combination of both may cause reproductive isolation and speciation, but their specific roles and potentially complex interactions in evolutionary radiations remain poorly understood, which defines a central knowledge gap at the interface of microevolution and macroevolution. Here I examine genome scans in combination with phenotypic and environmental data to disentangle how ecological and non-ecological processes contributed to population differentiation and speciation in an ongoing radiation of Lanistes gastropods from the Malawi Basin. I found a remarkable hierarchical structure of differentiation mechanisms in space and time: neutral and mutation-order processes are older and occur mainly between regions, whereas more recent adaptive processes are the main driver of genetic differentiation and reproductive isolation within regions. The strongest differentiation occurs between habitats and between regions, i.e. when ecological and non-ecological processes act synergistically. The structured occurrence of these processes based on the specific geographical setting and ecological opportunities strongly influenced the potential for evolutionary radiation. The results highlight the importance of interactions between various mechanisms of differentiation in evolutionary radiations, and suggest that non-ecological processes are important in adaptive radiations, including those of cichlids. Insight into such interactions is critical to understanding large-scale patterns of organismal diversity. [ABSTRACT FROM AUTHOR]

Evaluating the importance of coevolution for a wide range of evolutionary questions, such as the role parasites play in the evolution of sexual reproduction, requires that we understand the genetic basis of coevolutionary interactions. Despite its importance, little progress has been made identifying the genetic basis of coevolution, largely because we lack tools designed specifically for this purpose. Instead, coevolutionary studies are often forced to re-purpose single species techniques. Here, we propose a novel approach for identifying the genes mediating locally adapted coevolutionary interactions that relies on spatial correlations between genetic marker frequencies in the interacting species. Using individual-based multi-locus simulations, we quantify the performance of our approach across a range of coevolutionary genetic models. Our results show that when one species is strongly locally adapted to the other and a sufficient number of populations can be sampled, our approach accurately identifies functionally coupled host and parasite genes. Although not a panacea, the approach we outline here could help to focus the search for coevolving genes in a wide variety of well-studied systems for which substantial local adaptation has been demonstrated. [ABSTRACT FROM AUTHOR]

With the increasing availability of both molecular and topo-climatic data, the main challenges facing landscape genomics - that is the combination of landscape ecology with population genomics - include processing large numbers of models and distinguishing between selection and demographic processes (e.g. population structure). Several methods address the latter, either by estimating a null model of population history or by simultaneously inferring environmental and demographic effects. Here we present samβ ada, an approach designed to study signatures of local adaptation, with special emphasis on high performance computing of large-scale genetic and environmental data sets. samβ ada identifies candidate loci using genotype-environment associations while also incorporating multivariate analyses to assess the effect of many environmental predictor variables. This enables the inclusion of explanatory variables representing population structure into the models to lower the occurrences of spurious genotype-environment associations. In addition, samβ ada calculates local indicators of spatial association for candidate loci to provide information on whether similar genotypes tend to cluster in space, which constitutes a useful indication of the possible kinship between individuals. To test the usefulness of this approach, we carried out a simulation study and analysed a data set from Ugandan cattle to detect signatures of local adaptation with samβ ada, bayenv, lfmm and an FST outlier method ( FDIST approach in arlequin) and compare their results. samβ ada - an open source software for Windows, Linux and Mac OS X available at - outperforms other approaches and better suits whole-genome sequence data processing. [ABSTRACT FROM AUTHOR]

Analyses of genomewide polymorphism data have begun to shed light on speciation and adaptation. Genome scans to identify regions of the genome that are unusually different between populations or species, possibly due to divergent natural or sexual selection, are widespread in speciation genomics. Theoretical and empirical work suggests that such outlier regions may grow faster than linearly during speciation with gene flow due to a rapid transition between low and high reproductive isolation. We investigate whether this pattern could be attributed to neutral processes by simulating genomes under neutral evolution with varying amounts and timing of gene flow. Under both neutral evolution and divergent selection, simulations with little or no gene flow, or with a long allopatric period after its cessation, resulted in faster than linear growth of the proportion of the genome lying in outlier regions. Without selection, higher recent gene flow erased differentiation; with divergent selection, these same scenarios produced nonlinear growth to a plateau. Our results suggest that, given a history of gene flow, the growth of the divergent genome is informative about selection during divergence, but that in many scenarios, this pattern does not easily distinguish neutral and non-neutral processes during speciation with gene flow. [ABSTRACT FROM AUTHOR]

Selection has a deep impact on the distribution of genetic diversity and population differentiation along the genome (the genomic landscapes of diversity and differentiation), reducing diversity and elevating differentiation not only at the sites it targets, but also at linked neutral sites. Fuelled by the high-throughput sequencing revolution, these genomic footprints of selection have been extensively exploited over the past decade with the aim to identify genomic regions involved in adaptation and speciation. However, while this research has shown that the genomic landscapes of diversity and differentiation are usually highly heterogeneous, it has also led to the increasing realization that this heterogeneity may evolve under processes other than adaptation or speciation. In particular, instead of being an effect of selective sweeps or barriers to gene flow, accentuated differentiation can evolve by any process reducing genetic diversity locally within the genome (Charlesworth, ), including purifying selection at linked sites (background selection). In particular, in genomic regions where recombination is infrequent, accentuated differentiation can evolve as a by-product of diversity reductions unrelated to adaptation or speciation (Cruickshank & Hahn, ; Nachman & Payseur, ; Noor & Bennett, ). In such genomic regions, linkage extends over physically larger genome stretches, and selection affects a particularly high number of linked neutral sites. Even though the effects of selection on linked neutral diversity (linked selection) within populations are well documented (Cutter & Payseur, ), recent observations of diversity and differentiation landscapes that are highly correlated even among independent lineages suggest that the effects of long-term linked selection may have a deeper impact on the evolution of the genomic landscapes of diversity and differentiation than previously anticipated. The study on Saxicola stonechats by Van Doren et al. () reported in the current issue of Molecular Ecology lines in with a rapidly expanding body of evidence in this direction. Correlations of genomic landscapes extending from within stonechats to comparisons with Ficedula flycatchers add to recent insights into the timescales across which the effects of linked selection persist. Absent and inverted correlations of genomic landscapes in comparisons involving an island taxon, on the other hand, provide important empirical clues about the role of demographic constraints in the evolution of the genomic landscapes of diversity and differentiation. [ABSTRACT FROM AUTHOR]

Speciation, the evolution of reproductive isolation among populations, is continuous, complex, and involves multiple, interacting barriers. Until it is complete, the effects of this process vary along the genome and can lead to a heterogeneous genomic landscape with peaks and troughs of differentiation and divergence. When gene flow occurs during speciation, barriers restricting gene flow locally in the genome lead to patterns of heterogeneity. However, genomic heterogeneity can also be produced or modified by variation in factors such as background selection and selective sweeps, recombination and mutation rate variation, and heterogeneous gene density. Extracting the effects of gene flow, divergent selection and reproductive isolation from such modifying factors presents a major challenge to speciation genomics. We argue one of the principal aims of the field is to identify the barrier loci involved in limiting gene flow. We first summarize the expected signatures of selection at barrier loci, at the genomic regions linked to them and across the entire genome. We then discuss the modifying factors that complicate the interpretation of the observed genomic landscape. Finally, we end with a road map for future speciation research: a proposal for how to account for these modifying factors and to progress towards understanding the nature of barrier loci. Despite the difficulties of interpreting empirical data, we argue that the availability of promising technical and analytical methods will shed further light on the important roles that gene flow and divergent selection have in shaping the genomic landscape of speciation. [ABSTRACT FROM AUTHOR]

Identifying genomic regions underlying adaptation in extant lineages is key to understanding the trajectories along which biodiversity evolves. However, this task is complicated by evolutionary processes that obscure and mimic footprints of positive selection. Particularly, the long‐term effects of linked selection remain underappreciated and difficult to account for. Based on patterns emerging from recent research on the evolution of differentiation across the speciation continuum, I illustrate how long‐term linked selection affects the distribution of differentiation along genomes. I then argue that a comparative population genomics framework that exploits emergent features of long‐term linked selection can help overcome shortcomings of traditional genome scans for adaptive evolution, but needs to account for the temporal dynamics of differentiation landscapes. [ABSTRACT FROM AUTHOR]

Recent recognition that ecological and evolutionary processes can operate on similar timescales has led to a rapid increase in theoretical and empirical studies on eco-evolutionary dynamics. Progress in the fields of evolutionary biology, genomics and ecology is greatly enhancing our understanding of rapid adaptive processes, the predictability of adaptation and the genetics of ecologically important traits. However, progress in these fields has proceeded largely independently of one another. In an attempt to better integrate these fields, the centre for 'Adaptation to a Changing Environment' organized a conference entitled 'The genomic basis of eco-evolutionary change' and brought together experts in ecological genomics and eco-evolutionary dynamics. In this review, we use the work of the invited speakers to summarize eco-evolutionary dynamics and discuss how they are relevant for understanding and predicting responses to contemporary environmental change. Then, we show how recent advances in genomics are contributing to our understanding of eco-evolutionary dynamics. Finally, we highlight the gaps in our understanding of eco-evolutionary dynamics and recommend future avenues of research in eco-evolutionary dynamics. [ABSTRACT FROM AUTHOR]

Highly fecund marine species with dispersive life-history stages often display large population sizes and wide geographic distribution ranges. Consequently, they are expected to experience reduced genetic drift, efficient selection fueled by frequent adaptive mutations, and high migration loads. This has important consequences for understanding how local adaptation proceeds in the sea. A key issue in this regard, relates to the genetic architecture underlying fitness traits. Theory predicts that adaptation may involve many genes but with a high variance in effect size. Therefore, the effect of selection on allele frequencies may be substantial for the largest effect size loci, but insignificant for small effect genes. In such a context, the performance of population genomic methods to unravel the genetic basis of adaptation depends on the fraction of adaptive genetic variance explained by the cumulative effect of outlier loci. Here, we address some methodological challenges associated with the detection of local adaptation using molecular approaches. We provide an overview of genome scan methods to detect selection, including those assuming complex demographic models that better describe spatial population structure. We then focus on quantitative genetics approaches that search for genotype-phenotype associations at different genomic scales, including genome-wide methods evaluating the cumulative effect of variants. We argue that the limited power of single locus tests can be alleviated by the use of polygenic scores to estimate the joint contribution of candidate variants to phenotypic variation. [ABSTRACT FROM AUTHOR]

The spatial structure of the environment (e.g. the configuration of habitat patches) may play an important role in determining the strength of local adaptation. However, previous studies of habitat heterogeneity and local adaptation have largely been limited to simple landscapes, which poorly represent the multiscale habitat structure common in nature. Here, we use simulations to pursue two goals: (i) we explore how landscape heterogeneity, dispersal ability and selection affect the strength of local adaptation, and (ii) we evaluate the performance of several genotype-environment association (GEA) methods for detecting loci involved in local adaptation. We found that the strength of local adaptation increased in spatially aggregated selection regimes, but remained strong in patchy landscapes when selection was moderate to strong. Weak selection resulted in weak local adaptation that was relatively unaffected by landscape heterogeneity. In general, the power of detection methods closely reflected levels of local adaptation. False-positive rates (FPRs), however, showed distinct differences across GEA methods based on levels of population structure. The univariate GEA approach had high FPRs (up to 55%) under limited dispersal scenarios, due to strong isolation by distance. By contrast, multivariate, ordination-based methods had uniformly low FPRs (0-2%), suggesting these approaches can effectively control for population structure. Specifically, constrained ordinations had the best balance of high detection and low FPRs and will be a useful addition to the GEA toolkit. Our results provide both theoretical and practical insights into the conditions that shape local adaptation and how these conditions impact our ability to detect selection. [ABSTRACT FROM AUTHOR]