Repeated range expansion and glacial endurance of Potentilla glabra (Rosaceae) in the Qinghai-Tibetan plateau

Conserved DNA sequence analysis reveals the phylogeography and evolutionary events of Akebia trifoliata in the region across the eastern edge of the Tibetan Plateau and subtropical China

BACKGROUND

The eastern edge of the Qinghai‒Tibet Plateau (QTP) and subtropical China have various regions where plant species originate and thrive, but these regions have been the focus of very few integrative studies. Here, we elucidated the phylogeographic structure of a continuous and widespread Akebia trifoliata population across these two regions.

RESULTS

Sixty-one populations consisting of 391 genotypes were examined to assess population diversity and structure via network distribution analysis, maximum likelihood phylogenetic tree reconstruction, divergence time estimation, demographic history inference, and ancestral area reconstruction of both conserved internal transcribed spacer (ITS) and chloroplast (rps16) DNA sequences. The results showed that the ITS region was more variable than the rps16 region and could be suitable for studying intraspecific phylogeography. The A. trifoliata population displayed high genetic diversity, genetic differentiation and obvious phylogeographical structure, possibly originating on the eastern QTP, expanding during the last glacial-interglacial cycle, diverging in the early Pleistocene and middle Pleistocene, and extensively migrating thereafter. The migration route from west to east along rivers could be largely responsible for the long-distance dispersal of this species, while three main refuges (Qinba Mountains, Nanling Mountains and Yunnan-Guizhou Plateau) with multiple ice shelters facilitated its wide distribution.

CONCLUSIONS

Our results suggested that the from west to east long migration accompanying with the minor short reciprocal migration in the south-north direction, and the three main refuges (the Qinba Mountains, Nanling Mountains and Yunnan-Guizhou Plateau) contributed to the extant geographical distribution of A. trifoliata. In addition, this finding also strongly reduced the discrepancy between glacial contraction and postglacial expansion and the in situ survival hypothesis by simultaneously considering the existence of many similar climate-related ecological niches and migration influences.

Phylogeography and ecological niche modeling implicate multiple microrefugia of Swertia tetraptera during quaternary glaciations

BACKGROUND

Climate fluctuations during the Pleistocene and mountain uplift are vital driving forces affecting geographic distribution. Here, we ask how an annual plant responded to the Pleistocene glacial cycles.

METHODS

In this study, we analyzed the population demographic history of the annual herb Swertia tetraptera Maxim (Gentianaceae) endemic to Qinghai-Tibetan Plateau (QTP). A total of 301 individuals from 35 populations of S. tetraptera were analyzed based on two maternally inherited chloroplast fragments (trnL-trnF and trnS-trnG). Phylogeographic analysis was combined with species distribution modeling to detect the genetic variations in S. tetraptera.

RESULTS

The genetic diversity of S. tetraptera was high, likely due to its wide natural range, high proportion of endemic haplotypes and evolutionary history. Fifty-four haplotypes were identified in S. tetraptera. Only a few haplotypes were widespread (Hap_4, Hap_1, Hap_3), which were dispersed throughout the present geographical range of S. tetraptera, while many haplotypes were confined to single populations. The cpDNA dataset showed that phylogeographic structuring was lacking across the distribution range of S. tetraptera. Analyses of molecular variance showed that most genetic variation was found within populations (70.51%). In addition, the relationships of the haplotypes were almost completely unresolved by phylogenetic reconstruction. Both mismatch distribution analysis and neutrality tests showed a recent expansion across the distribution range of S. tetraptera. The MaxEnt analysis showed that S. tetraptera had a narrow distribution range during the Last Glacial Maximum (LGM) and a wide distribution range during the current time, with predictions into the future showing the distribution range of S. tetraptera expanding.

CONCLUSION

Our study implies that the current geographic and genetic distribution of S. tetraptera is likely to have been shaped by Quaternary periods. Multiple microrefugia of S. tetraptera existed during Quaternary glaciations. Rapid intraspecific diversification and hybridization and/or introgression may have played a vital role in shaping the current distribution patterns of S. tetraptera. The distribution range of S. tetraptera appeared to have experienced contraction during the LGM; in the future, when the global climate becomes warmer with rising carbon dioxide levels, the distribution of S. tetraptera will expand.

Phylogeography of the Plateau Pika (Ochotona curzoniae) in Response to the Uplift of the Qinghai-Tibet Plateau

The evolution and current distribution of species on the Qinghai-Tibet Plateau have been significantly impacted by historical occurrences, including the uplift of the plateau and the Quaternary climate upheaval. As a remnant species, the plateau pika (Ochotona curzoniae) is a great model for researching historical events. In this study, 302 samples from 42 sample sites were utilized to analyze the impact of historical events on the evolution and distribution pattern of plateau pikas. The genetic diversity, patterns of differentiation, and historical dynamics of the plateau pika were investigated using molecular markers that included four mitochondrial genes (COI, D-loop, Cytb, and 12S rRNA) and three nuclear genes (GHR, IRBP, and RAG1). The results showed that: (1) The genetic diversity of the plateau pika was high in the Tibetan Plateau (Hd = 0.9997, π = 0.01205), and the plateau pika evolved into five lineages that occupied different geographical areas, with lineage 1 (Group 1) in the south of the Yarlung Zangbo River, lineage 2 (Group 2) in the hinterland of the plateau, lineage 3 (Group 3) in the northeastern part of the plateau, lineage 4 (Group 4) in the Hengduan Mountains, and lineage 5 (Group 5) in the eastern part of the plateau. (2) The gene flow among the five lineages was low, and the differentiation level was high (Nm < 0.25; Fst > 0.25), indicating that the geographical barriers between the five lineages, such as the Yarlung Zangbo River, the Qaidam-Ghuong-Guide Basin, and the Lancang River, effectively promoted the population differentiation of the plateau pika. (3) The plateau pika first spread from the Hengduan Mountains to the entire Qinghai-Tibet Plateau and then conducted small-scale migration and dispersal in several refuges across the plateau in response to climate changes during the glacial and interglacial periods. (4) Except for Group 1 and Group 4, all the other populations exhibited a rapid expansion between 0.06 and 0.01 Mya, but the expansion was considerably delayed or halted by the effects of climate change during the last glacial maximum (0.02 Mya). Overall, the plateau pika on the Qinghai-Tibet Plateau exhibits high genetic diversity, and topographic obstacles, including mountains, valleys, and basins, created by the uplift of the plateau and climatic changes since the Quaternary period have played an important role in the differentiation and historical dynamics of the plateau pika population.

Genetic diversity and population dynamic of Ziziphus jujuba var. spinosa (Bunge) Hu ex H. F. Chow in Central China

Phylogeographic research concerning Central China has been rarely conducted. Population genetic and phylogeography of Ziziphus jujuba var. spinosa (also called sour jujube) were investigated to improve our understanding of plant phylogeographic patterns in Central China. Single-copy nuclear gene markers and complete chloroplast genome data were applied to 328 individuals collected from 21 natural populations of sour jujube in China. Nucleotide variation of sour jujube was relatively high (π = 0.00720, θ _w = 0.00925), which resulted from the mating system and complex population dynamics. Analysis of molecular variation analysis revealed that most of the total variation was attributed to variation within populations, and a high level of genetic differentiation among populations was detected (F _st = 0.197). Relatively low long-distance dispersal capability and vitality of pollen contributed to high genetic differentiation among populations. Differences in the environmental conditions and long distance among populations further restricted gene flow. Structure clustering analysis uncovered intraspecific divergence between central and marginal populations. Migrate analysis found a high level of gene flow between these two intraspecific groups. Bayesian skyline plot detected population expansion of these two intraspecific groups. Network and phylogeny analysis of chloroplast haplotypes also found intraspecific divergence, and the divergence time was estimated to occur at about 55.86 Ma. Haplotype native to the Loess Plateau was more ancient, and multiple glacial refugia of sour jujube were found to locate at the Loess Plateau, areas adjacent to the Qinling Mountains and Tianmu Mountains. Species distribution model analysis found a typical contraction-expansion model corresponding to the Quaternary climatic oscillations. In the future, the distribution of sour jujube may shift to high-latitude areas. This study provides new insights for phylogeographic research of temperate plant species distributed in Central China and sets a solid foundation for the application of the scientific management strategy of Z. jujuba var. spinosa.

Dispersal into the Qinghai-Tibet plateau: evidence from the genetic structure and demography of the alpine plant Triosteum pinnatifidum

Triosteum pinnatifidum Maxim., an alpine plant, is traditionally used for several medicinal purposes. Here, both chloroplast DNA sequences and nuclear low copy sequence markers were used to investigate the genetic diversity and population structure of T. pinnatifidum. Materials were collected from thirteen localities in the northeast Qinghai-Tibet Plateau (QTP) and adjacent highlands and advanced analytical toolkits were used to access their origin and range shifts. The results revealed a higher level of population differentiation based on chloroplast DNA (cpDNA) concatenated sequences compared with the nuclear DNA sequences (F _ST = 0.654 for cpDNA, F _ST = 0.398 for AT103), indicating that pollen flow was still extensive in T. pinnatifidum. A decline in haplotype variation was observed from the plateau edge and adjoining highlands toward the platform of the QTP. The hypothesis "dispersal into the QTP," proposing that T. pinnatifidum experienced migration from the plateau edge and adjacent highlands to the platform, was supported. These results were in line with the hypothesis that multiple refugia exist on the plateau edge and adjacent highlands rather than on the plateau platform. Our unimodal mismatch distribution, star-like network supported a recent expansion in T. pinnatifidum.

Phylogeography of Iris loczyi (Iridaceae) in Qinghai-Tibet Plateau revealed by chloroplast DNA and microsatellite markers

Quaternary climate oscillations and complex topography have tremendous effects on current distribution and genetic structure of species, and hence the Qinghai-Tibet Plateau (QTP), the largest plateau in the world, has become a hotspot for many phylogeographic studies. However, little is known about the phylogeographic pattern of herbaceous plants in QTP. Here, we investigate the genetic diversity, population structure and historical dynamics of Iris loczyi, using five chloroplast DNA (cpDNA) fragments and seven microsatellite markers. A total of 15 populations, and 149 individuals were sampled throughout the QTP. High genetic diversity was detected both in cpDNA (H _d = 0.820) and SSR (H _o = 0.689, H _e = 0.699). Ten cpDNA haplotypes and 163 alleles were identified. AMOVA and clustering analyses revealed obvious differentiation between regions. The N _st, G _st and Mantel test showed significant phylogeographic structure of I. loczyi. The neutrality test and mismatch distribution analyses indicated that I. loczyi could not have undergone a historical population expansion, but population XS from the Qilian Mountain area could have experienced a local expansion. Bottleneck analyses indicated that I. loczyi had not experienced bottleneck recently. Based on cpDNA and SSR results, the Qilian Mountain area was inferred as a potential glacial refuge, and the southern Tibet valley was considered as a 'microrefugia' for I. loczyi. These findings provided new insights into the location of glacial refuges for the species distributed in QTP, and supplemented more plant species data for the response of QTP species to the Quaternary climate.

Spatiotemporal maintenance of flora in the Himalaya biodiversity hotspot: Current knowledge and future perspectives

Mountain ecosystems support a significant one-third of all terrestrial biodiversity, but our understanding of the spatiotemporal maintenance of this high biodiversity remains poor, or at best controversial. The Himalaya hosts a complex mountain ecosystem with high topographic and climatic heterogeneity and harbors one of the world's richest floras. The high species endemism, together with increasing anthropogenic threats, has qualified the Himalaya as one of the most significant global biodiversity hotspots. The topographic and climatic complexity of the Himalaya makes it an ideal natural laboratory for studying the mechanisms of floral exchange, diversification, and spatiotemporal distributions. Here, we review literature pertaining to the Himalaya in order to generate a concise synthesis of the origin, distribution, and climate change responses of the Himalayan flora. We found that the Himalaya supports a rich biodiversity and that the Hengduan Mountains supplied the majority of the Himalayan floral elements, which subsequently diversified from the late Miocene onward, to create today's relatively high endemicity in the Himalaya. Further, we uncover links between this Miocene diversification and the joint effect of geological and climatic upheavals in the Himalaya. There is marked variance regarding species dispersal, elevational gradients, and impact of climate change among plant species in the Himalaya, and our review highlights some of the general trends and recent advances on these aspects. Finally, we provide some recommendations for conservation planning and future research. Our work could be useful in guiding future research in this important ecosystem and will also provide new insights into the maintenance mechanisms underpinning other mountain systems.

Effects of Mountain Uplift and Climatic Oscillations on Phylogeography and Species Divergence of Chamaesium (Apiaceae)

Exploring the effects of orographic events and climatic shifts on the geographic distribution of organisms in the Himalayas-Hengduan Mountains (HHM) region and Qinghai-Tibetan Plateau (QTP) is crucial to understand the impact of environmental changes on organism evolution. To gain further insight into these processes, we reconstructed the evolutionary history of nine Chamaesium species distributed across the HHM and QTP regions. In total, 525 individuals from 56 populations of the nine species were analyzed based on three maternally inherited chloroplast fragments (rpl16, trnT-trnL, and trnQ-rps16) and one nuclear DNA region (internal transcribed spacer, ITS). Fifty-two chloroplast DNA (cpDNA) and 47 ITS haplotypes were identified in nine species. All of the cpDNA and ITS haplotypes were species-specific. Phylogenetic analysis suggested that all nine species form a monophyletic clade with high support. Dating analysis and ancestral area reconstruction revealed that the ancestral group of Chamaesium originated in the southern Himalayan region at the beginning of the Paleogene (60.85 Ma). The nine species of Chamaesium then separated well during the last 25 million years started in Miocene. Our maxent modeling indicated the broad-scale distributions of all nine species remained fairly stable from LIG to the present and predicted that it will remain stable into the future. The initial split of Chamaesium was triggered by climate changes following the collision of the Indian plate with the Eurasia plate during the Eocene. Subsequently, divergences within Chamaesium may have been induced by the intense uplift of the QTP, the onset of the monsoon system, and Central Asian aridification. Long evolutionary history, sexual reproduction, and habitat fragmentation could contribute to the high level of genetic diversity of Chamaesium. The higher genetic differentiation among Chamaesium populations may be related to the drastic changes of the external environment in this region and limited seed/pollen dispersal capacity.

Phylogeographic analyses of poplar revealed potential glacial refugia and allopatric divergence in southwest China

The distribution pattern and genetic structure of plant species have been profoundly influenced by climate oscillations. Phylogeographic analyses have been numerously conducted in biodiversity hotspot regions and some general conclusions have been documented. However, other regions have received less attentions and these places may serve as potential glacial refugia for plant species to survive the Quaternary glaciation. Here, we used six nuclear and three cpDNA markers to estimate the phylogeographic pattern of Populus davidiana, a widespread species distributed in southwest China. As a widely distributed species in southwest China, the nucleotide diversity of P. davidiana was relatively high (N _a = 6.28, H _O = 0.534, and H _E = 0.658). Genetic differentiation (F _ST) between the two main distribution regions, Yunnan and Guizhou provinces, was 0.21221. According to the composition of chloroplast haplotypes and the result of structure in these populations, we clearly distinguished two distantly sublineages corresponding to two distribution regions. Results of the Mantel test showed that there was a significant correlation between genetic distance and geographical distance (R ² = 0.8252, p<.05). The topographically heterogeneous regions and the low dispersal ability of seed and pollen may lead to high genetic differentiation between these two regions. A potential glacial refugia for P. davidiana located in adjacent regions to the Hengduan range was revealed and allopatric divergence in separated glacial refugia may directly lead to the present phylogeographic pattern of this species.

Evolutionary history of endangered and relict tree species Dipteronia sinensis in response to geological and climatic events in the Qinling Mountains and adjacent areas

Geological and climatic events are considered to profoundly affect the evolution and lineage divergence of plant species. However, the evolutionary histories of tree species that have responded to past geological and climate oscillations in central China's mountainous areas remain mostly unknown. In this study, we assessed the evolutionary history of the endangered and relict tree species Dipteronia sinensis in the Qinling Mountains (QM) and adjacent areas in East Asia based on variations in the complete chloroplast genomes (cpDNA) and reduced-genomic scale single nucleotide polymorphisms (SNPs). Population structure and phylogenetic analysis based on the cpDNA variations suggested that D. sinensis could be divided into two intraspecific genetic lineages in the eastern and western sides of the QM (EQM and WQM, respectively) in East Asia. Molecular dating suggested that the intraspecific divergence of D. sinensis occurred approximately 39.2 million years ago during the later Paleogene. It was significantly correlated with the orogeny of the QM, where the formation of this significant geographic barrier in the region may have led to the divergence of independent lineages. Bayesian clustering and demographic analysis showed that intraspecific gene flow was restricted between the EQM and WQM lineages. Isolation-with-migration analysis indicated that the two genetic lineages experienced significant demographic expansions after the Pleistocene ice ages. However, the genetic admixture was determined in some populations between the two lineages by the large scale of SNP variations due to DNA incompatibility, the large significant population size, and rapid gene flow of nuclear DNA markers. Our results suggest that two different conservation and management units should be constructed for D. sinensis in the EQM and WQM areas. These findings provide novel insights into the unprecedented effects of tectonic changes and climatic oscillations on lineage divergence and plant population evolution in the QM and adjacent areas in East Asia.

Population Genetic Structure and Demographic History of Primula fasciculata in Southwest China

Understanding the factors that drive the genetic structure of a species and its responses to past climatic changes is an important first step in modern population management. The response to the last glacial maximum (LGM) has been well studied, however, the effect of previous glaciation periods on plant demographic history is still not well studied. Here we investigated the population structure and demographic history of Primula fasciculata that widely occurs in the Hengduan Mountains and Qinghai-Tibetan Plateau. We obtained genomic data for 234 samples of the species using restriction site-associated DNA (RAD) sequencing and combined approximate Bayesian computation (ABC) and species distribution modeling (SDM) to evaluate the effects of multiple glaciation periods by testing several population divergence models and demographic scenarios. The analyses of population structure showed that P. fasciculata displays a striking population structure with six groups that could be identified genetically. Our ABC modeling suggested that the current groups diverged from ancestral populations located in the eastern Hengduan Mountains after the largest glaciation occurred in the region (~ 0.8-0.5 million years ago), which is consistent with the result of SDMs. Each current group has survived in different glacial refugia during the LGM and experienced expansions and/or bottlenecks since their divergence during or across the following Quaternary glacial cycles. Our study demonstrates the usefulness of population genomics for evaluating the effects of past climatic changes in alpine plant species with shallow population structure.

Molecular Phylogeography and Evolutionary History of the Endemic Species Corydalis hendersonii (Papaveraceae) on the Tibetan Plateau Inferred From Chloroplast DNA and ITS Sequence Variation

In response to past climatic changes, the species with different habits or adaptive traits likely have experienced very different evolutionary histories, especially for species that restricted to high mountain areas. In order to trace how Quaternary climatic oscillations affected range distributions and intraspecific divergence of such alpine plants on the Tibetan Plateau, here, we investigated maternally inherited chloroplast DNA (cpDNA) markers and biparentally inherited nuclear ribosomal internal transcribed spacer (ITS) DNA variations and aimed to explore the phylogeographic history of the endemic alpine species Corydalis hendersonii Hemsl. (Papaveraceae). We sequenced four cpDNA fragments (trnS-trnG, trnT-trnL, atpH-atpI, and psbE-petL) and also the nuclear (ITS) region in 368 individuals from 30 populations across the species' range. The network and phylogenetic analysis based on cpDNA variations identified 15 chlorotypes that cluster into three distinct clades. However, our nuclear DNA results demonstrated that there were four genetic/geographical groups within C. hendersonii. Some common and highly divergent cpDNA and ITS haplotypes were distributed in the populations of central and northeastern Tibetan Plateau, and the highest nucleotide diversity and genetic differentiation were detected in the central region. Demographic tests further indicated that the populations of southwestern and western Tibet may have experienced recent range expansion, which most likely occurred during the last glacial maximum (LGM) and continued its expansion after the beginning of the Holocene. These two different groups of this species may have derived from potential refugia that existed in the central and/or northeastern regions of Tibet during recent interglacial periods. In addition, our AMOVA analyses detected high genetic differentiation along with the whole sampling range. Also, distinct phylogeographic structures were detected among populations of C. hendersonii based on both of cpDNA and ITS variation. These findings shed new light on the importance of climatic oscillations during Quaternary and complex local topography as causes of intraspecific diversification and demographic changes within cold-tolerant herbs in the Tibetan Plateau biodiversity hotspot.

Complex population evolutionary history of four cold-tolerant Notopterygium herb species in the Qinghai-Tibetan Plateau and adjacent areas

Historical geological and climatic events are the most important drivers of population expansions/contractions, range shifts, and interspecific divergence in plants. However, the species divergence and spatiotemporal population dynamics of alpine cold-tolerant herbal plants in the high-altitude Qinghai-Tibetan Plateau (QTP) and adjacent areas remain poorly understood. In this study, we investigated population evolutionary history of four endangered Notopterygium herb species in the QTP and adjacent regions. We sequenced 10 nuclear loci, 2 mitochondrial DNA regions, and 4 chloroplast DNA regions in a total of 72 natural populations from the 4 species, and tested the hypothesis that the population history of these alpine herbs was markedly affected by the Miocene-Pliocene QTP uplifts and Quaternary climatic oscillations. We found that the four Notopterygium species had generally low levels of nucleotide variability within populations. Molecular dating and isolation-with-migration analyses suggested that Notopterygium species diverged ~1.74-7.82 million years ago and their differentiation was significantly associated with recent uplifts of the eastern margin of the QTP. In addition, ecological niche modeling and population history analysis showed that N. incisum and N. franchetii underwent considerable demographic expansions during the last glacial period of the Pleistocene, whereas a demographic contraction and a expansion occurred for N. forrestii and N. oviforme during the antepenultimate interglacial period and penultimate glacial period, respectively. These findings highlight the importance of geological and climatic changes during the Miocene-Pliocene and Pleistocene as causes of species divergence and changes in population structure within cold-tolerant herbs in the QTP biodiversity hotspot.

Quaternary climate change and habitat preference shaped the genetic differentiation and phylogeography of Rhodiola sect. Prainia in the southern Qinghai-Tibetan Plateau

There are two long-standing biogeographic hypotheses regarding the glacial survival of plant species in the Qinghai-Tibetan Plateau (QTP): the in situ survival hypothesis and the tabula rasa hypothesis. We tested these two hypotheses in a phylogeographic study of Rhodiola sect. Prainia, a monophyletic section with ecologically divergent lineages. Molecular data from the nuclear internal transcribed spacer, six plastid markers and 13 nuclear microsatellite loci were analyzed for 240 individuals from 19 populations of this section. Environmental data were used to analyze the niches of major phylogenetic lineages within this section and to model changes in their distributions since the Last Glacial Maximum (LGM). We found that Rhodiola sect. Prainia consists of three evolutionary lineages: all populations of R. stapfii, R. prainii populations at the southern edge of the QTP, and R. prainii populations in the interior part of the QTP. During the LGM, the survival of R. prainii in the interior part of the QTP corresponded with the in situ survival hypothesis, while R. stapfii most probably survived the LGM in a manner corresponding with the tabula rasa hypothesis. The evolutionary history of different lineages of this section was shaped by topography, climate change, and lineage-specific habitat preferences.

Mountains as Evolutionary Arenas: Patterns, Emerging Approaches, Paradigm Shifts, and Their Implications for Plant Phylogeographic Research in the Tibeto-Himalayan Region

Recently, the "mountain-geobiodiversity hypothesis" (MGH) was proposed as a key concept for explaining the high levels of biodiversity found in mountain systems of the Tibeto-Himalayan region (THR), which comprises the Qinghai-Tibetan Plateau, the Himalayas, and the biodiversity hotspot known as the "Mountains of Southwest China" (Hengduan Mountains region). In addition to the MGH, which covers the entire life span of a mountain system, a complementary concept, the so-called "flickering connectivity system" (FCS), was recently proposed for the period of the Quaternary. The FCS focuses on connectivity dynamics in alpine ecosystems caused by the drastic climatic changes during the past ca. 2.6 million years, emphasizing that range fragmentation and allopatric speciation are not the sole factors for accelerated evolution of species richness and endemism in mountains. I here provide a review of the current state of knowledge concerning geological uplift, Quaternary glaciation, and the main phylogeographic patterns ("contraction/recolonization," "platform refugia/local expansion," and "microrefugia") of seed plant species in the THR. In addition, I make specific suggestions as to which factors future avenues of phylogeographic research should take into account based on the fundamentals presented by the MGH and FCS, and associated complementary paradigm shifts.

Nunataks or massif de refuge? A phylogeographic study of Rhodiola crenulata (Crassulaceae) on the world's highest sky islands

Background

Quaternary climatic oscillations had tremendous effects on the current distribution of species. Here, we aim to elucidate the glacial history of Rhodiola crenulata, a perennial herb almost exclusively restricted to rock crevices on mountain peaks, and to test whether the nunatak or massif de refuge hypotheses could explain its distribution pattern.

Results

Six haplotypes and six ribotypes were detected in the cpDNA data set and the ITS data set, respectively. The divergence of R. crenulata and its closest relatives was dated have occurred ca. 0.65 Mya, during the Naynayxungla glaciation on the QTP. Mismatch distribution analysis suggested that the species experienced a range expansion around 0.31 Mya. Populations with high genetic and haplotype diversity were found on the QTP platform as well in the Hengduan Mountains. The ecological niche modeling results showed that there were suitable habitats on both the QTP platform and in the Hengduan Mountains during the LGM.

Conclusion

Our results support a scenario that both nunataks and the massif de refuge hypotheses could explain the distribution of R. crenulata. We also confirmed that Quaternary climatic oscillations could promote plant speciation in some circumstances. This study adds to a growing body of evidence suggesting that the QTP plant lineages exhibited diverse reactions to the Quaternary climatic oscillations.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1270-6) contains supplementary material, which is available to authorized users.

Rapid Intraspecific Diversification of the Alpine Species Saxifraga sinomontana (Saxifragaceae) in the Qinghai-Tibetan Plateau and Himalayas

An increasing number of phylogeographic studies have been conducted for plant species in the Qinghai-Tibetan Plateau (QTP) and its flanking mountains. However, these studies have mainly focused on the determination of glacial refugia and routes of inter-/post-glacial expansions. Rapid intraspecific diversification of plants in this region have not been thoroughly discussed. Herein, we investigate the effects of the Quaternary climate changes on population genetic structure and diversifications of a herbaceous alpine species, Saxifraga sinomontana, which may have an evolutionary time scale <5 million years in the QTP and Himalayan regions. Using a total of 350 individuals from 29 populations, we studied the evolutionary history of S. sinomontana by analyzing cpDNA trnL-trnF, rpl16 and nrDNA ITS sequences. A total of 89 haplotypes and 158 genotypes were detected for cpDNA and ITS sequences, respectively. Only a few haplotypes/genotypes were widespread, while an extremely large number of haplotypes/genotypes were restricted to single populations, which were scattered throughout the current geographical range of S. sinomontana. This suggests the existence of microrefugia of this species during the Quaternary glaciations. In addition, the relationships of the haplotypes/genotypes were almost completely not resolved by phylogenetic reconstruction. Combining characteristics in terms of high haplotype richness, large proportion of private haplotypes, and shallow haplotype divergence, we speculate that recent intraspecific diversification has occurred in S. sinomontana. Molecular clock analysis estimated that the onset diversification within S. sinomontana to be 1.09 Ma (95% HPD = 0.80–1.45), coinciding with the extensive Quaternary glaciations on the QTP which started ca. 1.17 Ma. The Quaternary climatic oscillations may have triggered rapid intraspecific diversification in this QTP-Himalayan species. However, large niche breadth, as well as introgression/hybridization between the studied species and its closely related sympatric saxifrages, may also played a role to some extent on the current genetic structure of S. sinomontana, which need to be further studied.

The effect of Hengduan Mountains Region (HMR) uplift to environmental changes in the HMR and its eastern adjacent area: Tracing the evolutionary history of Allium section Sikkimensia (Amaryllidaceae)

Exploring the effects of orographic events and climatic shifts on geographic distribution of organism in the Hengduan Mountains Region (HMR) and its eastern adjacent area is crucial to the understanding of the environmental changes to organismal evolution. To gain further insight into these processes, we reconstruct evolutionary history of ten species in Allium section Sikkimensia, distributed across regions abovementioned. Using chloroplast and nuclear sequence variation of 79 populations of these ten Allium species with known morphological preferences, we elucidate the phylogenetic relationship, divergence time, ancestral area and genetic structures. Climatic variables analysis, Isolation by distance (IBD) and environment (IBE) and Species distribution modeling (SDM) were analyzed along different genetic clades. These analyses indicated that the initial split of Sikkimensia was triggered by climate changes following Qinghai-Tibet Plateau sensu lato (QTPsl) uplift during the late Miocene. Subsequently, divergences within lineage (lineage A)/among lineages (lineage C and D) in Sikkimensia may be induced by the intense uplift of the HMR around 3-4 Ma and abrupt intensifying of the Asian monsoon regimes. Furthermore, Sikkimensia populations exhibited lopsided demographic history in the Last Glacial Maximum (LGM), as was indicated by the expansion of their range in the QDM and contraction in the HMR. Our findings appear to suggest that the HMR uplift could have strengthened the orographic difference between the HMR and its eastern adjacent area and led to a colder climate in the HMR, while geological topography also played an important role for taxa to respond the climate change that had taken place in the HMR and its eastern adjacent area during the Pleistocene.

Westwards and northwards dispersal of Triosteum himalayanum (Caprifoliaceae) from the Hengduan Mountains region based on chloroplast DNA phylogeography

The varying topography and environment that resulted from paleoorogeny and climate fluctuations of the Himalaya–Hengduan Mountains (HHM) areas had a considerable impact on the evolution of biota during the Quaternary. To understand the phylogeographic pattern and historical dynamics of Triosteum himalayanum (Caprifoliaceae), we sequenced three chloroplast DNA fragments (rbcL-accD, rps15-ycf1, and trnH-psbA) from 238 individuals representing 20 populations. Nineteen haplotypes (H1–H19) were identified based on 23 single-site mutations and eight indels. Most haplotypes were restricted to a single population or neighboring populations. Analysis of molecular variance revealed that variations among populations were much higher than that within populations for the overall gene pool, as well as for the East Himalayan group (EH group) and the North Hengduan group (NHM group), but not for the Hengduan Mountains group (HM group). Ecoregions representing relatively high genetic diversity or high frequencies of private haplotypes were discovered, suggesting that this alpine herbaceous plant underwent enhanced allopatric divergence in isolated and fragmented locations during the Quaternary glaciations. The current phylogeographic structure of T. himalayanum might be due to heterogeneous habitats and Quaternary climatic oscillations. Based on the phylogeographic structure of T. himalayanum populations, the phylogenetic relationship of identified haplotypes and palaeodistributional reconstruction, we postulated both westwards and northwards expansion from the HM group for this species. The westwards dispersal corridor could be long, narrow mountain areas and/or the Yarlung Zangbo Valley, while the northwards movement path could be south–north oriented mountains and low-elevation valleys.

Climate Is Not All: Evidence From Phylogeography of Rhodiola fastigiata (Crassulaceae) and Comparison to Its Closest Relatives

How geological events and climate oscillations in the Pleistocene glaciation shaped the geographic distribution of genetic variation of species on the Qinghai-Tibetan Plateau (QTP) and its adjacent areas has been extensively studied. However, little studies have investigated whether closely related species in the same genus with similar physiological and life history traits responded similarly to the glacial climatic oscillations. If this is not the case, we would expect that the population histories of studied species were not driven by extrinsic environmental changes alone. Here we conducted a phylogeographic study of a succulent alpine plant Rhodiola fastigiata, using sequences from chloroplast genome and nrITS region, as well as ecological niche modeling. The results of R. fastigiata were compared to other congeneric species that have been studied, especially to R. alsia and R. crenulata. We found that for both markers, two geographic groups could be revealed, corresponding to the QTP plateau and the Hengduan Mountains, respectively, indicating isolated refugia in those two areas. The two groups diverged 1.23 Mya during the Pleistocene. We detected no significant population expansion by mismatch distribution analysis and Bayesian Skyline Plot. We found that even these similar species with similar physiological and life history traits have had different demographic histories in the Quaternary glacial periods. Our comparative phylogeographic study sheds new lights into phylogeographic research that extrinsic environmental changes are not the only factor that can drive population demography, and other factors, such as coevolved interactions between plants and their specialized pathogens, that probably played a role need to be examined with more case studies.

The demographic response of a deciduous shrub (the Indigofera bungeana complex, Fabaceae) to the Pleistocene climate changes in East Asia

East Asia harbors the highest level of floristic diversity among the world’s temperate regions. Despite the increase in phylogeographic studies of temperate plants in East Asia, far less attention has been paid to widely distributed deciduous shrubs that widespread across several floral regions. We sequenced two chloroplast DNA (cpDNA) fragments (ndhJ-trnF and trnD-trnT) and one nuclear DNA (Pgk1) of 472 individuals from 51 populations of such a group, the Indigofera bungeana complex. We used population genetic data as well as ecological niche modelling to examine the evolutionary history and glacial refugia during the Last Glacial Maximum (LGM) of this group. We recovered 133 cpDNA and 68 nuclear haplotypes. The star-phylogeny of the recovered cpDNA and nuclear haplotypes and demographic analyses suggested distinct range expansion of I. bungeana complex have occurred during the early and middle Pleistocene. The climate change of the LGM might have affected little on the distribution of this complex based on the niche modelling. However, these climate changes and geographic isolation probably resulted in fixtures of the private haplotypes and genetic differentiations between regions. Our results suggested that this arid-tolerant species complex may have different responses to the Quaternary climate changes with those climate-sensitive species.

Phylogeographic Insights into a Peripheral Refugium: The Importance of Cumulative Effect of Glaciation on the Genetic Structure of Two Endemic Plants

Quaternary glaciations and mostly last glacial maximum have shaped the contemporary distribution of many species in the Alps. However, in the Maritime and Ligurian Alps a more complex picture is suggested by the presence of many Tertiary paleoendemisms and by the divergence time between lineages in one endemic species predating the Late Pleistocene glaciation. The low number of endemic species studied limits the understanding of the processes that took place within this region. We used species distribution models and phylogeographical methods to infer glacial refugia and to reconstruct the phylogeographical pattern of Silene cordifolia All. and Viola argenteria Moraldo & Forneris. The predicted suitable area for last glacial maximum roughly fitted current known distribution. Our results suggest that separation of the major clades predates the last glacial maximum and the following repeated glacial and interglacial periods probably drove differentiations. The complex phylogeographical pattern observed in the study species suggests that both populations and genotypes extinction was minimal during the last glacial maximum, probably due to the low impact of glaciations and to topographic complexity in this area. This study underlines the importance of cumulative effect of previous glacial cycles in shaping the genetic structure of plant species in Maritime and Ligurian Alps, as expected for a Mediterranean mountain region more than for an Alpine region.

The Quaternary evolutionary history, potential distribution dynamics, and conservation implications for a Qinghai-Tibet Plateau endemic herbaceous perennial, Anisodus tanguticus (Solanaceae)

Various hypotheses have been proposed about the Quaternary evolutionary history of plant species on the Qinghai-Tibet Plateau (QTP), yet only a handful of studies have considered both population genetics and ecological niche context. In this study, we proposed and compared climate refugia hypotheses based on the phylogeographic pattern of Anisodus tanguticus (three plastid DNA fragments and nuclear internal transcribed spacer regions from 32 populations) and present and past species distribution models (SDMs). We detected six plastid haplotypes in two well-differentiated lineages. Although all haplotypes could be found in its western (sampling) area, only haplotypes from one lineage occurred in its eastern area. Meanwhile, most genetic variations existed between populations (F ST = 0.822). The SDMs during the last glacial maximum and last interglacial periods showed range fragmentation in the western area and significant range contraction in the eastern area, respectively, in comparison with current potential distribution. This species may have undergone intraspecific divergence during the early Quaternary, which may have been caused by survival in different refugia during the earliest known glacial in the QTP, rather than geological isolation due to orogenesis events. Subsequently, climate oscillations during the Quaternary resulted in a dynamic distribution range for this species as well as the distribution pattern of its plastid haplotypes and nuclear genotypes. The interglacial periods may have had a greater effect on A. tanguticus than the glacial periods. Most importantly, neither genetic data nor SDM alone can fully reveal the climate refugia history of this species. We also discuss the conservation implications for this important Tibetan folk medicine plant in light of these findings and SDMs under future climate models. Together, our results underline the necessity to combine phylogeographic and SDM approaches in future investigations of the Quaternary evolutionary history of species in topographically complex areas, such as the QTP.

Refugial isolation and range expansions drive the genetic structure of Oxyria sinensis (Polygonaceae) in the Himalaya-Hengduan Mountains

The formation of the Mekong-Salween Divide and climatic oscillations in Pleistocene were the main drivers for the contemporary diversity and genetic structure of plants in the Himalaya-Hengduan Mountains (HHM). To identify the relative roles of the two historical events in shaping population history of plants in HHM, we investigated the phylogeographic pattern of Oxyria sinensis, a perennial plant endemic to the HHM. Sixteen chloroplast haplotypes were identified and were clustered into three phylogenetic clades. The age of the major clades was estimated to be in the Pleistocene, falling into several Pleistocene glacial stages and postdating the formation of the Mekong-Salween Divide. Range expansions occurred at least twice in the early and middle Pleistocene, but the spatial genetic distribution rarely changed since the Last Glacial Maximum. Our results suggest that temporary mountain glaciers may act as barriers in promoting the lineage divergence in O. sinensis and that subsequential range expansions and secondary contacts might reshape the genetic distribution in geography and blur the boundary of population differentiation created in the earlier glacial stages. This study demonstrates that Pleistocene climatic change and mountain glaciers, rather than the Mekong-Salween Divide, play the primary role in shaping the spatial genetic structure of O. sinensis.

Phylogeographical structure inferred from cpDNA sequence variation of Zygophyllum xanthoxylon across north-west China

Zygophyllum xanthoxylon, a desert species, displaying a broad east-west continuous distribution pattern in arid Northwestern China, can be considered as a model species to investigate the biogeographical history of this region. We sequenced two chloroplast DNA spacers (psbK-psbI and rpl32-trnL) in 226 individuals from 31 populations to explore the phylogeographical structure. Median-joining network was constructed and analysis of AMOVA, SMOVA, neutrality tests and distribution analysis were used to examine genetic structure and potential range expansion. Using species distribution modeling, the geographical distribution of Z. xanthoxylon was modeled during the present and at the Last Glacial Maximum (LGM). Among 26 haplotypes, one was widely distributed, but most was restricted to either the eastern or western region. The populations with the highest levels of haplotype diversity were found in the Tianshan Mountains and its surroundings in the west, and the Helan Mountains and Alxa Plateau in the east. AMOVA and SAMOVA showed that over all populations, the species lacks phylogeographical structure, which is speculated to be the result of its specific biology. Neutrality tests and mismatch distribution analysis support past range expansions of the species. Comparing the current distribution to those cold and dry conditions in LGM, Z. xanthoxylon had a shrunken and more fragmented range during LGM. Based on the evidences from phylogeographical patterns, distribution of genetic variability, and paleodistribution modeling, Z. xanthoxylon is speculated most likely to have originated from the east and migrated westward via the Hexi Corridor.

Phylogeography of Rhodiola kirilowii (Crassulaceae): a story of Miocene divergence and quaternary expansion

The evolution and current distribution of the Sino-Tibetan flora have been greatly affected by historical geological events, such as the uplift of the Qinghai-Tibetan Plateau (QTP), and Quaternary climatic oscillations. Rhodiola kirilowii, a perennial herb with its distribution ranging from the southeastern QTP and the Hengduan Mountains (HM) to adjacent northern China and central Asia, provides an excellent model to examine and disentangle the effect of both geological orogeny and climatic oscillation on the evolutionary history of species with such distribution patterns. We here conducted a phylogeographic study using sequences of two chloroplast fragments (trnL-F and trnS-G) and internal transcribed spacers in 29 populations of R. kirilowii. A total of 25 plastid haplotypes and 12 ITS ribotypes were found. Molecular clock estimation revealed deep divergence between the central Asian populations and other populations from the HM and northern China; this split occurred ca. 2.84 million year ago. The majority of populations from the mountains of northern China were dominated by a single haplotype or ribotype, while populations of the HM harbored both high genetic diversity and high haplotype diversity. This distribution pattern indicates that HM was either a diversification center or a refugium for R. kirilowii during the Quaternary climatic oscillations. The present distribution of this species on mountains in northern China may have resulted from a rapid glacial population expansion from the HM. This expansion was confirmed by the mismatch distribution analysis and negative Tajima's D and Fu's FS values, and was dated to ca. 168 thousand years ago. High genetic diversity and population differentiation in both plastid and ITS sequences were revealed; these imply restricted gene flow between populations. A distinct isolation-by-distance pattern was suggested by the Mantel test. Our results show that in old lineages, populations may harbour divergent genetic forms that are sufficient to maintain or even increase overall genetic diversity despite fragmentation and low within-population variation.

Phytochemical profiles, antioxidant and antimicrobial activities of three Potentilla species

Background

Extracts from Potentilla species have been applied in traditional medicine and exhibit antioxidant, hypoglycemic, anti-inflammatory, antitumor and anti-ulcerogenic properties, but little has been known about the diversity of phytochemistry and pharmacology on this genus. This study investigated and compared the phytochemical profiles, antioxidant and antimicrobial activities of leaf extracts from three Potentilla species (Potentilla fruticosa, Potentilla glabra and Potentilla parvifolia) in order to discover new resources for lead structures and pharmaceutical products.

Methods

Chemical composition and content of six phenolic compounds were evaluated and determined by RP-HPLC; Total phenolic and total flavonoid content were determined using Folin-Ciocalteau colourimetric method and sodium borohydride/chloranil-based method (SBC); Antioxidant activities were determined using DPPH, ABTS and FRAP assays; Antimicrobial properties were investigated by agar dilution and mycelial growth rate method.

Results

The results showed hyperoside was the predominant phenolic compound in three Potentilla species by RP-HPLC assay, with the content of 8.86 (P. fruticosa), 2.56 (P. glabra) and 2.68 mg/g (P. parvifolia), respectively. The highest content of total identified phenolic compounds (hyperoside, (+)-catechin, caffeic acid, ferulic acid, rutin and ellagic acid) was observed in P. parvifolia (14.17 mg/g), follow by P. fruticosa (10.01 mg/g) and P. glabra (7.01 mg/g). P. fruticosa possessed the highest content of total phenolic (84.93 ± 0.50 mmol gallic acid equivalent/100 g) and total flavonoid (84.14 ± 0.03 mmol quercetin equivalent/100 g), which were in good correlation with its significant DPPH_IC50 (16.87 μg/mL), ABTS (2763.48 μmol Trolox equivalent/g) and FRAP (1398.70 μmol Trolox equivalent/g) capacities. Furthermore, the effective methodology to distinguish the different species of Potentilla was also established by chromatographic fingerprint analysis for the first time. The results of antimicrobial activities showed P. fruticosa exhibited the strongest inhibition aganist Gram-positive bacteria, Pseudomonas aeruginosa and Candida albicans with MIC values of 0.78–6.25 mg/mL. P. parvifolia possessed antibacterial and antifungal activities against all the microorganisms tested, with EC₅₀ and MIC values of 20.52–47.02 mg/mL and 0.78–50 mg/mL, respectively.

Conclusions

These results indicated that leaf extracts from three Potentilla species could become useful supplement for pharmaceutical products as a new antioxidant and antimicrobial agents.

Chloroplast DNA phylogeography reveals repeated range expansion in a widespread aquatic herb Hippuris vulgaris in the Qinghai-Tibetan Plateau and adjacent areas

Background

The Qinghai-Tibetan Plateau (QTP) is one of the most extensive habitats for alpine plants in the world. Climatic oscillations during the Quaternary ice age had a dramatic effect on species ranges on the QTP and the adjacent areas. However, how the distribution ranges of aquatic plant species shifted on the QTP in response to Quaternary climatic changes remains almost unknown.

Methodology and Principal Findings

We studied the phylogeography and demographic history of the widespread aquatic herb Hippuris vulgaris from the QTP and adjacent areas. Our sampling included 385 individuals from 47 natural populations of H. vulgaris. Using sequences from four chloroplast DNA (cpDNA) non-coding regions, we distinguished eight different cpDNA haplotypes. From the cpDNA variation in H. vulgaris, we found a very high level of population differentiation (G_ST = 0.819) but the phylogeographical structure remained obscure (N_ST = 0.853>G_ST = 0.819, P>0.05). Phylogenetic analyses revealed two main cpDNA haplotype lineages. The split between these two haplotype groups can be dated back to the mid-to-late Pleistocene (ca. 0.480 Myr). Mismatch distribution analyses showed that each of these had experienced a recent range expansion. These two expansions (ca. 0.12 and 0.17 Myr) might have begun from the different refugees before the Last Glacial Maximum (LGM).

Conclusions/Significance

This study initiates a research on the phylogeography of aquatic herbs in the QTP and for the first time sheds light on the response of an alpine aquatic seed plant species in the QTP to Quaternary climate oscillations.

Population genetic structure and geographic differentiation in the hot spring snake Thermophis baileyi (Serpentes, Colubridae): indications for glacial refuges in southern-central Tibet

The endangered hot spring snake (Thermophis baileyi) is endemic to the Tibetan Plateau and restricted to a few sites at high altitudes (>3500 m a.s.l.). Its strong preference for habitats with hot springs might be an early adaptation to the cooling climate during the plateau uplift. Some of these thermal sites may have been free of ice during the last glacial maximum (LGM) serving as refuges for the snake. To test this and to examine the population genetic structure data of 12 microsatellites and mitochondrial DNA were obtained from 153 individuals from 12 sites across the plateau. Four of the microsatellite loci were newly developed for this study. The results suggest that T. baileyi has at least two genetically diverse clades in Tibet, which developed during the Pleistocene and expanded after the LGM. Accordingly, the existence of separate glacial refuges on the central plateau can be assumed. Analyses of the genetic variation indicated a high level of geographic differentiation and population structure on a regional as well as on a range-wide scale. The study shows that, apart from the phylogeographic signatures, the diversification of current Thermophis populations is caused by (a) a limited dispersal due to mountain ranges, (b) a strong preference for hot springs and (c) the insular distribution of suitable habitats on the plateau. The present data should be taken into account in necessary monitoring and management plans for T. baileyi.

Phylogeography of the Sino-Himalayan fern Lepisorus clathratus on "the roof of the world"

BACKGROUND

The Qinghai-Tibetan Plateau (QTP) and its southern and southeastern mountain ranges, Himalaya-Hengduan Mountains (HHM), are one of the most extensive habitats for alpine plants in the world. How ferns occurring in QTP and HHM changed their distribution ranges in response to Quaternary climatic oscillations remains almost unknown.

METHODOLOGY AND RESULTS

We employed sequences of two chloroplast DNA regions, rps4-trnS and trnL-trnF, to reconstruct phylogeography of the Sino-Himalayan fern Lepisorus clathratus, occurring mainly in the QTP and HHM. Individuals of this species have either dehiscent or indehiscent sporangia with the latter evolved from the plesiomorphic dehiscent forms. Based on a range-wide sampling, we detected 27 cpDNA haplotypes that were divided into five groups by network analyses. Populations in the Hengduan Mountains possess the highest genetic diversity, while a single haplogroup is detected across the north-central region. A distinct phylogeographical subdivision was detected between the Hengduan Mountains and north-central region by AMOVA analysis. The haplogroup distribution pattern, coalescence and AMOVA analysis suggest that a long term survival area (refugia) of the species was located in the Hengduan Mountains during glaciations, with probable range expansions into north-central regions during interglacial periods. Populations with indehiscent sporangium can carry private haplotypes and are inclined to maintain genetic homogeneity. One group with indehiscent sporangia most likely survived in situ on the QTP during glaciations.

CONCLUSIONS/SIGNIFICANCE

This study for the first time sheds light on the response of alpine ferns in the QTP and HHM to the Quaternary climatic oscillations.

Phylogeography and allopatric divergence of cypress species (Cupressus L.) in the Qinghai-Tibetan Plateau and adjacent regions

BACKGROUND

Although allopatric speciation is viewed as the most common way in which species originate, allopatric divergence among a group of closely related species has rarely been examined at the population level through phylogeographic analysis. Here we report such a case study on eight putative cypress (Cupressus) species, which each have a mainly allopatric distribution in the Qinghai-Tibetan Plateau (QTP) and adjacent regions. The analysis involved sequencing three plastid DNA fragments (trnD-trnT, trnS-trnG and trnL-trnF) in 371 individuals sampled from populations at 66 localities.

RESULTS

Both phylogenetic and network analyses showed that most DNA haplotypes recovered or haplotype-clustered lineages resolved were largely species-specific. Across all species, significant phylogeographic structure (N(ST) > G(ST), P < 0.05) implied a high correlation between haplotypes/lineages and geographic distribution. Two species, C. duclouxiana and C. chengiana, which are distributed in the eastern QTP region, contained more haplotypes and higher diversity than five species with restricted distributions in the western highlands of the QTP. The remaining species, C. funebris, is widely cultivated and contained very little cpDNA diversity.

CONCLUSIONS

It is concluded that the formation of high mountain barriers separating deep valleys in the QTP and adjacent regions caused by various uplifts of the plateau since the early Miocene most likely promoted allopatric divergence in Cupressus by restricting gene flow and fixing local, species-specific haplotypes in geographically isolated populations. The low levels of intraspecific diversity present in most species might stem from population bottlenecks brought about by recurrent periods of unfavorable climate and more recently by the negative impacts of human activities on species' distributions. Our findings shed new light on the importance of geographical isolation caused by the uplift of the QTP on the development of high plant species diversity in the QTP biodiversity hotspot.