Using isolation-by-distance to jointly estimate effective population density and dispersal distance: a practical evaluation using bumble bees

Effective population density and intergenerational dispersal distance are key aspects of population biology, but obtaining empirical estimates of these parameters can be difficult. This is especially true for my study taxa, wild bees. In this paper, I apply and evaluate an existing but underutilized method to estimate effective density and dispersal distance of bumble bees (Bombus, Apidae). Specifically, using 10 datasets of bumble bees in North America, I use the relationship between genetic isolation-by-distance and Wright's neighborhood size to define a density-dispersal isocline-that is, a curve describing pairs of density and dispersal values consistent with observed rates of isolation-by-distance. These parameters are inversely related; as one increases the other decreases. I then use outside estimates of bumble bee dispersal distances to make more specific estimates of effective colony density. Compared to some prior estimates of census density (100s to 1000s colonies/km2), my estimated effective colony densities were very low (1-41 effective colonies/km2). I also hypothesize, however, that these estimates are affected by the spatial extent of sampling, due to scale-dependent patterns in the distribution of individuals. To test this hypothesis, I subsampled each dataset to simulate varying study extent, and repeated my analysis. Within populations, effective densities tended to decrease when measured across larger spatial extents. Altogether, I demonstrate a useful and under-appreciated tool for studying population biology, especially of small, mobile animals like bees, but also show that researchers must interpret their results carefully within the context of their study design.

Long-distance seasonal migration to the tropics promotes genetic diversity but not gene flow in boreal birds

Differences in life history can cause co-distributed species to evolve contrasting population genetic patterns, even as they occupy the same landscape. In high-latitude animals, evolutionary processes may be especially influenced by long-distance seasonal migration, a widespread adaptation to seasonality. Although migratory movements are intuitively linked to dispersal and therefore promotion of gene flow, their evolutionary genetic consequences remain poorly understood. Using ~1,700 genomes from 35 co-distributed boreal-breeding bird species that differ in non-breeding latitude and thus migration distance, we find that most long-distance migrants unexpectedly exhibit spatial genetic structure, despite their strong movement propensity. This result suggests evolutionary effects of philopatry-the tendency of many migrants to return to the same breeding site year after year, resulting in restricted dispersal. We further demonstrate that migration distance and genetic diversity are strongly positively correlated in our study species. This striking relationship suggests that the adaptive seasonal shifts in biogeography inherent to long-distance migration may enhance population stability, preserving genetic diversity in long-distance migrants relative to shorter-distance migrants that winter in harsher conditions at higher latitudes. Our results suggest that the major impact of long-distance seasonal migration on population genetic evolution occurs through promotion of demographic stability, rather than facilitation of dispersal.

Phosphorus acquisition capacity and size trait evolution in Achillea wilhelmsii reflect adaptation to environmental gradients

Understanding how plant traits evolve in response to environmental gradients is critical for elucidating mechanisms of local adaptation. This study investigated trait variation in Achillea wilhelmsii accessions from eight climatically diverse Iranian locations after cultivation under uniform conditions. Key traits-plant size, biomass, and phosphorus content-reflect adaptive divergence along an integrated environmental gradient (PCA1) derived from temperature, precipitation, and altitude. Warmer/drier, lower-altitude conditions corresponded with reduced phosphorus uptake and smaller plant size, while cooler/wetter, higher-altitude conditions favored increased phosphorus absorption and larger plant size. Principal component analyses revealed that 62.85% of observed trait variation arises from evolutionary responses via genetic divergence, driven by natural selection across environmental gradients. Populations from colder, high-altitude sites (positive PCA1 scores) evolved enhanced phosphorus uptake, likely due to selection on pre-existing genetic variation for cold tolerance, facilitating larger plant sizes. These patterns highlight how adaptation, as an inherent capacity of plants to respond to selection pressures, shapes trait divergence under environmental heterogeneity. Future studies should dissect the genetic architecture linking phosphorus metabolism, environmental gradients, and plant size evolution.

Population Genetic Structure of Three-Spined Sticklebacks in the St. Lawrence: A Gradient of Change

Understanding how environmental gradients shape population genetic structure is critical for elucidating evolutionary dynamics in heterogeneous landscapes. The St. Lawrence Estuary, spanning fluvial, middle, and marine zones, presents a steep salinity gradient that serves as an ideal setting to study such a question. Three-spined sticklebacks (Gasterosteus aculeatus) thrive across these zones, offering an ideal model system to investigate the interplay of gene flow and natural selection in shaping population structure. Using whole-genome resequencing of sticklebacks from 12 sites, this study aimed to resolve fine-scale population structure and investigate how genetic diversity and differentiation are influenced by selection and gene flow. By integrating single nucleotide polymorphisms (SNPs) and structural variants (SVs), we assessed differentiation patterns, examined clinal variation, and evaluated the relative roles of gene flow and selection in shaping population dynamics. Our findings reveal clear genetic differentiation between fluvial and saltwater populations, with Baie-Saint-Paul forming a potential third group. Salinity emerged as a key driver of genetic structure, with clinal variation in allele frequencies suggesting ongoing adaptation along the gradient. Demographic modeling indicated a history of secondary contact with recent and weak gene flow. Structural variants, particularly indels, complemented SNP-based analyses, underscoring their importance in detecting fine-scale population structure. These results highlight the complex interplay of evolutionary forces shaping biodiversity in transitional environments, providing a basis for exploring local adaptation in connected populations and contributing to broader efforts in conservation genomics.

Selection Over Small and Large Spatial Scales in the Face of High Gene Flow

Local adaptation represents the balance of selection and gene flow. Increasingly, studies find that adaptation can occur on spatial scales much smaller than the scale of dispersal, resulting in balanced polymorphisms within populations. However, microgeographic adaptation might be facilitated or hindered by large-scale environmental heterogeneity, such as across latitude. Marine systems present a special case, as many marine species have high dispersal capacity so that dispersal 'neighbourhoods' may encompass environmental heterogeneity over both small and large spatial scales. Here, we leverage fine-scale sampling across the California range of the Pacific purple sea urchin (Strongylocentrotus purpuratus), a species with previous evidence of both local adaptation and extremely high gene flow. We find that despite the complete absence of neutral population structure, satellite-based sea surface temperature and tidal zone are associated with subtle genetic differences among populations, suggesting that balanced polymorphisms can lead to adaptation across both large (latitudinal) and small (subtidal vs. intertidal) scales. In fact, some of the same genetic variants differentiate populations at both spatial scales, potentially because both environmental parameters are related to temperature. Further, we find that genes that are expressed at a single tissue or life history stage are more divergent than expected across both latitudinal and tidal zone comparisons, suggesting that these genes have specific functions that might generate phenotypic variation important for local adaptation. Together, these results suggest that even in species with little population structure, genetic variation can be sorted across varying spatial scales, potentially resulting in local adaptation across complex environmental mosaics.

Local Adaptation Is Highest in Populations With Stable Long-Term Growth

Theory suggests that the drivers of demographic variation and local adaptation are shared and may feedback on one other. Despite some evidence for these links in controlled settings, the relationship between local adaptation and demography remains largely unexplored in natural conditions. Using 10 years of demographic data and two reciprocal transplant experiments, we tested predictions about the relationship between the magnitude of local adaptation and demographic variation (population growth rates and their elasticities to vital rates) across 10 populations of a well-studied annual plant. In both years, we found a strong unimodal relationship between mean home-away local adaptation and stochastic population growth rates. Other predicted links were either weakly or not supported by our data. Our results suggest that declining and rapidly growing populations exhibit reduced local adaptation, potentially due to maladaptation and relaxed selection, respectively.

Intraspecific variation of thermal tolerance along elevational gradients: the case of a widespread diving beetle (Coleoptera: Dytiscidae)

Species distributed across wide elevational gradients are likely to experience local thermal adaptation and exhibit high thermal plasticity, as these gradients are characterised by steep environmental changes over short geographic distances (i.e., strong selection differentials). The prevalence of adaptive intraspecific variation in thermal tolerance with elevation remains unclear, however, particularly in freshwater taxa. We explored variation in upper and lower thermal limits and acclimation capacity among Iberian populations of adults of the widespread water beetle Agabus bipustulatus (Dytiscidae) across a 2000 m elevational gradient, from lowland to alpine areas. Since mean and extreme temperatures decline with elevation, we predicted that populations at higher elevations will show lower heat tolerance and higher cold tolerance than lowland ones. We also explored whether acclimation capacity is positively related with climatic variability across elevations. We found significant variation in thermal limits between populations of A. bipustulatus, but no evidence of local adaptation to different thermal conditions across the altitudinal gradient, as relationships between thermal limits and elevation or climatic variables were largely nonsignificant. Furthermore, plasticities of both upper and lower thermal limits were consistently low in all populations. These results suggest thermal niche conservatism in this species, likely due to gene flow counteracting the effects of divergent selection, or adaptations in other traits that buffer exposure to climate extremes. The limited adaptive potential and plasticity of thermal tolerance observed in A. bipustulatus suggest that even generalist species, distributed across wide environmental gradients, may have limited resilience to global warming.

Patterns of genetic variation and local adaptation of a native herbivore to a lethal invasive plant

Understanding the evolutionary processes that influence fitness is critical to predicting species' responses to selection. Interactions among evolutionary processes including gene flow, drift and the strength of selection can lead to either local adaptation or maladaptation, especially in heterogenous landscapes. Populations experiencing novel environments or resources are ideal for understanding the mechanisms underlying adaptation or maladaptation, specifically in locally co-evolved interactions. We used the interaction between a native herbivore that oviposits on a patchily distributed introduced plant that in turn causes significant mortality to the larvae to test for signatures of local adaptation in areas where the two co-occurred. We used whole-genome sequencing to explore population structure, patterns of gene flow and signatures of local adaptation. We found signatures of local adaptation in response to the introduced plant in the absence of strong population structure with no genetic differentiation and low genetic variation. Additionally, we found localized allele frequency differences within a single population between habitats with and without the lethal plant, highlighting the effects of strong selection. Finally, we identified that selection was acting on larval ability to feed on the plant rather than on females' ability to avoid oviposition, thus uncovering the specific ontogenetic target of selection. Our work highlights the potential for adaptation to occur in a fine-grained landscape in the presence of gene flow and low genetic variation.

Environmental gradients mediate dispersal evolution during biological invasions

Rapid evolution of increased dispersal at the edge of a range expansion can accelerate invasions. However, populations expanding across environmental gradients often face challenging environments that reduce fitness of dispersing individuals. We used an eco-evolutionary model to explore how environmental gradients influence dispersal evolution and, in turn, modulate the speed and predictability of invasion. Environmental gradients opposed evolution of increased dispersal during invasion, even leading to evolution of reduced dispersal along steeper gradients. Counterintuitively, reduced dispersal could allow for faster expansion by minimizing maladaptive gene flow and facilitating adaptation. While dispersal evolution across homogenous landscapes increased both the mean and variance of expansion speed, these increases were greatly dampened by environmental gradients. We illustrate our model's potential application to prediction and management of invasions by parameterizing it with data from a recent invertebrate range expansion. Overall, we find that environmental gradients strongly modulate the effect of dispersal evolution on invasion trajectories.

Weak evidence of provenance effects in spring phenology across Europe and North America

Forecasting the biological impacts of climate change requires understanding how species respond to warmer temperatures through interannual flexible variation vs through adaptation to local conditions. Yet, we often lack this information entirely or find conflicting evidence across studies, which is the case for spring phenology. We synthesized common garden studies across Europe and North America that reported spring event dates for a mix of angiosperm and gymnosperm tree species in the northern hemisphere, capturing data from 384 North American and 101 European provenances (i.e. populations) with observations from 1962 to 2019, alongside autumn event data when provided. Across continents, we found no evidence of provenance effects in spring phenology, but strong clines with latitude and mean annual temperature in autumn. These effects, however, appeared to diverge by continent and species type (gymnosperm vs angiosperm), with particularly pronounced clines in North America in autumn events. Our results suggest flexible, likely plastic responses, in spring phenology with warming, and potential limits - at least in the short term - due to provenance effects for autumn phenology. They also highlight that, after over 250 yr of common garden studies on tree phenology, we still lack a holistic predictive model of clines across species and phenological events.

Sensitivity of habitat network models to changes in maximum dispersal distance

Predicting the presence or absence (occurrence-state) of species in a certain area is highly important for conservation. Occurrence-state can be assessed by network models that take suitable habitat patches as nodes, connected by potential dispersal of species. To determine connections, a connectivity threshold is set at the species' maximum dispersal distance. However, this requires field observations prone to underestimation, so for most animal species there are no trustable maximum dispersal distance estimations. This limits the development of accurate network models to predict species occurrence-state. In this study, we performed a sensitivity analysis of the performance of network models to different settings of maximum dispersal distance. Our approach, applied on six amphibian species in Switzerland, used habitat suitability modelling to define habitat patches, which were linked within a dispersal distance threshold to form habitat networks. We used network topological measures, patch suitability, and patch size to explain species occurrence-state in habitat patches through boosted regression trees. These modelling steps were repeated on each species for different maximum dispersal distances, including a species-specific value from literature. We evaluated mainly the predictive performance and predictor importance among the network models. We found that predictive performance had a positive relation with the distance threshold, and that almost none of the species-specific values from literature yielded the best performance across tested thresholds. With increasing dispersal distance, the importance of the habitat-quality-related variable decreased, whereas that of the topology-related predictors increased. We conclude that the sensitivity of these models to the dispersal distance parameter stems from the very different topologies formed with different movement assumptions. Most reported maximum dispersal distances are underestimated, presumably due to leptokurtic dispersal distribution. Our results imply that caution should be taken when selecting a dispersal distance threshold, considering higher values than those derived from field reports, to account for long-distance dispersers.

Temporal Variation in Selection Influences Microgeographic Local Adaptation

AbstractEcological heterogeneity can lead to local adaptation when populations exhibit fitness trade-offs among habitats. However, the degree to which local adaptation is affected by the spatial and temporal scale of environmental variation is poorly understood. A multiyear reciprocal transplant experiment was performed with populations of the annual plant Leptosiphon parviflorus living on adjacent serpentine and nonserpentine soil. Local adaptation over this small geographic scale was observed, but there were differences in the temporal variability of selection across habitats. On serpentine soil, the local population had a consistently large survival advantage, presumably as a result of the temporal stability in selection imposed by soil cation content. In contrast, a fecundity advantage was observed for the sandstone population on its native soil type but only in the two study years with the highest rainfall. A manipulative greenhouse experiment demonstrated that the fitness advantage of the sandstone population in its native soil type depends critically on water availability. The temporal variability in local adaptation driven by variation in precipitation suggests that continued drought conditions have the potential to erode local adaptation in these populations. These results show how different selective factors can influence spatial and temporal patterns of variation in fitness trade-offs.

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.

Evidence for ecotone speciation across an African rainforest-savanna gradient

Accelerating climate change and habitat loss make it imperative that plans to conserve biodiversity consider species' ability to adapt to changing environments. However, in biomes where biodiversity is highest, the evolutionary mechanisms responsible for generating adaptative variation and, ultimately, new species are frequently poorly understood. African rainforests represent one such biome, as decadal debates continue concerning the mechanisms generating African rainforest biodiversity. These debates hinge on the relative importance of geographic isolation versus divergent natural selection across environmental gradients. Hindering progress is a lack of robust tests of these competing hypotheses. Because African rainforests are severely at-risk due to climate change and other anthropogenic activities, addressing this long-standing debate is critical for making informed conservation decisions. We use demographic inference and allele frequency-environment relationships to investigate mechanisms of diversification in an African rainforest skink, Trachylepis affinis, a species inhabiting the gradient between rainforest and rainforest-savanna mosaic (ecotone). We provide compelling evidence of ecotone speciation, in which gene flow has all but ceased between rainforest and ecotone populations, at a level consistent with infrequent hybridization between sister species. Parallel patterns of genomic, morphological, and physiological divergence across this environmental gradient and pronounced allele frequency-environment correlation indicate speciation is mostly probably driven by ecological divergence, supporting a central role for divergent natural selection. Our results provide strong evidence for the importance of ecological gradients in African rainforest speciation and inform conservation strategies that preserve the processes that produce and maintain biodiversity.

A biome-dependent distribution gradient of tree species range edges is strongly dictated by climate spatial heterogeneity

Understanding the causes of the arrest of species distributions has been a fundamental question in ecology and evolution. These questions are of particular interest for trees owing to their long lifespan and sessile nature. A surge in data availability evokes a macro-ecological analysis to determine the underlying forces limiting distributions. Here we analyse the spatial distribution of >3,600 major tree species to determine geographical areas of range-edge hotspots and find drivers for their arrest. We confirmed biome edges to be strong delineators of distributions. Importantly, we identified a stronger contribution of temperate than tropical biomes to range edges, adding strength to the notion that tropical areas are centres of radiation. We subsequently identified a strong association of range-edge hotspots with steep spatial climatic gradients. We linked spatial and temporal homogeneity and high potential evapotranspiration in the tropics as the strongest predictors of this phenomenon. We propose that the poleward migration of species in light of climate change might be hindered because of steep climatic gradients.

Repeated divergence of amphibians and reptiles across an elevational gradient in northern Madagascar

How environmental factors shape patterns of biotic diversity in tropical ecosystems is an active field of research, but studies examining the possibility of ecological speciation in terrestrial tropical ecosystems are scarce. We use the isolated rainforest herpetofauna on the Montagne d'Ambre (Amber Mountain) massif in northern Madagascar as a model to explore elevational divergence at the level of populations and communities. Based on intensive sampling and DNA barcoding of amphibians and reptiles along a transect ranging from ca. 470-1470 m above sea level (a.s.l.), we assessed a main peak in species richness at an elevation of ca. 1000 m a.s.l. with 41 species. The proportion of local endemics was highest (about 1/3) at elevations >1100 m a.s.l. Two species of chameleons (Brookesia tuberculata, Calumma linotum) and two species of frogs (Mantidactylus bellyi, M. ambony) studied in depth by newly developed microsatellite markers showed genetic divergence up the slope of the mountain, some quite strong, others very weak, but in each case with genetic breaks between 1100 and 1270 m a.s.l. Genetic clusters were found in transect sections significantly differing in bioclimate and herpetological community composition. A decrease in body size was detected in several species with increasing elevation. The studied rainforest amphibians and reptiles show concordant population genetic differentiation across elevation along with morphological and niche differentiation. Whether this parapatric or microallopatric differentiation will suffice for the completion of speciation is, however, unclear, and available phylogeographic evidence rather suggests that a complex interplay between ecological and allopatric divergence processes is involved in generating the extraordinary species diversity of Madagascar's biota. Our study reveals concordant patterns of diversification among main elevational bands, but suggests that these adaptational processes are only part of the complex of processes leading to species formation, among which geographical isolation is probably also important.

Phenology and plasticity can prevent adaptive clines in thermal tolerance across temperate mountains: The importance of the elevation-time axis

Critical thermal limits (CTmax and CTmin) decrease with elevation, with greater change in CTmin, and the risk to suffer heat and cold stress increasing at the gradient ends. A central prediction is that populations will adapt to the prevailing climatic conditions. Yet, reliable support for such expectation is scant because of the complexity of integrating phenotypic, molecular divergence and organism exposure. We examined intraspecific variation of CTmax and CTmin, neutral variation for 11 microsatellite loci, and micro- and macro-temperatures in larvae from 11 populations of the Galician common frog (Rana parvipalmata) across an elevational gradient, to assess (1) the existence of local adaptation through a PST-FST comparison, (2) the acclimation scope in both thermal limits, and (3) the vulnerability to suffer acute heat and cold thermal stress, measured at both macro- and microclimatic scales. Our study revealed significant microgeographic variation in CTmax and CTmin, and unexpected elevation gradients in pond temperatures. However, variation in CTmax and CTmin could not be attributed to selection because critical thermal limits were not correlated to elevation or temperatures. Differences in breeding phenology among populations resulted in exposure to higher and more variable temperatures at mid and high elevations. Accordingly, mid- and high-elevation populations had higher CTmax and CTmin plasticities than lowland populations, but not more extreme CTmax and CTmin. Thus, our results support the prediction that plasticity and phenological shifts may hinder local adaptation, promoting thermal niche conservatism. This may simply be a consequence of a coupled variation of reproductive timing with elevation (the "elevation-time axis" for temperature variation). Mid and high mountain populations of R. parvipalmata are more vulnerable to heat and cool impacts than lowland populations during the aquatic phase. All of this contradicts some of the existing predictions on adaptive thermal clines and vulnerability to climate change in elevational gradients.

Urban evolution of thermal physiology in a range-expanding, mycophagous fruit fly, Drosophila tripunctata

In Drosophila spp., their often high number of annual generations, large population sizes and large amounts of standing genetic variation should predispose them to undergo contemporary adaptation to climatic warming. However, a number of laboratory experimental evolution studies in this group of organisms suggest strong limits on the rate and magnitude of contemporary thermal adaptation. Here, we explore this discrepancy by examining the potential for rapid evolutionary divergence between wild populations of Drosophila tripunctata Loew, 1862 from rural and urban sites. We performed a multi-generation common garden study and found evidence for the evolution of higher heat tolerance (critical thermal maximum) in flies from urban populations. We also detected evolutionary divergence in cold resistance (chill coma recovery time), with diminished cold resistance in flies from urban populations, although the effect was weaker than the shift in heat tolerance. Our study provides evidence of contemporary urban thermal adaptation, although the magnitude of phenotypic change lagged the magnitude of environmental temperature change across the urbanization gradient, suggesting potential limits on the evolution of urban thermal physiology.

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.

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.

Plastic and quantitative genetic divergence mirror environmental gradients among wild, fragmented populations of Impatiens capensis

PREMISE

Habitat fragmentation generates molecular genetic divergence among isolated populations, but few studies have assessed phenotypic divergence and fitness in populations where the genetic consequences of habitat fragmentation are known. Phenotypic divergence could reflect plasticity, local adaptation, and/or genetic drift.

METHODS

We examined patterns and potential drivers of phenotypic divergence among 12 populations of jewelweed (Impatiens capensis) that show strong molecular genetic signals of isolation and drift among fragmented habitats. We measured morphological and reproductive traits in both maternal plants within natural populations and their self-fertilized progeny grown together in a common garden. We also quantified environmental divergence between home sites and the common garden.

RESULTS

Populations with less molecular genetic variation expressed less maternal phenotypic variation. Progeny in the common garden converged in phenotypes relative to their wild mothers but retained among-population differences in morphology, survival, and reproduction. Among-population phenotypic variance was 3-10× greater in home sites than in the common garden for 6 of 7 morphological traits measured. Patterns of phenotypic divergence paralleled environmental gradients in ways suggestive of adaptation. Progeny resembled their mothers less as the environmental distance between their home site and the common garden increased.

CONCLUSIONS

Despite strong molecular signatures of isolation and drift, phenotypic differences among these Impatiens populations appear to reflect both adaptive quantitative genetic divergence and plasticity. Quantifying the extent of local adaptation and plasticity and how these covary with molecular and phenotypic variation help us predict when populations may lose their adaptive capacity.

Evolution in Cities

Although research performed in cities will not uncover new evolutionary mechanisms, it could provide unprecedented opportunities to examine the interplay of evolutionary forces in new ways and new avenues to address classic questions. However, while the variation within and among cities affords many opportunities to advance evolutionary biology research, careful alignment between how cities are used and the research questions being asked is necessary to maximize the insights that can be gained. In this review, we develop a framework to help guide alignment between urban evolution research approaches and questions. Using this framework, we highlight what has been accomplished to date in the field of urban evolution and identify several up-and-coming research directions for further expansion. We conclude that urban environments can be used as evolutionary test beds to tackle both new and long-standing questions in evolutionary biology.

Selection on growth rate and local adaptation drive genomic adaptation during experimental range expansions in the protist Tetrahymena thermophila

Populations that expand their range can undergo rapid evolutionary adaptation of life‐history traits, dispersal behaviour and adaptation to the local environment. Such adaptation may be aided or hindered by sexual reproduction, depending on the context.

However, few empirical and experimental studies have investigated the genetic basis of adaptive evolution during range expansions. Even less attention has been given to the question how sexual reproduction may modulate such adaptive evolution during range expansions.

We here studied genomic adaptation during experimental range expansions of the protist Tetrahymena thermophila in landscapes with a uniform environment or a pH gradient. Specifically, we investigated two aspects of genomic adaptation during range expansion. First, we investigated adaptive genetic change in terms of the underlying numbers of allele frequency changes from standing genetic variation and de novo variants. We focused on how sexual reproduction may alter this adaptive genetic change. Second, we identified genes subject to selection caused by the expanding range itself, and directional selection due to the presence or absence of the pH gradient. We focused this analysis on alleles with large frequency changes that occurred in parallel in more than one population to identify the most likely candidate targets of selection.

We found that sexual reproduction altered adaptive genetic change both in terms of de novo variants and standing genetic variation. However, sexual reproduction affected allele frequency changes in standing genetic variation only in the absence of long‐distance gene flow. Adaptation to the range expansion affected genes involved in cell divisions and DNA repair, whereas adaptation to the pH gradient additionally affected genes involved in ion balance and oxidoreductase reactions. These genetic changes may result from selection on growth and adaptation to low pH.

In the absence of gene flow, sexual reproduction may have aided genetic adaptation. Gene flow may have swamped expanding populations with maladapted alleles, thus reducing the extent of evolutionary adaptation during range expansion. Sexual reproduction also altered the genetic basis of adaptation in our evolving populations via de novo variants, possibly by purging deleterious mutations or by revealing fitness benefits of rare genetic variants.

Evidence for genetic isolation and local adaptation in the field cricket Gryllus campestris

Understanding how species can thrive in a range of environments is a central challenge for evolutionary ecology. There is strong evidence for local adaptation along large-scale ecological clines in insects. However, potential adaptation among neighbouring populations differing in their environment has been studied much less. We used RAD sequencing to quantify genetic divergence and clustering of ten populations of the field cricket Gryllus campestris in the Cantabrian Mountains of northern Spain, and an outgroup on the inland plain. Our populations were chosen to represent replicate high and low altitude habitats. We identified genetic clusters that include both high and low altitude populations indicating that the two habitat types do not hold ancestrally distinct lineages. Using common-garden rearing experiments to remove environmental effects, we found evidence for differences between high and low altitude populations in physiological and life-history traits. As predicted by the local adaptation hypothesis, crickets with parents from cooler (high altitude) populations recovered from periods of extreme cooling more rapidly than those with parents from warmer (low altitude) populations. Growth rates also differed between offspring from high and low altitude populations. However, contrary to our prediction that crickets from high altitudes would grow faster, the most striking difference was that at high temperatures, growth was fastest in individuals from low altitudes. Our findings reveal that populations a few tens of kilometres apart have independently evolved adaptations to their environment. This suggests that local adaptation in a range of traits may be commonplace even in mobile invertebrates at scales of a small fraction of species' distributions.

Local Adaptation in Marine Foundation Species at Microgeographic Scales

AbstractNearshore foundation species in coastal and estuarine systems (e.g., salt marsh grasses, mangroves, seagrasses, corals) drive the ecological functions of ecosystems and entire biomes by creating physical structure that alters local abiotic conditions and influences species interactions and composition. The resilience of foundation species and the ecosystem functions they provide depends on their phenotypic and genetic responses to spatial and temporal shifts in environmental conditions. In this review, we explore what is known about the causes and consequences of adaptive genetic differentiation in marine foundation species over spatial scales shorter than dispersal capabilities (i.e., microgeographic scales). We describe the strength of coupling field and laboratory experiments with population genetic techniques to illuminate patterns of local adaptation, and we illustrate this approach by using several foundation species. Among the major themes that emerge from our review include (1) adaptive differentiation of marine foundation species repeatedly evolves along vertical (i.e., elevation or depth) gradients, and (2) mating system and phenology may facilitate this differentiation. Microgeographic adaptation is an understudied mechanism potentially underpinning the resilience of many sessile marine species, and this evolutionary mechanism likely has particularly important consequences for the ecosystem functions provided by foundation species.

Does Phenological Plasticity Help or Hinder Range Shifts Under Climate Change?

Climate warming is predicted to shift species’ ranges as previously uninhabitable environments just beyond the leading range edges become suitable habitat and trailing range edges become increasingly unsuitable. Understanding which aspects of the environment and species traits mediate these range shifts is critical for understanding species’ possible redistributions under global change, yet we have a limited understanding of the ecological and evolutionary responses underlying population spread or extinction at species’ range edges. Within plant populations, shifts in flowering phenology have been one of the strongest and most consistent responses to climate change, and are likely to play an important role in mediating population dynamics within and beyond species’ ranges. However, the role of phenological shifts, and particularly phenological plasticity, in species’ range shifts remains relatively unstudied. Here, we synthesize literature on phenology, plasticity, and adaptation to suggest ways in which phenological responses to climate may vary across species’ ranges and review the empirical evidence for and against these hypotheses. We then outline how phenological plasticity could facilitate or hinder persistence and potential consequences of phenological plasticity in range expansions, including phenological cues, shifts in correlated traits, altered species interactions, and effects on gene flow. Finally, we suggest future avenues for research, such as characterizing reaction norms for phenology across a species’ range and in beyond-the-range transplant experiments. Given the prevalence and magnitude of phenological shifts, future work should carefully dissect its costs and benefits for population persistence, and incorporate phenological plasticity into models predicting species’ persistence and geographic range shifts under climate change.

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.

The Effect of the Recombination Rate between Adaptive Loci on the Capacity of a Population to Expand Its Range

AbstractPrevious theoretical work on range expansions over heterogeneous environments showed that there is a critical environmental gradient where range expansion stops. For populations with freely recombining loci underlying the trait under selection (hereafter, "adaptive loci"), the critical gradient in one-dimensional habitats depends on the fitness cost of dispersal and the strength of selection relative to genetic drift. Here, we extend the previous work in two directions and ask, What is the role of the recombination rate between the adaptive loci during range expansions? And what effect does the ability of selfing as opposed to obligate outcrossing have on range expansions? To answer these questions, we use computer simulations. We demonstrate that while reduced recombination rates between adaptive loci slow down range expansions as a result of poor purging of locally deleterious alleles at the expansion front, they may also allow a species to occupy a greater range. In addition, we find that the allowance of selfing may improve the ability of populations to expand their ranges, for example, because selfing among potentially rare high-fitness individuals facilitates the establishment and maintenance of locally well-adapted genotypes. We conclude that during range expansions there is a trade-off between positive and negative effects of recombination within and between individuals.

Scale-dependent signatures of local adaptation in a foundation tree species

Understanding local adaptation is critical for conservation management under rapidly changing environmental conditions. Local adaptation inferred from genotype-environment associations may show different genomic patterns depending on the spatial scale of sampling, due to differences in the slope of environmental gradients and the level of gene flow. We compared signatures of local adaptation across the genome of mountain ash (Eucalyptus regnans) at two spatial scales: A species-wide data set and a topographically-complex subregional data set. We genotyped 367 individual trees at over 3700 single-nucleotide polymorphisms (SNPs), quantified patterns of spatial genetic structure among populations, and used two analytical methods to identify loci associated with at least one of three environmental variables at each spatial scale. Together, the analyses identified 549 potentially adaptive SNPs at the subregion scale, and 435 SNPs at the range-wide scale. A total of 39 genic or near-genic SNPs, associated with 28 genes, were identified at both spatial scales, although no SNP was identified by both methods at both scales. We observed that nongenic regions had significantly higher homozygote excess than genic regions, possibly due to selective elimination of inbred genotypes during stand development. Our results suggest that strong environmental selection occurs in mountain ash, and that the identification of putatively adaptive loci can differ substantially depending on the spatial scale of analyses. We also highlight the importance of multiple adaptive genetic architectures for understanding patterns of local adaptation across large heterogenous landscapes, with comparison of putatively adaptive loci among spatial scales providing crucial insights into the process of adaptation.

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

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

Microgeographical adaptation corresponds to elevational distributions of congeneric montane grasshoppers

Local adaptation can occur at small spatial scales relative to the dispersal capacity of species. Alpine ecosystems have sharp environmental clines that offer an opportunity to investigate the effects of fine-scale shifts in species' niche breadth on adaptive genetic processes. Here we examine two grasshopper species endemic to the Australian Alps (Kosciuscola spp.) that differ in elevational niche breadth: one broader, K. usitatus (1400-2200 m), and one narrower, K. tristis (1600-2000 m). We examine signatures of selection with respect to environmental and morphological variables in two mountain regions using F_ST outlier tests and environmental association analyses (EAAs) applied to single nucleotide polymorphism (SNP) data (K. usitatus: 9017 SNPs, n = 130; K. tristis: 7363 SNPs, n = 135). Stronger genetic structure was found in the more narrowly distributed K. tristis, which showed almost twice the number of SNPs under putative selection (10.8%) compared with K. usitatus (5.3%). When examining SNPs in common across species (n = 3058), 260 SNPs (8.5%) were outliers shared across species, and these were mostly associated with elevation, a proxy for temperature, suggesting parallel adaptive processes in response to climatic drivers. Additive polygenic scores (an estimate of the cumulative signal of selection across all candidate loci) were nonlinearly and positively correlated with elevation in both species. However, a steeper correlation in K. tristis indicated a stronger signal of spatially varying selection towards higher elevations. Our study illustrates that the niche breadth of co-occurring and related species distributed along the same environmental cline is associated with differences in patterns of microgeographical adaptation.

What Do We Really Know About Adaptation at Range Edges?

Recent theory and empirical evidence have provided new insights regarding how evolutionary forces interact to shape adaptation at stable and transient range margins. Predictions regarding trait divergence at leading edges are frequently supported. However, declines in fitness at and beyond edges show that trait divergence has sometimes been insufficient to maintain high fitness, so identifying constraints to adaptation at range edges remains a key challenge. Indirect evidence suggests that range expansion may be limited by adaptive genetic variation, but direct estimates of genetic constraints at and beyond range edges are still scarce. Sequence data suggest increased genetic load in edge populations in several systems, but its causes and fitness consequences are usually poorly understood. The balance between maladaptive and positive effects of gene flow on fitness at range edges deserves further study. It is becoming increasingly clear that characterizations about degree of adaptation based solely on geographical peripherality are unsupported.

Adaptation to elevation but limited local adaptation in an amphibian

We performed a reciprocal transplant experiment to estimate "parallel" adaptation to elevation and "unique" adaptation to local sites at the same elevation, using the frog Rana temporaria in the Swiss Alps. It is important to distinguish these two processes because they have different implications for population structure and ecological specialization. Larvae were reared from hatching to metamorphosis within enclosures installed in their pond of origin, in three foreign ponds at the same elevation, and in four ponds at different elevation (1500-2000 m higher or lower). There were two source populations from each elevation, and adults were held in a common environment for 1 year before they were crossed to produce offspring for the experiment. Fitness was a measure that integrated larval survival, development rate, and body size. Parallel adaptation to elevation was indicated by an advantage at the home elevation (11.5% fitness difference at low elevation and 47% at high elevation). This effect was stronger than that observed in most other studies, according to a survey of previous transplant experiments across elevation (N = 8 animal species and 71 plants). Unique local adaptation within elevational zones was only 0.3-0.7 times as strong as parallel adaptation, probably because gene flow is comparatively high among nearby wetlands at the same elevation. The home-elevation advantage may reduce gene flow across the elevational gradient and enable the evolution of habitat races specialized on elevation.

Small-scale population divergence is driven by local larval environment in a temperate amphibian

Genomic variation within and among populations is shaped by the interplay between natural selection and the effects of genetic drift and gene flow. Adaptive divergence can be found in small-scale natural systems even when population sizes are small, and the potential for gene flow is high, suggesting that local environments exert selection pressures strong enough to counteract the opposing effects of drift and gene flow. Here, we investigated genomic differentiation in nine moor frog (Rana arvalis) populations in a small-scale network of local wetlands using 16,707 ddRAD-seq SNPs, relating levels of differentiation with local environments, as well as with properties of the surrounding landscape. We characterized population structure and differentiation, and partitioned the effects of geographic distance, local larval environment, and landscape features on total genomic variation. We also conducted gene-environment association studies using univariate and multivariate approaches. We found small-scale population structure corresponding to 6-8 clusters. Local larval environment was the most influential component explaining 2.3% of the total genetic variation followed by landscape features (1.8%) and geographic distance (0.8%), indicative of isolation-by-environment, -by-landscape, and -by-distance, respectively. We identified 1000 potential candidate SNPs putatively under divergent selection mediated by the local larval environment. The candidate SNPs were involved in, among other biological functions, immune system function and development. Our results suggest that small-scale environmental differences can exert selection pressures strong enough to counteract homogenizing effects of gene flow and drift in this small-scale system, leading to observable population differentiation.

Switching from mesopredator to apex predator: how do responses vary in amphibians adapted to cave living?

Abstract

The effective detection of both prey and predators is pivotal for the survival of mesopredators. However, the condition of being a mesopredator is strongly context dependent. Here we focus on two aquatic caudate species that have colonised caves: the Pyrenean newt (Calotriton asper) and the olm (Proteus anguinus). The former maintains both surface and subterranean populations, while only cave-adapted populations of the latter exist. Both species are apex predators in underground waterbodies, while the Pyrenean newt is a mesopredator in surface waterbodies. Shifting to a higher level of the trophic web through colonising caves may promote the loss of anti-predator response against surface apex predators, and an increase in the ability to detect prey. To test these two non-exclusive hypotheses, we integrated classical behavioural characterisations with a novel approach: the assessment of lateralisation (i.e. preference for one body side exposure). Behavioural experiments were performed using laboratory-reared individuals. We performed 684 trials on 39 Pyrenean newts and eight olms. Under darkness and light conditions, we tested how exposure to different chemical cues (predatory fish, prey and unknown scent) affected individuals’ activity and lateralisation. Both cave and surface Pyrenean newts responded to predator cues, while olms did not. In Pyrenean newts, predator cues reduced the time spent in movement and time spent in lateralisation associated with hunting. Our results show that predator recognition is maintained in a species where recently separated populations inhabit environments lacking of higher predators, while such behaviour tends to be lost in populations with longer history of adaptation.

Significance statement

Predator recognition can be maintained in animals adapted to predator free habitats, but varies with their history of adaptation. Species that are not at the apex of the food web can become top predators if they colonise subterranean environments. We compared the behavioural responses of the olm, a strictly cave species with a long underground evolutionary history, and of the Pyrenean newt, a facultative cave species that also has stream-dwelling populations. Moreover, we integrated a classical behavioural characterisation, such as movement detection, with a novel approach: the assessment of lateralisation. While olms do not respond to external predators scent, cave-dwelling newts still recognise it. This clearly indicates that predator recognition is still maintained in species that have colonised predator-free environments more recently.