Patterns and effects of gene flow on adaptation across spatial scales: implications for management

Gene flow can have rapid effects on adaptation and is an important evolutionary tool available when undertaking biological conservation and restoration. This tool is underused partly because of the perceived risk of outbreeding depression and loss of mean fitness when different populations are crossed. In this article, we briefly review some theory and empirical findings on how genetic variation is distributed across species ranges, describe known patterns of gene flow in nature with respect to environmental gradients, and highlight the effects of gene flow on adaptation in small or stressed populations in challenging environments (e.g., at species range limits). We then present a case study involving crosses at varying spatial scales among mountain populations of a trigger plant (Stylidium armeria: Stylidiaceae) in the Australian Alps to highlight how some issues around gene flow effects can be evaluated. We found evidence of outbreeding depression in seed production at greater geographic distances. Nevertheless, we found no evidence of maladaptive gene flow effects in likelihood of germination, plant performance (size), and performance variance, suggesting that gene flow at all spatial scales produces offspring with high adaptive potential. This case study demonstrates a path to evaluating how increasing sources of gene flow in managed wild and restored populations could identify some offspring with high fitness that could bolster the ability of populations to adapt to future environmental changes. We suggest further ways in which managers and researchers can act to understand and consider adaptive gene flow in natural and conservation contexts under rapidly changing conditions.

Recruitment of European sea bass (Dicentrarchus labrax) in northerly UK estuaries indicates a mismatch between spawning and fisheries closure periods

European sea bass (Dicentrarchus labrax) is a species of high commercial and recreational value, but it exhibits highly variable recruitment rates and has been subject to recent declines. Emergency management measures put in place to protect spawning stocks include the annual closure of commercial and recreational fisheries over a 2-month, February-March, window. Whether this protection measure is having the desired outcome for this data-poor species remains unclear. Otolith microstructural analyses (counts and widths of daily growth rings and check marks indicative of settlement) were used to estimate (1) spawn timing, (2) pelagic larval duration and settlement timing, (3) growth rate and condition, and (4) the otolith-fish size relationship for juvenile European sea bass caught from two estuaries in Wales (Dwyryd, Y Foryd), located at the northern edge of the species range. We observed a significant mismatch between the timing of fisheries closures and the spawning, with 99.2% of recruits having been spawned after the fishery had reopened (back-calculated median spawn date = May 5 ± 17 days SD), suggesting that the closure may be too early to adequately protect this population. Further, we present the first empirically derived estimates of pelagic larval duration for sea bass recruits settling in UK habitats, which showed a strong negative relationship with spawn date. Finally, we found significant differences in fish condition between the two estuaries, suggesting local variation in habitat quality. The results suggest that the timing of current fisheries closures may not be adequately protecting the spawners supplying these northernmost estuaries, which are likely to become increasingly important as sea bass distributions shift northward in our climate future.

Plant kleptomaniacs: geographical genetic patterns in the amphi-apomictic Rubus ser. Glandulosi (Rosaceae) reveal complex reticulate evolution of Eurasian brambles

BACKGROUND AND AIMS

Rubus ser. Glandulosi provides a unique model of geographical parthenogenesis on a homoploid (2n = 4x) level. We aim to characterize evolutionary and phylogeographical patterns in this taxon and shed light on the geographical differentiation of apomicts and sexuals. Ultimately, we aim to evaluate the importance of phylogeography in the formation of geographical parthenogenesis.

METHODS

Rubus ser. Glandulosi was sampled across its Eurasian range together with other co-occurring Rubus taxa (587 individuals in total). Double-digest restriction site-associated DNA sequencing (ddRADseq) and modelling of suitable climate were used for evolutionary inferences.

KEY RESULTS

Six ancestral species were identified that contributed to the contemporary gene pool of R. ser. Glandulosi. Sexuals were introgressed from Rubus dolichocarpus and Rubus moschus in West Asia and from Rubus ulmifolius agg., Rubus canescens and Rubus incanescens in Europe, whereas apomicts were characterized by alleles of Rubus subsect. Rubus. Gene flow between sexuals and apomicts was also detected, as was occasional hybridization with other taxa.

CONCLUSIONS

We hypothesize that sexuals survived the last glacial period in several large southern refugia, whereas apomicts were mostly restricted to southern France, whence they quickly recolonized Central and Western Europe. The secondary contact of sexuals and apomicts was probably the principal factor that established geographical parthenogenesis in R. ser. Glandulosi. Sexual populations are not impoverished in genetic diversity along their borderline with apomicts, and maladaptive population genetic processes probably did not shape the geographical patterns.

Niche-based approach to explore the impacts of environmental disturbances on biodiversity

Globally, species are increasingly at risk from compounding threatening processes, an increasingly prominent driver of which is environmental disturbances. To facilitate effective conservation efforts following such events, methods that evaluate potential impacts across multiple species and provide landscape-scale information are needed to guide targeted responses. Often, the geographic overlap between a disturbance and species' distribution is calculated and then used as a proxy for potential impact. However, such methods do not account for the important influence of environmental heterogeneity throughout species' ranges. To address this shortcoming, we quantified the effects of environmental disturbances on species' environmental niche space. Using the Australian 2019 and 2020 Black Summer fires as a case study, we applied a niche-centric approach to examine the potential impacts of these fires on 387 vertebrate species. We examined the utility of established and novel niche metrics to assess the potential impacts of large-scale disturbance events on species by comparing the potential effects of the fires as determined by our various niche measures to those derived from geographic-based measures of impact. We examined the quality of environmental space affected by the disturbance by quantifying the position in niche space where the disturbance occurred (center or margin), the uniqueness of the environmental space that was burned, and the degree to which the remaining, unburned portion of the niche differed from a species' original prefire niche. There was limited congruence between the proportion of geographic and niche space affected, which showed that geographic-based approaches in isolation may have underestimated the impact of the fires for 56% of modeled species. For each species, when combined, these metrics provided a greater indication of postdisturbance recovery potential than geographic-based measures alone. Accordingly, the integration of niche-based analyses into conservation assessments following large-scale disturbance events will lead to a more nuanced understanding of potential impacts and guide more informed and effective conservation actions.

Native prey, not landscape change or novel prey, drive cougar (Puma concolor) distribution at a boreal forest range edge

Many large carnivores, despite widespread habitat alteration, are rebounding in parts of their former ranges after decades of persecution and exploitation. Cougars (Puma concolor) are apex predator with their remaining northern core range constricted to mountain landscapes and areas of western North America; however, cougar populations have recently started rebounding in several locations across North America, including northward in boreal forest landscapes. A camera-trap survey of multiple landscapes across Alberta, Canada, delineated a range edge; within this region, we deployed an array of 47 camera traps in a random stratified design across a landscape spanning a gradient of anthropogenic development relative to the predicted expansion front. We completed multiple hypotheses in an information-theoretic framework to determine if cougar occurrence is best explained by natural land cover features, anthropogenic development features, or competitor and prey activity. We predicted that anthropogenic development features from resource extraction and invading white-tailed deer (Odocoileus virgianius) explain cougar distribution at this boreal range edge. Counter to our predictions, the relative activity of native prey, predominantly snowshoe hare (Lepus americanus), was the best predictor of cougar occurrence at this range edge. Small-bodied prey items are particularly important for female and sub-adult cougars and may support breeding individuals in the northeast boreal forest. Also, counter to our predictions, there was not a strong relationship detected between cougar occurrence and gray wolf (Canis lupus) activity at this range edge. However, further investigation is recommended as the possibility of cougar expansion into areas of the multi-prey boreal system, where wolves have recently been controlled, could have negative consequences for conservation goals in this region (e.g. the recovery of woodland caribou [Rangifer tarandus caribou]). Our study highlights the need to monitor contemporary distributions to inform conservation management objectives as large carnivores recover across North America.

Recent range shifts of moths, butterflies, and birds are driven by the breadth of their climatic niche

Species are altering their ranges as a response to climate change, but the magnitude and direction of observed range shifts vary considerably among species. The ability to persist in current areas and colonize new areas plays a crucial role in determining which species will thrive and which decline as climate change progresses. Several studies have sought to identify characteristics, such as morphological and life-history traits, that could explain differences in the capability of species to shift their ranges together with a changing climate. These characteristics have explained variation in range shifts only sporadically, thus offering an uncertain tool for discerning responses among species. As long-term selection to past climates have shaped species' tolerances, metrics describing species' contemporary climatic niches may provide an alternative means for understanding responses to on-going climate change. Species that occur in a broader range of climatic conditions may hold greater tolerance to climatic variability and could therefore more readily maintain their historical ranges, while species with more narrow tolerances may only persist if they are able to shift in space to track their climatic niche. Here, we provide a first-filter test of the effect of climatic niche dimensions on shifts in the leading range edges in three relatively well-dispersing species groups. Based on the realized changes in the northern range edges of 383 moth, butterfly, and bird species across a boreal 1,100 km latitudinal gradient over c. 20 years, we show that while most morphological or life-history traits were not strongly connected with range shifts, moths and birds occupying a narrower thermal niche and butterflies occupying a broader moisture niche across their European distribution show stronger shifts towards the north. Our results indicate that the climatic niche may be important for predicting responses under climate change and as such warrants further investigation of potential mechanistic underpinnings.

Heat stress reveals a fertility debt owing to postcopulatory sexual selection

Climates are changing rapidly, demanding equally rapid adaptation of natural populations. Whether sexual selection can aid such adaptation is under debate; while sexual selection should promote adaptation when individuals with high mating success are also best adapted to their local surroundings, the expression of sexually selected traits can incur costs. Here we asked what the demographic consequences of such costs may be once climates change to become harsher and the strength of natural selection increases. We first adopted a classic life history theory framework, incorporating a trade-off between reproduction and maintenance, and applied it to the male germline to generate formalized predictions for how an evolutionary history of strong postcopulatory sexual selection (sperm competition) may affect male fertility under acute adult heat stress. We then tested these predictions by assessing the thermal sensitivity of fertility (TSF) in replicated lineages of seed beetles maintained for 68 generations under three alternative mating regimes manipulating the opportunity for sexual and natural selection. In line with the theoretical predictions, we find that males evolving under strong sexual selection suffer from increased TSF. Interestingly, females from the regime under strong sexual selection, who experienced relaxed selection on their own reproductive effort, had high fertility in benign settings but suffered increased TSF, like their brothers. This implies that female fertility and TSF evolved through genetic correlation with reproductive traits sexually selected in males. Paternal but not maternal heat stress reduced offspring fertility with no evidence for adaptive transgenerational plasticity among heat-exposed offspring, indicating that the observed effects may compound over generations. Our results suggest that trade-offs between fertility and traits increasing success in postcopulatory sexual selection can be revealed in harsh environments. This can put polyandrous species under immediate risk during extreme heat waves expected under future climate change.

Divergent mechanisms of reduced growth performance in Betula ermanii saplings from high-altitude and low-latitude range edges

The reduced growth performance of individuals from range edges is a common phenomenon in various taxa, and considered to be an evolutionary factor that limits the species' range. However, most studies did not distinguish between two mechanisms that can lead to this reduction: genetic load and adaptive selection to harsh conditions. To address this lack of understanding, we investigated the climatic and genetic factors underlying the growth performance of Betula ermanii saplings transplanted from 11 populations including high-altitude edge and low-latitude edge population. We estimated the climatic position of the populations within the overall B. ermanii's distribution, and the genetic composition and diversity using restriction-site associated DNA sequencing, and measured survival, growth rates and individual size of the saplings. The high-altitude edge population (APW) was located below the 95% significance interval for the mean annual temperature range, but did not show any distinctive genetic characteristics. In contrast, the low-latitude edge population (SHK) exhibited a high level of linkage disequilibrium, low genetic diversity, a distinct genetic composition from the other populations, and a high relatedness coefficient. Both APW and SHK saplings displayed lower survival rates, heights and diameters, while SHK saplings also exhibited lower growth rates than the other populations' saplings. The low heights and diameters of APW saplings was likely the result of adaptive selection to harsh conditions, while the low survival and growth rates of SHK saplings was likely the result of genetic load. Our findings shed light on the mechanisms underlying the reduced growth performance of range-edge populations.

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.

Bacterial motility can govern the dynamics of antibiotic resistance evolution

Spatial heterogeneity in antibiotic concentrations is thought to accelerate the evolution of antibiotic resistance, but current theory and experiments have overlooked the effect of cell motility on bacterial adaptation. Here, we study bacterial evolution in antibiotic landscapes with a quantitative model where bacteria evolve under the stochastic processes of proliferation, death, mutation and migration. Numerical and analytical results show that cell motility can both accelerate and decelerate bacterial adaptation by affecting the degree of genotypic mixing and ecological competition. Moreover, we find that for sufficiently high rates, cell motility can limit bacterial survival, and we derive conditions for all these regimes. Similar patterns are observed in more complex scenarios, namely where bacteria can bias their motion in chemical gradients (chemotaxis) or switch between motility phenotypes either stochastically or in a density-dependent manner. Overall, our work reveals limits to bacterial adaptation in antibiotic landscapes that are set by cell motility.

Biogeography and the evolution of acoustic communication in the polyploid North American grey treefrog complex

After polyploid species are formed, interactions between diploid and polyploid lineages may generate additional diversity in novel cytotypes and phenotypes. In anurans, mate choice by acoustic communication is the primary method by which individuals identify their own species and assess suitable mates. As such, the evolution of acoustic signals is an important mechanism for contributing to reproductive isolation and diversification in this group. Here, we estimate the biogeographical history of the North American grey treefrog complex, consisting of the diploid Hyla chrysoscelis and the tetraploid Hyla versicolor, focusing specifically on the geographical origin of whole genome duplication and the expansion of lineages out of glacial refugia. We then test for lineage-specific differences in mating signals by applying comparative methods to a large acoustic data set collected over 52 years that includes >1500 individual frogs. Along with describing the overall biogeographical history and call diversity, we found evidence that the geographical origin of H. versicolor and the formation of the midwestern polyploid lineage are both associated with glacial limits, and that the southwestern polyploid lineage is associated with a shift in acoustic phenotype relative to the diploid lineage with which they share a mitochondrial lineage. In H. chrysoscelis, we see that acoustic signals are largely split by Eastern and Western lineages, but that northward expansion along either side of the Appalachian Mountains is associated with further acoustic diversification. Overall, results of this study provide substantial clarity on the evolution of grey treefrogs as it relates to their biogeography and acoustic communication.

The evolution of plasticity at geographic range edges

Phenotypic plasticity enables rapid responses to environmental change, and could facilitate range shifts in response to climate change. What drives the evolution of plasticity at range edges, and the capacity of range-edge individuals to be plastic, remain unclear. Here, we propose that accurately predicting when plasticity itself evolves or mediates adaptive evolution at expanding range edges requires integrating knowledge on the demography and evolution of edge populations. Our synthesis shows that: (i) the demography of edge populations can amplify or attenuate responses to selection for plasticity through diverse pathways, and (ii) demographic effects on plasticity are modified by the stability of range edges. Our spatially explicit synthesis for plasticity has the potential to improve predictions for range shifts with climate change.

Gene flow drives genomic diversity in Asian Pikas distributed along the core and range-edge habitats in the Himalayas

Studying the genetic variation among different species distributed across their core and range-edge habitats can provide valuable insights into how genetic variation changes across the species' distribution range. This information can be important for understanding local adaptation, as well as for conservation and management efforts. In this study, we have carried out genomic characterization of six species of Asian Pikas distributed along their core and range-edge habitats in the Himalayas. We utilized a population genomics approach using ~28,000 genome-wide SNP markers obtained from restriction-site associated DNA sequencing. We identified low nucleotide diversity and high inbreeding coefficients in all six species across their core and range-edge habitats. We also identified evidence of gene flow among genetically diverse species. Our results provide evidence of reduced genetic diversity in Asian pikas distributed across the Himalayas and the neighboring regions and indicate that recurrent gene flow is possibly a key mechanism for maintaining genetic diversity and adaptive potential in these pikas. However, full-scale genomics studies that utilize whole-genome sequencing approaches will be needed to quantify the direction and timing of gene flow and functional changes associated with introgressed regions in the genome. Our results represent an important step toward understanding the patterns and consequences of gene flow in species, sampled at the least studied, yet climatically vulnerable part of their habitat that can be further used to inform conservation strategies that promote connectivity and gene flow between populations.

Acute cold stress and supercooling capacity of Mediterranean fruit fly populations across the Northern Hemisphere (Middle East and Europe)

The Mediterranean fruit fly, Ceratitis capitata (Diptera: Tephritidae), holds an impressive record of successful invasion events promoted by globalization in fruit trade and human mobility. In addition, C. capitata is gradually expanding its geographic distribution to cooler temperate areas of the Northern Hemisphere. Cold tolerance of C. capitata seems to be a crucial feature that promotes population establishment and hence invasion success. To elucidate the interplay between the invasion process in the northern hemisphere and cold tolerance of geographically isolated populations of C. capitata, we determined (a) the response to acute cold stress survival of adults, and (b) the supercooling capacity (SCP) of immature stages and adults. To assess the phenotypic plasticity in these populations, the effect of acclimation to low temperatures on acute cold stress survival in adults was also examined. The results revealed that survival after acute cold stress was positively related to low temperature acclimation, except for females originating from Thessaloniki (northern Greece). Adults from the warmer environment of South Arava (Israel) were less tolerant after acute cold stress compared with those from Heraklion (Crete, Greece) and Thessaloniki. Plastic responses to cold acclimation were population specific, with the South Arava population being more plastic compared to the two Greek populations. For SCP, the results revealed that there is little to no correlation between SCP and climate variables of the areas where C. capitata populations originated. SCP was much lower than the lowest temperature individuals are likely to experience in their respective habitats. These results set the stage for asking questions regarding the evolutionary adaptive processes that facilitate range expansions of C. capitata into cooler temperate areas of Europe.

Genome-environment associations along elevation gradients in two snowbed species of the North-Eastern Calcareous Alps

BACKGROUND

Anthropogenic climate change leads to increasing temperatures and altered precipitation and snowmelt patterns, especially in alpine ecosystems. To understand species' responses to climate change, assessment of genetic structure and diversity is crucial as the basis for the evaluation of migration patterns, genetic adaptation potential as well as the identification of adaptive alleles.

RESULTS

We studied genetic structure, diversity and genome-environment associations of two snowbed species endemic to the Eastern Alps with a large elevational range, Achillea clusiana Tausch and Campanula pulla L. Genotyping-by-sequencing was employed to assemble loci de novo, call variants and perform population genetic analyses. Populations of either species were distinguishable by mountain, and to some extent by elevation. We found evidence for gene flow between elevations. Results of genome-environment associations suggested similar selective pressures acting on both species, emanating mainly from precipitation and exposition rather than temperature.

CONCLUSIONS

Given their genetic structure and amount of gene flow among populations the two study species are suitable to serve as a model for genetic monitoring of climate change adaptation along an elevation gradient. Consequences of climate change will predominantly manifest via changes in precipitation and, thus, duration of snow cover in the snowbeds and indirectly via shrub encroachment accompanied by increasing shading of snowbeds at lower range margins. Assembling genomes of the study species and studying larger sample sizes and time series will be necessary to functionally characterize and validate the herein identified genomic loci putatively involved in adaptive processes.

Urban insect bioarks of the 21st century

Insects exhibit divergent biodiversity responses to cities. Many urban populations are not at equilibrium: biodiversity decline or recovery from environmental perturbation is often still in progress. Substantial variation in urban biodiversity patterns suggests the need to understand its mechanistic basis. In addition, current urban infrastructure decisions might profoundly influence future biodiversity trends. Although many nature-based solutions to urban climate problems also support urban insect biodiversity, trade-offs are possible and should be avoided to maximize biodiversity-climate cobenefits. Because insects are coping with the dual threats of urbanization and climate change, there is an urgent need to design cities that facilitate persistence within the city footprint or facilitate compensatory responses to global climate change as species transit through the city footprint.

Drivers of strong isolation and small effective population size at a leading range edge of a widespread plant

Climate change has influenced species distributions worldwide with upward elevational shifts observed in many systems. Leading range edge populations, like those at upper elevation limits, are crucial for climate change responses but can exhibit low genetic diversity due to founder effects, isolation, or limited outbreeding. These factors can hamper local adaptation at range limits. Using the widespread herb, Argentina anserina, we measured ecological attributes (population density on the landscape, area of population occupancy, and plant and flower density) spanning a 1000 m elevation gradient, with high elevation populations at the range limit. We measured vegetative clonal potential in the greenhouse for populations spanning the gradient. We combined these data with a ddRAD-seq dataset to test the hypotheses that high elevation populations would exhibit ecological and genomic signatures of leading range edge populations. We found that population density on the landscape declined towards the high elevation limit, as is expected towards range edges. However, plant density was elevated within edge populations. In the greenhouse, high elevation plants exhibited stronger clonal potential than low elevation plants, likely explaining increased plant density in the field. Phylogeographic analysis supported more recent colonization of high elevation populations which were also more genetically isolated, had more extreme heterozygote excess and had smaller effective population size than low. Results support that colonization of high elevations was likely accompanied by increased asexuality, contributing to a decline in effective population size. Despite high plant density in leading edge populations, their small effective size, isolation and clonality could constrain adaptive potential.

Understanding local plant extinctions before it is too late: bridging evolutionary genomics with global ecology

Understanding evolutionary genomic and population processes within a species range is key to anticipating the extinction of plant species before it is too late. However, most models of biodiversity risk under global change do not account for the genetic variation and local adaptation of different populations. Population diversity is critical to understanding extinction because different populations may be more or less susceptible to global change and, if lost, would reduce the total diversity within a species. Two new modeling frameworks advance our understanding of extinction from a population and evolutionary angle: Rapid climate change-driven disruptions in population adaptation are predicted from associations between genomes and local climates. Furthermore, losses of population diversity from global land-use transformations are estimated by scaling relationships of species' genomic diversity with habitat area. Overall, these global eco-evolutionary methods advance the predictability - and possibly the preventability - of the ongoing extinction of plant species.

Weak genetic differentiation but strong climate-induced selective pressure toward the rear edge of mountain pine in north-eastern Spain

Local differentiation at distribution limits may influence species' adaptive capacity to environmental changes. However, drivers, such gene flow and local selection, are still poorly understood. We focus on the role played by range limits in mountain forests to test the hypothesis that relict tree populations are subjected to genetic differentiation and local adaptation. Two alpine treelines of mountain pine (Pinus uncinata Ram. ex DC) were investigated in the Spanish Pyrenees. Further, an isolated relict population forming the species' southernmost distribution limit in north-eastern Spain was also investigated. Using genotyping by sequencing, a genetic matrix conformed by single nucleotide polymorphisms (SNPs) was obtained. This matrix was used to perform genotype-environment and genotype-phenotype associations, as well as to model risk of non-adaptedness. Increasing climate seasonality appears as an essential element in the interpretation of SNPs subjected to selective pressures. Genetic differentiations were overall weak. The differences in leaf mass area and radial growth rate, as well as the identification of several SNPs subjected to selective pressures, exceeded neutral predictions of differentiation among populations. Despite genetic drift might prevail in the isolated population, the Fst values (0.060 and 0.066) showed a moderate genetic drift and Nm values (3.939 and 3.555) indicate the presence of gene flow between the relict population and both treelines. Nonetheless, the SNPs subjected to selection pressures provide evidences of possible selection in treeline ecotones. Persistence in range boundaries seems to involve several selective pressures in species' traits, which were significantly related to enhanced drought seasonality at the limit of P. uncinata distribution range. We conclude that gene flow is unlikely to constrain adaptation in the P. uncinata rear edge, although this species shows vulnerability to future climate change scenarios involving warmer and drier conditions.

Genotype-environment associations to reveal the molecular basis of environmental adaptation

A fundamental goal in plant biology is to identify and understand the variation underlying plants' adaptation to their environment. Climate change has given new urgency to this goal, as society aims to accelerate adaptation of ecologically important plant species, endangered plant species, and crops to hotter, less predictable climates. In the pre-genomic era, identifying adaptive alleles was painstaking work, leveraging genetics, molecular biology, physiology, and ecology. Now, the rise of genomics and new computational approaches may facilitate this research. Genotype-environment associations (GEAs) use statistical associations between allele frequency and environment of origin to test the hypothesis that allelic variation at a given gene is adapted to local environments. Researchers may scan the genome for GEAs to generate hypotheses on adaptive genetic variants (environmental genome-wide association studies). Despite the rapid adoption of these methods, many important questions remain about the interpretation of GEA findings, which arise from fundamental unanswered questions on the genetic architecture of adaptation and limitations inherent to association-based analyses. We outline strategies to ground GEAs in the underlying hypotheses of genetic architecture and better test GEA-generated hypotheses using genetics and ecophysiology. We provide recommendations for new users who seek to learn about the molecular basis of adaptation. When combined with a rigorous hypothesis testing framework, GEAs may facilitate our understanding of the molecular basis of climate adaptation for plant improvement.

Using temporal genomics to understand contemporary climate change responses in wildlife

Monitoring the evolutionary responses of species to ongoing global climate change is critical for informing conservation. Population genomic studies that use samples from multiple time points (“temporal genomics”) are uniquely able to make direct observations of change over time. Consequently, only temporal studies can show genetic erosion or spatiotemporal changes in population structure. Temporal genomic studies directly examining climate change effects are currently rare but will likely increase in the coming years due to their high conservation value. Here, we highlight four key genetic indicators that can be monitored using temporal genomics to understand how species are responding to climate change. All indicators crucially rely on having a suitable baseline that accurately represents the past condition of the population, and we discuss aspects of study design that must be considered to achieve this.

The relative role of plasticity and demographic history in Capsella bursa-pastoris: a common garden experiment in Asia and Europe

The colonization success of a species depends on the interplay between its phenotypic plasticity, adaptive potential and demographic history. Assessing their relative contributions during the different phases of a species range expansion is challenging, and requires large-scale experiments. Here, we investigated the relative contributions of plasticity, performance and demographic history to the worldwide expansion of the shepherd's purse, Capsella bursa-pastoris. We installed two large common gardens of the shepherd's purse, a young, self-fertilizing, allopolyploid weed with a worldwide distribution. One common garden was located in Europe, the other in Asia. We used accessions from three distinct genetic clusters (Middle East, Europe and Asia) that reflect the demographic history of the species. Several life-history traits were measured. To explain the phenotypic variation between and within genetic clusters, we analysed the effects of (i) the genetic clusters, (ii) the phenotypic plasticity and its association to fitness and (iii) the distance in terms of bioclimatic variables between the sampling site of an accession and the common garden, i.e. the environmental distance. Our experiment showed that (i) the performance of C. bursa-pastoris is closely related to its high phenotypic plasticity; (ii) within a common garden, genetic cluster was a main determinant of phenotypic differences; and (iii) at the scale of the experiment, the effect of environmental distance to the common garden could not be distinguished from that of genetic clusters. Phenotypic plasticity and demographic history both play important role at different stages of range expansion. The success of the worldwide expansion of C. bursa-pastoris was undoubtedly influenced by its strong phenotypic plasticity.

The role of demographic compensation in stabilising marginal tree populations in North America

Demographic compensation-the opposing responses of vital rates along environmental gradients-potentially delays anticipated species' range contraction under climate change, but no consensus exists on its actual contribution. We calculated population growth rate (λ) and demographic compensation across the distributional ranges of 81 North American tree species and examined their responses to simulated warming and tree competition. We found that 43% of species showed stable population size at both northern and southern edges. Demographic compensation was detected in 25 species, yet 15 of them still showed a potential retraction from southern edges, indicating that compensation alone cannot maintain range stability. Simulated climatic warming caused larger decreases in λ for most species and weakened the effectiveness of demographic compensation in stabilising ranges. These findings suggest that climate stress may surpass the limited capacity of demographic compensation and pose a threat to the viability of North American tree populations.

Increasing temporal variance leads to stable species range limits

What prevents populations of a species from adapting to the novel environments outside the species' geographic distribution? Previous models highlighted how gene flow across spatial environmental gradients determines species expansion versus extinction and the location of species range limits. However, space is only one of two axes of environmental variation-environments also vary in time, and we know temporal environmental variation has important consequences for population demography and evolution. We used analytical and individual-based evolutionary models to explore how temporal variation in environmental conditions influences the spread of populations across a spatial environmental gradient. We find that temporal variation greatly alters our predictions for range dynamics compared to temporally static environments. When temporal variance is equal across the landscape, the fate of species (expansion versus extinction) is determined by the interaction between the degree of temporal autocorrelation in environmental fluctuations and the steepness of the spatial environmental gradient. When the magnitude of temporal variance changes across the landscape, stable range limits form where this variance increases maladaptation sufficiently to prevent local persistence. These results illustrate the pivotal influence of temporal variation on the likelihood of populations colonizing novel habitats and the location of species range limits.

Population expansion, divergence, and persistence in Western Fence Lizards (Sceloporus occidentalis) at the northern extreme of their distributional range

Population dynamics within species at the edge of their distributional range, including the formation of genetic structure during range expansion, are difficult to study when they have had limited time to evolve. Western Fence Lizards (Sceloporus occidentalis) have a patchy distribution at the northern edge of their range around the Puget Sound, Washington, where they almost exclusively occur on imperiled coastal habitats. The entire region was covered by Pleistocene glaciation as recently as 16,000 years ago, suggesting that populations must have colonized these habitats relatively recently. We tested for population differentiation across this landscape using genome-wide SNPs and morphological data. A time-calibrated species tree supports the hypothesis of a post-glacial establishment and subsequent population expansion into the region. Despite a strong signal for fine-scale population genetic structure across the Puget Sound with as many as 8-10 distinct subpopulations supported by the SNP data, there is minimal evidence for morphological differentiation at this same spatiotemporal scale. Historical demographic analyses suggest that populations expanded and diverged across the region as the Cordilleran Ice Sheet receded. Population isolation, lack of dispersal corridors, and strict habitat requirements are the key drivers of population divergence in this system. These same factors may prove detrimental to the future persistence of populations as they cope with increasing shoreline development associated with urbanization.

Effects of species interactions on the potential for evolution at species' range limits

Species’ ranges are limited by both ecological and evolutionary constraints. While there is a growing appreciation that ecological constraints include interactions among species, like competition, we know relatively little about how interactions contribute to evolutionary constraints at species' niche and range limits. Building on concepts from community ecology and evolutionary biology, we review how biotic interactions can influence adaptation at range limits by impeding the demographic conditions that facilitate evolution (which we term a ‘demographic pathway to adaptation’), and/or by imposing evolutionary trade-offs with the abiotic environment (a ‘trade-offs pathway’). While theory for the former is well-developed, theory for the trade-offs pathway is not, and empirical evidence is scarce for both. Therefore, we develop a model to illustrate how fitness trade-offs along biotic and abiotic gradients could affect the potential for range expansion and niche evolution following ecological release. The model shows that which genotypes are favoured at species' range edges can depend strongly on the biotic context and the nature of fitness trade-offs. Experiments that characterize trade-offs and properly account for biotic context are needed to predict which species will expand their niche or range in response to environmental change.

This article is part of the theme issue ‘Species’ ranges in the face of changing environments (Part II)’.

Environmental variation and biotic interactions limit adaptation at ecological margins: lessons from rainforest Drosophila and European butterflies

Models of local adaptation to spatially varying selection predict that maximum rates of evolution are determined by the interaction between increased adaptive potential owing to increased genetic variation, and the cost genetic variation brings by reducing population fitness. We discuss existing and new results from our laboratory assays and field transplants of rainforest Drosophila and UK butterflies along environmental gradients, which try to test these predictions in natural populations. Our data suggest that: (i) local adaptation along ecological gradients is not consistently observed in time and space, especially where biotic and abiotic interactions affect both gradient steepness and genetic variation in fitness; (ii) genetic variation in fitness observed in the laboratory is only sometimes visible to selection in the field, suggesting that demographic costs can remain high without increasing adaptive potential; and (iii) antagonistic interactions between species reduce local productivity, especially at ecological margins. Such antagonistic interactions steepen gradients and may increase the cost of adaptation by increasing its dimensionality. However, where biotic interactions do evolve, rapid range expansion can follow. Future research should test how the environmental sensitivity of genotypes determines their ecological exposure, and its effects on genetic variation in fitness, to predict the probability of evolutionary rescue at ecological margins. This article is part of the theme issue 'Species' ranges in the face of changing environments (Part II)'.

Understanding the biology of species' ranges: when and how does evolution change the rules of ecological engagement?

Understanding processes that limit species' ranges has been a core issue in ecology and evolutionary biology for many decades, and has become increasingly important given the need to predict the responses of biological communities to rapid environmental change. However, we still have a poor understanding of evolution at range limits and its capacity to change the ecological 'rules of engagement' that define these communities, as well as the time frame over which this occurs. Here we link papers in the current volume to some key concepts involved in the interactions between evolutionary and ecological processes at species' margins. In particular, we separate hypotheses about species' margins that focus on hard evolutionary limits, which determine how genotypes interact with their environment, from those concerned with soft evolutionary limits, which determine where and when local adaptation can persist in space and time. We show how theoretical models and empirical studies highlight conditions under which gene flow can expand local limits as well as contain them. In doing so, we emphasize the complex interplay between selection, demography and population structure throughout a species' geographical and ecological range that determines its persistence in biological communities. However, despite some impressively detailed studies on range limits, particularly in invertebrates and plants, few generalizations have emerged that can predict evolutionary responses at ecological margins. We outline some directions for future work such as considering the impact of structural genetic variants and metapopulation structure on limits, and the interaction between range limits and the evolution of mating systems and non-random dispersal. This article is part of the theme issue 'Species' ranges in the face of changing environments (Part II)'.

Temporal variation may have diverse impacts on range limits

Environmental fluctuations are pervasive in nature, but the influence of non-directional temporal variation on range limits has received scant attention. We synthesize insights from the literature and use simple models to make conceptual points about the potentially wide range of ecological and evolutionary effects of temporal variation on range limits. Because organisms respond nonlinearly to environmental conditions, temporal variation can directionally alter long-term growth rates, either to shrink or to expand ranges. We illustrate this diversity of outcomes with a model of competition along a mortality gradient. Temporal variation can permit transitions between alternative states, potentially facilitating range expansion. We show this for variation in dispersal, using simple source-sink population models (with strong Allee effects, or with gene flow hampering local adaptation). Temporal variation enhances extinction risk owing to demographic stochasticity, rare events, and loss of genetic variation, all tending to shrink ranges. However, specific adaptations to exploit variation (including dispersal) may permit larger ranges than in similar but constant environments. Grappling with temporal variation is essential both to understand eco-evolutionary dynamics at range limits and to guide conservation and management strategies. This article is part of the theme issue 'Species' ranges in the face of changing environments (Part II)'.

Population genomic signatures of the oriental fruit moth related to the Pleistocene climates

The Quaternary climatic oscillations are expected to have had strong impacts on the evolution of species. Although legacies of the Quaternary climates on population processes have been widely identified in diverse groups of species, adaptive genetic changes shaped during the Quaternary have been harder to decipher. Here, we assembled a chromosome-level genome of the oriental fruit moth and compared genomic variation among refugial and colonized populations of this species that diverged in the Pleistocene. High genomic diversity was maintained in refugial populations. Demographic analysis showed that the effective population size of refugial populations declined during the penultimate glacial maximum (PGM) but remained stable during the last glacial maximum (LGM), indicating a strong impact of the PGM rather than the LGM on this pest species. Genome scans identified one chromosomal inversion and a mutation of the circadian gene Clk on the neo-Z chromosome potentially related to the endemicity of a refugial population. In the colonized populations, genes in pathways of energy metabolism and wing development showed signatures of selection. These different genomic signatures of refugial and colonized populations point to multiple impacts of Quaternary climates on adaptation in an extant species.

The oriental fruit moth is a pest species native to East Asia with refugial and colonized populations throughout the region. Here, a chromosome-level assembly for the species is reported and used to identify genomic signatures related to Quaternary climate change.

Disentangling key species interactions in diverse and heterogeneous communities: A Bayesian sparse modelling approach

Modelling species interactions in diverse communities traditionally requires a prohibitively large number of species‐interaction coefficients, especially when considering environmental dependence of parameters. We implemented Bayesian variable selection via sparsity‐inducing priors on non‐linear species abundance models to determine which species interactions should be retained and which can be represented as an average heterospecific interaction term, reducing the number of model parameters. We evaluated model performance using simulated communities, computing out‐of‐sample predictive accuracy and parameter recovery across different input sample sizes. We applied our method to a diverse empirical community, allowing us to disentangle the direct role of environmental gradients on species’ intrinsic growth rates from indirect effects via competitive interactions. We also identified a few neighbouring species from the diverse community that had non‐generic interactions with our focal species. This sparse modelling approach facilitates exploration of species interactions in diverse communities while maintaining a manageable number of parameters.

Competition, Trait Variance Dynamics, and the Evolution of a Species' Range

Geographic ranges of communities of species evolve in response to environmental, ecological, and evolutionary forces. Understanding the effects of these forces on species' range dynamics is a major goal of spatial ecology. Previous mathematical models have jointly captured the dynamic changes in species' population distributions and the selective evolution of fitness-related phenotypic traits in the presence of an environmental gradient. These models inevitably include some unrealistic assumptions, and biologically reasonable ranges of values for their parameters are not easy to specify. As a result, simulations of the seminal models of this type can lead to markedly different conclusions about the behavior of such populations, including the possibility of maladaptation setting stable range boundaries. Here, we harmonize such results by developing and simulating a continuum model of range evolution in a community of species that interact competitively while diffusing over an environmental gradient. Our model extends existing models by incorporating both competition and freely changing intraspecific trait variance. Simulations of this model predict a spatial profile of species' trait variance that is consistent with experimental measurements available in the literature. Moreover, they reaffirm interspecific competition as an effective factor in limiting species' ranges, even when trait variance is not artificially constrained. These theoretical results can inform the design of, as yet rare, empirical studies to clarify the evolutionary causes of range stabilization.

Population transcriptomics reveals the effect of gene flow on the evolution of range limits

One of the most important questions in evolutionary biology is how the spatial distribution of species is limited. Asymmetric gene flow from core populations is suggested to increase the number of poorly adapted immigrants in the populations at the range edge. Genetic load due to migration, i.e., migration load, should prevent adaptation to the local habitat, leading to decreases in distribution range via local extinction or the limiting range expansion. However, few experimental studies have examined the effects of immigration on fitness and natural selection within recipient populations. To investigate the influence of migration load on the evolution of distribution range, we performed field and laboratory observations as well as population transcriptomics for the common river snail, Semisulcospira reiniana. This species meets the conditions that migration from source populations can prevent local adaptation in a sink population because they inhabit the broader range of environments, including middle/upper reaches of a river and estuaries within a single river and they may be more vulnerable to being swept away by water currents due to lowered spontaneous (upward) locomotion activity. We found that river steepness was related to the lower distribution limit of S. reiniana, with a narrower distribution range in the steeper river. Population transcriptomic analysis showed that gene flow was heavily asymmetric from the upstream populations to downstream ones in the steep river, suggesting a greater migration load in the steep river. The number of genes putatively involved in adaptation to the local habitat was lower in the steep river than in the gentle river. Gene expression profiles suggested that individuals achieve better local adaptation in the gentle river. Laboratory experiments suggested that evolutionary differences in salinity tolerance among local populations were only found in the gentle river. Our results consistent with the hypothesis that migration load owing to asymmetric gene flow disturbs local adaptation and restricts the distribution range of river snails.

Effects of habitat quality on abundance, size and growth of mussel recruits

Recruitment of mussels is a complex process with the successful arrival of individuals hinging on the availability of suitable habitats. We examined the effects of adult mussels as settlement habitat and the degree to which the suitability of habitat they offer is species-specific by comparing the recruitment success of intertidal mussels. We hypothesised that mussel recruitment and early growth are dictated by the quality of habitat offered by conspecifics adults. We used a unique experimental arena on the south coast of South Africa, where Mytilus galloprovincialis and two lineages of Perna perna co-exist. Treatments were based on the translocation of individuals of M. galloprovincialis, western- and eastern lineage of P. perna to a single site, where artificial beds were created and sampled monthly over one year. Recruit's number, their sizes and growth were greater within beds of the western lineage of Perna than eastern lineage or Mytilus beds. The results clearly demonstrate that the quality of settlement habitat offered by adult beds differs among adult lineages/species and affects rates of settlement and the early growth of recruits. This effect extends to the intraspecific level; we found greater differences in density and growth of recruits between lineages of Perna than between either lineage and M. galloprovincialis.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10750-022-04994-7.

Prevalence and Adaptive Impact of Introgression

Alleles that introgress between species can influence the evolutionary and ecological fate of species exposed to novel environments. Hybrid offspring of different species are often unfit, and yet it has long been argued that introgression can be a potent force in evolution, especially in plants. Over the last two decades, genomic data have increasingly provided evidence that introgression is a critically important source of genetic variation and that this additional variation can be useful in adaptive evolution of both animals and plants. Here, we review factors that influence the probability that foreign genetic variants provide long-term benefits (so-called adaptive introgression) and discuss their potential benefits. We find that introgression plays an important role in adaptive evolution, particularly when a species is far from its fitness optimum, such as when they expand their range or are subject to changing environments.

Why Are Species’ Traits Weak Predictors of Range Shifts?

We examine the evidence linking species’ traits to contemporary range shifts and find they are poor predictors of range shifts that have occurred over decades to a century. We then discuss reasons for the poor performance of traits for describing interspecific variation in range shifts from two perspectives: ( a) factors associated with species’ traits that degrade range-shift signals stemming from the measures used for species’ traits, traits that are typically not analyzed, and the influence of phylogeny on range-shift potential and ( b) issues in quantifying range shifts and relating them to species’ traits due to imperfect detection of species, differences in the responses of altitudinal and latitudinal ranges, and emphasis on testing linear relationships between traits and range shifts instead of nonlinear responses. Improving trait-based approaches requires a recognition that traits within individuals interact in unexpected ways and that different combinations of traits may be functionally equivalent.

Genome-wide analysis reveals associations between climate and regional patterns of adaptive divergence and dispersal in American pikas

Understanding the role of adaptation in species' responses to climate change is important for evaluating the evolutionary potential of populations and informing conservation efforts. Population genomics provides a useful approach for identifying putative signatures of selection and the underlying environmental factors or biological processes that may be involved. Here, we employed a population genomic approach within a space-for-time study design to investigate the genetic basis of local adaptation and reconstruct patterns of movement across rapidly changing environments in a thermally sensitive mammal, the American pika (Ochotona princeps). Using genotypic data at 49,074 single-nucleotide polymorphisms (SNPs), we analyzed patterns of genome-wide diversity, structure, and migration along three independent elevational transects located at the northern extent (Tweedsmuir South Provincial Park, British Columbia, Canada) and core (North Cascades National Park, Washington, USA) of the Cascades lineage. We identified 899 robust outlier SNPs within- and among-transects. Of those annotated to genes with known function, many were linked with cellular processes related to climate stress including ATP-binding, ATP citrate synthase activity, ATPase activity, hormone activity, metal ion-binding, and protein-binding. Moreover, we detected evidence for contrasting patterns of directional migration along transects across geographic regions that suggest an increased propensity for American pikas to disperse among lower elevation populations at higher latitudes where environments are generally cooler. Ultimately, our data indicate that fine-scale demographic patterns and adaptive processes may vary among populations of American pikas, providing an important context for evaluating biotic responses to climate change in this species and other alpine-adapted mammals.

Rules of Plant Species Ranges: Applications for Conservation Strategies

Earth is changing rapidly and so are many plant species’ ranges. Here, we synthesize eco-evolutionary patterns found in plant range studies and how knowledge of species ranges can inform our understanding of species conservation in the face of global change. We discuss whether general biogeographic “rules” are reliable and how they can be used to develop adaptive conservation strategies of native plant species across their ranges. Rules considered include (1) factors that set species range limits and promote range shifts; (2) the impact of biotic interactions on species range limits; (3) patterns of abundance and adaptive properties across species ranges; (4) patterns of gene flow and their implications for genetic rescue, and (5) the relationship between range size and conservation risk. We conclude by summarizing and evaluating potential species range rules to inform future conservation and management decisions. We also outline areas of research to better understand the adaptive capacity of plants under environmental change and the properties that govern species ranges. We advise conservationists to extend their work to specifically consider peripheral and novel populations, with a particular emphasis on small ranges. Finally, we call for a global effort to identify, synthesize, and analyze prevailing patterns or rules in ecology to help speed conservation efforts.

A genetically based ecological trade-off contributes to setting a geographic range limit

Understanding the ecological factors that shape geographic range limits and the evolutionary constraints that prevent populations from adaptively evolving beyond these limits is an unresolved question. Here, we investigated why the euryhaline fish, Poecila reticulata, is confined to freshwater within its native range, despite being tolerant of brackish water. We hypothesised that competitive interactions with a close relative, Poecilia picta, in brackish water prevents P. reticulata from colonising brackish water. Using a combination of field transplant, common garden breeding, and laboratory behaviour experiments, we find support for this hypothesis, as P. reticulata are behaviourally subordinate and have lower survival in brackish water with P. picta. We also found a negative genetic correlation between P. reticulata growth in brackish water versus freshwater in the presence of P. picta, suggesting a genetically based trade-off between salinity tolerance and competitive ability could constrain adaptive evolution at the range limit.

Dispersal evolution in temporally variable environments: implications for plant range dynamics

Dispersal-the movement of an individual from the site of birth to a different site for reproduction-is an ecological and evolutionary driver of species ranges that shapes patterns of colonization, connectivity, gene flow, and adaptation. In plants, the traits that influence dispersal often vary within and among species, are heritable, and evolve in response to the fitness consequences of moving through heterogeneous landscapes. Spatial and temporal variation in the quality and quantity of habitat are important sources of selection on dispersal strategies across species ranges. While recent reviews have evaluated the interactions between spatial variation in habitat and dispersal dynamics, the extent to which geographic variation in temporal variability can also shape range-wide patterns in dispersal traits has not been synthesized. In this paper, we summarize key predictions from metapopulation models that evaluate how dispersal evolves in response to spatial and temporal habitat variability. Next, we compile empirical data that quantify temporal variability in plant demography and patterns of dispersal trait variation across species ranges to evaluate the hypothesis that higher temporal variability favors increased dispersal at plant range limits. We found some suggestive evidence supporting this hypothesis while more generally identifying a major gap in empirical work evaluating plant metapopulation dynamics across species ranges and geographic variation in dispersal traits. To address this gap, we propose several future research directions that would advance our understanding of the interplay between spatiotemporal variability and dispersal trait variation in shaping the dynamics of current and future species ranges.

Environmental Heterogeneity Leads to Spatial Differences in Genetic Diversity and Demographic Structure of Acer caudatifolium

Under climate fluctuation, species dispersal may be disturbed by terrain and local climate, resulting in uneven spatial-genetic structure. In addition, organisms at different latitudes may be differentially susceptible to climate change. Here, we tracked the seed dispersal of Acer caudatifolium using chloroplast DNA to explore the relationships of terrain and local climate heterogeneity with range shifts and demography in Taiwan. Our results showed that the extant populations have shifted upward and northward to the mountains since the Last Glacial Maximum. The distributional upshift of A. caudatifolium is in contrast to the downward expansion of its closest relative in Taiwan, A. morrisonense. The northern populations of A. caudatifolium have acquired multiple-source chlorotypes and harbor high genetic diversity. However, effective gene flow between the north and south is interrupted by topography, geographic distance, north-south differences in October rainfall, and other climate heterogeneities, blocking southward genetic rescue. In addition, winter monsoon-driven rainfall may cause regional differences in the phenological schedule, resulting in adaptive effects on the timing of range shift and the genetic draft of chlorotype distribution. Terrain, distance, and local climate also differentiate the northernmost populations from the others, supporting the previous taxonomic treatment of Acer kawakamii var. taitonmontanum as an independent variety.

Range edges of North American marine species are tracking temperature over decades

Understanding the dynamics of species range edges in the modern era is key to addressing fundamental biogeographic questions about abiotic and biotic drivers of species distributions. Range edges are where colonization and extirpation processes unfold, and so these dynamics are also important to understand for effective natural resource management and conservation. However, few studies to date have analyzed time series of range edge positions in the context of climate change, in part because range edges are difficult to detect. We first quantified positions for 165 range edges of marine fishes and invertebrates from three U.S. continental shelf regions using up to five decades of survey data and a spatiotemporal model to account for sampling and measurement variability. We then analyzed whether those range edges maintained their edge thermal niche-the temperatures found at the range edge position-over time. A large majority of range edges (88%) maintained either summer or winter temperature extremes at the range edge over the study period, and most maintained both (76%), although not all of those range edges shifted in space. However, we also found numerous range edges-particularly poleward edges and edges in the region that experienced the most warming-that did not shift at all, shifted further than predicted by temperature alone, or shifted opposite the direction expected, underscoring the multiplicity of factors that drive changes in range edge positions. This study suggests that range edges of temperate marine species have largely maintained the same edge thermal niche during periods of rapid change and provides a blueprint for testing whether and to what degree species range edges track temperature in general.

How is adaptive potential distributed within species ranges?

Quantitative genetic variation (QGV) represents a major component of adaptive potential and, if reduced toward range-edge populations, could prevent a species' expansion or adaptive response to rapid ecological change. It has been hypothesized that QGV will be lower at the range edge due to small populations-often the result of poor habitat quality-and potentially decreased gene flow. However, whether central populations are higher in QGV is unknown. We used a meta-analytic approach to test for a general QGV-range position relationship, including geographic and climatic distance from range centers. We identified 35 studies meeting our criteria, yielding nearly 1000 estimates of QGV (including broad-sense heritability, narrow-sense heritability, and evolvability) from 34 species. The relationship between QGV and distance from the geographic range or climatic niche center depended on the focal trait and how QGV was estimated. We found some evidence that QGV declines from geographic centers but that it increases toward niche edges; niche and geographic distances were uncorrelated. Nevertheless, few studies have compared QGV in both central and marginal regions or environments within the same species. We call for more research in this area and discuss potential research avenues related to adaptive potential in the context of global change.

Insights into invasive species from whole-genome resequencing

Studies of invasive species can simultaneously inform management strategies and quantify rapid evolution in the wild. The role of genomics in invasion science is increasingly recognised, and the growing availability of reference genomes for invasive species is paving the way for whole-genome resequencing studies in a wide range of systems. Here, we survey the literature to assess the application of whole-genome resequencing data in invasion biology. For some applications, such as the reconstruction of invasion routes in time and space, sequencing the whole genome of many individuals can increase the accuracy of existing methods. In other cases, population genomic approaches such as haplotype analysis can permit entirely new questions to be addressed and new technologies applied. To date whole-genome resequencing has only been used in a handful of invasive systems, but these studies have confirmed the importance of processes such as balancing selection and hybridization in allowing invasive species to reuse existing adaptations and rapidly overcome the challenges of a foreign ecosystem. The use of genomic data does not constitute a paradigm shift per se, but by leveraging new theory, tools, and technologies, population genomics can provide unprecedented insight into basic and applied aspects of invasion science.

Adaptation across geographic ranges is consistent with strong selection in marginal climates and legacies of range expansion

Every species experiences limits to its geographic distribution. Some evolutionary models predict that populations at range edges are less well adapted to their local environments due to drift, expansion load, or swamping gene flow from the range interior. Alternatively, populations near range edges might be uniquely adapted to marginal environments. In this study, we use a database of transplant studies that quantify performance at broad geographic scales to test how local adaptation, site quality, and population quality change from spatial and climatic range centers toward edges. We find that populations from poleward edges perform relatively poorly, both on average across all sites (15% lower population quality) and when compared to other populations at home (31% relative fitness disadvantage), consistent with these populations harboring high genetic load. Populations from equatorial edges also perform poorly on average (18% lower population quality) but, in contrast, outperform foreign populations (16% relative fitness advantage), suggesting that populations from equatorial edges have strongly adapted to unique environments. Finally, we find that populations from sites that are thermally extreme relative to the species' niche demonstrate strong local adaptation, regardless of their geographic position. Our findings indicate that both nonadaptive processes and adaptive evolution contribute to variation in adaptation across species' ranges.