Genetically informed ecological niche models improve climate change predictions

We examined the hypothesis that ecological niche models (ENMs) more accurately predict species distributions when they incorporate information on population genetic structure, and concomitantly, local adaptation. Local adaptation is common in species that span a range of environmental gradients (e.g., soils and climate). Moreover, common garden studies have demonstrated a covariance between neutral markers and functional traits associated with a species’ ability to adapt to environmental change. We therefore predicted that genetically distinct populations would respond differently to climate change, resulting in predicted distributions with little overlap. To test whether genetic information improves our ability to predict a species’ niche space, we created genetically informed ecological niche models (gENMs) using Populus fremontii (Salicaceae), a widespread tree species in which prior common garden experiments demonstrate strong evidence for local adaptation. Four major findings emerged: (i) gENMs predicted population occurrences with up to 12‐fold greater accuracy than models without genetic information; (ii) tests of niche similarity revealed that three ecotypes, identified on the basis of neutral genetic markers and locally adapted populations, are associated with differences in climate; (iii) our forecasts indicate that ongoing climate change will likely shift these ecotypes further apart in geographic space, resulting in greater niche divergence; (iv) ecotypes that currently exhibit the largest geographic distribution and niche breadth appear to be buffered the most from climate change. As diverse agents of selection shape genetic variability and structure within species, we argue that gENMs will lead to more accurate predictions of species distributions under climate change.     

No‐analog climates and shifting realized niches during the late quaternary: implications for 21st‐century predictions by species distribution models

Empirically derived species distributions models (SDMs) are increasingly relied upon to forecast species vulnerabilities to future climate change. However, many of the assumptions of SDMs may be violated when they are used to project species distributions across significant climate change events. In particular, SDM's in theory assume stable fundamental niches, but in practice, they assume stable realized niches. The assumption of a fixed realized niche relative to climate variables remains unlikely for various reasons, particularly if novel future climates open up currently unavailable portions of species’ fundamental niches. To demonstrate this effect, we compare the climate distributions for fossil‐pollen data from 21 to 15 ka bp (relying on paleoclimate simulations) when communities and climates with no modern analog were common across North America to observed modern pollen assemblages. We test how well SDMs are able to project 20th century pollen‐based taxon distributions with models calibrated using data from 21 to 15 ka. We find that taxa which were abundant in areas with no‐analog late glacial climates, such as Fraxinus, Ostrya/Carpinus and Ulmus, substantially shifted their realized niches from the late glacial period to present. SDMs for these taxa had low predictive accuracy when projected to modern climates despite demonstrating high predictive accuracy for late glacial pollen distributions. For other taxa, e.g. Quercus, Picea, Pinus strobus, had relatively stable realized niches and models for these taxa tended to have higher predictive accuracy when projected to present. Our findings reinforce the point that a realized niche at any one time often represents only a subset of the climate conditions in which a taxon can persist. Projections from SDMs into future climate conditions that are based solely on contemporary realized distributions are potentially misleading for assessing the vulnerability of species to future climate change.     

Conservatism of ecological niche characteristics in North American plant species over the Pleistocene-to-Recent transition

Aim  To provide a test of the conservatism of a species' niche over the last 20,000 years by tracking the distribution of eight pollen taxa relative to climate type as they migrated across eastern North America following the Last Glacial Maximum (LGM).
Location  North America.
Methods  We drew taxon occurrence data from the North American pollen records in the Global Pollen Database, representing eight pollen types – all taxa for which ≥5 distinct geographic occurrences were available in both the present day and at the LGM (21,000 years ago ± 3000 years). These data were incorporated into ecological niche models based on present-day and LGM climatological summaries available from the Palaeoclimate Modelling Intercomparison Project to produce predicted potential geographic distributions for each species at present and at the LGM. The output for each time period was projected onto the 'other' time period, and tested using independent known occurrence information from that period.
Results  The result of our analyses was that all species tested showed general conservatism in ecological characteristics over the climate changes associated with the Pleistocene-to-Recent transition.
Main conclusions  This analysis constitutes a further demonstration of general and pervasive conservatism in ecological niche characteristics over moderate periods of time despite profound changes in climate and environmental conditions. As such, our results reinforce the application of ecological niche modelling techniques to the reconstruction of Pleistocene biodiversity distribution patterns, and to project the future potential distribution range of species in the face of global-scale climatic changes.     

Using ecological niche modelling to evaluate niche stability in deep time

Aim To investigate relative niche stability in species responses to various types of environmental pressure (biotic and abiotic) on geological time‐scales using the fossil record. Location The case study focuses on Late Ordovician articulate brachiopods of the Cincinnati Arch in eastern North America. Methods Species niches were modelled for a suite of fossil brachiopod species based on five environmental variables inferred from sedimentary parameters using GARP and Maxent . Niche stability was assessed by comparison of (1) the degree of overlap of species distribution models developed for a time‐slice and those generated by projecting niche models of the previous time‐slice onto environmental layers of a second time‐slice using GARP and Maxent , (2) Schoener’s D statistic, and (3) the similarity of the contribution of each environmental parameter within Maxent niche models between adjacent time‐slices. Results Late Ordovician brachiopod species conserved their niches with high fidelity during intervals of gradual environmental change but responded to inter‐basinal species invasions through niche evolution. Both native and invasive species exhibited similar levels of niche evolution in the invasion and post‐invasion intervals. Niche evolution was related mostly to decreased variance within the former ecological niche parameters rather than to shifts to new ecospace. Main conclusions Although the species examined exhibited morphological stasis during the study interval, high levels of niche conservatism were observed only during intervals of gradual environmental change. Rapid environmental change, notably inter‐basinal species invasions, resulted in high levels of niche evolution among the focal taxa. Both native and invasive species responded with similar levels of niche evolution during the invasion interval and subsequent environmental reorganization. The assumption of complete niche conservatism frequently employed in ecological niche modelling (ENM) analyses to forecast or hindcast species geographical distributions is more likely to be accurate for climate change studies than for invasive species analyses over geological time‐scales.     

Causal or spurious relationship? Climate and the distribution of Phelsuma geckos on Grand Comoro Island

The volcanic island of Grand Comoro, Malagasy biogeographic region, is inhabited by three species of Phelsuma day geckos; two island‐endemic taxa (Phelsuma comorensis and Phelsuma v‐nigra comoraegrandensis) and the introduced Phelsuma dubia. Phelsuma comorensis is restricted to elevations of greater than 150 m above sea level on the northern of the island's two volcanoes and is the only Phelsuma above 300 m. The other species are widespread at low elevations but also reach levels above 900 m at the southern volcano. To investigate these divergent distribution patterns, we used environmental niche models based on climate and habitat data and tested whether predicted climate change may influence species distributions. Analyses of niche overlap did not show significant differences between present‐day and predicted future potential distributions of any Phelsuma species studied, which could be seen as an indicator of resilience towards climate change. Climate models reflected the restricted distribution of P. comorensis with precipitation of the wettest month detected as most important variable, whereas habitat models predicted an island‐wide distribution. While climate appears to determine the distribution of P. comorensis, we propose isolation by migration barriers as an alternative and discuss the detection of causal versus spurious relationships in ecological niche models.     

Measuring ecological niche overlap from occurrence and spatial environmental data

Aim Concerns over how global change will influence species distributions, in conjunction with increased emphasis on understanding niche dynamics in evolutionary and community contexts, highlight the growing need for robust methods to quantify niche differences between or within taxa. We propose a statistical framework to describe and compare environmental niches from occurrence and spatial environmental data. Location Europe, North America and South America. Methods The framework applies kernel smoothers to densities of species occurrence in gridded environmental space to calculate metrics of niche overlap and test hypotheses regarding niche conservatism. We use this framework and simulated species with pre‐defined distributions and amounts of niche overlap to evaluate several ordination and species distribution modelling techniques for quantifying niche overlap. We illustrate the approach with data on two well‐studied invasive species. Results We show that niche overlap can be accurately detected with the framework when variables driving the distributions are known. The method is robust to known and previously undocumented biases related to the dependence of species occurrences on the frequency of environmental conditions that occur across geographical space. The use of a kernel smoother makes the process of moving from geographical space to multivariate environmental space independent of both sampling effort and arbitrary choice of resolution in environmental space. However, the use of ordination and species distribution model techniques for selecting, combining and weighting variables on which niche overlap is calculated provide contrasting results. Main conclusions The framework meets the increasing need for robust methods to quantify niche differences. It is appropriate for studying niche differences between species, subspecies or intra‐specific lineages that differ in their geographical distributions. Alternatively, it can be used to measure the degree to which the environmental niche of a species or intra‐specific lineage has changed over time.     

Habitat availability and gene flow influence diverging local population trajectories under scenarios of climate change: a place‐based approach

Ecological niche theory holds that species distributions are shaped by a large and complex suite of interacting factors. Species distribution models (SDMs) are increasingly used to describe species’ niches and predict the effects of future environmental change, including climate change. Currently, SDMs often fail to capture the complexity of species’ niches, resulting in predictions that are generally limited to climate‐occupancy interactions. Here, we explore the potential impact of climate change on the American pika using a replicated place‐based approach that incorporates climate, gene flow, habitat configuration, and microhabitat complexity into SDMs. Using contemporary presence–absence data from occupancy surveys, genetic data to infer connectivity between habitat patches, and 21 environmental niche variables, we built separate SDMs for pika populations inhabiting eight US National Park Service units representing the habitat and climatic breadth of the species across the western United States. We then predicted occurrence probability under current (1981–2010) and three future time periods (out to 2100). Occurrence probabilities and the relative importance of predictor variables varied widely among study areas, revealing important local‐scale differences in the realized niche of the American pika. This variation resulted in diverse and – in some cases – highly divergent future potential occupancy patterns for pikas, ranging from complete extirpation in some study areas to stable occupancy patterns in others. Habitat composition and connectivity, which are rarely incorporated in SDM projections, were influential in predicting pika occupancy in all study areas and frequently outranked climate variables. Our findings illustrate the importance of a place‐based approach to species distribution modeling that includes fine‐scale factors when assessing current and future climate impacts on species’ distributions, especially when predictions are intended to manage and conserve species of concern within individual protected areas.     

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.     

A synthesis of transplant experiments and ecological niche models suggests that range limits are often niche limits

Global change has made it important to understand the factors that shape species’ distributions. Central to this area of research is the question of whether species’ range limits primarily reflect the distribution of suitable habitat (i.e. niche limits) or arise as a result of dispersal limitation. Over‐the‐edge transplant experiments and ecological niche models are commonly used to address this question, yet few studies have taken advantage of a combined approach for inferring the causes of range limits. Here, we synthesise results from existing transplant experiments with new information on the predicted suitability of sites based on niche models. We found that individual performance and habitat suitability independently decline beyond range limits across multiple species. Furthermore, inferences from transplant experiments and niche models were generally concordant within species, with 31 out of 40 cases fully supporting the hypothesis that range limits are niche limits. These results suggest that range limits are often niche limits and that the factors constraining species’ ranges operate at scales detectable by both transplant experiments and niche models. In light of these findings, we outline an integrative framework for addressing the causes of range limits in individual species.     

Complex patterns of environmental niche evolution in Iberian columbines (genus Aquilegia,Ranunculaceae)

Aims This study explores the patterns of niche differentiation in a group of seven closely related columbines (genus Aquilegia, Ranunculaceae) from the Iberian Peninsula. Populations of these columbines are subject to complex patterns of divergent selection across environments, which partly explain the taxonomic structure of the group. This suggests the hypothesis that niche divergence must have occurred along the process of diversification of the group.Methods We used MaxEnt to build environmental niche models of seven subspecies belonging to the three species of Aquilegia present in the Iberian Peninsula. From these models, we compared the environmental niches through two different approaches: ENMtools and multivariate methods.Important findings MaxEnt distributions conformed closely to the actual distribution of the study taxa. ENMtools methods failed to uncover any clear patterns of niche differentiation or conservatism in Iberian columbines. Multivariate analyses indicate the existence of differentiation along altitudinal gradients and along a gradient of climatic conditions determined by the summer precipitation and temperatures. However, climatic conditions related to winter temperature and precipitation, as well as soil properties, were equally likely to show conservatism or divergence. The complex patterns of niche evolution we found suggest that Iberian Columbines have not been significantly constrained by forces of niche conservatism, so they could respond adaptively to the fast and profound climate changes in the Iberian Peninsula through the glacial cycles of the Pleistocene.     

Effects of climate change on species turnover and body mass frequency distributions of South African bird communities

Climate change will lead to substantial shifts in species distributions. Most of the predictions of shifting distributions rely on modelling future distributions with ecological niche models. We used these models to investigate (i) the expected species turnover, loss and gain within bird communities of four South African biomes and (ii) the expected changes in the body mass frequency distributions of these communities. We used distributional data of the Southern African Bird Atlas Project, current climate data and two scenarios of future climate change for 2050 to build ensemble models of bird distributions. Our results indicate that future species loss, gain and turnover within the four biomes will be considerable. Climate change will also have statistically significant effects on body mass frequency distributions, and these effects differ substantially depending on the severity of future climate change. We discuss the possible ecological effects of these predicted changes on ecosystem interactions and functions.     

Ecological niche conservatism in North American freshwater fishes

There are many hypotheses of relationships, and also of speciation processes, in North American freshwater fishes, although, to date, there have been no direct tests of whether there is evidence of ecological niche conservatism. In the present study, ecological niche modeling is used to look for evidence of ecological niche conservatism in six clades of freshwater fishes: the starheaded topminnows, sand darters, black basses, Notropis rubellus species group, Notropis longirostris species group, and the Hybopsis amblops species group. This is achieved by evaluating the reciprocal predictivity of distributional predictions based on ecological niche models developed for each individual taxon in a clade under the assumption that high reciprical predictivity between sister species can be taken as evidence of niche conservatism. Omission percentages, total and average commission, and the area under the curve in a receiver operating characteristic analysis, where calculated, are used to evaluate predictive ability. Occurrence data for each species were subset into a training and independent validation data set where possible. Across all clades and species, models predicted the validation data for a given species well. Ecological niche conservatism was found generally across the starheaded topminnows, the sand darters, and the N. longirostris species group. There was some inter-predictivity within the N. rubellus group, but almost no inter-predictivity within the black basses, indicating a lack of conservatism. These results demonstrate that ecological niches generally act as stable constraints on freshwater fish distributions in North America.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 282–295.     

Integrating coalescent and ecological niche modeling in comparative phylogeography

Understanding the factors that contribute to the formation of population genetic structure is a central goal of phylogeographic research, but achieving this goal can be complicated by the stochastic variance inherent to genetic processes. Statistical approaches to testing phylogeographic hypotheses accommodate this stochasticity by evaluating competing models of putative historical population structure, often by simulating null distributions of the expected variance. The effectiveness of these tests depends on the biological realism of the models. Information from the fossil record can aid in reconstructing the historical distributions of some taxa. However, for the majority of taxa, which lack sufficient fossils, paleodistributional modeling can provide valuable spatial-geographic data concerning ancestral distributions. Paleodistributional models are generated by projecting ecological niche models, which predict the current distribution of each species, onto a model of past climatic conditions. Here, we generate paleodistributional models describing the suitable habitat during the last glacial maximum for lineages from the mesic forests of the Pacific Northwest of North America, and use these models to generate alternative phylogeographic hypotheses. Coalescent simulations are then used to test these hypotheses to improve our understanding of the historical events that promoted the formation of population genetic structure in this ecosystem. Results from Pacific Northwest mesic forest organisms demonstrate the utility of these combined approaches. Paleodistribution models and population genetic structure are congruent across three amphibian lineages, suggesting that they have responded in a concerted manner to environmental change. Two other species, a willow and a water vole, despite being currently codistributed and having similar population genetic structure, were predicted by the paleodistributional model to have had markedly different distributions during the last glacial maximum. This suggests that congruent phylogeographic patterns can arise from incongruent ancestral distributions. Paleodistributional models introduce a much-needed spatial-geographic perspective to statistical phylogeography. In conjunction with coalescent models of population genetic structure, they have the potential to improve our understanding of the factors that promote population divergence and ultimately produce regional patterns of biodiversity.     

Climate‐driven geographic distribution of the desert locust during recession periods: Subspecies’ niche differentiation and relative risks under scenarios of climate change

The desert locust is an agricultural pest that is able to switch from a harmless solitarious stage, during recession periods, to swarms of gregarious individuals that disperse long distances and affect areas from western Africa to India during outbreak periods. Large outbreaks have been recorded through centuries, and the Food and Agriculture Organization keeps a long‐term, large‐scale monitoring survey database in the area. However, there is also a much less known subspecies that occupies a limited area in Southern Africa. We used large‐scale climatic and occurrence data of the solitarious phase of each subspecies during recession periods to understand whether both subspecies climatic niches differ from each other, what is the current potential geographical distribution of each subspecies, and how climate change is likely to shift their potential distribution with respect to current conditions. We evaluated whether subspecies are significantly specialized along available climate gradients by using null models of background climatic differences within and between southern and northern ranges and applying niche similarity and niche equivalency tests. The results point to climatic niche conservatism between the two clades. We complemented this analysis with species distribution modeling to characterize current solitarious distributions and forecast potential recession range shifts under two extreme climate change scenarios at the 2050 and 2090 time horizon. Projections suggest that, at a global scale, the northern clade could contract its solitarious recession range, while the southern clade is likely to expand its recession range. However, local expansions were also predicted in the northern clade, in particular in southern and northern margins of the current geographical distribution. In conclusion, monitoring and management practices should remain in place in northern Africa, while in Southern Africa the potential for the subspecies to pose a threat in the future should be investigated more closely.     

Latitudinal gradients in climatic‐niche evolution accelerate trait evolution at high latitudes

Despite the importance of divergent selection to the speed of evolution, it remains poorly understood if divergent selection is more prevalent in the tropics (where species richness is highest), or at high latitudes (where paleoclimate change has been most intense). We tested whether the rate of climatic‐niche evolution – one proxy for divergent selection – varies with latitude for 111 pairs of bird species. Using Brownian motion and Ornsetin–Ulhenbeck models, we show that evolutionary rates along two important axes of the climatic‐niche – temperature and seasonality – have been faster at higher latitudes. We then tested whether divergence of the climatic‐niche was associated with evolution in traits important in ecological differentiation (body mass) and reproductive isolation (song), and found that climatic divergence is associated with faster rates in both measures. These results highlight the importance of climate‐mediated divergent selection pressures in driving evolutionary divergence and reproductive isolation at high latitudes.     

Rates of projected climate change dramatically exceed past rates of climatic niche evolution among vertebrate species

A key question in predicting responses to anthropogenic climate change is: how quickly can species adapt to different climatic conditions? Here, we take a phylogenetic approach to this question. We use 17 time‐calibrated phylogenies representing the major tetrapod clades (amphibians, birds, crocodilians, mammals, squamates, turtles) and climatic data from distributions of > 500 extant species. We estimate rates of change based on differences in climatic variables between sister species and estimated times of their splitting. We compare these rates to predicted rates of climate change from 2000 to 2100. Our results are striking: matching projected changes for 2100 would require rates of niche evolution that are > 10 000 times faster than rates typically observed among species, for most variables and clades. Despite many caveats, our results suggest that adaptation to projected changes in the next 100 years would require rates that are largely unprecedented based on observed rates among vertebrate species.     

Predictability of stream insect distributions is dependent on niche position,but not on biological traits or taxonomic relatedness of species

Species distributions can be analysed under two perspectives: the niche‐based approach, which focuses on species–environment relationships; and the dispersal‐based approach, which focuses on metapopulation dynamics. The degree to which each of these two components affect species distributions may depend on habitat fragmentation, species traits and phylogenetic constraints. We analysed the distributions of 36 stream insect species across 60 stream sites in three drainage basins at high latitudes in Finland. We used binomial generalised linear models (GLMs) in which the predictor variables were environmental factors (E models), within‐basin spatial variables as defined by Moran's eigenvector maps (M models), among‐basin variability (B models), or a combination of the three (E + M + B models) sets of variables. Based on a comparative analysis, model performance was evaluated across all the species using Gaussian GLMs whereby the deviance accounted for by binomial GLMs was fitted on selected explanatory variables: niche position, niche breadth, site occupancy, biological traits and taxonomic relatedness. For each type of model, a reduced Gaussian GLM was eventually obtained after variable selection (Bayesian information criterion). We found that niche position was the only variable selected in all reduced models, implying that marginal species were better predicted than non‐marginal species. The influence of niche position was strongest in models based on environmental variables (E models) or a combination of all types of variables (E + M + B models), and weakest in spatial autocorrelation models (M models). This suggests that species–environment relationships prevail over dispersal processes in determining stream insect distributions at a regional scale. Our findings have clear implications for biodiversity conservation strategies, and they also emphasise the benefits of considering both the niche‐based and dispersal‐based approaches in species distribution modelling studies.     

Evaluating niche changes during invasion with seasonal models in Capsella bursa-pastoris

Premise

Researchers often use ecological niche models to predict where species might establish and persist under future or novel climate conditions. However, these predictive methods assume species have stable niches across time and space. Furthermore, ignoring the time of occurrence data can obscure important information about species reproduction and ultimately fitness. Here, we assess compare ecological niche models generated from full-year averages to seasonal models.

Methods

In this study, we generate full-year and monthly ecological niche models for Capsella bursa-pastoris in Europe and North America to see if we can detect changes in the seasonal niche of the species after long-distance dispersal.

Results

We find full-year ecological niche models have low transferability across continents and there are continental differences in the climate conditions that influence the distribution of C. bursa-pastoris. Monthly models have greater predictive accuracy than full-year models in cooler seasons, but no monthly models can predict North American summer occurrences very well.

Conclusions

The relative predictive ability of European monthly models compared to North American monthly models suggests a change in the seasonal timing between the native range to the non-native range. These results highlight the utility of ecological niche models at finer temporal scales in predicting species distributions and unmasking subtle patterns of evolution.