Species distribution models contribute to determine the effect of climate and interspecific interactions in moving hybrid zones

Climate is a major factor delimiting species’ distributions. However, biotic interactions may also be prominent in shaping geographical ranges, especially for parapatric species forming hybrid zones. Determining the relative effect of each factor and their interaction of the contact zone location has been difficult due to the lack of broad scale environmental data. Recent developments in species distribution modelling (SDM) now allow disentangling the relative contributions of climate and species’ interactions in hybrid zones and their responses to future climate change. We investigated the moving hybrid zone between the breeding ranges of two parapatric passerines in Europe. We conducted SDMs representing the climatic conditions during the breeding season. Our results show a large mismatch between the realized and potential distributions of the two species, suggesting that interspecific interactions, not climate, account for the present location of the contact zone. The SDM scenarios show that the southerly distributed species, Hippolais polyglotta, might lose large parts of its southern distribution under climate change, but a similar gain of novel habitat along the hybrid zone seems unlikely, because interactions with the other species (H. icterina) constrain its range expansion. Thus, whenever biotic interactions limit range expansion, species may become ‘trapped’ if range loss due to climate change is faster than the movement of the contact zone. An increasing number of moving hybrid zones are being reported, but the proximate causes of movement often remain unclear. In a global context of climate change, we call for more interest in their interactions with climate change.     

Projected climate changes threaten ancient refugia of kelp forests in the North Atlantic

Intraspecific genetic variability is critical for species adaptation and evolution and yet it is generally overlooked in projections of the biological consequences of climate change. We ask whether ongoing climate changes can cause the loss of important gene pools from North Atlantic relict kelp forests that persisted over glacial–interglacial cycles. We use ecological niche modelling to predict genetic diversity hotspots for eight species of large brown algae with different thermal tolerances (Arctic to warm temperate), estimated as regions of persistence throughout the Last Glacial Maximum (20,000 YBP), the warmer Mid‐Holocene (6,000 YBP), and the present. Changes in the genetic diversity within ancient refugia were projected for the future (year 2100) under two contrasting climate change scenarios (RCP2.6 and RCP8.5). Models predicted distributions that matched empirical distributions in cross‐validation, and identified distinct refugia at the low latitude ranges, which largely coincide among species with similar ecological niches. Transferred models into the future projected polewards expansions and substantial range losses in lower latitudes, where richer gene pools are expected (in Nova Scotia and Iberia for cold affinity species and Gibraltar, Alboran, and Morocco for warm‐temperate species). These effects were projected for both scenarios but were intensified under the extreme RCP8.5 scenario, with the complete borealization (circum‐Arctic colonization) of kelp forests, the redistribution of the biogeographical transitional zones of the North Atlantic, and the erosion of global gene pools across all species. As the geographic distribution of genetic variability is unknown for most marine species, our results represent a baseline for identification of locations potentially rich in unique phylogeographic lineages that are also climatic relics in threat of disappearing.     

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

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

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.     

Geography best explains global patterns of genetic diversity and postglacial co‐expansion in marine turtles

For many species, climate oscillations drove cycles of population contraction during cool glacial periods followed by expansion during interglacials. Some groups, however, show evidence of uniform and synchronous expansion, while others display differences in the timing and extent of demographic change. We compared demographic histories inferred from genetic data across marine turtle species to identify responses to postglacial warming shared across taxa and to examine drivers of past demographic change at the global scale. Using coalescent simulations and approximate Bayesian computation (ABC), we estimated demographic parameters, including the likelihood of past population expansion, from a mitochondrial data set encompassing 23 previously identified lineages from all seven marine turtle species. For lineages with a high posterior probability of expansion, we conducted a hierarchical ABC analysis to estimate the proportion of lineages expanding synchronously and the timing of synchronous expansion. We used Bayesian model averaging to identify variables associated with expansion and genetic diversity. Approximately 60% of extant marine turtle lineages showed evidence of expansion, with the rest mainly exhibiting patterns of genetic diversity most consistent with population stability. For lineages showing expansion, there was a strong signal of synchronous expansion after the Last Glacial Maximum. Expansion and genetic diversity were best explained by ocean basin and the degree of endemism for a given lineage. Geographic differences in sensitivity to climate change have implications for prioritizing conservation actions in marine turtles as well as for identifying areas of past demographic stability and potential resilience to future climate change for broadly distributed taxa.     

Genetic consequences of climate change for northern plants

Climate change will lead to loss of range for many species, and thus to loss of genetic diversity crucial for their long-term persistence. We analysed range-wide genetic diversity (amplified fragment length polymorphisms) in 9581 samples from 1200 populations of 27 northern plant species, to assess genetic consequences of range reduction and potential association with species traits. We used species distribution modelling (SDM, eight techniques, two global circulation models and two emission scenarios) to predict loss of range and genetic diversity by 2080. Loss of genetic diversity varied considerably among species, and this variation could be explained by dispersal adaptation (up to 57%) and by genetic differentiation among populations (F(ST); up to 61%). Herbs lacking adaptations for long-distance dispersal were estimated to lose genetic diversity at higher rate than dwarf shrubs adapted to long-distance dispersal. The expected range reduction in these 27 northern species was larger than reported for temperate plants, and all were predicted to lose genetic diversity according to at least one scenario. SDM combined with F(ST) estimates and/or with species trait information thus allows the prediction of species' vulnerability to climate change, aiding rational prioritization of conservation efforts.     

Geological habitat template overrides late Quaternary climate change as a determinant of range dynamics and phylogeography in some habitat‐specialist water beetles

Aim We investigated the roles of lithology and climate in constraining the ranges of four co‐distributed species of Iberian saline‐habitat specialist water beetles (Ochthebius glaber, Ochthebius notabilis, Enochrus falcarius and Nebrioporus baeticus) across the late Quaternary and in shaping their geographical genetic structure. The aim was to improve our understanding of the effects of past climate changes on the biota of arid Mediterranean environments and of the relative importance of history and landscape on phylogeographical patterns. Location Iberian Peninsula, Mediterranean. Methods We combined species distribution modelling (SDM) and comparative phylogeography. We used a multi‐model inference and model‐averaging approach both for assessment of range determinants (climate and lithology) and for provision of spatially explicit estimates of the species current and Last Glacial Maximum (LGM) potential ranges. Potential LGM distributions were then contrasted with the phylogeographical and population expansion patterns as assessed using mitochondrial DNA sequence data. We also evaluated the relative importance of geographical distance, habitat resistance and historical isolation for genetic structure in a causal modelling framework. Results Lithology poses a strong constraint on the distribution of Iberian saline‐habitat specialist water beetles, with a variable, but generally moderate, additional influence by climate. The degree to which potential LGM distributions were reduced and fragmented decreased with increasing importance of lithology. These SDM‐based suitability predictions were mostly congruent with phylogeographical and population genetic patterns across the study species, with stronger geographical structure in the genetic diversity of the more temperature‐sensitive species (O. glaber and E. falcarius). Furthermore, while historical isolation was the only factor explaining genetic structure in the more temperature‐sensitive species, lithology‐controlled landscape configuration also played an important role for those species with more lithology‐determined ranges (O. notabilis and N. baeticus). Main conclusions Our data show that lithology is an important constraint on the distribution and range dynamics of endemic Iberian saline‐habitat water beetles, in interaction with climate and long‐term climate change, and overrides the latter in importance for some species. Hence, geological landscape structure and long‐term history may codetermine the overall range and the distribution of genetic lineages in endemic species with specialized edaphic requirements.     

Ice‐age persistence and genetic isolation of the disjunct distribution of larch in Alaska

Larix laricina (eastern larch, tamarack) is a transcontinental North American conifer with a prominent disjunction in the Yukon isolating the Alaskan distribution from the rest of its range. We investigate whether in situ persistence during the last glacial maximum (LGM) or long‐distance postglacial migration from south of the ice sheets resulted in the modern‐day Alaskan distribution. We analyzed variation in three chloroplast DNA regions of 840 trees from a total of 69 populations (24 new sampling sites situated on both sides of the Yukon range disjunction pooled with 45 populations from a published source) and conducted ensemble species distribution modeling (SDM) throughout Canada and United States to hindcast the potential range of L. laricina during the LGM. We uncovered the genetic signature of a long‐term isolation of larch populations in Alaska, identifying three endemic chlorotypes and low levels of genetic diversity. Range‐wide analysis across North America revealed the presence of a distinct Alaskan lineage. Postglacial gene flow across the Yukon divide was unidirectional, from Alaska toward previously glaciated Canadian regions, and with no evidence of immigration into Alaska. Hindcast SDM indicates one of the broadest areas of past climate suitability for L. laricina existed in central Alaska, suggesting possible in situ persistence of larch in Alaska during the LGM. Our results provide the first unambiguous evidence for the long‐term isolation of L. laricina in Alaska that extends beyond the last glacial period and into the present interglacial period. The lack of gene flow into Alaska along with the overall probability of larch occurrence in Alaska being currently lower than during the LGM suggests that modern‐day Alaskan larch populations are isolated climate relicts of broader glacial distributions, and so are particularly vulnerable to current warming trends.     

Past and future demographic dynamics of alpine species: limited genetic consequences despite dramatic range contraction in a plant from the Spanish Sierra Nevada

Anthropogenic global climate change is expected to cause severe range contractions among alpine plants. Alpine areas in the Mediterranean region are of special concern because of the high abundance of endemic species with narrow ranges. This study combined species distribution models, population structure analyses and Bayesian skyline plots to trace the past and future distribution and diversity of Linaria glacialis, an endangered narrow endemic species that inhabits summits of Sierra Nevada (Spain). The results showed that: (i) the habitat of this alpine‐Mediterranean species in Sierra Nevada suffered little changes during glacial and interglacial stages of late Quaternary; (ii) climatic oscillations in the last millennium (Medieval Warm Period and Little Ice Age) moderately affected the demographic trends of L. glacialis; (iii) future warming conditions will cause severe range contractions; and (iv) genetic diversity will not diminish at the same pace as the distribution range. As a consequence of the low population structure of this species, genetic impoverishment in the alpine zones of Sierra Nevada should be limited during range contraction. We conclude that maintenance of large effective population sizes via high mutation rates and high levels of gene flow may promote the resilience of alpine plant species when confronted with global warming.     

Are Plant Species Able to Keep Pace with the Rapidly Changing Climate?

Future climate change is predicted to advance faster than the postglacial warming. Migration may therefore become a key driver for future development of biodiversity and ecosystem functioning. For 140 European plant species we computed past range shifts since the last glacial maximum and future range shifts for a variety of Intergovernmental Panel on Climate Change (IPCC) scenarios and global circulation models (GCMs). Range shift rates were estimated by means of species distribution modelling (SDM). With process-based seed dispersal models we estimated species-specific migration rates for 27 dispersal modes addressing dispersal by wind (anemochory) for different wind conditions, as well as dispersal by mammals (dispersal on animal''s coat – epizoochory and dispersal by animals after feeding and digestion – endozoochory) considering different animal species. Our process-based modelled migration rates generally exceeded the postglacial range shift rates indicating that the process-based models we used are capable of predicting migration rates that are in accordance with realized past migration. For most of the considered species, the modelled migration rates were considerably lower than the expected future climate change induced range shift rates. This implies that most plant species will not entirely be able to follow future climate-change-induced range shifts due to dispersal limitation. Animals with large day- and home-ranges are highly important for achieving high migration rates for many plant species, whereas anemochory is relevant for only few species.     

Investigating past range dynamics for a weed of cultivation,Silene vulgaris

Since the last glacial maximum (LGM), many plant and animal taxa have expanded their ranges by migration from glacial refugia. Weeds of cultivation may have followed this trend or spread globally following the expansion of agriculture or ruderal habitats associated with human‐mediated disturbance. We tested whether the range expansion of the weed Silene vulgaris across Europe fit the classical model of postglacial expansion from southern refugia, or followed known routes of the expansion of human agricultural practices. We used species distribution modeling to predict spatial patterns of postglacial expansion and contrasted these with the patterns of human agricultural expansion. A population genetic analysis using microsatellite loci was then used to test which scenario was better supported by spatial patterns of genetic diversity and structure. Genetic diversity was highest in southern Europe and declined with increasing latitude. Locations of ancestral demes from genetic cluster analysis were consistent with areas of predicted refugia. Species distribution models showed the most suitable habitat in the LGM on the southern coasts of Europe. These results support the typical postglacial northward colonization from southern refugia while refuting the east‐to‐west agricultural spread as the main mode of expansion for S. vulgaris. We know that S. vulgaris has recently colonized many regions (including North America and other continents) through human‐mediated dispersal, but there is no evidence for a direct link between the Neolithic expansion of agriculture and current patterns of genetic diversity of S. vulgaris in Europe. Therefore, the history of range expansion of S. vulgaris likely began with postglacial expansion after the LGM, followed by more recent global dispersal by humans.     

Incorporating intraspecific variation into species distribution models improves distribution predictions,but cannot predict species traits for a wide-spread plant species

The most common approach to predicting how species ranges and ecological functions will shift with climate change is to construct correlative species distribution models (SDMs). These models use a species’ climatic distribution to determine currently suitable areas for the species and project its potential distribution under future climate scenarios. A core, rarely tested, assumption of SDMs is that all populations will respond equivalently to climate. Few studies have examined this assumption, and those that have rarely dissect the reasons for intraspecific differences. Focusing on the arctic-alpine cushion plant Silene acaulis, we compared predictive accuracy from SDMs constructed using the species’ full global distribution with composite predictions from separate SDMs constructed using subpopulations defined either by genetic or habitat differences. This is one of the first studies to compare multiple ways of constructing intraspecific-level SDMs with a species-level SDM. We also examine the contested relationship between relative probability of occurrence and species performance or ecological function, testing if SDM output can predict individual performance (plant size) and biotic interactions (facilitation). We found that both genetic- and habitat-informed SDMs are considerably more accurate than a species-level SDM, and that the genetic model substantially differs from and outperforms the habitat model. While SDMs have been used to infer population performance and possibly even biotic interactions, in our system these relationships were extremely weak. Our results indicate that individual subpopulations may respond differently to climate, although we discuss and explore several alternative explanations for the superior performance of intraspecific-level SDMs. We emphasize the need to carefully examine how to best define intraspecific-level SDMs as well as how potential genetic, environmental, or sampling variation within species ranges can critically affect SDM predictions. We urge caution in inferring population performance or biotic interactions from SDM predictions, as these often-assumed relationships are not supported in our study.     

Chloroplast DNA variation and geographical structure of the Aristolochia kaempferi group (Aristolochiaceae)

The present study documents cpDNA variation in the Aristolochia kaempferi group (Aristolochiaceae), which consists of one Chinese and all Japanese and Taiwanese species of the subgenus Siphisia. In a phylogenetic analysis based on the nucleotide sequences of the matK gene, and the atpB-rbcL and trnS-trnG intergenic spacer regions, 38 haplotypes were recognized in the A. kaempferi group and as many as 24 within A. kaempferi. This is the most haplotypes reported for a single species to date. Although six highly significant major clades were identified in the phylogenetic analysis, they were not congruent with previous classifications. This might be attributed to the specific speciation process, such as convergent evolution, incomplete lineage sorting, and/or reticulate evolution. The six major clades had a clear geographical distribution pattern and were significantly associated with geographical distribution of haplotypes in a nested clade analysis and AMOVA. The results allow us to deduce a scenario in which multiple contractions and expansions of the geographical ranges brought about by Quaternary climatic oscillations affected the patterns of genetic diversity. The present geographic patterns of haplotype distribution within the A. kaempferi group can be explained by the last postglacial range expansion from different refugia, and the boundaries may be suture zones.     

Why is the choice of future climate scenarios for species distribution modelling important?

Species distribution models (SDMs) are common tools for assessing the potential impact of climate change on species ranges. Uncertainty in SDM output occurs due to differences among alternate models, species characteristics and scenarios of future climate. While considerable effort is being devoted to identifying and quantifying the first two sources of variation, a greater understanding of climate scenarios and how they affect SDM output is also needed. Climate models are complex tools: variability occurs among alternate simulations, and no single 'best' model exists. The selection of climate scenarios for impacts assessments should not be undertaken arbitrarily - strengths and weakness of different climate models should be considered. In this paper, we provide bioclimatic modellers with an overview of emissions scenarios and climate models, discuss uncertainty surrounding projections of future climate and suggest steps that can be taken to reduce and communicate climate scenario-related uncertainty in assessments of future species responses to climate change.     

Disregarding the edaphic dimension in species distribution models leads to the omission of crucial spatial information under climate change: the case of Quercus pubescens in France

Species distribution modelling is an easy, persuasive and useful tool for anticipating species distribution shifts under global change. Numerous studies have used only climate variables to predict future potential species range shifts and have omitted environmental factors important for determining species distribution. Here, we assessed the importance of the edaphic dimension in the niche‐space definition of Quercus pubescens and in future spatial projections under global change over the metropolitan French forest territory. We fitted two species distribution models (SDM) based on presence/absence data (111 013 plots), one calibrated from climate variables only (mean temperature of January and climatic water balance of July) and the other one from both climate and edaphic (soil pH inferred from plants) variables. Future predictions were conducted under two climate scenarios (PCM B2 and HadCM3 A2) and based on 100 simulations using a cellular automaton that accounted for seed dispersal distance, landscape barriers preventing migration and unsuitable land cover. Adding the edaphic dimension to the climate‐only SDM substantially improved the niche‐space definition of Q. pubescens, highlighting an increase in species tolerance in confronting climate constraints as the soil pH increased. Future predictions over the 21st century showed that disregarding the edaphic dimension in SDM led to an overestimation of the potential distribution area, an underestimation of the spatial fragmentation of this area, and prevented the identification of local refugia, leading to an underestimation of the northward shift capacity of Q. pubescens and its persistence in its current distribution area. Spatial discrepancies between climate‐only and climate‐plus‐edaphic models are strengthened when seed dispersal and forest fragmentation are accounted for in predicting a future species distribution area. These discrepancies highlight some imprecision in spatial predictions of potential distribution area of species under climate change scenarios and possibly wrong conclusions for conservation and management perspectives when climate‐only models are used.     

Phylogeography and paleodistribution models of a widespread birch (Betula platyphylla Suk.) across East Asia: Multiple refugia,multidirectional expansion,and heterogeneous genetic pattern

Widespread tree species cover large geographical areas and play important roles in various vegetation types. Understanding how these species responded to historical climatic changes is important for understanding community assembly mechanisms with evolutionary and conservation implications. However, the location of refugial areas and postglacial history of widespread trees in East Asia remain poorly known. We combined microsatellite data (63 populations, 1756 individuals) and ecological niche modeling to examine the range‐wide population diversity, genetic structure, and historical demography of a pioneer tree species, Asian white birch (Betula platyphylla Suk.) across East Asia. We found a north‐to‐south trend of declining genetic diversity and five clusters, corresponding to geographical regions. Different clusters were inferred to have diverged through Pleistocene climatic oscillations and have different expansion routes, leading to genetic admixture in some populations. Ecological niche models indicated that the distribution of B. platyphylla during the last glacial maximum still had a large latitude span with slight shifts toward southeast, and northern populations had more variable distribution ranges than those in the south during later climatic oscillations. Our results reflect the relatively stable distribution through the last glacial–interglacial cycles and recent multidirectional expansion of B. platyphylla, providing new hypotheses for the response pattern of widespread tree species to climate change. The gradual genetic pattern from northeast to southwest and alternative distribution dynamics possibly resulted from environmental differences caused by latitude and topographic heterogeneity.     

Northern glacial refugia for the pygmy shrew Sorex minutus in Europe revealed by phylogeographic analyses and species distribution modelling

The southern European peninsulas (Iberian, Italian and Balkan) are traditionally recognized as glacial refugia from where many species colonized central and northern Europe after the Last Glacial Maximum (LGM). However, evidence that some species had more northerly refugia is accumulating from phylogeographic, palaeontological and palynological studies, and more recently from species distribution modelling (SDM), but further studies are needed to test the idea of northern refugia in Europe. Here, we take a rarely implemented multidisciplinary approach to assess if the pygmy shrew Sorex minutus, a widespread Eurasian mammal species, had northern refugia during the LGM, and if these influenced its postglacial geographic distribution. First, we evaluated the phylogeographic and population expansion patterns using mtDNA sequence data from 123 pygmy shrews. Then, we used SDM to predict present and past (LGM) potential distributions using two different training data sets, two different algorithms (Maxent and GARP) and climate reconstructions for the LGM with two different general circulation models. An LGM distribution in the southern peninsulas was predicted by the SDM approaches, in line with the occurrence of lineages of S. minutus in these areas. The phylogeographic analyses also indicated a widespread and strictly northern‐central European lineage, not derived from southern peninsulas, and with a postglacial population expansion signature. This was consistent with the SDM predictions of suitable LGM conditions for S. minutus occurring across central and eastern Europe, from unglaciated parts of the British Isles to much of the eastern European Plain. Hence, S. minutus likely persisted in parts of central and eastern Europe during the LGM, from where it colonized other northern areas during the late‐glacial and postglacial periods. Our results provide new insights into the glacial and postglacial colonization history of the European mammal fauna, notably supporting glacial refugia further north than traditionally recognized.     

Rethinking long‐term vegetation dynamics: multiple glacial refugia and local expansion of a species complex

Evidence is accumulating that some arcto‐boreal plant taxa persisted through the last glacial maximum (LGM) in Alaska and adjacent Canada. However, the spatial patterns of glacial persistence and associated postglacial colonization remain largely unknown. In this study, we investigated the LGM refugia of an alder (Alnus) species complex (n = 3 taxa) and assess the spatiotemporal dynamics of Alnus in this vast region. Specifically, we conducted high‐throughput DNA sequencing (ddRADseq) on Alnus foliar samples collected from a dense population network to investigate patterns of genetic structure and infer the presence of glacial lineages. Species distribution modeling (SDM) was used to investigate the probability and possible locations of glacial persistence. These analyses were integrated and then compared with fossil pollen data to identify the locations of refugial populations and spatial patterns of postglacial colonization. Our genetic analyses revealed two glacial lineages with separate geographic origins for each Alnus taxon, suggesting that the genus persisted in multiple LGM refugia. Non‐overlapping hindcast distributions based on SDMs further support the presence of multiple, spatially distinct refugia. These ddRADseq and SDM results, in conjunction with reassessment of fossil pollen records, suggest that Alnus expanded from several population nuclei that existed during the LGM and coalesced during the Holocene to form its present range. These results challenge the unidirectional model for postglacial vegetation expansion, implying that climate buffering associated with landscape heterogeneity and adaptation to millennial‐scale environmental variability played important roles in driving late‐Quaternary population dynamics.     

Genetic signals of demographic expansion in Downy Woodpecker (Picoides pubescens) after the last North American glacial maximum

The glacial cycles of the Pleistocene have been recognized as important, large-scale historical processes that strongly influenced the demographic patterns and genetic structure of many species. Here we present evidence of a postglacial expansion for the Downy Woodpecker (Picoides pubescens), a common member of the forest bird communities in North America with a continental distribution. DNA sequences from the mitochondrial tRNA-Lys, and ATPase 6 and 8 genes, and microsatellite data from seven variable loci were combined with a species distribution model (SDM) to infer possible historical scenarios for this species after the last glacial maximum. Analyses of Downy Woodpeckers from 23 geographic areas suggested little differentiation, shallow genealogical relationships, and limited population structure across the species' range. Microsatellites, which have higher resolution and are able to detect recent differences, revealed two geographic groups where populations along the eastern edge of the Rocky Mountains (Montana, Utah, Colorado, and southern Alberta) were genetically isolated from the rest of the sampled populations. Mitochondrial DNA, an important marker to detect historical patterns, recovered only one group. However, populations in Idaho and southeast BC contained high haplotype diversity and, in general were characterized by the absence of the most common mtDNA haplotype. The SDM suggested several areas in the southern US as containing suitable Downy Woodpecker habitat during the LGM. The lack of considerable geographic structure and the starburst haplotype network, combined with several population genetic tests, suggest a scenario of demographic expansion during the last part of Pleistocene and early Holocene.