Effects of Mountain Uplift and Climatic Oscillations on Phylogeography and Species Divergence in Four Endangered Notopterygium Herbs

Abstract

Mountain uplift and climatic fluctuations are important driving forces that have affected the geographic distribution and population dynamics history of organisms. However, it is unclear how geological and climatic events might have affected the phylogeographic history and species divergence in high-alpine herbal plants. In this study, we analyzed the population demographic history and species differentiation of four endangered Notopterygium herbs on the high-altitude Qinghai–Tibetan Plateau (QTP) and adjacent areas. We combined phylogeographic analysis with species distribution modeling to detect the genetic variations in four Notopterygium species (N. incisum, N. franchetii, N. oviforme, and N. forrestii). In total, 559 individuals from 74 populations of the four species were analyzed based on three maternally inherited chloroplast fragments (matK, rbcL, and trnS-trnG) and one nuclear DNA region (internal transcribed spacer, ITS). Fifty-five chloroplast DNA (cpDNA) and 48 ITS haplotypes were identified in the four species. All of the cpDNA and ITS haplotypes were species-specific, except N. franchetii and N. oviforme shared one cpDNA haplotype, H32. Phylogenetic analysis suggested that all four species formed a monophyletic clade with high bootstrap support, where N. franchetii and N. oviforme were sisters. In addition, each Notopterygium species generated an individual clade that corresponded to their respective species in the ITS tree. Population dynamics analyses and species distribution modeling showed that the two widely distributed herbs N. incisum and N. franchetii exhibited obvious demographic expansions during the Pleistocene ice ages. Molecular dating suggested that the divergence of the four Notopterygium species occurred approximately between 3.6 and 1.2 Mya, and it was significantly associated with recent extensive uplifts of the QTP. Our results support the hypothesis that mountain uplift and Quaternary climatic oscillations profoundly shaped the population genetic divergence and demographic dynamics of Notopterygium species. The findings of this and previous studies provide important insights into the effects of QTP uplifts and climatic changes on phylogeography and species differentiation in high altitude mountainous areas. Our results may also facilitate the conservation of endangered herbaceous medicinal plants in the genus Notopterygium.

Introduction

Geological events and climatic fluctuations are considered to have profoundly shaped the distribution and population dynamics history of species in mountain areas (Hewitt, 2004; Hickerson et al., 2010). Thus, during glacial periods, most species experienced adverse weather conditions in high altitude mountains, where they contracted into refugia in low latitudes and then their ranges expanded again after the ice ages, thereby leading to species divergence or secondary contact evolution (Hewitt, 2004; Nybom, 2004; Ohsawa and Ide, 2008). In addition, the distribution patterns and population genetic structures of some species were reshaped due to geographic barriers and climatic oscillations. Moreover, the population size, mating system, and bio-characteristics of species had important effects on the divergence and evolutionary history of populations of species (Stewart et al., 2010; Zhang et al., 2015). For example, some studies have suggested that the locations of ice age refugia for plants were determined mainly by the adaptability of species to the external environment (Stewart et al., 2010; Liao et al., 2015; Wang et al., 2015).

During the Quaternary ice periods, many species experienced extinction events due to repeated bottlenecks and genetic drift, which led to further divergent evolution and genetic isolation within species (Soltis et al., 2006; Hickerson et al., 2010; Stewart et al., 2010; Keppel et al., 2012). Repeated environmental changes may also have promoted the fragmentation of habitats, as well as causing exotic distributions of different species or intraspecific genetic changes (Hickerson et al., 2010; Jia et al., 2012). Studies of alpine trees have shown that populations were crossed during mountain uplift processes, whereas the exchange of genes among populations was restricted due to climatic and geographic barriers (McLachlan et al., 2005; Soltis et al., 2006; Birks and Willis, 2008; Avise, 2009; Gao et al., 2012). Some alpine species experienced deep lineage divergence due to climatic changes and environmental isolations (Qiu et al., 2011; Liu et al., 2012; Xu et al., 2015).

In the high latitudes of Europe and North America, studies suggested that plant species could have survived in high elevation areas (“invisible refugia”) during the ice age periods (Anderson et al., 2006; Provan and Bennett, 2008; Parducci et al., 2012; Stewart and Stringer, 2012; Ortego et al., 2012, 2015; Guichoux et al., 2013; Allen et al., 2015; Cavender-Bares et al., 2015). The presence of Juniperus species in the Qinghai-Tibetan Plateau (QTP) region also supports the existence of invisible refugia (Opgenoorth et al., 2010). The QTP is the largest and highest plateau in the world, with a mean altitude of more than 4000 m. Various endangered species and high levels of global diversity are present on this plateau (Mittermeier et al., 2005). Studies suggest that extensive uplifts of the QTP occurred in the Miocene–Pliocene era between 3.6 and 1.7 Mya (Li and Fang, 1999; Zhang et al., 2000; Zhou et al., 2006). The lifting of mountains triggered species divergence and changed the genetic structure to affect the evolution of high-alpine plants (Liu et al., 2012; Wen et al., 2014; Ickert-Bond and Renner, 2016). In particular, the geological effects of the QTP on the genetic structure, geographic distribution, and species differentiation of plants have been clearly defined in this area (Wang et al., 2010; Xu et al., 2010; Li et al., 2013; Liu et al., 2013; Sun et al., 2014; Favre et al., 2015; Hughes and Atchison, 2015). However, most of these previous studies focused on the response patterns of tree or shrub species to mountain uplifts and climatic oscillations on the QTP (Liu et al., 2006; Wang et al., 2009; Mao et al., 2010; Xu et al., 2010; Tian et al., 2011; Qiu et al., 2011; Wen et al., 2014; Favre et al., 2015; Hughes and Atchison, 2015), whereas little is known about the effects of mountain uplifts and climate events on cold-tolerant herbal species in the high altitude QTP and adjacent regions.

The genus Notopterygium H. de Boissieu (Apiaceae) comprises perennial and endangered herbaceous medicinal plants, which are mainly distributed in the QTP and its surrounding high-altitude areas. According to records in the Flora of China, this genus comprises six species: N. incisum C. C. Ting ex H. T. Chang, N. oviforme R. H. Shan, N. franchetii H. de Boissieu, N. forrestii H. Wolff, N. tenuifolium M. L. Sheh and F. T. Pu, and N. pinnatiinvolucellum F. T. Pu and Y. P. Wang. N. incisum and N. franchetii have wide distribution ranges at altitudes of 3200–5100 m and 1700–4500 m, respectively. N. oviforme occurs in the eastern part of the QTP at altitudes of 1700–3200 m. The other three species, i.e., N. forrestii (4000–4300 m), N. tenuifolium (4300 m), and N. pinnatiinvolucellum (3400 m), have very limited distributions among the high-alpine shrubs and meadows in the west region of China. These herb species provide an excellent model for detecting the effects of the QTP uplifts and Quaternary climatic oscillations on the genetic structure and species divergence of plants. However, in recent years, due to high market demand, the wild resources of these Notopterygium species have decreased rapidly because of human over-exploitation (Zhou et al., 2010). The Notopterygium species are now listed as endangered herb species in the IUCN Red List, and their management and conservation are urgently required (Wu et al., 2005). Information regarding geographic distributions and genetic diversity is vital for formulating effective conservation strategies for wild plant resources. However, most of the previous studies of the Notopterygium species have focused mainly on their phylogenetic evolutionary relationships (Pu et al., 2000; Yang et al., 2017), morphological and physiological characteristics (She and Pu, 1996; Wang et al., 1996; Jiang et al., 2005), and comparative transcriptome analysis (Jia et al., 2017), whereas little is known about their genetic divergence and population demographic history.

In the current study, we sampled four species, i.e., N. incisum, N. oviforme, N. franchetii, and N. forrestii, across their entire geographic distributions in the high-altitude QTP and adjacent areas. We detected the genetic variations in three chloroplast DNA (cpDNA) markers and a nuclear DNA fragment in order to characterize the population histories and species divergence of these endangered herb plants. Our aims were: (1) to determine the genetic structure and population evolutionary history of four Notopterygium species; (2) to identify the phylogenetic relationships among these species and their phylogeographic history; (3) to explore the effects of QTP uplifts and climatic changes in the Quaternary on the divergence and phylogeography of these species; and (4) to propose reasonable conservation and management strategies for the endangered Notopterygium species.

Materials and Methods

Sample Collection

In this study, in order to obtain information about genetic variation over a wide area, 559 individuals from 74 populations were collected for the four Notopterygium species in Sichuan, Shaanxi, Gansu, Qinghai, and Shanxi provinces in the high-altitude QTP and adjacent areas. These samples covered the complete geographic distribution ranges of the four species in the QTP and surrounding areas. From 2 to 17 individuals were sampled from each population, where all of the samples collected were separated from each other by at least 100 m. Detailed information about the latitude, longitude, and altitude for all of the populations is provided in Figure 1 and Table 1. All of the materials and documents have been deposited in the College of Life Sciences, Northwest University. In addition, two species from the genus Pleurospermum, i.e., P. prattii and P. franchetianum, as well as Heracleum moellendorffii were used as outgroups.

FIGURE 1

DNA Extraction and Sequencing

Total DNA was extracted using the modified CTAB method (Doyle and Doyle, 1987) or with a plant DNA extraction kit (Tiangen, Beijing, China). We used 1% agarose gels to check the quality of the DNA extracted from the Notopterygium species. To screen for suitable primers, we first randomly selected 50 individuals (one individual from each population) to amplify the universal cpDNA primers and nDNA primers recommended by the Consortium for the Barcode of Life (CBOL) (CBOL Plant Working Group, 2009). Finally, three highly variable cpDNA primers, i.e., trnS-trnG, matK, and rbcL, and one nDNA internal transcribed spacer (ITS) primer were selected to determine the genetic variations in the genus Notopterygium after initial tests with six loci (Supplementary Table S1). Thus, two monomorphic cpDNA loci (trnL-trnF and rpl36-infA; Bai et al., 2010; Soumaya et al., 2014) were excluded from all of the subsequent analyses.

PCR amplification was performed in a volume of 25 μL containing 2 μL DNA template (10–50 ng/μL), 12.5 μL PCR MIX (Xi’an Runde, China), 0.75 μL of each primer (20 ng/μL), and 9 μL double-distilled H₂O. The PCR reaction conditions were as described in Supplementary Table S2. All of the high quality PCR products were sequenced using the amplified forward and reverse primers with an ABI 3730 XL genetic analyzer (Applied Biosystems, Foster City, CA, United States). All of the sequences were deposited in GenBank under accession numbers MF787385–MF787579.

Proofreading and Alignment of DNA, and Data Analysis

BioEdit v 7.0.9.0 (Hall, 1999) software was used for manual proofreading and checking the variable sites. MEGA v 7.0 (Sudhir et al., 2008) was used to remove low quality sequences and only high quality sequences were analyzed. For the ITS sequences, we visualized the possible color spectrum of the overlapping peaks at any one variable site. If a strong signal peak was more than half of a weak signal peak, then we used the strong peak for phrasing. If both the peaks overlapped, we used the following phrases for each variable site instead of both peaks as phrases: R: A+G, Y: C+T, M: A+C, K: G+T, S: G+C, W: A+T. DnaSP v 5.0 software was used for dividing the heterozygous loci into double sequence series (Librado and Rozas, 2009).

Genetic Variation and Genetic Structure Analysis

The genetic diversity of the cpDNA and ITS sequences were analyzed in all four Notopterygium species using PERMUT v 1.0 software, where we calculated the genetic diversity within the population of each species (h_S), total genetic diversity (h_T), and population genetic differentiation coefficients G_ST and N_ST (Grivet and Petit, 2002).

In addition, ARLEQUIN v 3.5 (Excoffier and Lischer, 2010) software was used to perform analysis of molecular variance (AMOVA) for the cpDNA and ITS sequences. AMOVA partitioned the genetic differentiation among the populations F_ST, within a population F_SC, and among species F_CT.

Phylogenetic Analysis

Phylogenetic analyses of the cpDNA and ITS sequences were performed with MEGA v 7.0. JModeltest v 3.06 (Posada and Crandall, 1998) was used to filter the best evolutionary model (GTR+G). One-thousand bootstrap replicates were performed for the maximum likelihood (ML) and maximum parsimony (MP) models to obtain the phylogenetic tree. MrBayes v 3.2.3 was also used to conduct phylogenetic analyses of the cpDNA and ITS sequences based on the Bayesian criterion (Huelsenbeck and Ronquist, 2001). We set the random tree rotation as 10,000,000 generation, where each 1000 generations were kept to construct a phylogenetic tree, with a burn-in of 2500.

NETWORK v 5.0.0 (Polzin and Daneshmand, 2003) was used to construct median-joining networks of the cpDNA and ITS sequences. ArcGIS v 10.2 (Bader, 2005) was employed to draw the haplotype distribution map. BEAST v 1.7.5 (Drummond and Rambaut, 2007) was used to estimate the divergence times of the cpDNA haplotypes where we used the cpDNA evolutionary rates (1.0–3.0 × 10^-9 s/s/y) recorded for other angiosperms to calibrate our datasets due to the lack of fossil evidence for Notopterygium plants (Wolfe et al., 1987). We employed the loose molecular clock method with an uncorrected log-normal distribution for the branch lengths. After a burn-in of 5,000,000 steps, all of the parameters were collected once every 1000 steps up to 50,000,000 Markov chain Monte Carlo (MCMC) algorithm steps. The convergence of the MCMC results was verified by using the Tracer v 1.5 program to check that the chain was balanced, where we then used the Tree Annotator v 1.7.5 program to obtain the best tree merging and Figtree v 1.3.1 (Rambaut, 2009) was employed to view the resulting tree.

Population Dynamics Analysis

DnaSP v 5.0 was used to analyze the genetic diversity parameters, including the haplotype diversity (H_d) (Nei and Tajima, 1981), nucleotide diversity (π) (Nei and Li, 1979), and number of haplotypes (H). We also used DnaSP v 5.0 to detect the mismatched distributions (Schneider and Excoffier, 1999) of cpDNA sequences in the four Notopterygium species. Population demographic expansions were tested using Arlequin v 3.5 (Excoffier and Lischer, 2010) and Tajima’s D (Tajima, 1989), Fu’s F_S (Fu, 1997), and Fu and Li’s F^∗ (Fu and Li, 1993) tests. We used the sum of the squared deviations between the observed and expected mismatches as well as Harpending’s raggedness index values (Rag) (Harpending, 1994) to determine the validity and significance level of the expansion model. According to the formula: τ =2ut (τ is the mismatch equilibrium expansion variable) (Rogers and Harpending, 1992), we calculated the expansion time t, where u is the mutation rate per generation calculated using the formula u = 2μkg, where μ is the mutation rate per nucleotide per year, k is the total length of a cpDNA sequence, and g is the generation time. According to our field investigations, the generation time for Notopterygium species was 3 years.

In order to further determine the signs of demographic growth in the four Notopterygium species, we used LAMARC v 2.1.8 (Kuhner, 2006) to calculate the population growth parameter g. The MCMC algorithm was run for 100,000 generations and sampled every 200,000 steps, where the first 25% of the sampled trees were discarded as the burn-in.

Species Distribution Modeling

We used MaxEnt v 3.3.3k (Phillips et al., 2006; Phillips and Dudík, 2008) to predict the current, last glacial maximum (LGM), last interglacial (LIG), and future distributions of two widespread Notopterygium species: N. incisum (148 distribution sites) and N. franchetii (80 distribution sites). The distribution sites of Notopterygium species were collected from previous studies as well as websites containing climate data and plant distributions. We also obtained some distribution sites based on field investigations. Bio-climatic environment data were downloaded from the WorldClim website¹ at a resolution of 2.5 arc-minutes. Six bioclimatic environmental variables (Supplementary Table S7) with significant effects on N. incisum and N. franchetii were used to detect changes in the distribution ranges of plants. We set the number of replicates to 10 and the maximum number of iterations to 500 for MaxEnt modeling. The accuracy of the model’s performance was assessed based on the area under the receiver operating characteristic curve (AUC) (Fawcett, 2006).

Results

cpDNA Variations and Haplotype Distributions

Three chloroplast fragments (matK, rbcL, and trnS-trnG) were used to analyze 559 individuals from 74 populations of the four Notopterygium species. The total length of the fragments was 1605 bp, and the lengths of the matK, rbcL, and trnS-trnG regions were 669, 668, and 268 bp, respectively, which included 21, seven, and eight nucleotide mutation sites (Supplementary Table S3). The cpDNA regions were uniparental inherited markers so we combined the three chloroplast fragments in the subsequent population genetics analysis.

In total, 55 cpDNA haplotypes were detected in the four Notopterygium species (Figures 1, 2). Most of the haplotypes were species-specific, except the H32 haplotype was shared by N. franchetii and N. oviforme. N. incisum contained 31 haplotypes, where haplotypes H1–H7, H12, H18, and H26 were shared among populations, and the remainder were unique to each population. Populations from the southeast part of the QTP (G, J, Q, U, Z, HB, and HC; see Table 1 for the site codes) had the highest haplotype diversity. N. franchetii had nine haplotypes with nine mutation sites. Populations from the west part of China (KG, KO, KZ, and YF) also had the highest haplotype diversity. N. oviforme had 12 haplotypes with 16 mutation sites, where the LE population had the highest haplotype diversity for this species. N. forrestii had three haplotypes with two mutation sites, where H55 was unique to the LCB population, and haplotypes H53 and H54 were shared by the other populations.

FIGURE 2

Notopterygium oviforme had the highest levels of genetic diversity and nucleotide diversity (H_d = 0.81, π = 0.0013), followed by N. incisum (H_d = 0.75, π = 0.00086) and N. forrestii (H_d = 0.39, π = 0.0002), whereas N. franchetii had the lowest level of diversity (H_d = 0.29, π = 0.00031) (Table 2).

ITS Sequence Variation

The total length of the sequenced ITS region was 593 bp and 48 haplotypes were identified with 66 nucleotide mutation sites (Figure 2 and Supplementary Table S4). All of the ITS haplotypes were species-specific in the four Notopterygium species. The total haplotype diversity (H_d) and π values for N. incisum, N. franchetii, and N. oviforme were 0.71 and 0.0025, 0.55 and 0.00364, and 0.69 and 0.0024, respectively. N. incisum populations from the southeast part of the QTP (E, G, H, J, U, W, X, HB and HC) had the highest haplotype diversity in this species, and haplotypes H1 and H7 had the highest distribution frequencies. N. franchetii populations from KG, KI, KJ, KK, KP, YC, YD, YE, and YF had the highest haplotype diversity in this species, and haplotype H22 had the highest frequency. In addition, N. oviforme and N. forrestii exhibited low haplotype diversity in terms of their ITS sequences (Table 2 and Supplementary Figure S1).

Genetic Diversity and Structure

The total genetic diversity (h_T) values based on the cpDNA datasets for N. incisum, N. franchetii, N. oviforme, and N. forrestii were 0.939, 0.766, 0.961, and 0.623, respectively, where N. incisum had the highest levels for h_S (0.404) and h_T, whereas N. forrestii had the lowest level of diversity (h_S = 0.167; h_T = 0.623) (Table 3). In addition, we calculated the genetic differentiation coefficients G_ST and N_ST for the four species. The U statistic (Gaussian test 1000 times) showed that N_ST was significantly larger than G_ST for N. incisum and N. oviforme (P < 0.05), thereby indicating that these two species exhibited significant phylogeographic structuring (Table 3).

AMOVA analysis of the cpDNA datasets detected genetic variations among the four species (F_CT = 0.5804) (Supplementary Table S5). In the four individual species, the genetic variations among populations (N. incisum, F_ST = 0.8196; N. franchetii, F_ST = 0.8391; N. oviforme, F_ST = 0.8474; and N. forrestii, F_ST = 0.7585) were significantly higher than those within populations (Supplementary Table S5). In addition, the AMOVA results obtained for the ITS sequences indicated similar genetic differentiation patterns to those based on the cpDNAs, where the differences among species in terms of the variation in the ITS were as high as 92% (F_CT = 0.9287) (Supplementary Table S6).

Phylogenetic Relationships

Phylogenetic trees of the cpDNA haplotypes were constructed based on the ML, MP, and Bayesian inference methods, which showed that the topological structures obtained were basically the same using the three methods (Figure 3). The four species of Notopterygium formed a larger monophyletic clade with high bootstrap support, where N. franchetii and N. oviforme were sisters. The median-joining network diagram produced using the cpDNA datasets was consistent with the phylogenetic analysis (Figure 2). The major haplotypes with the highest distribution frequencies (H1, H11, H32, H33, H41, and H53) were located in the central positions of the network. However, the phylogenetic relationships of the ITS sequences differed from those of the cpDNA sequences. No haplotypes were shared among the four Notopterygium species and each species formed its own individual branch (Figure 3) in the ITS tree. In addition, in order to confirm the phylogenetic positions of all four species considered in this study, we analyzed the other two species in the genus Notopterygium, i.e., N. tenuifolium and N. pinnatiinvolucellatum. Phylogenetic analyses based on variations in the chloroplast rbcL sequence showed that the four species considered in this study, i.e., N. incisum, N. franchetii, N. oviforme, and N. forrestii, were more closely related than N. tenuifolium and N. pinnatiinvolucellatum (Supplementary Figure S2).

FIGURE 3

Population Dynamics History and Divergence Time

Based on the cpDNA sequences, we performed various mathematical analyses to determine the population histories of the four Notopterygium species (Table 4 and Supplementary Figure S3). The mismatch distribution model had a single peak, with negative Tajima’s D and Fu’s F_S values for N. incisum and N. franchetii, which suggested that these two species had experienced rapid range expansions. The larger population growth indexes for N. incisum (g = 809) and N. franchetii (g = 2810.736) were also consistent with rapid population expansions. By contrast, N. oviforme and N. forrestii had bimodal mismatch distributions with positive Tajima’s D and Fu’s F_S values, where these results indicated that they did not experience expansion events. Therefore, we estimated the expansion times for N. incisum and N. franchetii as about 128–43 Kya and 51–17 Kya in the Pleistocene, respectively (Table 5).

We estimated the divergence times between the four species of Notopterygium based on a range of mutation rates (1.0–3.0 × 10^-9 s/s/y). The first divergence among the four species occurred between approximately 3.6 Mya (95% highest posterior density (HPD), 2.1–5.3 Mya) and 1.2 Mya (95% HPD, 0.67–1.8 Mya), whereas the estimated divergence between N. forrestii and N. incisum occurred between 2.24 Mya (95% HPD, 1.14–3.4 Mya) and 0.75 Mya (95% HPD, 0.4–1.14 Mya). In addition, the divergence between the major lineages of N. franchetii and N. oviforme occurred between 1.3 Mya (95% HPD, 0.54–2.2 Mya) and 0.42 Mya (95% HPD, 0.18–0.73 Mya) (Figure 4).

FIGURE 4

Species Distribution Modeling

In this study, MaxEnt modeling had the highest predictive capacity (AUC > 0.9) for the two widely distributed Notopterygium species (N. franchetii and N. incisum). The distribution ranges predicted for these two species were consistent with the current geographic distributions in the QTP and adjacent areas (Figures 5, 6 and Supplementary Table S7). Species distribution modeling also showed that the range of N. incisum was limited in the LIG period whereas it expanded very rapidly in the LGM period. However, there were no significant changes in the distribution range from the LGM until the current period for this species. For N. franchetii, MaxEnt modeling suggested that the distribution range of this species increased very rapidly from the LIG until the LGM period. However, it was interesting that the distribution range of N. franchetii did not change greatly from the LGM until the current period.

FIGURE 5

FIGURE 6

Discussion

Genetic Diversity and Structure

In the current study, our analysis of the cpDNA sequences showed that N. oviforme had the highest level of genetic diversity (H_d = 0.81, π = 0.0013), followed by N. incisum (H_d = 0.75, π = 0.00086) and N. forrestii (H_d = 0.39, π = 0.0002), whereas N. franchetii had the lowest level of diversity (H_d = 0.29, π = 0.00031) (Table 2). However, the results were different according to the ITS sequence analysis, where N. incisum had the highest diversity (H_d = 0.71, π = 0.0025), followed by N. oviforme (H_d = 0.69, π = 0.0024) and N. franchetii (H_d = 0.55, π = 0.0036), whereas N. forrestii exhibited no variation (Table 2). In addition, all four Notopterygium species had high levels of genetic differentiation, where the genetic variations in the cpDNA and ITS sequences mainly occurred among the populations within each species (Table 3 and Supplementary Tables S5, S6).

In general, N. franchetii has the most extensive natural geographic distribution range, but we found that its cpDNA and ITS sequences had low diversity. We consider that this low diversity may be due to harvesting and climatic changes, where many of the natural populations of N. franchetii have become extinct because of habitat destruction, thereby causing low diversity and high genetic differentiation (Lowe et al., 2005). N. incisum is another widely distributed species but we found that it had a high level of genetic diversity compared with other three species, which may be explained by the less extensive destruction of the wild populations of this species. According to the field investigations, we found that this species generally occurs in higher altitude areas (≥3000 m) compared with other Notopterygium species, and thus its less frequent harvesting might explain the high genetic variation (Wang G.N. et al., 2014). In addition, the high level of genetic diversity in N. oviforme according to this study might be explained by the lower altitude range of this species (1700–3200 m), which is consistent with a previous report of high species diversity at low altitudes (Lu et al., 2012). The lower genetic diversity of N. forrestii may be due to its narrow geographical distribution, where the smaller localized populations can interbreed and the gene flow is greater, thereby leading to a low level of diversity.

These Notopterygium species may also have been affected by adverse environmental changes in the high altitude QTP and adjacent areas. Thus, repeated climatic oscillations and geological events may have led to genetic drift and the fragmentation of habitats, thereby reducing their diversity (N. oviforme had slightly higher diversity compared with the other three species) and causing a high level of genetic differentiation among the populations of Notopterygium species (Hamrick and Loveless, 1989).

Relationships among Species

Phylogenetic analysis based on the cpDNA and ITS haplotypes showed that all four Notopterygium species formed a monophyletic clade with high bootstrap support (Figure 3 and Supplementary Figure S2). N. franchetii and N. oviforme shared a common branch in the phylogenetic tree based on the cpDNA sequences, and they also shared cpDNA haplotype H32. However, there were no shared ITS haplotypes among the four Notopterygium species, where each species formed an individual clade in the ITS tree. Thus, the ITS marker could identify the species at greater resolution than the cpDNA fragments. In general, cpDNA is a uniparentally inherited region whereas nuclear ITS fragments are biparentally inherited markers in most angiosperms, so ITS markers are superior for discriminating lineages and species than cpDNA fragments (Fan et al., 2011; Wang et al., 2012).

In addition, the shared cpDNA H32 haplotype was found in the parapatric populations (LF and YA) of N. franchetii and N. oviforme. These parapatric geographic distributions may have provided the opportunity for interspecific gene flow and hybridization among the two species. According to the field observations, we found that these two species have overlapping flowering times, which may have facilitated genetic introgression among these species. Previous studies have also suggested the occurrence of hybridization among species distributed in the same geographic regions and subsequent backcrosses with one of the parental species, where these processes resulted in high levels of shared plastid genotypes (Hamzeh et al., 2006; Li et al., 2013; Wang Z. et al., 2014). However, it is also possible that incomplete lineage sorting could have lead to the sharing of cpDNA haplotypes among species. The perennial herb Notopterygium species have large population sizes and long generation times, which may have led to the sharing of ancestral polymorphisms among species.

Species Divergence and Population Dynamics History

Mountain barriers may play a key role in speciation and diversification because their topographic complexity can lead to ecological stratification and environmental heterogeneity (Fjeldså et al., 2012). In the present study, we estimated the divergence time of the four Notopterygium species based on three cpDNA fragments, which showed that their divergence occurred between about 3.6 Mya (95% HPD, 2.1–5.3 Mya) and 1.2 Mya (95% HPD, 0.67–1.8 Mya) in the Pliocene and Pleistocene periods. The divergence of N. forrestii and N. incisum was estimated as occurring between 2.24 Mya (95% HPD, 1.14–3.4 Mya) and 0.75 Mya (95% HPD, 0.4–1.14 Mya) in the early to middle Pleistocene period (Figure 4). We suggest that the divergence of the four Notopterygium species was significantly related to the uplift of the QTP. Previous studies and geological data indicated that the uplift of the QTP started in the Oligocene to Miocene (25–17 Mya), middle of the Miocene (15–13 Mya), late Miocene (8–7 Mya), or in the Pliocene to early Pleistocene period (3.6–1.8 Mya) (Harrison et al., 1992; Coleman and Hodges, 1995; Shi et al., 1998; Spicer et al., 2003). During the uplift of the QTP and adjacent Himalayan mountains, long-term geological events generated great environmental differences, which might have triggered the diversification of species in the genus Notopterygium. Other studies have also shown that the recent extensive uplift of the QTP and adjacent mountains triggered the lineage divergence and evolution of many herb species due to geographical isolation and climatic changes (Lei et al., 2007, 2015).

In addition, dramatic variations in the environment and climate might have affected the genetic structure and geographic distributions of the Notopterygium species. MaxEnt modeling showed that the two cold-tolerant species comprising N. incisum and N. franchetii exhibited significant range expansions from the LIG to the LGM period (Figures 5, 6). The mismatch analysis, neutrality test, and population growth index results also supported similar expansions by these two species. Therefore, we estimated the expansion times for N. incisum and N. franchetii as about 128–43 and 51–17 Kya, respectively, during the late Pleistocene (Table 5). We showed that the geographic ranges of these two species increased significantly during the ice ages in the Pleistocene. Demographic expansions of cold-tolerant tree species during the glacial periods have also been reported in high altitude areas of the QTP (Qiu et al., 2011; Wen et al., 2014). Moreover, repeated founder and bottleneck effects during the expansion processes may explain the low genetic variation in the two species. By contrast, we found that the populations of N. incisum and N. franchetii had high levels of genetic diversity in the southeast part of the QTP. For example, some N. incisum populations in Qinghai (population G), Sichuan (populations J, Q, U, and X), and Gansu (population Z) had high diversity and many more unique haplotypes. The N. franchetii populations in Sichuan (populations KG and YF) and Gansu (populations KO and KZ) also had high genetic diversity and a rich abundance of haplotypes (Table 2). In addition, the LE population of N. oviforme and the LCA population of N. forrestii had high levels of haplotype diversity. These areas may have provided important glacial refugia for these endemic perennial herb species. Mountain areas at low latitudes can also provide relatively stable environmental conditions according to the “ecological stability hypothesis” (Qu et al., 2014; Lei et al., 2015), which implies that these populations should have high genetic diversity and a rich diversity of haplotypes (Tzedakis et al., 2002; Abbott and Brochmann, 2003; Petit et al., 2003; Schonswetter et al., 2005). Similar results have been obtained for other organisms, such as birds (Lei et al., 2007), mammals (Wu et al., 2005), spiders (Meng et al., 2008), and aphids (Huang et al., 2010).

Conservation Strategies for Endangered Notopterygium Species

The genus Notopterygium comprises unique perennial herbaceous plants with medicinal applications in China (Zhou et al., 2010). These species have high economic value so the market demand is great, especially for N. incisum and N. franchetii. However, in recent years, due to their continuous harvesting, slow growth rate, and low reproductive capacity, the natural populations of these Notopterygium species have been greatly depleted (Zhou et al., 2010). According to field investigations, we found that many of the previously recorded natural populations of species in the genus Notopterygium were extinct, and thus these important resources require urgent conservation and management.

According to the results of our population genetics analysis, we propose that the natural populations of wild Notopterygium species should be protected in situ, especially in the natural refugia areas (i.e., populations J, Q, U, and X of N. incisum; populations KO and KZ of N. franchetii; population LE of N. oviforme; and population LCA of N. forrestii). In addition, it is necessary to control all activities that deplete the sizes of the populations (e.g., illegal harvesting) and genetic fragmentations (e.g., habitat loss) (Hamilton et al., 2012).

In order to conserve the populations of a species, it is necessary to understand the genetic diversity and population structure of the natural populations (Schmitt, 2007). In this study, we found that the genetic variability and haplotype diversity were low for the widely distributed species, thereby indicating that the habitats have been destroyed or fragmented for these species, where interbreeding has occurred with nearby populations of individuals, thereby reducing the haplotype diversity. It is also necessary to protect the populations in different regions in order to increase the genetic links among populations. Finally, the mature seeds from each population should be collected and artificially planted with other populations in order to improve the habitats and to strengthen the gene exchange among populations (Cabrera-Toledo et al., 2012).

References

• 1

  AbbottR. J.BrochmannC. (2003). History and evolution of the arctic flora: in the footsteps of Eric Hultén.Mol. Ecol.12299–313. 10.1046/j.1365-294X.2003.01731.x

  • CrossRef
  • Google Scholar

• 2

  AllenG. A.MarrK. L.McCormickL. J.HebdaR. J. (2015). Geographical origins, migration patterns and refugia of Sibbaldia procumbens, an arctic-alpine plant with a fragmented range.J. Biogeogr.421665–1676. 10.1111/jbi.12543

  • CrossRef
  • Google Scholar

• 3

  AndersonL. L.HuF. S.NelsonD. M.PetitR. J.PaigeK. N. (2006). Ice-age endurance: DNA evidence of a white spruce refugium in Alaska.Proc. Natl. Acad. Sci. U.S.A.10312447–12450. 10.1073/pnas.0605310103

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  AviseJ. C. (2009). Phylogeography: retrospect and prospect.J. Biogeogr.363–15. 10.1111/j.1365-2699.2008.02032.x

  • CrossRef
  • Google Scholar

• 5

  BaderE. (2005). ArcGIS 9: ArcGIS Server Administrator and Developer Guide.Redlands, CA: Esri Press.

  • Google Scholar

• 6

  BaiW. N.LiaoW. J.ZhangD. Y. (2010). Nuclear and chloroplast DNA phylogeography reveal two refuge areas with asymmetrical gene flow in a temperate walnut tree from East Asia.New Phytol.188892–901. 10.1111/j.1469-8137.2010.03407.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  BirksH. J. B.WillisK. J. (2008). Alpines, trees, and refugia in Europe.Plant Ecol. Divers.1147–160. 10.1080/17550870802349146

  • CrossRef
  • Google Scholar

• 8

  Cabrera-ToledoD.González-AstorgaJ.Flores-VázquezJ. C. (2012). Fine-scale spatial genetic structure in two Mexican Cycads species Dioon caputoi and Dioon merolae (Zamiaceae, Cycadales): implications for conservation.Biochem. Syst. Ecol.4043–48. 10.1016/j.bse.2011.09.004

  • CrossRef
  • Google Scholar

• 9

  Cavender-BaresJ.González-RodríguezA.EatonD. A.HippA. A.BeulkeA.ManosP. S. (2015). Phylogeny and biogeography of the American live oaks (Quercus subsection Virentes): a genomic and population genetics approach.Mol. Ecol.243668–3687. 10.1111/mec.13269

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10

  CBOL Plant Working Group (2009). A DNA barcode for land plants.Proc. Natl. Acad. Sci. U.S.A.10612794–12797. 10.1073/pnas.0905845106

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11

  ColemanM.HodgesK. (1995). Evidence for Tibetan plateau uplift before 14 Myr ago from a new minimum age for east-west extension.Nature37449–52. 10.1073/pnas.0605310103

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12

  DoyleJ. J.DoyleJ. L. (1987). A rapid DNA isolation procedure for small quantities of fresh leaf material.Phytochem. Bull.1911–15.

  • Google Scholar

• 13

  DrummondA. J.RambautA. (2007). BEAST: Bayesian evolutionary analysis by sampling Trees.BMC Evol. Biol.7:214. 10.1186/1471-2148-7-214

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14

  ExcoffierL.LischerH. E. (2010). Arlequin suite v 3.5: a new series of programs to perform population genetics analyses under Linux and Windows.Mol. Ecol. Resour.10564–567. 10.1111/j.1755-0998.2010.02847.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15

  FanW. Y.ZhouC. J.WangY. C. (2011). Confirmation of natural hybrids between two species of Nitraria in Inner Mongolia based on molecular evidence.Acta Bot. Boreali Occident. Sin.5913–919.

  • Google Scholar

• 16

  FavreA.PäckertM.PaulsS. U.JähnigS. C.UhlD.MichalakI.et al (2015). The role of the uplift of the Qinghai-Tibetan plateau for the evolution of Tibetan biotas.Biol. Rev.90236–253. 10.1111/brv.12107

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17

  FawcettT. (2006). An introduction to ROC analysis.Pattern Recogn. Lett.27861–874. 10.1016/j.patrec.2005.10.010

  • CrossRef
  • Google Scholar

• 18

  FjeldsåJ.BowieR. C. K.RahbekC. (2012). The role of mountain ranges in the diversification of birds.Annu. Rev. Ecol. Evol. Syst.43249–256. 10.1146/annurev-ecolsys-102710-145113

  • CrossRef
  • Google Scholar

• 19

  FuY. X. (1997). Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection.Genetics147915–925.

  • Pubmed Abstract
  • Google Scholar

• 20

  FuY. X.LiW. H. (1993). Statistical tests of neutrality of mutations.Genetics133693–709.

  • Google Scholar

• 21

  GaoJ.WangB.MaoJ. F.IngvarssonP.ZengQ. Y.WangX. R. (2012). Demography and speciation history of the homoploid hybrid pine Pinus densata on the Tibetan Plateau.Mol. Ecol.214811–4827. 10.1111/j.1365-294X.2012.05712.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22

  GrivetD.PetitR. J. (2002). Phylogeography of the common ivy (Hedera sp.) in Europe: genetic differentiation through space and time.Mol. Ecol.111351–1362. 10.1046/j.1365-294X.2002.01522.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23

  GuichouxE.Garnier-GéréP.LagacheL.LangT.BouryC.PetitR. J. (2013). Outlier loci highlight the direction of introgression in Oaks.Mol. Ecol.22450–462. 10.1111/mec.12125

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  HallT. A. (1999). BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT.Nucleic Acids Symp. Ser.4195–98.

  • Google Scholar

• 25

  HamiltonA.PeiS.HuaiH.AndersonS. (2012). Why and how to make plant conservation ecosystem-based.Sustain. Agric. Res.148–54. 10.5539/sar.v1n1p48

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  HamrickJ. L.LovelessM. D. (1989). “The genetic structure of tropical tree populations: associations with reproductive biology,” in The Evolutionary Ecology of Plants, edsBockJ.LinhartY. B. (Boulder, CO: Westview Press), 129–146.

  • Google Scholar

• 27

  HamzehM.PérinetP.DayanandanS. (2006). Genetic relationships among species of Populus (Salicaceae) based on nuclear genomic data 1.J. Torrey Bot. Soc.133519–527. 10.3159/1095-5674(2006)133[519:GRASOP]2.0.CO;2

  • CrossRef
  • Google Scholar

• 28

  HarpendingH. C. (1994). Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution.Hum. Biol.66591–600.

  • Pubmed Abstract
  • Google Scholar

• 29

  HarrisonT. M.CopelandP.KiddW. S. F.YinA. (1992). Raising Tibet.Science2551663–1670. 10.1126/science.255.5052.1663

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  HewittG. (2004). Genetic consequences of climatic oscillations in the Quaternary.Philos. Trans. R. Soc. Lond. Biol. Sci.359183–195. 10.1098/rstb.2003.1388

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  HickersonM. J.CarstensB. C.Cavender-BaresJ.CrandallK. A.GrahamC. H.JohnsonJ. B.et al (2010). Phylogeography’s past, present, and future: 10 years after.Mol. Phylogenet. Evol.54291–301. 10.1016/j.ympev.2009.09.016

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  HuangX. L.QiaoG. X.LeiF. M. (2010). Use of parsimony analysis to identify areas of endemism of Chinese birds: implications for conservation and biogeography.Int. J. Mol. Sci.112097–2108. 10.3390/ijms11052097

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  HuelsenbeckJ. P.RonquistF. (2001). MRBAYES: Bayesian inference of phylogenetic tree.Bioinformatics17754–755. 10.1093/bioinformatics/17.8.754

  • CrossRef
  • Google Scholar

• 34

  HughesC. E.AtchisonG. W. (2015). The ubiquity of alpine plant radiations: from the Andes to the Hengduan Mountains.New Phytol.207275–282. 10.1111/nph.13230

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  Ickert-BondS. M.RennerS. S. (2016). The Gnetales: recent insights on their morphology, reproductive biology, chromosome numbers, biogeography, and divergence times.J. Syst. Evol.541–16. 10.1111/jse.12190

  • CrossRef
  • Google Scholar

• 36

  JiaD. R.AbbottR. J.LiuT. L.MaoK. S.BartishI. V.LiuJ. Q. (2012). Out of the Qinghai-Tibet Plateau: evidence for the origin and dispersal of Eurasian temperate plants from a phylogeographic study of Hippophaë rhamnoides (Elaeagnaceae).New Phytol.1941123–1133. 10.1111/j.1469-8137.2012.04115.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  JiaY.LiuM. L.YueM.ZhaoZ.ZhaoG. F.LiZ. H. (2017). Comparative transcriptome analysis reveals adaptive evolution of Notopterygium incisum and Notopterygium franchetii, two high-alpine herbal species endemic to China.Molecules22:1158. 10.3390/molecules22071158

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38

  JiangS. Y.SunH.HuangX. J.ZhouY.MaX. J.YangZ. R. (2005). Environmental pedology of Notopterygium incisum and N. forbesii.Chin. Tradit. Herb. Drugs36918–921.

  • Google Scholar

• 39

  KeppelG.Van NielK. P.Wardell-JohnsonG. W.YatesC. J.ByrneM.MucinaL.et al (2012). Refugia: identifying and understanding safe havens for biodiversity under climate change.Glob. Ecol. Biogeogr.21393–404. 10.1111/j.1466-8238.2011.00686.x

  • CrossRef
  • Google Scholar

• 40

  KuhnerM. K. (2006). LAMARC 2.0: maximum likelihood and Bayesian estimation of population parameters.Bioinformatics22768–770. 10.1093/bioinformatics/btk051

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  LeiF. M.QuY. H.SongG.AlströmP.FjeldsåJ. (2015). The potential drivers in forming avian biodiversity hotpots in the east Himalaya-Mountains of southwest China.Integr. Zool.10171–181. 10.1111/1749-4877.12121

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  LeiF. M.WeiG. A.ZhaoH. F.YinZ. H.LuJ. L. (2007). China subregional avian endemism and biodiversity conservation.Biodivers. Conserv.161119–1130. 10.1007/s10531-006-9080-3

  • CrossRef
  • Google Scholar

• 43

  LiJ. J.FangX. M. (1999). Uplift of the Tibetan Plateau and environmental changes.China Sci. Bull.442117–2124. 10.1007/BF03182692

  • CrossRef
  • Google Scholar

• 44

  LiL.AbbottR. J.LiuB.SunY.LiL.ZouJ. B.et al (2013). Pliocene intraspecific divergence and Plio-Pleistocene range expansions within Picea likiangensis (Lijiang spruce), a dominant forest tree of the Qinghai-Tibet Plateau.Mol. Ecol.225237–5255. 10.1111/mec.12466

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45

  LiaoY. Y.GichiraA. W.WangQ. F.ChenJ. M. (2015). Molecular phylogeography of four endemic Sagittaria species (Alismataceae) in the Sino Japanese Floristic Region of East Asia.Bot. J. Linn. Soc.1806–20. 10.1111/boj.12351

  • CrossRef
  • Google Scholar

• 46

  LibradoP.RozasJ. (2009). DnaSPv5: a software for comprehensive analysis of DNA polymorphism data.Bioinformatics251451–1452. 10.1093/bioinformatics/btp187

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47

  LiuJ.M€ollerM.ProvanJ.GaoL. M.PoudelR. C.LiD. Z. (2013). Geological and ecological factors drive cryptic speciation of yews in a biodiversity hotspot.New Phytol.1991093–1108. 10.1111/nph.12336

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  LiuJ. Q.SunY. S.GeX. J.GaoL. M.QiuY. X. (2012). Phylogeographic studies of plants in China: advances in the past and directions in the future.J. Syst. Evol.50267–275. 10.1111/j.1759-6831.2012.00214.x

  • CrossRef
  • Google Scholar

• 49

  LiuJ. Q.WangY. J.WangA. L.HideakiO.AbbottR. J. (2006). Radiation and diversification within the Ligularia-Cremanthodium-Parasenecio complex (Asteraceae) triggered by uplift of the Qinghai-Tibetan Plateau.Mol. Phylogenet. Evol.3831–49. 10.1016/j.ympev.2005.09.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  LoweA. J.BoshierD.WardM.BaclesC. F.NavarroC. (2005). Genetic resource impacts of habitat loss and degradation; reconciling empirical evidence and predicted theory for neotropical trees.Heredity95255–273. 10.1038/sj.hdy.6800725

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51

  LuB.ZhengY.MurphyR. W.ZengX. (2012). Coalescence patterns of endemic Tibetan species of stream salamanders (Hynobiidae: Batrachuperus).Mol. Ecol.213308–3324. 10.1111/j.1365-294X.2012.05606.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  MaoK. S.HaoG.LiuJ. Q.AdamsR. P.MilneR. I. (2010). Diversification and biogeography of Juniperus (Cupressaceae): variable diversification rates and multiple intercontinental dispersals.New Phytol.188254–272. 10.1111/j.1469-8137.2010.03351.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53

  McLachlanJ. S.ClarkJ. S.ManosP. S. (2005). Molecular indicators of tree migration capacity under rapid climate change.Ecology862088–2098. 10.1890/04-1036

  • CrossRef
  • Google Scholar

• 54

  MengK.LiS. Q.MurphyR. W. (2008). Biogeographical patterns of Chinese spiders (Arachnida: Araneae) based on a parsimony analysis of endemicity.J. Biogeogr.351241–1249. 10.1111/j.1365-2699.2007.01843.x

  • CrossRef
  • Google Scholar

• 55

  MittermeierR. A.GilP. R.HoffmanM.PilgrimJ.BrooksT.MittermeierC. G.et al (2005). Hotspots Revisited: Earth’s Biologically Richest and Most Endangered Terrestrial Ecoregions.Chicago, IL: University of Chicago Press.

  • Google Scholar

• 56

  NeiM.LiW. H. (1979). Mathematical model for studying genetic variation in terms of restriction endonucleases.Proc. Natl. Acad. Sci. U.S.A.765269–5273. 10.1073/pnas.76.10.5269

  • CrossRef
  • Google Scholar

• 57

  NeiM.TajimaF. (1981). DNA polymorphism detectable by restriction endonucleases.Genetics97145–163.

  • Google Scholar

• 58

  NybomH. (2004). Comparison of different nuclear DNA markers for estimating intraspecific genetic diversity in plants.Mol. Ecol.131143–1155. 10.1111/j.1365-294X.2004.02141.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59

  OhsawaT.IdeY. (2008). Global patterns of genetic variation in plant species along vertical and horizontal gradients on mountains.Glob. Ecol. Biogeogr.17152–163. 10.1111/j.1466-8238.2007.00357.x

  • CrossRef
  • Google Scholar

• 60

  OpgenoorthL.VendraminG. G.MaoK.MieheG.MieheS.LiepeltS.et al (2010). Tree endurance on the Tibetan Plateau marks the world’s highest known tree line of the Last Glacial Maximum.New Phytol.185332–342. 10.1111/j.1469-8137.2009.03007.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61

  OrtegoJ.NogueralesV.GuggerP. F.SorkV. L. (2015). Evolutionary and demographic history of the Californian scrub white Oak species complex: an integrative approach.Mol. Ecol.246188–6208. 10.1111/mec.13457

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 62

  OrtegoJ.RiordanE. C.GuggerP. F.SorkV. L. (2012). Influence of environmental heterogeneity on genetic diversity and structure in an endemic southern Californian Oak.Mol. Ecol.213210–3223. 10.1111/j.1365-294X.2012.05591.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 63

  ParducciL.JørgensenT.TollefsrudM. M.ElverlandE.AlmT.FontanaS. L.et al (2012). Glacial survival of boreal trees in northern Scandinavia.Science3351083–1086. 10.1126/science.1216043

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64

  PetitR. J.AguinagaldeI.de BeaulieuJ. L.BittkauC.BrewerS.CheddadiR.et al (2003). Glacial refugia: hotspots but not melting pots of genetic diversity.Science3001563–1565. 10.1126/science.1083264

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 65

  PhillipsS. J.AndersonR. P.SchapireR. E. (2006). Maximum entropy modeling of species geographic distributions.Ecol. Model.190231–259. 10.1016/j.ecolmodel.2005.03.026

  • CrossRef
  • Google Scholar

• 66

  PhillipsS. J.DudíkM. (2008). Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation.Ecography31161–175. 10.1111/j.0906-7590.2008.5203.x

  • CrossRef
  • Google Scholar

• 67

  PolzinT.DaneshmandS. V. (2003). On Steiner trees and minimum spanning trees in hypergraphs.Oper. Res. Lett.3112–20. 10.1016/S0167-6377(02)00185-2

  • CrossRef
  • Google Scholar

• 68

  PosadaD.CrandallK. A. (1998). Modeltest: testing the model of DNA substitution.Bioinformatics14817–818. 10.1093/bioinformatics/14.9.817

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 69

  ProvanJ.BennettK. D. (2008). Phylogeographic insights into cryptic glacial refugia.Trends Ecol. Evol.23564–571. 10.1016/j.tree.2008.06.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 70

  PuF.WangP.ZhengZ.WangY. P. (2000). A reclassification of Notopterygium boissieu (Umbelliferae).Acta Phytotaxon. Sin.38430–436.

  • Google Scholar

• 71

  QiuY. X.FuC. X.ComesH. P. (2011). Plant molecular phylogeography in China and adjacent regions: tracing the genetic imprints of Quaternary climate and environmental change in the world’s most diverse temperate flora.Mol. Phylogenet. Evol.59225–244. 10.1016/j.ympev.2011.01.012

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 72

  QuY. H.EricsonP. G.QuanQ.SongG.ZhangR. Y.GaoB.et al (2014). Long-term isolation and stability explain high genetic diversity in the Eastern Himalaya.Mol. Ecol.23705–720. 10.1111/mec.12619

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73

  RambautA. (2009). FigTree v.1.3.1. Available at: http://tree.bio.ed.ac.uk/software/figtree/

  • Google Scholar

• 74

  RogersA. R.HarpendingH. (1992). Population growth makes waves in the distribution of pairwise genetic differences.Mol. Biol. Evol.9552–569.

  • Pubmed Abstract
  • Google Scholar

• 75

  SchmittT. (2007). Molecular biogeography of Europe: Pleistocene cycles and postglacial trends.Front. Zool.4:11. 10.1186/1742-9994-4-11

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 76

  SchneiderS.ExcoffierL. (1999). Estimation of past demographic parameters from the distribution of pairwise differences when the mutation rates vary among sites: application to human mitochondrial DNA.Genetics1521079–1089.

  • Pubmed Abstract
  • Google Scholar

• 77

  SchonswetterP.StehlikI.HoldereggerR.TribschA. (2005). Molecular evidence for glacial refugia of mountain plants in the European Alps.Mol. Ecol.143547–3555. 10.1111/j.1365-294X.2005.02683.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 78

  SheM.PuF. (1996). A new species of Notopterygium de Boiss. from China.J. Plant Resour. Environ.641–42.

  • Google Scholar

• 79

  ShiY. F.LiJ. J.LiB. Y. (1998). Uplift and Environmental Changes on Qinghai-Xizang (Tibetan) Plateau in the Late Cenozoic.Guangzhou: Guangdong Science and Technology Press, 1–289.

  • Google Scholar

• 80

  SoltisD. E.MorrisA. B.McLachlanJ. S.ManosP. S.SoltisP. S. (2006). Comparative phylogeography of unglaciated eastern North America.Mol. Ecol.154261–4293. 10.1111/j.1365-294X.2006.03061.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 81

  SoumayaR. C.SarraC.SalwaZ. A.KhaledC.KhaledS. (2014). Molecular polymorphism and phylogenetic relationships within Tunisian date palm (Phoenix dactylifera L.): evidence of non-coding trnL-trnF regions of chloroplast DNAs.Sci. Hortic.17032–38. 10.1016/j.scienta.2014.02.027

  • CrossRef
  • Google Scholar

• 82

  SpicerR. A.HarrisN. B. W.WiddowsonM.HermanA. B.GuoS.ValdesP. J.et al (2003). Constant elevation of southern Tibet over the past 15 million years.Nature421622–624. 10.1038/nature01356

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 83

  StewartJ. R.ListerA. M.BarnesI.DalénL. (2010). Refugia revisited: individualistic responses of species in space and time.Proc. Biol. Sci.277661–671. 10.1098/rspb.2009.1272

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 84

  StewartJ. R.StringerC. B. (2012). Human evolution out of Africa: the role of refugia and climate change.Science3351317–1321. 10.1126/science.1215627

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 85

  SudhirK.MasatoshiN.JoelD.TamuraK. (2008). MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequence.Brief. Bioinform.9299–306. 10.1093/bib/bbn017

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 86

  SunY.AbbottR. J.LiL.LiL.ZhouJ. B.LiuJ. Q. (2014). Evolutionary history of Purple cone spruce (Picea purpurea) in the Qinghai-Tibet Plateau: homoploid hybrid origin and Pleistocene expansion.Mol. Ecol.23343–359. 10.1111/mec.12599

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 87

  TajimaF. (1989). Statistical method for testing the neutral mutation hypothesis by DNA polymorphism.Genetics123585–595.

  • Pubmed Abstract
  • Google Scholar

• 88

  TianX. M.LuoJ.WangA. L.MaoK. S.LiuJ. Q. (2011). On the origin of the woody buckwheat Fagopyrum tibeticum (Parapteropyrum tibeticum) in the Qinghai-Tibetan Plateau.Mol. Phylogenet. Evol.61515–520. 10.1016/j.ympev.2011.07.001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 89

  TzedakisP. C.FrogleyM. R.HeatonT. H. E. (2002). Duration of last interglacial conditions in northwestern Greece.Quat. Res.5853–55. 10.1006/qres.2002.2328

  • CrossRef
  • Google Scholar

• 90

  WangG. N.HeX. Y.MieheG.MaoK. S. (2014). Phylogeography of the Qinghai-Tibet Plateau endemic alpine herb Pomatosace filicula (Primulaceae).J. Syst. Evol.52289–302. 10.1111/jse.12089

  • CrossRef
  • Google Scholar

• 91

  WangH.QiongL. A.SunK.LuF.WangY.SongZ. P.et al (2010). Phylogeographic structure of Hippophae tibetana (Elaeagnaceae) highlights the highest microrefugia and the rapid uplift of the Qinghai-Tibetan Plateau.Mol. Ecol.192964–2979. 10.1111/j.1365-294X.2010.04729.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 92

  WangS. X.LiQ. H.XuJ.GaoT. T.XinZ. M. (2012). Experimental research on the pollination characteristics of 4 plant species of genus Nitraria L.J. Biol.2949–51.

  • Google Scholar

• 93

  WangY. J.SusannaA.Von Raab-StraubeE.MilneR.LiuJ. Q. (2009). Island-like radiation of Saussurea (Asteraceae: Cardueae) triggered by uplifts of the Qinghai-Tibetan Plateau.Biol. J. Linn. Soc.97893–903. 10.1111/j.1095-8312.2009.01225.x

  • CrossRef
  • Google Scholar

• 94

  WangY. P.PuF.WangP.HeX. J. (1996). Studies on the systematics of the Chinese endemic genus Notopterygium.Acta Bot. Yunnanica18424–430.

  • Google Scholar

• 95

  WangZ.DuS.DayanandanS.WangD.ZengY.ZhangJ. (2014). Phylogeny reconstruction and hybrid analysis of Populus (Salicaceae) based on nucleotide sequences of multiple single-copy nuclear genes and plastid fragments.PLOS ONE9:e103645. 10.1371/journal.pone.0103645

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 96

  WangZ. W.ChenS. T.NieZ. L.ZhangJ. W.ZhouZ.DengT.et al (2015). Climatic factors drive population divergence and demography: insights based on the phylogeography of a riparian plant species endemic to the Hengduan Mountains and adjacent regions.PLOS ONE10:e0145014. 10.1371/journal.pone.0145014

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 97

  WenJ.ZhangJ. Q.NieZ. L.ZhongY.SunH. (2014). Evolutionary diversifications of plants on the Qinghai-Tibetan Plateau.Front. Genet.5:4. 10.3389/fgene.2014.00004

  • CrossRef
  • Google Scholar

• 98

  WolfeK. H.LiW. H.SharpP. M. (1987). Rates of nucleotide substitution vary greatly among plant mitochondrial, chloroplast, and nuclear DNAs.Proc. Natl. Acad. Sci. U.S.A.849054–9058. 10.1073/pnas.84.24.9054

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 99

  WuZ. Y.RavenP. H.HongD. Y. (eds). (2005). “Apiaceae through Ericaceae,” inFlora of China (Beijing: Science Press).

  • Google Scholar

• 100

  XuJ.DengM.JiangX. L.WestwoodM.SongY. G.TurkingtonR. (2015). Phylogeography of Quercus glauca (Fagaceae), a dominant tree of East Asian subtropical evergreen forests, based on three chloroplast DNA interspace sequences.Tree Genet. Genomes11:805. 10.1007/s11295-014-0805-2

  • CrossRef
  • Google Scholar

• 101

  XuT. T.AbbottR. J.MilneR. I.MaoK. S.DuF. K.WuG. L.et al (2010). Phylogeography and allopatric divergence of Cypress species (Cupressus L.) in the Qinghai-Tibetan Plateau and adjacent regions.BMC Evol. Biol.10:194. 10.1186/1471-2148-10-194

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 102

  YangJ.YueM.NiuC.MaX. F.LiZ. H. (2017). Comparative analysis of the complete chloroplast genome of four endangered herbals of Notopterygium.Genes8:E124. 10.3390/genes8040124

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 103

  ZhangD. F.FengquanL.JianminB. (2000). Eco-environmental effects of the Qinghai-Tibet Plateau uplift during the Quaternary in China.Environ. Geol.391352–1358. 10.1016/j.ympev.2015.10.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 104

  ZhangX. W.LiY.LiuC. Y.XiaT.ZhangQ.FangY. M. (2015). Phylogeography of the temperate tree species Quercus acutissima in China: inferences from chloroplast DNA variations.Biochem. Syst. Ecol.63190–197. 10.1016/j.bse.2015.10.010

  • CrossRef
  • Google Scholar

• 105

  ZhouG.YangL.LiC.XuW.ChenG. (2010). Genetic diversity in endangered Notopterygium forbesii Boissieu based on intraspecies sequence variation of chloroplast DNA and implications for conservation.Biochem. Syst. Ecol.38911–916. 10.1016/j.bse.2010.09.012

  • CrossRef
  • Google Scholar

• 106

  ZhouS. Z.WangX. L.WangJ.XuL. B. (2006). A preliminary study on timing of the oldest Pleistocene glaciation in Qinghai-Tibetan Plateau.Quat. Int.15444–51. 10.1016/j.quaint.2006.02.002

  • CrossRef
  • Google Scholar