The shaping of genetic variation in edge‐of‐range populations under past and future climate change

With rates of climate change exceeding the rate at which many species are able to shift their range or adapt, it is important to understand how future changes are likely to affect biodiversity at all levels of organisation. Understanding past responses and extent of niche conservatism in climatic tolerance can help predict future consequences. We use an integrated approach to determine the genetic consequences of past and future climate changes on a bat species, Plecotus austriacus. Glacial refugia predicted by palaeo‐modelling match those identified from analyses of extant genetic diversity and model‐based inference of demographic history. Former refugial populations currently contain disproportionately high genetic diversity, but niche conservatism, shifts in suitable areas and barriers to migration mean that these hotspots of genetic diversity are under threat from future climate change. Evidence of population decline despite recent northward migration highlights the need to conserve leading‐edge populations for spearheading future range shifts.     

Contrasting climate niches among co‐occurring subdominant forbs of the sagebrush steppe

Aim

Abiotic conditions are key components that determine the distribution of species. However, co‐occurring species can respond differently to the same factors, and determining which climate components are most predictive of geographic distributions is important for understanding community response to climate change. Here, we estimate and compare climate niches of ten subdominant, herbaceous forb species common in sagebrush steppe systems, asking how niches differ among co‐occurring species and whether more closely related species exhibit higher niche overlap.

Location

Western United States.

Methods

We used herbarium records and ecological niche modelling to estimate area of occupancy, niche breadth and overlap, and describe characteristics of suitable climate. We compared mean values and variability in summer precipitation and minimum temperatures at occurrence locations among species, plant families, and growth forms, and related estimated phylogenetic distances to niche overlap.

Results

Species varied in the size and spatial distribution of suitable climate and in niche breadth. Species also differed in the variables contributing to their suitable climate and in mean values, spatial variation and interannual variation in highly predictive climate variables. Only two of ten species shared comparable climate niches. We found family‐level differences associated with variation in summer precipitation and minimum temperatures, as well as in mean minimum temperatures. Growth forms differed in their association with variability in summer precipitation and minimum temperatures. We found no relationship between phylogenetic distance and niche overlap among our species.

Main conclusions

We identified contrasting climate niches for ten Great Basin understorey forbs, including differences in both mean values and climate variability. These estimates can guide species selection for restoration by identifying species with a high tolerance for climate variability and large climatic niches. They can also help conservationists to understand which species may be least tolerant of climate variability, and potentially most vulnerable to climate change.

     

Links between plant species’ spatial and temporal responses to a warming climate

To generate realistic projections of species’ responses to climate change, we need to understand the factors that limit their ability to respond. Although climatic niche conservatism, the maintenance of a species’s climatic niche over time, is a critical assumption in niche-based species distribution models, little is known about how universal it is and how it operates. In particular, few studies have tested the role of climatic niche conservatism via phenological changes in explaining the reported wide variance in the extent of range shifts among species. Using historical records of the phenology and spatial distribution of British plants under a warming climate, we revealed that: (i) perennial species, as well as those with weaker or lagged phenological responses to temperature, experienced a greater increase in temperature during flowering (i.e. failed to maintain climatic niche via phenological changes); (ii) species that failed to maintain climatic niche via phenological changes showed greater northward range shifts; and (iii) there was a complementary relationship between the levels of climatic niche conservatism via phenological changes and range shifts. These results indicate that even species with high climatic niche conservatism might not show range shifts as instead they track warming temperatures during flowering by advancing their phenology.     

Influence of environmental heterogeneity on genetic diversity and structure in an endemic southern Californian oak

Understanding how specific environmental factors shape gene flow while disentangling their importance relative to the effects of geographical isolation is a major question in evolutionary biology and a specific goal of landscape genetics. Here, we combine information from nuclear microsatellite markers and ecological niche modelling to study the association between climate and spatial genetic structure and variability in Engelmann oak (Quercus engelmannii), a wind-pollinated species with high potential for gene flow. We first test whether genetic diversity is associated with climatic niche suitability and stability since the Last Glacial Maximum (LGM). Second, we use causal modelling to analyse the potential influence of climatic factors (current and LGM niche suitability) and altitude in the observed patterns of genetic structure. We found that genetic diversity is negatively associated with local climatic stability since the LGM, which may be due to higher immigration rates in unstable patches during favourable climatic periods and/or temporally varying selection. Analyses of spatial genetic structure revealed the presence of three main genetic clusters, a pattern that is mainly driven by two highly differentiated populations located in the northern edge of the species distribution range. After controlling for geographic distance, causal modelling analyses showed that genetic relatedness decreases with the environmental divergence among sampling sites estimated as altitude and current and LGM niche suitability. Natural selection against nonlocal genotypes and/or asynchrony in reproductive phenology may explain this pattern. Overall, this study suggests that local environmental conditions can shape patterns of genetic structure and variability even in species with high potential for gene flow and relatively small distribution ranges.     

Genetic structure and ecological niche segregation of Indian gray mongoose (Urva edwardsii) in Iran

Combining genetic data with ecological niche models is an effective approach for exploring climatic and nonclimatic environmental variables affecting spatial patterns of intraspecific genetic variation. Here, we adopted this combined approach to evaluate genetic structure and ecological niche of the Indian gray mongoose (Urva edwardsii) in Iran, as the most western part of the species range. Using mtDNA, we confirmed the presence of two highly differentiated clades. Then, we incorporated ensemble of small models (ESMs) using climatic and nonclimatic variables with genetic data to assess whether genetic differentiation among clades was coupled with their ecological niche. Climate niche divergence was also examined based on a principal component analysis on climatic factors only. The relative habitat suitability values predicted by the ESMs for both clades revealed their niche separation. Between‐clade climate only niche comparison revealed that climate space occupied by clades is similar to some extent, but the niches that they utilize differ between the distribution ranges of clades. We found that in the absence of evidence for recent genetic exchanges, distribution models suggest the species occurs in different niches and that there are apparent areas of disconnection across the species range. The estimated divergence time between the two Iranian clades (4.9 Mya) coincides with the uplifting of the Zagros Mountains during the Early Pliocene. The Zagros mountain‐building event seems to have prevented the distribution of U. edwardsii populations between the western and eastern parts of the mountains as a result of vicariance events. Our findings indicated that the two U. edwardsii genetic clades in Iran can be considered as two conservation units and can be utilized to develop habitat‐specific and climate change‐integrated management strategies.     

Intraspecific trait variation across elevation predicts a widespread tree species' climate niche and range limits

Global change is widely altering environmental conditions which makes accurately predicting species range limits across natural landscapes critical for conservation and management decisions. If climate pressures along elevation gradients influence the distribution of phenotypic and genetic variation of plant functional traits, then such trait variation may be informative of the selective mechanisms and adaptations that help define climatic niche limits. Using extensive field surveys along 16 elevation transects and a large common garden experiment, we tested whether functional trait variation could predict the climatic niche of a widespread tree species (Populus angustifolia) with a double quantile regression approach. We show that intraspecific variation in plant size, growth, and leaf morphology corresponds with the species' total climate range and certain climatic limits related to temperature and moisture extremes. Moreover, we find evidence of genetic clines and phenotypic plasticity at environmental boundaries, which we use to create geographic predictions of trait variation and maximum values due to climatic constraints across the western US. Overall, our findings show the utility of double quantile regressions for connecting species distributions and climate gradients through trait‐based mechanisms. We highlight how new approaches like ours that incorporate genetic variation in functional traits and their response to climate gradients will lead to a better understanding of plant distributions as well as identifying populations anticipated to be maladapted to future environments.     

Tests of species‐specific models reveal the importance of drought in postglacial range shifts of a Mediterranean‐climate tree: insights from integrative distributional,demographic and coalescent modelling and ABC model selection

Past climate change has caused shifts in species distributions and undoubtedly impacted patterns of genetic variation, but the biological processes mediating responses to climate change, and their genetic signatures, are often poorly understood. We test six species‐specific biologically informed hypotheses about such processes in canyon live oak (Quercus chrysolepis) from the California Floristic Province. These hypotheses encompass the potential roles of climatic niche, niche multidimensionality, physiological trade‐offs in functional traits, and local‐scale factors (microsites and local adaptation within ecoregions) in structuring genetic variation. Specifically, we use ecological niche models (ENMs) to construct temporally dynamic landscapes where the processes invoked by each hypothesis are reflected by differences in local habitat suitabilities. These landscapes are used to simulate expected patterns of genetic variation under each model and evaluate the fit of empirical data from 13 microsatellite loci genotyped in 226 individuals from across the species range. Using approximate Bayesian computation (ABC), we obtain very strong support for two statistically indistinguishable models: a trade‐off model in which growth rate and drought tolerance drive habitat suitability and genetic structure, and a model based on the climatic niche estimated from a generic ENM, in which the variables found to make the most important contribution to the ENM have strong conceptual links to drought stress. The two most probable models for explaining the patterns of genetic variation thus share a common component, highlighting the potential importance of seasonal drought in driving historical range shifts in a temperate tree from a Mediterranean climate where summer drought is common.     

Model-based uncertainty in species range prediction

Aim Many attempts to predict the potential range of species rely on environmental niche (or ‘bioclimate envelope’) modelling, yet the effects of using different niche‐based methodologies require further investigation. Here we investigate the impact that the choice of model can have on predictions, identify key reasons why model output may differ and discuss the implications that model uncertainty has for policy‐guiding applications. Location The Western Cape of South Africa. Methods We applied nine of the most widely used modelling techniques to model potential distributions under current and predicted future climate for four species (including two subspecies) of Proteaceae. Each model was built using an identical set of five input variables and distribution data for 3996 sampled sites. We compare model predictions by testing agreement between observed and simulated distributions for the present day (using the area under the receiver operating characteristic curve (AUC) and kappa statistics) and by assessing consistency in predictions of range size changes under future climate (using cluster analysis). Results Our analyses show significant differences between predictions from different models, with predicted changes in range size by 2030 differing in both magnitude and direction (e.g. from 92% loss to 322% gain). We explain differences with reference to two characteristics of the modelling techniques: data input requirements (presence/absence vs. presence‐only approaches) and assumptions made by each algorithm when extrapolating beyond the range of data used to build the model. The effects of these factors should be carefully considered when using this modelling approach to predict species ranges. Main conclusions We highlight an important source of uncertainty in assessments of the impacts of climate change on biodiversity and emphasize that model predictions should be interpreted in policy‐guiding applications along with a full appreciation of uncertainty.     

HOW IS THE RATE OF CLIMATIC‐NICHE EVOLUTION RELATED TO CLIMATIC‐NICHE BREADTH?

The rate of climatic‐niche evolution is important to many research areas in ecology, evolution, and conservation biology, including responses of species to global climate change, spread of invasive species, speciation, biogeography, and patterns of species richness. Previous studies have implied that clades with higher rates of climatic‐niche evolution among species should have species with narrower niche breadths, but there is also evidence suggesting the opposite pattern. However, the relationships between rate and breadth have not been explicitly analyzed. Here, we examine the relationships between the rate of climatic‐niche evolution and climatic‐niche breadth using phylogenetic and climatic data for 250 species in the salamander family Plethodontidae, a group showing considerable variation in both rates of climatic‐niche evolution and climatic‐niche breadths. Contrary to some expectations, we find no general relationship between climatic‐niche breadth and the rate of climatic‐niche evolution. Climatic‐niche breadths for some ecologically important climatic variables considered separately (temperature seasonality and annual precipitation) do show significant relationships with the rate of climatic‐niche evolution, but rates are faster in clades in which species have broader (not narrower) niche breadths. In summary, our results show that narrower niche breadths are not necessarily associated with faster rates of niche evolution.     

Influence of late Quaternary climate change on present patterns of genetic variation in valley oak,Quercus lobata Née

Phylogeography and ecological niche models (ENMs) suggest that late Quaternary glacial cycles have played a prominent role in shaping present population genetic structure and diversity, but have not applied quantitative methods to dissect the relative contribution of past and present climate vs. other forces. We integrate multilocus phylogeography, climate‐based ENMs and multivariate statistical approaches to infer the effects of late Quaternary climate change on contemporary genetic variation of valley oak (Quercus lobata Née). ENMs indicated that valley oak maintained a stable distribution with local migration from the last interglacial period (~120 ka) to the Last Glacial Maximum (~21 ka, LGM) to the present compared with large‐scale range shifts for an eastern North American white oak (Quercus alba L.). Coast Range and Sierra Nevada foothill populations diverged in the late Pleistocene before the LGM [104 ka (28–1622)] and have occupied somewhat distinct climate niches, according to ENMs and coalescent analyses of divergence time. In accordance with neutral expectations for stable populations, nuclear microsatellite diversity positively correlated with niche stability from the LGM to present. Most strikingly, nuclear and chloroplast microsatellite variation significantly correlated with LGM climate, even after controlling for associations with geographic location and present climate using partial redundancy analyses. Variance partitioning showed that LGM climate uniquely explains a similar proportion of genetic variance as present climate (16% vs. 11–18%), and together, past and present climate explains more than geography (19%). Climate can influence local expansion–contraction dynamics, flowering phenology and thus gene flow, and/or impose selective pressures. These results highlight the lingering effect of past climate on genetic variation in species with stable distributions.     

The Gray Toad-headed Agama,Phrynocephalus scutellatus,on the Iranian Plateau: The degree of niche overlap depends on the phylogenetic distance

The Gray Toad-headed Agama (Phrynocephalus scutellatus) occurs in Iran, Pakistan and Afghanistan and is represented in Iran by four distinctive genetic clades. We built distribution models for three of these clades (one clade was not included due to a low number of distribution records) using Maximum Entropy Algorithm in order to determine the contribution of ecological factors to the distribution pattern. The degree of spatial niche overlap between every pair of clades were measured using Schoener's D niche overlap metric. The results showed that at species-level climate variables (annual precipitation, annual mean temperature) were the most influential parameters determining the boundaries of the distribution in Iran. Temperature seasonality was found to be the most influential factor in the distribution of both Clade I and Clade II. However, this variable was replaced by the annual mean temperature for Clade VI. Based on the results of Schoener's D metric, Clades I and II had the lowest, and Clades II and VI the highest level of ecological niche overlap. Comparing the result of niche overlap with genetic distance between the clades, it was found that the ecologically least similar clades were those with the longer history of genetic segregation.     

Reconstructing the demise of Tethyan plants: climate‐driven range dynamics of Laurus since the Pliocene

Aim Climate changes are thought to be responsible for the retreat and eventual extinction of subtropical lauroid species that covered much of Europe and North Africa during the Palaeogene and early Neogene; little is known, however, of the spatial and temporal patterns of this demise. Herein we calibrate ecological niche models to assess the climatic requirements of Laurus L. (Lauraceae), an emblematic relic from the Tethyan subtropical flora, subsequently using these models to infer how the range dynamics of Laurus were affected by Plio‐Pleistocene climate changes. We also provide predictions of likely range changes resulting from future climatic scenarios. Location The Mediterranean Basin and Macaronesian islands (Canaries, Madeira, Azores). Methods We used a maximum‐entropy algorithm (Maxent) to model the relationship between climate and Laurus distribution over time. The models were fitted both to the present and to the middle Pliocene, based on fossil records. We employed climatic reconstructions for the mid‐Pliocene (3 Ma), the Last Glacial Maximum (21 ka) and a CO₂‐doubling future scenario to project putative species distribution in each period. We validated the model projections with Laurus fossil and present occurrences. Results Laurus preferentially occupied warm and moist areas with low seasonality, showing a marked stasis of its climatic niche. Models fitted to Pliocene conditions successfully predicted the current species distribution. Large suitable areas existed during the Pliocene, which were strongly reduced during the Pleistocene, but humid refugia within the Mediterranean Basin and Macaronesian islands enabled long‐term persistence. Future climate conditions are likely to re‐open areas suitable for colonization north of the current range. Main conclusions The climatic requirements of Laurus remained virtually unchanged over the last 3 Myr. This marked niche conservatism imposed largely deterministic range dynamics driven by climate conditions. This species's relatively high drought tolerance might account for the survival of Laurus in continental Europe throughout the Quaternary whilst other Lauraceae became extinct. Climatic scenarios for the end of this century would favour an expansion of the species's range towards northern latitudes, while severely limiting southern populations due to increased water stress.     

Evolution under changing climates: climatic niche stasis despite rapid evolution in a non-native plant

A topic of great current interest is the capacity of populations to adapt genetically to rapidly changing climates, for example by evolving the timing of life-history events, but this is challenging to address experimentally. I use a plant invasion as a model system to tackle this question by combining molecular markers, a common garden experiment and climatic niche modelling. This approach reveals that non-native Lactuca serriola originates primarily from Europe, a climatic subset of its native range, with low rates of admixture from Asia. It has rapidly refilled its climatic niche in the new range, associated with the evolution of flowering phenology to produce clines along climate gradients that mirror those across the native range. Consequently, some non-native plants have evolved development times and grow under climates more extreme than those found in Europe, but not among populations from the native range as a whole. This suggests that many plant populations can adapt rapidly to changed climatic conditions that are already within the climatic niche space occupied by the species elsewhere in its range, but that evolution to conditions outside of this range is more difficult. These findings can also help to explain the prevalence of niche conservatism among non-native species.     

Phenological plasticity will not help all species adapt to climate change

Concerns are rising about the capacity of species to adapt quickly enough to climate change. In long‐lived organisms such as trees, genetic adaptation is slow, and how much phenotypic plasticity can help them cope with climate change remains largely unknown. Here, we assess whether, where and when phenological plasticity is and will be adaptive in three major European tree species. We use a process‐based species distribution model, parameterized with extensive ecological data, and manipulate plasticity to suppress phenological variations due to interannual, geographical and trend climate variability, under current and projected climatic conditions. We show that phenological plasticity is not always adaptive and mostly affects fitness at the margins of the species' distribution and climatic niche. Under current climatic conditions, phenological plasticity constrains the northern range limit of oak and beech and the southern range limit of pine. Under future climatic conditions, phenological plasticity becomes strongly adaptive towards the trailing edges of beech and oak, but severely constrains the range and niche of pine. Our results call for caution when interpreting geographical variation in trait means as adaptive, and strongly point towards species distribution models explicitly taking phenotypic plasticity into account when forecasting species distribution under climate change scenarios.     

Identifying refugia from climate change using coupled ecological and genetic data in a transitional Mediterranean‐temperate tree species

Populations occurring in areas of overlap between the current and future distribution of a species are particularly important because they can represent “refugia from climate change”. We coupled ecological and range‐wide genetic variation data to detect such areas and to evaluate the impacts of habitat suitability changes on the genetic diversity of the transitional Mediterranean‐temperate tree Fraxinus angustifolia. We sampled and genotyped 38 natural populations comprising 1006 individuals from across Europe. We found the highest genetic diversity in western and northern Mediterranean populations, as well as a significant west to east decline in genetic diversity. Areas of potential refugia that correspond to approximately 70% of the suitable habitat may support the persistence of more than 90% of the total number of alleles in the future. Moreover, based on correlations between Bayesian genetic assignment and climate, climate change may favour the westward spread of the Black Sea gene pool in the long term. Overall, our results suggest that the northerly core areas of the current distribution contain the most important part of the genetic variation for this species and may serve as in situ macrorefugia from ongoing climate change. However, rear‐edge populations of the southern Mediterranean may be exposed to a potential loss of unique genetic diversity owing to habitat suitability changes unless populations can persist in microrefugia that have facilitated such persistence in the past.     

Population trends influence species ability to track climate change

Shifts of distributions have been attributed to species tracking their fundamental climate niches through space. However, several studies have now demonstrated that niche tracking is imperfect, that species’ climate niches may vary with population trends, and that geographic distributions may lag behind rapid climate change. These reports of imperfect niche tracking imply shifts in species’ realized climate niches. We argue that quantifying climate niche shifts and analyzing them for a suite of species reveal general patterns of niche shifts and the factors affecting species’ ability to track climate change. We analyzed changes in realized climate niche between 1984 and 2012 for 46 species of North American birds in relation to population trends in an effort to determine whether species differ in the ability to track climate change and whether differences in niche tracking are related to population trends. We found that increasingly abundant species tended to show greater levels of niche expansion (climate space occupied in 2012 but not in 1980) compared to declining species. Declining species had significantly greater niche unfilling (climate space occupied in 1980 but not in 2012) compared to increasing species due to an inability to colonize new sites beyond their range peripheries after climate had changed at sites of occurrence. Increasing species, conversely, were better able to colonize new sites and therefore showed very little niche unfilling. Our results indicate that species with increasing trends are better able to geographically track climate change compared to declining species, which exhibited lags relative to changes in climate. These findings have important implications for understanding past changes in distribution, as well as modeling dynamic species distributions in the face of climate change.     

Mechanistic niche modelling: combining physiological and spatial data to predict species' ranges

Species distribution models (SDMs) use spatial environmental data to make inferences on species' range limits and habitat suitability. Conceptually, these models aim to determine and map components of a species' ecological niche through space and time, and they have become important tools in pure and applied ecology and evolutionary biology. Most approaches are correlative in that they statistically link spatial data to species distribution records. An alternative strategy is to explicitly incorporate the mechanistic links between the functional traits of organisms and their environments into SDMs. Here, we review how the principles of biophysical ecology can be used to link spatial data to the physiological responses and constraints of organisms. This provides a mechanistic view of the fundamental niche which can then be mapped to the landscape to infer range constraints. We show how physiologically based SDMs can be developed for different organisms in different environmental contexts. Mechanistic SDMs have different strengths and weaknesses to correlative approaches, and there are many exciting and unexplored prospects for integrating the two approaches. As physiological knowledge becomes better integrated into SDMs, we will make more robust predictions of range shifts in novel or non-equilibrium contexts such as invasions, translocations, climate change and evolutionary shifts.     

Changing climate mediates sapsucker (Aves: Sphyrapicus) hybrid zone movement

Hybrid zones, where two divergent taxa meet and interbreed, offer unique opportunities to investigate how climate contributes to reproductive isolation between closely related taxa and how these taxa may respond to climatic changes. Red‐naped (Sphyrapicus nuchalis) and Red‐breasted (Sphyrapicus ruber) sapsuckers (Aves: Picidae) hybridize along a narrow contact zone that stretches from northern California to British Columbia. The hybrid zone between these species has been studied extensively for more than 100 years and represents an excellent system for investigations of the evolution of reproductive isolation. Shifts in the proportions of phenotypes at hybrid localities since 1910 that were inferred using specimens from museum collections were confirmed using species distribution models. We predicted the historical, current, and future distributions of parental and hybrid sapsuckers using Random Forests models to quantify how climate change is affecting hybrid zone movement in the Pacific Northwest. We found observed distribution shifts of parental sapsuckers were likely the result of climate change over the past 100 years, with these shifts predicted to continue for both sapsuckers over the next 80 years. We found Red‐breasted Sapsuckers are predicted to continue to expand, while Red‐naped Sapsuckers are predicted to contract substantially under future climate scenarios. As a result of the predicted changes, the amount of overlap in the distribution of these sapsuckers may decrease. Using hybrid phenotypes, we found the climate niche occupied by the hybrid zone is predicted to disappear under future conditions. The disappearance of this climate niche where the two parental species come into contact and hybridize may lead to a substantial reduction in genetic introgression. Understanding the impacts of global climate change on hybrid zones may help us to better understand how speciation has been shaped by climate in the past, as well as how evolution may respond to climate change in the future.     

Unravelling the evolutionary history and future prospects of endemic species restricted to former glacial refugia

The contemporary distribution and genetic composition of biodiversity bear a signature of species’ evolutionary histories and the effects of past climatic oscillations. For many European species, the Mediterranean peninsulas of Iberia, Italy and the Balkans acted as glacial refugia and the source of range recolonization, and as a result, they contain disproportionately high levels of diversity. As these areas are particularly threatened by future climate change, it is important to understand how past climatic changes affected their biodiversity. We use an integrated approach, combining markers with different evolutionary rates and combining phylogenetic analysis with approximate Bayesian computation and species distribution modelling across temporal scales. We relate phylogeographic processes to patterns of genetic variation in Myotis escalerai, a bat species endemic to the Iberian Peninsula. We found a distinct population structure at the mitochondrial level with a strong geographic signature, indicating lineage divergence into separate glacial refugia within the Iberian refugium. However, microsatellite markers suggest higher levels of gene flow resulting in more limited structure at recent time frames. The evolutionary history of M. escalerai was shaped by the effects of climatic oscillations and changes in forest cover and composition, while its future is threatened by climatically induced range contractions and the role of ecological barriers due to competition interactions in restricting its distribution. This study warns that Mediterranean peninsulas, which provided refuge for European biodiversity during past glaciation events, may become a trap for limited dispersal and ecologically limited endemic species under future climate change, resulting in loss of entire lineages.