Ephemeral Pleistocene woodlands connect the dots for highland rattlesnakes of the Crotalus intermedius group

Aim To test how Pleistocene climatic changes affected diversification of the Crotalus intermedius species complex. Location Highlands of Mexico and the south‐western United States (Arizona). Methods We synthesize the matrilineal genealogy based on 2406 base pairs of mitochondrial DNA sequences, fossil‐calibrated molecular dating, reconstruction of ancestral geographic ranges, and climate‐based modelling of species distributions to evaluate the history of female dispersion. Results The presently fragmented distribution of the C. intermedius group is the result of both Neogene vicariance and Pleistocene pine–oak habitat fragmentation. Most lineages appear to have a Quaternary origin. The Sierra Madre del Sur and northern Sierra Madre Oriental are likely to have been colonized during this time. Species distribution models for the Last Glacial Maximum predict expansions of suitable habitat for taxa in the southern Sierra Madre Occidental and northern Sierra Madre Oriental. Main conclusions Lineage diversification in the C. intermedius group is a consequence of Pleistocene climate cycling. Distribution models for two sister taxa in the northern and southern Sierra Madre Occidental and northern Sierra Madre Oriental during the Last Glacial Maximum provide evidence for the expansion of pine–oak habitat across the Central Mexican Plateau. Downward displacement and subsequent expansions of highland vegetation across Mexico during cooler glacial cycles may have allowed dispersal between highlands, which resulted in contact between previously isolated taxa and the colonization of new habitats.     

Identification of glacial refugia in south-eastern North America by phylogeographical analyses of a forest understorey plant, Trillium cuneatum

Aim We examine several hypotheses emerging from biogeographical and fossil records regarding glacial refugia of a southern thermophilic plant species. Specifically, we investigated the glacial history and post‐glacial colonization of a forest understorey species, Trillium cuneatum. We focused on the following questions: (1) Did T. cuneatum survive the Last Glacial Maximum (LGM) in multiple refugia, and (if so) where were they located, and is the modern genetic structure congruent with the fossil record‐based reconstruction of refugia for mesic deciduous forests? (2) What are the post‐glacial colonization patterns in the present geographical range? Location South‐eastern North America. Methods We sampled 45 populations of T. cuneatum throughout its current range. We conducted phylogeographical analyses based on maternally inherited chloroplast DNA (cpDNA haplotypes) and used TCS software to reconstruct intraspecific phylogeny. Results We detected six cpDNA haplotypes, geographically highly structured into non‐overlapping areas. With one exception, none of the populations had mixed haplotype composition. TCS analysis resulted in two intraspecific cpDNA lineages, with one clade subdivided further by shallower diversification. Main conclusions Our investigation revealed that T. cuneatum survived the LGM in multiple refugia, belonging to two (western, eastern) genealogical lineages geographically structured across south‐eastern North America. The western clade is confined to the south‐western corner of T. cuneatum’s modern range along the Lower Mississippi Valley, where fossil records document a major refugium of mesic deciduous forest. For the eastern clade, modern patterns of cpDNA haplotype distribution suggest cryptic vicariance, in the form of forest contractions and subsequent expansions associated with Pleistocene glacial cycles, rather than simple southern survival and subsequent northward colonization. The north–south partitioning of cpDNA haplotypes was unexpected, suggesting that populations of this rather southern thermophilic species may have survived in more northern locations than initially expected based on LGM climate reconstruction, and that the Appalachian Mountains functioned as a barrier to the dispersal of propagules originating in more southern refugia. Furthermore, our results reveal south‐west to north‐east directionality in historical migration through the Valley and Ridge region of north‐west Georgia.     

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.     

Interglacial genetic diversification of Moussonia deppeana (Gesneriaceae), a hummingbird‐pollinated,cloud forest shrub in northern Mesoamerica

Recent empirical work on cloud forest‐adapted species supports the role of both old divergences across major geographical areas and more recent divergences attributed to Pleistocene climate changes. The shrub Moussonia deppeana is distributed in northern Mesoamerica, with geographically disjunct populations. Based on sampling throughout the species range and employing plastid and nuclear markers, we (i) test whether the fragmented distribution is correlated with main evolutionary lineages, (ii) reconstruct its phylogeographical history to infer the history of cloud forest in northern Mesoamerica and (iii) evaluate a set of refugia/vicariance scenarios for the region and demographic patterns of the populations whose ranges expanded and tracked cloud forest conditions during the Last Glacial Maximum. We found a deep evolutionary split in M. deppeana about 6–3 Ma, which could be consistent with a Pliocene divergence. Comparison of variation in plastid and nuclear markers revealed several lineages mostly congruent with their isolated geographical distribution and restricted gene flow among groups. Results of species distribution modelling and coalescent simulations fit a model of multiple refugia diverging during interglacial cycles. The demographic history of M. deppeana is not consistent with an expanding–contracting cloud forest archipelago model during the Last Glacial Maximum. Instead, our data suggest that populations persisted across the geographical range throughout the glacial cycles, and experienced isolation and divergence during interglacial periods.     

Mitochondrial DNA variation and Quaternary range dynamics in the endangered Yellow Spotted Mountain Newt,Neurergus derjugini (Caudata,Salamandridae)

In temperate regions of the Earth Pleistocene, climatic fluctuations significantly influenced distribution of species. However, little is known on how glacial and interglacial cycles affected range dynamics of the species occupying lower latitudes. In this study, we investigated mitochondrial DNA (mtDNA) variation and reconstructed the potential current and past (during the mid‐Holocene, 6 ka BP, and the Last Glacial Maximum, LGM, 21 ka BP) distribution of Neurergus derjugini, an endangered amphibian species endemic to the mid‐Zagros Mountains in Iran and Iraq. Six haplotypes identified in the control region (D‐loop) form a well‐supported monophyletic clade, distinct from other Neurergus species and revealing a sister relationship to Neurergus kaiseri. Nucleotide diversity quantifying mean divergence between the sequences is low and does not support the recognition of distinct evolutionary lineages in Neurergus derjugini. The landscape connectivity analysis and the haplotype parsimony network reveal higher gene flow rate between the breeding streams in the southern part of the range, while the northern populations are more isolated. The potential distribution of Neurergus derjugini is restricted to valleys close to mountain tops, wherein very high elevations and dry habitats appear to be unsuitable. During the mid‐Holocene and LGM conditions, the range of the species may have been more extended and shifted to lower elevations. These findings show retraction of the Neurergus derjugini range during the Quaternary and indicate that range dynamics of the species occupying lower latitudes may not follow a scenario of glacial retraction and postglacial expansion.     

Pleistocene climate changes shaped the divergence and demography of Asian populations of the great tit Parus major: evidence from phylogeographic analysis and ecological niche models

Different scales and frequencies of glaciations developed in Europe and Asia during the Pleistocene. Because species’ responses to climate change are influenced by interactive factors including ecology and local topography, the pattern and tempo of species diversification may vary significantly across regions. The great tit Parus major is a widespread Eurasian passerine with a range that encircles the central Asian desert and high‐altitude areas of the Tibetan Plateau. A number of genetic studies have assessed the effect of paleo‐climate changes on the distribution of the European population. However, none have comprehensively addressed how paleo‐climate change affected the distribution of the great tit in China, an apparent hotspot of P. major subspecific diversity. Here, we describe likely paleo‐climatic effects on P. major populations in China based on a combination of phylogeography and ecological niche models (ENMs). We sequenced three mitochondrial DNA markers from 28 populations (213 individuals), and downloaded 112 sequences from outside its Chinese range. As the first step in clarifying the intra‐specific relationships among haplotypes, we attempted to clarify the divergence and demography of populations in China. Phylogeographic analysis revealed that P. major is comprised of five highly divergent clades with geographic breaks corresponding to steep mountains and dry deserts. A previously undescribed monophyletic clade with high genetic diversity, stable niches and a long and independent evolutionary history was detected in the mountainous areas of southwest China. The estimated times at which these clades diverged was traced back to the Early‐Middle Pleistocene (2.19–0.61 mya). Contrary to the post‐LGM (the Last Glacial Maximum) expansion of European populations, demographic history indicates that Asian populations expanded before the LGM after which they remained relatively stable or grew slowly through the LGM. ENMs support this conclusion and predict a similar distribution in the present and the LGM. Our genetic and ecological results demonstrate that Pleistocene climate changes shaped the divergence and demography of P. major in China.     

Evolutionary history of Apocheima cinerarius (Lepidoptera: Geometridae), a female flightless moth in northern China

The alterations in the phylogeographical structures of insects in response to the uplift of the Qinghai–Tibet Plateau and the Quaternary glaciations in eastern Asia, particularly in northern China, remain largely unknown. In this study, we selected Apocheima cinerarius, a moth with flightless females, using molecular data (complete mitochondrial genomes and nuclear data) and ecological niche modelling (ENM) to investigate the effects of paleoclimatic changes on the evolutionary history of insects in the area of northern China. The phylogenetic tree of complete mitochondrial genomes indicated that there were two lineages, the western and eastern lineages. The nuclear gene analyses also detected unique haplotypes in each lineage. Time of the most recent common ancestor (TMRCA) of the two lineages was approximately in Early–Middle Pleistocene. Bayesian skyline plots revealed that the western lineage underwent a population expansion event after the Last Glacial Maximum, whereas the eastern lineage underwent expansion between the Last Interglacial and the Last Glacial Maximum. Our results suggest that A. cinerarius expanded eastward from western sites until the moth was distributed across the entire region of northern China. Then, A. cinerarius underwent contraction into isolated glacial refugia followed by subsequent expansion driven by Pleistocene climate changes, which established a narrow sympatric area. Our results indicate that the Quaternary environmental fluctuations had profound influences on the diversification and demography of an insect in northern China, and the same species in north‐western China and north‐eastern China have different demographic histories.     

REVISITING WALLACE'S HAUNT: COALESCENT SIMULATIONS AND COMPARATIVE NICHE MODELING REVEAL HISTORICAL MECHANISMS THAT PROMOTED AVIAN POPULATION DIVERGENCE IN THE MALAY ARCHIPELAGO

Sundaland, a biogeographic region of Southeast Asia, is a major biodiversity hotspot. However, little is known about the relative importance of Pleistocene habitat barriers and rivers in structuring populations and promoting diversification here. We sampled 16 lowland rainforest bird species primarily from peninsular Malaysia and Borneo to test the long‐standing hypothesis that animals on different Sundaic landmasses intermixed extensively when lower sea‐levels during the Last Glacial Maximum (LGM) exposed land‐bridges. This hypothesis was rejected in all but five species through coalescent simulations. Furthermore, we detected a range of phylogeographic patterns; Bornean populations are often genetically distinct from each other, despite their current habitat connectivity. Environmental niche modeling showed that the presence of unsuitable habitats between western and eastern Sundaland during the LGM coincided with deeper interpopulation genetic divergences. The location of this habitat barrier had been hypothesized previously based on other evidence. Paleo‐riverine barriers are unlikely to have produced such a pattern, but we cannot rule out that they acted with habitat changes to impede population exchanges across the Sunda shelf. The distinctiveness of northeastern Borneo populations may be underlied by a combination of factors such as rivers, LGM expansion of montane forests and other aspects of regional physiography.     

Effects of Pleistocene climate changes on species ranges and evolutionary processes in the Neotropical Atlantic Forest

The effects of global glaciations on the distribution of organisms is an essential element of many diversification models. However, the empirical evidence supporting this idea is mixed, in particular with respect to explaining tropical forest evolution. In the present study, we evaluated the impacts of range shifts associated with Pleistocene global glacial cycles on the evolution of tropical forests. In particular, we tested the predictions: (1) that population genetic structure increases with fragmentation variation between the present and the Last Glacial Maximum (LGM) and also (2) with geographical range instability; and (3) that genetic diversity increases with range stability and (4) decreases with fragmentation variation between periods. To address our predictions, we studied population genetic structures and modelled present and past distributions of 15 Atlantic Forest (AF) endemic birds. Afterwards, we evaluated the relationship of population genetic parameters with metrics of species range shifts between the present and the LGM. We found that geographical ranges of AF birds changed in concert with Pleistocene glacial cycles but, unexpectedly, our findings suggest the novel idea that ranges during glacial maxima were slightly larger on average, as well as equally fragmented and displaced from the interglacial ranges. Our findings suggest that range shifts over the late Pleistocene impacted on the diversification of forest organisms, although they did not show that those range shifts had a strong effect. We found that a combination of fragmentation variation across time, small current range size, and range stability increased population genetic structure. However, neither fragmentation, nor range stability affected genetic diversity. Our study showed that evolutionary responses to range shifts across AF birds have a high variance, which could explain the mixed support given by single‐species studies to the action of Pleistocene range shifts on population evolution.     

The critical role of local refugia in postglacial colonization of Chinese pine: joint inferences from DNA analyses,pollen records,and species distribution modeling

The importance of long‐distance migration from low to high latitudes relative to local spread from northern refugia after the Last Glacial Maximum (LGM) remains a focus of debate for many temperate tree species. We assessed the dynamics of Chinese pine Pinus tabulaeformis, a widespread species endemic to northern China, since the LGM by integrating cytoplasmic DNA data, mapped pollen records and ecological niche modeling. Genetic variation among 544 individuals from 50 populations spanning the entire natural species range revealed eight genetic clusters with distinct geographic distribution, indicating glacial lineages likely originating from multiple local microrefugia. Palynological evidence suggested that the northernmost part of the natural distribution originated from local postglacial spread. Niche modeling indicated high probability of the species being present in the area of the Loess Plateau and coastal areas north of the Yangtze River during the LGM. The three lines of evidence jointly suggest that the species persisted through the last glaciation in the mountains surrounding the Loess Plateau of northern China and that the current distribution of the species originated primarily from the spread of local refugial populations, instead of long‐distance migration. These results cast doubt on the notion that Chinese pine migrated from areas south of the Yangtze River and underscore the importance of northern refugia.     

Distributional and demographic consequences of Pleistocene climate fluctuations for a marine demersal fish in the north-eastern Atlantic

Aim The Pleistocene glaciations were the most significant historical event during the evolutionary life span of most extant species. However, little is known about the consequences of these climate changes for the distribution and demography of marine animals of the north‐eastern Atlantic. The present study focuses on the phylogeographic and demographic patterns of the sand goby, Pomatoschistus minutus (Teleostei: Gobiidae), a small marine demersal fish. Location North‐eastern Atlantic, Mediterranean, Irish, North and Baltic seas. Methods Analysis was carried out by sequencing the mtDNA cytochrome b gene of sand gobies from 12 localities throughout the species’ range, and using this information in combination with published data of allozyme markers and mtDNA control region sequences. Several phylogenetic methods and a network analysis were used to explore the phylogeographic pattern. The historical demography of P. minutus was studied through a mismatch analysis and a Bayesian skyline plot. Results Reciprocal monophyly was found between a Mediterranean Sea (MS) clade and an Atlantic Ocean (AO) clade, both with a Middle Pleistocene origin. The AO Clade contains two evolutionary significant units (ESUs): the Iberian Peninsula (IB) Group and the North Atlantic (NA) Group. These two groups diverged during Middle Pleistocene glacial cycles. For the NA Group there is evidence for geographic sorting of the ancestral haplotypes with recent radiations in the Baltic Sea, Irish Sea, North Sea and Bay of Biscay. The demographic histories of the Mediterranean Clade and the two Atlantic ESUs were influenced mainly by expansions dated as occurring during the Middle Pleistocene glaciations and post‐Eem, respectively. Main conclusions The pre‐LGM (Last Glacial Maximum) subdivision signals were not erased for P. minutus during the LGM. Middle Pleistocene glaciations yielded isolated and differently evolving sets of populations. In contrast to the case for most other taxa, only the northern Atlantic group contributed to the post‐glacial recolonization. The historical demography of Mediterranean sand gobies was influenced mainly by Middle Pleistocene glaciations, in contrast to that of the Atlantic populations, which was shaped by Late Pleistocene expansions.     

History and evolution of alpine plants endemic to the Qinghai-Tibetan Plateau: Aconitum gymnandrum (Ranunculaceae)

How Quaternary climatic oscillations affected range distributions and intraspecific divergence of alpine plants on the Qinghai‐Tibetan Plateau (QTP) remains largely unknown. Here, we report a survey of chloroplast DNA (cpDNA) and nuclear ribosomal internal transcribed spacer (ITS) DNA variation aimed at exploring the phylogeographical history of the QTP alpine endemic Aconitum gymnandrum. We sequenced three cpDNA fragments (rpl20–rps12 intergenic spacer, the trnV intron and psbA‐trnH spacer) and also the nuclear (ITS) region in 245 individuals from 23 populations sampled throughout the species’ range. Two distinct lineages, with eastern and western geographical distributions respectively, were identified from a phylogenetic analysis of ITS sequence variation. Based on a fast substitution rate, these were estimated to have diverged from each other in the early Pleistocene approximately 1.45 Ma. The analysis of cpDNA variation identified nine chlorotypes that clustered into two major clades that were broadly congruent in geographical distribution with the two ITS lineages. The east–west split of cpDNA divergence was supported by an amova which partitioned approximately half of the total variance between these two groups of populations. Analysis of the spatial distribution of chlorotypes showed that each clade was subdivided into two groups of populations such that a total of four population groups existed in the species. It is suggested that these different groups derive from four independent glacial refugia that existed during the Last Glacial Maximum (LGM), and that three of these refugia were located at high altitude on the QTP platform itself at that time. Coalescent simulation of chlorotype genealogies supported both an early Pleistocene origin of the two main cpDNA clades and also the ‘four‐refugia’ hypothesis during the LGM. Two previous phylogeographical studies of QTP alpine plants indicated that such plants retreated to refugia at the eastern/south‐eastern plateau edge during the LGM and/or previous glacial maxima. However, the results for A. gymnandrum suggest that at least some of these cold‐tolerant species may have also survived centrally on the QTP platform throughout the Quaternary.     

HOLARCTIC PHYLOGEOGRAPHY OF AN ASEXUAL SPECIES COMPLEX I. MITOCHONDRIAL DNA VARIATION IN ARCTIC DAPHNIA

Pleistocene glacial cycles undoubtedly altered the evolutionary trajectories of many taxa, yet few studies have examined the impact of such events on genetic differentiation and phylogeography at large geographic scales. Here we present the results of a circumarctic survey of mitochondrial DNA diversity in members of the Daphnia pulex complex. The analysis involved the survey of restriction site polymorphisms in a 2100-bp fragment of the NADH-4 (ND4) and NADH-5 (ND5) genes for 276 populations representing the two major groups (tenebrosa and pulicaria) in this complex across their Holarctic range. A comparison of the distribution patterns for seven clades in this complex revealed very clear phylogeographic structuring. Most notably, pulicaria group lineages were restricted primarily to the Nearctic, with some colonization of formerly glaciated portions of northern Europe. This group was not detected from vast expanses of northern Eurasia, including the Beringian glacial refuge. In contrast, tenebrosa group haplotypes showed considerable intercontinental divergence between Eurasian and North American lineages, but were absent from Greenland and Iceland, as well as the Canadian arctic archipelago. Dispersal in Eurasia was primarily in a westerly direction from Beringia, whereas dispersal in the Nearctic followed proglacial drainage patterns. Long-distance dispersal of certain lineages was observed in both groups, and variation in haplotype richness and nucleotide diversity allowed us to make inferences about the positioning of putative glacial refugia. Overall, the phylogeographic pattern of diversification in this arctic complex is characterized by the apparently unique postglacial histories for each clade, indicating that even closely allied taxa can respond independently to the allopatric effects of glacial cycles. This is in sharp contrast to other phylogeographic studies of species assemblages from more southern (unglaciated) latitudes, which are often characterized by concordant patterns.     

Beringia and the global dispersal of modern humans

Until recently, the settlement of the Americas seemed largely divorced from the out‐of‐Africa dispersal of anatomically modern humans, which began at least 50,000 years ago. Native Americans were thought to represent a small subset of the Eurasian population that migrated to the Western Hemisphere less than 15,000 years ago. Archeological discoveries since 2000 reveal, however, that Homo sapiens occupied the high‐latitude region between Northeast Asia and northwest North America (that is, Beringia) before 30,000 years ago and the Last Glacial Maximum (LGM). The settlement of Beringia now appears to have been part of modern human dispersal in northern Eurasia. A 2007 model, the Beringian Standstill Hypothesis, which is based on analysis of mitochondrial DNA (mtDNA) in living people, derives Native Americans from a population that occupied Beringia during the LGM. The model suggests a parallel between ancestral Native Americans and modern human populations that retreated to refugia in other parts of the world during the arid LGM. It is supported by evidence of comparatively mild climates and rich biota in south‐central Beringia at this time (30,000‐15,000 years ago). These and other developments suggest that the settlement of the Americas may be integrated with the global dispersal of modern humans.     

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.     

Multiple Pleistocene refugia and recent diversification for Streptocarpus ionanthus (Gesneriaceae) complex: Insights from multiple molecular sources

Streptocarpus ionanthus (Gesneriaceae) is endemic to Tanzania and Kenya, distributed in Tanga, Morogoro, and Kilifi regions. The species houses nine subspecies characterized by complex morphotypes and poorly understood evolutionary relationships, and thus is an ideal model for investigating evolutionary dynamics over time. Using multiple methods, we sought to test our hypothesis that the infraspecific taxa in Str. ionanthus are slightly variable and evolving populations. We first examined the genetic diversity, population differentiation, and phylogeographic structure among the populations of Str. ionanthus using both chloroplast and nuclear markers. We then estimated the divergence time of Str. ionanthus lineages and modeled past and future distribution. Despite Str. ionanthus exhibiting bottleneck events across its range, the populations maintain relatively high genetic diversity attributed to historical population admixture or local adaptation arising from habitat heterogeneity. The phylogeographic and genetic structure revealed a high connection among the Usambara mountains populations, while molecular dating suggested most diversification of haplotypes began ~1.32–0.18 million years ago and intensified toward the present, a conclusion of recent diversification. Phylogenetic relationship of Str. ionanthus cpDNA haplotypes revealed five main lineages with unique haplotypes that could be suggestive of past isolated refugia during the Pleistocene climate shifts. According to niche modeling, the stability of suitable areas during the Last Glacial Maximum (LGM) offered protective micro-habitats that have preserved the genetic diversity of Str. ionanthus to date. In conclusion, our findings suggest a complex Str. ionanthus with slightly variable lineages or populations attributed to multiple refugia and on the verge of divergence.     

On the origin of the Norwegian lemming

The Pleistocene glacial cycles resulted in significant changes in species distributions, and it has been discussed whether this caused increased rates of population divergence and speciation. One species that is likely to have evolved during the Pleistocene is the Norwegian lemming (Lemmus lemmus). However, the origin of this species, both in terms of when and from what ancestral taxon it evolved, has been difficult to ascertain. Here, we use ancient DNA recovered from lemming remains from a series of Late Pleistocene and Holocene sites to explore the species' evolutionary history. The results revealed considerable genetic differentiation between glacial and contemporary samples. Moreover, the analyses provided strong support for a divergence time prior to the Last Glacial Maximum (LGM), therefore likely ruling out a postglacial colonization of Scandinavia. Consequently, it appears that the Norwegian lemming evolved from a small population that survived the LGM in an ice‐free Scandinavian refugium.     

Lineage diversification and historical demography of a sky island salamander,Plethodon ouachitae,from the Interior Highlands

Sky islands provide ideal opportunities for understanding how climatic changes associated with Pleistocene glacial cycles influenced species distributions, genetic diversification, and demography. The salamander Plethodon ouachitae is largely restricted to high‐elevation, mesic forest on six major mountains in the Ouachita Mountains. Because these mountains are separated by more xeric, low‐elevation valleys, the salamanders appear to be isolated on sky islands where gene flow among populations on different mountains may be restricted. We used DNA sequence data along with ecological niche modelling and coalescent simulations to test several hypotheses related to diversifications in sky island habitats. Our results revealed that P. ouachitae is composed of seven well‐supported lineages structured across six major mountains. The species originated during the Late Pliocene, and lineage diversification occurred during the Middle Pleistocene in a stepping stone fashion with a cyclical pattern of dispersal to a new mountain followed by isolation and divergence. Diversification occurred primarily on an east–west axis, which is likely related to the east–west orientation of the Ouachita Mountains and the more favourable cooler and wetter environmental conditions on north slopes compared to south‐facing slopes and valleys. All non‐genealogical coalescent methods failed to detect significant population expansion in any lineages. Bayesian skyline plots showed relatively stable population sizes over time, but indicated a slight to moderate amount of population growth in all lineages starting approximately 10 000–12 000 years ago. Our results provide new insight into sky island diversifications from a previously unstudied region, and further demonstrate that climatic changes during the Pleistocene had profound effects on lineage diversification and demography, especially in species from environmentally sensitive habitats in montane regions.     

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.