Weixue Luo1,2‡
Shunxin Han1‡
Ting Yu1
Peng Wang1
Yuxuan Ma1
Maji Wan1
Jinchun Liu1,2
Zongfeng Li1,2
Jianping Tao1,2*†- 1Key Laboratory of Eco-environments in Three Gorges Reservoir Region (Ministry of Education), Chongqing Key Laboratory of Plant Ecology and Resources Research in Three Gorges Reservoir Region, School of Life Sciences, Southwest University, Chongqing, China
- 2Chongqing Jinfo Mountain Karst Ecosystem National Observation and Research Station, Southwest University, Chongqing, China
Climate change exerts profound influences on the ecological environments on a global scale, leading to habitat destruction and altering distribution patterns for numerous plant species. Traditional Chinese medicinal plants, such as those belonging to the Sambucus genus, have been extensively utilized for several centuries to treat fractures, rheumatism, and inflammation. However, our understanding of their geographic distribution and climatic adaptation within China still needs to be improved. In this study, we screened the optimal predictive model (random forest model) to predict the potential suitable distribution of three Sambucus species (Sambucus adnata, Sambucus javanica, and Sambucus williamsii) across China under both current and future climate scenarios. Moreover, we identified key climate factors that influence their potential distributions. Our findings revealed that S. adnata and S. javanica are predominantly shaped by temperature seasonality and mean diurnal range, respectively, whereas S. williamsii is significantly affected by the precipitation of the wettest month. Currently, S. williamsii is primarily distributed in north and central south China (covering 9.57 × 105 km2), S. javanica is prevalent in the south and east regions (covering 6.41×105 km2), and S. adnata predominantly thrives in the southwest China (covering 1.99×105 km2). Under future climate change scenarios, it is anticipated that S. adnata may migrate to higher latitudes while S. javanica may shift to lower latitudes. However, potentially suitable areas for S. williamsii may contract under certain scenarios for the years 2050 and 2090, with an expansion trend under the SSP585 scenario for the year 2090. Our study emphasizes the importance of climatic variables in influencing the potential geographic distribution of Sambucus species. These findings provide valuable theoretical insights for the preservation, cultivation, and utilization of Sambucus medicinal plant resources in the context of ongoing climate change.
1 Introduction
Sambucus genus plants hold significant medicinal value in China and have been extensively utilized for centuries owing to their multifaceted pharmacological properties (Wieland et al., 2021; Waswa et al., 2022). The initial classification of the Sambucus genus under the Caprifoliaceae family was later revised to Adoxaceae, incorporating genetic evidence and morphological comparisons, including nucleotide sequences and preliminary morphology. The genus Sambucus comprises approximately 20 species distributed in tropical and subtropical regions (Li et al., 2021). In China, the primary Sambucus species include Sambucus adnata, Sambucus javanica, and Sambucus williamsii (Liao et al., 2005; Xiao et al., 2016; Yuan et al., 2020). S. adnata, a perennial tall herb or subshrub, is commonly employed to treat fractures, traumatic injuries, and chronic nephritis diseases (Yuan et al., 2020; Li et al., 2021; Waswa et al., 2022). S. javanica, a perennial herb with multiple pharmacological functions, including anti-inflammatory and liver protection, is generally used to treat rheumatism, chronic inflammation of the airways, and viral hepatitis (Liao et al., 2005; Zhang et al., 2010; Chen et al., 2020). In addition, S. williamsii, another well-known herb, exhibits notable anti-inflammatory and fracture-healing properties and is frequently used to treat fractures and joint illness (Lv et al., 2015; Xiao et al., 2016). In recent years, herbal remedies for treating respiratory ailments have gained increased attention, particularly in the context of global respiratory infectious diseases such as COVID-19 (Wieland et al., 2021). Nevertheless, despite the recognized pharmacological worth of Sambucus plants (Waswa et al., 2022), there remains a lack of clarity regarding their geographical distribution and conservation status.
Climate change has been identified as the primary threat to biodiversity in the next century, with an estimated one-quarter of plant species at risk of extinction (Outlook, 2010; Dawson et al., 2011; Roman-Palacios and Wiens, 2020). According to the IPCC fifth assessment report, the global average temperature is projected to increase by 0.3 to 4.8°C by the end of this century (Change, 2013). Under such rapidly changing climate conditions, numerous terrestrial plants will modify their physiological, genetic, and community structure characteristics to shift their spatial distribution towards higher elevations and latitudes (Chen et al., 2011; Shen and Wang, 2013; Becklin et al., 2016; Yi et al., 2018). While Sambucus plants offer substantial economic and pharmacological value, their wild populations are limited, making it challenging to meet future demand, particularly in light of rapid climate change (Wang et al., 2021). Therefore, investigating the suitable geographical distribution of Sambucus species under climate change can serve the dual purpose of enhancing their commercial value and health benefits while preserving their natural habitat.
Global warming primarily results from the increased anthropogenic emissions of atmospheric greenhouse gases (Kim and Bae, 2021). To better understand the correlation between different socioeconomic development patterns and climate change risks, the Coupled Model Intercomparison Project Phase 6 (CMIP6) proposed four shared socioeconomic pathways (SSPs) based on radiative forcing levels and current national circumstances (O'Neill et al., 2014). Numerous studies have demonstrated that climate factors play a crucial role in shaping species distribution and can be used to model habitat changes in response to future climate scenarios (Kaky and Gilbert, 2017; Shrestha et al., 2018). These findings have significant implications for the conservation of medicinal plant resources, including Coptis (Li et al., 2020), Artemisia annua L. (Wang et al., 2022), Homonoia riparia (Yi et al., 2016), and Scutellaria baicalensis (Xu et al., 2020). In addition to climate factors, soil nutrients, and topographic conditions have long been recognized as key drivers of changes in species distribution that limit resources access to wild plant populations (Chytrý et al., 2003; Soberon and Peterson, 2005). Recent studies have shown that topographic conditions control the redistribution of solar radiation and precipitation (Cantón et al., 2004), while soil nutrients indicate the availability of belowground resources (John et al., 2007). Unfortunately, the critical limiting factors determining the potential distribution of Sambucus plants remain unclear.
In recent years, the rapid development of species distribution models (SDMs), based on the ecological and distributional characteristics of species, has provided scientific guidance for prioritizing species conservation (Zhang et al., 2018; Li et al., 2020; Rousseau and Betts, 2022; Zhao et al., 2023). Using the occurrence of target species and environment variables, SDMs can simulate the potential distribution of suitable habitats for species under climate change (Yi et al., 2016; Wang et al., 2022). However, it is important to note that the accuracy of SDMs relies on selecting appropriate species distribution models, which is critical for formulating effective conservation strategies. Over the past two decades, several noteworthy species distribution models have been formulated, including the generalized linear model (GLM, Lehmann et al., 2002), maximum entropy model ((MaxEnt, Phillips et al., 2006), random forest model (RF, Evans et al., 2011), genetic algorithms for rule set generation (GRAP, Stockwell and Peters, 1999), and ecological niche factor analysis (ENFA, Hirzel et al., 2002). The GLM model is frequently employed by scholars due to its simplicity and convenience, effectively handling the relatively simple relationship between response variables and environmental variables (Maria and Udo, 2017; Gábor et al., 2022). Similarly, the MaxEnt model, based on maximum entropy theory, requires only a small sample size to generate reliable results and has become one of the most popular species distribution models (Hernandez et al., 2006; Yi et al., 2016). Meanwhile, the RF model is known for being one of the best performing classification regression tree models, which can capture the complex nonlinear relationships among multiple ecological variables (Mi et al., 2017; Mohapatra et al., 2019; Rather et al., 2020). Although these species distribution models have been well evaluated in most studies, identifying the optimal model for accurately predicting the potential distribution of different plant species remains controversial.
This study focuses on investigating the implications of climate change on the potential distribution of three Sambucus species (S. adnata, S. javanica, and S. williamsii) in China. To achieve this, our study has three primary objectives: (1) to conduct a comparative analysis of the performance of three models (GLM, MaxEnt, and RF); (2) to simulate potential suitable distribution of the three species under current and future climate scenarios (SSP 2-4.5 and SSP 5-8.5) at two time points (2050 and 2090) in China; (3) to evaluate the relative importance of bioclimate, soil, and topographic factors in the spatial distributions of three Sambucus species. We hope this study can provide a basis for the scientific conservation and effective utilization of these important medicinal plants in China.
2 Materials and methods
2.1 Species occurrence data
We collected the specimen records of S. adnata, S. javanica, and S. williamsii dating back to the year 1955 through online databases, including the Global Biodiversity Information Facility (GBIF, https://doi.org/10.15468/dl.xyt3gh) and Chinese Virtual Herbarium (CVH, https://www.cvh.ac.cn/). To ensure the accuracy of geographical locations, we excluded problematic records with duplicate information (i.e., same sampling time and location), ambiguous geo-location and time, and artificially planted specimen information. For records lacking precise geographic coordinates, we used the Baidu Map Picking System to obtain the corresponding longitude and latitude. To minimize the spatial autocorrelation of species distribution points, we kept only one record closest to the centroid within each 5 km×5 km grid cell, using the “spThin” package in R (Aiello-Lammens et al., 2015). In the end, we obtained a total of 239 records for S. adnata, 863 records for S. javanica, and 509 records for S. williamsii in China (Figure 1A). Moreover, considering the spatial heterogeneity of administrative units and climatic conditions, we divided all provinces in China into six subregions: northeast, northwest, north, east, southwest, and central south (Figure 1B).
Figure 1 Original occurrence records of three Sambucus species in China. Occurrence records of S. adnata, S. javanica and S. williamsii in China (A). Red points indicate S. adnata, blue points indicate S. javanica and golden points indicate S. williamsii. (B) The subregions division of China, including Northeast China (NE), Northwest China (NW), North China (N), East China (E), Southwest China (SW), and Central south China (CS).
2.2 Selection and preprocessing of environmental variables
In this study, we initially selected 28 environmental factors (as detailed in Table S1) that may affect the spatial distribution of S. adnata, S. javanica, and S. williamsii. We categorized these factors into three distinct categories: bioclimatic, topographic, and soil factors. The bioclimatic dataset consisted of nineteen variables obtained from the WorldClim database, each with a spatial resolution of 5 km × 5 km (Hijmans et al., 2005). These variables encompassed essential climate parameters such as annual mean precipitation, annual mean temperature, minimum monthly precipitation, and maximum monthly precipitation, spanning the period from 1970 to 2000. This dataset effectively captured both annual and seasonal variations in temperature and precipitation, providing a comprehensive depiction of the prevailing climatic conditions. To project the gradual shifts in the potential distribution of Sambucus species in response to climate change, we utilized both current climate data and future climate projections for the 2050s and 2090s. The 2050s represent the mean values for the period 2041-2060, while the 2090s correspond to the average conditions projected for the years 2081-2100. This implies that the selected variables could effectively reflect the climatic conditions at the specified occurrence coordinates. Shared socioeconomic pathways are often used to model and represent the likely economic patterns of future societies and the resulting emissions of greenhouse gases and air pollutants. To evaluate the potential impact of climate change on these three Sambucus species, we incorporated future climate data estimated by eight global climate models: ACCESS-CM2, BCC-CSM2-MR, CMCC-ESM2, GISS-E2-1-G, IPSL-CM6A-LR, MIROC6, MPI-ESM1-2-HR and MRI-ESM2-0. Specifically, we integrated two shared socioeconomic pathways climate scenarios, SSP245 (moderate scenario) and SSP585 (pessimistic scenario), for the years 2050 and 2090 to accurately predict the potential distribution of these three species. For each future climate scenario, we calculated the mean values of 19 climate variables across the eight climate models.
In addition to bioclimatic factors, we recognized the importance of topographic and soil factors as key drivers of species distribution, as underscored in previous studies (Svenning et al., 2009; Buri et al., 2017). Specifically, we incorporated the latest digital elevation model (DEM) from the Shuttle Radar Topography Mission (SRTM) (Jarvis et al., 2008). Elevation, aspect, and slope were extracted from the SRTM with a spatial resolution of 90 m and subsequently resampled to a coarser resolution of 5 km to maintain consistency with other data layers. Furthermore, we considered six soil parameters crucial for plant growth, including soil pH and total nitrogen (Table S1). These soil variables were sourced from the SoilGrids database (Hengl et al., 2017) at a sampling depth from 5 to 15 cm with a spatial resolution of 1 km and were resampled to a spatial resolution of 5 km.
The selection of environmental variables has a significant impact on the accurate evaluation of species’ habitat suitability. High collinearity between environmental variables can reduce the accuracy of the species distribution models (SDMs) (Cruz-Cárdenas et al., 2014). Firstly, we performed Spearman correlation analysis using the ‘ggcorr’ function in the corrplot R package to exclude predictors with correlation coefficients greater than 0.8 (Pearson et al., 2007). Furthermore, we automatically selected a subset of predictors for each species by iteratively excluding variables with the variance inflation factor (VIF) greater than 10 using the ‘vifstep’ function of the usdm R package (Naimi et al., 2014). In the end, we retained 10 environmental variables for modeling S. adnata (Figure S1A), 12 for modeling S. javanica (Figure S1B), and 11 for modeling S. williamsii (Figure S1C).
2.3 Model implementation and accuracy assessment
To identify the most accurate species distribution models for S. adnata, S. javanica, and S. williamsii in China, we evaluated three widely used SDMs: generalized linear models (GLM; Nelder and Wedderburn, 1972), maximum entropy models (MaxEnt; Phillips et al., 2006), and random forest models (RF; Breiman, 2001). In addition to the existing presence points of the three species, we randomly generated 5000 background points for each species in China. Then, we randomly divided the dataset into a 30% test set and a 70% training set with 10 repetitions using bootstrap validation to reduce the uncertainty. The sdm package in R 4.2.1 was used to simulate the three SDMs models (Naimi and Araujo, 2016).
The accuracy of the three species distribution models was assessed mainly by four complementary metrics, including the area under the receiver operating characteristic (ROC) curve (AUC), true skill statistics (TSS), Sensitivity and Specificity (Allouche et al., 2006; Lobo et al., 2008; Liu et al., 2010). The AUC is a threshold-independent evaluation criterion (Hanley and McNeil, 1982), which is used to assess the ability of the predictive model (Phillips, 2008). In this study, three levels were used for model evaluation, as follows: AUC values below 0.7 were considered poor, values between 0.7 and 0.9 were considered moderate, and values greater than 0.9 were considered excellent (Li et al., 2020). In addition, we used another model evaluation metric, the TSS, which is calculated as the true positive rate minus the false positive rate (Allouche et al., 2006). TSS values between 0.4 and 0.6 were categorized as “fair,” while those between 0.6 and 1 were categorized as “good” (Xiao et al., 2019). Sensitivity is the probability that correctly predicts a presence, while Specificity is the probability that correctly predicts an absence (Liu et al., 2010). It is worth noting that the TSS value equals to Sensitivity value plus Specificity minus 1.
Finally, the model with the highest AUC, TSS, Sensitivity, and Specificity values was selected as the best predictive model. To further elucidate the effects of environmental factors on species distribution, we calculated the relative contribution of each predictor variable to the distribution model. Subsequently, response curves were generated to investigate the intrinsic relationships between the key factors affecting the potential distribution of the three species.
2.4 Evaluation of the dynamics of species habitat suitability
We used the best model to predict the spatial distribution of S. adnata, S. javanica, and S. williamsii under current and multiple future climate scenarios (SSP245-2050, SSP245-2090, SSP585-2050, and SSP585-2090). Using the species distribution model outputs, we transformed the continuous suitability scores (ranging from 0-1) into habitat suitability maps in ArcMap 10.8. To assess changes in the predicted distribution of each species, we constructed binary distribution maps for both current and future climate scenarios using the SDMtoolbox. The areas of suitable and unsuitable zones of three species were calculated. Firstly, we employed the sdm package in R to calculate the minimum suitability score threshold of each species. Below the score threshold represents unsuitable for the survival of the species (Absence), and greater than or equal to the threshold represents suitable for the survival (Presence). Secondly, we calculated the suitable and unsuitable survival areas in the SDMtoolbox of ArcGIS. Specifically, we quantified trends in future species suitability areas by calculating species distribution under future climate scenarios relative to current distributions. To evaluate changes in suitability, we classified the results into three categories: contraction trend, no change trend, and expansion trend.
3 Results
3.1 Model assessment
All three species (S. adnata, S. javanica, and S. williamsii) performed well under current climate conditions in all species distribution models (GLM, MaxEnt, and RF; AUC > 0.9, TSS > 0.6), except for the S. williamsii in the GLM model (AUC = 0.84, TSS = 0.56) (Table 1). The RF models had the best predictive performance for all three species, with AUC ranging from 0.96 to 0.98 and TSS ranging from 0.8 to 0.86, while the GLM models had lower suitability, with AUC ranging from 0.84 to 0.92 and TSS ranging from 0.56 to 0.75 (Table 1). In addition, the RF models had the highest Specificity and Sensitivity values, while the GLM models had the lowest Specificity and Sensitivity values. Therefore, we used the RF models to predict the current and potential future distribution of all three species.
Table 1 Accuracy evaluation of GLM, MaxEnt and RF models in predicting the potential distribution of S. adnata, S. javanica and S. williamsii in China.
3.2 Influence of predictor variables
Overall, the geographic distribution of S. adnata, S. javanica, and S. williamsii was primarily influenced by climatic factors, followed by soil and topography factors (Figure 2). Specifically, for S. adnata, the four most important predictor variables were temperature seasonality (bio4, 47.6%), soil total nitrogen (soil_TN, 17.4%), mean temperature of the wettest quarter (bio8, 10%), and slope (topo_slope, 9%) (Figure 2A). To understand how each key environmental variable influenced potential distribution of all three species, we plotted variable response curves (Figure 3). The probability of occurrence of S. adnata was negatively correlated with bio4, with an exception of a slight initial positive correlation (Figure 3A). In contrast, within a specific range, the probability of occurrence of S. adnata was positively correlated with the soil_TN, topo_slope, and bio8 (Figure 3A).
Figure 2 The relative contribution of environmental variables used in this study in predicting the current distribution of S. adnata (A), S. javanica(B), and S. williamsii (C) in China.
Figure 3 The response curves for the four most important contributing environment variables for S. adnata (A), S. javanica (B), and S. williamsii (C) under the current climate conditions in China, which reflect the relationship between environmental variables and probability of occurrence of species.
For S. javanica, the most important factors for potential distribution were the mean diurnal range (bio2, 10.4%), mean temperature of the driest quarter (bio9, 9.1%), slope (topo_slope, 7.7%), and temperature seasonality (bio4, 7.2%) (Figure 3B). There was a positive relationship between the occurrence probability for S. javanica and bio9 and topo_slope, while a negative relationship between bio2 and bio4 (Figure 3B).
For S. williamsii, the four most dominant contributing environmental variables were the precipitation of wettest month (bio13, 23.9%), temperature seasonality (bio4, 13.5%), isothermality (bio3, 11.1%), and slope (topo_slope, 11.1%). Moreover, we found that there was a unimodal relationship bio13 and the probability of occurrence of S. williamsii (Figure 3C). Meanwhile, the probability of occurrence of S. williamsii was positively associated with the topo_slope, but negatively correlated with bio3 and bio4 (Figure 3C).
3.3 Predicting distribution under current and future conditions
The distribution patterns of S. adnata, S. javanica, and S. williamsii varied widely across China under current climatic conditions (Figure 4). Overall, all three species have high habitat suitability scores, with values greater than or equal to 0.8. Specifically, S. adnata was mainly distributed in northwest and southwest China, with only small areas in central south China (Figures 1B, 4A). The highly suitable areas for S. adnata were located throughout the southwest China, extending north to the southern part of northwest China (Figures 1B, 4A). Moreover, S. javanica had suitable areas distributed widely in almost all southern regions of China. Notably, the highly suitable areas for S. javanica occurred at the boundaries of three subregions: the northwest and southwest edges of China, the southwest and central south edges of China, and the entire Taiwan Province (Figures 1B, 4B). In addition, S. williamsii was mainly located in north, east, southwest and central south China, with only a small portion in northwest China (Figures 1B, 4C). The highly suitable areas were in the border of southwest and central south China, and a small southern region of northwest China.
Figure 4 The potential distributions of S. adnata (A), S. javanica (B), and S. williamsii (C) under current climate conditions in China. The legend colors indicate the differences of presence probability of species in different area.
The potential distributions of the three species under the four future climate scenarios (SSP245-2050, SSP245-2090, SSP585-2050, and SSP585-2090) were highly divergent compared to their current distributions (Figure 5). Specifically, the potential suitable distribution of S. adnata was expected to largely expand under future climate scenarios, with significant increases in relatively high suitability areas in northeast and north China (Figures 1B, 5A). Similarly, except for northwest China, the future potential suitability of S. javanica would be expanded in other regions of China (Figures 1B, 5B). However, the highly potentially suitable areas of S. javanica would basically be maintained in the current distribution areas (Figure 5B). Furthermore, the potential future suitability of S. williamsii would be expanded in northeast, north, and east China. Notably, the potentially suitable areas of S. williamsii increased most significantly under the SSP585-2090 scenario, especially in northeast China (Figures 1B, 5C).
Figure 5 The potential distributions of S. adnata (A), S. javanica (B), and S. williamsii (C) under SSP245-2050, SSP585-2050, SSP245-2090 and SSP585-2090 scenarios in China. The legend colors indicate the differences of the presence probability of species in a different area.
To further quantify the areas of suitable regions, we mapped the binary distribution of the three species under current and future climate scenarios (SSP245-2050, SSP245-2090, SSP585-2050, and SSP585-2090) (Figure 6). The areas of suitable zones for S. adnata under current and four future climate scenarios were calculated: 1.99×105, 7.14×105, 5.59×105, 5.59×105, and 8.75×105 km2, respectively (Figure 6A, Table S2). Similarly, the areas of suitable zones for S. javanica under current and future climate scenarios were 6.41×105, 9.46×105, 13.54×105, 13.54×105, and 13.54×105 km2, respectively (Figure 6B, Table S2). Finally, the areas of suitable zones for S. williamsii under current and future climate scenarios were 9.57×105, 10.79×105, 10.79×105, 3.55×105, 21.68×105 km2, respectively (Figure 6C, Table S2).
Figure 6 The presence and absence distributions of S. adnata (A), S. javanica (B), and S. williamsii (C) under current and future climatic conditions in China. The future climate conditions including SSP245-2050, SSP585-2050, SSP245-2090 and SSP585-2090 scenarios.
3.4 Future changes in suitable habitats
Our analysis of four future climate scenarios revealed significant differences in the expansion and contraction trends of all three species (Figure 7, Table S3). The expansion of S. adnata was observed in northeast China, while its contraction occurred in current high suitable zone (Figure 7A). Similarly, S. javanica’s expansion was seen in the southern part of China, while its contraction occurred in the current high suitability area. (Figure 7B). In addition, the potentially suitable distribution areas of S. williamsii varied significantly under different future climate scenarios (Figure 7C). We observed a contraction trend in north and parts of northwest China, and the borders of southwest and central south China under SSP245-2050, SSP245-2090, and SSP585-2050 scenarios. Conversely, northeast China showed a large expansion trend under the SSP585-2090 scenario.
Figure 7 Changes in the distribution of S. adnata (A), S. javanica (B), and S. williamsii (C) under SSP245-2050, SSP585-2050, SSP245-2090, and SSP585-2090 scenarios compared to current distribution in China. The legend colors indicate as follows: indicolite green denotes unchanged areas, red indicates expansion areas, and blue represents contraction areas.
4 Discussion
The genus Sambucus has significant roles in traditional Chinese medicine due to its pharmacological properties, such as the treatment of respiratory illness (Wieland et al., 2021), bone fractures (Liao et al., 2005; Lin et al., 2009; Wang et al., 2020) and diabetes (Waswa et al., 2022). In this study, using the best predictive model (RF model), we mapped, for the first time, the potential suitable distribution of three Sambucus species, S. adnata, S. javanica, and S. williamsii, under the current and future climate scenarios in China. Moreover, we further explored potential drivers influencing the habitat distribution of these three Sambucus species.
4.1 Model predictive performance
According to the evaluation of AUC, TSS, Specificity, and Sensitivity values, our results indicated that the GLM model was less predictive than the MaxEnt and RF models for the three Sambucus species (Table 1). This suggests that GLM models are more appropriate for analyzing linear correlations, but cannot adequately capture nonlinear relationships between high-dimensional variables (Elith and Graham, 2009; Mainali et al., 2015). In contrast, MaxEnt models offer several advantages for species distribution prediction, such as an easily interpreted modeling mechanism (Zhang and Wang, 2022) and not requiring background data (Elith et al., 2011). However, it is worth noting that the parameter settings of the MaxEnt model need to be optimized for accurate prediction results (Radosavljevic and Anderson, 2014; Morales et al., 2017; Lissovsky and Dudov, 2021). Compared to GLM and MaxEnt models, our findings indicated that that RF performed the best in predicting the potentially suitable distribution of Sambucus, which is consistent with previous studies (González-Irusta et al., 2015; Mi et al., 2017; Čengić et al., 2020). Notably, RF models offer the benefits of computational efficiency and reliable prediction results without parameter adjustment, which are applicable to both large and small samples (Luan et al., 2020; Zhang et al., 2020; Zhao et al., 2022). Overall, our study emphasizes the importance of selecting appropriate modeling techniques for predicting the potential distribution of species.
4.2 Habitat distribution and key environmental factors under current environment
Environmental heterogeneity can influence distribution and abundance of species by altering conditions in temperature, precipitation, light, and nutrient availability (Palmer and Dixon, 1990; Wang et al., 2012; Ehrlén and Morris, 2015). Our findings indicated that climate drivers were more influential in determining the habitat distribution of S. adnata, S. javanica, and S. williamsii than topographic and soil drivers, which is consistent with recent studies (Dyderski et al., 2018; Yue et al., 2019; Naudiyal et al., 2021; Wang et al., 2022).
Specifically, we observed a strong negative effect of the temperature seasonality range (bio4) on the geographical distribution of all three species, particularly S. adnata (Figure 3). The observed distribution of Sambucus species in central, and southwestern China underscores their preference for stable temperature habitats (Figure 4). Furthermore, Sambucus plants exhibit a propensity for thriving within forested environments, benefiting from the temperature-buffering effect the forest canopy offers (De Lombaerde et al., 2022). Moreover, our findings revealed that the mean temperature of the wettest quarter (bio8) and the mean temperature of the driest quarter (bio9) played pivotal roles in shaping the distribution of S. adnata and S. javanica, respectively. This finding aligns with previous studies, which have consistently highlighted the positive influence of temperature on plant growth by extending the growing season and accelerating the phenological process (Zhang et al., 2015; Ma and Sun, 2018; Zhang et al., 2018). In addition, recent studies have also emphasized the importance of precipitation as a major factor influencing the distribution and survival of plant species (Petrie et al., 2016; Wang et al., 2018). In contrast to the responses of S. adnata and S. javanica to temperature, S. williamsii was more sensitive to variations in precipitation than temperature (Figure 3). Similar to a recent study (Naudiyal et al., 2021), the precipitation of wettest month (bio13) was identified as the most influential variable in determining the stable habitat distribution of S. williamsii, spanning a range from 0 to 270 mm. This underscores the notion that moderate precipitation levels confer benefits to plant growth and seedling survival (Petrie et al., 2016), while excessive soil water can create anaerobic conditions, detrimentally impacting plant metabolism and constraining growth (Zhang et al., 2019). Collectively, these findings underscore the pivotal role of temperature and precipitation as prerequisites for suitable distributions of Sambucus species, concurrently highlighting variations in ecological suitability among the three species.
In addition to climatic factors, topographic variables, especially slope, emerged as pivotal determinants shaping the habitat distributions of S. adnata, S. javanica and S. williamsii in the central and southwest China (Figure 3), where the genus Sambucus was mainly distributed (Figure 4). This aligns with previous research that underscores the provision of ample light to understory vegetation on steep slopes, facilitating their growth and regeneration (Austin and Van Niel, 2011; Chen et al., 2018). Additionally, our findings highlighted the varying impact of soil organic nitrogen (soil_TN) on the distribution patterns of the three species. Notably, S. adnata displayed a stronger positive correlation with soil_TN than S. javanica and S. williamsii, supporting prior studies that have elucidated the differential influence of soil nutrient availability on species distributions (Swaine, 1996; Lévesque et al., 2016; Walthert and Meier, 2017). Therefore, regions characterized by high soil nutrient availability are more suitable for the distribution of S. adnata.
It is worth noting that global climate change profoundly impacts soil properties across various aspects (Certini and Scalenghe, 2023). For instance, recent studies have demonstrated effects on soil pH patterns across different depths and grassland types, as well as on soil nutrient cycling processes (Gao et al., 2014; Sun et al., 2021b). Moreover, the regulation of moisture and temperature by climate change has a direct impact on soil biogeochemical responses, which in turn affect plant growth and distribution (Classen et al., 2015; Certini and Scalenghe, 2023).
4.3 Future distribution of Sambucus under climate change
Several recent studies have provided evidence indicating that climate change would significantly alter species ranges (Chen et al., 2011; Scheffers et al., 2016; Li et al., 2019), resulting in species migrating northward or to higher latitudes (Hickling et al., 2006; Bertrand et al., 2011). To simulate the maximum potential impact of greenhouse gas emissions on the Sambucus species, their potential distributions were predicted by our most robust prediction model (i.e., RF) across multiple future climate scenarios. Our findings indicate that the potentially suitable habitats of S. adnata, S. javanica, and S. williamsii are differ significantly under the impacts of climate change.
The potential suitable areas of S. adnata and S. javanica were predicted to expand significantly under future climate scenarios, but the directions of expansion were different. Specifically, S. adnata would shift its potential habitats to higher latitudes, while S. javanica would shift to lower latitudes. Contrary to previous studies that found negative effects of the habitat degeneration or population loss due to climate warming (Brecka et al., 2018; Dyderski et al., 2018; Sun et al., 2021a), both S. adnata and S. javanica are expected to benefit from global warming, possibly due to their broad thermal tolerance ranges (Calosi et al., 2008). Studies have shown that species with narrow thermal tolerances are more vulnerable to climate change than those with broader thermal ranges (Jiguet et al., 2006; Khaliq et al., 2014).
With global warming occurs, animals and plants have been observed to shift their distribution areas to higher latitudes or elevations to adapt to the changing climate (Parmesan, 2006; Chen et al., 2011; Feeley, 2012; Lenoir and Svenning, 2013). However, the direction and magnitude of changes in the geographic distribution of different species vary considerably (Lenoir et al., 2008; Fei et al., 2017). The distributions of species are not regulated by temperature alone (Wisz et al., 2013; Grytnes et al., 2014), but are driven by a combination of ecological factors, particularly the combined effects of temperature and moisture conditions (Zu et al., 2021). Recent studies revealed that the distinct roles played by temperature and precipitation in shaping the geographic distribution of species determine the direction of species shift under climate change (Zu et al., 2021; Chen et al., 2022).
However, the potential suitable areas of S. williamsii are expected to exhibit significantly different contraction and expansion trends under various future climate scenarios. There will be a noticeable trend of contraction in parts of northwest China, and the borders of southwest and central south China under SSP245-2050, SSP585-2050, and SSP245-2090 scenarios, but a trend of expansion in northeast China under the SSP585-2090 scenario. Compared to the SSP245 scenario, future climate change under the SSP585 scenario will become increasingly complicated over time (Liu et al., 2022). By the end of this century, the potential suitable habitats of S. williamsii are predicted to contract under the SSP245 scenario, but to expand under the SSP585 scenario. This illustrates the dual impact of climate change on suitable habitats of organisms (Liu et al., 2022). However, the complexity of climate change also contributes to uncertainty in species distribution projections (Millar et al., 2007; Yousefpour et al., 2012; Lindner et al., 2014).
4.4 Prospect
Our study provides a framework to predict the current and future geographic distribution of plant species based on climate, soil, and topography factors. However, there are still some uncertainties in our findings. First, the accuracy of the model is critical to the prediction results and can be improved through parameter adjustments (Radosavljevic and Anderson, 2014). Second, while climate, soil, and topography factors play key roles in species distribution, other factors such as human activities, land use change (Yu et al., 2019), biological interaction (Boulangeat et al., 2012), dispersal distances and limitations (Peterson et al., 2008), also impact species distribution patterns. Due to the inherent challenges in quantifying certain aspects and the potential issue of multicollinearity arising from incorporating all variables in the model, our study primarily focused on examining the effects of key variables. Therefore, our results may exhibit slight variations from the actual distribution pattern of the three species, and subsequent studies can supplement our findings.
Notably, the impact of human activities on plant species distribution cannot be overlooked, as they have become increasingly influential in shaping suitable habitats for animals and plants (Xu et al., 2019). Specifically, the impact of human activities on the diversity and distribution of plant communities can be either positive or negative (Wang et al., 2022; Huang and Fu, 2023). Moreover, the structure and function of different plant communities show considerable variation in response to human activities (Sun et al., 2021b; Zha et al., 2022). Therefore, it is crucial to consider the effects of human activities on species distributions in future studies.
Despite these uncertainties, we have confidence in the predictions of the RF model, which provides scientific theoretical foundations for the conservation and exploitation of the three medicinal plants from Sambucus. Therefore, we suggest that the contraction areas of S. adnata (southwest China), S. javanica (southwest China), and S. williamsii (north, central south, and southwest China) should be prioritized for germplasm resource collection. The expansion areas of S. adnata (northeast China), S. javanica (east and central south China), and S. williamsii (northeast China) may be suitable for introduction and cultivation.
5 Conclusions
This study demonstrated the significant ecological habitat differences among S. adnata, S. javanica, and S. williamsii in China. Furthermore, our findings underscore the pivotal roles of temperature and precipitation in shaping the potentially suitable habitats for these species. Specifically, S. adnata displayed heightened sensitivity to temperature seasonality, while S. javanica exhibited a strong dependence on mean diurnal range, and S. williamsii demonstrated a pronounced reliance on precipitation of the wettest month. Currently, S. williamsii primarily occupies north and central south China (9.57 × 105 km2), S. javanica is prevalent in the south and east regions (6.41×105 km2), and S. adnata predominantly occupies southwest areas (1.99×105 km2). Under future climate scenarios, we anticipated that S. adnata may shift to higher latitudes, while S. javanica may shift to lower latitudes. Our findings provide valuable insights into the potential impacts of climate change on the distributions of three Sambucus species in China. These insights can inform the development of effective conservation strategies, ensuring the preservation and sustainable management of these valuable medicinal plant resources.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
WL: Conceptualization, Methodology, Formal analysis, Writing – review & editing. SH: Investigation, Formal analysis, Data curation, Writing – original draft. TY: Investigation, Writing – review & editing. PW: Methodology, Writing – review & editing. YM: Investigation, Writing – review & editing. MW: Investigation, Writing – review & editing. JL: Supervision, Writing – review & editing. ZL: Supervision, Writing – review & editing. JT: Conceptualization, Validation, Funding acquisition, Writing – review & editing.
Funding
This research is supported by the Program of the National Natural Science Foundation of China (32201312), the Research Innovation Program for graduate students of Chongqing (5250101036), and the Fundamental Research Funds for the Central Universities (SWU-KQ22009).
Acknowledgments
We would like to thank all the students who participated in the data collection.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1194444/full#supplementary-material
References
Aiello-Lammens, M. E., Boria, R. A., Radosavljevic, A., Vilela, B., Anderson, R. P. (2015). spThin: an R package for spatial thinning of species occurrence records for use in ecological niche models. Ecography 38 (5), 541–545. doi: 10.1111/ecog.01132
Allouche, O., Tsoar, A., Kadmon, R. (2006). Assessing the accuracy of species distribution models: prevalence, kappa and the true skill statistic (TSS). J. Appl. Ecol. 43 (6), 1223–1232. doi: 10.1111/j.1365-2664.2006.01214.x
Austin, M. P., Van Niel, K. P. (2011). Improving species distribution models for climate change studies: variable selection and scale. J. Biogeogr. 38 (1), 1–8. doi: 10.1111/j.1365-2699.2010.02416.x
Becklin, K. M., Anderson, J. T., Gerhart, L. M., Wadgymar, S. M., Wessinger, C. A., Ward, J. K. (2016). Examining plant physiological responses to climate change through an evolutionary lens. Plant Physiol. 172 (2), 635–649. doi: 10.1104/pp.16.00793
Bertrand, R., Lenoir, J., Piedallu, C., Riofrio-Dillon, G., de Ruffray, P., Vidal, C., et al. (2011). Changes in plant community composition lag behind climate warming in lowland forests. Nature 479 (7374), 517–520. doi: 10.1038/nature10548
Boulangeat, I., Gravel, D., Thuiller, W. (2012). Accounting for dispersal and biotic interactions to disentangle the drivers of species distributions and their abundances. Ecol. Lett. 15 (6), 584–593. doi: 10.1111/j.1461-0248.2012.01772.x
Brecka, A. F. J., Shahi, C., Chen, H. Y. H. (2018). Climate change impacts on boreal forest timber supply. For. Policy Econ. 92, 11–21. doi: 10.1016/j.forpol.2018.03.010
Buri, A., Cianfrani, C., Pinto-Figueroa, E., Yashiro, E., Spangenberg, J. E., Adatte, T., et al. (2017). Soil factors improve predictions of plant species distribution in a mountain environment. Prog. Phys. Geogr. 41 (6), 703–722. doi: 10.1177/0309133317738162
Calosi, P., Bilton, D. T., Spicer, J. I. (2008). Thermal tolerance, acclimatory capacity and vulnerability to global climate change. Biol. Lett. 4 (1), 99–102. doi: 10.1098/rsbl.2007.0408
Cantón, Y., Del Barrio, G., Solé-Benet, A., Lázaro, R. (2004). Topographic controls on the spatial distribution of ground cover in the Tabernas badlands of SE Spain. Catena 55 (3), 341–365. doi: 10.1016/s0341-8162(03)00108-5
Čengić, M., Rost, J., Remenska, D., Janse, J. H., Huijbregts, M. A., Schipper, A. M. (2020). On the importance of predictor choice, modelling technique, and number of pseudo-absences for bioclimatic envelope model performance. Ecol. Evol. 10 (21), 12307–12317. doi: 10.1002/ece3.6859
Certini, G., Scalenghe, R. (2023). The crucial interactions between climate and soil. Sci. Total Environ. 856, 159169. doi: 10.1016/j.scitotenv.2022.159169
Change, I. C. (2013). The physical science basis. Contribution working group I to fifth Assess. Rep. intergovernmental panel Climate Change 1535, 2013.
Chen, I. C., Hill, J. K., Ohlemuller, R., Roy, D. B., Thomas, C. D. (2011). Rapid range shifts of species associated with high levels of climate warming. Science 333 (6045), 1024–1026. doi: 10.1126/science.1206432
Chen, Y. G., Le, X. G., Chen, Y. H., Cheng, W. X., Du, J. G., Zhong, Q. L., et al. (2022). [Identification of the potential distribution area of Cunninghamia lanceolata in China under climate change based on the MaxEnt model]. Ying Yong Sheng Tai Xue Bao 33 (5), 1207–1214. doi: 10.13287/j.1001-9332.202205.024
Chen, F., Liu, D. L., Wang, W., Lv, X. M., Li, W., Shao, L. D., et al. (2020). Bioactive triterpenoids from Sambucus javanica Blume. Nat. Prod. Res. 34 (19), 2816–2821. doi: 10.1080/14786419.2019.1596092
Chen, Y., Yuan, Z., Bi, S., Wang, X., Ye, Y., Svenning, J. C. (2018). Macrofungal species distributions depend on habitat partitioning of topography, light, and vegetation in a temperate mountain forest. Sci. Rep. 8 (1), 13589. doi: 10.1038/s41598-018-31795-7
Chytrý, M., Tichý, L., Roleček, J. (2003). Local and regional patterns of species richness in Central European vegetation types along the pH/calcium gradient. Folia Geobot. 38 (4), 429–442. doi: 10.1007/bf02803250
Classen, A. T., Sundqvist, M. K., Henning, J. A., Newman, G. S., Moore, J. A. M., Cregger, M. A., et al. (2015). Direct and indirect effects of climate change on soil microbial and soil microbial-plant interactions: What lies ahead? Ecosphere. 6 (8), 1–21. doi: 10.1890/es15-00217.1
Cruz-Cárdenas, G., López-Mata, L., Villaseñor, J. L., Ortiz, E. (2014). Potential species distribution modeling and the use of principal component analysis as predictor variables. Rev. Mex. Biodivers. 85 (1), 189–199. doi: 10.7550/rmb.36723
Dawson, T. P., Jackson, S. T., House, J. I., Prentice, I. C., Mace, G. M. (2011). Beyond predictions: biodiversity conservation in a changing climate. Science 332 (6025), 53–58. doi: 10.1126/science.1200303
De Lombaerde, E., Vangansbeke, P., Lenoir, J., Van Meerbeek, K., Lembrechts, J., Rodriguez-Sanchez, F., et al. (2022). Maintaining forest cover to enhance temperature buffering under future climate change. Sci. Total Environ. 810, 151338. doi: 10.1016/j.scitotenv.2021.151338
Dyderski, M. K., Paz, S., Frelich, L. E., Jagodzinski, A. M. (2018). How much does climate change threaten European forest tree species distributions? Glob. Change Biol. 24 (3), 1150–1163. doi: 10.1111/gcb.13925
Ehrlén, J., Morris, W. F. (2015). Predicting changes in the distribution and abundance of species under environmental change. Ecol. Lett. 18 (3), 303–314. doi: 10.1111/ele.12410
Elith, J., Graham, C. H. (2009). Do they? How do they? WHY do they differ? On finding reasons for differing performances of species distribution models. Ecography 32 (1), 66–77. doi: 10.1111/j.1600-0587.2008.05505.x
Elith, J., Phillips, S. J., Hastie, T., Dudik, M., Chee, Y. E., Yates, C. J. (2011). A statistical explanation of MaxEnt for ecologists. Divers. Distrib. 17 (1), 43–57. doi: 10.1111/j.1472-4642.2010.00725.x
Evans, J. S., Murphy, M. A., Holden, Z. A., Cushman, S. A. (2011). “Modeling species distribution and change using random forest,” in Predictive species and habitat modeling in landscape ecology (New York: Springer), 139–159.
Feeley, K. J. (2012). Distributional migrations, expansions, and contractions of tropical plant species as revealed in dated herbarium records. Glob. Change Biol. 18 (4), 1335–1341. doi: 10.1111/j.1365-2486.2011.02602.x
Fei, S., Desprez, J. M., Potter, K. M., Jo, I., Knott, J. A., Oswalt, C. M. (2017). Divergence of species responses to climate change. Sci. Adv. 3 (5), e1603055. doi: 10.1126/sciadv.1603055
Gábor, L., Šímová, P., Keil, P., Zarzo-Arias, A., Marsh, C. J., Rocchini, D., et al. (2022). Habitats as predictors in species distribution models: Shall we use continuous or binary data? Ecography 2022 (7), e06022. doi: 10.1111/ecog.06022
Gao, Q.-z., Li, Y., Xu, H.-m., Wan, Y.-f., Jiangcun, W.-z. (2014). Adaptation strategies of climate variability impacts on alpine grassland ecosystems in Tibetan Plateau. Mitig. Adapt. Strateg. Glob. Change. 19, 199–209. doi: 10.1007/s11027-012-9434-y
González-Irusta, J. M., González-Porto, M., Sarralde, R., Arrese, B., Almón, B., Martín-Sosa, P. (2015). Comparing species distribution models: a case study of four deep sea urchin species. Hydrobiologia 745 (1), 43–57. doi: 10.1007/s10750-014-2090-3
Grytnes, J. A., Kapfer, J., Jurasinski, G., Birks, H. H., Henriksen, H., Klanderud, K., et al. (2014). Identifying the driving factors behind observed elevational range shifts on E uropean mountains. Glob. Ecol. Biogeogr. 23 (8), 876–884. doi: 10.1111/geb.12170
Hanley, J. A., McNeil, B. J. (1982). The meaning and use of the area under a receiver operating characteristic (ROC) curve. Radiology 143 (1), 29–36. doi: 10.1148/radiology.143.1.7063747
Hengl, T., Mendes de Jesus, J., Heuvelink, G. B., Ruiperez Gonzalez, M., Kilibarda, M., Blagotić, A., et al. (2017). SoilGrids250m: Global gridded soil information based on machine learning. PloS One 12 (2), e0169748. doi: 10.1371/journal.pone.0169748
Hernandez, P. A., Graham, C. H., Master, L. L., Albert, D. L. (2006). The effect of sample size and species characteristics on performance of different species distribution modeling methods. Ecography 29 (5), 773–785. doi: 10.1111/j.0906-7590.2006.04700.x
Hickling, R., Roy, D. B., Hill, J. K., Fox, R., Thomas, C. D. (2006). The distributions of a wide range of taxonomic groups are expanding polewards. Glob. Change Biol. 12 (3), 450–455. doi: 10.1111/j.1365-2486.2006.01116.x
Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G., Jarvis, A. (2005). Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 25 (15), 1965–1978. doi: 10.1002/joc.1276
Hirzel, A. H., Hausser, J., Chessel, D., Perrin, N. (2002). Ecological-niche factor analysis: how to compute habitat-suitability maps without absence data? Ecology 83 (7), 2027–2036. doi: 10.1890/0012-9658(2002)083[2027:enfaht]2.0.co;2
Huang, S. L., Fu, G. (2023). Impacts of climate change and human activities on plant species α-diversity across the tibetan grasslands. Remote Sens. 15 (11), 2947. doi: 10.3390/rs15112947
Jarvis, A., Reuter, H. I., Nelson, A., Guevara, E. (2008). Hole-filled SRTM for the globe Version 4, available from the CGIAR-CSI SRTM 90m Database. 15 (25-54), 5.
Jiguet, F., Julliard, R., Thomas, C. D., Dehorter, O., Newson, S. E., Couvet, D. (2006). Thermal range predicts bird population resilience to extreme high temperatures. Ecol. Lett. 9 (12), 1321–1330. doi: 10.1111/j.1461-0248.2006.00986.x
John, R., Dalling, J. W., Harms, K. E., Yavitt, J. B., Stallard, R. F., Mirabello, M., et al. (2007). Soil nutrients influence spatial distributions of tropical tree species. Proc. Natl. Acad. Sci. U.S.A. 104 (3), 864–869. doi: 10.1073/pnas.0604666104
Kaky, E., Gilbert, F. (2017). Predicting the distributions of Egypt's medicinal plants and their potential shifts under future climate change. PloS One 12 (11), e0187714. doi: 10.1371/journal.pone.0187714
Khaliq, I., Hof, C., Prinzinger, R., Bohning-Gaese, K., Pfenninger, M. (2014). Global variation in thermal tolerances and vulnerability of endotherms to climate change. Proc. Natl. Acad. Sci. U. S. A. 281 (1789), 20141097. doi: 10.1098/rspb.2014.1097
Kim, J. B., Bae, D. H. (2021). The impacts of global warming on climate zone changes over Asia based on CMIP6 projections. Earth Space Sci. 8 (8), 1–17. doi: 10.1029/2021ea001701
Lehmann, A., Overton, J. M., Leathwick, J. R. (2002). GRASP: generalized regression analysis and spatial prediction. Ecol. Model. 157 (2-3), 189–207. doi: 10.1016/s0304-3800(02)00195-3
Lenoir, J., Gegout, J. C., Marquet, P. A., de Ruffray, P., Brisse, H. (2008). A significant upward shift in plant species optimum elevation during the 20th century. Science 320 (5884), 1768–1771. doi: 10.1126/science.1156831
Lenoir, J., Svenning, J.-C. (2013). Latitudinal and elevational range shifts under contemporary climate change. Encyclopedia Biodiversity (Second Edition) 599-611. doi: 10.1016/B978-0-12-384719-5.00375-0
Lévesque, M., Walthert, L., Weber, P. (2016). Soil nutrients influence growth response of temperate tree species to drought. J. Ecol. 104 (2), 377–387. doi: 10.1111/1365-2745.12519
Li, J. Y., Chang, H., Liu, T., Zhang, C. (2019). The potential geographical distribution of Haloxylon across Central Asia under climate change in the 21st century. Agric. For. Meteorol. 275, 243–254. doi: 10.1016/j.agrformet.2019.05.027
Li, J., Fan, G., He, Y. (2020). Predicting the current and future distribution of three Coptis herbs in China under climate change conditions, using the MaxEnt model and chemical analysis. Sci. Total Environ. 698, 134141. doi: 10.1016/j.scitotenv.2019.134141
Li, Q. Y., Wang, W., Li, L. H., Wang, W. J. (2021). Chemical components from Sambucus adnata wall. Biochem. Syst. Ecol. 96, 0305–1978. doi: 10.1016/j.bse.2021.104266
Liao, Q., Yang, W., Jia, Y., Chen, X., Gao, Q., Bi, K. (2005). LC-MS determination and pharmacokinetic studies of ursolic acid in rat plasma after administration of the traditional chinese medicinal preparation Lu-Ying extract. Yakugaku Zasshi 125 (6), 509–515. doi: 10.1248/yakushi.125.509
Lin, M., Mei, J., Liu, X. (2009). The clinical observation of the effects of traditional chinese medicine on healing fracture. J. Mod Med. Health 25 (6), 901–902.
Lindner, M., Fitzgerald, J. B., Zimmermann, N. E., Reyer, C., Delzon, S., van der Maaten, E., et al. (2014). Climate change and European forests: what do we know, what are the uncertainties, and what are the implications for forest management? J. Environ. Manage. 146, 69–83. doi: 10.1016/j.jenvman.2014.07.030
Lissovsky, A., Dudov, S. (2021). Species-distribution modeling: advantages and limitations of its application. 2. MaxEnt. Biol. Bull. Rev. 11 (3), 265–275. doi: 10.1134/s2079086421030087
Liu, C., White, M., Newell, G. (2010). Measuring and comparing the accuracy of species distribution models with presence-absence data. Ecography 34 (2), 232–243. doi: 10.1111/j.1600-0587.2010.06354.x
Liu, J. M., Xu, Y. N., Sun, C. W., Wang, X., Zheng, Y. L., Shi, S. L., et al. (2022). Distinct ecological habits and habitat responses to future climate change in three east and southeast Asian Sapindus species. For. Ecol. Manage. 507, 119982. doi: 10.1016/j.foreco.2021.119982
Lobo, J. M., Jimenez-Valverde, A., Real, R. (2008). AUC: a misleading measure of the performance of predictive distribution models. Glob. Ecol. Biogeogr. 17 (2), 145–151. doi: 10.1111/j.1466-8238.2007.00358.x
Luan, J., Zhang, C. L., Xu, B. D., Xue, Y., Ren, Y. P. (2020). The predictive performances of random forest models with limited sample size and different species traits. Fisheries Res. 227, 105534. doi: 10.1016/j.fishres.2020.105534
Lv, H., Chen, S. S., Xu, X. L., Zhu, M. M., Zhao, W. F., Liu, K. W., et al. (2015). Isolation of linoleic acid from sambucus williamsii seed oil extracted by high pressure fluid and its antioxidant, antiglycemic, hypolipidemic activities. Int. J. Food Eng. 11 (3), 383–391. doi: 10.1515/ijfe-2014-0234
Ma, B., Sun, J. (2018). Predicting the distribution of Stipa purpurea across the Tibetan Plateau via the MaxEnt model. BMC Ecol. 18 (1), 10. doi: 10.1186/s12898-018-0165-0
Mainali, K. P., Warren, D. L., Dhileepan, K., McConnachie, A., Strathie, L., Hassan, G., et al. (2015). Projecting future expansion of invasive species: comparing and improving methodologies for species distribution modeling. Glob. Change Biol. 21 (12), 4464–4480. doi: 10.1111/gcb.13038
Maria, B., Udo, S. (2017). Why input matters: Selection of climate data sets for modelling the potential distribution of a treeline species in the Himalayan region. Ecol. Model. 359, 92–102. doi: 10.1016/j.ecolmodel.2017.05.021
Mi, C., Huettmann, F., Guo, Y., Han, X., Wen, L. (2017). Why choose Random Forest to predict rare species distribution with few samples in large undersampled areas? Three Asian crane species models provide supporting evidence. PeerJ 5, e2849. doi: 10.7717/peerj.2849
Millar, C. I., Stephenson, N. L., Stephens, S. L. (2007). Climate change and forests of the future: Managing in the face of uncertainty. Ecol. Appl. 17 (8), 2145–2151. doi: 10.1890/06-1715.1
Mohapatra, J., Singh, C. P., Hamid, M., Verma, A., Semwal, S. C., Gajmer, B., et al. (2019). Modelling Betula utilis distribution in response to climate-warming scenarios in Hindu-Kush Himalaya using random forest. Biodivers. Conserv. 28 (8-9), 2295–2317. doi: 10.1007/s10531-019-01731-w
Morales, N. S., Fernández, I. C., Baca-González, V. (2017). MaxEnt’s parameter configuration and small samples: are we paying attention to recommendations? A systematic review. PeerJ 5, e3093. doi: 10.7717/peerj.3093
Naimi, B., Araujo, M. B. (2016). sdm: a reproducible and extensible R platform for species distribution modelling. Ecography 39 (4), 368–375. doi: 10.1111/ecog.01881
Naimi, B., Hamm, N. A. S., Groen, T. A., Skidmore, A. K., Toxopeus, A. G. (2014). Where is positional uncertainty a problem for species distribution modelling? Ecography 37 (2), 191–203. doi: 10.1111/j.1600-0587.2013.00205.x
Naudiyal, N., Wang, J. N., Ning, W., Gaire, N. P., Peili, S., Wei, Y. Q., et al. (2021). Potential distribution of Abies, Picea, and Juniperus species in the sub-alpine forest of Minjiang headwater region under current and future climate scenarios and its implications on ecosystem services supply. Ecol. Indic. 121, 107131. doi: 10.1016/j.ecolind.2020.107131
Nelder, J. A., Wedderburn, R. W. (1972). Generalized linear models. J. R. Stat. Soc Ser. A-Stat. Soc 135 (3), 370–384. doi: 10.2307/2344614
O'Neill, B. C., Kriegler, E., Riahi, K., Ebi, K. L., Hallegatte, S., Carter, T. R., et al. (2014). A new scenario framework for climate change research: the concept of shared socioeconomic pathways. Clim. Change 122 (3), 387–400. doi: 10.1007/s10584-013-0905-2
Outlook, G. B. (2010). “Global biodiversity outlook 3,” in Montréal, Canada: Secretariat of the Convention on Biological Diversity (Phil. Trans. R. Soc. B). Available at: http://gbo3.cbd.int/.
Palmer, M. W., Dixon, P. M. (1990). Small-scale environmental heterogeneity and the analysis of species distributions along gradients. J. Veg. Sci. 1 (1), 57–65. doi: 10.2307/3236053
Parmesan, C. (2006). Ecological and evolutionary responses to recent climate change. Annu. Rev. Ecol. Evol. Syst. 37, 637–669. doi: 10.1146/annurev.ecolsys.37.091305.110100
Pearson, R. G., Raxworthy, C. J., Nakamura, M., Peterson, A. T. (2007). Predicting species distributions from small numbers of occurrence records: a test case using cryptic geckos in Madagascar. J. Biogeogr. 34 (1), 102–117. doi: 10.1111/j.1365-2699.2006.01594.x
Peterson, A. T., Papes, M., Soberon, J. (2008). Rethinking receiver operating characteristic analysis applications in ecological niche modeling. Ecol. Model. 213 (1), 63–72. doi: 10.1016/j.ecolmodel.2007.11.008
Petrie, M. D., Wildeman, A. M., Bradford, J. B., Hubbard, R. M., Lauenroth, W. K. (2016). A review of precipitation and temperature control on seedling emergence and establishment for ponderosa and lodgepole pine forest regeneration. For. Ecol. Manage. 361, 328–338. doi: 10.1016/j.foreco.2015.11.028
Phillips, S. J. (2008). Transferability, sample selection bias and background data in presence-only modelling: a response to Peterson et al (2007). Ecography 31 (2), 272–278. doi: 10.1111/j.0906-7590.2008.5378.x
Phillips, S. J., Anderson, R. P., Schapire, R. E. (2006). Maximum entropy modeling of species geographic distributions. Ecol. Model. 190 (3-4), 231–259. doi: 10.1016/j.ecolmodel.2005.03.026
Radosavljevic, A., Anderson, R. P. (2014). Making better MAXENT models of species distributions: complexity, overfitting and evaluation. J. Biogeogr. 41 (4), 629–643. doi: 10.1111/jbi.12227
Rather, T. A., Kumar, S., Khan, J. A. (2020). Multi-scale habitat modelling and predicting change in the distribution of tiger and leopard using random forest algorithm. Sci. Rep. 10 (1), 11473. doi: 10.1038/s41598-020-68167-z
Roman-Palacios, C., Wiens, J. J. (2020). Recent responses to climate change reveal the drivers of species extinction and survival. Proc. Natl. Acad. Sci. U. S. A. 117 (8), 4211–4217. doi: 10.1073/pnas.1913007117
Rousseau, J. S., Betts, M. G. (2022). Factors influencing transferability in species distribution models. Ecography 2022 (7), e06060. doi: 10.1111/ecog.06060
Scheffers, B. R., De Meester, L., Bridge, T. C., Hoffmann, A. A., Pandolfi, J. M., Corlett, R. T., et al. (2016). The broad footprint of climate change from genes to biomes to people. Science 354 (6313), aaf7671. doi: 10.1126/science.aaf7671
Shen, Y., Wang, G. (2013). Key findings and assessment results of IPCC WGI fifth assessment report. J. Glaciology Geocryology 35 (5), 1068–1076.
Shrestha, U. B., Sharma, K. P., Devkota, A., Siwakoti, M., Shrestha, B. B. (2018). Potential impact of climate change on the distribution of six invasive alien plants in Nepal. Ecol. Indic. 95, 99–107. doi: 10.1016/j.ecolind.2018.07.009
Soberon, J., Peterson, A. T. (2005). Interpretation of models of fundamental ecological niches and species’ Distributional areas. Biodivers. Inform. 2 (0), 1–10. doi: 10.17161/bi.v2i0.4
Stockwell, D., Peters, D. (1999). The GARP modelling system: problems and solutions to automated spatial prediction. Int. J. Geogr. Inf. Sci. 13 (2), 143–158. doi: 10.1080/136588199241391
Sun, J. J., Jiao, W. X., Wang, Q., Wang, T. L., Yang, H. Q., Jin, J. X., et al. (2021a). Potential habitat and productivity loss of Populus deltoides industrial forest plantations due to global warming. For. Ecol. Manage. 496, 119474. doi: 10.1016/j.foreco.2021.119474
Sun, W., Li, S., Wang, J., Fu, G. (2021b). Effects of grazing on plant species and phylogenetic diversity in alpine grasslands, Northern Tibet. Ecol. Eng. 170, 106331. doi: 10.1016/j.ecoleng.2021.106331
Svenning, J. C., Harlev, D., Sorensen, M., Balslev, H. (2009). Topographic and spatial controls of palm species distributions in a montane rain forest, southern Ecuador. Biodivers. Conserv. 18 (1), 219–228. doi: 10.1007/s10531-008-9468-3
Swaine, M. D. (1996). Rainfall and soil fertility as factors limiting forest species distributions in Ghana. J. Ecol. 84 (3), 419–428. doi: 10.2307/2261203
Walthert, L., Meier, E. S. (2017). Tree species distribution in temperate forests is more influenced by soil than by climate. Ecol. Evol. 7 (22), 9473–9484. doi: 10.1002/ece3.3436
Wang, Z. H., Fang, J. Y., Tang, Z. Y., Lin, X. (2012). Relative role of contemporary environment versus history in shaping diversity patterns of China's woody plants. Ecography 35 (12), 1124–1133. doi: 10.1111/j.1600-0587.2011.06781.x
Wang, J., Seyler, B. C., Ticktin, T., Zeng, Y., Ayu, K. (2020). An ethnobotanical survey of wild edible plants used by the Yi people of Liangshan Prefecture, Sichuan Province, China. J. Ethnobiol. Ethnomed. 16 (1), 10. doi: 10.1186/s13002-019-0349-5
Wang, D., Shi, C., Alamgir, K., Kwon, S., Pan, L., Zhu, Y., et al. (2022). Global assessment of the distribution and conservation status of a key medicinal plant (Artemisia annua L.): The roles of climate and anthropogenic activities. Sci. Total Environ. 821, 153378. doi: 10.1016/j.scitotenv.2022.153378
Wang, J., Sun, Z. Y., Hui, D. F., Yang, L., Wang, F. M., Liu, N., et al. (2018). Responses of seedling performance to altered seasonal precipitation in a secondary tropical forest, southern China. For. Ecol. Manage. 410, 27–34. doi: 10.1016/j.foreco.2017.12.035
Wang, S., Yu, Y., Cui, M., Liu, K., Liu, K. (2021). Seed oil quality and cultivation of sambucus williamsii hance as a new oil crop. Front. Nutr. 8. doi: 10.3389/fnut.2021.796175
Waswa, E. N., Li, J., Mkala, E. M., Wanga, V. O., Mutinda, E. S., Nanjala, C., et al. (2022). Ethnobotany, phytochemistry, pharmacology, and toxicology of the genus Sambucus L. (Viburnaceae). J. Ethnopharmacol. 292, 115102. doi: 10.1016/j.jep.2022.115102
Wieland, L. S., Piechotta, V., Feinberg, T., Ludeman, E., Hutton, B., Kanji, S., et al. (2021). Elderberry for prevention and treatment of viral respiratory illnesses: a systematic review. BMC Complement. Altern. Med. 21 (1), 112. doi: 10.1186/s12906-021-03283-5
Wisz, M. S., Pottier, J., Kissling, W. D., Pellissier, L., Lenoir, J., Damgaard, C. F., et al. (2013). The role of biotic interactions in shaping distributions and realised assemblages of species: implications for species distribution modelling. Biol. Rev. 88 (1), 15–30. doi: 10.1111/j.1469-185X.2012.00235.x
Xiao, J., Eziz, A., Zhang, H., Wang, Z., Tang, Z., Fang, J. (2019). Responses of four dominant dryland plant species to climate change in the Junggar Basin, northwest China. Ecol. Evol. 9 (23), 13596–13607. doi: 10.1002/ece3.5817
Xiao, H. H., Zhang, Y., Cooper, R., Yao, X. S., Wong, M. S. (2016). Phytochemicals and potential health effects of Sambucus williamsii Hance (Jiegumu). Chin. Med. 11 (1), 36. doi: 10.1186/s13020-016-0106-9
Xu, N., Meng, F. Y., Zhou, G. F., Li, Y. F., Wang, B., Lu, H. (2020). Assessing the suitable cultivation areas for Scutellaria baicalensis in China using the Maxent model and multiple linear regression. Biochem. Syst. Ecol. 90, 104052. doi: 10.1016/j.bse.2020.104052
Xu, W. B., Svenning, J. C., Chen, G. K., Zhang, M. G., Huang, J. H., Chen, B., et al. (2019). Human activities have opposing effects on distributions of narrow-ranged and widespread plant species in China. Proc. Natl. Acad. Sci. 116 (52), 26674–26681. doi: 10.1073/pnas.1911851116
Yi, Y. J., Cheng, X., Yang, Z. F., Zhang, S. H. (2016). Maxent modeling for predicting the potential distribution of endangered medicinal plant (H. riparia Lour) in Yunnan, China. Ecol. Eng. 92, 260–269. doi: 10.1016/j.ecoleng.2016.04.010
Yi, Y. J., Zhou, Y., Cai, Y. P., Yang, W., Li, Z. W., Zhao, X. (2018). The influence of climate change on an endangered riparian plant species: The root of riparian Homonoia. Ecol. Indic. 92, 40–50. doi: 10.1016/j.ecolind.2017.05.004
Yousefpour, R., Jacobsen, J. B., Thorsen, B. J., Meilby, H., Hanewinkel, M., Oehler, K. (2012). A review of decision-making approaches to handle uncertainty and risk in adaptive forest management under climate change. Ann. For. Sci. 69 (1), 1–15. doi: 10.1007/s13595-011-0153-4
Yu, F., Wang, T., Groen, T. A., Skidmore, A. K., Yang, X., Ma, K., et al. (2019). Climate and land use changes will degrade the distribution of Rhododendrons in China. Sci. Total Environ. 659, 515–528. doi: 10.1016/j.scitotenv.2018.12.223
Yuan, L., Zhong, Z. C., Liu, Y. (2020). Structural characterisation and immunomodulatory activity of a neutral polysaccharide from Sambucus adnata Wall. Int. J. Biol. Macromolecules 154, 1400–1407. doi: 10.1016/j.ijbiomac.2019.11.021
Yue, Y., Zhang, P., Shang, Y. (2019). The potential global distribution and dynamics of wheat under multiple climate change scenarios. Sci. Total Environ. 688, 1308–1318. doi: 10.1016/j.scitotenv.2019.06.153
Zha, X., Tian, Y., Ouzhu, N., Fu, G. (2022). Response of forage nutrient storages to grazing in alpine grasslands. Front. Plant Sci. 13, 991287. doi: 10.3389/fpls.2022.991287
Zhang, B., Wang, H. J. (2022). Exploring the advantages of the maximum entropy model in calibrating cellular automata for urban growth simulation: a comparative study of four methods. GISci. Remote Sens. 59 (1), 71–95. doi: 10.1080/15481603.2021.2016240
Zhang, H., Wang, X., Zhang, Y., Hu, R., Pan, Y., Chen, N. (2015). Responses of plant growth of different life forms to rainfall amount changes in an arid desert area. Chin. J. Ecol. 34 (7), 1847.
Zhang, X. H., Wei, H. Y., Zhao, Z. F., Liu, J., Zhang, Q. Z., Zhang, X. Y., et al. (2020). The Global Potential Distribution of Invasive Plants: Anredera cordifolia under Climate Change and Human Activity Based on Random Forest Models. Sustainability 12 (4), 1491. doi: 10.3390/su12041491
Zhang, K., Yao, L., Meng, J., Tao, J. (2018). Maxent modeling for predicting the potential geographical distribution of two peony species under climate change. Sci. Total Environ. 634, 1326–1334. doi: 10.1016/j.scitotenv.2018.04.112
Zhang, K. L., Zhang, Y., Zhou, C., Meng, J. S., Sun, J., Zhou, T. H., et al. (2019). Impact of climate factors on future distributions of Paeonia ostii across China estimated by MaxEnt. Ecol. Inform. 50, 62–67. doi: 10.1016/j.ecoinf.2019.01.004
Zhang, T. H., Zhu, M., Chen, X. H., Bi, K. S. (2010). Simultaneous analysis of seven bioactive compounds in sambucus chinensis lindl by HPLC. Anal. Lett. 43 (16), 2525–2533. doi: 10.1080/00032711003731399
Zhao, K., Li, X., Yang, J., Huang, Z., Li, C., Yao, L., et al. (2023). Effects of climate change on the geographical distribution and potential distribution areas of 35 Millettia Species in China. Environ. Sci. pollut. Res. 30 (7), 18535–18545. doi: 10.1007/s11356-022-23515-6
Zhao, Z. Y., Xiao, N. W., Shen, M., Li, J. S. (2022). Comparison between optimized MaxEnt and random forest modeling in predicting potential distribution: A case study with Quasipaa boulengeri in China. Sci. Total Environ. 842, 156867. doi: 10.1016/j.scitotenv.2022.156867
Keywords: climate change, Sambucus, medicinal plants, habitat distribution, random forest, species distribution model
Citation: Luo W, Han S, Yu T, Wang P, Ma Y, Wan M, Liu J, Li Z and Tao J (2023) Assessing the suitability and dynamics of three medicinal Sambucus species in China under current and future climate scenarios. Front. Plant Sci. 14:1194444. doi: 10.3389/fpls.2023.1194444
Received: 27 March 2023; Accepted: 02 October 2023;
Published: 19 October 2023.
Edited by:
Sailesh Ranjitkar, N. Gene Solution of Natural Innovation, NepalReviewed by:
Hamid Moazzeni, Ferdowsi University of Mashhad, IranGang Fu, Chinese Academy of Sciences (CAS), China
Khurram Shahzad, Chinese Academy of Sciences (CAS), China
Copyright © 2023 Luo, Han, Yu, Wang, Ma, Wan, Liu, Li and Tao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jianping Tao, taojp@swu.edu.cn
†Present address: Jianping Tao, School of Life Sciences, Southwest University, Chongqing, China
‡These authors share first authorship