The search for loci under selection: trends, biases and progress

Adaptation to the abiotic environment in insects: the influence of variability on ecophysiology and evolutionary genomics

Advances in tools to gather environmental, phenotypic, and molecular data have accelerated our ability to detect abiotic drivers of variation across the genome-to-phenome spectrum in model and non-model insects. However, differences in the spatial and temporal resolution of these data sets may create gaps in our understanding of linkages between environment, genotype, and phenotype that yield missed or misleading results about adaptive variation. In this review we highlight sources of variability that might impact studies of phenotypic and 'omic environmental adaptation, challenges to collecting data at relevant scales, and possible solutions that link intensive fine-scale reductionist studies of mechanisms to large-scale biogeographic patterns.

Sex matters: Otolith shape and genomic variation in deacon rockfish (Sebastes diaconus)

Little is known about intraspecific variation within the deacon rockfish (Sebastes diaconus), a recently described species found in the northeast Pacific Ocean. We investigated population structure among fish sampled from two nearshore reefs (Siletz Reef and Seal Rock) and one offshore site (Stonewall Bank) within a <50-km2 area off the Oregon coast. Fish from the three sample sites exhibited small but statistically significant differences based on genetic variation at >15,000 neutral loci, whether analyzed independently or classified into nearshore and offshore groups. Male and females were readily distinguished using genetic data and 92 outlier loci were associated with sex, potentially indicating differential selection between males and females. Morphometric results indicated that there was significant secondary sexual dimorphism in otolith shape, but further sampling is required to disentangle potential confounding influence of age. This study is the first step toward understanding intraspecific variation within the deacon rockfish and the potential management implications. Since differentiation among the three sample sites was small, we consider the results to be suggestive of a single stock. However, future studies should evaluate how the stock is affected by differences in sex, age, and gene flow between the nearshore and offshore environments.

Landscape drivers of genomic diversity and divergence in woodland Eucalyptus

Spatial genetic patterns are influenced by numerous factors, and they can vary even among coexisting, closely related species due to differences in dispersal and selection. Eucalyptus (L'Héritier 1789; the "eucalypts") are foundation tree species that provide essential habitat and modulate ecosystem services throughout Australia. Here we present a study of landscape genomic variation in two woodland eucalypt species, using whole-genome sequencing of 388 individuals of Eucalyptus albens and Eucalyptus sideroxylon. We found exceptionally high genetic diversity (π ≈ 0.05) and low genome-wide, interspecific differentiation (FST  = 0.15) and intraspecific differentiation between localities (FST  ≈ 0.01-0.02). We found no support for strong, discrete population structure, but found substantial support for isolation by geographic distance (IBD) in both species. Using generalized dissimilarity modelling, we identified additional isolation by environment (IBE). Eucalyptus albens showed moderate IBD, and environmental variables have a small but significant amount of additional predictive power (i.e. IBE). Eucalyptus sideroxylon showed much stronger IBD and moderate IBE. These results highlight the vast adaptive potential of these species and set the stage for testing evolutionary hypotheses of interspecific adaptive differentiation across environments.

Predicting the success of an invader: Niche shift versus niche conservatism

Invasive species can encounter environments different from their source populations, which may trigger rapid adaptive changes after introduction (niche shift hypothesis). To test this hypothesis, we investigated whether postintroduction evolution is correlated with contrasting environmental conditions between the European invasive and source ranges in the Asian tiger mosquito Aedes albopictus. The comparison of environmental niches occupied in European and source population ranges revealed more than 96% overlap between invasive and source niches, supporting niche conservatism. However, we found evidence for postintroduction genetic evolution by reanalyzing a published ddRADseq genomic dataset from 90 European invasive populations using genotype-environment association (GEA) methods and generalized dissimilarity modeling (GDM). Three loci, among which a putative heat-shock protein, exhibited significant allelic turnover along the gradient of winter precipitation that could be associated with ongoing range expansion. Wing morphometric traits weakly correlated with environmental gradients within Europe, but wing size differed between invasive and source populations located in different climatic areas. Niche similarities between source and invasive ranges might have facilitated the establishment of populations. Nonetheless, we found evidence for environmental-induced adaptive changes after introduction. The ability to rapidly evolve observed in invasive populations (genetic shift) together with a large proportion of unfilled potential suitable areas (80%) pave the way to further spread of Ae. albopictus in Europe.

Impact of merging commercial breeding lines on the genetic diversity of Landrace pigs

BACKGROUND

The pig breeding industry has undergone a large number of mergers in the past decades. Various commercial lines were merged or discontinued, which is expected to reduce the genetic diversity of the pig species. The objective of the current study was to investigate the genetic diversity of different former Dutch Landrace breeding lines and quantify their relationship with the current Dutch Landrace breed that originated from these lines.

RESULTS

Principal component analysis clearly divided the former Landrace lines into two main clusters, which are represented by Norwegian/Finnish Landrace lines and Dutch Landrace lines. Structure analysis revealed that each of the lines that are present in the Dutch Gene bank has a unique genetic identity. The current Dutch Landrace breed shows a high level of admixture and is closely related to the six former lines. The Dumeco N-line, which is conserved in the Dutch Gene bank, is poorly represented in the current Dutch Landrace. All seven lines (the six former and the current line) contribute almost equally to the genetic diversity of the Dutch Landrace breed. As expected, the current Dutch Landrace breed comprises only a small proportion of unique genetic diversity that was not present in the other lines. The genetic diversity level, as measured by Eding's core set method, was equal to 0.89 for the current Dutch Landrace breed, whereas total genetic diversity across the seven lines, measured by the same method, was equal to 0.99.

CONCLUSIONS

The current Dutch Landrace breed shows a high level of admixture and is closely related to the six former Dutch Landrace lines. Merging of commercial Landrace lines has reduced the genetic diversity of the Landrace population in the Netherlands, although a large proportion of the original variation is maintained. Thus, our recommendation is to conserve breeding lines in a gene bank before they are merged.

The active spread of adaptive variation for reef resilience

The speed at which species adapt depends partly on the rates of beneficial adaptation generation and how quickly they spread within and among populations. Natural rates of adaptation of corals may not be able to keep pace with climate warming. Several interventions have been proposed to fast-track thermal adaptation, including the intentional translocation of warm-adapted adults or their offspring (assisted gene flow, AGF) and the ex situ crossing of warm-adapted corals with conspecifics from cooler reefs (hybridization or selective breeding) and field deployment of those offspring. The introgression of temperature tolerance loci into the genomic background of cooler-environment corals aims to facilitate adaptation to warming while maintaining fitness under local conditions. Here we use research on selective sweeps and connectivity to understand the spread of adaptive variants as it applies to AGF on the Great Barrier Reef (GBR), focusing on the genus Acropora. Using larval biophysical dispersal modeling, we estimate levels of natural connectivity in warm-adapted northern corals. We then model the spread of adaptive variants from single and multiple reefs and assess if the natural and assisted spread of adaptive variants will occur fast enough to prepare receiving central and southern populations given current rates of warming. We also estimate fixation rates and spatial extent of fixation under multiple release scenarios to inform intervention design. Our results suggest that thermal tolerance is unlikely to spread beyond northern reefs to the central and southern GBR without intervention, and if it does, 30+ generations are needed for adaptive gene variants to reach fixation even under multiple release scenarios. We argue that if translocation, breeding, and reseeding risks are managed, AGF using multiple release reefs can be beneficial for the restoration of coral populations. These interventions should be considered in addition to conventional management and accompanied by strong mitigation of CO2 emissions.

Population Genomic Signatures of Genetic Structure and Environmental Selection in the Catadromous Roughskin Sculpin Trachidermus fasciatus

Understanding the patterns of genetic diversity and adaptation across species' range is crucial to assess its long-term persistence and determine appropriate conservation measures. The impacts of human activities on the genetic diversity and genetic adaptation to heterogeneous environments remain poorly understood in the marine realm. The roughskin sculpin (Trachidermus fasciatus) is a small catadromous fish, and has been listed as a second-class state protected aquatic animal since 1988 in China. To elucidate the underlying mechanism of population genetic structuring and genetic adaptations to local environments, RAD tags were sequenced for 202 individuals in nine populations across the range of T. fasciatus in China. The pairwise FST values over 9,271 filtered SNPs were significant except that between Dongying and Weifang. All the genetic clustering analysis revealed significant population structure with high support for eight distinct genetic clusters. Both the minor allele frequency spectra and Ne estimations suggested extremely small Ne in some populations (e.g., Qinhuangdao, Rongcheng, Wendeng, and Qingdao), which might result from recent population bottleneck. The strong genetic structure can be partly attributed to genetic drift and habitat fragmentation, likely due to the anthropogenic activities. Annotations of candidate adaptive loci suggested that genes involved in metabolism, development, and osmoregulation were critical for adaptation to spatially heterogenous environment of local populations. In the context of anthropogenic activities and environmental change, results of the present population genomic work provided important contributions to the understanding of genetic differentiation and adaptation to changing environments.

What can be learned by scanning the genome for molecular convergence in wild populations?

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.

Landscape genomics to the rescue of a tropical bee threatened by habitat loss and climate change

Habitat degradation and climate change are currently threatening wild pollinators, compromising their ability to provide pollination services to wild and cultivated plants. Landscape genomics offers powerful tools to assess the influence of landscape modifications on genetic diversity and functional connectivity, and to identify adaptations to local environmental conditions that could facilitate future bee survival. Here, we assessed range-wide patterns of genetic structure, genetic diversity, gene flow, and local adaptation in the stingless bee Melipona subnitida, a tropical pollinator of key biological and economic importance inhabiting one of the driest and hottest regions of South America. Our results reveal four genetic clusters across the species' full distribution range. All populations were found to be under a mutation-drift equilibrium, and genetic diversity was not influenced by the amount of reminiscent natural habitats. However, genetic relatedness was spatially autocorrelated and isolation by landscape resistance explained range-wide relatedness patterns better than isolation by geographic distance, contradicting earlier findings for stingless bees. Specifically, gene flow was enhanced by increased thermal stability, higher forest cover, lower elevations, and less corrugated terrains. Finally, we detected genomic signatures of adaptation to temperature, precipitation, and forest cover, spatially distributed in latitudinal and altitudinal patterns. Taken together, our findings shed important light on the life history of M. subnitida and highlight the role of regions with large thermal fluctuations, deforested areas, and mountain ranges as dispersal barriers. Conservation actions such as restricting long-distance colony transportation, preserving local adaptations, and improving the connectivity between highlands and lowlands are likely to assure future pollination services.

Standing genomic variation within coding and regulatory regions contributes to the adaptive capacity to climate in a foundation tree species

Global climate is rapidly changing, and the ability for tree species to adapt is dependent on standing genomic variation; however, the distribution and abundance of functional and adaptive variants are poorly understood in natural systems. We test key hypotheses regarding the genetics of adaptive variation in a foundation tree: genomic variation is associated with climate, and genomic variation is more likely to be associated with temperature than precipitation or aridity. To test these hypotheses, we used 9,593 independent, genomic single-nucleotide polymorphisms (SNPs) from 270 individuals sampled from Corymbia calophylla's entire distribution in south-western Western Australia, spanning orthogonal temperature and precipitation gradients. Environmental association analyses returned 537 unique SNPs putatively adaptive to climate. We identified SNPs associated with climatic variation (i.e., temperature [458], precipitation [75] and aridity [78]) across the landscape. Of these, 78 SNPs were nonsynonymous (NS), while 26 SNPs were found within gene regulatory regions. The NS and regulatory candidate SNPs associated with temperature explained more deviance (27.35%) than precipitation (5.93%) and aridity (4.77%), suggesting that temperature provides stronger adaptive signals than precipitation. Genes associated with adaptive variants include functions important in stress responses to temperature and precipitation. Patterns of allelic turnover of NS and regulatory SNPs show small patterns of change through climate space with the exception of an aldehyde dehydrogenase gene variant with 80% allelic turnover with temperature. Together, these findings provide evidence for the presence of adaptive variation to climate in a foundation species and provide critical information to guide adaptive management practices.

Identification of Arabis alpina genomic regions associated with climatic variables along an elevation gradient through whole genome scan

We performed a pooled whole-genome sequencing on samples of the alpine plant Arabis alpina, harvested in ten populations along an elevation gradient in the French Alps. A large dataset of genetic variations was produced as single nucleotide polymorphisms (SNPs). A combined genome scan approach enabled detecting genomic regions associated with a synthetic environmental variable characterizing the climate at each sampling location. Positive loci detected by two methods were retained and belong to 19 regions in the Arabis alpina genome. The most significant region harbors an ortholog of the AtNAC062 gene, encoding a membrane-bound transcription factor described as linking the cold response and pathogen resistance that may confer protection to plants under extended snow coverage at high elevations. Other genes involved in the stress response or in flowering regulation were also detected. Altogether, our results indicated that Arabis alpina represent a suitable model for studying genomic adaptation in alpine perennial plants.

The Effects of Epistasis and Pleiotropy on Genome-Wide Scans for Adaptive Outlier Loci

With the advent of next-generation sequencing approaches, the search for individual loci underlying local adaptation has become a major enterprise in evolutionary biology. One promising method to identify such loci is to examine genome-wide patterns of differentiation, using an FST-outlier approach. The effects of pleiotropy and epistasis on this approach are not yet known. Here, we model 2 populations of a sexually reproducing, diploid organism with 2 quantitative traits, one of which is involved in local adaptation. We consider genetic architectures with and without pleiotropy and epistasis. We also model neutral marker loci on an explicit genetic map as the 2 populations diverge and apply FST outlier approaches to determine the extent to which quantitative trait loci (QTL) are detectable. Our results show, under a wide range of conditions, that only a small number of QTL are typically responsible for most of the trait divergence between populations, even when inheritance is highly polygenic. We find that the loci making the largest contributions to trait divergence tend to be detectable outliers. These loci also make the largest contributions to within-population genetic variance. The addition of pleiotropy reduces the extent to which quantitative traits can evolve independently but does not reduce the efficacy of outlier scans. The addition of epistasis, however, reduces the mean FST values for causative QTL, making these loci more difficult, but not impossible, to detect in outlier scans.

Stress response, behavior, and development are shaped by transposable element-induced mutations in Drosophila

Most of the current knowledge on the genetic basis of adaptive evolution is based on the analysis of single nucleotide polymorphisms (SNPs). Despite increasing evidence for their causal role, the contribution of structural variants to adaptive evolution remains largely unexplored. In this work, we analyzed the population frequencies of 1,615 Transposable Element (TE) insertions annotated in the reference genome of Drosophila melanogaster, in 91 samples from 60 worldwide natural populations. We identified a set of 300 polymorphic TEs that are present at high population frequencies, and located in genomic regions with high recombination rate, where the efficiency of natural selection is high. The age and the length of these 300 TEs are consistent with relatively young and long insertions reaching high frequencies due to the action of positive selection. Besides, we identified a set of 21 fixed TEs also likely to be adaptive. Indeed, we, and others, found evidence of selection for 84 of these reference TE insertions. The analysis of the genes located nearby these 84 candidate adaptive insertions suggested that the functional response to selection is related with the GO categories of response to stimulus, behavior, and development. We further showed that a subset of the candidate adaptive TEs affects expression of nearby genes, and five of them have already been linked to an ecologically relevant phenotypic effect. Our results provide a more complete understanding of the genetic variation and the fitness-related traits relevant for adaptive evolution. Similar studies should help uncover the importance of TE-induced adaptive mutations in other species as well.

Characterizing neutral and adaptive genomic differentiation in a changing climate: The most northerly freshwater fish as a model

Arctic freshwater ecosystems have been profoundly affected by climate change. Given that the Arctic charr (Salvelinus alpinus) is often the only fish species inhabiting these ecosystems, it represents a valuable model for studying the impacts of climate change on species life-history diversity and adaptability. Using a genotyping-by-sequencing approach, we identified 5,976 neutral single nucleotide polymorphisms and found evidence for reduced gene flow between allopatric morphs from two high Arctic lakes, Linne'vatn (Anadromous, Normal, and Dwarf) and Ellasjøen (Littoral and Pelagic). Within each lake, the degree of genetic differentiation ranged from low (Pelagic vs. Littoral) to moderate (Anadromous and Normal vs. Dwarf). We identified 17 highly diagnostic, putatively adaptive SNPs that differentiated the allopatric morphs. Although we found no evidence for adaptive differences between morphs within Ellasjøen, we found evidence for moderate (Anadromous vs. Normal) to high genetic differentiation (Anadromous and Normal vs. Dwarf) among morphs within Linne'vatn based on two adaptive loci. As these freshwater ecosystems become more productive, the frequency of sympatric morphs in Ellasjøen will likely shift based on foraging opportunities, whereas the propensity to migrate may decrease in Linne'vatn, increasing the frequency of the Normal morph. The Dwarf charr was the most genetically distinct group. Identifying the biological basis for small body size should elucidate the potential for increased growth and subsequent interbreeding with sympatric morphs. Overall, neutral and adaptive genomic differentiation between allopatric and some sympatric morphs suggests that the response of Arctic charr to climate change will be variable across freshwater ecosystems.

Genome-wide signals of drift and local adaptation during rapid lineage divergence in a songbird

The formation of independent evolutionary lineages involves neutral and selective factors, and understanding their relative roles in population divergence is a fundamental goal of speciation research. Correlations between allele frequencies and environmental variability can reveal the role of selection, yet the relative contribution of drift can be difficult to establish. Recently diversified taxa like the Oregon junco (Aves, Passerellidae, Junco hyemalis oreganus) of western North America provide ideal scenarios to apply genetic-environment association analyses (GEA) while controlling for population structure. Analysis of genome-wide SNP loci revealed marked genetic structure consisting of differentiated populations in isolated, dry southern mountain ranges, and less divergent, recently expanded populations in humid northern latitudes. We used correlations between genomic and environmental variance to test for three specific modes of evolutionary divergence: (a) drift in geographic isolation, (b) differentiation along continuous selective gradients and (c) isolation-by-adaptation. We found evidence of strong drift in southern mountains, but also signals of local adaptation driven by temperature, precipitation, elevation and vegetation, especially when controlling for population history. We identified numerous variants under selection scattered across the genome, suggesting that local adaptation can promote rapid differentiation when acting over multiple independent loci.

Local adaptation in mainland anole lizards: Integrating population history and genome-environment associations

Environmental gradients constrain physiological performance and thus species' ranges, suggesting that species occurrence in diverse environments may be associated with local adaptation. Genome-environment association analyses (GEAA) have become central for studies of local adaptation, yet they are sensitive to the spatial orientation of historical range expansions relative to landscape gradients. To test whether potentially adaptive genotypes occur in varied climates in wide-ranged species, we implemented GEAA on the basis of genomewide data from the anole lizards Anolis ortonii and Anolis punctatus, which expanded from Amazonia, presently dominated by warm and wet settings, into the cooler and less rainy Atlantic Forest. To examine whether local adaptation has been constrained by population structure and history, we estimated effective population sizes, divergence times, and gene flow under a coalescent framework. In both species, divergence between Amazonian and Atlantic Forest populations dates back to the mid-Pleistocene, with subsequent gene flow. We recovered eleven candidate genes involved with metabolism, immunity, development, and cell signaling in A. punctatus and found no loci whose frequency is associated with environmental gradients in A. ortonii. Distinct signatures of adaptation between these species are not associated with historical constraints or distinct climatic space occupancies. Similar patterns of spatial structure between selected and neutral SNPs along the climatic gradient, as supported by patterns of genetic clustering in A. punctatus, may have led to conservative GEAA performance. This study illustrates how tests of local adaptation can benefit from knowledge about species histories to support hypothesis formulation, sampling design, and landscape gradient characterization.

Comparison of Migratory and Resident Populations of Brown Trout Reveals Candidate Genes for Migration Tendency

Candidate genes associated with migration have been identified in multiple taxa: including salmonids, many of whom perform migrations requiring a series of physiological changes associated with the freshwater–saltwater transition. We screened over 5,500 SNPs for signatures of selection related to migratory behavior of brown trout Salmo trutta by focusing on ten differentially migrating freshwater populations from two watersheds (the Koutajoki and the Oulujoki). We found eight outlier SNPs potentially associated with migratory versus resident life history using multiple (≥3) outlier detection approaches. Comparison of three migratory versus resident population pairs in the Koutajoki watershed revealed seven outlier SNPs, of which three mapped close to genes ZNF665-like, GRM4-like, and PCDH8-like that have been previously associated with migration and smoltification in salmonids. Two outlier SNPs mapped to genes involved in mucus secretion (ST3GAL1-like) and osmoregulation (C14orf37-like). The last two strongly supported outlier SNPs mapped to thermally induced genes (FNTA1-like, FAM134C-like). Within the Oulujoki, the only consistent outlier SNP mapped close to a gene (EZH2) that is associated with compensatory growth in fasted trout. Our results suggest that a relatively small yet common set of genes responsible for physiological functions associated with resident and migratory life histories is evolutionarily conserved.

Molecular signatures of divergence and selection in closely related pine taxa

Efforts to detect loci under selection in plants have mostly focussed on single species. However, assuming that intraspecific divergence may lead to speciation, comparisons of genetic variation within and among recently diverged taxa can help to locate such genes. In this study, coalescent and outlier detection methods were used to assess nucleotide polymorphism and divergence at 79 nuclear gene fragments (1212 SNPs) in 16 populations (153 individuals) of the closely related, but phenotypically and ecologically distinct, pine taxa Pinus mugo, P. uliginosa and P. uncinata across their European distributions. Simultaneously, mitochondrial DNA markers, which are maternally inherited in pines and distributed by seeds at short geographic distance, were used to assess genetic relationships of the focal populations and taxa. The majority of nuclear loci showed homogenous patterns of variation between the taxa due to a high number of shared SNPs and haplotypes, similar levels of polymorphism, and low net divergence. However, against this common genetic background and an overall low population structure within taxa at mitochondrial markers, we identified several genes showing signatures of selection, accompanied by significant intra- and interspecific divergence. Our results indicate that loci involved in species divergence may be involved in intraspecific local adaptation.

Identifying genomic hotspots of differentiation and candidate genes involved in the adaptive divergence of pea aphid host races

Identifying the genomic bases of adaptation to novel environments is a long-term objective in evolutionary biology. Because genetic differentiation is expected to increase between locally adapted populations at the genes targeted by selection, scanning the genome for elevated levels of differentiation is a first step towards deciphering the genomic architecture underlying adaptive divergence. The pea aphid Acyrthosiphon pisum is a model of choice to address this question, as it forms a large complex of plant-specialized races and cryptic species, resulting from recent adaptive radiation. Here, we characterized genomewide polymorphisms in three pea aphid races specialized on alfalfa, clover and pea crops, respectively, which we sequenced in pools (poolseq). Using a model-based approach that explicitly accounts for selection, we identified 392 genomic hotspots of differentiation spanning 47.3 Mb and 2,484 genes (respectively, 9.12% of the genome size and 8.10% of its genes). Most of these highly differentiated regions were located on the autosomes, and overall differentiation was weaker on the X chromosome. Within these hotspots, high levels of absolute divergence between races suggest that these regions experienced less gene flow than the rest of the genome, most likely by contributing to reproductive isolation. Moreover, population-specific analyses showed evidence of selection in every host race, depending on the hotspot considered. These hotspots were significantly enriched for candidate gene categories that control host-plant selection and use. These genes encode 48 salivary proteins, 14 gustatory receptors, 10 odorant receptors, five P450 cytochromes and one chemosensory protein, which represent promising candidates for the genetic basis of host-plant specialization and ecological isolation in the pea aphid complex. Altogether, our findings open new research directions towards functional studies, for validating the role of these genes on adaptive phenotypes.