Figures
Abstract
Global warming has created opportunities and challenges for the survival and development of species. Determining how climate change may impact multiple ecosystem levels and lead to various species adaptations is necessary for both biodiversity conservation and sustainable biological resource utilization. In this study, we employed Maxent to predict changes in the habitat range and altitude of Polygala tenuifolia Willd. under current and future climate scenarios in China. Four representative concentration pathways (RCP2.6, RCP4.5, RCP6.0, and RCP8.5) were modeled for two time periods (2050 and 2070). The model inputs included 732 presence points and nine sets of environmental variables under the current conditions and the four RCPs in 2050 and 2070. The area under the receiver-operating characteristic (ROC) curve (AUC) was used to evaluate model performance. All of the AUCs were greater than 0.80, thereby placing these models in the “very good” category. Using a jackknife analysis, the precipitation in the warmest quarter, annual mean temperature, and altitude were found to be the top three variables that affect the range of P. tenuifolia. Additionally, we found that the predicted highly suitable habitat was in reasonable agreement with its actual distribution. Furthermore, the highly suitable habitat area was slowly reduced over time.
Citation: Jiang H, Liu T, Li L, Zhao Y, Pei L, Zhao J (2016) Predicting the Potential Distribution of Polygala tenuifolia Willd. under Climate Change in China. PLoS ONE 11(9): e0163718. https://doi.org/10.1371/journal.pone.0163718
Editor: Daniel E. Naya, Universidad de la Republica Uruguay, URUGUAY
Received: April 17, 2016; Accepted: September 13, 2016; Published: September 23, 2016
Copyright: © 2016 Jiang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All occurrence files are available from the figshare database (https://dx.doi.org/10.6084/m9.figshare.3504812).
Funding: Our work was supported by three funds. The first one was Graduate Innovation Ability Training Project of Hebei Province, award no. was sj2015008, recipient was Hongjun Jiang. The second one was Special Fund on Public Health from the Central Finance "Resource survey and monitoring on original plant of essential medicine", award no. was Cai She[2011]76, recipient was Jiancheng Zhao. The last one was Special Fund of the National Traditional Chinese Medicine, award no. was 201207002, recipient was Jiancheng Zhao. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Radix Polygalae, called Yuanzhi in Chinese, is a traditional Chinese herb officially listed in the medical pharmacopoeia [1,2] and is used as a mucolytic, tonic, sedative, antipsychotic, and expectorant [1,3–5]. Various medicinal components are extracted from this herb. In addition to the various saponins, xanthones, and oligosaccharides recorded in the Pharmacopoeia Commission of the People's Republic of China 2005 [2], Radix polygalae contains flavonoids, coumarins, hydroxycinnamic acid conjugates, and lignans [6]. It is the dry root of Polygala tenuifolia Willd. or P. sibiric L., which belong to the Polygalaceae family [1], P. tenuifolia Willd. is the most important origin plant of Radix Polygalae.
It is a perennial herb mainly distributed in Gansu, Hebei, Heilongjiang, Henan, Jiangsu, Jiangxi, Liaoning, Inner Mongolia, Ningxia, Qinghai, Shaanxi, Shanxi, and Sichuan in China, and in Mongolia, the Russian Federation, and South Korea [7]. As the understanding of Traditional Chinese Medicine has improved, sales and consumption of Radix Polygalae have been growing in last few years. But its wild resource has been getting worsen and worsen by habitat destruction and over-exploitation[8]. And this slow-growing herb has a low yield and is often harvested by dredging before seed maturation. Thus, it has been included in the List of Grade III Key State-Protected Wild Medicinal Species in China[9]. Conflict between high consumption and yield decrement has gradually forced people to exploit more wild resources and cultivate them on a large scale in recent years. However, overexploitation and promotion of blind cultivation result in loss of a valuable wild resource and a decline in the quantity and quality of Radix Polygalae. Thus, determining the optimal habitat and adaptability of wild P. tenuifolia has great implications for its protection, recovery, and utilization.
Over the period from 1880 to 2012, average global temperatures have increased by 0.85°C. The period from 1983 to 2012 was very likely the warmest 30-year period in the last 800 years [10]. A great many shifts in the distribution and abundances of species occurred in this period [11,12]. Approximately 20% of all of the world’s plant species are on the brink of extinction [13]. Thomas et al. stated that 15–37% of species in their sample of regions and taxa will be extinct based on mid-range climate-warming scenarios by 2050 [14]. Survival of many species is threatened by many issues caused by global warming, such as rising sea levels, increased rainfall, increased drought, extreme weather, and climate events [10]. Particularly, those species with high economic values face over-harvesting and environmental disruption. Although evolution of species has proven that they can adapt to the changing environment [12], enormous challenges for wild species and their ecosystems have emerged with this rapid rate of global warming. Currently, species shift, contract, expand, or fragment their range in response to the changing climate [15]. These range shifts that track suitable climatic conditions have been observed in more than 1000 species, and habitat shifts appear to be one of the most important and effective methods by which species have adapted to changes in the environment [16]. Thus, understanding the impacts of climate change on range shift of different species is important for biodiversity conservation and biological resource protection [17,18]. The increasing work on species distribution models (SDMs) has provided a new approach for predicting those changes.
SDMs were developed in the mid-1980s [19]. These models involve the comprehensive application of geographic information systems (GISs), statistics, and ecology [20–22]. In recent years, these systems have been widely used in many aspects of resource management and conservation planning [23]. By associating species occurrence with environmental variables, SDMs can model or extrapolate the occurrence of a species and yield maps of habitat suitability or a predicted species distribution [20,24]. These results can be interpreted as the probability of species presence, habitat suitability, or species richness [23,25]. These maps and distributions are imperative for biodiversity assessment [26], reserve design [27], habitat management and restoration [28], distribution modeling [29], ecological restoration [30], invasive species risk assessment [31,32], agricultural disease and insect pest forecasting [33], and impact evaluation of environmental change [14,26,34]. A number of alternative SDMs are available to predict potential species distribution, including BIOCLIM [35], BIOMAPPER [36], BRT [37], CLIMEX [38], DOMAIN [39], Favorability Function [22], GAM [40], GLM [41], GARP [42], HABITAT [43], Maxent [44], and MARS [45]. However, these models have their advantages and disadvantages. Maxent has shown higher predictive accuracy than many other methods when applied to “presence-only” species occurrence data [20,25,44,46 47].
Maxent is a general purpose machine learning algorithm [44,48]. Currently, it is the most widely used SDM [49]. Since 2006, this algorithm has been used in several studies on plants and vertebrates (e.g. birds, reptiles) [26,48,50–55]. Maxent has been used in many reports on medicinal plants. Yang et al predicted the potential distribution of Justicia adhatoda L. in the Himalayan foothills [56]. Remya et al. found that the habitat suitable for the distribution of Myristica dactyloides on the Kolli Hill in India would significantly decrease by the years 2050 and 2070 [46]. Yi et al. showed that seven variables are dominant factors in determining the suitable habitat of Homonoia riparia Lour [57]. Zhang et al. predicted potential suitable cultivation regions and explored the key environmental factors that affect the content of active ingredients in Scutellaria baicalensis Georgi in China [58]. Maxent has been demonstrated to provide the most accurate predictions using presence-only data and exhibit a better overall performance than other methods that use both presence and absence data [48,59], and it is less sensitive to overfitting [60–62].
To seek suitable habitats and evaluate the impact of climate change on P. tenuifolia, our study used Maxent to model its potential distribution based on occurrence records and environmental variables (including soil, land cover, vegetation coverage, topographical variables, and bioclimatic variables). First, dominant environmental variables were selected to build a model; second, current bioclimatic variables were used to estimate the current climatic suitable habitat; finally, potential distributions under different future climate scenarios were predicted.
Materials and Methods
Species records
In this study, the total number of P. tenuifolia occurrences was 2332. Of these occurrence records, 2030 were collected from National Specimen Information Infrastructure (NSII; http://www.nsii.org.cn/) and Chinese Virtual Herbarium (CVH; http://www.cvh.org.cn/), 167 were collected from GBIF (http://www.gbif.org), and 135 were obtained from the fourth national survey on Chinese material medical resources. First, we georeferenced the occurrence records with detailed location information using GPS and Coordinate pick up system of Baidu Map (http://api.map.baidu.com/lbsapi/getpoint/index.html). Then, we removed 885 incomplete records as well as 657 duplicated entries. Finally, 784 accurate presence points remained. (Fig 1).
A 1 km grid [63,64] was used to reduce any negative effects that may have been caused by spatial autocorrelations in the high-density collections from populated areas [65,66]. With this grid in place, the occurrences were then filtered using SDM Tools [67]; one record was randomly picked from each grid. After filtering, 732 records remained.
General circulation models and environmental variables.
A geographical base map of China was downloaded from DIVA-GIS (http://www.diva-gis.org). The DEM data with a 90 m spatial resolution was obtained from the CGIAR Consortium for Spatial Information (CGIAR-CSI, http://srtm.csi.cgiar.org/)[68]. Using ArcGIS 10.0 (Esri, Redlands, CA, USA), the slope and aspect were calculated.
The soil data were downloaded from the Harmonized World Soil Database (Version 1.2) [69], which was the result of a collaboration between the FAO, Chinese Academy of Sciences (ISSCAS), and several other organizations. Eighteen variables of topsoil were included in this dataset.
Global land cover and vegetation data were acquired from the International Steering Committee for Global Mapping (ISCGM, http://www.iscgm.org/gmd/).
The climate data consisted of a set of 19 bioclimatic variables that were originally derived for BIOCLIM [35], which is currently widely used in various SDMs. The current (1950–2000) climate data were obtained from the WorldClim database [70] (http://worldclim.org/), and the future (2041–2060, 2061–2080) climate data under different scenarios were obtained from Climate Change, Agriculture, and Food Security (http://www.ccafs-climate.org). The most recent climate projections were used in the Fifth Assessment IPCC report. These data were generated from general circulation models (GCMs) for four representative concentration pathways (RCP2.6, RCP4.5, RCP6, and RCP8.5). To reduce the bias in certain areas from one GCM, the multi-model ensemble (MME) average was used. This method has been shown to provide superior results to those obtained from one model [71–74], and it could maintain its essential characteristics while the number of GCMs was reduced from 25 to five [72]. Therefore, we chose five models for use in the Inter-Sectoral Impact Model Inter comparison Project (ISI-MIP) [75], including GFDL-ESM2M, HadGEM2-ES, IPSL-CM5A-LR, MIROC-SEM-CHEM, and NorESM1-M.
A total of 42 environmental variables pertaining to climate, soil, topography, land cover, and vegetation coverage were collected (S1 Table). And all environmental variables were converted into 30 arc second (~1 km) spatial resolution. High correlations among these variables decreased the accuracy of the SDM [76]; therefore, these variables were filtered to reduce this effect. The 19 bioclimatic variables were clustered into four groups that were highly correlated with the mean annual value of temperature and precipitation and their intra-annual precipitation fluctuations [77,78]. According the results in the pretest (S2 Table) and the correlations (|R| > 0.7) between variables, 20 variables were selected from the 42 variables (S2 Table) using the SDM Tools [33,67,79].
Modeling procedure
We predicted the potential distribution of P. tenuifolia under different climatic scenarios using Maxent (http://www.cs.princeton.edu/~schapire/maxent/, version 3.3.3k, [44]). The 20 environmental variables and 732 occurrences were used as inputs for Maxent. Of these records, 25% were randomly selected to be test data, and the remainder were used as training data. The “10th percentile training presence” was chosen as the threshold rule to remove biologically irrelevant noise from the model prediction maps [64,80–82] and resulted in reliable species distributions [83].
All other parameters used the default setting. The model was replicated 10 times.
Using the parameters above, the current potential distribution of P. tenuifolia was projected, and then the future potential distributions from four different RCP scenarios (RCPs 2.6, 4.5, 6.0 and 8.5) were projected in two future periods (2050 and 2070) based on the assumptions that soil, topography, land cover, and vegetation coverage would not change under different climatic scenarios.
The area under the receiver-operating characteristic (ROC) curve (AUC) was employed in the model performance evaluation [44]. The AUC values were between 0 and 1, and higher values indicated better model performance [84]. When the AUC was below 0.5, the model performed was worse than chance. Model performance was classified as poor (0.5–0.6), fair (0.6–0.7), good (0.7–0.8), very good (0.8–0.9), or excellent (0.9–1) [84].
The jackknife results and response curves were used to evaluate the importance of each environmental variable to the distribution of P. tenuifolia. According to the presence probability, suitable regions for P. tenuifolia were divided into three levels: low suitable regions (0-T), medium suitable regions (T-0.5), and highly suitable regions (0.5–1.0). The term “T” represents the 10th percentile threshold [82]. Finally, the suitable regions and their altitudes under different climate conditions were analyzed using ArcGIS 10. Difference in the mean pixel freauency and mean altitude of highly suitable region for current and future climatic scenarios were tested using one-way ANOVAs and Tukey’s HSD comparisons in IBM SPSS Statistics 22 (IBM, Armonk, NY, USA).
Results
Prediction accuracy
Table 1 lists the AUC values of the Maxent predictions for the potential P. tenuifolia distributions based on environmental variables. The AUC values in all scenarios exceeded 0.8.
Dominant environmental variables
According to the results of the jackknife analyses (Fig 2) and the estimations of the relative contributions (Table 2) of environmental variables to the Maxent models under the current scenario, the top three variables that affected the distribution of P. tenuifolia were as follows: precipitation in the warmest quarter (BIO18), annual mean temperature (BIO1), and altitude (ALT). The total contribution of the three variables was 73.91%, suggesting that P. tenuifolia distributions are strongly influenced by these three environmental variables. The total contribution of five climatic variables was 61.82%. The total contribution of three topographic variables was 19.96%. However, the contribution of ten soil variables totaled 12.20%. This showed that climatic variables had a much greater impact than topography, soil, land cover, and vegetation coverage on the habitat distribution of P. tenuifolia.
The abbreviations of variables could be looked up from S1 Table.
The graph shows the result of the Jackknife test of variable importance using regularized training gain, test gain and AUC on test data respectively. The red bars indicate the gain using solo environmental variable, the hollow bars indicate the gain excluding the single variable from the full model, and the blue bars indicate the gain considering all variables.
Relationship between the species distribution and the dominant environmental variables
The relationship between presence probability of P. tenuifolia and precipitation in the warmest quarter (Fig 3a) showed that when the precipitation was below 23 mm, the presence probability of finding P. tenuifolia was less than 5%. When the precipitation in the warmest quarter ranged between 148 mm and 512 mm, the presence probability was higher than 50%. Precipitation in the warmest quarter from 148 mm and 512 mm might be effective predictor for modeling the potential distribution of P. tenuifolia.
(a) Precipitation in the Warmest Quarter, (b) Annual Mean Temperature and (c) Altitude.
The annual mean temperature was another important variable that affected the distribution of P. tenuifolia (Fig 3b). When the annual mean temperature fell below -2.5°C, the presence probability of finding P. tenuifolia was less than 5%. When this temperature was in the range of 8.4°C to 15.4°C, the presence probability was higher than 50%. In other words, regions with such a range of annual temperature may be better candidates for planting P. tenuifolia.
Altitude also played an important role in forecasting the potential distribution of P. tenuifolia, which displayed a presence threshold of 3740 m (Fig 3c). The presence probability gradually declined with increasing elevation from 100 m to 3087 m.
Potential suitable distribution areas for P. tenuifolia
Using ArcGIS 10.0, the potential distribution in the current situation predicted by Maxent was analyzed (Fig 4). It showed that the highly suitable regions (the presence probability was greater than or equal to 50%) for P. tenuifolia are primarily located in the north of China, including the following: West Liaoning, North and West Hebei, most of Beijing, Central and South Shanxi, the northwest corner and west of Henan, North and Central Shaanxi, East and South Gansu, Northwest Hubei, Central Shandong, and some fragmented plots in Heilongjiang, Jilin, Inner Mongolia, Ningxia, Anhui, Zhejiang, Jiangxi, Hunan, Chongqing, Sichuan, and Guizhou. The total area of this region is approximately 0.68 × 106 km2, which accounts for 7.07% of the total land in China.
Blue to red colors show the increase of presence probability. 0 indicates completely unsuitable and 1 indicates optimal.
The medium suitable region was located around the edges of the highly suitable region as follows: Southwest Heilongjiang, most of Jilin, Central and West Liaoning, Northwest and Southeast Hebei, West and East Shandong, East Henan, West Hubei, South Shaanxi, Central Gansu, the east corner of Qinghai, most of Ningxia, Central and West Inner Mongolia, the northeast corner of Yunnan, North Guizhou, Central Hunan, West and Southeast Jiangxi, West Zhejiang, South Anhui, and other scattered regions. The total area of this region was approximately 1.21 × 106 km2, which accounts for 12.62% of the total land in China. The low suitable region had the widest distribution, which included Xinjiang, Qinghai, Xizang, West Gansu, most of Sichuan, most of Yunnan, Guangxi, Guangdong, Hainan, Taiwan, West and North Inner Mongolia, and North Heilongjiang. The total area of this region was approximately 7.71 × 106 km2, which accounts for approximately 80.31% of the total land in China.
Furthermore, the altitudinal pattern of the highly suitable habitat under the current situation was analyzed. Thirty percent of the highly suitable habitat was below 500 m, 70% was below 1200 m, and almost 90% was below 1700 m. The highly suitable habitat above 1700 m consisted of slightly less than 10%. In the medium suitable region, the region below 200 m accounted for nearly 30%, the region below 600 m accounted for more than 50%, and 90% of this region was not higher than 2000 m (Fig 5).
The cells indicate the proportion of each hundred altitude to the whole.
Changes of the suitable climatic conditions according to climate warming scenarios
The potential distributions of P. tenuifolia under the four RCP scenarios (RCP2.6, RCP4.5, RCP6.0 and RCP8.5) in 2050 and 2070 were compared and analyzed in Fig 6.
Blue to red colors indicate the presence probability of the areas from 0 to 1.
Compared to the current distribution, in 2050, the total area of the medium suitable region for P. tenuifolia under the four RCPs (RCP2.6, RCP4.5, RCP6.0 and RCP8.5) will increase by 9.43%, 0.03%, 5.00% and 4.81%, respectively. The total area of the highly suitable district will decrease by 4.09%, 1.00%, 8.37% and 8.55% (Fig 7). Under scenario RCP2.6 2050 and RCP6.0 2050, the areas of the low suitable regions will decline by 1.12% and 0.05%, but almost no change will occur under RCP4.5 2050 and RCP8.5 in 2050. In Northwest Fujian, East Jiangxi, West Zhejiang, Northeast Yunnan, West Guizhou, North Inner Mongolia, East Henan, and East Liaoning, low suitability regions will be converted into medium suitable regions. The region in which medium suitable regions will become low suitable regions is mainly distributed in East Jilin, Hunan, Guizhou, and North Inner Mongolia. The highly suitable region that will be transformed from medium suitable regions includes Southwest Liaoning, North Hebei, Northwest Shanxi, West Henan, Southeast Shaanxi, South Ningxia, and South Gansu. In Central Shandong, Southwest Shanxi, Central Shaanxi, Central Gansu, Central Ningxia, and West Jilin, highly suitable regions will be transformed into medium suitability regions. Altogether, the potential distribution in 2050 will contract in the northeast and expand in southeast, southwest and north (Fig 8).
The mean pixel frequency of highly suitable habitat in the current scenario was set as 1.0 and the others were expressed relative to it. Bars indicate means ± sd (n = 10 circle), and significant differences at P < 0.01 within in each group are indicated by different letters (a, b and c).
In 2070, the total areas of the medium suitable region will increase by 9.59%, 0.04%, 1.05% and 2.23%. The total highly suitable habitat areas will decrease by 8.19%, 7.39%, 6.31% and 8.83%. The low suitable regions will increase by 0.64%, 0.39% and 0.43% under RCP4.5, RCP 6.0 and RCP8.5, respectively, but will decline by 0.79% under RCP2.6. The distribution pattern will be similar to that in 2050 and shows an expansion in the north, southeast, and southwest and a contraction in the east, west, and southwest (Fig 8).
A separate analysis that focused on the altitude patterns among the highly suitable districts (Fig 9) is shown below. Under all scenarios, except for RCP 4.5 2050, the highly suitable habitat area below 100 m will increase in the range of 4.37% to 113.18%, and decrease in the range of -7.34% to -42.79% between 100 m and 200 m. In the range of 200 m to 1500 m, the highly suitable habitat area will fluctuate between -18.41% and 9.63%. The highly suitable habitat area in the range of 1500 m to 2400 m, will vary between -34.29% and 19.92%. Above 2400 m, under all the future scenarios except for RCP2.6 in 2050, highly suitable habitat area will decline from -5.30% to -100% as the increasing elevation.
The blue columns indicate the future highly suitable habitat area increased relative to the current in the corresponding altitudinal interval, and the red decreased. The following numbers were the magnitude of changes.
No significant difference was found between the current mean elevation and that predicted in 2050 and 2070. However, a significant difference was found between the mean elevations of highly suitable districts in 2050 and that in 2070 (p < 0.01); the mean elevation of the highly suitable regions in 2050 was higher than that predicted in 2070 (Fig 10).
Bars indicate means ± sd (n = 10 circle), and significant differences at P < 0.01 within each group are indicated by different letters (a, b and c).
Discussion
The rapid increase in the global mean temperature threatens the growth and survival of many wild species[14]. The habitats of many plants have been altered [13], and P. tenuifolia is no exception. The potential distributions of several traditional Chinese herbs are predicted to change under global warming [29,46,57], and P. tenuifolia would be impacted as well. Therefore, we modeled the potential distribution of this medicinal plant and predicted how its distribution might be affected under current and eight different future climatic scenarios.
Prediction Performance of Maxent
In our research, the AUC was adopted to evaluate the discrimination performance of Maxent. The AUC is one of the best model evaluation indexes [48,85]. All of the AUC values for the model predicting the potential distribution of P. tenuifolia under different climate conditions were above 0.8, which is considered to be very good for Maxent projections of geographic distribution [44,84].
Dominant environmental variables
The top three variables with high contributions toward determining the distribution of P. tenuifolia were precipitation in the warmest quarter, annual mean temperature, and altitude. The precipitation in the warmest quarter showed the highest contribution (28.71%) and ranged from 148 mm to 512 mm in highly suitable areas. The annual mean temperature provided an almost identical contribution. The highly suitable range was 8.4°C to 15.4°C. Altitude was a significant variable that could explain the distribution of P. tenuifolia. Its contribution accounted for 17.13%, and the highly suitable range was 100 m to 2000 m. These results were consistent with the current habit requirements of P. tenuifolia, which grows well in cool temperatures, is drought resistant, and is commonly found on sunny slopes, forest edges, roadsides, and ridges of fields [7].
Current suitable habitat
The Maxent predictions showed that the highly suitable habitat of P. tenuifolia was mainly located in Heilongjiang, Jilin, Liaoning, Inner Mongolia, Hebei, Henan, Beijing, Shanxi, Shaanxi, Ningxia, Gansu, Qinghai, Sichuan, Hubei, Hunan, Shandong, Anhui, Zhejiang and Jiangxi; these regions were consistent with the description in Flora of China [7]. Our research also found 90% of the highly suitable habitat in the current scenario was below 1700 m, and less than 10% of medium suitable habitat located above 2000 m.
Impacts of climate change on the potential distribution
Previous studies have concluded that global warming will greatly influence species distributions by causing expansions, shifts, or contractions in the species ranges [14,29,86,87]. Our research has revealed that the area of highly suitable habitat for P. tenuifolia in all eight future scenarios will decrease, and the medium suitable habitat will increase relative to the current area. In four scenarios (RCP2.6 2050, RCP6.0 2050, RCP8.5 2050 and RCP2.6 2070), the low and highly suitable habitat will decline, and the medium suitable habitat will increase. In the other four scenarios (RCP4.5 2050, RCP4.5 2070, RCP6.0 2070 and RCP8.5 2070), the low and medium suitable habitats will increase, and the highly suitable habitat will decline.
Not only the range but also the altitudinal limits of species distribution were affected by global warming. The growth line for plants and butterflies has declined in the Alpine mountain regions[88]. Shrestha found that the Chinese Caterpillar Fungus (Ophiocordyceps sinensis) expanded its range at both high and low altitudes [89]. However, in our study, the current mean altitude of highly suitable habitats was lower than that in scenarios RCP2.6 and RCP 4.5 in 2050, almost equal to the value in scenarios RCP6.0 and RCP 8.5 in 2050, but was higher than that in all scenarios in 2070. The centroids of the highly suitable habitat for P. tenuifolia in all four scenarios in 2050 and RCP2.6 in 2070 moved to lower latitudes, but they moved to higher latitudes in scenarios RCP45, RCP6.0 and RCP8.5 in 2070. Analyzing of the variation in every hundred meters elevation relative to the current suggested that the highly suitable habitat of P. tenuifolia was easier to be affected by climatic change at high elevation.
The result under different scenarios did not show same tendency, which may be due to several reasons. First, the forecasted ecological niche was wider than the actual niche, which resulted in a larger projected distribution. Second, other environmental variables may affect the distribution range in addition to the selected variables. Third, the MME has been shown to be better than any one model alone [73], but the result has not been validated by observations, and some deviations may occur.
Conclusions
This study projected the potential distributions of P. tenuifolia under current and future climate change scenarios and determined the dominant environmental variables that affect changes in the distribution. This result could play an important role in location selection for P. tenuifolia cultivation and wild P. tenuifolia reserve design. In addition, the changes in the range and altitude of the highly suitable habitats for P. tenuifolia were analyzed and compared. The results indicate that the highly suitable habitat of P. tenuifolia will obviously decrease under future climate change scenarios in 2050 and 2070. We also found that the area of highly suitable habitats below 500 m and above 2300 m will decline in the future. Therefore, the impact of climate change on plant resource protection and sustainable development must be thoroughly investigated, particularly for estimating the extinction risks of environmentally sensitive species.
Supporting Information
S1 Table. Environmental variables used for modeling the potential distribution of P. tenuifolia.
The bold represented the selected variables used to develop the models, and the others were not used to develop the model.
https://doi.org/10.1371/journal.pone.0163718.s001
(DOCX)
S2 Table. Percent contribution and permutation importance of each environmental variable to model performance under current climatic scenarios in pretest.
https://doi.org/10.1371/journal.pone.0163718.s002
(DOCX)
Acknowledgments
We sincerely thank Prof. Ran Haichao of Hebei University of Science & Technology in China for his valuable comments with regard to data processing. Furthermore, we thank Min Li, Shuo Shi, Pingping Wan, and all others who offered their assistance during our work.
Author Contributions
- Conceptualization: HJ TL JZ.
- Formal analysis: HJ TL.
- Funding acquisition: HJ JZ.
- Investigation: HJ LL YZ.
- Methodology: HJ TL JZ.
- Project administration: HJ JZ.
- Resources: LP.
- Supervision: JZ.
- Visualization: HJ.
- Writing – original draft: HJ TL.
- Writing – review & editing: HJ TL.
References
- 1.
Tang W, Eisenbrand G. Chinese drugs of plant origin. Berling, Heidelberg, New York, London, Paris, Tokyo, Hong Kong, Barcelona, Budapest: Springer-Verlag; 1992.
- 2.
State Pharmacopooeia Commission of PRC. Pharmacopoeia of the People's Republic of China. 1st ed. Beijing: China Medical Science Press; 2005.
- 3. Zhang F, Li X, Li Z, Xu X, Peng B, Qin X, et al. UPLC/Q-TOF MS-based metabolomics and qRT-PCR in enzyme gene screening with key role in triterpenoid saponin biosynthesis of Polygala tenuifolia. PLoS ONE. 2014;9: e105765. pmid:25148032
- 4. Zhang F, Song X, Li L, Wang J, Lin L, Li C, et al. Polygala tenuifolia polysaccharide PTP induced apoptosis in ovarian cancer cells via a mitochondrial pathway. Tumor Biol. 2015;36: 2913–2919.
- 5. Dong XZ, Huang CL, Yu BY, Hu Y, Mu LH, Liu P. Effect of Tenuifoliside A isolated from Polygala tenuifolia on the ERK and PI3K pathways in C6 glioma cells. Phytomedicine. 2014;21: 1178–1188. pmid:24877714
- 6. Tian H, Xu X, Zhang F, Wang Y, Guo S, Qin X, et al. Analysis of Polygala tenuifolia transcriptome and description of secondary metabolite biosynthetic pathways by illumina sequencing. International Journal of Genomics. 2015;2015: 782635. pmid:26543847
- 7.
Flora of China Editorial Committee. Flora of China (checklist & addendum). Beijing & St. Louis: Science Press & Missouri Botanical Garden Press; 1988. doi: https://doi.org/10.3897/phytokeys.59.6261 26929706
- 8. Fang MF, Li J, Zhou TH, Yang J, Zhao GF. Genetic diversity in natural populations of the medicinal herb Polygala tenuifolia Willd. and its implications for conservation. Biochemical Systematics & Ecology. 2012;44: 400–406.
- 9.
Endangered Species Scientific Commission PRC. List of key state-protected wild medicinal species. 2012. Available: http://www.cites.org.cn/article/show.php?itemid=589.
- 10.
IPCC. Climate change 2014: synthesis report. Contribution of working groups I, II and III to the fifth assessment report of the intergovernmental panel on climate change. Geneva, Switzerland: IPCC; 2014.
- 11. Camille P, Gary Y. A globally coherent fingerprint of climate change impacts across natural systems. Nature. 2003;421: 37–42. pmid:12511946
- 12. Root TL, Price JT, Hall KR, Schneider SH, Rosenzweig C, Pounds JA. Fingerprints of global warming on wild animals and plants. Nature. 2003;421: 57–60. pmid:12511952
- 13.
Brummitt N, Bachman S. Plants under pressure a global assessment. The first report of the IUCN sampled red list index for plants. London, UK: Natural History Museum Press; 2010.
- 14. Thomas CD, Cameron A, Green RE, Bakkenes M, Beaumont LJ, Collingham YC, et al. Extinction risk from climate change. Nature. 2004;427: 145–148. pmid:14712274
- 15. Chen IC, Hill JK, Ohlemuller R, Roy DB, Thomas CD. Rapid range shifts of species associated with high levels of climate warming. Science. 2011;333: 1024–1026. pmid:21852500
- 16. Bellard C, Bertelsmeier C, Leadley P, Thuiller W, Courchamp F. Impacts of climate change on the future of biodiversity. Ecology Letters. 2012;15: 365–377. pmid:22257223
- 17. Summers DM, Bryan BA, Crossman ND, Meyer WS. Species vulnerability to climate change: impacts on spatial conservation priorities and species representation. Global Change Biology. 2012;18: 2335–2348.
- 18. Pressey RL, Mar C, Watts ME, Cowling RM, Wilson KA. Conservation planning in a changing world. Trends in Ecology & Evolution. 2007;22: 583–592.
- 19. Booth TH, Nix HA, Busby JR, Hutchinson MF. BIOCLIM: the first species distribution modelling package, its early applications and relevance to most current MAXENT studies. Diversity and Distributions. 2014;20: 1–9.
- 20.
Franklin J, Miller JA. Mapping species distributions: spatial inference and prediction. New York: Cambridge University Press; 2009.
- 21. Miller J. Species distribution modeling. Geography Compass. 2010;4: 490–509.
- 22. Acevedo P, Real R. Favourability: concept, distinctive characteristics and potential usefulness. Naturwissenschaften. 2012;99: 515–522. pmid:22660474
- 23. Franklin J. Species distribution models in conservation biogeography: developments and challenges. Diversity and Distributions. 2013;19: 1217–1223.
- 24. Guisan A, Thuiller W. Predicting species distribution: offering more than simple habitat models. Ecology letters. 2005;8: 993–1009.
- 25.
Elith J, Franklin J. Species distribution modeling. In: Levin SA, editor. Encyclopedia of Biodiversity Second Edition ed. Oxford: Academic Press; 2013. pp. 692–705.
- 26. Ihlow F, Dambach J, Engler JO, Flecks M, Hartmann T, Nekum S, et al. On the brink of extinction? How climate change may affect global chelonian species richness and distribution. Global Change Biology. 2012;18: 1520–1530.
- 27. Guisan A, Tingley R, Baumgartner JB, Naujokaitis-Lewis I, Sutcliffe PR, Tulloch AIT, et al. Predicting species distributions for conservation decisions. Ecology Letters. 2013;16: 1424–1435. pmid:24134332
- 28. Nakao K, Higa M, Tsuyama I, Matsui T, Horikawa M, Tanaka N. Spatial conservation planning under climate change: Using species distribution modeling to assess priority for adaptive management of Fagus crenata in Japan. Journal for Nature Conservation. 2013;21: 406–413.
- 29. Yuan H, Wei Y, Wang X. Maxent modeling for predicting the potential distribution of Sanghuang, an important group of medicinal fungi in China. Fungal Ecology. 2015;17: 140–145.
- 30. Angelieri CC, Adams-Hosking C, Ferraz KM, de Souza MP, McAlpine CA. Using species distribution models to predict potential landscape restoration effects on puma conservation. PLoS ONE. 2016;11: e0145232. pmid:26735128
- 31. Murray JV, Stokes KE, van Klinken RD. Predicting the potential distribution of a riparian invasive plant: the effects of changing climate, flood regimes and land-use patterns. Global Change Biology. 2012;18: 1738–1753.
- 32. Sheppard CS, Burns BR, Stanley MC. Predicting plant invasions under climate change: are species distribution models validated by field trials? Global Change Biology. 2014;20: 2800–2814. pmid:24446429
- 33. Murray KA, Retallick RWR, Puschendorf R, Skerratt LF, Dan R, McCallum HI, et al. Assessing spatial patterns of disease risk to biodiversity: implications for the management of the amphibian pathogen, Batrachochytrium dendrobatidis. Journal of Applied Ecology. 2011;48: 163–173.
- 34. Keith DA, Elith J, Simpson CC. Predicting distribution changes of a mire ecosystem under future climates. Diversity and Distributions. 2014;20: 440–454.
- 35. Busby JR. BIOCLIM-a bioclimate analysis and prediction system. Plant Prot Q. 1991
- 36. Hirzel AH, Hausser J, Chessel D, Perrin N. Ecological-niche factor analysis: How to compute habitat-suitability maps without absence data? Ecology. 2002;83: 2027–2036.
- 37. Friedman J, Hastie T, Tibshirani R. Additive logistic regression: a statistical view of boosting. Annals of Statistics. 2000;28: 374–376.
- 38.
Sutherst RW. Predicting insect distributions in a changed climate. In: Harrington R, Stork NE, editors. Insects in a changing environment. London, UK: Academic Press; 1995. pp. 59–91.
- 39. Carpenter G, Gillison AN, Winter J. DOMAIN: a flexible modelling procedure for mapping potential distributions of plants and animals. Biodiversity & Conservation. 1993;2: 667–680.
- 40. Yee TW, Mitchell ND. Generalized additive models in plant ecology. Journal of Vegetation Science. 1991;2: 587–602.
- 41. Lehmann A, Overton JM, Leathwick JR. GRASP: generalized regression analysis and spatial prediction. Ecological Modelling. 2003;160: 165–183.
- 42. Stockwell D, Peters D. The GARP modelling system: problems and solutions to automated spatial prediction. International Journal of Geographical Information Science. 1999;13: 143–158.
- 43. Walker PA, Cocks KD. HABITAT: a proceduce for modelling a disjoint environmental envelope for a plant or animal species. Global Ecology & Biogeography Letters. 1991;1: 108–118.
- 44. Phillips SJ, Anderson RP, Schapire RE. Maximum entropy modeling of species geographic distributions. Ecological modelling. 2006;190: 231–259.
- 45. Thuiller W, Lafourcade B, Engler R, Araújo MB. BIOMOD–a platform for ensemble forecasting of species distributions. Ecography. 2009;32: 369–373.
- 46. Remya K, Ramachandran A, Jayakumar S. Predicting the current and future suitable habitat distribution of Myristica dactyloides Gaertn. using MaxEnt model in the Eastern Ghats, India. Ecological Engineering. 2015;82: 184–188.
- 47.
Peterson AT, Soberón J, Pearson RG, Anderson RP, Martínez-Meyer E, Nakamura M, et al. Ecological niches and geographic distributions (MPB-49). Princeton: Princeton University Press; 2011.
- 48. Elith J, Graham CH, Anderson RP, Dudík M, Ferrier S, Guisan A, et al. Novel methods improve prediction of species' distributions from occurrence data. Ecography. 2006;29: 129–151.
- 49. Merow C, Smith MJ, Silander JA. A practical guide to MaxEnt for modeling species' distributions: what it does, and why inputs and settings matter. Ecography. 2013;36: 1058–1069.
- 50. Barnes MA, Jerde CL, Wittmann ME, Chadderton WL, Ding JQ, Zhang JL, et al. Geographic selection bias of occurrence data influences transferability of invasive Hydrilla verticillata distribution models. Ecology and Evolution. 2014;4: 2584–2593. pmid:25360288
- 51. Pinkernell S, Beszteri B. Potential effects of climate change on the distribution range of the main silicate sinker of the Southern Ocean. Ecology and Evolution. 2014;4: 3147–3161. pmid:25473469
- 52. Barve N, Martin C, Brunsell NA, Peterson AT. The role of physiological optima in shaping the geographic distribution of Spanish moss. Global Ecology and Biogeography. 2014;23: 633–645.
- 53. Yu J, Ma YH, Guo SL. Modeling the geographic distribution of the epiphytic moss Macromitrium japonicum in China. Annales Botanici Fennici. 2013;50: 35–42.
- 54. Barbet-Massin M, Jetz W. The effect of range changes on the functional turnover, structure and diversity of bird assemblages under future climate scenarios. Global Change Biology. 2015;21: 2917–2928. pmid:25931153
- 55. Mendoza-Gonzalez G, Martinez ML, Rojas-Soto OR, Vazquez G, Gallego-Fernandez JB. Ecological niche modeling of coastal dune plants and future potential distribution in response to climate change and sea level rise. Global Change Biology. 2013;19: 2524–2535. pmid:23625760
- 56. Yang XQ, Kushwaha SPS, Saran S, Xu J, Roy PS. Maxent modeling for predicting the potential distribution of medicinal plant, Justicia adhatoda L. in Lesser Himalayan foothills. Ecological Engineering. 2013;51: 83–87.
- 57. Yi Y, Cheng X, Yang Z, Zhang S. Maxent modeling for predicting the potential distribution of endangered medicinal plant (H. riparia Lour) in Yunnan, China. Ecological Engineering. 2016;92: 260–269.
- 58. Zhang L, Cao B, Bai C, Li G, Mao M. Predicting suitable cultivation regions of medicinal plants with Maxent modeling and fuzzy logics: a case study of Scutellaria baicalensis in China. Environmental Earth Sciences. 2016;75: 361.
- 59. Pineda E, Lobo JM. Assessing the accuracy of species distribution models to predict amphibian species richness patterns. Journal of Animal Ecology. 2009;78: 182–190. pmid:18771504
- 60. Elith J, Graham CH. Do they? How do they? WHY do they differ? On finding reasons for differing performances of species distribution models. Ecography. 2009;32: 66–77.
- 61. Hijmans RJ, Graham CH. The ability of climate envelope models to predict the effect of climate change on species distributions. Global Change Biology. 2006;12: 2272–2281.
- 62. Synes NW, Osborne PE. Choice of predictor variables as a source of uncertainty in continental-scale species distribution modelling under climate change. Global Ecology and Biogeography. 2011;20: 904–914.
- 63. Anderson RP, Raza A. The effect of the extent of the study region on GIS models of species geographic distributions and estimates of niche evolution: preliminary tests with montane rodents (genus Nephelomys) in Venezuela. Journal of Biogeography. 2010;37: 1378–1393.
- 64. Radosavljevic A, Anderson RP. Making better MAXENT models of species distributions: complexity, overfitting and evaluation. Journal of Biogeography. 2014;41: 629–643.
- 65. Adams-Hosking C, Grantham HS, Rhodes JR, McAlpine C, Moss PT. Modelling climate-change-induced shifts in the distribution of the koala. Wildlife Research. 2011;38: 122–130.
- 66. Mckenney DW, Pedlar JH, Lawrence K, Campbell K, Hutchinson MF. Beyond traditional hardiness zones: using climate envelopes to map plant range limits. Bioscience. 2007;57: 929–937.
- 67. Brown JL. SDMtoolbox: a python-based GIS toolkit for landscape genetic, biogeographic and species distribution model analyses. Methods in Ecology & Evolution. 2014;5: 694–700.
- 68.
Jarvis A, Reuter HI, Nelson A, Guevara E. Hole-filled SRTM for the globe Version. 2008. Available: the CGIAR-CSI SRTM 90m Database (http://srtm.csi.cgiar.org).
- 69.
Nachtergaele. F, Velthuizen. Hv, Verelst. L, Wiberg. D, Batjes. N, Dijkshoorn. K, et al. Harmonized world soil database (Version 1.2). Laxenburg, Austria: Food and Agriculture Organization of the UN, International Institute for Applied Systems Analysis, ISRIC—World Soil Information, Institute of Soil Science—Chinese Academy of Sciences, Joint Research Centre of the EC; 2012.
- 70. Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A. Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology. 2005;25: 1965–1978.
- 71. Moss RH, Edmonds JA, Hibbard KA, Manning MR, Rose SK, Vuuren DPV, et al. The next generation of scenarios for climate change research and assessment. Nature. 2010;463: 747–756. pmid:20148028
- 72. Mendlik T, Gobiet A. Selecting climate simulations for impact studies based on multivariate patterns of climate change. Climatic Change. 2016;135: 381–393. pmid:27429499
- 73. Pierce DW, Barnett TP, Santer BD, Gleckler PJ. Selecting global climate models for regional climate change studies. Proceedings of the National Academy of Sciences of the United States of America. 2009;106: 8441–8446. pmid:19439652
- 74. Wang L, Chen W. A CMIP5 multimodel projection of future temperature, precipitation, and climatological drought in China. International Journal of Climatology. 2014;34: 2059–2078.
- 75. Warszawski L, Frieler K, Huber V, Piontek F, Serdeczny O, Schewe J. The inter-sectoral impact model intercomparison project (ISI-MIP): project framework. Proceedings of the National Academy of Sciences. 2014;111: 3228–3232.
- 76. Elith J, Kearney M, Phillips S. The art of modelling range-shifting species. Methods in Ecology & Evolution. 2010;1: 330–342.
- 77. van Zonneveld M, Castaneda N, Scheldeman X, van Etten J, Van Damme P. Application of consensus theory to formalize expert evaluations of plant species distribution models. Applied Vegetation Science. 2014;17: 528–542.
- 78. Zonneveld MV, Jarvis A, Dvorak W, Lema G, Leibing C. Climate change impact predictions on Pinus patula and Pinus tecunumanii populations in Mexico and Central America. Forest Ecology & Management. 2009;257: 1566–1576.
- 79. Pickles RSA, Thornton D, Feldman R, Marques A, Murray DL. Predicting shifts in parasite distribution with climate change: a multitrophic level approach. Global Change Biology. 2013;19: 2645–2654. pmid:23666800
- 80. Liu C, Berry PM, Dawson TP, Pearson RG. Selecting thresholds of occurrence in the prediction of species distributions. Ecography. 2005;28: 385–393.
- 81. Bean WT, Stafford R, Brashares JS. The effects of small sample size and sample bias on threshold selection and accuracy assessment of species distribution models. Ecography. 2012;35: 250–258.
- 82. Engler JO, Rodder D, Elle O, Hochkirch A, Secondi J. Species distribution models contribute to determine the effect of climate and interspecific interactions in moving hybrid zones. Journal of Evolutionary Biology. 2013;26: 2487–2496. pmid:24016292
- 83. Escalante T, Rodríguez-Tapia G, Linaje M, Illoldi-Rangel P, González-López R. Identification of areas of endemism from species distribution models: threshold selection and nearctic mammals. TIP Revista Especializada en Ciencias Químico-Biológicas. 2013;16: 5–17.
- 84. Swets JA. Measuring the accuracy of diagnostic systems. Science. 1988;240: 1285–1293. pmid:3287615
- 85. Mcpherson JM, Walter J, Rogers DJ. The effects of species’ range sizes on the accuracy of distribution models: ecological phenomenon or statistical artefact? Journal of Applied Ecology. 2004;41: 811–823.
- 86. Li R, Xu M, Wong MHG, Qiu S, Sheng Q, Li X, et al. Climate change-induced decline in bamboo habitats and species diversity: implications for giant panda conservation. Diversity and Distributions. 2015;21: 379–391.
- 87. Lenoir J, Gégout JC, Marquet PA, De Ruffray P, Brisse H. A significant upward shift in plant species optimum elevation during the 20th century. science. 2008;320: 1768–1771. pmid:18583610
- 88. Tobias R, Matthias P, Valentin A. Plants, birds and butterflies: short-term responses of species communities to climate warming vary by taxon and with altitude. PLoS ONE. 2014;9: e82490. pmid:24416144
- 89. Shrestha UB, Bawa KS. Impact of climate change on potential distribution of Chinese caterpillar fungus (Ophiocordyceps sinensis) in Nepal Himalaya. PLoS ONE. 2014;9: e106405. pmid:25180515