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.     

Phylogeography and ecological niche modelling of the New Zealand stick insect Clitarchus hookeri (White) support survival in multiple coastal refugia

Aim Increasing our understanding of the effects of the Last Glacial Maximum (LGM) and determining the location of refugia requires studies on widely distributed species with dense sampling of populations. We have reconstructed the biogeographic history of Clitarchus hookeri (White), a widespread species of New Zealand stick insect that exhibits geographic parthenogenesis, using phylogeographic analysis and ecological niche modelling. Location New Zealand. Methods We used DNA sequence data from the mitochondrial cytochrome c oxidase subunit I gene to reconstruct phylogenetic relationships among haplotypes from C. hookeri and two undescribed Clitarchus species. We also used distribution data from our own field surveys and museum records to reconstruct the geographic distribution of C. hookeri during the present and the LGM, using ecological niche modelling. Results The ecological niche models showed that the geographic distribution of C. hookeri has expanded dramatically since the LGM. Our model predicted large areas of suitable LGM habitat in upper North Island, and small patches along the east coast of South Island. The phylogeographic analysis shows that populations in the northern half of North Island contain much higher levels of genetic variation than those from southern North Island and South Island, and is congruent with the ecological niche model. The distribution of bisexual populations is also non-random, with males completely absent from South Island and very rare in southern North Island. Main conclusions During the LGM C. hookeri was most likely restricted to several refugia in upper North Island and one or more smaller refugia along the east coast of South Island. The unisexual populations predominate in post-glacial landscapes and are clearly favoured in the recolonization of such areas. Our study exemplifies the utility of integrating ecological niche modelling and phylogeographic analysis.     

Multiple glacial refugia for cool‐temperate deciduous trees in northern East Asia: the Mongolian oak as a case study

In East Asia, temperate forests are predicted to have retracted southward to c. 30° N during the last glacial maximum (LGM) based on fossil pollen data, whereas phylogeographic studies have often suggested glacial in situ survival of cool‐temperate deciduous trees in their modern northern ranges. Here we report a study of the genetic diversity and structure of 29 natural Mongolian oak (Quercus mongolica) populations using 19 nuclear simple sequence repeat (nSSR) loci and four chloroplast DNA fragments. Bayesian clustering analysis with nSSRs revealed five groups, which were inferred by approximate Bayesian computation (ABC) to have diverged in multiple refugia through multiple glacial–interglacial cycles. Analysis of chloroplast DNA variation revealed four lineages that were largely but incompletely geographically disjunct. Ecological niche modelling (ENMs) indicated a southward range shift of the oak's distribution at the LGM, although high suitability scores were also evident in the Changbai Mts. (Northeast China), the Korean Peninsula, areas surrounding the Bohai Sea, and along the coast of the Russian Far East. In addition, endemic chloroplast DNA haplotypes and nuclear lineages occurred in high‐latitude northern areas where the ENM predicted no suitable habitat. The combined evidence from nuclear and chloroplast DNA, and the results of the ENM clearly demonstrate that multiple northern refugia, including cryptic ones, were maintained across the current distributional range of the Mongolian oak during the LGM or earlier glacial periods. Though spatially limited, postglacial expansions from these refugia have led to a pattern of decreased genetic diversity with increasing latitude.     

Contrasting genetic patterns between two coexisting Eleutherococcus species in northern China

Climate oscillations are the key factors to understand the patterns in modern biodiversity. East Asia harbors the most diverse temperate flora, largely because an extensive terrestrial ice cap was absent during repeated Pleistocene glaciation–interglacial cycles. Comparing the demographic histories of species that are codistributed and are close relatives may provide insight into how the process of climate change influences species ranges. In this study, we compared the spatial genetic structure and demographic histories of two coexisting Eleutherococcus species, Eleutherococcus senticosus and E. sessiliflorus. Both species are distributed in northern China, regions that are generally considered to be sensitive to climatic fluctuations. These regions once hosted temperate forest, but this temperate forest was replaced by tundra and taiga forest during the Last Glacial Maximum (LGM), according to pollen records. Using three chloroplast DNA fragments, we assessed the genetic structure of 20 and 9 natural populations of E. senticosus and E. sessiliflorus, respectively. Extremely contrasting genetic patterns were found between the two species; E. sessiliflorus had little genetic variation, whereas E. senticosus had considerably higher levels of genetic variation (15 haplotypes). We speculated that a recent severe bottleneck may have resulted in the extremely low genetic diversity in E. sessiliflorus. In E. senticosus, populations in Northeast China (NEC) harbored all of the haplotypes found in this species and included private haplotypes. The populations in NEC had higher levels of genetic diversity than did those from North China (NC). Therefore, we suggest that both the NC and NEC regions can sustain LGM refugia and that lineage admixture from multiple refugia took place after the LGM elevated the local genetic diversity in NEC. In NEC, multiple genetic hot spots were found in the Changbai Mountains and the Xiaoxing'an Range, which implied that multiple locations in NEC may sustain LGM refugia, even in the Xiaoxing'an Range.     

River islands,refugia and genetic structuring in the endemic brown frog Rana kukunoris (Anura,Ranidae) of the Qinghai‐Tibetan Plateau

Frequently, Pleistocene climatic cycling has been found to be the diver of genetic structuring in populations, even in areas that did not have continental ice sheets, such as on the Qinghai‐Tibetan Plateau (QTP). Typically, species distributed on the plateau have been hypothesized to re‐treat to south‐eastern refugia, especially during the Last Glacial Maximum (LGM). We evaluated sequence variation in the mitochondrial DNA gene Cytb and the nuclear DNA gene RAG‐1 in Rana kukunoris, a species endemic to the QTP. Two major lineages, N and S, were identified, and lineage N was further subdivided into N1 and N2. The geographical distribution and genealogical divergences supported the hypothesis of multiple refugia. However, major lineages and sublineages diverged prior to the LGM. Demographical expansion was detected only in lineage S and sublineage N2. Sublineage N1 might have survived several glacial cycles in situ and did not expand after the LGM because of the absence of suitable habitat; it survived in river islands. Genetic analysis and environment modelling suggested that the north‐eastern edge of QTP contained a major refugium for R. kukunoris. From here, lineage S dispersed southwards after the LGM. Two microrefugia in northern Qilian Mountains greatly contributed to current level of intraspecific genetic diversity. These results were found to have important implications for the habitat conservation in Northwest China.     

Combining ecological niche modelling and morphology to assess the range‐wide population genetic structure of bobcats (Lynx rufus)

Despite a broad distribution, general habitat requirements, and a large dispersal potential, bobcats (Lynx rufus) exhibit a genetic division that longitudinally transects central North America. We investigated (1) whether the climate of the Last Glacial Maximum (LGM; 21 kya) isolated bobcats into refugia and also whether the current climate influences gene flow between the segregate populations and (2) whether the geographical patterns in cranial morphology reflect population identity. We created ecological niche models (ENMs) to evaluate climatic suitability and to estimate distributions of the disparate populations under both historical (LGM) and contemporary conditions. We used two‐dimensional geometric morphometric methods to evaluate variations in the cranium and mandible. These variations were then regressed across geographical variables to assess morphological differences throughout the range of the bobcat. ENMs projected onto LGM climate provided evidence of refugia during the LGM via increased suitability in the north‐west and south‐east portions of this species' range. Contemporarily, our models suggest that the Great Plains may be restricting bobcat migration and gene flow, effectively maintaining disparate populations. Morphological analyses identified a significant linear trend in shape variation across latitudinal and longitudinal gradients rather than distinct morphological divergence between lineages. Similar shape variations, however, did converge in approximate locations of assumed refugia. The findings of the present study provide a robust assessment of the biogeographical considerations for the population genetic structure of bobcats.     

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.     

Identifying climate refugia and its potential impact on Tibetan brown bear (Ursus arctos pruinosus) in Sanjiangyuan National Park,China

Climate change has direct impacts on wildlife and future biodiversity protection efforts. Vulnerability assessment and habitat connectivity analyses are necessary for drafting effective conservation strategies for threatened species such as the Tibetan brown bear (Ursus arctos pruinosus). We used the maximum entropy (MaxEnt) model to assess the current (1950–2000) and future (2041–2060) habitat suitability by combining bioclimatic and environmental variables, and identified potential climate refugia for Tibetan brown bears in Sanjiangyuan National Park, China. Next, we selected Circuit model to simulate potential migration paths based on current and future climatically suitable habitat. Results indicate a total area of potential suitable habitat under the current climate scenario of approximately 31,649.46 km², of which 28,778.29 km² would be unsuitable by the 2050s. Potentially suitable habitat under the future climate scenario was projected to cover an area of 23,738.6 km². Climate refugia occupied 2,871.17 km², primarily in the midwestern and northeastern regions of Yangtze River Zone, as well as the northern region of Yellow River Zone. The altitude of climate refugia ranged from 4,307 to 5,524 m, with 52.93% lying at altitudes between 4,300 and 4,600 m. Refugia were mainly distributed on bare rock, alpine steppe, and alpine meadow. Corridors linking areas of potentially suitable brown bear habitat and a substantial portion of paths with low‐resistance value were distributed in climate refugia. We recommend various actions to ameliorate the impact of climate change on brown bears, such as protecting climatically suitable habitat, establishing habitat corridors, restructuring conservation areas, and strengthening monitoring efforts.     

Palaeodistribution modelling and genetic evidence highlight differential post‐glacial range shifts of a rain forest conifer distributed across a latitudinal gradient

Aim We examine the range expansion/contraction dynamics during the last glacial cycle of the late‐successional tropical rain forest conifer Podocarpus elatus using a combination of modelling and molecular marker analyses. Specifically, we test whether distributional changes predicted by environmental niche modelling are in agreement with (1) the glacial maximum contractions inferred from the southern fossil record, and (2) population genetic‐based estimates of range disjunctions and demographic dynamics. In addition, we test whether northern and southern ranges are likely to have experienced similar expansion/contraction dynamics. Location Eastern Australian tropical and subtropical rain forests. Methods Environmental niche modelling was completed for three time periods during the last glacial cycle and was interpreted in light of the known palynology. We collected 109 samples from 32 populations across the entire range of P. elatus. Six microsatellite loci and Bayesian coalescence analysis were used to infer population expansion/contraction dynamics, and five sequenced loci (one plastid and four nuclear) were used to quantify genetic structure/diversity. Results Environmental niche modelling suggested that the northern and southern ranges of P. elatus experienced different expansion/contraction dynamics. In the northern range, the habitat suitable for P. elatus persisted in a small refugial area during the Last Glacial Maximum (LGM, 21 ka) and then expanded during the post‐glacial period. Conversely, in the south suitable habitat was widespread during the LGM but subsequently contracted. These differential dynamics were supported by Bayesian analyses of the population genetic data (northern dispersal) and are consistent with the greater genetic diversity in the south compared with the north. A contact zone between the two genetically divergent groups (corresponding to the Macleay Overlap Zone) was supported by environmental niche modelling and molecular analyses. Main conclusions The climatic fluctuations of the Quaternary have differentially impacted the northern and southern ranges of a broadly distributed rain forest tree in Australia. Recurrent contraction/expansion cycles contributed to the genetic distinction between northern and southern distributions of P. elatus. By combining molecular and environmental niche modelling evidence, this unique study undermines the general assumption that broadly distributed species respond in a uniform way to climate change.     

Quaternary climate instability is correlated with patterns of population genetic variability in Bombus huntii

Climate oscillations have left a significant impact on the patterns of genetic diversity observed in numerous taxa. In this study, we examine the effect of Quaternary climate instability on population genetic variability of a bumble bee pollinator species, Bombus huntii in western North America. Pleistocene and contemporary B. huntii habitat suitability (HS) was estimated with an environmental niche model (ENM) by associating 1,035 locality records with 10 bioclimatic variables. To estimate genetic variability, we genotyped 380 individuals from 33 localities at 13 microsatellite loci. Bayesian inference was used to examine population structure with and without a priori specification of geographic locality. We compared isolation by distance (IBD) and isolation by resistance (IBR) models to examine population differentiation within and among the Bayesian inferred genetic clusters. Furthermore, we tested for the effect of environmental niche stability (ENS) on population genetic diversity with linear regression. As predicted, high‐latitude B. huntii habitats exhibit low ENS when compared to low‐latitude habitats. Two major genetic clusters of B. huntii inhabit western North America: (a) a north genetic cluster predominantly distributed north of 28°N and (b) a south genetic cluster distributed south of 28°N. In the south genetic cluser, both IBD and IBR models are significant. However, in the north genetic cluster, IBD is significant but not IBR. Furthermore, the IBR models suggest that low‐latitude montane populations are surrounded by habitat with low HS, possibly limiting dispersal, and ultimately gene flow between populations. Finally, we detected high genetic diversity across populations in regions that have been climatically unstable since the last glacial maximum (LGM), and low genetic diversity across populations in regions that have been climatically stable since the LGM. Understanding how species have responded to climate change has the potential to inform management and conservation decisions of both ecological and economic concerns.     

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

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

Disentangling the effects of geographic peripherality and habitat suitability on neutral and adaptive genetic variation in Swiss stone pine

It is generally accepted that the spatial distribution of neutral genetic diversity within a species’ native range mostly depends on effective population size, demographic history, and geographic position. However, it is unclear how genetic diversity at adaptive loci correlates with geographic peripherality or with habitat suitability within the ecological niche. Using exome‐wide genomic data and distribution maps of the Alpine range, we first tested whether geographic peripherality correlates with four measures of population genetic diversity at > 17,000 SNP loci in 24 Alpine populations (480 individuals) of Swiss stone pine (Pinus cembra) from Switzerland. To distinguish between neutral and adaptive SNP sets, we used four approaches (two gene diversity estimates, F_ST outlier test, and environmental association analysis) that search for signatures of selection. Second, we established ecological niche models for P. cembra in the study range and investigated how habitat suitability correlates with genetic diversity at neutral and adaptive loci. All estimates of neutral genetic diversity decreased with geographic peripherality, but were uncorrelated with habitat suitability. However, heterozygosity (H_e) at adaptive loci based on Tajima's D declined significantly with increasingly suitable conditions. No other diversity estimates at adaptive loci were correlated with habitat suitability. Our findings suggest that populations at the edge of a species' geographic distribution harbour limited neutral genetic diversity due to demographic properties. Moreover, we argue that populations from suitable habitats went through strong selection processes, are thus well adapted to local conditions, and therefore exhibit reduced genetic diversity at adaptive loci compared to populations at niche margins.     

Evolutionary history of the scorpionfly Dicerapanorpa magna (Mecoptera,Panorpidae)

Climate changes can have fundamental impacts on the distributional patterns of montane species, and range shifts frequently lead to allopatric divergence followed by the establishment of secondary contact zones. Many European and North American organisms have retreated to southern refugia during glacial periods and colonized northward during postglacial periods, but little is known about the evolutionary response of cold‐adapted insects to Pleistocene climate changes in eastern Asia. The scorpionfly Dicerapanorpa magna (Chou), with cold temperate habitat preference and weak dispersal ability, provides a good model system to explore how climate changes have influenced the distribution and divergence of cold‐adapted insects in eastern Asia. This study reconstructed the demographic dynamics and evolutionary history of D. magna with phylogeographic approaches, and predicted the species’ suitable areas under the Last Glacial Maximum (LGM) and current scenarios with the ecological niche modelling analysis. The mitochondrial cytochrome c oxidase subunit I resolved three phylogenetic lineages in D. magna dating back to Pleistocene, corresponding well with the geographically isolated Qinling, Bashan and Minshan Mountains. The ecological niche modelling recovered the suitable habitats for D. magna were the Qinling and Bashan Mountains under LGM and current conditions. The three lineages of D. magna might be in a process of incipient speciation, and likely derived their current distribution from separate glacial origins, followed by vicariance and divergence.     

Genetic turnovers and northern survival during the last glacial maximum in European brown bears

The current phylogeographic pattern of European brown bears (Ursus arctos) has commonly been explained by postglacial recolonization out of geographically distinct refugia in southern Europe, a pattern well in accordance with the expansion/contraction model. Studies of ancient DNA from brown bear remains have questioned this pattern, but have failed to explain the glacial distribution of mitochondrial brown bear clades and their subsequent expansion across the European continent. We here present 136 new mitochondrial sequences generated from 346 remains from Europe, ranging in age between the Late Pleistocene and historical times. The genetic data show a high Late Pleistocene diversity across the continent and challenge the strict confinement of bears to traditional southern refugia during the last glacial maximum (LGM). The mitochondrial data further suggest a genetic turnover just before this time, as well as a steep demographic decline starting in the mid‐Holocene. Levels of stable nitrogen isotopes from the remains confirm a previously proposed shift toward increasing herbivory around the LGM in Europe. Overall, these results suggest that in addition to climate, anthropogenic impact and inter‐specific competition may have had more important effects on the brown bear's ecology, demography, and genetic structure than previously thought.     

Phylogeographic analysis of North American populations of the parasitic herbaceous plant Monotropa hypopitys L. reveals a complex history of range expansion from multiple late glacial refugia

Aim We carried out a phylogeographic study across the range of the herbaceous plant species Monotropa hypopitys L. in North America to determine whether its current disjunct distribution is due to recolonization from separate eastern and western refugia after the Last Glacial Maximum (LGM). Location North America: Pacific Northwest and north‐eastern USA/south‐eastern Canada. Methods Palaeodistribution modelling was carried out to determine suitable climatic regions for M. hypopitys at the LGM. We analysed between 155 and 176 individuals from 39 locations spanning the species’ entire range in North America. Sequence data were obtained for the chloroplast rps2 gene (n = 168) and for the nuclear ITS region (n = 158). Individuals were also genotyped for eight microsatellite loci (n = 176). Interpolation of diversity values was used to visualize the range‐wide distribution of genetic diversity for each of the three marker classes. Minimum spanning networks were constructed showing the relationships between the rps2 and ITS haplotypes, and the geographical distributions of these haplotypes were plotted. The numbers of genetic clusters based on the microsatellite data were estimated using Bayesian clustering approaches. Results The palaeodistribution modelling indicated suitable climate envelopes for M. hypopitys at the LGM in both the Pacific Northwest and south‐eastern USA. High levels of genetic diversity and endemic haplotypes were found in Oregon, the Alexander Archipelago, Wisconsin, and in the south‐eastern part of the species’ distribution range. Main conclusions Our results suggest a complex recolonization history for M. hypopitys in North America, involving persistence in separate eastern and western refugia. A generally high degree of congruence between the different marker classes analysed indicated the presence of multiple refugia, with at least two refugia in each area. In the west, putative refugia were identified in Oregon and the Alexander Archipelago, whereas eastern refugia may have been located in the southern part of the species’ current distribution, as well as in the ‘Driftless Area’. These findings are in contrast to a previous study on the related species Orthilia secunda, which has a similar disjunct distribution to M. hypopitys, but which appears to have recolonized solely from western refugia.     

Postglacial colonization of Europe by the barbastelle bat: agreement between molecular data and past predictive modelling

The barbastelle (Barbastella barbastellus) is a rare forest bat with a wide distribution in Europe. Here, we combine results from the analysis of two mtDNA fragments with species distribution modelling to determine glacial refugia and postglacial colonization routes. We also investigated whether niche conservatism occurs in this species. Glacial refugia were identified in the three southern European peninsulas: Iberia, Italy and the Balkans. These latter two refugia played a major role in the postglacial colonization process, with their populations expanding to England and central Europe, respectively. Palaeo‐distribution models predicted that suitable climatic conditions existed in the inferred refugia during the last glacial maximum (LGM). Nevertheless, the overlap between the current and the LGM distributions was almost inexistent in Italy and in the Balkans, meaning that B. barbastellus populations were forced to shift range between glacial and interglacial periods, a process that probably caused some local extinctions. In contrast, Iberian populations showed a ‘refugia within refugium’ pattern, with two unconnected areas containing stable populations (populations that subsisted during both glacial and interglacial phases). Moreover, the match between LGM models and the refugial areas determined by molecular analysis supported the hypothesis of niche conservatism in B. barbastellus. We argue that geographic patterns of genetic structuring, altogether with the modelling results, indicate the existence of four management units for conservation: Morocco, Iberia, Italy and UK, and Balkans and central Europe. In addition, all countries sampled possessed unique gene pools, thus stressing the need for the conservation of local populations.     

Palaeodistribution modelling does not support disjunct Pleistocene refugia in several Central American plant taxa

Aim We combine evidence from palaeoniche modelling studies of several tree species to estimate the extent of Central American forest during the Last Glacial Maximum (LGM). In particular, we ask whether the distributions of these species are likely to have changed since the LGM, and whether LGM distributions coincide with previously proposed Pleistocene refugia in this area. Location Central American wet and seasonally dry forests. Methods We developed ecological niche models using two simulations of Pleistocene climate and occurrence data for 15 Neotropical plant species. We focused on palaeodistribution models of three ‘focal’ tree species that occur in wet and seasonally dry Central American forests, where recent phylogeographic data suggest Pleistocene differentiation coincident with previously proposed refugia. We added predictions from six wet‐forest and six seasonally dry‐forest obligate plant species to gauge whether Pleistocene range shifts were specific to habitat type. Correlation analyses were performed between projected LGM and present distributions, LGM distributions and previously proposed refugia. We also asked whether modelled palaeodistributions were smaller than their current extents. Results According to our models, the ranges of the study species were not reduced during the LGM, and did not correlate with refugial models, regardless of habitat type. Relative range sizes between present and LGM distributions did not indicate significant range changes since the LGM. However, relative range sizes differed overall between the two palaeoclimate models. Main conclusions Many of the modelled palaeodistributions of study species were not restricted to refugia during the LGM, regardless of forest type. While constrained from higher elevations, most species found suitable habitat at coastal margins and on newly exposed land due to lowered sea levels during the LGM. These results offer no corroboration for Pleistocene climate change as a driver of genetic differentiation in the ‘focal’ species. We offer alternative explanations for genetic differentiation found in plant species in this area.     

Distance,elevation and environment as drivers of diversity and divergence in bumble bees across latitude and altitude

Identifying drivers of dispersal limitation and genetic differentiation is a key goal in biogeography. We examine patterns of population connectivity and genetic diversity using restriction site‐associated DNA sequencing (RADseq) in two bumble bee species, Bombus vosnesenskii and Bombus bifarius, across latitude and altitude in mountain ranges from California, Oregon and Washington, U.S.A. Bombus vosnesenskii, which occurs across a broader elevational range at most latitudes, exhibits little population structure while B. bifarius, which occupies a relatively narrow higher elevation niche across most latitudes, exhibits much stronger population differentiation, although gene flow in both species is best explained by isolation with environmental niche resistance. A relationship between elevational habitat breadth and genetic diversity is also apparent, with B. vosnesenskii exhibiting relatively consistent levels of genetic diversity across its range, while B. bifarius has reduced genetic diversity at low latitudes, where it is restricted to high‐elevation habitat. The results of this study highlight the importance of the intersect between elevational range and habitat suitability in influencing population connectivity and suggest that future climate warming will have a fragmenting effect even on populations that are presently well connected, as they track their thermal niches upward in montane systems.     

Phylogeography,population genetics and distribution modelling reveal vulnerability of Scirpus longii (Cyperaceae) and the Atlantic Coastal Plain Flora to climate change

A proactive approach to conservation must be predictive, anticipating how habitats will change and which species are likely to decline or prosper. We use composite species distribution modelling to identify suitable habitats for 18 members of the North American Atlantic Coastal Plain Flora (ACPF) since the Last Glacial Maximum and project these into the future. We then use Scirpus longii (Cyperaceae), a globally imperiled ACPF sedge with many of the characteristics of extinction vulnerability, as a case study. We integrate phylogeographical and population genetic analyses and species distribution modelling to develop a broad view of its current condition and prognosis for conservation. We use genotyping‐by‐sequencing to characterize the genomes of 142 S. longii individuals from 20 populations distributed throughout its range (New Jersey to Nova Scotia). We measure the distribution of genetic diversity in the species and reconstruct its phylogeographical history using the snapp and rase models. Extant populations of S. longii originated from a single refugium south of the Laurentide ice sheet around 25 ka. The genetic diversity of S. longii is exceedingly low, populations exhibit little genetic structure and the species is slightly inbred. Projected climate scenarios indicate that nearly half of extant populations of S. longii will be exposed to unsuitable climate by 2070. Similar changes in suitable habitat will occur for many other northern ACPF species—centres of diversity will shift northward and Nova Scotia may become the last refuges for those species not extinguished.     

Phylogeography of the Northern Alligator Lizard (Squamata,Anguidae): Hidden diversity in a western endemic

Western North America includes the California Floristic Province and the Pacific Northwest, biologically diverse regions highlighted by a complex topography, geology, climate and history. A number of animals span these regions and show distinctive patterns of dispersal, vicariance and lineage diversification. Examining phylogeographic patterns in the fauna of this area aids in our understanding of the forces that have contributed to the generation and maintenance of regional biodiversity. Here, we investigate the biogeography and population structure of the Northern Alligator Lizard (Elgaria coerulea), a wide‐ranging anguid endemic to western North America. We sequenced two mtDNA fragments (ND2 and ND4) for 181 individuals across the range of the species and analysed these data with phylogenetic approaches to infer population and biogeographic history, and date major divergences within the taxon. We further used Bayesian clustering methods to assess major patterns of population structure and performed ecological niche modelling (ENM) to aid in our interpretation of geographic structure and diversification of E. coerulea lineages. Our phylogeographic examination of E. coerulea uncovered surprising diversity and structure, recovering 10 major lineages, each with substantial geographic substructure. While some divergences within the species are relatively old (Pliocene, 5.3–2.6 mya), most intraspecific variation appears to be of more recent origin (Pleistocene, 2.6 mya‐11,700 ya). Current diversity appears to have arisen in the Sierra Nevada Mountains and spread west and north since the Pliocene. Finally, our ENMs suggest that much of the Coast Ranges in California provided ideal habitat during the Last Glacial Maxima (LGM) that has since contracted dramatically and shifted northwards, whereas significant portions of the Sierra Nevada were unsuitable during the LGM and have since become more suitable. Interestingly, E. coerulea shares a number of genetic boundaries with other sympatric taxa, suggesting common historical events and geomorphological features have shaped the biota of this region.