Phylogeographic Patterns of mtDNA Variation Revealed Multiple Glacial Refugia for the Frog Species Feirana taihangnica Endemic to the Qinling Mountains

Diversification patterns and demography of montane species are affected by Pleistocene climate fluctuations. Empirical cases from the Qinling Mountains (QM) region, which is a major biogeographic divider of East Asia, are few. We used DNA sequence data of the complete mitochondrial ND2 gene to detect effects of the Pleistocene glaciations on phylogeographic profiles of a frog species, Feirana taihangnica, which is endemic to the QM. Four distinct lineages consisting of seven sublineages were revealed. The strongest signal of biogeographical structure (F _ct = 0.971, P < 0.01) was found when populations were grouped according to these seven sublineages. One narrow secondary contact zone was detected in the middle QM between the lineage from middle QM and the lineage from eastern QM. Coalescent simulations indicated that this species colonized the QM region by a stepping-stone model. Divergences among lineages had likely been influenced by the uplift of the Tibetan Plateau during the late Miocene-to-late Pleistocene, as well as by the Pleistocene climatic cycles. Coalescent simulations also suggested that F. taihangnica populations have persisted through the Pleistocene glacial periods in multiple refugia across the QM region. Demographic analyses indicated that all lineages, except the lineage in the Funiu Mountains, have been experienced postglacial expansion of population size and distribution range. In conclusion, Pleistocene climate fluctuations and tectonic changes during the late Miocene-late Pleistocene have profoundly influenced the phylogeography and historical demography of F. taihangnica.     

Time of diversification in the Cape fauna endemisms,inferred by phylogenetic studies of the genus Iselma (Coleoptera: Meloidae: Eleticinae)

The endemic diversification of the Cape zone fauna and the phylogenetic relationships among the 30 species of the blister beetle genus Iselma are investigated. We analysed morphological, molecular (mtDNA 16S) and combined datasets of characters using a number of approaches (maximum parsimony, maximum likelihood, Bayesian inference analyses). We propose hypotheses of diversification times among taxa from molecular clock analyses. Morphological and molecular analyses produced similar results. According to maximum likelihood molecular studies, radiation within the genus Iselma occurred during the Pliocene and Pleistocene, roughly contemporaneously with the shift in African climate, vegetation and faunal assemblages. Two main lineages, one Namibian and one South African, separated c. 4.9 Ma. Within the South African lineage we identify two groups of endemic species, one in Little Karoo and one that extend towards southern Namaqualand. The remaining South African groups of species are spread through Namaqualand and the southwestern Cape area. These include several endemisms, with different times of diversification during the Pleistocene, probably related in part to glacial cycles. These endemisms are distributed in small areas of the following ecosystems: coastal Strandveld, Succulent Karoo and particularly Fynbos.     

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.     

Phylogeographic and demographic effects of Quaternary climate oscillations in Hexinia polydichotoma (Asteraceae) in Tarim Basin and adjacent areas

In order to investigate the genetic diversity and influence of climate oscillations on evolutionary processes of organisms in Northwest China, we selected Hexinia polydichotoma, a species endemic to China, and examined the phylogeographic structure and historical factors that influenced the evolutionary history of this species in its entire cover range, Tarim Basin and adjacent areas. In the study, 17 haplotypes were identified in H. polydichotoma on the basis of two chloroplast DNA sequences (trnH–psbA and ycf6–psbM). Shown in the network, the two common haplotypes, A and D, respectively, mainly distribute along the northern and southern rims of the basin. The analyses of molecular variance analysis suggest that genetic variation primarily occurs among populations, and all populations were subdivided into five groups by SAMOVA. Geographic range expansion along the southern and northern rims of the basin was supported by the significant value for Tajima’s D and by the unimodal mismatch distribution. It is possible that during the interglacial period of the middle Pleistocene, a large amount of snow and glacial ice melted from the mountains surrounding Tarim Basin. This increased water, the expanding desert, and the dispersal ability of H. polydichotoma were important factors driving not only geographic range expansion, but also the current phylogeographic structure of this species. It is possible that during the middle Pleistocene, the climatic fluctuations resulted in expansion and contraction cycles of river systems and oases, and may consequently have caused population fragmentation.     

Climate impacts on transocean dispersal and habitat in gray whales from the Pleistocene to 2100

Arctic animals face dramatic habitat alteration due to ongoing climate change. Understanding how such species have responded to past glacial cycles can help us forecast their response to today's changing climate. Gray whales are among those marine species likely to be strongly affected by Arctic climate change, but a thorough analysis of past climate impacts on this species has been complicated by lack of information about an extinct population in the Atlantic. While little is known about the history of Atlantic gray whales or their relationship to the extant Pacific population, the extirpation of the Atlantic population during historical times has been attributed to whaling. We used a combination of ancient and modern DNA, radiocarbon dating and predictive habitat modelling to better understand the distribution of gray whales during the Pleistocene and Holocene. Our results reveal that dispersal between the Pacific and Atlantic was climate dependent and occurred both during the Pleistocene prior to the last glacial period and the early Holocene immediately following the opening of the Bering Strait. Genetic diversity in the Atlantic declined over an extended interval that predates the period of intensive commercial whaling, indicating this decline may have been precipitated by Holocene climate or other ecological causes. These first genetic data for Atlantic gray whales, particularly when combined with predictive habitat models for the year 2100, suggest that two recent sightings of gray whales in the Atlantic may represent the beginning of the expansion of this species' habitat beyond its currently realized range.     

Glacial History of a Modern Invader: Phylogeography and Species Distribution Modelling of the Asian Tiger Mosquito Aedes albopictus

Background

The tiger mosquito, Aedes albopictus, is one of the 100 most invasive species in the world and a vector of human diseases. In the last 30 years, it has spread from its native range in East Asia to Africa, Europe, and the Americas. Although this modern invasion has been the focus of many studies, the history of the species’ native populations remains poorly understood. Here, we aimed to assess the role of Pleistocene climatic changes in shaping the current distribution of the species in its native range.

Methodology/Principal Findings

We investigated the phylogeography, historical demography, and species distribution of Ae. albopictus native populations at the Last Glacial Maximum (LGM). Individuals from 16 localities from East Asia were analyzed for sequence variation at two mitochondrial genes. No phylogeographic structure was observed across the study area. Demographic analyses showed a signature of population expansion that started roughly 70,000 years BP. The occurrence of a continuous and climatically suitable area comprising Southeast China, Indochinese Peninsula, and Sundaland during LGM was indicated by species distribution modelling.

Conclusions/Significance

Our results suggest an evolutionary scenario in which, during the last glacial phase, Ae. albopictus did not experience a fragmentation phase but rather persisted in interconnected populations and experienced demographic growth. The wide ecological flexibility of the species probably played a crucial role in its response to glacial-induced environmental changes. Currently, there is little information on the impact of Pleistocene climatic changes on animal species in East Asia. Most of the studies focused on forest-associated species and suggested cycles of glacial fragmentation and post-glacial expansion. The case of Ae. albopictus, which exhibits a pattern not previously observed in the study area, adds an important piece to our understanding of the Pleistocene history of East Asian biota.     

Phylogeography and historical demography of Sicydium salvini in the eastern Pacific

Amphidromy, characterized by freshwater adult and marine larval stages, has been shown to be an important influence on the genetic structure of aquatic populations. Sicydium salvini is a widespread goby species (Teleostei: Gobiidae: Sicydiinae) in the eastern Pacific with a continuous distribution from Mexico to Panama. Here, we use mitochondrial data to infer population genetic and historical demography of this species. Sequences were collected for the mitochondrial gene cytochrome b for 162 specimens sampled across the range of S. salvini with a concentration on rivers in Costa Rica. No genetic structure was detected between regions or rivers in the AMOVA analysis; the phylogeny for this species showed no geographic affinities and very little resolution. Historical demographic analyses indicated a population expansion during the late Pleistocene. These results are consistent with a panmictic population with expansion influenced strongly by Pleistocene glacial cycles and geologic uplift.     

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.     

Multilocus phylogeography of the European ground squirrel: cryptic interglacial refugia of continental climate in Europe

The theory of classical and cryptic Pleistocene refugia is based mainly on historical changes in temperature, and the refugia are usually defined within a latitudinal gradient. However, the gradient of oceanic–continental climate (i.e. longitudinal) was also significantly variable during glacial cycles with important biotic consequences. Range‐wide phylogeography of the European ground squirrel (EGS) was used to interpret the evolutionary and palaeogeographical history of the species in Europe and to shed light on its glacial–interglacial dynamic. The EGS is a steppe‐inhabiting species and the westernmost member of the genus in the Palaearctic region. We have analysed 915 specimens throughout the present natural range by employing mitochondrial DNA sequences (cytochrome b gene) and 12 nuclear microsatellite markers. The reconstructed phylogeography divides the species into two main geographical groups, with deep substructuring within both groups. Bulgaria is the centre of the ancestral area, and it also has the highest genetic diversity within the species. The northernmost group of the EGS survived in the southern part of Pannonia throughout several glacial–interglacial cycles. Animals from this population probably repeatedly colonized areas further to the north and west during the glacial periods, while in the interglacial periods, the EGS distribution contracted back to this Pannonian refugium. The EGS thus represents a species with a glacial expansion/interglacial contraction palaeogeographical dynamics, and the Pannonian and southeastern Balkanian steppes are supported as cryptic refugia of continental climate during Pleistocene interglacials.     

Independent Demographic Responses to Climate Change among Temperate and Tropical Milksnakes (Colubridae: Genus Lampropeltis)

The effects of Late Quaternary climate change have been examined for many temperate New World taxa, but the impact of Pleistocene glacial cycles on Neotropical taxa is less well understood, specifically with respect to changes in population demography. Here, we examine historical demographic trends for six species of milksnake with representatives in both the temperate and tropical Americas to determine if species share responses to climate change as a taxon or by area (i.e., temperate versus tropical environments). Using a multilocus dataset, we test for the demographic signature of population expansion and decline using non-genealogical summary statistics, as well as coalescent-based methods. In addition, we determine whether range sizes are correlated with effective population sizes for milksnakes. Results indicate that there are no identifiable trends with respect to demographic response based on location, and that species responded to changing climates independently, with tropical taxa showing greater instability. There is also no correlation between range size and effective population size, with the largest population size belonging to the species with the smallest geographic distribution. Our study highlights the importance of not generalizing the demographic histories of taxa by region and further illustrates that the New World tropics may not have been a stable refuge during the Pleistocene.     

Isolation by instability: Historical climate change shapes population structure and genomic divergence of treefrogs in the Neotropical Cerrado savanna

Although the impact of Pleistocene glacial cycles on the diversification of the tropical biota was once dismissed, increasing evidence suggests that Pleistocene climatic fluctuations greatly affected the distribution and population divergence of tropical organisms. Landscape genomic analyses coupled with paleoclimatic distribution models provide a powerful way to understand the consequences of past climate changes on the present‐day tropical biota. Using genome‐wide SNP data and mitochondrial DNA, combined with projections of the species distribution across the late Quaternary until the present, we evaluate the effect of paleoclimatic shifts on the genetic structure and population differentiation of Hypsiboas lundii, a treefrog endemic to the South American Cerrado savanna. Our results show a recent and strong genetic divergence in H. lundii across the Cerrado landscape, yielding four genetic clusters that do not seem congruent with any current physical barrier to gene flow. Isolation by distance (IBD) explains some of the population differentiation, but we also find strong support for past climate changes promoting range shifts and structuring populations even in the presence of IBD. Post‐Pleistocene population persistence in four main areas of historical stable climate in the Cerrado seems to have played a major role establishing the present genetic structure of this treefrog. This pattern is consistent with a model of reduced gene flow in areas with high climatic instability promoting isolation of populations, defined here as “isolation by instability,” highlighting the effects of Pleistocene climatic fluctuations structuring populations in tropical savannas.     

Phylogeographic and population genetic analyses reveal Pleistocene isolation followed by high gene flow in a wide ranging,but endangered,freshwater mussel

Freshwater organisms of North America have had their contemporary genetic structure shaped by vicariant events, especially Pleistocene glaciations. Life history traits promoting dispersal and gene flow continue to shape population genetic structure. Cumberlandia monodonta, a widespread but imperiled (IUCN listed as endangered) freshwater mussel, was examined to determine genetic diversity and population genetic structure throughout its range. Mitochondrial DNA sequences and microsatellite loci were used to measure genetic diversity and simulate demographic events during the Pleistocene using approximate Bayesian computation (ABC) to test explicit hypotheses explaining the evolutionary history of current populations. A phylogeny and molecular clock suggested past isolation created two mtDNA lineages during the Pleistocene that are now widespread. Two distinct groups were also detected with microsatellites. ABC simulations indicated the presence of two glacial refugia and post-glacial admixture of them followed by simultaneous dispersal throughout the current range of the species. The Ouachita population is distinct from others and has the lowest genetic diversity, indicating that this is a peripheral population of the species. Gene flow within this species has maintained high levels of genetic diversity in most populations; however, all populations have experienced fragmentation. Extirpation from the center of its range likely has isolated remaining populations due to the geographic distances among them.     

Genetic consequences of interglacial isolation in a steppe bird

In response to climate changes that have occurred during Pleistocene glacial cycles, taxa associated to steppe vegetation might have followed a pattern of historical evolution in which isolation and fragmentation of populations occurred during the short interglacials and expansion events occurred during the long glacial periods, in contrast to the pattern described for temperate species. Here, we use molecular genetic data to evaluate this idea in a steppe bird with Palaearctic distribution, the little bustard (Tetrax tetrax). Overall, extremely low genetic diversity and differentiation was observed among eight little bustard populations distributed in Spain and France. Mismatch distribution analyses showed that most little bustard populations expanded during cooling periods previous to, and just after, the last interglacial period (127,000-111,000 years before present), when steppe habitats were widespread across Europe. Coalescent-based methods suggested that glacial expansions have resulted in substantial admixture in Western Europe due to the existence of different interglacial refugia. Our results are consistent with a model of evolution and genetic consequences of Pleistocene cycles with low between-population genetic differentiation as a result of short-term isolation periods during interglacials and long-term exchange during glacial periods.     

Sailing through the Late Pleistocene: unusual historical demography of an East Asian endemic, the Chinese Hwamei (Leucodioptron canorum canorum), during the last glacial period

Pleistocene climate fluctuations shaped the patterns of genetic diversity observed in extant species. In contrast to Europe and North America where the effects of recent glacial cycles on genetic diversity have been well studied, the genetic legacy of the Late Pleistocene for East Asia, a region of great topographical complexity and presumably milder historical climate, remains poorly understood. We analysed 3.86 kb of the mitochondrial genome of 186 Chinese Hwamei birds, Leucodioptron canorum canorum , and found that contrary to the conventional expectation of population decline during cold periods (stadials), the demographic history of this species shows continuous population growth since the penultimate glacial period (about 170 000 years ago). Refugia were identified in the south, coastal regions, and northern inland areas, implying that topographic complexity played a substantial role in providing suitable habitats for the Chinese Hwamei during cold periods. Intermittent gene flow between these refugia during the warmer periods (interstadials) might have resulted in a large effective population of this bird through the last glacial period.     

Pleistocene Niche Stability and Lineage Diversification in the Subtropical Spider Araneus omnicolor (Araneidae)

The influence of Quaternary climate oscillations on the diversification of the South American fauna is being increasingly explored. However, most of these studies have focused on taxa that are endemic to tropical environments, and relatively few have treated organisms restricted to subtropical biomes. Here we used an integrative phylogeographical framework to investigate the effects of these climate events on the ecological niche and genetic patterns of the subtropical orb-weaver spider Araneus omnicolor (Araneidae). We analyzed the mitochondrial (Cytochrome Oxidase I, COI) and nuclear (Internal Transcribed Subunit II, ITS2) DNA of 130 individuals throughout the species’ range, and generated distribution models in three different climate scenarios [present, Last Glacial Maximum (LGM), and Last Interglacial Maximum (LIG)]. Additionally, we used an Approximate Bayesian Computation (ABC) approach to compare possible demographic scenarios and select the hypothesis that better explains the genetic patterns of A. omnicolor. We obtained high haplotype diversity but low nucleotide variation among sequences. The population structure and demographic analyses showed discrepancies between markers, suggesting male-biased dispersal in the species. The time-calibrated COI phylogenetic inference showed a recent diversification of lineages (Middle/Late Pleistocene), while the paleoclimate modeling indicated niche stability since ~120 Kya. The ABC results agreed with the niche models, supporting a panmictic population as the most likely historical scenario for the species. These results indicate that A. omnicolor experienced no niche or population reductions during the Late Pleistocene, despite the intense landscape modifications that occurred in the subtropical region, and that other factors beside LGM and LIG climate oscillations might have contributed to the demographic history of this species. This pattern may be related to the high dispersal ability and wide environmental tolerance of A. omnicolor, highlighting the need for more phylogeographical studies with invertebrates and other generalist taxa, in order to understand the effects of Quaternary climate changes on Neotropical biodiversity.     

High levels of population subdivision in a morphologically conserved Mediterranean toad (Alytes cisternasii) result from recent,multiple refugia: evidence from mtDNA,microsatellites and nuclear genealogies

Pleistocene glaciations often resulted in differentiation of taxa in southern European peninsulas, producing the high levels of endemism characteristic of these regions (e.g. the Iberian Peninsula). Despite their small ranges, endemic species often exhibit high levels of intraspecific differentiation as a result of a complex evolutionary history dominated by successive cycles of fragmentation, expansion and subsequent admixture of populations. Most evidence so far has come from the study of species with an Atlantic distribution in northwestern Iberia, and taxa restricted to Mediterranean‐type habitats remain poorly studied. The Iberian Midwife toad (Alytes cisternasii) is a morphologically conserved species endemic to southwestern and central Iberia and a typical inhabitant of Mediterranean habitats. Applying highly variable genetic markers from both mitochondrial and nuclear genomes to samples collected across the species’ range, we found evidence of high population subdivision within A. cisternasii. Mitochondrial haplotypes and microsatellites show geographically concordant patterns of genetic diversity, suggesting population fragmentation into several refugia during Pleistocene glaciations followed by subsequent events of geographical and demographic expansions with secondary contact. In addition, the absence of variation at the nuclear β‐fibint7 and Ppp3caint4 gene fragments suggests that populations of A. cisternasii have been recurrently affected by episodes of extinction and recolonization, and that documented patterns of population subdivision are the outcome of recent and multiple refugia. We discuss the evolutionary history of the species with particular interest in the increasing relevance of Mediterranean refugia for the survival of genetically differentiated populations during the Pleistocene glaciations as revealed by studies in co‐distributed taxa.     

Multiple lines of evidence for demographic and range expansion of a temperate species (Hyla sarda) during the last glaciation

Many temperate species experienced demographic and range contractions in response to climatic changes during Pleistocene glaciations. In this study, we investigate the evolutionary history of the Tyrrhenian tree frog Hyla sarda, a species inhabiting the Corsica-Sardinia island system (Western Mediterranean basin). We used sequence analysis of two mitochondrial (overall 1229 bp) and three nuclear (overall 1692 bp) gene fragments to assess the phylogeography and demographic history of this species, and species distribution modelling (SDM) to predict its range variation over time. Phylogeographic, historical demographic and SDM analyses consistently indicate that H. sarda does not conform to the scenario generally expected for temperate species but rather underwent demographic and range expansion mostly during the last glacial phase. Palaeogeographic data and SDM analyses suggest that such expansion was driven by the glaciation-induced increase in lowland areas during marine regression. This unusual scenario suggests that at least some temperate species may not have suffered the adverse effects of glacial climate on their population size and range extent, owing to the mitigating effects of other glaciations-induced palaeoenvironmental changes. We discuss previous clues for the occurrence of such a scenario in other species and some possible challenges with its identification. Early phylogeographic literature suggested that responses to the Pleistocene glacial-interglacial cycles were expected to vary among species and regions. Our results point out that such variation may have been greater than previously thought.     

Population status and distribution modelling of the critically endangered riverine rabbit (Bunolagus monticularis)

The riverine rabbit (Bunolagus monticularis) is one of the most endangered mammals in the world due to fragmentation of its habitat in the semi‐arid Karoo region of South Africa, to which it is endemic. It is an umbrella conservation species for the riparian shrubland associated with the seasonal drainage system of the Karoo, where its presence is an indicator of ecosystem health. In this study, we analysed historical survey data to derive an improved assessment of the current B. monticularis population status and distribution. Geospatial analysis was conducted using geographical information systems, and distribution modelling was performed using Maxent. Extent of occurrence for the species is 54,227 km², and area of occupancy is 2943 km². Estimates of 157–207 mature individuals confirm an alarmingly small species population size, and it appears that no subpopulation has >50 mature individuals. Our findings thus support the continued classification of this species as ‘critically endangered’ under IUCN Red List criteria. However, with most remaining habitat occurring outside of protected areas, and with habitat loss being exacerbated by climate change, a viable conservation plan remains elusive.     

A range wide geographic pattern of genetic diversity and population structure of Hexinia polydichotoma (Asteraceae) in Tarim Basin and adjacent areas

In order to investigate the genetic diversity and influence of climate oscillations on evolutionary processes of organisms in northwest China, we selected Hexinia polydichotoma, a species endemic to China, and examined the geographic pattern of genetic variation in its entire cover range, Tarim Basin and adjacent areas. In the study, 22 chloroplast (cp) haplotypes were identified based on two cpDNA sequences (trnH–psbA and ycf6–psbM), and five ITS sequence variations were found. Shown in the cp haplotype network, the two common cp haplotypes mainly distributed along the northern and southern rims of the basin respectively and intersected in the center of the basin, whereas in the ITS network, ITS genotype 1 was widespread across the whole distribution area, and rare genotypes were concentrated in the western rim of the basin. Genetic variation primarily occurred among populations and SAMOVA groups. Fragmented desert habitat may have caused gene flow barrier among populations or groups far from each other, leading to significant genetic differentiation at these levels. It appears that expansion and contraction cycles of river systems and oases during the middle Pleistocene was the source of the observed fragmentation. Geographic range expansion in H. polydichotoma was supported by the significant value for Tajima's D, Fu's F_S, and by a unimodal mismatch distribution. It was possible that the enlargement of the Taklimakan Desert during the middle Pleistocene may have provided appropriate conditions for the range dispersal. We identified western rim of the basin as the center of genetic diversity of H. polydichotoma based on the present dataset.     

The effect of past climatic oscillations on spatial genetic structure of Atraphaxis manshurica (Polygonoideae) in the Horqin sandlands,northern China

To investigate the effect of past climate oscillations on the genetic diversity and demographic history of organisms in the Horqin sandlands, at the margin of the monsoonal zone in northern China. We selected Atraphaxis manshurica as a model and examined the spatial genetic structure in populations of this species across most of its covered range. Five haplotypes were identified on the basis of two chloroplast DNA sequences (psbK-psbI and psbB-psbH) in 127 individuals from 11 natural populations. The total genetic diversity (H_T = 0.437) of A. manshurica is low compared with species to the west in arid northwestern China, outside of monsoonal influence. Analysis of molecular variance (AMOVA) indicates that genetic differentiation primarily occurred among two geographical groups distinguished by spatial analysis of molecular variance (SAMOVA). The BEAST molecular clock approach revealed that genetic divergence of the species mostly occurred in middle-late Pleistocene, in accordance with the periodic glacial periods and accompanying monsoonal changes. From species distribution modeling (SDM), we found that the species experienced range contraction and southward retreat during the Last Glacial Maximum (LGM). The current genetic pattern and demographic history may have thus been shaped by glacial–interglacial cycles and changes of habitat since the middle Pleistocene.