Better forecasts of range dynamics using genetic data

Drivers of genomic differentiation landscapes in populations of disparate ecological and geographical settings within mainland Apis cerana

Elucidating the evolutionary processes that drive population divergence can enhance our understanding of the early stages of speciation and inform conservation management decisions. The honeybee Apis cerana displays extensive population divergence, providing an informative natural system for exploring these processes. The mainland lineage A. cerana includes several peripheral subspecies with disparate ecological and geographical settings radiated from a central ancestor. Under this evolutionary framework, we can explore the patterns of genome differentiation and the evolutionary models that explain them. We can also elucidate the contribution of non-genomic spatiotemporal mechanisms (extrinsic features) and genomic mechanisms (intrinsic features) that influence these genomic differentiation landscapes. Based on 293 whole genomes, a small part of the genome is highly differentiated between central-peripheral subspecies pairs, while low and partial parallelism partly reflects idiosyncratic responses to environmental differences. Combined elements of recurrent selection and speciation-with-gene-flow models generate the heterogeneous genome landscapes. These elements weight differently between central-island and other central-peripheral subspecies pairs, influenced by glacial cycles superimposed on different geomorphologies. Although local recombination rates exert a significant influence on patterns of genomic differentiation, it is unlikely that low-recombination rates regions were generated by structural variation. In conclusion, complex factors including geographical isolation, divergent ecological selection and non-uniform genome features have acted concertedly in the evolution of reproductive barriers that could reduce gene flow in part of the genome and facilitate the persistence of distinct populations within mainland lineage of A. cerana.

Reconstructing colonization dynamics to establish how human activities transformed island biodiversity

Drivers and dynamics of initial human migrations across individual islands and archipelagos are poorly understood, hampering assessments of subsequent modification of island biodiversity. We developed and tested a new statistical-simulation approach for reconstructing the pattern and pace of human migration across islands at high spatiotemporal resolutions. Using Polynesian colonisation of New Zealand as an example, we show that process-explicit models, informed by archaeological records and spatiotemporal reconstructions of past climates and environments, can provide new and important insights into the patterns and mechanisms of arrival and establishment of people on islands. We find that colonisation of New Zealand required there to have been a single founding population of approximately 500 people, arriving between 1233 and 1257 AD, settling multiple areas, and expanding rapidly over both North and South Islands. These verified spatiotemporal reconstructions of colonisation dynamics provide new opportunities to explore more extensively the potential ecological impacts of human colonisation on New Zealand's native biota and ecosystems.

Millennial processes of population decline, range contraction and near extinction of the European bison

European bison (Bison bonasus) were widespread throughout Europe during the late Pleistocene. However, the contributions of environmental change and humans to their near extinction have never been resolved. Using process-explicit models, fossils and ancient DNA, we disentangle the combinations of threatening processes that drove population declines and regional extinctions of European bison through space and across time. We show that the population size of European bison declined abruptly at the termination of the Pleistocene in response to rapid environmental change, hunting by humans and their interaction. Human activities prevented populations of European bison from rebounding in the Holocene, despite improved environmental conditions. Hunting caused range loss in the north and east of its distribution, while land use change was responsible for losses in the west and south. Advances in hunting technologies from 1500 CE were needed to simulate low abundances observed in 1870 CE. While our findings show that humans were an important driver of the extinction of the European bison in the wild, vast areas of its range vanished during the Pleistocene-Holocene transition because of post-glacial environmental change. These areas of its former range have been climatically unsuitable for millennia and should not be considered in reintroduction efforts.

Using landscape genomics to delineate future adaptive potential for climate change in the Yosemite toad (Anaxyrus canorus)

An essential goal in conservation biology is delineating population units that maximize the probability of species persisting into the future and adapting to future environmental change. However, future-facing conservation concerns are often addressed using retrospective patterns that could be irrelevant. We recommend a novel landscape genomics framework for delineating future "Geminate Evolutionary Units" (GEUs) in a focal species: (1) identify loci under environmental selection, (2) model and map adaptive conservation units that may spawn future lineages, (3) forecast relative selection pressures on each future lineage, and (4) estimate their fitness and likelihood of persistence using geo-genomic simulations. Using this process, we delineated conservation units for the Yosemite toad (Anaxyrus canorus), a U.S. federally threatened species that is highly vulnerable to climate change. We used a genome-wide dataset, redundancy analysis, and Bayesian association methods to identify 24 candidate loci responding to climatic selection (R 2 ranging from 0.09 to 0.52), after controlling for demographic structure. Candidate loci included genes such as MAP3K5, involved in cellular response to environmental change. We then forecasted future genomic response to climate change using the multivariate machine learning algorithm Gradient Forests. Based on all available evidence, we found three GEUs in Yosemite National Park, reflecting contrasting adaptive optima: YF-North (high winter snowpack with moderate summer rainfall), YF-East (low to moderate snowpack with high summer rainfall), and YF-Low-Elevation (low snowpack and rainfall). Simulations under the RCP 8.5 climate change scenario suggest that the species will decline by 29% over 90 years, but the highly diverse YF-East lineage will be least impacted for two reasons: (1) geographically it will be sheltered from the largest climatic selection pressures, and (2) its standing genetic diversity will promote a faster adaptive response. Our approach provides a comprehensive strategy for protecting imperiled non-model species with genomic data alone and has wide applicability to other declining species.

New insights on patterns of genetic admixture and phylogeographic history in Iberian high mountain populations of midwife toads

Multiple Quaternary glacial refugia in the Iberian Peninsula, commonly known as "refugia within refugia", allowed diverging populations to come into contact and admix, potentially boosting substantial mito-nuclear discordances. In this study, we employ a comprehensive set of mitochondrial and nuclear markers to shed light onto the drivers of geographical differentiation in Iberian high mountain populations of the midwife toads Alytes obstetricans and A. almogavarii from the Pyrenees, Picos de Europa and Guadarrama Mountains. In the three analysed mountain regions, we detected evidence of extensive mito-nuclear discordances and/or admixture between taxa. Clustering analyses identified three major divergent lineages in the Pyrenees (corresponding to the eastern, central and central-western Pyrenees), which possibly recurrently expanded and admixed during the succession of glacial-interglacial periods that characterised the Late Pleistocene, and that currently follow a ring-shaped diversification pattern. On the other hand, populations from the Picos de Europa mountains (NW Iberian Peninsula) showed a mitochondrial affinity to central-western Pyrenean populations and a nuclear affinity to populations from the central Iberian Peninsula, suggesting a likely admixed origin for Picos de Europa populations. Finally, populations from the Guadarrama Mountain Range (central Iberian Peninsula) were depleted of genetic diversity, possibly as a consequence of a recent epidemic of chytridiomycosis. This work highlights the complex evolutionary history that shaped the current genetic composition of high mountain populations, and underscores the importance of using a multilocus approach to better infer the dynamics of population divergence.

Spatiotemporal influences of climate and humans on muskox range dynamics over multiple millennia

Processes leading to range contractions and population declines of Arctic megafauna during the late Pleistocene and early Holocene are uncertain, with intense debate on the roles of human hunting, climatic change, and their synergy. Obstacles to a resolution have included an overreliance on correlative rather than process-explicit approaches for inferring drivers of distributional and demographic change. Here, we disentangle the ecological mechanisms and threats that were integral in the decline and extinction of the muskox (Ovibos moschatus) in Eurasia and in its expansion in North America using process-explicit macroecological models. The approach integrates modern and fossil occurrence records, ancient DNA, spatiotemporal reconstructions of past climatic change, species-specific population ecology, and the growth and spread of anatomically modern humans. We show that accurately reconstructing inferences of past demographic changes for muskox over the last 21,000 years require high dispersal abilities, large maximum densities, and a small Allee effect. Analyses of validated process-explicit projections indicate that climatic change was the primary driver of muskox distribution shifts and demographic changes across its previously extensive (circumpolar) range, with populations responding negatively to rapid warming events. Regional analyses show that the range collapse and extinction of the muskox in Europe (~13,000 years ago) was likely caused by humans operating in synergy with climatic warming. In Canada and Greenland, climatic change and human activities probably combined to drive recent population sizes. The impact of past climatic change on the range and extinction dynamics of muskox during the Pleistocene-Holocene transition signals a vulnerability of this species to future increased warming. By better establishing the ecological processes that shaped the distribution of the muskox through space and time, we show that process-explicit macroecological models have important applications for the future conservation and management of this iconic species in a warming Arctic.

The antique genetic plight of the Mediterranean monk seal (Monachus monachus)

Disentangling the impact of Late Quaternary climate change from human activities can have crucial implications on the conservation of endangered species. We investigated the population genetics and demography of the Mediterranean monk seal (Monachus monachus), one of the world's most endangered marine mammals, through an unprecedented dataset encompassing historical (extinct) and extant populations from the eastern North Atlantic to the entire Mediterranean Basin. We show that Cabo Blanco (Western Sahara/Mauritania), Madeira, Western Mediterranean (historical range) and Eastern Mediterranean regions segregate into four populations. This structure is probably the consequence of recent drift, combined with long-term isolation by distance (R2 = 0.7), resulting from prevailing short-distance (less than 500 km) and infrequent long-distance dispersal (less than 1500 km). All populations (Madeira especially), show high levels of inbreeding and low levels of genetic diversity, seemingly declining since historical time, but surprisingly not being impacted by the 1997 massive die-off in Cabo Blanco. Approximate Bayesian Computation analyses support scenarios combining local extinctions and a major effective population size decline in all populations during Antiquity. Our results suggest that the early densification of human populations around the Mediterranean Basin coupled with the development of seafaring techniques were the main drivers of the decline of Mediterranean monk seals.

Process-explicit models reveal the structure and dynamics of biodiversity patterns

With ever-growing data availability and computational power at our disposal, we now have the capacity to use process-explicit models more widely to reveal the ecological and evolutionary mechanisms responsible for spatiotemporal patterns of biodiversity. Most research questions focused on the distribution of diversity cannot be answered experimentally, because many important environmental drivers and biological constraints operate at large spatiotemporal scales. However, we can encode proposed mechanisms into models, observe the patterns they produce in virtual environments, and validate these patterns against real-world data or theoretical expectations. This approach can advance understanding of generalizable mechanisms responsible for the distributions of organisms, communities, and ecosystems in space and time, advancing basic and applied science. We review recent developments in process-explicit models and how they have improved knowledge of the distribution and dynamics of life on Earth, enabling biodiversity to be better understood and managed through a deeper recognition of the processes that shape genetic, species, and ecosystem diversity.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Climate change effects on animal ecology: butterflies and moths as a case study

Butterflies and moths (Lepidoptera) are one of the most studied, diverse, and widespread animal groups, making them an ideal model for climate change research. They are a particularly informative model for studying the effects of climate change on species ecology because they are ectotherms that thermoregulate with a suite of physiological, behavioural, and phenotypic traits. While some species have been negatively impacted by climatic disturbances, others have prospered, largely in accordance with their diversity in life-history traits. Here we take advantage of a large repertoire of studies on butterflies and moths to provide a review of the many ways in which climate change is impacting insects, animals, and ecosystems. By studying these climate-based impacts on ecological processes of Lepidoptera, we propose appropriate strategies for species conservation and habitat management broadly across animals.

StableClim, continuous projections of climate stability from 21000 BP to 2100 CE at multiple spatial scales

Paleoclimatic data are used in eco-evolutionary models to improve knowledge of biogeographical processes that drive patterns of biodiversity through time, opening windows into past climate–biodiversity dynamics. Applying these models to harmonised simulations of past and future climatic change can strengthen forecasts of biodiversity change. StableClim provides continuous estimates of climate stability from 21,000 years ago to 2100 C.E. for ocean and terrestrial realms at spatial scales that include biogeographic regions and climate zones. Climate stability is quantified using annual trends and variabilities in air temperature and precipitation, and associated signal-to-noise ratios. Thresholds of natural variability in trends in regional- and global-mean temperature allow periods in Earth’s history when climatic conditions were warming and cooling rapidly (or slowly) to be identified and climate stability to be estimated locally (grid-cell) during these periods of accelerated change. Model simulations are validated against independent paleoclimate and observational data. Projections of climatic stability, accessed through StableClim, will improve understanding of the roles of climate in shaping past, present-day and future patterns of biodiversity.

Measurement(s)	climate change • climate • temperature of air • volume of hydrological precipitation
Technology Type(s)	computational modeling technique • digital curation
Factor Type(s)	timing of temperature and precipitation estimates
Sample Characteristic - Environment	climate system
Sample Characteristic - Location	Earth (planet)

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12831935

Using paleo-archives to safeguard biodiversity under climate change

Strategies for 21st-century environmental management and conservation under global change require a strong understanding of the biological mechanisms that mediate responses to climate- and human-driven change to successfully mitigate range contractions, extinctions, and the degradation of ecosystem services. Biodiversity responses to past rapid warming events can be followed in situ and over extended periods, using cross-disciplinary approaches that provide cost-effective and scalable information for species' conservation and the maintenance of resilient ecosystems in many bioregions. Beyond the intrinsic knowledge gain such integrative research will increasingly provide the context, tools, and relevant case studies to assist in mitigating climate-driven biodiversity losses in the 21st century and beyond.

Harnessing paleo-environmental modeling and genetic data to predict intraspecific genetic structure

Spatially explicit simulations of gene flow within complex landscapes could help forecast the responses of populations to global and anthropological changes. Simulating how past climate change shaped intraspecific genetic variation can provide a validation of models in anticipation of their use to predict future changes. We review simulation models that provide inferences on population genetic structure. Existing simulation models generally integrate complex demographic and genetic processes but are less focused on the landscape dynamics. In contrast to previous approaches integrating detailed demographic and genetic processes and only secondarily landscape dynamics, we present a model based on parsimonious biological mechanisms combining habitat suitability and cellular processes, applicable to complex landscapes. The simulation model takes as input (a) the species dispersal capacities as the main biological parameter, (b) the species habitat suitability, and (c) the landscape structure, modulating dispersal. Our model emphasizes the role of landscape features and their temporal dynamics in generating genetic differentiation among populations within species. We illustrate our model on caribou/reindeer populations sampled across the entire species distribution range in the Northern Hemisphere. We show that simulations over the past 21 kyr predict a population genetic structure that matches empirical data. This approach looking at the impact of historical landscape dynamics on intraspecific structure can be used to forecast population structure under climate change scenarios and evaluate how species range shifts might induce erosion of genetic variation within species.

Insights into the neutral and adaptive processes shaping the spatial distribution of genomic variation in the economically important Moroccan locust (Dociostaurus maroccanus)

Understanding the processes that shape neutral and adaptive genomic variation is a fundamental step to determine the demographic and evolutionary dynamics of pest species. Here, we use genomic data obtained via restriction site-associated DNA sequencing to investigate the genetic structure of Moroccan locust (Dociostaurus maroccanus) populations from the westernmost portion of the species distribution (Iberian Peninsula and Canary Islands), infer demographic trends, and determine the role of neutral versus selective processes in shaping spatial patterns of genomic variation in this pest species of great economic importance. Our analyses showed that Iberian populations are characterized by high gene flow, whereas the highly isolated Canarian populations have experienced strong genetic drift and loss of genetic diversity. Historical demographic reconstructions revealed that all populations have passed through a substantial genetic bottleneck around the last glacial maximum (~21 ka BP) followed by a sharp demographic expansion at the onset of the Holocene, indicating increased effective population sizes during warm periods as expected from the thermophilic nature of the species. Genome scans and environmental association analyses identified several loci putatively under selection, suggesting that local adaptation processes in certain populations might not be impeded by widespread gene flow. Finally, all analyses showed few differences between outbreak and nonoutbreak populations. Integrated pest management practices should consider high population connectivity and the potential importance of local adaptation processes on population persistence.

Climate connectivity of the bobcat in the Great Lakes region

The Great Lakes and the St. Lawrence River are imposing barriers for wildlife, and the additive effect of urban and agricultural development that dominates the lower Great Lakes region likely further reduces functional connectivity for many terrestrial species. As the climate warms, species will need to track climate across these barriers. It is important therefore to investigate land cover and bioclimatic hypotheses that may explain the northward expansion of species through the Great Lakes. We investigated the functional connectivity of a vagile generalist, the bobcat, as a representative generalist forest species common to the region. We genotyped tissue samples collected across the region at 14 microsatellite loci and compared different landscape hypotheses that might explain the observed gene flow or functional connectivity. We found that the Great Lakes and the additive influence of forest stands with either low or high canopy cover and deep lake-effect snow have disrupted gene flow, whereas intermediate forest cover has facilitated gene flow. Functional connectivity in southern Ontario is relatively low and was limited in part by the low amount of forest cover. Pathways across the Great Lakes were through the Niagara region and through the Lower Peninsula of Michigan over the Straits of Mackinac and the St. Marys River. These pathways are important routes for bobcat range expansion north of the Great Lakes and are also likely pathways that many other mobile habitat generalists must navigate to track the changing climate. The extent to which species can navigate these routes will be important for determining the future biodiversity of areas north of the Great Lakes.

Genes on the edge: A framework to detect genetic diversity imperiled by climate change

Ongoing global warming is disrupting several ecological and evolutionary processes, spanning different levels of biological organization. Species are expected to shift their ranges as a response to climate change, with relevant implications to peripheral populations at the trailing and leading edges. Several studies have analyzed the exposure of species to climate change but few have explored exposure at the intraspecific level. We introduce a framework to forecast exposure to climate change at the intraspecific level. We build on existing methods by combining correlative species distribution models, a model of species range dynamics, and a model of phylogeographic interpolation. We demonstrate the framework by applying it to 20 Iberian amphibian and reptile species. Our aims were to: (a) identify which species and intraspecific lineages will be most exposed to future climate change; (b) test if nucleotide diversity at the edges of species ranges are significantly higher or lower than on the overall range; and (c) analyze if areas of higher species gain, loss, and turnover coincide with those predicted for lineages richness and nucleotide diversity. We found that about 80% of the studied species are predicted to contract their range. Within each species, some lineages were predicted to contract their range, while others were predicted to maintain or expand it. Therefore, estimating the impacts of climate change at the species level only can underestimate losses at the intraspecific level. Some species had significant high amount of nucleotide at the trailing or leading edge, or both, but we did not find a consistent pattern across species. Spatial patterns of species richness, gain, loss, and turnover were fairly concurrent with lineages richness and nucleotide diversity. Our results support the need for increased attention to intraspecific diversity regarding monitoring and conservation strategies under climate change.

Implications of existing local (mal)adaptations for ecological forecasting under environmental change

Standing genetic variation represents a genetic load on population fitness but can also support a rapid response to short-term environmental change, and the greatest potential source of such standing genetic variation typically exists among locally adapted populations living along an environmental gradient. Here, we develop a spatially explicit simulation model to quantify the contribution of existing genetic variation arising from migration-mutation-selection-drift balance to time to extinction under environmental change. Simulations reveal that local adaptation across a species range associated with an underlying environmental gradient could extend time to extinction by nearly threefold irrespective of the rate of environmental change. The potential for preadapted alleles to increase the rate of adaptation changes the relative importance of established extinction risk factors; in particular, it reduced the importance of the breadth of environmental tolerance and it increased the relative importance of fecundity. Although migration of preadapted alleles generally increased persistence time, it decreased it at rates of environmental change close to the critical rate of change by creating a population bottleneck, which ultimately limited the rate at which de novo mutations could arise. An analysis of the extinction dynamics further revealed that one consequence of gene flow is the potential to maximize population growth rate in at least part of the species range, which is likely to have consequences for forecasting the consequences of ecological interactions. Our study shows that predictions of persistence time change fundamentally when existing local adaptations are explicitly taken into account, underscoring the need to preserve and manage genetic diversity.

Phylogeography and niche modelling: reciprocal enlightenment

Phylogeography examines the spatial genetic structure of species. Environmental niche modelling (or ecological niche modelling; ENM) examines the environmental limits of a species’ ecological niche. These two fields have great potential to be used together. ENM can shed light on how phylogeographical patterns develop and help identify possible drivers of spatial structure that need to be further investigated. Specifically, ENM can be used to test for niche differentiation among clades, identify factors limiting individual clades and identify barriers and contact zones. It can also be used to test hypotheses regarding the effects of historical and future climate change on spatial genetic patterns by projecting niches using palaeoclimate or future climate data. Conversely, phylogeographical information can populate ENM with within-species genetic diversity. Where adaptive variation exists among clades within a species, modelling their niches separately can improve predictions of historical distribution patterns and future responses to climate change. Awareness of patterns of genetic diversity in niche modelling can also alert conservationists to the potential loss of genetically diverse areas in a species’ range. Here, we provide a simplistic overview of both fields, and focus on their potential for integration, encouraging researchers on both sides to take advantage of the opportunities available.

Nunataks or massif de refuge? A phylogeographic study of Rhodiola crenulata (Crassulaceae) on the world's highest sky islands

Background

Quaternary climatic oscillations had tremendous effects on the current distribution of species. Here, we aim to elucidate the glacial history of Rhodiola crenulata, a perennial herb almost exclusively restricted to rock crevices on mountain peaks, and to test whether the nunatak or massif de refuge hypotheses could explain its distribution pattern.

Results

Six haplotypes and six ribotypes were detected in the cpDNA data set and the ITS data set, respectively. The divergence of R. crenulata and its closest relatives was dated have occurred ca. 0.65 Mya, during the Naynayxungla glaciation on the QTP. Mismatch distribution analysis suggested that the species experienced a range expansion around 0.31 Mya. Populations with high genetic and haplotype diversity were found on the QTP platform as well in the Hengduan Mountains. The ecological niche modeling results showed that there were suitable habitats on both the QTP platform and in the Hengduan Mountains during the LGM.

Conclusion

Our results support a scenario that both nunataks and the massif de refuge hypotheses could explain the distribution of R. crenulata. We also confirmed that Quaternary climatic oscillations could promote plant speciation in some circumstances. This study adds to a growing body of evidence suggesting that the QTP plant lineages exhibited diverse reactions to the Quaternary climatic oscillations.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1270-6) contains supplementary material, which is available to authorized users.

Cracking the Code of Biodiversity Responses to Past Climate Change

How individual species and entire ecosystems will respond to future climate change are among the most pressing questions facing ecologists. Past biodiversity dynamics recorded in the paleoecological archives show a broad array of responses, yet significant knowledge gaps remain. In particular, the relative roles of evolutionary adaptation, phenotypic plasticity, and dispersal in promoting survival during times of climate change have yet to be clarified. Investigating the paleo-archives offers great opportunities to understand biodiversity responses to future climate change. In this review we discuss the mechanisms by which biodiversity responds to environmental change, and identify gaps of knowledge on the role of range shifts and tolerance. We also outline approaches at the intersection of paleoecology, genomics, experiments, and predictive models that will elucidate the processes by which species have survived past climatic changes and enhance predictions of future changes in biological diversity.

Why decadal to century timescale palaeoclimate data are needed to explain present-day patterns of biological diversity and change

The current distribution of species, environmental conditions and their interactions represent only one snapshot of a planet that is continuously changing, in part due to human influences. To distinguish human impacts from natural factors, the magnitude and pace of climate shifts, since the Last Glacial Maximum, are often used to determine whether patterns of diversity today are artefacts of past climate change. In the absence of high-temporal resolution palaeoclimate reconstructions, this is generally done by assuming that past climate change occurred at a linear pace between widely spaced (usually, ≥1,000 years) climate snapshots. We show here that this is a flawed assumption because regional climates have changed significantly across decades and centuries during glacial-interglacial cycles, likely causing rapid regional replacement of biota. We demonstrate how recent atmosphere-ocean general circulation model (AOGCM) simulations of the climate of the past 21,000 years can provide credible estimates of the details of climate change on decadal to centennial timescales, showing that these details differ radically from what might be inferred from longer timescale information. High-temporal resolution information can provide more meaningful estimates of the magnitude and pace of climate shifts, the location and timing of drivers of physiological stress, and the extent of novel climates. They also produce new opportunities to directly investigate whether short-term climate variability is more important in shaping biodiversity patterns rather than gradual changes in long-term climatic means. Together, these more accurate measures of past climate instability are likely to bring about a better understanding of the role of palaeoclimatic change and variability in shaping current macroecological patterns in many regions of the world.

How complex should models be? Comparing correlative and mechanistic range dynamics models

Criticism has been levelled at climate-change-induced forecasts of species range shifts that do not account explicitly for complex population dynamics. The relative importance of such dynamics under climate change is, however, undetermined because direct tests comparing the performance of demographic models vs. simpler ecological niche models are still lacking owing to difficulties in evaluating forecasts using real-world data. We provide the first comparison of the skill of coupled ecological-niche-population models and ecological niche models in predicting documented shifts in the ranges of 20 British breeding bird species across a 40-year period. Forecasts from models calibrated with data centred on 1970 were evaluated using data centred on 2010. We found that more complex coupled ecological-niche-population models (that account for dispersal and metapopulation dynamics) tend to have higher predictive accuracy in forecasting species range shifts than structurally simpler models that only account for variation in climate. However, these better forecasts are achieved only if ecological responses to climate change are simulated without static snapshots of historic land use, taken at a single point in time. In contrast, including both static land use and dynamic climate variables in simpler ecological niche models improve forecasts of observed range shifts. Despite being less skilful at predicting range changes at the grid-cell level, ecological niche models do as well, or better, than more complex models at predicting the magnitude of relative change in range size. Therefore, ecological niche models can provide a reasonable first approximation of the magnitude of species' potential range shifts, especially when more detailed data are lacking on dispersal dynamics, demographic processes underpinning population performance, and change in land cover.

Genetic variation during range expansion: effects of habitat novelty and hybridization

How species' ranges evolve remains an enduring problem in ecology and evolutionary biology. Species' range limits are potentially set by the inability of peripheral populations to adapt to range-edge habitat. Indeed, peripheral populations are often assumed to have reduced genetic diversity and population sizes, which limit evolvability. However, support for this assumption is mixed, possibly because the genetic effects of range expansion depend on two factors: the extent that habitat into which expansion occurs is novel and sources of gene flow. Here, we used spadefoot toads, Spea bombifrons, to contrast the population genetic effects of expansion into novel versus non-novel habitat. We further evaluated gene flow from conspecifics and from heterospecifics via hybridization with a resident species. We found that range expansion into novel habitat, relative to non-novel habitat, resulted in higher genetic differentiation, lower conspecific gene flow and bottlenecks. Moreover, we found that hybridizing with a resident species introduced genetic diversity in the novel habitat. Our results suggest the evolution of species' ranges can depend on the extent of differences in habitat between ancestral and newly occupied ranges. Furthermore, our results highlight the potential for hybridization with a resident species to enhance genetic diversity during expansions into novel habitat.

Bioclimatic transect networks: Powerful observatories of ecological change

Transects that traverse substantial climate gradients are important tools for climate change research and allow questions on the extent to which phenotypic variation associates with climate, the link between climate and species distributions, and variation in sensitivity to climate change among biomes to be addressed. However, the potential limitations of individual transect studies have recently been highlighted. Here, we argue that replicating and networking transects, along with the introduction of experimental treatments, addresses these concerns. Transect networks provide cost‐effective and robust insights into ecological and evolutionary adaptation and improve forecasting of ecosystem change. We draw on the experience and research facilitated by the Australian Transect Network to demonstrate our case, with examples, to clarify how population‐ and community‐level studies can be integrated with observations from multiple transects, manipulative experiments, genomics, and ecological modeling to gain novel insights into how species and systems respond to climate change. This integration can provide a spatiotemporal understanding of past and future climate‐induced changes, which will inform effective management actions for promoting biodiversity resilience.

Testing the role of ancient and contemporary landscapes on structuring genetic variation in a specialist grasshopper

Understanding the processes underlying spatial patterns of genetic diversity and structure of natural populations is a central topic in evolutionary biogeography. In this study, we combine data on ancient and contemporary landscape composition to get a comprehensive view of the factors shaping genetic variation across the populations of the scrub-legume grasshopper (Chorthippus binotatus binotatus) from the biogeographically complex region of southeast Iberia. First, we examined geographical patterns of genetic structure and employed an approximate Bayesian computation (ABC) approach to compare different plausible scenarios of population divergence. Second, we used a landscape genetic framework to test for the effects of (1) Late Miocene paleogeography, (2) Pleistocene climate fluctuations, and (3) contemporary topographic complexity on the spatial patterns of population genetic differentiation. Genetic structure and ABC analyses supported the presence of three genetic clusters and a sequential west-to-east splitting model that predated the last glacial maximum (LGM, c. 21 Kya). Landscape genetic analyses revealed that population genetic differentiation was primarily shaped by contemporary topographic complexity, but was not explained by any paleogeographic scenario or resistance distances based on climate suitability in the present or during the LGM. Overall, this study emphasizes the need of integrating information on ancient and contemporary landscape composition to get a comprehensive view of their relative importance to explain spatial patterns of genetic variation in organisms inhabiting regions with complex biogeographical histories.

Life on the rocks: Multilocus phylogeography of rock hyrax (Procavia capensis) from southern Africa

Understanding the role of geography and climatic cycles in determining patterns of biodiversity is important in comparative and evolutionary biology and conservation. We studied the phylogeographic pattern and historical demography of a rock-dwelling small mammal species from southern Africa, the rock hyrax Procavia capensis capensis. Using a multilocus coalescent approach, we assessed the influence of strong habitat dependence and fluctuating regional climates on genetic diversity. We sequenced a mitochondrial gene (cytochrome b) and two nuclear introns (AP5, PRKC1) supplemented with microsatellite genotyping, in order to assess evolutionary processes over multiple temporal scales. In addition, distribution modelling was used to investigate the current and predicted distribution of the species under different climatic scenarios. Collectively, the data reveal a complex history of isolation followed by secondary contact shaping the current intraspecific diversity. The cyt b sequences confirmed the presence of two previously proposed geographically and genetically distinct lineages distributed across the southern African Great Escarpment and north-western mountain ranges. Molecular dating suggests Miocene divergence of the lineages, yet there are no discernible extrinsic barriers to gene flow. The nuclear markers reveal incomplete lineage sorting or ongoing mixing of the two lineages. Although the microsatellite data lend some support to the presence of two subpopulations, there is weak structuring within and between lineages. These data indicate the presence of gene flow from the northern into the southern parts of the southern African sub-region likely following the secondary contact. The distribution modelling predictably reveal the species' preference for rocky areas, with stable refugia through time in the northern mountain ranges, the Great Escarpment, as well as restricted areas of the Northern Cape Province and the Cape Fold Mountains of South Africa. Different microclimatic variables appear to determine the distributional range of the species. Despite strong habitat preference, the micro-habitat offered by rocky crevices and unique life history traits likely promoted the adaptability of P. capensis, resulting in the widespread distribution and persistence of the species over a long evolutionary period. Spatio-temporal comparison of the evolutionary histories of other co-distributed species across the rocky landscapes of southern Africa will improve our understanding of the regional patterns of biodiversity and local endemism.

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.

Advances using molecular data in insect systematics

The size of molecular datasets has been growing exponentially since the mid 1980s, and new technologies have now dramatically increased the slope of this increase. New datasets include genomes, transcriptomes, and hybrid capture data, producing hundreds or thousands of loci. With these datasets, we are approaching a consensus on the higher level insect phylogeny. Huge datasets can produce new challenges in interpreting branch support, and new opportunities in developing better models and more sophisticated partitioning schemes. Dating analyses are improving as we recognize the importance of careful fossil calibration selection. With thousands of genes now available, coalescent methods have come of age. Barcode libraries continue to expand, and new methods are being developed for incorporating them into phylogenies with tens of thousands of individuals.

Towards Process-based Range Modeling of Many Species

Understanding and forecasting species' geographic distributions in the face of global change is a central priority in biodiversity science. The existing view is that one must choose between correlative models for many species versus process-based models for few species. We suggest that opportunities exist to produce process-based range models for many species, by using hierarchical and inverse modeling to borrow strength across species, fill data gaps, fuse diverse data sets, and model across biological and spatial scales. We review the statistical ecology and population and range modeling literature, illustrating these modeling strategies in action. A variety of large, coordinated ecological datasets that can feed into these modeling solutions already exist, and we highlight organisms that seem ripe for the challenge.

Extracting DNA from 'jaws': high yield and quality from archived tiger shark (Galeocerdo cuvier) skeletal material

Archived specimens are highly valuable sources of DNA for retrospective genetic/genomic analysis. However, often limited effort has been made to evaluate and optimize extraction methods, which may be crucial for downstream applications. Here, we assessed and optimized the usefulness of abundant archived skeletal material from sharks as a source of DNA for temporal genomic studies. Six different methods for DNA extraction, encompassing two different commercial kits and three different protocols, were applied to material, so-called bio-swarf, from contemporary and archived jaws and vertebrae of tiger sharks (Galeocerdo cuvier). Protocols were compared for DNA yield and quality using a qPCR approach. For jaw swarf, all methods provided relatively high DNA yield and quality, while large differences in yield between protocols were observed for vertebrae. Similar results were obtained from samples of white shark (Carcharodon carcharias). Application of the optimized methods to 38 museum and private angler trophy specimens dating back to 1912 yielded sufficient DNA for downstream genomic analysis for 68% of the samples. No clear relationships between age of samples, DNA quality and quantity were observed, likely reflecting different preparation and storage methods for the trophies. Trial sequencing of DNA capture genomic libraries using 20 000 baits revealed that a significant proportion of captured sequences were derived from tiger sharks. This study demonstrates that archived shark jaws and vertebrae are potential high-yield sources of DNA for genomic-scale analysis. It also highlights that even for similar tissue types, a careful evaluation of extraction protocols can vastly improve DNA yield.

Downscaled and debiased climate simulations for North America from 21,000 years ago to 2100AD

Increasingly, ecological modellers are integrating paleodata with future projections to understand climate-driven biodiversity dynamics from the past through the current century. Climate simulations from earth system models are necessary to this effort, but must be debiased and downscaled before they can be used by ecological models. Downscaling methods and observational baselines vary among researchers, which produces confounding biases among downscaled climate simulations. We present unified datasets of debiased and downscaled climate simulations for North America from 21 ka BP to 2100AD, at 0.5° spatial resolution. Temporal resolution is decadal averages of monthly data until 1950AD, average climates for 1950–2005 AD, and monthly data from 2010 to 2100AD, with decadal averages also provided. This downscaling includes two transient paleoclimatic simulations and 12 climate models for the IPCC AR5 (CMIP5) historical (1850–2005), RCP4.5, and RCP8.5 21st-century scenarios. Climate variables include primary variables and derived bioclimatic variables. These datasets provide a common set of climate simulations suitable for seamlessly modelling the effects of past and future climate change on species distributions and diversity.

Contemporary genetic structure and postglacial demographic history of the black scorpionfish, Scorpaena porcus, in the Mediterranean and the Black Seas

Understanding the distribution of genetic diversity in the light of past demographic events linked with climatic shifts will help to forecast evolutionary trajectories of ecosystems within the current context of climate change. In this study, mitochondrial sequences and microsatellite loci were analysed using traditional population genetic approaches together with Bayesian dating and the more recent approximate Bayesian computation scenario testing. The genetic structure and demographic history of a commercial fish, the black scorpionfish, Scorpaena porcus, was investigated throughout the Mediterranean and Black Seas. The results suggest that the species recently underwent population expansions, in both seas, likely concomitant with the warming period following the Last Glacial Maximum, 20 000 years ago. A weak contemporaneous genetic differentiation was identified between the Black Sea and the Mediterranean Sea. However, the genetic diversity was similar for populations of the two seas, suggesting a high number of colonizers entered the Black Sea during the interglacial period and/or the presence of a refugial population in the Black Sea during the glacial period. Finally, within seas, an east/west genetic differentiation in the Adriatic seems to prevail, whereas the Black Sea does not show any structured spatial genetic pattern of its population. Overall, these results suggest that the Black Sea is not that isolated from the Mediterranean, and both seas revealed similar evolutionary patterns related to climate change and changes in sea level.

The Quaternary evolutionary history, potential distribution dynamics, and conservation implications for a Qinghai-Tibet Plateau endemic herbaceous perennial, Anisodus tanguticus (Solanaceae)

Various hypotheses have been proposed about the Quaternary evolutionary history of plant species on the Qinghai-Tibet Plateau (QTP), yet only a handful of studies have considered both population genetics and ecological niche context. In this study, we proposed and compared climate refugia hypotheses based on the phylogeographic pattern of Anisodus tanguticus (three plastid DNA fragments and nuclear internal transcribed spacer regions from 32 populations) and present and past species distribution models (SDMs). We detected six plastid haplotypes in two well-differentiated lineages. Although all haplotypes could be found in its western (sampling) area, only haplotypes from one lineage occurred in its eastern area. Meanwhile, most genetic variations existed between populations (F ST = 0.822). The SDMs during the last glacial maximum and last interglacial periods showed range fragmentation in the western area and significant range contraction in the eastern area, respectively, in comparison with current potential distribution. This species may have undergone intraspecific divergence during the early Quaternary, which may have been caused by survival in different refugia during the earliest known glacial in the QTP, rather than geological isolation due to orogenesis events. Subsequently, climate oscillations during the Quaternary resulted in a dynamic distribution range for this species as well as the distribution pattern of its plastid haplotypes and nuclear genotypes. The interglacial periods may have had a greater effect on A. tanguticus than the glacial periods. Most importantly, neither genetic data nor SDM alone can fully reveal the climate refugia history of this species. We also discuss the conservation implications for this important Tibetan folk medicine plant in light of these findings and SDMs under future climate models. Together, our results underline the necessity to combine phylogeographic and SDM approaches in future investigations of the Quaternary evolutionary history of species in topographically complex areas, such as the QTP.

Where to Dig for Fossils: Combining Climate-Envelope, Taphonomy and Discovery Models

Fossils represent invaluable data to reconstruct the past history of life, yet fossil-rich sites are often rare and difficult to find. The traditional fossil-hunting approach focuses on small areas and has not yet taken advantage of modelling techniques commonly used in ecology to account for an organism's past distributions. We propose a new method to assist finding fossils at continental scales based on modelling the past distribution of species, the geological suitability of fossil preservation and the likelihood of fossil discovery in the field, and apply it to several genera of Australian megafauna that went extinct in the Late Quaternary. Our models predicted higher fossil potentials for independent sites than for randomly selected locations (mean Kolmogorov-Smirnov statistic = 0.66). We demonstrate the utility of accounting for the distribution history of fossil taxa when trying to find the most suitable areas to look for fossils. For some genera, the probability of finding fossils based on simple climate-envelope models was higher than the probability based on models incorporating current conditions associated with fossil preservation and discovery as predictors. However, combining the outputs from climate-envelope, preservation, and discovery models resulted in the most accurate predictions of potential fossil sites at a continental scale. We proposed potential areas to discover new fossils of Diprotodon, Zygomaturus, Protemnodon, Thylacoleo, and Genyornis, and provide guidelines on how to apply our approach to assist fossil hunting in other continents and geological settings.

Tackling intraspecific genetic structure in distribution models better reflects species geographical range

Genetic diversity provides insight into heterogeneous demographic and adaptive history across organisms' distribution ranges. For this reason, decomposing single species into genetic units may represent a powerful tool to better understand biogeographical patterns as well as improve predictions of the effects of GCC (global climate change) on biodiversity loss. Using 279 georeferenced Iberian accessions, we used classes of three intraspecific genetic units of the annual plant Arabidopsis thaliana obtained from the genetic analyses of nuclear SNPs (single nucleotide polymorphisms), chloroplast SNPs, and the vernalization requirement for flowering. We used SDM (species distribution models), including climate, vegetation, and soil data, at the whole-species and genetic-unit levels. We compared model outputs for present environmental conditions and with a particularly severe GCC scenario. SDM accuracy was high for genetic units with smaller distribution ranges. Kernel density plots identified the environmental variables underpinning potential distribution ranges of genetic units. Combinations of environmental variables accounted for potential distribution ranges of genetic units, which shrank dramatically with GCC at almost all levels. Only two genetic clusters increased their potential distribution ranges with GCC. The application of SDM to intraspecific genetic units provides a detailed picture on the biogeographical patterns of distinct genetic groups based on different genetic criteria. Our approach also allowed us to pinpoint the genetic changes, in terms of genetic background and physiological requirements for flowering, that Iberian A. thaliana may experience with a GCC scenario applying SDM to intraspecific genetic units.

Predicting the genetic consequences of future climate change: The power of coupling spatial demography, the coalescent, and historical landscape changes

PREMISE OF THE STUDY

Climate change is a widely accepted threat to biodiversity. Species distribution models (SDMs) are used to forecast whether and how species distributions may track these changes. Yet, SDMs generally fail to account for genetic and demographic processes, limiting population-level inferences. We still do not understand how predicted environmental shifts will impact the spatial distribution of genetic diversity within taxa.

METHODS

We propose a novel method that predicts spatially explicit genetic and demographic landscapes of populations under future climatic conditions. We use carefully parameterized SDMs as estimates of the spatial distribution of suitable habitats and landscape dispersal permeability under present-day, past, and future conditions. We use empirical genetic data and approximate Bayesian computation to estimate unknown demographic parameters. Finally, we employ these parameters to simulate realistic and complex models of responses to future environmental shifts. We contrast parameterized models under current and future landscapes to quantify the expected magnitude of change.

KEY RESULTS

We implement this framework on neutral genetic data available from Penstemon deustus. Our results predict that future climate change will result in geographically widespread declines in genetic diversity in this species. The extent of reduction will heavily depend on the continuity of population networks and deme sizes.

CONCLUSIONS

To our knowledge, this is the first study to provide spatially explicit predictions of within-species genetic diversity using climatic, demographic, and genetic data. Our approach accounts for climatic, geographic, and biological complexity. This framework is promising for understanding evolutionary consequences of climate change, and guiding conservation planning.

Mesocosms Reveal Ecological Surprises from Climate Change

Understanding, predicting, and mitigating the impacts of climate change on biodiversity poses one of the most crucial challenges this century. Currently, we know more about how future climates are likely to shift across the globe than about how species will respond to these changes. Two recent studies show how mesocosm experiments can hasten understanding of the ecological consequences of climate change on species' extinction risk, community structure, and ecosystem functions. Using a large-scale terrestrial warming experiment, Bestion et al. provide the first direct evidence that future global warming can increase extinction risk for temperate ectotherms. Using aquatic mesocosms, Yvon-Durocher et al. show that human-induced climate change could, in some cases, actually enhance the diversity of local communities, increasing productivity. Blending these theoretical and empirical results with computational models will improve forecasts of biodiversity loss and altered ecosystem processes due to climate change.