The Evolutionary History of Nebraska Deer Mice: Local Adaptation in the Face of Strong Gene Flow

Genetic mechanisms of animal camouflage: an interdisciplinary perspective

Camouflage is a classic example of a trait wherein animals respond to natural selection to avoid predation or attract prey. This unique phenomenon has attracted significant recent attention and the rapid development of integrative research methods is facilitating advances in our understanding of the in-depth genetic mechanisms of camouflage. In this review article, we revisit camouflage definitions and strategies and then we examine the underlying mechanisms of the two most common forms of camouflage, crypsis and masquerade, that have recently been elucidated using multiple approaches. We also discuss unresolved questions related to camouflage. Ultimately, we highlight the implications of camouflage for informing various key issues in ecology and evolution.

Demographic history and genomic signatures of selection in a widespread vertebrate ectotherm

Environmental conditions vary greatly across large geographic ranges, and yet certain species inhabit entire continents. In such species, genomic sequencing can inform our understanding of colonization history and the impact of selection on the genome as populations experience diverse local environments. As ectothermic vertebrates are among the most vulnerable to environmental change, it is critical to understand the contributions of local adaptation to population survival. Widespread ectotherms offer an opportunity to explore how species can successfully inhabit such differing environments and how future climatic shifts will impact species' survival. In this study, we investigated the widespread painted turtle (Chrysemys picta) to assess population genomic structure, demographic history, and genomic signatures of selection in the western extent of the range. We found support for a substantial role of serial founder effects in shaping population genomic structure: demographic analysis and runs of homozygosity were consistent with bottlenecks of increasing severity from eastern to western populations during and following the Last Glacial Maximum, and edge populations were more strongly diverged and had less genetic diversity than those from the centre of the range. We also detected outlier loci, but allelic patterns in many loci could be explained by either genetic surfing or selection. While range expansion complicates the identification of loci under selection, we provide candidates for future study of local adaptation in a long-lived, widespread ectotherm that faces an uncertain future as the global climate continues to rapidly change.

Local adaptation with gene flow in a highly dispersive shark

Adaptive divergence in response to environmental clines are expected to be common in species occupying heterogeneous environments. Despite numerous advances in techniques appropriate for non-model species, gene-environment association studies in elasmobranchs are still scarce. The bronze whaler or copper shark (Carcharhinus brachyurus) is a large coastal shark with a wide distribution and one of the most exploited elasmobranchs in southern Africa. Here, we assessed the distribution of neutral and adaptive genomic diversity in C. brachyurus across a highly heterogeneous environment in southern Africa based on genome-wide SNPs obtained through a restriction site-associated DNA method (3RAD). A combination of differentiation-based genome-scan (outflank) and genotype-environment analyses (redundancy analysis, latent factor mixed models) identified a total of 234 differentiation-based outlier and candidate SNPs associated with bioclimatic variables. Analysis of 26,299 putatively neutral SNPs revealed moderate and evenly distributed levels of genomic diversity across sites from the east coast of South Africa to Angola. Multivariate and clustering analyses demonstrated a high degree of gene flow with no significant population structuring among or within ocean basins. In contrast, the putatively adaptive SNPs demonstrated the presence of two clusters and deep divergence between Angola and all other individuals from Namibia and South Africa. These results provide evidence for adaptive divergence in response to a heterogeneous seascape in a large, mobile shark despite high levels of gene flow. These results are expected to inform management strategies and policy at the national and regional level for conservation of C. brachyurus populations.

Genomic signals of local adaptation across climatically heterogenous habitats in an invasive tropical fruit fly (Bactrocera tryoni)

Local adaptation plays a key role in the successful establishment of pest populations in new environments by enabling them to tolerate novel biotic and abiotic conditions experienced outside their native range. However, the genomic underpinnings of such adaptive responses remain unclear, especially for agriculturally important pests. We investigated population genomic signatures in the tropical/subtropical Queensland fruit fly, Bactrocera tryoni, which has an expanded range encompassing temperate and arid zones in Australia, and tropical zones in the Pacific Islands. Using reduced representation sequencing data from 28 populations, we detected allele frequency shifts associated with the native/invasive status of populations and identified environmental factors that have likely driven population differentiation. We also determined that precipitation, temperature, and geographic variables explain allelic shifts across the distribution range of B. tryoni. We found spatial heterogeneity in signatures of local adaptation across various climatic conditions in invaded areas. Specifically, disjunct invasive populations in the tropical Pacific Islands and arid zones of Australia were characterised by multiple significantly differentiated single nucleotide polymorphisms (SNPs), some of which were associated with genes with well-understood function in environmental stress (e.g., heat and desiccation) response. However, invasive populations in southeast Australian temperate zones showed higher gene flow with the native range and lacked a strong local adaptive signal. These results suggest that population connectivity with the native range has differentially affected local adaptive patterns in different invasive populations. Overall, our findings provide insights into the evolutionary underpinnings of invasion success of an important horticultural pest in climatically distinct environments.

Candidate genes under selection in song sparrows co-vary with climate and body mass in support of Bergmann's Rule

Ecogeographic rules denote spatial patterns in phenotype and environment that may reflect local adaptation as well as a species' capacity to adapt to change. To identify genes underlying Bergmann's Rule, which posits that spatial correlations of body mass and temperature reflect natural selection and local adaptation in endotherms, we compare 79 genomes from nine song sparrow (Melospiza melodia) subspecies that vary ~300% in body mass (17 - 50 g). Comparing large- and smaller-bodied subspecies revealed 9 candidate genes in three genomic regions associated with body mass. Further comparisons to the five smallest subspecies endemic to California revealed eight SNPs within four of the candidate genes (GARNL3, RALGPS1, ANGPTL2, and COL15A1) associated with body mass and varying as predicted by Bergmann's Rule. Our results support the hypothesis that co-variation in environment, body mass and genotype reflect the influence of natural selection on local adaptation and a capacity for contemporary evolution in this diverse species.

Rapid evolution of novel biotic interactions in the UK Brown Argus butterfly uses genomic variation from across its geographical range

Understanding the rate and extent to which populations can adapt to novel environments at their ecological margins is fundamental to predicting the persistence of biological communities during ongoing and rapid global change. Recent range expansion in response to climate change in the UK butterfly Aricia agestis is associated with the evolution of novel interactions with a larval food plant, and the loss of its ability to use an ancestral host species. Using ddRAD analysis of 61,210 variable SNPs from 261 females from throughout the UK range of this species, we identify genomic regions at multiple chromosomes that are associated with evolutionary responses, and their association with demographic history and ecological variation. Gene flow appears widespread throughout the range, despite the apparently fragmented nature of the habitats used by this species. Patterns of haplotype variation between selected and neutral genomic regions suggest that evolution associated with climate adaptation is polygenic, resulting from the independent spread of alleles throughout the established range of this species, rather than the colonization of pre-adapted genotypes from coastal populations. These data suggest that rapid responses to climate change do not depend on the availability of pre-adapted genotypes. Instead, the evolution of novel forms of biotic interaction in A. agestis has occurred during range expansion, through the assembly of novel genotypes from alleles from multiple localities.

An evolutionarily distinct ringed seal in the Ilulissat Icefjord

The Earth's polar regions are low rates of inter- and intraspecific diversification. An extreme mammalian example is the Arctic ringed seal (Pusa hispida hispida), which is assumed to be panmictic across its circumpolar Arctic range. Yet, local Inuit communities in Greenland and Canada recognize several regional variants; a finding supported by scientific studies of body size variation. It is however unclear whether this phenotypic variation reflects plasticity, morphs or distinct ecotypes. Here, we combine genomic, biologging and survey data, to document the existence of a unique ringed seal ecotype in the Ilulissat Icefjord (locally 'Kangia'), Greenland; a UNESCO World Heritage site, which is home to the most productive marine-terminating glacier in the Arctic. Genomic analyses reveal a divergence of Kangia ringed seals from other Arctic ringed seals about 240 kya, followed by secondary contact since the Last Glacial Maximum. Despite ongoing gene flow, multiple genomic regions appear under strong selection in Kangia ringed seals, including candidate genes associated with pelage coloration, growth and osmoregulation, potentially explaining the Kangia seal's phenotypic and behavioural uniqueness. The description of 'hidden' diversity and adaptations in yet another Arctic species merits a reassessment of the evolutionary processes that have shaped Arctic diversity and the traditional view of this region as an evolutionary freezer. Our study highlights the value of indigenous knowledge in guiding science and calls for efforts to identify distinct populations or ecotypes to understand how these might respond differently to environmental change.

Analysis of the genetic loci of pigment pattern evolution in vertebrates

Vertebrate pigmentation patterns are amongst the best characterised model systems for studying the genetic basis of adaptive evolution. The wealth of available data on the genetic basis for pigmentation evolution allows for analysis of trends and quantitative testing of evolutionary hypotheses. We employed Gephebase, a database of genetic variants associated with natural and domesticated trait variation, to examine trends in how cis-regulatory and coding mutations contribute to vertebrate pigmentation phenotypes, as well as factors that favour one mutation type over the other. We found that studies with lower ascertainment bias identified higher proportions of cis-regulatory mutations, and that cis-regulatory mutations were more common amongst animals harbouring a higher number of pigment cell classes. We classified pigmentation traits firstly according to their physiological basis and secondly according to whether they affect colour or pattern, and identified that carotenoid-based pigmentation and variation in pattern boundaries are preferentially associated with cis-regulatory change. We also classified genes according to their developmental, cellular, and molecular functions. We found a greater proportion of cis-regulatory mutations in genes implicated in upstream developmental processes compared to those involved in downstream cellular functions, and that ligands were associated with a higher proportion of cis-regulatory mutations than their respective receptors. Based on these trends, we discuss future directions for research in vertebrate pigmentation evolution.

Genomic architecture of migration timing in a long-distance migratory songbird

The impact of climate change on spring phenology poses risks to migratory birds, as migration timing is controlled predominantly by endogenous mechanisms. Despite recent advances in our understanding of the underlying genetic basis of migration timing, the ways that migration timing phenotypes in wild individuals may map to specific genomic regions requires further investigation. We examined the genetic architecture of migration timing in a long-distance migratory songbird (purple martin, Progne subis subis) by integrating genomic data with an extensive dataset of direct migratory tracks. A moderate to large amount of variance in spring migration arrival timing was explained by genomics (proportion of phenotypic variation explained by genomics = 0.74; polygenic score R² = 0.24). On chromosome 1, a region that was differentiated between migration timing phenotypes contained genes that could facilitate nocturnal flights and act as epigenetic modifiers. Overall, these results advance our understanding of the genomic underpinnings of migration timing.

Life on a beach leads to phenotypic divergence despite gene flow for an island lizard

Limited spatial separation within small islands suggests that observed population divergence may occur due to habitat differences without interruption to gene flow but strong evidence of this is scarce. The wall lizard Teira dugesii lives in starkly contrasting shingle beach and inland habitats on the island of Madeira. We used a matched pairs sampling design to examine morphological and genomic divergence between four beach and adjacent (<1 km) inland areas. Beach populations are significantly darker than corresponding inland populations. Geometric morphometric analyses reveal divergence in head morphology: beach lizards have generally wider snouts. Genotyping-by-sequencing allows the rejection of the hypothesis that beach populations form a distinct lineage. Bayesian analyses provide strong support for models that incorporate gene flow, relative to those that do not, replicated at all pairs of matched sites. Madeiran lizards show morphological divergence between habitats in the face of gene flow, revealing how divergence may originate within small islands.

Rapid and predictable genome evolution across three hybrid ant populations

Hybridization is frequent in the wild but it is unclear when admixture events lead to predictable outcomes and if so, at what timescale. We show that selection led to correlated sorting of genetic variation rapidly after admixture in 3 hybrid Formica aquilonia × F. polyctena ant populations. Removal of ancestry from the species with the lowest effective population size happened in all populations, consistent with purging of deleterious load. This process was modulated by recombination rate variation and the density of functional sites. Moreover, haplotypes with signatures of positive selection in either species were more likely to fix in hybrids. These mechanisms led to mosaic genomes with comparable ancestry proportions. Our work demonstrates predictable evolution over short timescales after admixture in nature.

Local adaptation: Causal agents of selection and adaptive trait divergence

Divergent selection across the landscape can favor the evolution of local adaptation in populations experiencing contrasting conditions. Local adaptation is widely observed in a diversity of taxa, yet we have a surprisingly limited understanding of the mechanisms that give rise to it. For instance, few have experimentally confirmed the biotic and abiotic variables that promote local adaptation, and fewer yet have identified the phenotypic targets of selection that mediate local adaptation. Here, we highlight critical gaps in our understanding of the process of local adaptation and discuss insights emerging from in-depth investigations of the agents of selection that drive local adaptation, the phenotypes they target, and the genetic basis of these phenotypes. We review historical and contemporary methods for assessing local adaptation, explore whether local adaptation manifests differently across life history, and evaluate constraints on local adaptation.

Genetic basis of growth reaction to drought stress differs in contrasting high-latitude treeline ecotones of a widespread conifer

Climate change is increasing the frequency and intensity of drought events in many boreal forests. Trees are sessile organisms with a long generation time, which makes them vulnerable to fast climate change and hinders fast adaptations. Therefore, it is important to know how forests cope with drought stress and to explore the genetic basis of these reactions. We investigated three natural populations of white spruce (Picea glauca) in Alaska, located at one drought-limited and two cold-limited treelines with a paired plot design of one forest and one treeline plot. We obtained individual increment cores from 458 trees and climate data to assess dendrophenotypes, in particular the growth reaction to drought stress. To explore the genetic basis of these dendrophenotypes, we genotyped the individual trees at 3000 single nucleotide polymorphisms in candidate genes and performed genotype-phenotype association analysis using linear mixed models and Bayesian sparse linear mixed models. Growth reaction to drought stress differed in contrasting treeline populations. Therefore, the populations are likely to be unevenly affected by climate change. We identified 40 genes associated with dendrophenotypic traits that differed among the treeline populations. Most genes were identified in the drought-limited site, indicating comparatively strong selection pressure of drought-tolerant phenotypes. Contrasting patterns of drought-associated genes among sampled sites and in comparison to Canadian populations in a previous study suggest that drought adaptation acts on a local scale. Our results highlight genes that are associated with wood traits which in turn are critical for the establishment and persistence of future forests under climate change.

Evolutionary Genomics of a Subdivided Species

The ways in which genetic variation is distributed within and among populations is a key determinant of the evolutionary features of a species. However, most comprehensive studies of these features have been restricted to studies of subdivision in settings known to have been driven by local adaptation, leaving our understanding of the natural dispersion of allelic variation less than ideal. Here, we present a geographic population-genomic analysis of 10 populations of the freshwater microcrustacean Daphnia pulex, an emerging model system in evolutionary genomics. These populations exhibit a pattern of moderate isolation-by-distance, with an average migration rate of 0.6 individuals per generation, and average effective population sizes of ∼650,000 individuals. Most populations contain numerous private alleles, and genomic scans highlight the presence of islands of excessively high population subdivision for more common alleles. A large fraction of such islands of population divergence likely reflect historical neutral changes, including rare stochastic migration and hybridization events. The data do point to local adaptive divergence, although the precise nature of the relevant variation is diffuse and cannot be associated with particular loci, despite the very large sample sizes involved in this study. In contrast, an analysis of between-species divergence highlights positive selection operating on a large set of genes with functions nearly nonoverlapping with those involved in local adaptation, in particular ribosome structure, mitochondrial bioenergetics, light reception and response, detoxification, and gene regulation. These results set the stage for using D. pulex as a model for understanding the relationship between molecular and cellular evolution in the context of natural environments.

An enhancer of Agouti contributes to parallel evolution of cryptically colored beach mice

Identifying the genetic basis of repeatedly evolved traits provides a way to reconstruct their evolutionary history and ultimately investigate the predictability of evolution. Here, we focus on the oldfield mouse (Peromyscus polionotus), which occurs in the southeastern United States, where it exhibits considerable color variation. Dorsal coats range from dark brown in mainland mice to near white in mice inhabiting sandy beaches; this light pelage has evolved independently on Florida's Gulf and Atlantic coasts as camouflage from predators. To facilitate genomic analyses, we first generated a chromosome-level genome assembly of Peromyscus polionotus subgriseus. Next, in a uniquely variable mainland population (Peromyscus polionotus albifrons), we scored 23 pigment traits and performed targeted resequencing in 168 mice. We find that pigment variation is strongly associated with an ∼2-kb region ∼5 kb upstream of the Agouti signaling protein coding region. Using a reporter-gene assay, we demonstrate that this regulatory region contains an enhancer that drives expression in the dermis of mouse embryos during the establishment of pigment prepatterns. Moreover, extended tracts of homozygosity in this Agouti region indicate that the light allele experienced recent and strong positive selection. Notably, this same light allele appears fixed in both Gulf and Atlantic coast beach mice, despite these populations being separated by >1,000 km. Together, our results suggest that this identified Agouti enhancer allele has been maintained in mainland populations as standing genetic variation and from there, has spread to and been selected in two independent beach mouse lineages, thereby facilitating their rapid and parallel evolution.

Whole-genome analysis of multiple wood ant population pairs supports similar speciation histories, but different degrees of gene flow, across their European ranges

The application of demographic history modelling and inference to the study of divergence between species has become a cornerstone of speciation genomics. Speciation histories are usually reconstructed by analysing single populations from each species, assuming that the inferred population history represents the actual speciation history. However, this assumption may not be met when species diverge with gene flow, for example, when secondary contact may be confined to specific geographic regions. Here, we tested whether divergence histories inferred from heterospecific populations may vary depending on their geographic locations, using the two wood ant species Formica polyctena and F. aquilonia. We performed whole-genome resequencing of 20 individuals sampled in multiple locations across the European ranges of both species. Then, we reconstructed the histories of distinct heterospecific population pairs using a coalescent-based approach. Our analyses always supported a scenario of divergence with gene flow, suggesting that divergence started in the Pleistocene (c. 500 kya) and occurred with continuous asymmetrical gene flow from F. aquilonia to F. polyctena until a recent time, when migration became negligible (2-19 kya). However, we found support for contemporary gene flow in a sympatric pair from Finland, where the species hybridise, but no signature of recent bidirectional gene flow elsewhere. Overall, our results suggest that divergence histories reconstructed from a few individuals may be applicable at the species level. Nonetheless, the geographical context of populations chosen to represent their species should be taken into account, as it may affect estimates of migration rates between species when gene flow is spatially heterogeneous.

The response of a metapopulation to a changing environment

A species distributed across diverse environments may adapt to local conditions. We ask how quickly such a species changes its range in response to changed conditions. Szép et al. (Szép E, Sachdeva H, Barton NH. 2021 Polygenic local adaptation in metapopulations: a stochastic eco-evolutionary model. Evolution75, 1030-1045 (doi:10.1111/evo.14210)) used the infinite island model to find the stationary distribution of allele frequencies and deme sizes. We extend this to find how a metapopulation responds to changes in carrying capacity, selection strength, or migration rate when deme sizes are fixed. We further develop a 'fixed-state' approximation. Under this approximation, polymorphism is only possible for a narrow range of habitat proportions when selection is weak compared to drift, but for a much wider range otherwise. When rates of selection or migration relative to drift change in a single deme of the metapopulation, the population takes a time of order m-1 to reach the new equilibrium. However, even with many loci, there can be substantial fluctuations in net adaptation, because at each locus, alleles randomly get lost or fixed. Thus, in a finite metapopulation, variation may gradually be lost by chance, even if it would persist in an infinite metapopulation. When conditions change across the whole metapopulation, there can be rapid change, which is predicted well by the fixed-state approximation. This work helps towards an understanding of how metapopulations extend their range across diverse environments. This article is part of the theme issue 'Species' ranges in the face of changing environments (Part II)'.

Habitat-linked genetic variation supports microgeographic adaptive divergence in an island-endemic bird species

We investigated the potential mechanisms driving habitat-linked genetic divergence within a bird species endemic to a single 250 km² island. The island scrub-jay (Aphelocoma insularis) exhibits microgeographic divergence in bill morphology across pine-oak ecotones on Santa Cruz Island, California (USA) similar to adaptive differences described in mainland congeners over much larger geographic scales. To test whether individuals exhibit genetic differentiation related to habitat type and divergence in bill length, we genotyped over 3,000 single nucleotide polymorphisms (SNPs) in 123 adult island scrub-jay males from across Santa Cruz Island using restriction site-associated DNA sequencing (RADseq). Neutral landscape genomic analyses revealed that genome-wide genetic differentiation was primarily related to geographic distance and differences in habitat composition. We also found 168 putatively adaptive loci associated with habitat type using multivariate redundancy analysis (RDA) while controlling for spatial effects. Finally, two genome-wide association analyses revealed a polygenic basis to variation in bill length with multiple loci detected in or near genes known to affect bill morphology in other birds. Our findings support the hypothesis that divergent selection at microgeographic scales can cause adaptive divergence in the presence of ongoing gene flow.

Genetic isolation by distance underlies colour pattern divergence in red-eyed treefrogs (Agalychnis callidryas)

Investigating the spatial distribution of genetic and phenotypic variation can provide insights into the evolutionary processes that shape diversity in natural systems. We characterized patterns of genetic and phenotypic diversity to learn about drivers of colour-pattern diversification in red-eyed treefrogs (Agalychnis callidryas) in Costa Rica. Along the Pacific coast, red-eyed treefrogs have conspicuous leg colour patterning that transitions from orange in the north to purple in the south. We measured phenotypic variation of frogs, with increased sampling at sites where the orange-to-purple transition occurs. At the transition zone, we discovered the co-occurrence of multiple colour-pattern morphs. To explore possible causes of this variation, we generated a single nucleotide polymorphism data set to analyse population genetic structure, measure genetic diversity and infer the processes that mediate genotype-phenotype dynamics. We investigated how patterns of genetic relatedness correspond to individual measures of colour pattern along the coast, including testing for the role of hybridization in geographic regions where orange and purple phenotypic groups co-occur. We found no evidence that colour-pattern polymorphism in the transition zone arose through recent hybridization. Instead, a strong pattern of genetic isolation by distance indicates that colour-pattern variation was either retained through other processes such as ancestral colour polymorphisms or ancient secondary contact, or else it was generated by novel mutations. We found that phenotype changes along the Pacific coast more than would be expected based on genetic divergence and geographic distance alone. Combined, our results suggest the possibility of selective pressures acting on colour pattern at a small geographic scale.

Local adaptation to climate anomalies relates to species phylogeny

Climatic anomalies are increasing in intensity and frequency due to rapid rates of global change, leading to increased extinction risk for many species. The impacts of anomalies are likely to vary between species due to different degrees of sensitivity and extents of local adaptation. Here, we used long-term butterfly monitoring data of 143 species across six European bioclimatic regions to show how species’ population dynamics have responded to local or globally-calculated climatic anomalies, and how species attributes mediate these responses. Contrary to expectations, degree of apparent local adaptation, estimated from the relative population sensitivity to local versus global anomalies, showed no associations with species mobility or reproductive rate but did contain a strong phylogenetic signal. The existence of phylogenetically-patterned local adaptation to climate has important implications for forecasting species responses to current and future climatic conditions and for developing appropriate conservation practices.

Melero et al. investigate butterfly responses to climatic anomalies from long-term monitoring observations in the field. They found the degree of adaptation to local fluctuations in climate had a strong phylogenetic signal but was not associated with mobility or reproductive rate of a species.

Postglacial colonization in the Great Lakes Region by the white-footed mouse (Peromyscus leucopus): conflicts between genomic and field data

For regions that were covered by ice during the Pleistocene glaciations, species must have emigrated from unglaciated regions. However, it can be difficult to discern when and from what ancestral source populations such expansions took place, especially since warming climates introduce the possibility of very recent expansions. For example, in the Great Lakes region, pronounced climatic change includes past glaciations as well as recent, rapid warming. Here we evaluate different expansion hypotheses with a genomic study of the white-footed mouse (Peromyscus leucopus noveboracensis), which is one of the most common mammals throughout the Great Lakes region. Ecological surveys coupled with historical museum records suggest a recent range expansion of P. leucopus associated with the warming climate over the last decades. These detailed records have yet to be complemented by genomic data that provide the requisite resolution for detecting recent expansion, although some mitochondrial DNA (mtDNA) sequences have suggested possible hypotheses about the geography of expansion. With more than 7,000 loci generated using RADseq, we evaluate support for multiple hypotheses of a geographic expansion in the Upper Peninsula of Michigan (UP). Analysis of a single random single-nucleotide polymorphism per locus revealed a fine-scale population structure separating the Lower Peninsula (LP) population from all other populations in the UP. We also detected a genetic structure that reflects an evolutionary history of postglacial colonization from two different origins into the UP, one coming from the LP and one coming from the west. Instead of supporting a climate-driven range expansion, as suggested by field surveys, our results support more ancient postglacial colonization of the UP from two different ancestral sources. With these results, we offer new insights about P. leucopus geographic expansion history, as well as a more general phylogeographic framework for testing range shifts in the Great Lakes region.

Genetic guidelines for translocations: Maintaining intraspecific diversity in the lion (Panthera leo)

Conservation translocations have become an important management tool, particularly for large wildlife species such as the lion (Panthera leo). When planning translocations, the genetic background of populations needs to be taken into account; failure to do so risks disrupting existing patterns of genetic variation, ultimately leading to genetic homogenization, and thereby reducing resilience and adaptability of the species. We urge wildlife managers to include knowledge of the genetic background of source/target populations, as well as species-wide patterns, in any management intervention. We present a hierarchical decision-making tool in which we list 132 lion populations/lion conservation units and provide information on genetic assignment, uncertainty and suitability for translocation for each source/target combination. By including four levels of suitability, from 'first choice' to 'no option', we provide managers with a range of options. To illustrate the extent of international trade of lions, and the potential disruption of natural patterns of intraspecific diversity, we mined the CITES Trade Database for estimated trade quantities of live individuals imported into lion range states during the past 4 decades. We identified 1056 recorded individuals with a potential risk of interbreeding with wild lions, 772 being captive-sourced. Scoring each of the records with our decision-making tool illustrates that only 7% of the translocated individuals were 'first choice' and 73% were 'no option'. We acknowledge that other, nongenetic factors are important in the decision-making process, and hence a pragmatic approach is needed. A framework in which source/target populations are scored based on suitability is not only relevant to lion, but also to other species of wildlife that are frequently translocated. We hope that the presented overview supports managers to include genetics in future management decisions and contributes towards conservation of the lion in its full diversity.

Exome sequencing of deer mice on two California Channel Islands identifies potential adaptation to strongly contrasting ecological conditions

Understanding the forces that drive genotypic and phenotypic change in wild populations is a central goal of evolutionary biology. We examined exome variation in populations of deer mice from two of the California Channel Islands: Peromyscus maniculatus elusus from Santa Barbara Island and P. m. santacruzae from Santa Cruz Island exhibit significant differences in olfactory predator recognition, activity timing, aggressive behavior, morphology, prevalence of Sin Nombre virus, and population densities. We characterized variation in protein-coding regions using exome capture and sequencing of 25 mice from Santa Barbara Island and 22 mice from Santa Cruz Island. We identified and examined 386,256 SNPs using three complementary methods (BayeScan, pcadapt, and LFMM). We found strong differences in molecular variation between the two populations and 710 outlier SNPs in protein-coding genes that were detected by all three methods. We identified 35 candidate genes from this outlier set that were related to differences in phenotypes between island populations. Enrichment analyses demonstrated that patterns of molecular variation were associated with biological processes related to response to chemical stimuli and regulation of immune processes. Candidate genes associated with olfaction (Gfy, Tlr2, Vmn13r2, numerous olfactory receptor genes), circadian activity (Cry1), anxiety (Brca1), immunity (Cd28, Eif2ak4, Il12a, Syne1), aggression (Cyp19a, Lama2), and body size (Bc16, Syne1) exhibited non-synonymous mutations predicted to have moderate to large effects. Variation in olfaction-related genes, including a stop codon in the Santa Barbara Island population, suggests loss of predator-recognition traits at the molecular level, consistent with a lack of behavioral aversion to fox feces. These findings also suggest that divergent pathogen prevalence and population density may have influenced adaptive immunity and behavioral phenotypes, such as reduced aggression. Overall, our study indicates that ecological differences between islands are associated with signatures of selection in protein-coding genes underlying phenotypes that promote success in those environments.

Phenotypic and genomic differences between biomes of the South America marsh rat, Holochilus brasiliensis

Abiotic factors can influence genetic and phenotypic divergence in several ways, and identifying the mechanisms responsible for generating this variation is challenging. However, when evaluated in combination, ecological characteristics and genetic and phenotypic information can help us to understand how habitat preferences can influence morphological and genetic patterns exhibited by taxa distributed between distinct biomes, such as the Atlantic Forest and Pampas biomes in South America. By combining distributional, environmental, phenotypic and genomic information from a habitat-specialist semi-aquatic rodent (Holochilus brasiliensis), we quantified the relationship between ecological niche differences and the phenotypic and genetic variation. The results demonstrate notable segregation among the ecological niches of H. brasiliensis within each biome, although we could not refute the hypothesis of niche similarity or equivalency. Such differences are consistent with a solid morphometric variation associated with the size of these rodents. However, the ecological and morphometric differentiation is not accompanied by the same pattern of genetic variation. Despite differences in the connectivity patterns in both biomes, the genetic differences corroborate a consistent level of migration history between biomes. Additionally, the association tests show that the environment explains a small and non- significant part of the genetic variation but a significant portion of the morphometric variation.

Rapid Evolutionary Adaptation in Response to Selection on Quantitative Traits

Evolutionary adaptation after sudden environmental changes can occur very rapidly. The mechanisms facilitating rapid adaptation range from strong positive directional selection leading to large shifts in the allele frequencies at a few loci (selective sweeps) to polygenic selection causing small changes in allele frequencies at many loci. In addition, combinations of these two extreme mechanisms may also result in fast evolution. In recent years, following reports of new case studies of rapid adaptation, population genetic models have been proposed to explain these observations. In these models, the role of the major selective forces (positive directional and stabilizing selection) is highlighted as well as the genetic architecture of quantitative traits. Furthermore, the factors limiting the speed of adaptation are analyzed, in particular, the effects of random genetic drift and demography due to finite population size.

Using replicate hybrid zones to understand the genomic basis of adaptive divergence

Combining hybrid zone analysis with genomic data is a promising approach to understanding the genomic basis of adaptive divergence. It allows for the identification of genomic regions underlying barriers to gene flow. It also provides insights into spatial patterns of allele frequency change, informing about the interplay between environmental factors, dispersal and selection. However, when only a single hybrid zone is analysed, it is difficult to separate patterns generated by selection from those resulting from chance. Therefore, it is beneficial to look for repeatable patterns across replicate hybrid zones in the same system. We applied this approach to the marine snail Littorina saxatilis, which contains two ecotypes, adapted to wave-exposed rocks vs. high-predation boulder fields. The existence of numerous hybrid zones between ecotypes offered the opportunity to test for the repeatability of genomic architectures and spatial patterns of divergence. We sampled and phenotyped snails from seven replicate hybrid zones on the Swedish west coast and genotyped them for thousands of single nucleotide polymorphisms. Shell shape and size showed parallel clines across all zones. Many genomic regions showing steep clines and/or high differentiation were shared among hybrid zones, consistent with a common evolutionary history and extensive gene flow between zones, and supporting the importance of these regions for divergence. In particular, we found that several large putative inversions contribute to divergence in all locations. Additionally, we found evidence for consistent displacement of clines from the boulder-rock transition. Our results demonstrate patterns of spatial variation that would not be accessible without continuous spatial sampling, a large genomic data set and replicate hybrid zones.

Multigenerational backcrossing and introgression between two woodrat species at an abrupt ecological transition

When organisms experience secondary contact after allopatric divergence, genomic regions can introgress differentially depending on their relationships with adaptation, reproductive isolation, recombination, and drift. Analyses of genome-wide patterns of divergence and introgression could provide insight into the outcomes of hybridization and the potential relationship between allopatric divergence and reproductive isolation. Here, we generate population genetic data (26,262 SNPs; 353 individuals) using a reduced-representation sequencing approach to quantify patterns of ancestry, differentiation, and introgression between a pair of ecologically distinct mammals-the desert woodrat (N. lepida) and Bryant's woodrat (N. bryanti)-that hybridize at a sharp ecotone in southern California. Individual ancestry estimates confirmed that hybrids were rare in this bimodal hybrid zone, and entirely consisted of a few F₁ individuals and a broad range of multigenerational backcrosses. Genomic cline analyses indicated more than half of loci had elevated introgression from one genomic background into the other. However, introgression was not associated with relative or absolute measures of divergence, and loci with extreme values for both were not typically found near detoxification enzymes previously implicated in dietary specialization for woodrats. The decoupling of differentiation and introgression suggests that processes other than adaptation, such as drift, may underlie the extreme clines at this contact zone.

Genomic Adaptations to Salinity Resist Gene Flow in the Evolution of Floridian Watersnakes

The migration-selection balance often governs the evolution of lineages, and speciation with gene flow is now considered common across the tree of life. Ecological speciation is a process that can facilitate divergence despite gene flow due to strong selective pressures caused by ecological differences; however, the exact traits under selection are often unknown. The transition from freshwater to saltwater habitats provides strong selection targeting traits with osmoregulatory function. Several lineages of North American watersnakes (Nerodia spp.) are known to occur in saltwater habitat and represent a useful system for studying speciation by providing an opportunity to investigate gene flow and evaluate how species boundaries are maintained or degraded. We use double digest restriction-site associated DNA sequencing to characterize the migration-selection balance and test for evidence of ecological divergence within the Nerodia fasciata-clarkii complex in Florida. We find evidence of high intraspecific gene flow with a pattern of isolation-by-distance underlying subspecific lineages. However, we identify genetic structure indicative of reduced gene flow between inland and coastal lineages suggesting divergence due to isolation-by-environment. This pattern is consistent with observed environmental differences where the amount of admixture decreases with increased salinity. Furthermore, we identify significantly enriched terms related to osmoregulatory function among a set of candidate loci, including several genes that have been previously implicated in adaptation to salinity stress. Collectively, our results demonstrate that ecological differences, likely driven by salinity, cause strong divergent selection which promotes divergence in the N. fasciata-clarkii complex despite significant gene flow.

Adaptive Divergence under Gene Flow along an Environmental Gradient in Two Coexisting Stickleback Species

There is a general and solid theoretical framework to explain how the interplay between natural selection and gene flow affects local adaptation. Yet, to what extent coexisting closely related species evolve collectively or show distinctive evolutionary responses remains a fundamental question. To address this, we studied the population genetic structure and morphological differentiation of sympatric three-spined and nine-spined stickleback. We conducted genotyping-by-sequencing and morphological trait characterisation using 24 individuals of each species from four lowland brackish water (LBW), four lowland freshwater (LFW) and three upland freshwater (UFW) sites in Belgium and the Netherlands. This combination of sites allowed us to contrast populations from isolated but environmentally similar locations (LFW vs. UFW), isolated but environmentally heterogeneous locations (LBW vs. UFW), and well-connected but environmentally heterogenous locations (LBW vs. LFW). Overall, both species showed comparable levels of genetic diversity and neutral genetic differentiation. However, for all three spatial scales, signatures of morphological and genomic adaptive divergence were substantially stronger among populations of the three-spined stickleback than among populations of the nine-spined stickleback. Furthermore, most outlier SNPs in the two species were associated with local freshwater sites. The few outlier SNPs that were associated with the split between brackish water and freshwater populations were located on one linkage group in three-spined stickleback and two linkage groups in nine-spined stickleback. We conclude that while both species show congruent evolutionary and genomic patterns of divergent selection, both species differ in the magnitude of their response to selection regardless of the geographical and environmental context.

Connecting tree-ring phenotypes, genetic associations and transcriptomics to decipher the genomic architecture of drought adaptation in a widespread conifer

As boreal forests face significant threats from climate change, understanding evolutionary trajectories of coniferous species has become fundamental to adapting management and conservation to a drying climate. We examined the genomic architecture underlying adaptive variation related to drought tolerance in 43 populations of a widespread boreal conifer, white spruce (Piceaglauca [Moench] Voss), by combining genotype–environment associations, genotype–phenotype associations, and transcriptomics. Adaptive genetic variation was identified by correlating allele frequencies for 6,153 single nucleotide polymorphisms from 2,606 candidate genes with temperature, precipitation and aridity gradients, and testing for significant associations between genotypes and 11 dendrometric and drought‐related traits (i.e., anatomical, growth response and climate‐sensitivity traits) using a polygenic model. We identified a set of 285 genes significantly associated with a climatic factor or a phenotypic trait, including 110 that were differentially expressed in response to drought under greenhouse‐controlled conditions. The interlinked phenotype–genotype–environment network revealed eight high‐confidence genes involved in white spruce adaptation to drought, of which four were drought‐responsive in the expression analysis. Our findings represent a significant step toward the characterization of the genomic basis of drought tolerance and adaptation to climate in conifers, which is essential to enable the establishment of resilient forests in view of new climate conditions.

see also the Perspective by Lars Opgenoorth and Christian Rellstab

Environmental correlates of genetic variation in the invasive European starling in North America

Populations of invasive species that colonize and spread in novel environments may differentiate both through demographic processes and local selection. European starlings (Sturnus vulgaris) were introduced to New York in 1890 and subsequently spread throughout North America, becoming one of the most widespread and numerous bird species on the continent. Genome-wide comparisons across starling individuals and populations can identify demographic and/or selective factors that facilitated this rapid and successful expansion. We investigated patterns of genomic diversity and differentiation using reduced-representation genome sequencing of 17 winter-season sampling sites. Consistent with this species' high dispersal rate and rapid expansion history, we found low geographical differentiation and few F_ST outliers even at a continental scale. Despite starting from a founding population of ~180 individuals, North American starlings show only a moderate genetic bottleneck, and models suggest a dramatic increase in effective population size since introduction. In genotype-environment associations we found that ~200 single-nucleotide polymorphisms are correlated with temperature and/or precipitation against a background of negligible genome- and range-wide divergence. Given this evidence, we suggest that local adaptation in North American starlings may have evolved rapidly even in this wide-ranging and evolutionarily young system. This survey of genomic signatures of expansion in North American starlings is the most comprehensive to date and complements ongoing studies of world-wide local adaptation in these highly dispersive and invasive birds.

Evolutionary fine-tuning of background-matching camouflage among geographical populations in the sandy beach tiger beetle

Background-matching camouflage is a widespread adaptation in animals; however, few studies have thoroughly examined its evolutionary process and consequences. The tiger beetle Chaetodera laetescripta exhibits pronounced variation in elytral colour pattern among sandy habitats of different colour in the Japanese Archipelago. In this study, we performed digital image analysis with avian vision modelling to demonstrate that elytral luminance, which is attributed to proportions of elytral colour components, is fine-tuned to match local backgrounds. Field predation experiments with model beetles showed that better luminance matching resulted in a lower attack rate and corresponding lower mortality. Using restriction site-associated DNA (RAD) sequence data, we analysed the dispersal and evolution of colour pattern across geographical locations. We found that sand colour matching occurred irrespective of genetic and geographical distances between populations, suggesting that locally adapted colour patterns evolved after the colonization of these habitats. Given that beetle elytral colour patterns presumably have a quantitative genetic basis, our findings demonstrate that fine-tuning of background-matching camouflage to local habitat conditions can be attained through selection by visual predators, as predicted by the earliest proponent of natural selection.

Polygenic Adaptation in a Population of Finite Size

Polygenic adaptation in response to selection on quantitative traits has become an important topic in evolutionary biology. Here we review the recent literature on models of polygenic adaptation. In particular, we focus on a model that includes mutation and both directional and stabilizing selection on a highly polygenic trait in a population of finite size (thus experiencing random genetic drift). Assuming that a sudden environmental shift of the fitness optimum occurs while the population is in a stochastic equilibrium, we analyze the adaptation of the trait to the new optimum. When the shift is not too large relative to the equilibrium genetic variance and this variance is determined by loci with mostly small effects, the approach of the mean phenotype to the optimum can be approximated by a rapid exponential process (whose rate is proportional to the genetic variance). During this rapid phase the underlying changes to allele frequencies, however, may depend strongly on genetic drift. While trait-increasing alleles with intermediate equilibrium frequencies are dominated by selection and contribute positively to changes of the trait mean (i.e., are aligned with the direction of the optimum shift), alleles with low or high equilibrium frequencies show more of a random dynamics, which is expected when drift is dominating. A strong effect of drift is also predicted for population size bottlenecks. Our simulations show that the presence of a bottleneck results in a larger deviation of the population mean of the trait from the fitness optimum, which suggests that more loci experience the influence of drift.

Genome sequencing and population genomics modeling provide insights into the local adaptation of weeping forsythia

Understanding the genetic basis underlying the local adaptation of nonmodel species is a fundamental goal in evolutionary biology. In this study, we explored the genetic mechanisms of the local adaptation of Forsythia suspensa using genome sequence and population genomics data obtained from specific-locus amplified fragment sequencing. We assembled a high-quality reference genome of weeping forsythia (Scaffold N50 = 7.3 Mb) using ultralong Nanopore reads. Then, genome-wide comparative analysis was performed for 15 natural populations of weeping forsythia across its current distribution range. Our results revealed that candidate genes associated with local adaptation are functionally correlated with solar radiation, temperature and water variables across heterogeneous environmental scenarios. In particular, solar radiation during the period of fruit development and seed drying after ripening, cold, and drought significantly contributed to the adaptive differentiation of F. suspensa. Natural selection exerted by environmental factors contributed substantially to the population genetic structure of F. suspensa. Our results supported the hypothesis that adaptive differentiation should be highly pronounced in the genes involved in signal crosstalk between different environmental variables. Our population genomics study of F. suspensa provides insights into the fundamental genetic mechanisms of the local adaptation of plant species to climatic gradients.

Rapid local adaptation linked with phenotypic plasticity

Models of "plasticity-first" evolution are attractive because they explain the rapid evolution of new complex adaptations. Nevertheless, it is unclear whether plasticity can facilitate rapid microevolutionary change between diverging populations. Here, we show how plasticity may have generated adaptive differences in fecundity between neighboring wild populations of burying beetles Nicrophorus vespilloides. These populations occupy distinct Cambridgeshire woodlands that are just 2.5 km apart and that probably originated from a common ancestral population about 1000-4000 years ago. We find that populations are divergently adapted to breed on differently sized carrion. Adaptive differences in clutch size and egg size are associated with divergence at loci connected with oogenesis. The populations differ specifically in the elevation of the reaction norm linking clutch size to carrion size (i.e., genetic accommodation), and in the likelihood that surplus offspring will be lost after hatching. We suggest that these two processes may have facilitated rapid local adaptation on a fine-grained spatial scale.

Recent population genomic insights into the genetic basis of arsenic tolerance in humans: the difficulties of identifying positively selected loci in strongly bottlenecked populations

Recent advances in genomics have enabled researchers to shed light on the evolutionary processes driving human adaptation, by revealing the genetic architectures underlying traits ranging from lactase persistence, to skin pigmentation, to hypoxic response, to arsenic tolerance. Complicating the identification of targets of positive selection in modern human populations is their complex demographic history, characterized by population bottlenecks and expansions, population structure, migration, and admixture. In particular, founder effects and recent strong population size reductions, such as those experienced by the indigenous peoples of the Americas, have severe impacts on genetic variation that can lead to the accumulation of large allele frequency differences between populations due to genetic drift rather than natural selection. While distinguishing the effects of demographic history from selection remains challenging, neglecting neutral processes can lead to the incorrect identification of candidate loci. We here review the recent population genomic insights into the genetic basis of arsenic tolerance in Andean populations, and utilize this example to highlight both the difficulties pertaining to the identification of local adaptations in strongly bottlenecked populations, as well as the importance of controlling for demographic history in selection scans.

Selective sweeps on novel and introgressed variation shape mimicry loci in a butterfly adaptive radiation

Natural selection leaves distinct signatures in the genome that can reveal the targets and history of adaptive evolution. By analysing high-coverage genome sequence data from 4 major colour pattern loci sampled from nearly 600 individuals in 53 populations, we show pervasive selection on wing patterns in the Heliconius adaptive radiation. The strongest signatures correspond to loci with the greatest phenotypic effects, consistent with visual selection by predators, and are found in colour patterns with geographically restricted distributions. These recent sweeps are similar between co-mimics and indicate colour pattern turn-over events despite strong stabilising selection. Using simulations, we compare sweep signatures expected under classic hard sweeps with those resulting from adaptive introgression, an important aspect of mimicry evolution in Heliconius butterflies. Simulated recipient populations show a distinct 'volcano' pattern with peaks of increased genetic diversity around the selected target, characteristic of sweeps of introgressed variation and consistent with diversity patterns found in some populations. Our genomic data reveal a surprisingly dynamic history of colour pattern selection and co-evolution in this adaptive radiation.

Do different rates of gene flow underlie variation in phenotypic and phenological clines in a montane grasshopper community?

Species responses to environmental change are likely to depend on existing genetic and phenotypic variation, as well as evolutionary potential. A key challenge is to determine whether gene flow might facilitate or impede genomic divergence among populations responding to environmental change, and if emergent phenotypic variation is dependent on gene flow rates. A general expectation is that patterns of genetic differentiation in a set of codistributed species reflect differences in dispersal ability. In less dispersive species, we predict greater genetic divergence and reduced gene flow. This could lead to covariation in life-history traits due to local adaptation, although plasticity or drift could mirror these patterns. We compare genome-wide patterns of genetic structure in four phenotypically variable grasshopper species along a steep elevation gradient near Boulder, Colorado, and test the hypothesis that genomic differentiation is greater in short-winged grasshopper species, and statistically associated with variation in growth, reproductive, and physiological traits along this gradient. In addition, we estimate rates of gene flow under competing demographic models, as well as potential gene flow through surveys of phenological overlap among populations within a species. All species exhibit genetic structure along the elevation gradient and limited gene flow. The most pronounced genetic divergence appears in short-winged (less dispersive) species, which also exhibit less phenological overlap among populations. A high-elevation population of the most widespread species, Melanoplus sanguinipes, appears to be a sink population derived from low elevation populations. While dispersal ability has a clear connection to the genetic structure in different species, genetic distance does not predict growth, reproductive, or physiological trait variation in any species, requiring further investigation to clearly link phenotypic divergence to local adaptation.

Parasite-mediated selection in a natural metapopulation of Daphnia magna

Parasite-mediated selection varying across time and space in metapopulations is expected to result in host local adaptation and the maintenance of genetic diversity in disease-related traits. However, nonadaptive processes like migration and extinction-(re)colonization dynamics might interfere with adaptive evolution. Understanding how adaptive and nonadaptive processes interact to shape genetic variability in life-history and disease-related traits can provide important insights into their evolution in subdivided populations. Here we investigate signatures of spatially fluctuating, parasite-mediated selection in a natural metapopulation of Daphnia magna. Host genotypes from infected and uninfected populations were genotyped at microsatellite markers, and phenotyped for life-history and disease traits in common garden experiments. Combining phenotypic and genotypic data a Q_ST -F_ST -like analysis was conducted to test for signatures of parasite mediated selection. We observed high variation within and among populations for phenotypic traits, but neither an indication of host local adaptation nor a cost of resistance. Infected populations have a higher gene diversity (Hs) than uninfected populations and Hs is strongly positively correlated with fitness. These results suggest a strong parasite effect on reducing population level inbreeding. We discuss how stochastic processes related to frequent extinction-(re)colonization dynamics as well as host and parasite migration impede the evolution of resistance in the infected populations. We suggest that the genetic and phenotypic patterns of variation are a product of dynamic changes in the host gene pool caused by the interaction of colonization bottlenecks, inbreeding, immigration, hybrid vigor, rare host genotype advantage and parasitism. Our study highlights the effect of the parasite in ameliorating the negative fitness consequences caused by the high drift load in this metapopulation.

The population genetics of crypsis in vertebrates: recent insights from mice, hares, and lizards

By combining well-established population genetic theory with high-throughput sequencing data from natural populations, major strides have recently been made in understanding how, why, and when vertebrate populations evolve crypsis. Here, we focus on background matching, a particular facet of crypsis that involves the ability of an organism to conceal itself through matching its color to the surrounding environment. While interesting in and of itself, the study of this phenotype has also provided fruitful population genetic insights into the interplay of strong positive selection with other evolutionary processes. Specifically, and predicated upon the findings of previous candidate gene association studies, a primary focus of this recent literature involves the realization that the inference of selection from DNA sequence data first requires a robust model of population demography in order to identify genomic regions which do not conform to neutral expectations. Moreover, these demographic estimates provide crucial information about the origin and timing of the onset of selective pressures associated with, for example, the colonization of a novel environment. Furthermore, such inference has revealed crypsis to be a particularly useful phenotype for investigating the interplay of migration and selection—with examples of gene flow constraining rates of adaptation, or alternatively providing the genetic variants that may ultimately sweep through the population. Here, we evaluate the underlying evidence, review the strengths and weaknesses of the many population genetic methodologies used in these studies, and discuss how these insights have aided our general understanding of the evolutionary process.

Genomic Patterns of Local Adaptation under Gene Flow in Arabidopsis lyrata

Short-scale local adaptation is a complex process involving selection, migration, and drift. The expected effects on the genome are well grounded in theory but examining these on an empirical level has proven difficult, as it requires information about local selection, demographic history, and recombination rate variation. Here, we use locally adapted and phenotypically differentiated Arabidopsis lyrata populations from two altitudinal gradients in Norway to test these expectations at the whole-genome level. Demography modeling indicates that populations within the gradients diverged <2 kya and that the sites are connected by gene flow. The gene flow estimates are, however, highly asymmetric with migration from high to low altitudes being several times more frequent than vice versa. To detect signatures of selection for local adaptation, we estimate patterns of lineage-specific differentiation among these populations. Theory predicts that gene flow leads to concentration of adaptive loci in areas of low recombination; a pattern we observe in both lowland-alpine comparisons. Although most selected loci display patterns of conditional neutrality, we found indications of genetic trade-offs, with one locus particularly showing high differentiation and signs of selection in both populations. Our results further suggest that resistance to solar radiation is an important adaptation to alpine environments, while vegetative growth and bacterial defense are indicated as selected traits in the lowland habitats. These results provide insights into genetic architectures and evolutionary processes driving local adaptation under gene flow. We also contribute to understanding of traits and biological processes underlying alpine adaptation in northern latitudes.

Tipping the Scales: The Migration-Selection Balance Leans toward Selection in Snake Venoms

The migration-selection interaction is the strongest determinant of whether a beneficial allele increases in frequency within a population. Results of empirical studies examining the role of gene flow in an adaptive context, however, have largely been equivocal, with examples of opposing outcomes being repeatedly documented (e.g., local adaptation with high levels of gene flow vs. gene swamping). We compared neutral genomic and venom expression divergence for three sympatric pit vipers with differing ecologies to determine if and how migration-selection disequilibria result in local adaptation. We specifically tested whether neutral differentiation predicted phenotypic differentiation within an isolation-by-distance framework. The decoupling of neutral and phenotypic differentiation would indicate selection led to adaptive divergence irrespective of migration, whereas a significant relationship between neutral and venom expression differentiation would provide evidence in favor of the constraining force of gene flow. Neutral differentiation and geographic distance predicted phenotypic differentiation only in the generalist species, indicating that selection was the predominant mechanism in the migration-selection balance underlying adaptive venom evolution in both specialists. Dispersal is thought to be a stronger influence on genetic differentiation than specialization, but our results suggest the opposite. Greater specialization may lead to greater diversification rates, and extreme spatial and temporal variation in selective pressures can favor generalist phenotypes evolving under strong stabilizing selection. Our results are consistent with these expectations, suggesting that the equivocal findings of studies examining the role of gene flow in an adaptive context may be explained by ecological specialization theory.

Sampling related individuals within ponds biases estimates of population structure in a pond-breeding amphibian

Effective conservation and management of pond-breeding amphibians depends on the accurate estimation of population structure, demographic parameters, and the influence of landscape features on breeding-site connectivity. Population-level studies of pond-breeding amphibians typically sample larval life stages because they are easily captured and can be sampled nondestructively. These studies often identify high levels of relatedness between individuals from the same pond, which can be exacerbated by sampling the larval stage. Yet, the effect of these related individuals on population genetic studies using genomic data is not yet fully understood. Here, we assess the effect of within-pond relatedness on population and landscape genetic analyses by focusing on the barred tiger salamanders (Ambystoma mavortium) from the Nebraska Sandhills. Utilizing genome-wide SNPs generated using a double-digest RADseq approach, we conducted standard population and landscape genetic analyses using datasets with and without siblings. We found that reduced sample sizes influenced parameter estimates more than the inclusion of siblings, but that within-pond relatedness led to the inference of spurious population structure when analyses depended on allele frequencies. Our landscape genetic analyses also supported different models across datasets depending on the spatial resolution analyzed. We recommend that future studies not only test for relatedness among larval samples but also remove siblings before conducting population or landscape genetic analyses. We also recommend alternative sampling strategies to reduce sampling siblings before sequencing takes place. Biases introduced by unknowingly including siblings can have significant implications for population and landscape genetic analyses, and in turn, for species conservation strategies and outcomes.