Losing your edge: climate change and the conservation value of range-edge populations

Heterogeneity of Locked-Pasture Snow Conditions Modulate Habitat and Movement Choices of a Facultative Migrant

Habitat selection and movement are key mechanisms by which animals can respond to and potentially cope with highly variable environmental conditions. Optimal responses likely vary, however, depending on the severity and scope of conditions. We tested this hypothesis using a facultative migrant species, the Great Gray Owl (Strix nebulosa), which exhibits high inter- and intra-individual variation in the timing, direction, and distance of winter movements. Specifically, we evaluated whether episodic, spatiotemporally variable "locked-pasture" snow conditions, which restrict access to subnivean food, prompted shifts in habitat selection or long-distance movements by owls. We quantified the movement of 42 owls using global positioning system (GPS) data within the Greater Yellowstone Ecosystem, USA, during 2017-2022. We used a novel ecological application of SnowModel, a snow evolution modeling system, to estimate fine-scale, physical snow properties likely to influence access to prey. Variables included snow depth, snow crusts produced by wind, and ice crusts produced by melt-freeze and rain-on-snow events. Owls avoided heterogeneously distributed wind crusts via local shifts in habitat selection. More homogenous ice crusts elicited long-distance movements away from affected home ranges. Finally, owls employed both proximate shifts in habitat selection and long-distance movements to avoid deeper snow. Ultimately, owls exhibited behavioral flexibility in response to limiting snow conditions that can vary in terms of severity, spatial extent, and duration. Such behavioral responses determine species distribution, with implications for population and community dynamics in spatiotemporally variable systems. Understanding the effects of, and responses to, environmental controls is increasingly important given the scope of on-going global change.

Ecotypic Variation in Leaf Thermoregulation and Heat Tolerance but Not Thermal Safety Margins in Tropical Trees

To avoid reaching lethal temperatures during periods of heat stress, plants may acclimate either their biochemical thermal tolerance or leaf morphological and physiological characteristics to reduce leaf temperature (Tleaf). While plants from warmer environments may have a greater capacity to regulate Tleaf, the extent of intraspecific variation and contribution of provenance is relatively unexplored. We tested whether upland and lowland provenances of four tropical tree species grown in a common garden differed in their thermal safety margins by measuring leaf thermal traits, midday leaf-to-air temperature differences (∆Tleaf) and critical leaf temperatures defined by chlorophyll fluorescence (Tcrit). Provenance variation was species- and trait-specific. Higher ∆Tleaf and Tcrit were observed in the lowland provenance for Terminalia microcarpa, and in the upland provenance for Castanospermum australe, with no provenance effects in the other two species. Within-species covariation of Tcrit and ∆Tleaf led to a convergence of thermal safety margins across provenances. While future studies should expand the number of provenances and species investigated, our findings suggest that lowland and upland provenances may not differ substantially in their vulnerability to heat stress, as determined by thermal safety margins, despite differences in operating temperatures and Tcrit.

The importance of peripheral populations in the face of novel environmental change

Anthropogenically driven environmental change has imposed substantial threats on biodiversity, including the emergence of infectious diseases that have resulted in declines of wildlife globally. In response to pathogen invasion, maintaining diversity within host populations across heterogenous environments is essential to facilitating species persistence. White-nose syndrome is an emerging fungal pathogen that has caused mass mortalities of hibernating bats across North America. However, in the northeast, peripheral island populations of the endangered northern myotis (Myotis septentrionalis) appear to be persisting despite infection while mainland populations in the core of the species range have experienced sharp declines. Thus, this study investigated host and environmental factors that may contribute to divergent population responses. We compared patterns of pathogen exposure and infection intensity between populations and documented the environmental conditions and host activity patterns that may promote survival despite disease invasion. For island populations, we found lower prevalence and less severe infections, possibly due to a shorter hibernation duration compared to the mainland, which may reduce the time for disease progression. The coastal region of the northern myotis range may serve as habitat refugia that enables this species to persist despite pathogen exposure; however, conservation efforts could be critical to supporting species survival in the long term.

Critical reproductive behaviors in Scaled Quail and Northern Bobwhite are affected by thermal variability and mean temperature

Animals can respond differently to shifting thermal variability versus thermal averages, both of which are changing due to climate warming. How these thermal variables affect parental care behaviors can reveal the ability of parents to modify their behaviors to meet the competing demands of their offspring's thermal needs and self-maintenance, which becomes critical in suboptimal thermal conditions. Further, the time frame used to examine the interplay between temperature and behavioral shifts (e.g., seasonal patterns in care vs. drivers of individual care decisions) can provide different information about the plasticity of parental care behavior. We investigated the relationship between thermal means, thermal variability, and incubation behaviors across multiple timescales in Scaled Quail and Northern Bobwhite. Both species decreased off-bout length during periods of high thermal variability, a novel finding among studies of avian parental behavior. Further relationships between thermal endpoints (mean vs. variation) and behavior differed depending on the temporal scale. For instance, total daily time spent off the nest was not influenced by daily average temperature, yet individual off-bout duration increased with increasing average temperature in the 2 h prior to the off-bout. These results provide evidence that thermal-behavioral relationships differ across scales and likely represent a bird's ability to modify their incubation strategy to rapidly respond to the immediate thermal environment (altering individual off-bout length based on temperature) to meet self-maintenance needs while resulting in a similar outcome for their nest (total daily off-bout time). However, longer off-bout durations during high temperature events can come with reproductive costs, sometimes resulting in acute offspring mortality when eggs or chicks experience lethal temperatures.

A Review of Theory: Comparing Invasion Ecology and Climate Change-Induced Range Shifting

Human actions have led to large-scale shifts in the distributions of species, which have accelerated over recent decades. Two contributing factors include the introduction of non-native species to new regions, and more recently, the shift of species into new ranges to track a human-altered climate. While the context of these species redistributions is different, we argue here that scientists studying the effects of either invasion or range shifting are interested in the same underlying mechanistic questions: (1) What traits make a species likely to survive in a new location? and (2) Which recipient ecosystems are likely to support a newly arrived species? A survey of the theoretical literature surrounding these topics demonstrates the usefulness of this comparison and highlights key differences between range shifting and invasion in factors including genetic diversity, climatic tolerance, local adaptation, and coevolutionary history. This review does not support the blanket application of an invasion framework to climate change-induced range shifts. However, we suggest the use of modified invasion theories, experimental designs, and risk assessments could aid in predicting outcomes and prioritizing management resources for climate-threatened species.

Drivers of Intraspecific Variation in Thermal Traits and Their Importance for Resilience to Global Change in Amphibians

Intraspecific variation can be as great as variation across species, but the role of intraspecific variation in driving local and large-scale patterns is often overlooked, particularly in the field of thermal biology. In amphibians, which depend on environmental conditions and behavior to regulate body temperature, recognizing intraspecific thermal trait variation is essential to comprehensively understanding how global change impacts populations. Here, we examine the drivers of micro- and macrogeographical intraspecific thermal trait variation in amphibians. At the local scale, intraspecific variation can arise via changes in ontogeny, body size, and between the sexes, and developmental plasticity, acclimation, and maternal effects may modulate predictions of amphibian performance under future climate scenarios. At the macrogeographic scale, local adaptation in thermal traits may occur along latitudinal and elevational gradients, with seasonality and range-edge dynamics likely playing important roles in patterns that may impact future persistence. We also discuss the importance of considering disease as a factor affecting intraspecific variation in thermal traits and population resilience to climate change, given the impact of pathogens on thermal preferences and critical thermal limits of hosts. Finally, we make recommendations for future work in this area. Ultimately, our goal is to demonstrate why it is important for researchers to consider intraspecific variation to determine the resilience of amphibians to global change.

The consequences of heatwaves for the reproductive success and physiology of the wingless sub-Antarctic fly Anatalanta aptera

Sub-lethal effects of warming temperatures are an important, yet sometimes overlooked impact of climate change that may threaten the long-term survival of numerous species. This, like many other effects of climate change, is especially concerning for cold-adapted ectotherms living in rapidly warming polar regions. This study examines the effects of warmer temperatures on cold-adapted Diptera, using the long-lived sub-Antarctic sphaerocerid fly, Anatalanta aptera, as a focal species. We conducted two experiments to assess heat stress in adult flies, one varying the intensity of the heat stress (daily heating from 4 °C to 8 °C, 20 °C, or 24 °C) and one varying the frequency of heat stress exposure (heating from 4 °C to 12 °C every one, two, or three days) and examined consequences for reproductive success and metabolic responses. We found that more heat stress reduced reproductive output, but not timing of reproduction. Surprisingly, individuals sampled at different times during heat stress exposure were undifferentiable when all metabolite concentrations were analysed with redundancy analysis, however some individual metabolites did exhibit significant differences. Overall, our findings suggest that warmer temperatures in the sub-Antarctic may put this species at greater risk, especially when combined with other concurrent threats from biological invasions.

"A comparison of thermal stress response between Drosophila melanogaster and Drosophila pseudoobscura reveals differences between species and sexes"

The environment is changing faster than anticipated due to climate change, making species more vulnerable to its impacts. The level of vulnerability of species is influenced by factors such as the degree and duration of exposure, as well as the physiological sensitivity of organisms to changes in their environments, which has been shown to vary among species, populations, and individuals. Here, we compared physiological changes in fecundity, critical thermalmaximum (CTmax), respiratory quotient (RQ), and DNA damage in ovaries in response to temperature stress in two species of fruit fly, Drosophila melanogaster (25 vs. 29.5 °C) and Drosophila pseudoobscura (20.5 vs. 25 °C). The fecundity of D. melanogaster was more affected by high temperatures when exposed during egg through adult development, while D. pseudoobscura was most significantly affected when exposed to high temperatures exclusively during egg through pupal development. Additionally, D. melanogaster males exhibited a decrease of CTmax under high temperatures, while females showed an increase of CTmax when exposed to high temperatures during egg through adult development. while D. pseudoobscura females and males showed an increased CTmax only when reared at high temperatures during egg through pupae development. Moreover, both species showed an acceleration in oogenesis and an increase in apoptosis due to heat stress. These changes can likely be attributed to key differences in the geographic range, thermal range, development time, and other different factors between these two systems. Through this comparison of variation in physiology and developmental response to thermal stress, we found important differences between species and sexes that suggest future work needs to account for these factors separately in understanding the effects of constant increased temperatures.

The role of genotypic and climatic variation at the range edge: A case study in winegrapes

PREMISE

Changes in habitat suitability due to climate change are causing range shifts, with new habitat potentially available at cold range edges. We must predict these range shifts, but forecasters have limited knowledge of how genetic differences in plant physiological tolerances influence range shifts. Here, we focus on a major determinant of species ranges-physiological tolerance to extreme cold-to ask how warming over recent decades and genetic variation shape expansion across complex landscapes.

METHODS

We examined how genotypes vary in maximum cold tolerance from 9 years of cold hardiness data across 18 genotypes from 13 sites, using winegrapes (Vitis vinifera subsp. vinifera) as a case study. Combining a Bayesian hierarchical dose-response model with gridded climate data, we then project changes in climatic suitability near winegrapes' current cold range-edge between 1949 and 2016.

RESULTS

Plants increased maximum cold hardiness non-linearly with decreasing air temperature (maximum cold hardiness: -23.6°C), but with substantial (by 2°C) variation across genotypes. Our results suggest, since the 1980s, decreasing freeze injury risk has made conditions more favorable for all genotypes at the cold range-edge, but conditions remained more favorable for more cold hardy genotypes and in warmer areas. There was substantial spatial variation in habitat suitability, with the majority of suitably warm habitat located in a narrow north-south oriented strip.

CONCLUSIONS

We highlight the importance of genotypic differences in physiological tolerances when assessing range shift potential with climate change. Habitat improvements were unevenly distributed over the spatially complex landscape, though, emphasizing the importance of dispersal in range expansion.

Climate change and carnivores: shifts in the distribution and effectiveness of protected areas in the Amazon

Background

Carnivore mammals are animals vulnerable to human interference, such as climate change and deforestation. Their distribution and persistence are affected by such impacts, mainly in tropical regions such as the Amazon. Due to the importance of carnivores in the maintenance and functioning of the ecosystem, they are extremely important animals for conservation. We evaluated the impact of climate change on the geographic distribution of carnivores in the Amazon using Species Distribution Models (SDMs). Do we seek to answer the following questions: (1) What is the effect of climate change on the distribution of carnivores in the Amazon? (2) Will carnivore species lose or gain representation within the Protected Areas (PAs) of the Amazon in the future?

Methods

We evaluated the distribution area of 16 species of carnivores mammals in the Amazon, based on two future climate scenarios (RCP 4.5 and RCP 8.5) for the year 2070. For the construction of the SDMs we used bioclimatic and vegetation cover variables (land type). Based on these models, we calculated the area loss and climate suitability of the species, as well as the effectiveness of the protected areas inserted in the Amazon. We estimated the effectiveness of PAs on the individual persistence of carnivores in the future, for this, we used the SDMs to perform the gap analysis. Finally, we analyze the effectiveness of PAs in protecting taxonomic richness in future scenarios.

Results

The SDMs showed satisfactory predictive performance, with Jaccard values above 0.85 and AUC above 0.91 for all species. In the present and for the future climate scenarios, we observe a reduction of potencial distribution in both future scenarios (RCP4.5 and RCP8.5), where five species will be negatively affected by climate change in the RCP 4.5 future scenario and eight in the RCP 8.5 scenario. The remaining species stay stable in terms of total area. All species in the study showed a loss of climatic suitability. Some species lost almost all climatic suitability in the RCP 8.5 scenario. According to the GAP analysis, all species are protected within the PAs both in the current scenario and in both future climate scenarios. From the null models, we found that in all climate scenarios, the PAs are not efficient in protecting species richness.

Combined threats of climate change and contaminant exposure through the lens of bioenergetics

Organisms face energetic challenges of climate change in combination with suites of natural and anthropogenic stressors. In particular, chemical contaminant exposure has neurotoxic, endocrine-disrupting, and behavioral effects which may additively or interactively combine with challenges associated with climate change. We used a literature review across animal taxa and contaminant classes, but focused on Arctic endotherms and contaminants important in Arctic ecosystems, to demonstrate potential for interactive effects across five bioenergetic domains: (1) energy supply, (2) energy demand, (3) energy storage, (4) energy allocation tradeoffs, and (5) energy management strategies; and involving four climate change-sensitive environmental stressors: changes in resource availability, temperature, predation risk, and parasitism. Identified examples included relatively equal numbers of synergistic and antagonistic interactions. Synergies are often suggested to be particularly problematic, since they magnify biological effects. However, we emphasize that antagonistic effects on bioenergetic traits can be equally problematic, since they can reflect dampening of beneficial responses and result in negative synergistic effects on fitness. Our review also highlights that empirical demonstrations remain limited, especially in endotherms. Elucidating the nature of climate change-by-contaminant interactive effects on bioenergetic traits will build toward determining overall outcomes for energy balance and fitness. Progressing to determine critical species, life stages, and target areas in which transformative effects arise will aid in forecasting broad-scale bioenergetic outcomes under global change scenarios.

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

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

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

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

Jack pine of all trades: Deciphering intraspecific variability of a key adaptive trait at the rear edge of a widespread fire-embracing North American conifer

PREMISE

Understanding mechanisms fostering long-term persistence of marginal populations should provide key insights about species resilience facing climate change. Cone serotiny is a key adaptive trait in Pinus banksiana (jack pine), which shows phenotypic variation according to the fire regime. Compared to range-core populations within the fire-prone boreal forest, low and variable serotiny in rear-edge populations suggest local adaptation to uncommon and unpredictable wildfire regime. We assessed environmental/physiological factors that might modulate intraspecific variation in cone serotiny.

METHODS

We experimentally subjected closed cones to incrementing temperatures, then tested seed germination to determine whether and how various ecological factors (cone age, branch height, tree size, tree age) are related to cone dehiscence and seed viability in jack pines from rear-edge and range-core populations in eastern Canada.

RESULTS

Cones from rear-edge populations dehisce at a lower opening temperature, which increases with cone age. Cones from range-core stands open at a more constant, yet higher temperature. Cones from rear-edge stands take between 13 and 27 years to reach the level of serotiny achieved at the range core. At the rear edge, seed viability is steady (51%), whereas it decreases from 70% to 30% in 20 years at the range core.

CONCLUSIONS

We inferred the mechanisms of a bet-hedging strategy in rear-edge populations, which ensures steady recruitment during fire-free intervals and successful postfire regeneration. This capacity to cope with infrequent and unpredictable fire regime should increase the resilience of jack pine populations as global changes alter fire dynamics of the boreal forest.

Exploring Old Data with New Tricks: Long-Term Monitoring Indicates Spatial and Temporal Changes in Populations of Sympatric Prairie Grouse in the Nebraska Sandhills

The contiguous grasslands of the Sandhills region in Nebraska, USA, provide habitat for two sympatric, grassland-obligate species of grouse, the greater prairie-chicken (Tympanuchus cupido pinnatus) and the plains sharp-tailed grouse (Tympanuchus phasianellus jamesi). Collectively referred to as prairie grouse, these birds are monitored and managed jointly by wildlife practitioners who face the novel challenge of conserving historically allopatric species in shared range. We reconstructed region-wide and route-specific prairie grouse population trends in the Sandhills, using a 63-year timeseries of breeding ground counts aggregated from old reports and paper archives. Our objective was to repurpose historical data collected for harvest management to address questions pertinent to the conservation of prairie grouse, species whose populations have declined precipitously throughout their respective ranges. Because we cannot change the sampling protocol of historical data to answer new questions, we applied 3 different methods of data analysis—traditional regional mean counts used to adjust harvest regulations, spatially implicit, site-specific counts, and spatially explicit trends. Prairie-chicken populations have increased since the 1950s, whereas sharp-tailed grouse populations have remained stable or slightly declined. However, each species exhibited unique shifts in abundance and distribution over time, and regional indices masked important aspects of population change. Our findings indicate that legacy data have the capacity to tell new stories apart from the questions they were collected to answer. By integrating concepts from landscape ecology—a discipline that emerged decades after the collection of our count data began—we demonstrate the potential of historical data to address questions of modern-day conservation concern, using prairie grouse as a case study.

Once upon a time in Mexico: Holocene biogeography of the spotted bat (Euderma maculatum)

Holocene-era range expansions are relevant to understanding how a species might respond to the warming and drying climates of today. The harsh conditions of North American deserts have phylogenetically structured desert bat communities but differences in flight capabilities are expected to affect their ability to compete, locate, and use habitat in the face of modern climate change. A highly vagile but data-deficient bat species, the spotted bat (Euderma maculatum), is thought to have expanded its range from central Mexico to western Canada during the Holocene. With specimens spanning this latitudinal extent, we examined historical demography, and used ecological niche modeling (ENM) and phylogeography (mitochondrial DNA), to investigate historic biogeography from the rear to leading edges of the species' range. The ENM supported the notion that Mexico was largely the Pleistocene-era range, whereas haplotype pattern and Skyline plots indicated that populations expanded from the southwestern US throughout the Holocene. This era provided substantial gains in suitable climate space and likely facilitated access to roosting habitat throughout the US Intermountain West. Incongruent phylogenies among different methods prevented a precise understanding of colonization history. However, isolation at the southern-most margin of the range suggests a population was left behind in Mexico as climate space contracted and are currently of unknown status. The species appears historically suited to follow shifts in climate space but differences in flight behaviors between leading edge and core-range haplogroups suggest range expansions could be influenced by differences in habitat quality or climate (e.g., drought). Although its vagility could facilitate response to environmental change and thereby avoid extinction, anthropogenic pressures at the core range could still threaten the ability for beneficial alleles to expand into the leading edge.

Past climatic refugia and landscape resistance explain spatial genetic structure in Oriental beech in the South Caucasus

Predicting species-level effects of climatic changes requires unraveling the factors affecting the spatial genetic composition. However, disentangling the relative contribution of historical and contemporary drivers is challenging. By applying landscape genetics and species distribution modeling, we investigated processes that shaped the neutral genetic structure of Oriental beech (Fagus orientalis), aiming to assess the potential risks involved due to possible future distribution changes in the species. Using nuclear microsatellites, we analyze 32 natural populations from the Georgia and Azerbaijan (South Caucasus). We found that the species colonization history is the most important driver of the genetic pattern. The detected west-east gradient of genetic differentiation corresponds strictly to the Colchis and Hyrcanian glacial refugia. A significant signal of associations to environmental variables suggests that the distinct genetic composition of the Azerbaijan and Hyrcanian stands might also be structured by the local climate. Oriental beech retains an overall high diversity; however, in the context of projected habitat loss, its genetic resources might be greatly impoverished. The most affected are the Azerbaijan and Hyrcanian populations, for which the detected genetic impoverishment may enhance their vulnerability to environmental change. Given the adaptive potential of range-edge populations, the loss of these populations may ultimately affect the specie's adaptation, and thus the stability and resilience of forest ecosystems in the Caucasus ecoregion. Our study is the first approximation of the potential risks involved, inducing far-reaching conclusions about the need of maintaining the genetic resources of Oriental beech for a species' capacity to cope with environmental change.

Evolution of thermal tolerance and phenotypic plasticity under rapid and slow temperature fluctuations

Global warming is associated with an increase in sea surface temperature and its variability. The consequences of evolving in variable, fluctuating environments are explored by a large body of theory: when populations evolve in fluctuating environments the frequency of fluctuations determines the shapes of tolerance curves (indicative of habitats that organisms can inhabit) and trait reaction norms (the phenotypes that organisms display across these environments). Despite this well-established theoretical backbone, predicting how trait and tolerance curves will evolve in organisms at the foundation of marine ecosystems remains a challenge. Here, we used a globally distributed phytoplankton, Thalassiosira pseudonana, and show that fluctuations in temperature on scales of 3–4 generations rapidly selected for populations with enhanced trait plasticity and elevated thermal tolerance. Fluctuations spanning 30–40 generations selected for the formation of two stable, genetically and physiologically distinct populations, one evolving high trait plasticity and enhanced thermal tolerance, and the other, akin to samples evolved under constant warming, with lower trait plasticity and a smaller increase in thermal tolerance.

Genetic Structure of Natural Northern Range-Margin Mainland, Peninsular, and Island Populations of Northern Red Oak (Quercus rubra L.)

Plant populations at the leading edge of the species’ native range often exhibit genetic structure as a result of genetic drift and adaptation to harsh environmental conditions. Hence, they are likely to harbour rare genetic adaptations to local environmental conditions and therefore are of particular interest to understand climate adaptation. We examined genetic structure of nine northern marginal mainland, peninsular and isolated island natural populations of northern red oak (Quercus rubraL.), a valuable long-lived North American hardwood tree species, covering a wide climatic range, using 17 nuclear microsatellites. We found pronounced genetic differentiation of a disjunct isolated island population from all mainland and peninsular populations. Furthermore, we observed remarkably strong fine-scale spatial genetic structure (SGS) in all investigated populations. Such high SGS values are uncommon and were previously solely observed in extreme range-edge marginal oak populations in one other study. We found a significant correlation between major climate parameters and SGS formation in northern range-edge red oak populations, with more pronounced SGS in colder and drier regions. Most likely, the harsh environment in leading edge populations influences the density of reproducing trees within the populations and therefore leads to restricted overlapping of seed shadows when compared to more central populations. Accordingly, SGS was negatively correlated with effective population size and increased with latitude of the population locations. The significant positive association between genetic distances and precipitation differences between populations may be indicative of isolation by adaptation in the observed range-edge populations. However, this association was not confirmed by a multiple regression analysis including geographic distances and precipitation distances, simultaneously. Our study provides new insights in the genetic structure of long-lived tree species at their leading distribution edge.

Do marginal plant populations enhance the fitness of larger core units under ongoing climate change? Empirical insights from a rare carnation

Assisted gene flow (AGF) can restore fitness in small plant populations. Due to climate change, current fitness patterns could vary in the future ecological scenario, as highly performant lineages can undergo maladaptation under the new climatic contexts. Peripheral populations have been argued to represent a potential source of species adaptation against climate change, but experimental evidence is poor. This paper considers the consequences of within- and between-population mating between a large core population and the southernmost population, the rare Dianthus guliae, to evaluate optimal AGF design under current and future conditions. We performed experimental self-pollinations and within- and between-population cross-pollinations to generate seed material and test its adaptive value to aridity. Seed germination, seedling growth and survival were measured under current and expected aridity. Effects of population type, pollination treatment and stress treatment on fitness components were analysed by generalized linear models. Relative measures of inbreeding depression and heterosis were taken under different stress treatments. Self-pollination reduced fitness for all the considered traits compared to within- and between-population cross-pollination. Under current aridity regime, the core population expressed higher fitness, and a larger magnitude of inbreeding depression. This indicated the core unit is close to its fitness optimum and could allow for restoring the fitness of the small peripheral population. Contrarily, under increased aridity, the fitness of outbred core lineages decreased, suggesting the rise of maladaptation. In this scenario, AGF from the small peripheral population enhanced the fitness of the core unit, whereas AGF from the core population promoted a fitness loss in the peripheral population. Hence, the small peripheral population could improve fitness of large core units versus climate change, while the contrary could be not true. Integrating reciprocal breeding programmes and fitness analyses under current and predicted ecological conditions can support optimal AGF design in a long-term perspective.

Genetic purity of a rear-edge population of Carex podogyna Franch. et Sav. (Cyperaceae) maintained under interspecific hybridization

Interspecific hybridization is a critical issue in conservation biology because it may drive small populations to extinction through direct or indirect processes. In this study, to develop a conservation strategy for an endangered rear-edge population of Carex podogyna in Ashiu, Kyoto, Japan, we performed a molecular genetic analysis of the wild population and an ex-situ population established from wild seeds. Microsatellite genotypic data revealed a complete loss of genetic diversity in the wild population, suggesting that it has long been prone to genetic drift due to isolation as a small population. In contrast, microsatellite analysis of 13 ex-situ individuals detected multiple alleles that are not harbored in the wild C. podogyna population. Sequence analysis revealed that these individuals are likely natural hybrids between C. podogyna and a co-occurring species, C. curvicollis, although established hybrids have never been found in the natural habitat. Based on our observation of variegated leaves in hybrid individuals, we propose that hybrids have been excluded by natural selection and/or interspecific competition caused by insufficient productivity of photosynthesis, although other genetic and ecological factors may also be influential. Overall, this study indicates that natural mechanisms selectively removing the hybrids have maintained the genetic purity of this rear-edge population of C. podogyna, and also emphasizes the importance of genetic assessment in ex-situ conservation programs.

Adult cold tolerance and potential North American distribution of Myllocerus undecimpustulatus undatus (Coleoptera: Curculionidae)

Cold tolerance and potential distribution of Myllocerus undecimpustulatus undatus Marshall, a polyphagous pest in the United States, were investigated. Adult survivorship after 2 days at 0 °C and − 5 °C averaged 60% and 18%, respectively. Four days of exposure resulted in survivorship of 11% at 0 °C and 4% at − 5 °C, respectively. Summer-collected weevils at − 5 °C through repeated cold exposure of 2 h survived 3 times longer than those subjected to sustained cold period of 10 h. Leaf consumption did not differ among summer-collected weevils at constant 20 °C and repeated cold exposure treatments; weevils under sustained cold exposure consumed less than weevils in repeated cold exposure treatments. Leaf area consumed after cold exposure was 2–4 times greater in winter-collected weevils compared to summer-collected weevils. Leaf consumption by winter-collected weevils decreased as the number of repeated cold exposure periods increased. Locality data from collections in Florida during 2000–2012 were used to produce a correlative model complemented by a mechanistic model from the cold tolerance data to project the potential distribution of M. undecimpustulatus undatus in North America. The models support the hypothesis that M. undecimpustulatus undatus could spread to areas of the southeastern and western United States. The predicted northern distribution followed an isothermal line about 33° North. The niche model defined an area along the western Gulf Coast as unsuitable for the weevil, possibly because the area receives greater annual rainfall than other areas of the southeastern United States and has aquic or udic soil unlike the well-drained sandy soil of peninsular Florida.

Seasonal fluctuations in symbiotic bacteria and their role in environmental adaptation of the scleractinian coral Acropora pruinosa in high-latitude coral reef area of the South China Sea

Coral-associated bacterial communities are paramount for coral ecosystems and holobiont health. However, the role of symbiotic bacteria in the adaptation of high-latitude corals to seasonal fluctuations remains underexplored. Therefore, we used 16S rRNA-based high-throughput sequencing to analyze the symbiotic bacterial diversity, composition, and core bacterial community in high-latitude coral and explored the seasonal fluctuation characteristics of symbiotic bacterial communities. We found that bacterial richness and α-diversity changed significantly across different seasons. Additionally, the community structure recombined seasonally, with different dominant bacterial phyla and genera in different seasons. However, the symbiotic bacterial community structures of Acropora pruinosa in winter and spring were similar. Proteobacteria were the dominant bacteria in spring, autumn, and winter. In summer, the dominant bacterial taxa were Bacteroidota and Proteobacteria. Ralstonia was the dominant bacterial genus in spring and winter, whereas in autumn, BD1-7_clade was dominant. Linear discriminant analysis effect size identified 20 abundant genera between the different groups. Core microbiome analysis revealed that 12 core bacterial operational taxonomic units were associated with A. pruinosa in all seasons, seven of which varied with the seasons, changing between dominant and rare. Distance-based redundancy and variation partitioning analyses revealed that sea surface temperature was the major contributor of variation in the microbial community structure. We hypothesized that the high diversity and abundance of symbiotic bacteria and the increase in Prosthecochloris abundance in coral in summer can help A. pruinosa maintain its physiological functions, ameliorating the negative physiological effects of the decrease in Symbiodiniaceae density under high-temperature stress. Thus, the rapid reorganization of the symbiotic bacterial community structure and core microflora in different seasons may allow the corals to adapt to large seasonal environmental fluctuations. In conclusion, seasonal variation of bacteria plays an important role in coral adaptation to large environmental fluctuations.

Factors influencing distributional shifts and abundance at the range core of a climate-sensitive mammal

Species are frequently responding to contemporary climate change by shifting to higher elevations and poleward to track suitable climate space. However, depending on local conditions and species' sensitivity, the nature of these shifts can be highly variable and difficult to predict. Here, we examine how the American pika (Ochotona princeps), a philopatric, montane lagomorph, responds to climatic gradients at three spatial scales. Using mixed-effects modeling in an information-theoretic approach, we evaluated a priori model suites regarding predictors of site occupancy, relative abundance, and elevational-range retraction across 760 talus patches, nested within 64 watersheds across the Northern Rocky Mountains of North America, during 2017-2020. The top environmental predictors differed across these response metrics. Warmer temperatures in summer and winter were associated with lower occupancy, lower relative abundances, and greater elevational retraction across watersheds. Occupancy was also strongly influenced by habitat patch size, but only when combined with climate metrics such as actual evapotranspiration. Using a second analytical approach, acute heat stress and summer precipitation best explained retraction residuals (i.e., the relative extent of retraction given the original elevational range of occupancy). Despite the study domain occurring near the species' geographic-range center, where populations might have higher abundances and be at lower risk of climate-related stress, 33.9% of patches showed evidence of recent extirpations. Pika-extirpated sites averaged 1.44℃ warmer in summer than did occupied sites. Additionally, the minimum elevation of pika occupancy has retracted upslope in 69% of watersheds (mean: 281 m). Our results emphasize the nuance associated with evaluating species' range dynamics in response to climate gradients, variability, and temperature exceedances, especially in regions where species occupy gradients of conditions that may constitute multiple range edges. Furthermore, this study highlights the importance of evaluating diverse drivers across response metrics to improve the predictive accuracy of widely used, correlative models.

The evolutionary genomics of species' responses to climate change

Climate change is a threat to biodiversity. One way that this threat manifests is through pronounced shifts in the geographical range of species over time. To predict these shifts, researchers have primarily used species distribution models. However, these models are based on assumptions of niche conservatism and do not consider evolutionary processes, potentially limiting their accuracy and value. To incorporate evolution into the prediction of species' responses to climate change, researchers have turned to landscape genomic data and examined information about local genetic adaptation using climate models. Although this is an important advancement, this approach currently does not include other evolutionary processes-such as gene flow, population dispersal and genomic load-that are critical for predicting the fate of species across the landscape. Here, we briefly review the current practices for the use of species distribution models and for incorporating local adaptation. We next discuss the rationale and theory for considering additional processes, reviewing how they can be incorporated into studies of species' responses to climate change. We summarize with a conceptual framework of how manifold layers of information can be combined to predict the potential response of specific populations to climate change. We illustrate all of the topics using an exemplar dataset and provide the source code as potential tutorials. This Perspective is intended to be a step towards a more comprehensive integration of population genomics with climate change science.

Demographic modeling informs functional connectivity and management interventions in Graham’s beardtongue

Functional connectivity (i.e., the movement of individuals across a landscape) is essential for the maintenance of genetic variation and persistence of rare species. However, illuminating the processes influencing functional connectivity and ultimately translating this knowledge into management practice remains a fundamental challenge. Here, we combine various population structure analyses with pairwise, population-specific demographic modeling to investigate historical functional connectivity in Graham’s beardtongue (Penstemon grahamii), a rare plant narrowly distributed across a dryland region of the western US. While principal component and population structure analyses indicated an isolation-by-distance pattern of differentiation across the species’ range, spatial inferences of effective migration exposed an abrupt shift in population ancestry near the range center. To understand these seemingly conflicting patterns, we tested various models of historical gene flow and found evidence for recent admixture (~ 3400 generations ago) between populations near the range center. This historical perspective reconciles population structure patterns and suggests management efforts should focus on maintaining connectivity between these previously isolated lineages to promote the ongoing transfer of genetic variation. Beyond providing species-specific knowledge to inform management options, our study highlights how understanding demographic history may be critical to guide conservation efforts when interpreting population genetic patterns and inferring functional connectivity.

Climate warming threatens the persistence of a community of disturbance-adapted native annual plants

With ongoing climate change, populations are expected to exhibit shifts in demographic performance that will alter where a species can persist. This presents unique challenges for managing plant populations and may require ongoing interventions, including in situ management or introduction into new locations. However, few studies have examined how climate change may affect plant demographic performance for a suite of species, or how effective management actions could be in mitigating climate change effects. Over the course of two experiments spanning 6 yr and four sites across a latitudinal gradient in the Pacific Northwest, United States, we manipulated temperature, precipitation, and disturbance intensity, and quantified effects on the demography of eight native annual prairie species. Each year we planted seeds and monitored germination, survival, and reproduction. We found that disturbance strongly influenced demographic performance and that seven of the eight species had increasingly poor performance with warmer conditions. Across species and sites, we observed 11% recruitment (the proportion of seeds planted that survived to reproduction) following high disturbance, but just 3.9% and 2.3% under intermediate and low disturbance, respectively. Moreover, mean seed production following high disturbance was often more than tenfold greater than under intermediate and low disturbance. Importantly, most species exhibited precipitous declines in their population growth rates (λ) under warmer-than-ambient experimental conditions and may require more frequent disturbance intervention to sustain populations. Aristida oligantha, a C4 grass, was the only species to have λ increase with warmer conditions. These results suggest that rising temperatures may cause many native annual plant species to decline, highlighting the urgency for adaptive management practices that facilitate their restoration or introduction to newly suitable locations. Frequent and intense disturbances are critical to reduce competitors and promote native annuals' persistence, but even such efforts may prove futile under future climate regimes.

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.

Contrasting multitaxon responses to climate change in Mediterranean mountains

We explored the influence of climatic factors on diversity patterns of multiple taxa (lichens, bryophytes, and vascular plants) along a steep elevational gradient to predict communities' dynamics under future climate change scenarios in Mediterranean regions. We analysed (1) species richness patterns in terms of heat-adapted, intermediate, and cold-adapted species; (2) pairwise beta-diversity patterns, also accounting for its two different components, species replacement and richness difference; (3) the influence of climatic variables on species functional traits. Species richness is influenced by different factors between three taxonomic groups, while beta diversity differs mainly between plants and cryptogams. Functional traits are influenced by different factors in each taxonomic group. On the basis of our observations, poikilohydric cryptogams could be more impacted by climate change than vascular plants. However, contrasting species-climate and traits-climate relationships were also found between lichens and bryophytes suggesting that each group may be sensitive to different components of climate change. Our study supports the usefulness of a multi-taxon approach coupled with a species traits analysis to better unravel the response of terrestrial communities to climate change. This would be especially relevant for lichens and bryophytes, whose response to climate change is still poorly explored.

The influence of species life history and distribution characteristics on species responses to habitat fragmentation in an urban landscape

Fragmentation within urbanized environments often leads to a loss of native species diversity; however, variation exists in responses among-species and among-populations within species. We aimed to identify patterns in species biogeography in an urbanized landscape to understand anthropogenic effects on vertebrate communities and identify species that are more sensitive or resilient to landscape change. We investigated patterns in species richness and species responses to fragmentation in southern Californian small vertebrate communities using multispecies occupancy models and determined factors associated with overall commonness and sensitivity to patch size for 45 small vertebrate species both among and within remaining non-developed patches. In general, smaller patches had fewer species, with amphibian species richness being particularly sensitive to patch size effects. Mammals were generally more common, occurring both in a greater proportion of patches and a higher proportion of the sites within occupied patches. Alternatively, amphibians were generally restricted to larger patches but were more ubiquitous within smaller patches when occupied. Species range size was positively correlated with how common a species was across and within patches, even when controlling for only patches that fell within a species' range. We found sensitivity to patch size was greater for more fecund species and depended on where the patch occurred within a species' range. While all taxa were more likely to occur in patches in the warmer portions of their ranges, amphibians and mammals were more sensitive to fragmentation in these warmer areas as compared to the rest of their ranges. Similarly, amphibians occurred at a smaller proportion of sites within patches in drier portions of their ranges. Mammals occurred at a higher proportion of sites that were also in drier portions of their range while reptiles did not differ in their sensitivity to patch size by range position. We demonstrate that taxonomy, life history, range size and range position can predict commonness and sensitivity of species across this highly fragmented yet biodiverse landscape. The impacts of fragmentation on species communities within an urban landscape depend on scale, with differences emerging among and within species and populations.

Ecological Diversity within Rear-Edge: A Case Study from Mediterranean Quercus pyrenaica Willd

Understanding the ecology of populations located in the rear edge of their distribution is key to assessing the response of the species to changing environmental conditions. Here, we focus on rear-edge populations of Quercus pyrenaica in Sierra Nevada (southern Iberian Peninsula) to analyze their ecological and floristic diversity. We perform multivariate analyses using high-resolution environmental information and forest inventories to determine how environmental variables differ among oak populations, and to identify population groups based on environmental and floristic composition. We find that water availability is a key variable in explaining the distribution of Q. pyrenaica and the floristic diversity of their accompanying communities within its rear edge. Three cluster of oak populations were identified based on environmental variables. We found differences among these clusters regarding plant diversity, but not for forest attributes. A remarkable match between the populations clustering derived from analysis of environmental variables and the ordination of the populations according to species composition was found. The diversity of ecological behaviors for Q. pyrenaica populations in this rear edge are consistent with the high genetic diversity shown by populations of this oak in the Sierra Nevada. The identification of differences between oak populations within the rear-edge with respect to environmental variables can aid with planning the forest management and restoration actions, particularly considering the importance of some environmental factors in key ecological aspects.

Hurricanes overcome migration lag and shape intraspecific genetic variation beyond a poleward mangrove range limit

Expansion of many tree species lags behind climate change projections. Extreme storms can rapidly overcome this lag, especially for coastal species, but how will storm-driven expansion shape intraspecific genetic variation? Do storms provide recruits only from the nearest sources, or from more distant sources? Answers to these questions have ecological and evolutionary implications, but empirical evidence is absent from the literature. In 2017, Hurricane Irma provided an opportunity to address this knowledge gap at the northern range limit of the neotropical black mangrove (Avicennia germinans) on the Atlantic coast of Florida, USA. We observed massive post-hurricane increases in beach-stranded A. germinans propagules at, and past, this species' present day range margin when compared to a previously surveyed nonhurricane year. Yet, propagule dispersal does not guarantee subsequent establishment and reproductive success (i.e., effective dispersal). We also evaluated prior effective dispersal along this coastline with isolated A. germinans trees identified beyond the most northern established population. We used 12 nuclear microsatellite loci to genotype 896 hurricane-driven drift propagules from nine sites and 10 isolated trees from four sites, determined their sources of origin, and estimated dispersal distances. Almost all drift propagules and all isolated trees came from the nearest sources. This research suggests that hurricanes are a prerequisite for poleward range expansion of a coastal tree species and that storms can shape the expanding gene pool by providing almost exclusively range-margin genotypes. These insights and empirical estimates of hurricane-driven dispersal distances should improve our ability to forecast distributional shifts of coastal species.

Geographic variation in adult and embryonic desiccation tolerance in a terrestrial-breeding frog

Intraspecific variation in the ability of individuals to tolerate environmental perturbations is often neglected when considering the impacts of climate change. Yet this information is potentially crucial for mitigating deleterious effects of climate change on threatened species. Here we assessed patterns of intraspecific variation in desiccation tolerance in the frog Pseudophryne guentheri, a terrestrial-breeding species experiencing a drying climate. Adult frogs were collected from six populations across a rainfall gradient and their dehydration and rehydration rates were assessed. We also compared desiccation tolerance of embryos and hatchlings originating from within-population parental crosses from four of the populations. Embryos were reared on soil at three soil-water potentials and their desiccation tolerance was assessed across a range of traits. We found significant and strong patterns of intraspecific variation in almost all traits, both in adults and first-generation offspring. Adult frogs exhibited clinal variation in their water balance responses, with populations from drier sites both dehydrating and rehydrating more slowly compared to frogs from more mesic sites. Similarly, desiccation tolerance of first-generation offspring was significantly greater in populations from xeric sites. Our findings suggest that populations within this species will respond differently to the regional reduction in rainfall predicted by climate change models.

Integrating intraseasonal grassland dynamics in cross-scale distribution modeling to support waterbird recovery plans

Despite much discussion about the utility of remote sensing for effective conservation, the inclusion of these technologies in species recovery plans remains largely anecdotal. We developed a modeling approach for the integration of local, spatially measured ecosystem functional dynamics into a species distribution modeling (SDM) framework in which other ecologically relevant factors are modeled separately at broad scales. To illustrate the approach, we incorporated intraseasonal water-vegetation dynamics into a cross-scale SDM for the Common Snipe (Gallinago gallinago), which is highly dependent on water and vegetation dynamics. The Common Snipe is an Iberian grassland waterbird characteristic of European agricultural meadows and a member of one of the most threatened bird guilds. The intraseasonal dynamics of water content of vegetation were measured using the standard deviation of the normalized difference water index time series computed from bimonthly images of the Sentinel-2 satellite. The recovery plan for the Common Snipe in Galicia (northwestern Iberian Peninsula) provided an opportunity to apply our modeling framework. Model accuracy in predicting the species' distribution at a regional scale (resulting from integration of downscaled climate projections with regional habitat-topographic suitability models) was very high (area under the curve [AUC] of 0.981 and Boyce's index of 0.971). Local water-vegetation dynamic models, based exclusively on Sentinel-2 imagery, were good predictors (AUC of 0.849 and Boyce's index of 0.976). The predictive power improved (AUC of 0.92 and Boyce's index of 0.98) when local model predictions were restricted to areas identified by the continental and regional models as priorities for conservation. Our models also performed well (AUC of 0.90 and Boyce's index of 0.93) when projected to updated water-vegetation conditions. Our modeling framework enabled incorporation of key ecosystem processes closely related to water and carbon cycles while accounting for other factors ecologically relevant to endangered grassland waterbirds across different scales, allowed identification of priority areas for conservation, and provided an opportunity for cost-effective recovery planning by monitoring management effectiveness from space.

Range edges in heterogeneous landscapes: Integrating geographic scale and climate complexity into range dynamics

The impacts of climate change have re-energized interest in understanding the role of climate in setting species geographic range edges. Despite the strong focus on species' distributions in ecology and evolution, defining a species range edge is theoretically and empirically difficult. The challenge of determining a range edge and its relationship to climate is in part driven by the nested nature of geography and the multidimensionality of climate, which together generate complex patterns of both climate and biotic distributions across landscapes. Because range-limiting processes occur in both geographic and climate space, the relationship between these two spaces plays a critical role in setting range limits. With both conceptual and empirical support, we argue that three factors-climate heterogeneity, collinearity among climate variables, and spatial scale-interact to shape the spatial structure of range edges along climate gradients, and we discuss several ways that these factors influence the stability of species range edges with a changing climate. We demonstrate that geographic and climate edges are often not concordant across species ranges. Furthermore, high climate heterogeneity and low climate collinearity across landscapes increase the spectrum of possible relationships between geographic and climatic space, suggesting that geographic range edges and climatic niche limits correspond less frequently than we may expect. More empirical explorations of how the complexity of real landscapes shapes the ecological and evolutionary processes that determine species range edges will advance the development of range limit theory and its applications to biodiversity conservation in the context of changing climate.

Prediction of present and future distribution of the Schlegel’s Japanese gecko (Gekko japonicus) using MaxEnt modeling

Background

Understanding the geographical distribution of a species is a key component of studying its ecology, evolution, and conservation. Although Schlegel’s Japanese gecko (Gekko japonicus) is widely distributed in Northeast Asia, its distribution has not been studied in detail. We predicted the present and future distribution of G. japonicus across China, Japan, and Korea based on 19 climatic and 5 environmental variables using the maximum entropy (MaxEnt) species distribution model.

Results

Present time major suitable habitats for G. japonicus, having greater than 0.55 probability of presence (threshold based on the average predicted probability of the presence records), are located at coastal and inland cities of China; western, southern, and northern coasts of Kyushu and Honshu in Japan; and southern coastal cities of Korea. Japan contained 69.3% of the suitable habitats, followed by China (27.1%) and Korea (4.2%). Temperature seasonality (66.5% of permutation importance) was the most important predictor of the distribution. Future distributions according to two climate change scenarios predicted that by 2070, and overall suitable habitats would decrease compared to the present habitats by 18.4% (scenario RCP 4.5) and 10.4% (scenario RCP 8.5). In contrast to these overall trends, range expansions are expected in inland areas of China and southern parts of Korea.

Conclusions

Suitable habitats predicted for G. japonicus are currently located in coastal cities of Japan, China, and Korea, as well as in isolated patches of inland China. Due to climate change, suitable habitats are expected to shrink along coastlines, particularly at the coastal-edge of climate change zones. Overall, our results provide essential distribution range information for future ecological studies of G. japonicus across its distribution range.