Limits to evolution at range margins: when and why does adaptation fail?

Heatwave-Induced Paternal Effects Have Limited Adaptive Benefits in Offspring

As the threat of climate change and associated heatwaves grows, we need to understand how natural populations will respond. Inter-generational non-genetic inheritance may play a key role in rapid adaptation, but whether such mechanisms are truly adaptive and sufficient to protect wild populations is unclear. The contribution of paternal effects in particular is not fully understood, even though the male reproductive system may be highly sensitive to heatwaves. We used the zebrafish Danio rerio to investigate the effects of heatwaves on male fertility and assess potential adaptive benefits to their offspring in a number of large-scale heatwave experiments. Heatwave conditions had negative effects on male fertility by reducing gamete quality and fertilisation success, and we found indications of an adaptive effect on hatching in offspring produced by heatwave-exposed males. Our findings highlight the importance of including male and female fertility when determining species ability to cope with extreme conditions and suggest that parental effects provide limited adaptive benefits.

Pollinator shifts despite hybridisation in the Cape's hyperdiverse heathers (Erica, Ericaceae)

Interrogating the ecological and geographic factors that influence population divergence dynamics can reveal why some groups of organisms diversify more prolifically than others. One such group is the heathers (Erica, Ericaceae), the largest plant genus in the Cape Floristic Region. We study Erica abietina, a highly variable species complex with four subspecies differing in geographic range, habitat and pollination syndrome. We test for population differentiation, hybridisation, introgression and pollinator-driven divergence using genotyping-by-sequencing on samples across the entire distribution. We find five variably distinct genetic groups, with one subspecies comprising two independent lineages that are geographically isolated and occur on different soil types. Phylogenetic analysis suggests two independent shifts between bird and insect pollination, with accompanying genetic divergence. However, for one pair of populations with different pollinators, we uncover several individuals of hybrid origin at a site of sympatry. These results suggest that floral differentiation driven by divergent selection acts in concert with geographic isolation to maintain reproductive isolation and promote speciation. Our investigations reveal a highly dynamic system whose diversity has been shaped by a variety of interacting forces. We suggest that such a system could be a model for much of the diversification of the Cape flora.

Effects of extreme temperatures and recovery potential of Gongolaria barbata from a coastal lagoon in the northern Adriatic Sea: an ex situ approach

BACKGROUND AND AIMS

Globally, rising seawater temperatures contribute to the regression of marine macroalgal forests. Along the Istrian coastline (northern Adriatic), an isolated population of Gongolaria barbata persists in a coastal lagoon, representing one of the last marine macroalgal forests in the region. Our objective was to examine the impact of extreme temperatures on the morphology and physiology of G. barbata and test its potential for recovery after simulating marine heatwave (MHW) conditions.

METHODS

We explored the occurrence of marine heatwaves in southern Istria, adjacent to the study area, in addition to extreme temperatures inside the area itself. Subsequently, we performed a thermotolerance experiment, consisting of a stress and recovery phase, in which we exposed G. barbata thalli to four extreme (28, 30, 32 and 34 °C) and one favourable (18 °C) temperature. We monitored morphological and physiological responses.

KEY RESULTS

Our findings indicate a significant rise in frequency, duration and intensity of MHWs over decades on the southern Istrian coast. Experimental results show that G. barbata demonstrates potential for both morphological and physiological recovery after exposure to temperatures as high as 32 °C. However, exposure to 34 °C led to thallus decay, with limited ability to regenerate.

CONCLUSIONS

Our results show that G. barbata has a remarkable resilience to long-term exposure to extreme temperatures ≤32 °C and suggest that short-term exposure to temperatures beyond this, as currently recorded inside the lagoon, do not notably affect the physiology or morphology of local G. barbata. With more MHWs expected in the future, such an adapted population might represent an important donor suitable for future restoration activities along the Istrian coast. These results emphasize the resilience of this unique population, but also warn of the vulnerability of marine macroalgal forests to rising seawater temperatures in rapidly changing climatic conditions.

Mate-choice copying accelerates species range expansion

Mate-choice copying is a type of social learning in which females can change their mate preference after observing the choice of others. This behaviour can potentially affect population evolution and ecology, namely through increased dispersal and reduced local adaptation. Here, we simulated the effects of mate-choice copying in populations expanding across an environmental gradient to understand whether it can accelerate or retard the expansion process. Two mate-choice copying strategies were used: (i) when females target a single individual and (ii) when females target similar individuals. We also simulated cases where the male trait singled out by females with mate choice maps perfectly onto his genotype or is influenced by genotype-by-environment interactions. These rules have different effects on the results. When a trait is determined by genotype alone, populations where copier females target all similar males expand faster and the number of potential copiers increased. However, when preference is determined by genotype-by-environment interactions, populations where copier females target a single male had higher dispersal and also expand faster, but the potential number of copiers decreases. The results show that mate-choice copying can accelerate the expansion process, although its adaptiveness depends on the information animals use in different contexts.

Identifying the abiotic factors that determine the inland range limits of a mesic-adapted lizard species

For most species, the factors that determine geographical range limits are unknown. In mesic-adapted species, populations occurring near the edge of the species' distribution provide ideal study systems to investigate what limits distributional ranges. We aimed to identify the abiotic constraints that preclude an east-Australian mesic-adapted lizard (Lampropholis delicata) from occupying arid environments. We performed lizard surveys at sites spanning an elevation/aridity gradient (380-1070 m) and measured the prevalence of habitat features (logs, rocks, leaf litter, bare ground, solar radiation) in addition to hourly temperatures in a variety of microhabitats available to lizards. Species distribution models (SDM) were used to identify the macroclimatic variables limiting the species' distribution. At its inland range limit, L. delicata is associated with mesic high-elevation forests with complex microhabitat structures, which gradually decline in availability toward lower (and more arid) elevations where the species is absent. Moreover, L. delicata is absent from sites with a shallow leaf litter layer, in which daily temperatures exceed the species' thermal preference range, which we determined in a laboratory thermal gradient. In regards to macroclimate, SDM revealed that temperature seasonality is the primary variable predicting the species' distribution, suggesting that L. delicata avoids inland areas owing to their high annual thermal variability. By combining multiple lines of evidence, this research highlights that habitat and microclimate suitability-not solely macroclimate suitability-are important range-limiting factors for mesic ectotherms and should be incorporated in studies addressing range-limiting hypotheses.

Heatwave-like events affect drone production and brood-care behaviour in bumblebees

Climate change is currently considered one of the major threats to biodiversity and is associated with an increase in the frequency and intensity of extreme weather events, such as heatwaves. Heatwaves create acutely stressful conditions that may lead to disruption in the performance and survival of ecologically and economically important organisms, such as insect pollinators. In this study, we investigated the impact of simulated heatwaves on the performance of queenless microcolonies of Bombus terrestris audax under laboratory conditions. Our results indicate that heatwaves can have significant impacts on bumblebee performance. However, contrary to our expectations, exposure to heatwaves did not affect survival. Exposure to a mild 5-day heatwave (30-32 °C) resulted in increased offspring production compared to those exposed to an extreme heatwave (34-36 °C) and to the control group (24 °C). We also found that brood-care behaviours were impacted by the magnitude of the heatwave. Wing fanning occurred occasionally at temperatures of 30-32 °C, whereas at 34-36 °C the proportion of workers engaged in this thermoregulatory behaviour increased significantly. Our results provide insights into the effects of heatwaves on bumblebee colony performance and underscore the use of microcolonies as a valuable tool for studying the effects of extreme weather events. Future research, especially field-based studies replicating natural foraging conditions, is crucial to complement laboratory-based studies to comprehend how heatwaves compromise the performance of pollinators. Such studies may potentially help to identify those species more resilient to climate change, as well as those that are most vulnerable.

Effects of Climate Change on the Distribution of Prosthechea mariae (Orchidaceae) and within Protected Areas in Mexico

The impact of climate change on the distribution of native species in the Neotropics remains uncertain for most species. Prosthechea mariae is an endemic epiphytic orchid in Mexico, categorized as threatened. The objective of this study was to assess the effect of climate change on the natural distribution of P. mariae and the capacity of protected areas (PAs) to safeguard optimal environmental conditions for the species in the future. Historical records were obtained from herbaria collections and through field surveys. We utilized climate variables from WorldClim for the baseline scenario and for the 2050 period, using the general circulation models CCSM4 and CNRM-CM5 (RCP 4.5). Three sets of climate data were created for the distribution models, and multiple models were evaluated using the kuenm package. We found that the species is restricted to the eastern region of the country. The projections of future scenarios predict not only a substantial reduction in habitat but also an increase in habitat fragmentation. Ten PAs were found within the current distribution area of the species; in the future, the species could lose between 36% and 48% of its available habitat within these PAs. The results allowed for the identification of locations where climate change will have the most severe effects, and proposals for long-term conservation are addressed.

Rare, long-distance dispersal underpins genetic connectivity in the pink sea fan, Eunicella verrucosa

Characterizing patterns of genetic connectivity in marine species is of critical importance given the anthropogenic pressures placed on the marine environment. For sessile species, population connectivity can be shaped by many processes, such as pelagic larval duration, oceanographic boundaries and currents. This study combines restriction-site associated DNA sequencing (RADseq) and passive particle dispersal modelling to delineate patterns of population connectivity in the pink sea fan, Eunicella verrucosa, a temperate octocoral. Individuals were sampled from 20 sites covering most of the species' northeast Atlantic range, and a site in the northwest Mediterranean Sea to inform on connectivity across the Atlantic-Mediterranean transition. Using 7510 neutral SNPs, a geographic cline of genetic clusters was detected, partitioning into Ireland, Britain, France, Spain (Atlantic), and Portugal and Spain (Mediterranean). Evidence of significant inbreeding was detected at all sites, a finding not detected in a previous study of this species based on microsatellite loci. Genetic connectivity was characterized by an isolation by distance pattern (IBD) (r 2 = 0.78, p < 0.001), which persisted across the Mediterranean-Atlantic boundary. In contrast, exploration of ancestral population assignment using the program ADMIXTURE indicated genetic partitioning across the Bay of Biscay, which we suggest represents a natural break in the species' range, possibly linked to a lack of suitable habitat. As the pelagic larval duration (PLD) is unknown, passive particle dispersal simulations were run for 14 and 21 days. For both modelled PLDs, inter-annual variations in particle trajectories suggested that in a long-lived, sessile species, range-wide IBD is driven by rare, longer dispersal events that act to maintain gene flow. These results suggest that oceanographic patterns may facilitate range-wide stepping-stone genetic connectivity in E. verrucosa and highlight that both oceanography and natural breaks in a species' range should be considered in the designation of ecologically coherent MPA networks.

Enhanced and asymmetric signatures of hybridization at climatic margins: Evidence from closely related dioecious fig species

Hybridization plays a significant role in biological evolution. However, it is not clear whether ecological contingency differentially influences likelihood of hybridization, particularly at ecological margins where parental species may exhibit reduced fitnesses. Moreover, it is unknown whether future ecosystem change will increase the prevalence of hybridization. Ficus heterostyla and F. squamosa are closely related species co-distributed from southern Thailand to southwest China where hybridization, yielding viable seeds, has been documented. As a robust test of ecological factors driving hybridization, we investigated spatial hybridization signatures based on nuclear microsatellites from extensive population sampling across a widespread contact range. Both species showed high population differentiation and strong patterns of isolation by distance. Admixture estimates exposed asymmetric interspecific gene flow. Signatures of hybridization increase significantly towards higher latitude zones, peaking at the northern climatic margins. Geographic variation in reproductive phenology combined with ecologically challenging marginal habitats may promote this phenomenon. Our work is a first systematic evaluation of such patterns in a comprehensive, latitudinally-based clinal context, and indicates that tendency to hybridize appears strongly influenced by environmental conditions. Moreover, that future climate change scenarios will likely alter and possibly augment cases of hybridization at ecosystem scales.

Out of the ancient lake: Multiple riverine colonizations and diversification of the freshwater snails in the genus Semisulcospira around Lake Biwa

Ancient lakes are a hotspot of biodiversity. Freshwater species often experience spectacular species radiation after colonizing lakes from riverine habitats. Therefore, the relationship between the fauna of the ancient lakes and the surrounding riverine system has a special significance in understanding their origin and evolutionary history. The study of ancient lake species often focused on the lake colonization of riverine species. In contrast, far less attention has been placed on the reverse direction: the riverine colonization of the lake species, despite its importance in disentangling their complex evolutionary history. The freshwater snails in the genus Semisulcospira involve endemic groups that radiated in the ancient Lake Biwa. Using genetics and fossil records, we inferred that the ancestors of these lake-endemic Semisulcospira snails historically colonized the riverine habitats at least three times during the Middle Pleistocene. Each colonization resulted in the formation of a new lineage that was genetically and morphologically distinct from other lineages. Further, one of these colonizations was followed by hybridization with a cosmopolitan riverine species, which potentially facilitated the population persistence of the colonizers in the new environment. Despite their complex histories, all these colonizers were currently grouped within a single species, Semisulcospira kurodai, suggesting cryptic diversity in this species. This study highlights the significance of the riverine colonizations of the lake species to fully understand the diversification history of freshwater fauna in and around the ancient lakes.

Monitoring of species' genetic diversity in Europe varies greatly and overlooks potential climate change impacts

Genetic monitoring of populations currently attracts interest in the context of the Convention on Biological Diversity but needs long-term planning and investments. However, genetic diversity has been largely neglected in biodiversity monitoring, and when addressed, it is treated separately, detached from other conservation issues, such as habitat alteration due to climate change. We report an accounting of efforts to monitor population genetic diversity in Europe (genetic monitoring effort, GME), the evaluation of which can help guide future capacity building and collaboration towards areas most in need of expanded monitoring. Overlaying GME with areas where the ranges of selected species of conservation interest approach current and future climate niche limits helps identify whether GME coincides with anticipated climate change effects on biodiversity. Our analysis suggests that country area, financial resources and conservation policy influence GME, high values of which only partially match species' joint patterns of limits to suitable climatic conditions. Populations at trailing climatic niche margins probably hold genetic diversity that is important for adaptation to changing climate. Our results illuminate the need in Europe for expanded investment in genetic monitoring across climate gradients occupied by focal species, a need arguably greatest in southeastern European countries. This need could be met in part by expanding the European Union's Birds and Habitats Directives to fully address the conservation and monitoring of genetic diversity.

Adaptive constraints at the range edge of a widespread and expanding invasive plant

Identifying the factors that facilitate and limit invasive species' range expansion has both practical and theoretical importance, especially at the range edges. Here, we used reciprocal common garden experiments spanning the North/South and East/West range that include the North American core, intermediate and range edges of the globally invasive plant, Johnsongrass (Sorghum halepense) to investigate the interplay of climate, biotic interactions (i.e. competition) and patterns of adaptation. Our results suggest that the rapid range expansion of Johnsongrass into diverse environments across wide geographies occurred largely without local adaptation, but that further range expansion may be restricted by a fitness trade-off that limits population growth at the range edge. Interestingly, plant competition strongly dampened Johnsongrass growth but did not change the rank order performance of populations within a garden, though this varied among gardens (climates). Our findings highlight the importance of including the range edge when studying the range dynamics of invasive species, especially as we try to understand how invasive species will respond to accelerating global changes.

Local niches explain coexistence in environmentally-distinct contact zones between Western Mediterranean vipers

Species' ecological niches are frequently analysed to gain insights into how anthropogenic changes affect biodiversity. Coping with these changes often involves shifts in niche expression, which can disrupt local biotic interactions. Secondary contact zones, where competition and ecological segregation commonly occur, are ideal for studying the ecological factors influencing species' niches. In this study, we investigated the effect of climate and landscape factors on the ecological niches of two viper species, Vipera aspis and Vipera latastei, across three contact zones in northern Iberia, characterized by varying levels of landscape alteration. Using niche overlap tests, ecological niche models and spatial analyses we observed local variation in the expression of the species' niches across the three contact zones, resulting from the different abiotic and biotic conditions of each area. Rather than spatial niche segregation, we observed high niche overlap, suggesting niche convergence. A pattern of asymmetrical niche variation was identified in all contact zones, driven by species' climatic tolerances and the environmental conditions of each area. V. aspis generally exhibited a wider niche, except in the southernmost zone where the pure Mediterranean climate favored V. latastei. Human-induced landscape changes intensified niche asymmetry, by favoring the most generalist V. aspis over the specialist V. latastei, increasing habitat overlap, and likely competition. This study presents a comprehensive analysis of niche expression at range margins, anticipating a heightened impact of landscape changes in V. latastei. The methodological framework implemented here, and our findings, hold significant relevance for biodiversity management and conservation in human-impacted areas.

Landscape genomics analysis provides insights into future climate change-driven risk in rhesus macaque

Climate change significantly affects the suitability of wildlife habitats. Thus, understanding how animals adapt ecologically and genetically to climate change is important for targeted species protection. Rhesus macaques (Macaca mulatta) are widely distributed and multi-climatically adapted primates. This study explored how rhesus macaques adapt to climate change by integrating ecological and genetic methods and applying species distribution models (SDMs) and a gradient forest (GF) model. The findings suggested that temperature seasonality primarily affects habitat suitability and indicated that climate change will have a dramatic impact on macaque populations in the future. We also applied genotype-environment association (GEA) analyses and selection signature analyses to identify genes associated with climate change and provide possible explanations for the adaptation of rhesus macaques to climatic environments. The population genomics analyses suggested that the Taihang population has the highest genomic vulnerability with inbreeding and low heterozygosity. Combined with the higher ecological vulnerability, additional conservation strategies are required for this population under higher risk of climate change. Our work measured the impact of climate change and enabled the identification of populations that exhibit high vulnerability to severe climate change. Such information is useful for selecting populations of rhesus macaques as subject of long-term monitoring or evolutionary rescue under future climate change.

Integrative analysis of physiology and genomics provides insights into freeze tolerance adaptations of Acacia koa along an elevational cline

Natural selection for plant species in heterogeneous environments creates genetic variation for traits such as cold tolerance. While physiological or molecular analyses have been used to evaluate stress tolerance adaptations, combining these approaches may provide deeper insight. Acacia koa (koa) occurs from sea level to 2300 m in Hawai'i, USA. At high elevations, natural koa populations have declined due to deforestation, and freeze tolerance is a limiting factor for tree regeneration. We used physiology and molecular analyses to evaluate cold tolerance of koa populations from low (300-750 m), middle (750-1500 m), and high elevations (1500-2100 m). Half of the seedlings were cold acclimated by exposure to progressively lowered air temperatures for eight weeks (from 25.6/22.2°C to 8/4°C, day/night). Using the whole plant physiology-freezing test and koa C-repeat Binding Factor CBF genes, our results indicated that koa can be cold-acclimated when exposed to low, non-freezing temperatures. Seedlings from high elevations had consistently higher expression of Koa CBF genes associated with cold tolerance, helping to explain variation in cold-hardy phenotypes. Evaluation of the genetic background of 22 koa families across the elevations with low coverage RNA sequencing indicated that high elevation koa had relatively low values of heterozygosity, suggesting that adaptation is more likely to arise in the middle and low elevation sources. This physiology and molecular data for cold tolerance of koa across the elevation gradient of the Hawaiian Islands provides insights into natural selection processes and may help to support guidelines for conservation and seed transfer in forest restoration efforts.

Proteomic differences in seminal fluid of social insects whose sperm differ in heat tolerance

In the coming years, climate change is likely to increase the frequency and intensity of heatwaves. In many organisms, heat stress provokes physiological perturbations and can lead to decreased male fertility. Bumblebees are endo-heterothermic but display interspecific differences in thermotolerance that could have conservation implications. For the species of concern Bombus magnus, exposure to high temperatures can severely reduce sperm quality and, consequently, reproductive success. Such is not the case for B. terrestris, a ubiquitous species. To decipher the mechanisms at play, we characterized the seminal fluid proteomes of the two species. We quantified 1121 proteins, of which 522 were differentially expressed between B. terrestris and B. magnus. Several proteins with protective functions, such as proteases, antioxidant proteins and various heat-shock proteins, were present at higher levels in B. terrestris than in B. magnus under both control and heat-stress conditions. The same was true for proteins involved in cellular homeostasis, immunity, lipid/sugar metabolism and thermotolerance. Furthermore, proteins involved in the capture and elimination of reactive oxygen species also occurred at much high levels in B. terrestris. Overall, these results clearly indicate differences in the seminal proteome of the more thermotolerant B. terrestris versus B. magnus. The differences may contribute to explaining interspecific differences in sperm survival.

Genomic Signatures of Local Adaptation under High Gene Flow in Lumpfish-Implications for Broodstock Provenance Sourcing and Larval Production

Aquaculture of the lumpfish (Cyclopterus lumpus L.) has become a large, lucrative industry owing to the escalating demand for "cleaner fish" to minimise sea lice infestations in Atlantic salmon mariculture farms. We used over 10K genome-wide single nucleotide polymorphisms (SNPs) to investigate the spatial patterns of genomic variation in the lumpfish along the coast of Norway and across the North Atlantic. Moreover, we applied three genome scans for outliers and two genotype-environment association tests to assess the signatures and patterns of local adaptation under extensive gene flow. With our 'global' sampling regime, we found two major genetic groups of lumpfish, i.e., the western and eastern Atlantic. Regionally in Norway, we found marginal evidence of population structure, where the population genomic analysis revealed a small portion of individuals with a different genetic ancestry. Nevertheless, we found strong support for local adaption under high gene flow in the Norwegian lumpfish and identified over 380 high-confidence environment-associated loci linked to gene sets with a key role in biological processes associated with environmental pressures and embryonic development. Our results bridge population genetic/genomics studies with seascape genomics studies and will facilitate genome-enabled monitoring of the genetic impacts of escapees and allow for genetic-informed broodstock selection and management in Norway.

Widespread learned predator recognition to an alien predator across populations in an amphibian species

Alien predators are a major cause of decline and extinction of species worldwide, since native organisms are rarely equipped with specific antipredatory strategies to cope with them. However, phenotypic plasticity and learned predator recognition may help prey populations to survive novel predators. Here we examine geographical variation in the learning ability of larval spadefoot toads (Pelobates cultripes) to recognize invasive predatory crayfish (Procambarus clarkii). We compare the learning-mediated behavioural responses of tadpoles from six populations across two regions in Spain (central and southern), with different histories of exposure to the presence of the invasive species. Two of the populations showed innate recognition of chemical cues from the invasive crayfish, whereas three of them learned to recognize such cues as a threat after conditioning with conspecific alarm cues. Learning abilities did not differ among southern populations, but they did among central populations. We assessed patterns of genetic variation within and among these two regions through microsatellite markers and found low genetic divergence among the southern populations but greater differentiation among the central ones. We hypothesize that similar responses to the invasive crayfish in southern populations may have arisen from a combination of extended historical exposure to this introduced predator (~ 50 y) and higher levels of gene flow, as they inhabit a highly interconnected pond network. In contrast, populations from central Spain show lower connectivity, have been exposed to the invasive crayfish for a shorter period of time, and are more divergent in their plastic responses.

Hidden genetic variation in plasticity provides the potential for rapid adaptation to novel environments

Rapid environmental change is forcing populations into environments where plasticity will no longer maintain fitness. When populations are exposed to novel environments, evolutionary theory predicts that genetic variation in fitness will increase and should be associated with genetic differences in plasticity. If true, then genetic variation in plasticity can increase adaptive potential in novel environments, and population persistence via evolutionary rescue is more likely. To test whether genetic variation in fitness increases in novel environments and is associated with plasticity, we transplanted 8149 clones of 314 genotypes of a Sicilian daisy (Senecio chrysanthemifolius) within and outside its native range, and quantified genetic variation in fitness, and plasticity in leaf traits and gene expression. Although mean fitness declined by 87% in the novel environment, genetic variance in fitness increased threefold and was correlated with plasticity in leaf traits. High fitness genotypes showed greater plasticity in gene expression, but lower plasticity in most leaf traits. Interestingly, genotypes with the highest fitness in the novel environment had the lowest fitness at the native site. These results suggest that standing genetic variation in plasticity could help populations to persist and adapt to novel environments, despite remaining hidden in native environments.

Corals at the edge of environmental limits: A new conceptual framework to re-define marginal and extreme coral communities

The worldwide decline of coral reefs has renewed interest in coral communities at the edge of environmental limits because they have the potential to serve as resilience hotspots and climate change refugia, and can provide insights into how coral reefs might function in future ocean conditions. These coral communities are often referred to as marginal or extreme but few definitions exist and usage of these terms has therefore been inconsistent. This creates significant challenges for categorising these often poorly studied communities and synthesising data across locations. Furthermore, this impedes our understanding of how coral communities can persist at the edge of their environmental limits and the lessons they provide for future coral reef survival. Here, we propose that marginal and extreme coral communities are related but distinct and provide a novel conceptual framework to redefine them. Specifically, we define coral reef extremeness solely based on environmental conditions (i.e., large deviations from optimal conditions in terms of mean and/or variance) and marginality solely based on ecological criteria (i.e., altered community composition and/or ecosystem functioning). This joint but independent assessment of environmental and ecological criteria is critical to avoid common pitfalls where coral communities existing outside the presumed optimal conditions for coral reef development are automatically considered inferior to coral reefs in more traditional settings. We further evaluate the differential potential of marginal and extreme coral communities to serve as natural laboratories, resilience hotspots and climate change refugia, and discuss strategies for their conservation and management as well as priorities for future research. Our new classification framework provides an important tool to improve our understanding of how corals can persist at the edge of their environmental limits and how we can leverage this knowledge to optimise strategies for coral reef conservation, restoration and management in a rapidly changing ocean.

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.

Changes in Community Composition and Functional Traits of Bumblebees in an Alpine Ecosystem Relate to Climate Warming

Climate warming has been observed as the main cause of changes in diversity, community composition, and spatial distribution of different plant and invertebrate species. Due to even stronger warming compared to the global mean, bumblebees in alpine ecosystems are particularly exposed to these changes. To investigate the effects of climate warming, we sampled bumblebees along an elevational gradient, compared the records with data from 1935 and 1936, and related our results to climate models. We found that bumblebee community composition differed significantly between sampling periods and that increasing temperatures in spring were the most plausible factor explaining these range shifts. In addition, species diversity estimates were significantly lower compared to historical records. The number of socio-parasitic species was significantly higher in the historical communities, while recent communities showed increases in climate generalists and forest species at lower elevations. Nevertheless, no significant changes in community-weighted means of a species temperature index (STI) or the number of cold-adapted species were detected, likely due to the historical data resolution. We conclude that the composition and functionality of bumblebee communities in the study area have been significantly affected by climate warming, with changes in land use and vegetation cover likely playing an additional important role.

Are populations less genetically diverse and more differentiated at the species range edges? Analysis of the quantitative trait and molecular variation in wild oat Avena sterilis

Although understanding factors determining the genetic makeup of natural populations has long been an important goal of evolutionary biology, the effect of population position within the species range (i.e., interior vs. edge) on species genetic variation is still unclear. According to the ‘abundant center’ hypothesis, the range edge populations are expected to exhibit lower genetic variation and higher differentiation than core populations because of their greater spatial isolation and smaller size. We tested these predictions by comparing the extent and structure of neutral (SSR) and quantitative trait genetic variation in 20 populations of an annual grass Avena sterilis hierarchically sampled at the species range core and two opposite edges. Within-population genetic diversity was reduced at the desert range edge compared to the range core as assessed by six SSR markers but not by eight quantitative traits; no reduction was detected at the mountain edge. Genetic differentiation among populations was higher at the desert range edge than the range core in both molecular markers and quantitative traits, but not at the mountain edge. Our results imply that the pattern of population genetic variation at the species range edges largely depends on the steepness of the environmental cline that has a major effect on species fitness. The more gradual the environmental cline from the species interior towards the edge, the higher the probability of detecting reduced genetic diversity and increased differentiation of peripheral populations as predicted by the ‘abundant center’ hypothesis.

Parallel adaptation prompted core-periphery divergence of Ammopiptanthus mongolicus

Range expansion requires peripheral populations to shift adaptive optima to breach range boundaries. Opportunities for range expansion can be assessed by investigating the associations of core-periphery environmental and genetic differences. This study investigates differences in the core-periphery adaptation of Ammopiptanthus mongolicus, a broad-leaved evergreen shrub species in a relatively homogeneous temperate Asian desert environment, to explore the environmental factors that limit the expansion of desert plants. Temperate deserts are characterized by severe drought, a large diurnal temperature range, and seasonality. Long-standing adaptation to the harsh desert environment may confine the genetic diversity of A. mongolicus, despite its distribution over a wide range of longitude, latitude, and altitude. Since range edges defined by climate niches may have different genetic responses to environmental extremes, we compared genome-wide polymorphisms between nine environmental core populations and ten fragmented peripheral populations to determine the "adaptive peripheral" populations. At least four adaptive peripheral populations had similar genetic-environmental association patterns. High elevations, summer drought, and winter cold were the three main determinants of converging these four adaptive peripheral populations. Elevation mainly caused similar local climates among different geographic regions. Altitudinal adaptation resulting from integrated environmental-genetic responses was a breakthrough in breaching niche boundaries. These peripheral populations are also located in relatively humid and warmer environments. Relaxation of the drought and cold constraints facilitated the genetic divergence of these peripheral populations from the core population's adaptive legacy. We conclude that pleiotropic selection synchronized adaptative divergence to cold and drought vs. warm and humid environments between the core and peripheral populations. Such parallel adaptation of peripheral populations relies on selection under a background of abundant new variants derived from the core population's standing genetic variation, i.e., integration of genetic surfing and local adaptation.

Genetic variability and the ecology of geographic range: A test of the central-marginal hypothesis in Australian scincid lizards

For many species, both local abundance and regional occupancy are highest near the centre of their geographic distributions. One hypothesis for this pattern is that niche suitability declines with increasing distance from a species geographic centre, such that populations near range margins are characterized by reduced density and increased patchiness. In these smaller edge populations, genetic drift is more powerful, leading to the loss of genetic diversity. This simple verbal model has been formalized as the central-marginal hypothesis, which predicts that core populations should have greater genetic diversity than edge populations. Here, we tested the central-marginal hypothesis using a genomic data set of 25 species-level taxa of Australian scincid lizards in the genera Ctenotus and Lerista. A majority of taxa in our data set showed range-wide patterns of genetic variation consistent with central-marginal hypothesis, and eight of 25 taxa showed significantly greater genetic diversity in the centre of their range. We then explored biological, historical, and methodological factors that might predict which taxa support the central-marginal hypothesis. We found that taxa with the strongest evidence for range expansion were the least likely to follow predictions of the central-marginal hypothesis. The majority of these taxa had range expansions that originated at the range edge, which led to a gradient of decreasing genetic diversity from the range edge to the core, contrary to the central-marginal hypothesis.

Mutation in phcA Enhanced the Adaptation of Ralstonia solanacearum to Long-Term Acid Stress

Bacterial wilt, caused by the plant pathogen Ralstonia solanacearum, occurs more severely in acidified soil according to previous reports. However, R. solanacearum cannot grow well in acidic environments under barren nutrient culture conditions, especially when the pH is lower than 5. With the worsening acidification of farmland, further determination of how R. solanacearum adapts to the long-term acidic environment is worthwhile. In this study, experimental evolution was applied to evaluate the adaptability and mechanism of the R. solanacearum experimental population responding to long-term acid stress. We chose the CQPS-1 strain as the ancestor, and minimal medium (MM medium) with different pH values as the culture environment to simulate poor soil. After 1500 generations of serial passage experiments in pH 4.9 MM, acid-adapted experimental strains (denoted as C49 strains) were obtained, showing significantly higher growth rates than the growth rates of control experimental strains (serial passage experiment in pH 6.5 MM, denoted as C65 strains). Competition experiments showed that the competitive indices (CIs) of all selected clones from C49 strains were superior to the ancestor in acidic environment competitiveness. Based on the genome variation analysis and functional verification, we confirmed that loss of function in the phcA gene was associated with the acid fitness gain of R. solanacearum, which meant that the inactivation of the PhcA regulator caused by gene mutation mediated the population expansion of R. solanacearum when growing in an acidic stress environment. Moreover, the swimming motility of acid evolution strains and the phcA deletion mutant was significantly enhanced compared to CQPS-1. This work provided evidence for understanding the adaptive strategy of R. solanacearum to the long-term acidic environment.

Macroevolutionary dynamics of climatic niche space

How and why lineages evolve along with niche space as they diversify and adapt to different environments is fundamental to evolution. Progress has been hampered by the difficulties of linking a robust empirical characterization of species niches with flexible evolutionary models that describe their evolution. Consequently, the relative influence of abiotic and biotic factors remains poorly understood. Here, we characterize species' two-dimensional temperature and precipitation niche space occupied (i.e. species niche envelope) as complex geometries and assess their evolution across all Aves using a model that captures heterogeneous evolutionary rates on time-calibrated phylogenies. We find that extant birds coevolved from warm, mesic climatic niches into colder and drier environments and responded to the Cretaceous-Palaeogene (K-Pg) boundary with a dramatic increase in disparity. Contrary to expectations of subsiding rates of niche evolution, our results show that overall rates have increased steadily, with some lineages experiencing exceptionally high evolutionary rates, associated with the colonization of novel niche spaces, and others showing niche stasis. Both competition- and environmental change-driven niche evolution transpire and result in highly heterogeneous rates near the present. Our findings highlight the growing ecological and conservation insights arising from the model-based integration of comprehensive environmental and phylogenetic information.

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)'.

Sympatry leads to reduced body condition in chickadees that occasionally hybridize

Both abiotic and biotic drivers influence species distributions. Abiotic drivers such as climate have received considerable attention, even though biotic drivers such as hybridization often interact with abiotic drivers. We sought to explore the (1) costs of co-occurrence for ecologically similar species that hybridize and (2) associations between ecological factors and condition to understand how abiotic and biotic factors influence species distributions. For two closely related and ecologically similar songbirds, black-capped and mountain chickadees, we characterized body condition, as a proxy for fitness, using a 1358-individual range-wide dataset. We compared body condition in sympatry and allopatry with several abiotic and biotic factors using species-specific generalized linear mixed models. We generated genomic data for a subset of 217 individuals to determine the extent of hybridization-driven admixture in our dataset. Within this data subset, we found that ~11% of the chickadees had hybrid ancestry, and all hybrid individuals had typical black-capped chickadee plumage. In the full dataset, we found that birds of both species, independent of demographic and abiotic factors, had significantly lower body condition when occurring in sympatry than birds in allopatry. This could be driven by either the inclusion of cryptic, likely poor condition, hybrids in our full dataset, competitive interactions in sympatry, or range edge effects. We are currently unable to discriminate between these mechanisms. Our findings have implications for mountain chickadees in particular, which will encounter more black-capped chickadees as black-capped chickadee ranges shift upslope and could lead to local declines in mountain chickadee populations.

Variation in intraspecific demography drives localised concordance but species-wide discordance in response to past climatic change

BACKGROUND

Understanding how species biology may facilitate resilience to climate change remains a critical factor in detecting and protecting species at risk of extinction. Many studies have focused on the role of particular ecological traits in driving species responses, but less so on demographic history and levels of standing genetic variation. Additionally, spatial variation in the interaction of demographic and adaptive factors may further complicate prediction of species responses to environmental change. We used environmental and genomic datasets to reconstruct the phylogeographic histories of two ecologically similar and largely co-distributed freshwater fishes, the southern (Nannoperca australis) and Yarra (N. obscura) pygmy perches, to assess the degree of concordance in their responses to Plio-Pleistocene climatic changes. We described contemporary genetic diversity, phylogenetic histories, demographic histories, and historical species distributions across both species, and statistically evaluated the degree of concordance in co-occurring populations.

RESULTS

Marked differences in contemporary genetic diversity, historical distribution changes and historical migration were observed across the species, with a distinct lack of genetic diversity and historical range expansion suggested for N. obscura. Although several co-occurring populations within a shared climatic refugium demonstrated concordant demographic histories, idiosyncratic population size changes were found at the range edges of the more spatially restricted species. Discordant responses between species were associated with low standing genetic variation in peripheral populations. This might have hindered adaptive potential, as documented in recent demographic declines and population extinctions for the two species.

CONCLUSION

Our results highlight both the role of spatial scale in the degree of concordance in species responses to climate change, and the importance of standing genetic variation in facilitating range shifts. Even when ecological traits are similar between species, long-term genetic diversity and historical population demography may lead to discordant responses to ongoing and future climate change.

Genetic load and extinction in peripheral populations: the roles of migration, drift and demographic stochasticity

We analyse how migration from a large mainland influences genetic load and population numbers on an island, in a scenario where fitness-affecting variants are unconditionally deleterious, and where numbers decline with increasing load. Our analysis shows that migration can have qualitatively different effects, depending on the total mutation target and fitness effects of deleterious variants. In particular, we find that populations exhibit a genetic Allee effect across a wide range of parameter combinations, when variants are partially recessive, cycling between low-load (large-population) and high-load (sink) states. Increased migration reduces load in the sink state (by increasing heterozygosity) but further inflates load in the large-population state (by hindering purging). We identify various critical parameter thresholds at which one or other stable state collapses, and discuss how these thresholds are influenced by the genetic versus demographic effects of migration. Our analysis is based on a 'semi-deterministic' analysis, which accounts for genetic drift but neglects demographic stochasticity. We also compare against simulations which account for both demographic stochasticity and drift. Our results clarify the importance of gene flow as a key determinant of extinction risk in peripheral populations, even in the absence of ecological gradients. This article is part of the theme issue 'Species' ranges in the face of changing environments (part I)'.

The role of phenotypic plasticity in the establishment of range margins

It has been argued that adaptive phenotypic plasticity may facilitate range expansions over spatially and temporally variable environments. However, plasticity may induce fitness costs. This may hinder the evolution of plasticity. Earlier modelling studies examined the role of plasticity during range expansions of populations with fixed genetic variance. However, genetic variance evolves in natural populations. This may critically alter model outcomes. We ask: how does the capacity for plasticity in populations with evolving genetic variance alter range margins that populations without the capacity for plasticity are expected to attain? We answered this question using computer simulations and analytical approximations. We found a critical plasticity cost above which the capacity for plasticity has no impact on the expected range of the population. Below the critical cost, by contrast, plasticity facilitates range expansion, extending the range in comparison to that expected for populations without plasticity. We further found that populations may evolve plasticity to buffer temporal environmental fluctuations, but only when the plasticity cost is below the critical cost. Thus, the cost of plasticity is a key factor involved in range expansions of populations with the potential to express plastic response in the adaptive trait. This article is part of the theme issue 'Species' ranges in the face of changing environments (part I)'.

Variation in insulative feather structure in songbirds replacing each other along a tropical elevation gradient

High-elevation organisms are expected to evolve physiological adaptations to cope with harsh environmental conditions. Yet, evidence for such adaptive differences, especially compared to closely related lowland taxa occurring along the same elevational gradient, is rare. Revisiting an anecdotal natural history observation by O. Bangs from 1899 and based on new measurements of museum specimens, we confirmed that the high-elevation hermit wood wren (Henicorhina anachoreta) from the Sierra Nevada de Santa Marta, Colombia, has longer, more insulative feathers on the chest and back, than its lower-elevation counterpart the grey-breasted wood wren (H. leucophrys). However, we did not find evidence for the same specializations in subspecies of H. leucophrys that live at high elevations on other elevational gradients in the Colombian Andes, although similar adaptive solutions have arisen in separate mountain systems like the Himalayas. Adaptations in plumage may be associated with the recurrence of elevational species replacements throughout the tropics.

New records of Darevskia praticola at the northern limit of its distribution range in Romania

In the summer of 2021 we identified three new distribution localities of Darevskia praticola north of the Mureș River, and one locality south of the river. The habitats populated by D. praticola (broad-leaved forest with wet areas) and the altitude (175–245 m) of the new records are typical for this species. Nevertheless, D. praticola had not been recorded in 12 other localities with similar conditions from an area previously considered suitable for this species. Thus, D. praticola may be slowly expanding from a bridgehead north of the Mureș River, occupying new favorable habitats. Probably, D. praticola recently crossed the Mureș River, possibly on a bridge, or with the timber trucks which exploit the woods from both sides of the river.

Interactions between microenvironment, selection and genetic architecture drive multiscale adaptation in a simulation experiment

When environmental conditions differ both within and among populations, multiscale adaptation results from processes at both scales and interference across scales. We hypothesize that within-population environmental heterogeneity influences the chance of success of migration events, both within and among populations, and maintains within-population adaptive differentiation. We used a simulation approach to analyze the joint effects of environmental heterogeneity patterns, selection intensity and number of QTL controlling a selected trait on local adaptation in a hierarchical metapopulation design. We show the general effects of within-population environmental heterogeneity: (i) it increases occupancy rate at the margins of distribution ranges, under extreme environments and high levels of selection; (ii) it increases the adaptation lag in all environments; (iii) it impacts the genetic variance in each environment, depending on the ratio of within- to between-populations environmental heterogeneity; (iv) it reduces the selection-induced erosion of adaptive gene diversity. Most often, the smaller the number of QTL involved, the stronger are these effects. We also show that both within- and between-populations phenotypic differentiation (Q_ST ) mainly results from covariance of QTL effects rather than QTL differentiation (F_STq ), that within-population QTL differentiation is negligible, and that stronger divergent selection is required to produce adaptive differentiation within populations than among populations. With a high number of QTL, when the difference between environments within populations exceeds the smallest difference between environments across populations, high levels of within-population differentiation can be reached, reducing differentiation among populations. Our study stresses the need to account for within-population environmental heterogeneity when investigating local adaptation.

Into the range: a latitudinal gradient or a center-margins differentiation of ecological strategies in Arabidopsis thaliana?

BACKGROUND AND AIMS

Determining within-species large-scale variation in phenotypic traits is central to elucidate the drivers of species' ranges. Intraspecific comparisons offer the opportunity to understand how trade-offs and biogeographical history constrain adaptation to contrasted environmental conditions. Here we test whether functional traits, ecological strategies from the CSR scheme and phenotypic plasticity in response to abiotic stress vary along a latitudinal or a center- margins gradient within the native range of Arabidopsis thaliana.

METHODS

We experimentally examined the phenotypic outcomes of plant adaptation at the center and margins of its geographic range using 30 accessions from southern, central and northern Europe. We characterized the variation of traits related to stress tolerance, resource use, colonization ability, CSR strategy scores, survival and fecundity in response to high temperature (34 °C) or frost (- 6 °C), combined with a water deficit treatment.

KEY RESULTS

We found evidence for both a latitudinal and a center-margins differentiation for the traits under scrutiny. Age at maturity, leaf dry matter content, specific leaf area and leaf nitrogen content varied along a latitudinal gradient. Northern accessions presented a greater survival to stress than central and southern accessions. Leaf area, C-scores, R-scores and fruit number followed a center-margins differentiation. Central accessions displayed a higher phenotypic plasticity than northern and southern accessions for most studied traits.

CONCLUSIONS

Traits related to an acquisitive/conservative resource-use trade-off followed a latitudinal gradient. Traits associated with a competition/colonization trade-off differentiated along the historic colonization of the distribution range and then followed a center-margins differentiation. Our findings pinpoint the need to consider the joint effect of evolutionary history and environmental factors when examining phenotypic variation across the distribution range of a species.

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.

Indications for rapid evolution of trait means and thermal plasticity in range-expanding populations of a butterfly

Currently, poleward range expansions are observed in many taxa, often in response to anthropogenic climate change. At the expanding front, populations likely face cooler and more variable temperature conditions, imposing thermal selection. This may result in changes in trait means or plasticity, the relative contribution of which is not well understood. We, here, investigate evolutionary change in range-expanding populations of the butterfly Pieris mannii, by comparing populations from the core and the newly established northern range under laboratory conditions. We observed both changes in trait means and in thermal reaction norms. Range-expanding populations showed a more rapid development, potentially indicative of counter-gradient variation and an increased cold tolerance compared with core populations. Genotype-environment interactions prevailed in all associated traits, such that the above differences were restricted to cooler environmental conditions. In range-expanding populations, plasticity was decreased in developmental traits enabling relatively rapid growth even under cooler conditions but increased in cold tolerance arguably promoting higher activity under thermally challenging conditions. Notably, these changes must have occurred within a time period of ca. 10 years only. Our results suggest, in line with contemporary theory, that the evolution of plasticity may play a hitherto underestimated role for adaptation to climatic variation. However, rather than generally increased or decreased levels of plasticity, our results indicate fine-tuned, trait-specific evolutionary responses to increase fitness in novel environments.

Body and wing size, but not wing shape, vary along a large-scale latitudinal gradient in a damselfly

Large-scale latitudinal studies that include both north and south edge populations and address sex differences are needed to understand how selection has shaped trait variation. We quantified the variation of flight-related morphological traits (body size, wing size, ratio between wing size and body size, and wing shape) along the whole latitudinal distribution of the damselfly Lestes sponsa, spanning over 2700 km. We tested predictions of geographic variation in the flight-related traits as a signature of: (1) stronger natural selection to improve dispersal in males and females at edge populations; (2) stronger sexual selection to improve reproduction (fecundity in females and sexual behaviors in males) at edge populations. We found that body size and wing size showed a U-shaped latitudinal pattern, while wing ratio showed the inverse shape. However, wing shape varied very little along the latitudinal gradient. We also detected sex-differences in the latitudinal patterns of variation. We discuss how latitudinal differences in natural and sexual selection regimes can lead to the observed quadratic patterns of variation in body and wing morphology via direct or indirect selection. We also discuss the lack of latitudinal variation in wing shape, possibly due to aerodynamic constraints.

Balanced polymorphism fuels rapid selection in an invasive crab despite high gene flow and low genetic diversity

Adaptation across environmental gradients has been demonstrated in numerous systems with extensive dispersal, despite high gene flow and consequently low genetic structure. The speed and mechanisms by which such adaptation occurs remain poorly resolved, but are critical to understanding species spread and persistence in a changing world. Here, we investigate these mechanisms in the European green crab Carcinus maenas, a globally distributed invader. We focus on a northwestern Pacific population that spread across >12 degrees of latitude in 10 years from a single source, following its introduction <35 years ago. Using six locations spanning >1500 km, we examine genetic structure using 9376 single nucleotide polymorphisms (SNPs). We find high connectivity among five locations, with significant structure between these locations and an enclosed lagoon with limited connectivity to the coast. Among the five highly connected locations, the only structure observed was a cline driven by a handful of SNPs strongly associated with latitude and winter temperature. These SNPs are almost exclusively found in a large cluster of genes in strong linkage disequilibrium that was previously identified as a candidate for cold tolerance adaptation in this species. This region may represent a balanced polymorphism that evolved to promote rapid adaptation in variable environments despite high gene flow, and which now contributes to successful invasion and spread in a novel environment. This research suggests an answer to the paradox of genetically depauperate yet successful invaders: populations may be able to adapt via a few variants of large effect despite low overall diversity.

Demographic history shapes genomic ancestry in hybrid zones

Demographic factors such as migration rate and population size can impede or facilitate speciation. In hybrid zones, reproductive boundaries between species are tested and demography mediates the opportunity for admixture between lineages that are partially isolated. Genomic ancestry is a powerful tool for revealing the history of admixed populations, but models and methods based on local ancestry are rarely applied to structured hybrid zones. To understand the effects of demography on ancestry in hybrids zones, we performed individual-based simulations under a stepping-stone model, treating migration rate, deme size, and hybrid zone age as parameters. We find that the number of ancestry junctions (the transition points between genomic regions with different ancestries) and heterogenicity (the genomic proportion heterozygous for ancestry) are often closely connected to demographic history. Reducing deme size reduces junction number and heterogenicity. Elevating migration rate increases heterogenicity, but migration affects junction number in more complex ways. We highlight the junction frequency spectrum as a novel and informative summary of ancestry that responds to demographic history. A substantial proportion of junctions are expected to fix when migration is limited or deme size is small, changing the shape of the spectrum. Our findings suggest that genomic patterns of ancestry could be used to infer demographic history in hybrid zones.