Phylogeographic patterns of the desert poplar in Northwest China shaped by both geology and climatic oscillations

Phylogeographic analysis reveals multiple origins of the desert shrub Reaumuria songarica in northern Xinjiang, involving homoploid and tetraploid hybrids

Hybrid speciation plays an important role in species diversification. The establishment of reproductive isolation is crucial for hybrid speciation, and the identification of diverse types of hybrids, particularly homoploid hybrid species, contributes to a comprehensive understanding of this process. Reaumuria songarica is a constructive shrub widespread in arid Central Asia. Previous studies have inferred that the R. songarica populations in the Gurbantunggut Desert (GuD) originated from homoploid hybridizations between its eastern and western lineages and may have evolved into an incipient species. To further elucidate the genetic composition of different hybrid populations and to determine the species boundary of this hybrid lineage, we investigated the overall phylogeographic structure of R. songarica based on variation patterns of five cpDNA and one nrITS sequences across 32 populations. Phylogenetic analyses demonstrated that within the GuD lineage, the Wuerhe population evolved directly from ancestral lineages, whereas the others originated from hybridizations between the eastern and western lineages. PCoA and genetic barrier analysis supported the subdivision of the GuD lineage into the southern (GuD-S) and northern (GuD-N) groups. Populations in the GuD-S group had a consistent genetic composition and the same ancestral female parent, indicating that they belonged to a homoploid hybrid lineage. However, the GuD-N group experienced genetic admixture of the eastern and western lineages on nrITS and cpDNA, with some populations inferred to be allopolyploid based on ploidy data. Based on cpDNA haplotypes, BEAST analyses showed that the GuD-S and GuD-N groups originated after 0.5 Ma. Our results suggest that multiple expansions and contractions of GuD, driven by Quaternary climatic oscillations and the Kunlun-Yellow River tectonic movement, are important causes of the complex origins of R. songarica populations in northern Xinjiang. This study highlights the complex origins of the Junggar Basin flora and the underappreciated role of hybridization in increasing its species diversity.

The Phylogeography of Deciduous Tree Ulmus macrocarpa (Ulmaceae) in Northern China

Disentangling how climate oscillations and geographical events significantly influence plants' genetic architecture and demographic history is a central topic in phytogeography. The deciduous ancient tree species Ulmus macrocarpa is primarily distributed throughout Northern China and has timber and horticultural value. In the current study, we studied the phylogenic architecture and demographical history of U. macrocarpa using chloroplast DNA with ecological niche modeling. The results indicated that the populations' genetic differentiation coefficient (NST) value was significantly greater than the haplotype frequency (GST) (p < 0.05), suggesting that U. macrocarpa had a clear phylogeographical structure. Phylogenetic inference showed that the putative chloroplast haplotypes could be divided into three groups, in which the group Ⅰ was considered to be ancestral. Despite significant genetic differentiation among these groups, gene flow was detected. The common ancestor of all haplotypes was inferred to originate in the middle-late Miocene, followed by the haplotype overwhelming diversification that occurred in the Quaternary. Combined with demography pattern and ecological niche modeling, we speculated that the surrounding areas of Shanxi and Inner Mongolia were potential refugia for U. macrocarpa during the glacial period in Northern China. Our results illuminated the demography pattern of U. macrocarpa and provided clues and references for further population genetics investigations of precious tree species distributed in Northern China.

Diversification and historical demography of Haloxylon ammodendron in relation to Pleistocene climatic oscillations in northwestern China

The influence of aridification and climatic oscillations on the genetic diversity and evolutionary processes of organisms during the Quaternary in northwestern China is examined using Haloxylon ammodendron. Based on the variation of two cpDNA regions (trnS-trnG and trnV) and one nDNA sequence (ITS1-ITS4) in 420 individuals from 36 populations, the spatial genetic structure and demographic history of H. ammodendron in arid China is examined. Median-joining network and Bayesian inference trees enabled the identification of three diverged lineages within H. ammodendron from 24 different haplotypes and 16 ribotypes, distributed across western (Xinjiang), eastern (Gansu and Inner Mongolia) and southern (Qinghai) regions. AMOVA analysis demonstrated that more than 80% of observed genetic variation related to lineage split was based on cpDNA and nDNA variation. Allopatric divergence among the three groups was mainly triggered by geographical isolation due to Xingxingxia rock and uplift of the Qilian Mountains during the Quaternary. Local adaptive differentiation among western, eastern and southern groups occurred due to gene flow obstruction resulting from arid landscape fragmentation accompanied by local environmental heterogeneity of different geographical populations. The southern margin of the Junggar Basin and the Tengger Desert possibly served as two independent glacial refugia for H. ammodendron. The distribution of genetic variation, coupled with SDMs and LCP results, indicated that H. ammodendron probably moved northward along the Junggar Basin and westward along Tengger Desert at the end of the last glacial maximum; postglacial re-colonization was probably westward and southward along the Hexi Corridor.

Phylogeography and Population Genetics Analyses Reveal Evolutionary History of the Desert Resource Plant Lycium ruthenicum (Solanaceae)

Climactic oscillations during the Quaternary played a significant role in the formation of genetic diversity and historical demography of numerous plant species in northwestern China. In this study, we used 11 simple sequence repeats derived from expressed sequence tag (EST-SSR), two chloroplast DNA (cpDNA) fragments, and ecological niche modeling (ENM) to investigate the population structure and the phylogeographic history of Lycium ruthenicum, a plant species adapted to the climate in northwestern China. We identified 20 chloroplast haplotypes of which two were dominant and widely distributed in almost all populations. The species has high haplotype diversity and low nucleotide diversity based on the cpDNA data. The EST-SSR results showed a high percentage of total genetic variation within populations. Both the cpDNA and EST-SSR results indicated no significant differentiation among populations. By combining the evidence from ENM and demographic analysis, we confirmed that both the last interglacial (LIG) and late-glacial maximum (LGM) climatic fluctuations, aridification might have substantially narrowed the distribution range of this desert species, the southern parts of the Junggar Basin, the Tarim Basin, and the eastern Pamir Plateau were the potential glacial refugia for L. ruthenicum during the late middle Pleistocene to late Pleistocene Period. During the early Holocene, the warm, and humid climate promoted its demographic expansion in northwestern China. This work may provide new insights into the mechanism of formation of plant diversity in this arid region.

Molecular Phylogeography and Intraspecific Divergences in Siberian Wildrye (Elymus sibiricus L.) Wild Populations in China, Inferred From Chloroplast DNA Sequence and cpSSR Markers

A detailed understanding of the distribution and degree of genetic variation within a species is important for determining their evolutionary potential, which in return facilitates the development of efficient conservation strategies aimed at preserving adaptive genetic variation. As an important perennial, cool-season grass in temperate Eurasia, increasing attention has been paid to Siberian wildrye (Elymus sibiricus) due to its excellent ecological utilization value and forage production potential in China, particularly in the Qinghai-Tibet Plateau (QTP) regions. In this study, we applied two chloroplast (cp) genes (matK and rbcL), three cp spacer regions (trnY-GUA∼trnD-GUC, atpH∼atpF, and rps4∼trnT-UGU), and six cpSSR markers to the genetic and phylogenetic analysis of 137 wild E. sibiricus accessions from 23 natural populations that represent the main distribution regions in China. The results show the highest genetic diversity (h = 0.913) and haplotype richness (10 haplotypes) for the QTP population, which indicates QTP as the probable diversity center and geographic origin of E. sibiricus in China. Population divergence was high, indicating a significant phylogeographic structure together with a significantly higher N_st value (N_st > G_st, P < 0.05) at the species level, QTP+XJ (combined populations from QTP and Xinjiang), QTP+NC (combined populations from QTP and North China), and XJ+NC (combined populations from Xinjiang and North China) group levels, respectively. An expansion was revealed in the distributional range of E. sibiricus in China from paleo times up to the recent past, while a dramatic range of contraction was predicted for the near future. The predicted main limiting factor for the further spread of E. sibiricus is an increasing global mean temperature. We recommend that the combination of Es-cpDNA1 and Es-cpDNA3+4+5 can be used as effective markers for phylogenetic analysis and phylogeographical history analysis of E. sibiricus. These findings shed new light on the historical population dynamics of cold-season herbs in the QTP region and the north of China and are of great significance for the future establishment of protection and collection strategies for wild E. sibiricus germplasm.

Genetic Diversity Study on Geographical Populations of the Multipurpose Species Elsholtzia stauntonii Using Transferable Microsatellite Markers

Elsholtzia stauntonii Benth. (Lamiaceae) is an economically important ornamental, medicinal and aromatic plant species. To meet the increasing market demand for E. stauntonii, it is necessary to assess genetic diversity within the species to accelerate the process of genetic improvement. Analysis of the transferability of simple sequence repeat (SSR) markers from related species or genera is a fast and economical method to evaluate diversity, and can ensure the availability of molecular markers in crops with limited genomic resources. In this study, the cross-genera transferability of 497 SSR markers selected from other members of the Lamiaceae (Salvia L., Perilla L., Mentha L., Hyptis Jacq., Leonurus L., Pogostemon Desf., Rosmarinus L., and Scutella L.) to E. stauntonii was 9.05% (45 primers). Among the 45 transferable markers, 10 markers revealed relatively high polymorphism in E. stauntonii. The genetic variation among 825 individuals from 18 natural populations of E. stauntonii in Hebei Province of China was analyzed using the 10 polymorphic SSR markers. On the basis of the SSR data, the average number of alleles (N _A), expected heterozygosity (H _E), and Shannon's information index (I) of the 10 primers pairs were 7.000, 0.478, and 0.688, respectively. Lower gene flow (N _m = 1.252) and high genetic differentiation (F _st = 0.181) were detected in the populations. Analysis of molecular variance (AMOVA) revealed that most of the variation (81.47%) was within the populations. Integrating the results of STRUCTURE, UPGMA (Unweighted Pair Group Method with Arithmetic Mean) clustering, and principal coordinate analysis, the 825 samples were grouped into two clusters associated with geographical provenance (southwestern and northeastern regions), which was consistent with the results of a Mantel test (r = 0.56, p < 0.001). Overall, SSR markers developed in related genera were effective to study the genetic structure and genetic diversity in geographical populations of E. stauntonii. The results provide a theoretical basis for conservation of genetic resources, genetic improvement, and construction of a core collection for E. stauntonii.

Insights into plant biodiversity conservation in large river valleys in China: A spatial analysis of species and phylogenetic diversity

Large river valleys (LRVs) are heterogeneous in habitat and rich in biodiversity, but they are largely overlooked in policies that prioritize conservation. Here, we aimed to identify plant diversity hotspots along LRVs based on species richness and spatial phylogenetics, evaluate current conservation effectiveness, determine gaps in the conservation networks, and offer suggestions for prioritizing conservation. We divided the study region into 50 km × 50 km grid cells and determined the distribution patterns of seed plants by studying 124,927 occurrence points belonging to 14,481 species, using different algorithms. We generated phylogenies for the plants using the "V. PhyloMaker" R package, determined spatial phylogenetics, and conducted correlation analyses between different distribution patterns and spatial phylogenetics. We evaluated the effectiveness of current conservation practices and discovered gaps of hotspots within the conservation networks. In the process, we identified 36 grid cells as hotspots (covering 10% of the total area) that contained 83.4% of the species. Fifty-eight percent of the hotspot area falls under the protection of national nature reserves (NNRs) and 83% falls under national and provincial nature reserves (NRs), with 42% of the area identified as conservation gaps of NNRs and 17% of the area as gaps of NRs. The hotspots contained high proportions of endemic and threatened species, as did conservation gaps. Therefore, it is necessary to optimize the layout of current conservation networks, establish micro-nature reserves, conduct targeted conservation priority planning focused on specific plant groups, and promote conservation awareness. Our results show that the conservation of three hotspots in Southwest China, in particular, is likely to positively affect the protection of biodiversity in the LRVs, especially with the participation of the neighboring countries, India, Myanmar, and Laos.

Phylogeography of Amygdalus mongolica in relation to Quaternary climatic aridification and oscillations in northwestern China

Quaternary period geological events and climatic oscillations significantly affect the geographic structure and genetic diversity of species distribution in arid northwestern China. Amygdalus mongolica is a relict and endangered shrub that occurs primarily in arid areas of northwestern China. Based on variation patterns present at three cpDNA regions (psbK-psbI, trnL-trnF and trnV) and in one nDNA sequence (ITS1-ITS4) in 174 individuals representing 15 populations, the spatial genetic structure and demographic history of A. mongolica was examined across its entire geographic range. The 17 different haplotypes and 10 ribotypes showed two lineages, distributed across the Western (Mazong Mountains, Hexi Corridor, and Alxa Left Banner) and Eastern regions (Urad Houqi, Yinshan Mountains, Urad Zhongqi, and Daqing Mountains) according to the median-joining network and the BI (Bayesian inference) and ML (Maximum likelihood) trees. AMOVA analysis demonstrated that over 65% of the observed genetic variation was related to this lineage split. The expansions of the Ulanbuhe and Tengger deserts and the eastward extension of the Yinshan Mountains since the Quaternary period likely interrupted gene flow and triggered the observed divergence in the two allopatric regions; arid landscape fragmentation accompanied by local environmental heterogeneity further increased local adaptive differentiation between the Western and Eastern groups. Based on the evidence from phylogeographical patterns and the distribution of genetic variation, A. mongolica distributed in the eastern and western regions are speculated to have experienced eastward migration along the southern slopes of the Lang Mountains and westward migration along the margins of the Ulanbuhe and Tengger deserts to the Hexi Corridor, respectively. For setting a conservation management plan, it is recommended that the south slopes of the Lang Mountains and northern Helan Mountains be identified as the two primary conservation areas, as they have high genetic variation and habitats that are more suitable.

Evolutionary history of a desert perennial Arnebia szechenyi (Boraginaceae): Intraspecific divergence, regional expansion and asymmetric gene flow

The complex interactions of historical, geological and climatic events on plant evolution have been an important research focus for many years. However, the role of desert formation and expansion in shaping the genetic structures and demographic histories of plants occurring in arid areas has not been well explored. In the present study, we investigated the phylogeography of Arnebia szechenyi, a desert herb showing a near-circular distribution surrounding the Tengger Desert in Northwest China. We measured genetic diversity of populations using three maternally inherited chloroplast DNA (cpDNA) fragments and seven bi-paternally inherited nuclear DNA (nDNA) loci that were sequenced from individuals collected from 16 natural populations across its range and modelled current and historical potential habitats of the species. Our data indicated a considerably high level of genetic variation within A. szechenyi and noteworthy asymmetry in historical migration from the east to the west. Moreover, two nuclear genetic groups of populations were revealed, corresponding to the two geographic regions separated by the Tengger Desert. However, analysis of cpDNA data did not show significant geographic structure. The most plausible explanation for the discrepancy between our findings based on cpDNA and nDNA data is that A. szechenyi populations experienced long periods of geographic isolation followed by range expansion, which would have promoted generalized recombination of the nuclear genome. Our findings further highlight the important role that the Tengger Desert, together with the Helan Mountains, has played in the evolution of desert plants and the preservation of biodiversity in arid Northwest China.

Phylogeography Reveals Geographic and Environmental Factors Driving Genetic Differentiation of Populus sect. Turanga in Northwest China

Populus sect. Turanga (hereafter referred to as "Populus"), including Populus euphratica and Populus pruinosa, are the predominant tree species in desert riparian forests in northwestern China. These trees play key roles in maintaining ecosystem balance, curbing desertification, and protecting biodiversity. However, the distribution area of Populus forests has been severely diminished and degraded in recent years due to increased habitat destruction and human activity. Understanding the genetic diversity among Populus individuals and populations is essential for designing conservation strategies, but comprehensive studies of their genetic diversity in northwest China are lacking. Here, we assessed the population structures and genetic diversity of 1,620 samples from 85 natural populations of Populus (59 P. euphratica and 26 P. pruinosa populations) covering all of northwestern China using 120 single nucleotide polymorphism (SNP) markers. Analysis of population structure revealed significant differentiation between these two sister species and indicated that strong geographical distribution patterns, a geographical barrier, and environmental heterogeneity shaped the extant genetic patterns of Populus. Both P. euphratica and P. pruinosa populations in southern Xinjiang had higher genetic diversity than populations in other clades, perhaps contributing to local geographic structure and strong gene flow. Analysis of molecular variance (AMOVA) identified 15% variance among and 85% variance within subpopulations. Mantel tests suggested that the genetic variation among P. euphratica and P. pruinosa populations could be explained by both geographical and environmental distance. The genetic diversity of P. euphratica showed a significant negative correlation with latitude and longitude and a positive correlation with various environmental factors, such as precipitation of warmest quarter and driest month, temperature seasonality, and annual mean temperature. These findings provide insights into how the genetic differentiation of endangered Populus species was driven by geographical and environmental factors, which should be helpful for designing strategies to protect these genetic resources in the future.

Demographic history and genetic differentiation of an endemic and endangered Ulmus lamellosa (Ulmus)

BACKGROUND

Ulmus lamellosa (one of the ancient species of Ulmus) is an endemic and endangered plant that has undergone climatic oscillations and geographical changes. The elucidation of its demographic history and genetic differentiation is critical for understanding the evolutionary process and ecological adaption to forests in Northern China.

RESULTS

Polymorphic haplotypes were detected in most populations of U. lamellosa via DNA sequencing. All haplotypes were divided into three phylogeographic clades fundamentally corresponding to their geographical distribution, namely THM (Taihang Mountains), YM (Yinshan Mountains), and YSM (Yanshan Mountains) groups. The YSM group, which is regarded as ancestral, possessed higher genetic diversity and significant genetic variability in contrast to the YSM and YM groups. Meanwhile, the divergence time of intraspecies haplotypes occurred during the Miocene-Pliocene, which was associated with major Tertiary geological and/or climatic events. Different degrees of gene exchanges were identified between the three groups. During glaciation, the YSM and THM regions might have served as refugia for U. lamellosa. Based on ITS data, range expansion was not expected through evolutionary processes, except for the THM group. A series of mountain uplifts (e.g., Yanshan Mountains and Taihang Mountains) following the Miocene-Pliocene, and subsequently quaternary climatic oscillations in Northern China, further promoted divergence between U. lamellosa populations.

CONCLUSIONS

Geographical topology and climate change in Northern China played a critical role in establishing the current phylogeographic structural patterns of U. lamellosa. These results provide important data and clues that facilitate the demographic study of tree species in Northern China.

Enhancing the expression of ARK1 genes in poplar leads to multiple branches and transcriptomic changes

The ARBORKNOX1 (ARK1) gene is an important gene for regulating plant growth and development; however, transcriptomic responses of enhancing expression of ARK1 gene in poplar are poorly investigated. To provide insight into the gene function of the ARK1 gene in poplar, the ARK1 transgenic poplar '717' and '84 K' lines were obtained, the morphology of transgenic plants was observed, and transcriptome profiles were compared. The results showed that there were multiple branches in ARK1 transgenic seedlings compared with non-transgenic seedlings. The results of transcriptome analysis showed that there were significant differences in transcriptome profiles between the transgenic lines of '717' and '84 K', and between non-transgenic lines (CK) and transgenic plants. The real-time quantitative polymerase chain reaction (RT-qPCR) analysis confirmed the expression levels of the genes involved in the pathway of zeatin biosynthesis and brassinosteroid biosynthesis. The increase in expression levels of AHP and CYCD3 was related to multiple branches. Enhancing the expression of the ARK1 gene in poplar seedlings leads to multiple branches and transcriptomic changes.

Genomic Phylogeography of Gymnocarpos przewalskii (Caryophyllaceae): Insights into Habitat Fragmentation in Arid Northwestern China

Extensive range of deserts and gobis (rocks) had promoted habitat fragmentation of species in arid northwestern China. Distribution of endangered Gymnocarpos przewalskii Maxim. covers most of gobis (rocks) and desert terrain across arid regions of northwestern China. In the present study, we had employed genomic phylogeographical analysis to investigate population structure of G. przewalskii and test the effect of environmental conditions on spatial pattern of genetic diversity. Results showed four groups were identified from east to west: Edge of the Alxa Desert, Hexi Corridor, Hami Basin, and North edge of the Tarim Basin. Genetic diversity was at an equal level among four groups. General linear model (GLM) analysis showed spatial pattern of genetic diversity was significant correlated with three habitat variables including habitat suitability at present (Npre) and last glacial maximum (LGM) (NLGM) periods, and locality habitat stability (NStab). It concluded that habitat fragmentation had triggered lineage divergences of G. przewalskii in response to long-term aridification. Genome-wide single nucleotide polymorphisms (SNPs) could increase the ability of clarifying population structures in comparison with traditional molecular markers. Spatial pattern of genetic diversity was determined by fragmented habitats with high habitat suitability (Npre and NLGM) and stability (NStab). At last, we propose to establish four conservation units which are in consistent with the population grouping to maintain the genetic integrity of this endangered species.

The East Asian Winter Monsoon Acts as a Major Selective Factor in the Intraspecific Differentiation of Drought-Tolerant Nitraria tangutorum in Northwest China

The influence of Quaternary climate fluctuation on the geographical structure and genetic diversity of species distributed in the regions of the Qinghai–Tibet Plateau (QTP) has been well established. However, the underlying role of the East Asian monsoon system (EAMS) in shaping the genetic structure of the population and the demography of plants located in the arid northwest of China has not been explored. In the present study, Nitraria tangutorum, a drought-tolerant desert shrub that is distributed in the EAMS zone and has substantial ecological and economic value, was profiled to better understand the influence of EAMS evolution on its biogeographical patterns and demographic history. Thus, the phylogeographical structure and historical dynamics of this plant species were elucidated using its five chloroplast DNA (cpDNA) fragments. Hierarchical structure analysis revealed three distinct, divergent lineages: West, East-A, and East-B. The molecular dating was carried out using a Bayesian approach to estimate the time of intraspecies divergence. Notably, the eastern region, which included East-A and East-B lineages, was revealed to be the original center of distribution and was characterized by a high level of genetic diversity, with the intraspecific divergence time dated to be around 2.53 million years ago (Ma). These findings, combined with the data obtained by ecological niche modeling analysis, indicated that the East lineages have undergone population expansion and differentiation, which were closely correlated with the development of the EAMS, especially the East Asian winter monsoon (EAWM). The West lineage appears to have originated from the migration of N. tangutorum across the Hexi corridor at around 1.85 Ma, and subsequent colonization of the western region. These results suggest that the EAWM accelerated the population expansion of N. tangutorum and subsequent intraspecific differentiation. These findings collectively provide new information on the impact of the evolution of the EAMS on intraspecific diversification and population demography of drought-tolerant plant species in northwest China.

Genome resequencing reveals demographic history and genetic architecture of seed salinity tolerance in Populus euphratica

Populus euphratica is a dominant tree species in desert riparian forests and possesses extraordinary adaptation to salinity stress. Exploration of its genomic variation and molecular underpinning of salinity tolerance is important for elucidating population evolution and identifying stress-related genes. Here, we identify approximately 3.15 million single nucleotide polymorphisms using whole-genome resequencing. The natural populations of P. euphratica in northwest China are divided into four distinct clades that exhibit strong geographical distribution patterns. Pleistocene climatic fluctuations and tectonic deformation jointly shaped the extant genetic patterns. A seed germination rate-based salinity tolerance index was used to evaluate seed salinity tolerance of P. euphratica and a genome-wide association study was implemented. A total of 38 single nucleotide polymorphisms were associated with seed salinity tolerance and were located within or near 82 genes. Expression profiles showed that most of these genes were regulated under salt stress, revealing the genetic complexity of seed salinity tolerance. Furthermore, DEAD-box ATP-dependent RNA helicase 57 and one undescribed gene (CCG029559) were demonstrated to improve the seed salinity tolerance in transgenic Arabidopsis. These results provide new insights into the demographic history and genetic architecture of seed salinity tolerance in desert poplar.

Phylogeography and ecological niche modeling unravel the evolutionary history of the Yarkand hare, Lepus yarkandensis (Mammalia: Leporidae), through the Quaternary

Background

The Taklimakan Desert in China is characterized by unique geological and historical dynamics and endemic flora and fauna, but the influence of historical climate oscillations on the evolutionary history of endemic animals is poorly understood. Lepus yarkandensis is an oases-dependent Near Threatened species that lives in fragmented oasis habitats in the Taklimakan Desert, China. We investigated the geological and climatic impacts on its geographical differentiation, demographic history and influence of Pleistocene glacial-interglacial cycles on the evolutionary history of L. yarkandensis. Further, studied the impact of climatic oscillation based modification on phylogeography, distribution and diversification pattern of Yarkand hare by using Cytb (1140 bp), MGF (592 bp) and SPTBN1 (619 bp) markers. Ecological niche modeling (ENM) revealed the evolutionary history of this species in response to climate change during the Quaternary. Paleodistribution modeling was used to identify putative refugia and estimate their historical distributions.

Results

Both historical demographic analyses and climatic niche modeling revealed strong effects of glacial climate changes, suggesting recurrent range contractions and expansions. The EBSP results indicated clear population expansion of L. yarkandensis since the Pleistocene. In the “early Pleistocene”, the demographic expansion continued from 0.83 MYA to the last glacial period. The ENM analysis supported a wide distribution of Lepus yarkandensis at high altitudes during the last interglacial (LIG) period. During the last glacial maximum (LGM), the suitable climate was reduced and restricted to the western part of the Taklimakan Desert.

Conclusions

Inland aridification, oasis evolution and river flow played major roles in the population differentiation and demographic history of Yarkand hares. Historically, the large, continuous oases in the Taklimakan Desert contained a viable and unique population of L. yarkandensis. The fragmented desert environment might have caused low gene flow between individuals or groups, thus leading to predominant genetic differentiation. The Pleistocene climatic cycles triggered the diversification and expansion of this species during cold and warm periods, respectively, leading to multiple colonization events within the Taklimakan Desert. These events might be due to the expansion of the Taklimakan Desert during the Middle Pleistocene. Yarkand hare previously occupied vast areas at low and intermediate altitudes in Xinjiang, Gansu, Shanxi, Henan and Shaanxi Provinces in China. The past aridification, climate change-induced oasis modifications, changes in river volumes and flow directions, and human activities all affected the population demography and phylogeography of the Yarkand hare.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1426-z) contains supplementary material, which is available to authorized users.