- 1National Engineering Laboratory for Resource Development of Endangered Crude Drugs in Northwest China, The Key Laboratory of Medicinal Resources and Natural Pharmaceutical Chemistry, The Ministry of Education, College of Life Sciences, Shaanxi Normal University, Xi’an, China
- 2Henan Sesame Research Center, Henan Academy of Agricultural Sciences, Zhengzhou, Henan, China
- 3Shaanxi Panlong Pharmaceutical Group Limited by Share Ltd, Shangluo, Shaanxi, China
Schisandrae Sphenantherae Fructus (SSF), the dry ripe fruit of Schisandra sphenanthera Rehd. et Wils., is a traditional Chinese medicine with wide application potential. The quality of SSF indicated by the composition and contents of secondary metabolites is closely related to environmental factors, such as regional climate and soil conditions. The aims of this study were to predict the distribution patterns of potentially suitable areas for S. sphenanthera in China and pinpoint the major environmental factors influencing its accumulation of medicinal components. An optimized maximum entropy model was developed and applied under current and future climate scenarios (SSP1-RCP2.6, SSP3-RCP7, and SSP5-RCP8.5). Results show that the total suitable areas for S. sphenanthera (179.58×104 km2) cover 18.71% of China’s territory under the current climatic conditions (1981–2010). Poorly, moderately, and highly suitable areas are 119.00×104 km2, 49.61×104 km2, and 10.98×104 km2, respectively. The potentially suitable areas for S. sphenanthera are predicted to shrink and shift westward under the future climatic conditions (2041–2070 and 2071–2100). The areas of low climate impact are located in southern Shaanxi, northwestern Guizhou, southeastern Chongqing, and western Hubei Provinces (or Municipality), which exhibit stable and high suitability under different climate scenarios. The contents of volatile oils, lignans, and polysaccharides in SSF are correlated with various environmental factors. The accumulation of major secondary metabolites is primarily influenced by temperature variation, seasonal precipitation, and annual precipitation. This study depicts the potential distribution of S. sphenanthera in China and its spatial change in the future. Our findings decipher the influence of habitat environment on the geographical distribution and medicinal quality of S. sphenanthera, which could have great implications for natural resource conservation and artificial cultivation.
1 Introduction
Medicinal plants are a bioresource that provides pharmaceutically active components. The distribution and growth of medicinal plants are dependent on a range of environmental factors, such as geographical location, climate, and soil conditions (Shan et al., 2022; Feng et al., 2023). Various medicinal plant species show distinct habitat preferences. For example, Scutellaria baicalensis Georgi is tolerant to low temperature and drought, with hindered growth in hot and humid regions (Ji et al., 2021). Atractylodes lancea (Thunb.) DC. is also widely distributed in arid regions (Han et al., 2020), in contrast to Paris yunnanensis Franch. that generally favors humid climate zones with concentrated precipitation (Luo et al., 2017). Even cultivars of the same medicinal plant species (e.g., Anemarrhena asphodeloides Bunge.) vary in yield and quality across regions (Li et al., 2015). Therefore, the formulation of effective strategies for the exploitation and utilization of medicinal plants necessitates accurate prediction of their potentially suitable habitats.
Ecological niche modeling (ENM) has been widely used to predict the potential distribution of various species and assess their habitat suitability in the natural environment (Blank and Blaustein, 2012). Owing to the development of geographic information systems (GIS), the potentially suitable habitats for a given species can be predicted by analyzing its regional environmental conditions based on the correlation between environmental factors and species occurrence data (Title and Bemmels, 2018). Among the ENM methods, maximum entropy (MaxEnt) modeling has an outstanding predictive power. In particular, the MaxEnt method can predict the distribution of potentially suitable habitats for a specific species in the past and future based on climate data. This would help identify the areas of low climate impact (i.e., low impact areas) for the target species (Xu et al., 2019). MaxEnt model can handle large-capacity data, and can still show excellent prediction ability in the case of small sample size and spatial deviation of distribution point data (Graham et al., 2008). As one of the most popular modeling approaches, the MaxEnt method has been successfully used in the prediction of potentially suitable areas for various species and the conservation planning of medicinal plants (Kaky et al., 2020; Chen et al., 2021).
Medicinal plants synthesize a variety of secondary metabolites (SMs) in response to their habitat environments. Natural SMs with diverse functions enable the plants to defend against pathogenic attacks and environmental stresses (Zhou et al., 2010). Given their remarkable biological activities, many SMs are used as active components in pharmaceutics, cosmetics, fine chemicals, nutraceuticals, and food supplements (Yang et al., 2013). The composition and contents of SMs are quality indicators that determine the therapeutic effects of medicinal plants (Yuan et al., 2020). Farmers pursue high yield and quality (measured by the levels of major SMs) when growing medicinal plants. Importantly, the accumulation profile of SMs in plants varies with specific local environment, in addition to plant genotype. For example, temperature and sunshine duration respectively affect the accumulation of polysaccharides and total alkaloids in the traditional Chinese orchid herb (Dendrobium officinale Kimura & Migo) (Yuan et al., 2020). Volatile organic compounds accumulate in lavender (Lavandula angustifolia Mill.) mainly depending on the latitudinal gradient (Demasi et al., 2018). In contrast, the SM profile of the bulk phloem of lodgepole pine (Pinus contorta Douglas ex Loudon) shifts with elevation (Mullin et al., 2021).
Schisandra sphenanthera Rehd. et Wils. is a perennial medicinal plant that is primarily found in western and southern China. It prefers shade and grows well under forest canopies. S. sphenanthera develops a solid lateral root and plays a positive role in soil and water conservation (Huang et al., 2021). Its dry ripe fruit, known as Schisandrae Sphenantherae Fructus (SSF), is one of the most famous traditional Chinese medicines and has been used to treat hepatitis, Alzheimer’s disease, osteoporosis, and insomnia (Li et al., 2022). SSF also has wide application in food and cosmetics (Tao et al., 2020). Lignans are the major bioactive components in SSF, and the content of schisantherin A is considered to be a standard indicator for quality evaluation (The State Pharmacopoeia Commission of China, 2020). Lignan components, such as schisantherin and schisanhenol, show pharmacological effects exemplified by hepatoprotective and antiviral activities (Yang et al., 2010; Sowndhararajan et al., 2018). Water-soluble polysaccharides from SSF also exhibit beneficial effects, including antioxidant, immunomodulatory, hepatoprotective, antitumor, and hypolipidemic activities (Zhao et al., 2013). Additionally, volatile oils are the major pharmacological components of SSF that could protect islet B cells and lower blood sugar (Song et al., 2007).
Nowadays, the herbal medicine SSF is mainly harvested from the wild resources of S. sphenanthera. In recent decades, it has been challenging to meet the growing market demand for SSF with wild resources only. The development of sustainable conservation and management strategies for S. sphenanthera calls for a holistic understanding of its potential geographical distribution and habitat preferences (Azareh et al., 2019). Moreover, climate change reportedly has significant implications for the habitat requirements of various species (Peng and Guo, 2017; Li et al., 2020a). Hence, ascertaining whether climate change affects the habitat suitability for S. sphenanthera is another critical issue related to its economic value and ecological significance (Zhang et al., 2020). So far, the ecological distribution of S. sphenanthera has been insufficiently investigated, and it is necessary to identify the priority areas where this medicinal plant should be cultivated under climate change.
In the present study, an optimized MaxEnt model was adopted to predict the distribution patterns of potentially suitable areas for S. sphenanthera under current and future climatic conditions. The purposes of the present study were to: (1) predict the potential geographical distribution of S. sphenanthera and identify its low impact areas under different climate scenarios and (2) determine the major environmental factors influencing the accumulation of medicinal components in SSF across regions. Results of this study could provide guidance for the exploitation and conservation of wild S. sphenanthera resources, as well as for the promotion of artificial cultivation over large scales.
2 Materials and methods
2.1 Prediction of potentially suitable habitats
2.1.1 Study area and data sources
The study area covered the territory of China with geographic coordinates between 73.1–135.5° E and 17.7–53.9° N. A total of 262 distribution records of S. sphenanthera were retrieved from publicly accessible databases, including the China Virtual Herbarium (CVH, https://www.cvh.ac.cn/; accessed on 08/12/2022) and the China National Knowledge Infrastructure (CNKI, https://www.cnki.net/; accessed on 08/12/2022). All the retrieved records were filtered to remove dubiously unreliable and redundant data using the method of Ye et al. (2020). To reduce errors caused by clustering effects, we retained only one distribution point in a 5×5 km grid using the ENMtools v1.4.1 software (Warren et al., 2010). A final set of 222 species occurrence points were used for modeling analysis (Figures S1 and S2, Table S1).
The raw environmental data used in this study were downloaded from network databases, such as CHELSA (https://chelsa-climate.org/), EarthEnv (http://www.earthenv.org/topography/), and ORNL DAAC (https://doi.org/10.3334/ORNLDAAC/1304/). The environmental variables included climate (Karger et al., 2017; Title and Bemmels, 2018), topography/landform (Amatulli et al., 2018), and soil conditions (Pelletier et al., 2016). Some variables were generated by R code provided in the original document.
A total of 47 environmental variables (Table S2) were considered for model prediction of the current (1981–2010) and future periods (2041–2070 and 2071–2100), and the same 47 variables were included in the prediction of each period. All the environmental data had a spatial resolution of 2.5 arc-minute. The geographical scope of the modeling analysis covered the territory of China. To limit the geographical scope, we used the Administrative Map of China and the Map of National Border–Territorial Sea as base maps for masking.
Considering the influence of climate scenario selection on model prediction, we selected the following five atmospheric general circulation models for climate data of future decades: GFDL-ESM4, IPSL-CM6A-LR, MPI-ESM1-2-HR, MRI-ESM2-0, and UKESM1-0-LL. Each model corresponded to a combination of three shared socioeconomic pathways (SSPs) and representative concentration pathways (RCPs): SSP1-RCP2.6, SSP3-RCP7, and SSP5-RCP8.5. The climate severity of the three climate scenarios increased from SSP1-RCP2.6 to SSP5-RCP8.5. To reduce the prediction uncertainty caused by a single atmospheric circulation model, the output layers corresponding to the three climate scenarios were arithmetically averaged. Therefore, 31 sets of climate data were used in this study, with a spatial resolution of 2.5 arc-minutes.
Significant environmental variables were selected by the MaxEnt v3.4.1 (Phillips et al., 2017) and ENMTools v1.4.1 software (Warren et al., 2010) to construct the MaxEnt model. Based on the data of species occurrence points and 47 environmental variables, we identified the most prominent bioclimatic factors that influenced the potential distribution of S. sphenanthera (Shcheglovitova and Anderson, 2013). First, the Jackknife method was employed to evaluate the importance of each variable. Second, Pearson correlation coefficients were used to evaluate the correlations among different variables. For each pair of significantly correlated variables (|r| ≥ 0.8), only one with a greater significant contribution to habitat suitability was retained (Wei et al., 2021; Yang et al., 2021).
2.1.2 Parameter optimization, model construction, and reliability test
The MaxEnt v3.4.1 software (Phillips et al., 2017) was used to construct a MaxEnt model for S. sphenanthera. To ensure that the potential distribution of S. sphenanthera was close to normal probability, 75% of the data were selected for model training, and the remaining data for model testing. The maximum number of iterations was 5000, where each process was repeated 20 times. Other parameter settings were set by default (Sun et al., 2020; Yan et al., 2020).
The ‘kuenm’ package in the R v3.4.1 software (Cobos et al., 2019) was used to optimize the feature class (FC) and regularization multiplier (RM) of the MaxEnt model. The RM was initially set to 0.1–4.0 at 0.1 intervals, resulting in 40 values. Then, four FCs (linear, L; quadratic, Q; hinge, H; and product, P) were arranged and combined to create 15 FC combinations (L, P, Q, H, LP, LQ, LH, LPQ, LPH, LQH, LPQH, PQH, PQ, PH, and QH). Thus, 600 parameter combinations were obtained with FC and RM. Based on the optimization results, we selected a statistically significant model that had an omission rate < 0.05 and a small-sample corrected delta Akaike’s information criterion ≤ 2 (Akpoti et al., 2020; Arenas-Castro et al., 2020).
The area under the receiver operating characteristic curve (AUC) was calculated to evaluate the accuracy of model prediction. The AUC range is (0,1), and a greater AUC value indicates a higher model reliability of distinguishing between suitable and unsuitable habitats. AUC ≥ 0.9 indicates high accuracy of model prediction (Guo et al., 2019; Liu et al., 2019).
2.1.3 Prediction of potentially suitable areas, low impact areas, and spatial pattern change
The habitat suitability for a given species is generally scored 0–1, and a higher suitability value indicates a better adaptive ability of the species. The maximum test sensitivity plus specificity (MTSPS) threshold is superior to other threshold options for the classification of potentially suitable areas (Tang et al., 2018). Therefore, we used the MTSPS threshold of 0.1463 in this study, and those areas with suitability values below the threshold were considered unsuitable for S. sphenanthera. The suitability values greater than the MTSPS threshold were equally divided into three groups, corresponding to poorly suitable areas (0.1463–0.4309), moderately suitable areas (0.4309–0.7154), and highly suitable areas (0.7154–1), respectively (Ye et al., 2018; Wei et al., 2021).
Low impact areas were defined as those areas where climate change had low impact on S. sphenanthera growth (Pan et al., 2020). These areas were identified by overlaying the binary prediction maps of potentially suitable areas in different periods and extracting the completely overlapped portions. First, the maps of potentially suitable areas in different periods were superimposed using the DIVA-GIS v7.5 software (http://www.diva-gis.org/). Then, the spatial units with distribution probability greater or lower than the MTSPS threshold were re-defined as suitable and unsuitable areas, respectively. This led to the establishment of the matrices of unsuitable and suitable areas, and the completely overlapped portions in the overlay layer were extracted as low impact areas. We predicted the low impact areas under three climate scenarios in the current (1981–2010) and future periods (2041–2070 and 2071–2100).
The spatial pattern change of potentially suitable areas across different periods was analyzed using the DIVA-GIS v7.5 software (Zurell et al., 2020; Santos-Hernandez et al., 2021; Wu et al., 2021). Based on the matrix table of suitable and unsuitable areas, we analyzed changes in the percentage of potentially suitable areas for S. sphenanthera under current and future climatic conditions.
2.2 Determination of secondary metabolites
2.2.1 Sample collection
A total of 32 fresh fruit samples were collected from different SSF-producing regions (cities or counties) throughout China in August 2020 (Table 1). The samples were oven-dried at 40°C to constant weight, crushed into powder, and passed through a 200-mesh sieve before analysis.
2.2.2 Chromatographic analysis
Lignans were extracted from SSF samples with methanol under ultrasonic irradiation and analyzed by ultra-high-performance liquid chromatography (Huang et al, 2008). We analyzed the extracts on an Ultimate UHPLC PolarRP C18 (100 mm × 2.1 mm, 1.8 mm) column connected to a Waters UPLCH-Class ultra-high-performance liquid chromatography system. There were six replicates for each sample. Nine standard substances (schisandrol A, B; schisantherin A, B; schisandrin A, B, C; gomisin J; and schisanhenol) were purchased from Solarbio Science & Technology Co. Ltd. (Beijing, China). Major peak identification, standard solution preparation, and methodology validation were conducted following the protocol of Wang et al. (2018).
Volatile oils were isolated by ethyl ether extraction and their chemical composition was analyzed by gas chromatography–mass spectrometry (Wang et al., 2017). Nine replicates were prepared for each sample.
Polysaccharide extraction was performed in ethyl alcohol, and the absorbance of sample extracts was measured by spectrophotometry at 490 nm with glucose as the substrate (Wang et al., 2011). The polysaccharide content in the samples was calculated as follows: polysaccharide content (%) = C*D*F/W × 100, where C is the concentration of glucose in the test solution (mg/mL), D is the dilution factor of the test solution, F is the conversion factor (F = w/c × d, where w is the polysaccharide content (mg), c is the glucose concentration in the polysaccharide diluent (mg/mL), and d is the dilution factor of polysaccharide), and W is the weight of the test sample (mg). Three replicates of each sample were tested.
2.3 Data analysis
All tests, including the prediction of potentially suitable habitats and the determination of SM contents, were carried out at least three times. The data were statistically analyzed using the R v4.0.5 software with different packages (yyeasy 0.0.4.0, https://gitee.com/yanpd01/yyeasy; psych 2.1.9, https://CRAN.R-project.org/package=psych; ggplot2 3.3.5, https://CRAN.R-project.org/package=ggplot2) (Robinson et al., 2010; Ginestet, 2011; Love et al., 2014). One-way analysis of variance was conducted to compare group means, and if differences existed, Tukey’s test was performed to assess their significance (p < 0.05). Pearson’s correlation analysis and mapping between the 47 environmental variables and the contents of major SMs in SSF were conducted. Pearson’s correlation analysis and mapping were performed among 47 environmental variables. Heatmaps of major SMs in SSF were created using the TBtools v1.113 software (Chen et al., 2020).
3 Results
3.1 Model parameter optimization and accuracy evaluation
Among the 47 environmental variables, eight were ultimately selected and used to construct the MaxEnt model (Table 2). These variables were average soil and sedimentary deposit thicknesses across upland hillslopes and valley bottoms (AvgSoilSedimDeposThick), isothermality (BIO03), annual range of temperature (BIO07), elevation, maximum temperature of the coldest month (maxTempColdest), net primary productivity (NPP), monthly variability in potential evapotranspiration (PETseasonality), and vector ruggedness measure (VRM). NPP (26.19%) was the largest contributor to habitat suitability, followed by BIO07 (18.22%), BIO03 (17.03%), and AvgSoilSedimDeposThick (13.66%). The relative contributions of PETseasonality (9.86%), elevation (8.33%), maxTempColdest (3.75%), and VRM (2.96%) were relatively low (Figures S3 and S4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1302417/fpls-14-1302417-HTML/image_m/fpls-14-1302417-t002.jpg)
Table 2 Environmental variables selected for model construction and their contributions to habitat suitability.
Based on the results of model optimization, the optimal FC and RM were LQ and 0.2, respectively. The AUCTRAIN and AUCTEST values of the optimized model were 0.9630 ± 0.0021 and 0.9615 ± 0.0085, respectively, with an |AUCDIFF| value of 0.0015. The partial AUC value was 0.9415 ± 0.0081 (Figure S5). All these estimates were indicative of high prediction accuracy.
Based on the MTSPS threshold (0.1463), the spatial units were subdivided as follows: 0-0.1463 Unsuitable areas; 0.1218-0.4309 Poorly suitable areas; 0.4309-0.7069 Moderately suitable areas; 0.7154-1 Highly suitable areas.
3.2 Projection of potentially suitable areas and their spatial pattern change
The potentially suitable areas for S. sphenanthera were 179.58×104 km2 in total under the current climatic conditions, accounting for 18.71% of China’s territory. These areas were primarily distributed across Shaanxi, Chongqing, Guizhou, Sichuan, and Hubei Provinces (or Municipality; Figure 1). Poorly, moderately, and highly suitable areas were 119.00×104 km2, 49.61×104 km2, and 10.98×104 km2, respectively. The highly suitable areas were mainly located in southern Shaanxi, northern Guizhou, southeastern Chongqing, western Hubei, and northern Hunan Provinces (or Municipality; Table 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1302417/fpls-14-1302417-HTML/image_m/fpls-14-1302417-g001.jpg)
Figure 1 Predicted distribution of S. sphenanthera in China under current climatic conditions (1981–2010). 0-0.1463 Unsuitable areas, represented by grey; 0.1218-0.4309 Poorly suitable areas, represented by yellow; 0.4309-0.7069 Moderately suitable areas, represented by orange; 0.7154-1 Highly suitable areas, represented by red.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1302417/fpls-14-1302417-HTML/image_m/fpls-14-1302417-t003.jpg)
Table 3 Predicted suitable areas for S. sphenanthera in China under current and future climatic conditions.
The potentially suitable areas for S. sphenanthera were predicted to decrease from the current to the future, and this trend would be more prominent with increasing severity of climate change (SSP1-RCP2.6 → SSP5-RCP8.5) (Figures 2, 3; Table 3). Specifically, the model prediction revealed a reduction in the suitable areas from 1981–2010 to 2071–2100 under the SSP1-RCP2.6 scenario, involving southern Shaanxi, northern Guizhou, southeastern Chongqing, and western Hubei Provinces (or Municipality). Under the SSP3-RCP7 scenario, the total suitable areas would shrink to varying degrees across the two future periods, accounting for 31.46% (2041–2070) and 28.51% (2071–2100) of the current suitable areas. Under the SSP5-RCP8.5 scenario, the future suitable areas would be remarkably smaller than the current suitable areas. The potentially suitable areas would shift westward from the current to the future, and the highly suitable areas would disappear completely by 2100 under different climate scenarios.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1302417/fpls-14-1302417-HTML/image_m/fpls-14-1302417-g002.jpg)
Figure 2 Predicted distribution of S. sphenanthera in China under future climatic scenarios (2041–2100).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1302417/fpls-14-1302417-HTML/image_m/fpls-14-1302417-g003.jpg)
Figure 3 Spatial change of potentially suitable areas for S. sphenanthera from the current to future periods.
3.3 Identification of low impact areas
The total low impact areas for S. sphenanthera would decrease from 39.57×104 to 30.49×104 km2 as a consequence of increased severity of climate change (SSP1-RCP2.6 → SSP5-RCP8.5). The corresponding proportions of low impact areas in the current suitable areas would also change from 22.03% to 16.98% (Table 4). However, specific areas represented by southern Shaanxi, northwestern Guizhou, southeastern Chongqing, and western Hubei Provinces (or Municipality) were identified as stable suitable areas despite different climate scenarios (Figure 4). Furthermore, the low impact areas under various combinations of climate scenarios were predicted. The largest low impact areas were obtained under the SSP1-RCP2.6 scenario, accounting for 22.03% of the current suitable areas. The low impact areas were smallest under the SSP1-RCP2.6 + SSP3-RCP7 + SSP5-RCP8.5 scenario combination, accounting for 16.51% of the current suitable areas (Figures 4 and S6).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1302417/fpls-14-1302417-HTML/image_m/fpls-14-1302417-g004.jpg)
Figure 4 Composite prediction of low impact areas for S. sphenanthera under various combinations of climate scenarios.
3.4 Accumulation patterns of major secondary metabolites across regions
The lignan content of SSF samples varied across different producing regions. Among the nine lignan components identified, schisantherin A was detected at the highest level in all regions, with gomisin J and schizandrol B at the lowest levels (Figure S7). The content of schisantherin A was exceptionally high at sites S25, S26, S28, and S32. The highest content of schisanhenol was observed at sites S13, S20, and S21. The volatile oil content of SSF samples also varied distinctly across regions and occurred at the highest levels at sites S31 and S32 (Figure S8). Although the composition of SSF was dominated by terpenoids (Table S3), there was considerable variation in the contents of volatile oil components (Figure S9). The polysaccharide content reached its highest level at sites S23 and S29, in contrast to the lowest levels at sites S1 and S2 (Figure S10). The accumulation of major SMs in SSF samples showed distinct patterns in various regions (Figure 5). The content of schisantherin A was highest at site S3, whereas the total content of lignans peaked at site S11. The total contents of both volatile oils and polysaccharides reached their maximum levels at site S23 (Table S4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1302417/fpls-14-1302417-HTML/image_m/fpls-14-1302417-g005.jpg)
Figure 5 Heatmap showing the contents of major secondary metabolites in Schisandrae Sphenantherae Fructus across 32 regions.
3.5 Relationship between metabolite accumulation and environmental factors
The content of schisantherin A in SSF showed high correlations with 14 environmental variables (Figure 6). For example, it was positively correlated with BIO07, temperature seasonality (BIO04), and mean temperature difference between the hottest and coldest months (continentality) (p < 0.01), and there was also a high positive correlation between these environmental factors (Figure S12). Conversely, negative correlation was observed between the content of schisantherin A and NPP, mean monthly precipitation of the wettest quarter (BIO16), and mean monthly precipitation of the warmest quarter (BIO18) (p < 0.05), and there is a high correlation between these environmental factors (Figure S12). The content of schisandrin A was positively correlated with BIO07 and BIO04 (p < 0.01). Furthermore, the content of total lignans was positively correlated with precipitation seasonality (BIO15) and elevation (p < 0.05), negatively correlated with mean daily temperature of the wettest quarter (BIO08) and mean monthly potential evapotranspiration of the wettest quarter (PETWettestQuarter) (p < 0.05). We also found that there was a high positive correlation between the environmental variables that had the same effect on the accumulation of the SMs of SSF, and the environment variables that had the opposite effect on the accumulation of SMs of SSF were negatively correlated (Figure S11). The correlations between major metabolite contents and the 47 environmental variables are shown in Figure S12.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1302417/fpls-14-1302417-HTML/image_m/fpls-14-1302417-g006.jpg)
Figure 6 Correlation heatmap between the contents of major secondary metabolites in Schisandrae Sphenantherae Fructus and a subset of environmental variables. * p < 0.05 and ** p < 0.01.
4 Discussion
4.1 Comprehensiveness and rationality of MaxEnt model prediction
MaxEnt modeling has become a popular tool in global change biology and biogeography, especially for predicting the potential distribution of animal and plant species (Araujo et al., 2005; Araujo and Luoto, 2007). Compared with other ENM models, the MaxEnt model uses relatively simple continuous or categorical variables and performs well with limited or incomplete datasets (Elith et al., 2011). However, previous studies have often neglected the optimization of model parameters, or the proposed models usually require further optimization, which may affect the prediction accuracy (Garcia-Valdes et al., 2013; Li et al., 2020a; Yan et al., 2020). Therefore, we developed an optimized MaxEnt model to predict potentially suitable areas for S. sphenanthera in China. The proposed model achieved high prediction accuracy based on a large AUC value.
In this study, a total of 47 environmental variables were considered in MaxEnt model prediction. The environmental dataset included land cover, soil regolith and sedimentary deposit thickness, cloud cover, topography, temperature, and precipitation. This dataset was more comprehensive than previously used in other studies. Many ENM-based studies have developed prediction models only using environmental data layers from WorldClim (Aidoo et al., 2023; Akyol et al., 2023). For example, Parveen et al. (2022) selected 11 temperature and eight precipitation-derived variables to build a MaxEnt model, and then used it to predict the distribution of potentially suitable habitats for the genus Nymphaea in India.
Furthermore, we took into consideration various factors affecting plant growth to optimize the MaxEnt model parameters systematically. These factors included topography, land cover, and cloud cover (Tuanmu and Jetz, 2014; Wilson and Jetz, 2016; Amatulli et al., 2018). The MaxEnt model was run for a maximum of 5000 iterations, with each process repeating 20 times, which could reduce prediction errors. Collectively, our analysis procedures ensured the accuracy, comprehensiveness, and rationality of the optimized MaxEnt model.
4.2 Current and future suitable areas for S. sphenanthera
Based on MaxEnt modeling, the current suitable areas for S. sphenanthera are mainly located in southern Shaanxi, northern Guizhou, southeastern Chongqing, western Hubei, and northern Hunan Provinces, in agreement with previous studies (Heyu et al., 2016). Temperature, precipitation, and soil conditions are major environmental factors driving medicinal plant distribution (Li et al., 2020b). However, the primary factors that determine the range of potentially suitable areas for diverse medicinal plants are variable. For example, Hu et al. (2012) identified mean annual and extreme temperature (minimum and maximum) as the most important limiting factors for Schisandra chinensis (Turcz.) Baill during the critical growth period. This is inconsistent with our finding that the growth of S. sphenanthera is principally influenced by the carbon sequestration capacity of plants (NPP), temperature-dependent climatic factors (BIO07 and BIO03), and soil physical quality (AvgSoilSedimDeposThick).
NPP is primarily determined by photosynthetic energy and mainly constrained by light, water, and nutrient resources, which are necessary for plant yield formation (Briggs et al., 2018). Temperature and soil physical conditions are also essential factors for plant growth. Field surveys and research reports indicate that S. sphenanthera prefers to grow in warm and humid climate zones with a long, sunny, and warm summer (Zhang et al., 2014). Specifically, the mean annual temperature of 15–20°C, the mean annual precipitation of 1,000–2,000 mm, and the soil thickness of 30 cm are favorable conditions for S. sphenanthera growth (Lu and Chen, 2009). These exciting findings corroborate the accuracy of the MaxEnt model proposed in this study.
There will be a remarkable trend of decreasing suitable areas for S. sphenanthera in the future decades. Predictions under the SSP5-RCP8.5 scenario produced the worst outcome, followed by SSP3-RCP7 and SSP1-RCP2.6. This means that the wild resources of S. sphenanthera in China will be depleted in the 21st century as a consequence of decreased growing areas under possible climate change. Under all three climate scenarios, the potentially suitable areas will shift westward from the current to the future and eventually disappear in specific regions of southeastern China. Various plant species respond differently to climate change in terms of geographical distribution. For example, the suitable habitat distribution of Paeonia veitchii was predicted to shift toward higher elevations and latitudes under the RCP8.5 scenario (Zhang et al., 2019a). A northward shift and an eastward expansion of planting boundaries were predicted for rice and maize cultivars under the SRES A1B and A2 scenarios similar to RCP8.5 (Yuan et al., 2012; Zhang et al., 2017).
4.3 Influence of environmental factors on medicinal plant quality
As expected, the habitat environment had a profound influence on the accumulation of SMs in SSF. This was demonstrated by the correlations between the contents of various SMs in 32 samples of SSF collected so far and different environmental factors. Temperature- and precipitation-dependent climatic factors were the major drivers of SM accumulation patterns in SSF. Specifically, the content of schisantherin A exhibited an upward trend in response to increasing annual temperature extremes (i.e., difference between the annual maximum and minimum temperatures) or seasonal temperature variation. The opposite trend was observed for schisantherin A with increasing annual or seasonal precipitation. Additionally, the total content of lignans exhibited an upward trend with increasing seasonal precipitation and elevation, whereas a downward trend emerged with increasing temperature and evapotranspiration of the wettest quarter. All these findings lead us to posit that specific environmental stress could enhance the accumulation of major medicinal components in SSF. A greater change in temperature or precipitation is likely to enable S. sphenanthera plants to accumulate more lignans (especially schisantherin A) in the fruit. In the analysis of environmental factors’ autocorrelation, it was found that there is correlation among the environmental factors that affect the content of SMs in SSF. For example, elevation and BIO15 (precipitation seasonality) were positively correlated with the total lignan content of SSF, and there was also a high positive correlation between these two environmental factors. On the other hand, environmental factors such as BIO08 (mean Temperature of wettest quarter) and BIO05 (max temperature of warmest month) were negatively correlated with the total lignan content. However, these environmental factors show a negative correlation with elevation and BIO15 (Figures 6 and S11). Such correlation aligns with the fact that as altitude increases, temperature gradually decreases (Kad et al., 2022). Therefore, when selecting planting areas for high-quality S. sphenanthera, it is possible to summarize these environmental factors that affect the content of one medicinal component (such as total lignans) and make a comprehensive decision on the planting areas for S. sphenanthera. Within the suitable habitat range, S. sphenanthera plants would produce SSF with higher levels of schisantherin A when growing in areas with higher temperature, less intense light, and lower precipitation. More lignans would accumulate in SSF when planting S. sphenanthera in low-temperature areas.
In addition to their medicinal functions, SMs play an essential role in plant adaptation to abiotic stresses, including temperature, drought, nutrient deficiency, carbon dioxide elevation, salinity, and ultraviolet light (Renault et al., 2010; Ramakrishna and Ravishankar, 2011; Rao et al., 2021). For example, terpenoids accumulate in plants in response to abiotic stresses, providing antioxidant protection against drought, temperature, light, and salt stresses (Zhang et al., 2019b). Additionally, nitrogen-containing compounds, such as alkaloids and glucosinolates, occur at high levels in plants with stress responses, which can alleviate oxidative stress (Zhu et al., 2015). Our study also found that the accumulation of some SMs in SSF was stimulated under certain environmental stress. Therefore, the selection of planting areas for S. sphenanthera should take into account not only the habitat suitability, but also the accumulation of medicinal components. In this study, samples were collected from 32 S. sphenanthera growing areas. Results indicate that the content of lignans in Huxian County, Shaanxi Province, and Pingwu County, Sichuan Province, were relatively high, while polysaccharides and volatile oil content were higher in Zhashui County, Shaanxi Province. These locations are also considered moderately to highly suitable for the growth of S. sphenanthera. Therefore, priority could be given to the cultivation of S. sphenanthera in these areas.
Like environmental factors, human activities play a non-negligible role in shaping the geographical distribution of S. sphenanthera (Zhang et al., 2014). In fact, no single environmental and human factors are expected to independently control the quality of medicinal plants. A combination of multiple factors can provide a better explanation for the variation in medicinal plant quality (Shen et al., 2010). This underscores the necessity of taking into account both environmental factors (e.g., climate, soil, terrain) and human activities (e.g., fertilization, irrigation) when planting S. sphenanthera for high-yield and high-quality production of SSF.
5 Conclusions
This study developed an optimized maximum entropy model and successfully used it to predict the distribution patterns and spatial changes of potentially suitable areas for S. sphenanthera in China. We found that the potentially suitable areas would shrink from the current (1981–2010) to the future (2041–2100), depending on the severity of climate change. The areas of low climate impact were mainly found in southern Shaanxi, northern Guizhou, northwestern Guizhou, southeast Chongqing, and western Hubei Provinces. These low impact areas should be prioritized to establish nature reserves for the conservation of wild S. sphenanthera resources or develop farmland for the cultivation of this medicinal plant. The contents of secondary metabolites in Schisandrae Sphenantherae Fructus varied across regions as a consequence of different climatic environments. Diurnal temperature variation, together with seasonal and annual precipitation, prominently influenced the accumulation of secondary metabolites as medicinal components and stress response regulators. Our findings indicate that specific environmental stress is conducive to the improvement of medicinal plant quality. This research can provide scientific guidance for the conservation planning of wild plant resources, large-scale cultivation of S. sphenanthera, and quality improvement of finished medicinal materials.
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 authors.
Author contributions
JS: Formal Analysis, Methodology, Visualization, Writing – original draft, Writing – review & editing. QZ: Methodology, Software, Visualization, Writing – review & editing. PY: Methodology, Software, Visualization, Writing – review & editing. MS: Software, Visualization, Writing – review & editing. HS: Project administration, Validation, Writing – review & editing. HL: Conceptualization, Resources, Writing – review & editing. DZ: Funding acquisition, Investigation, Resources, Writing – review & editing. ZQ: Conceptualization, Data curation, Writing – review & editing. LC: Conceptualization, Data curation, Resources, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Key Research and Development Program (2022YFD1602000, 2019YFC1712600), the Central Guidance on Local Science and Technology Development Fund (2021ZY2 -CG-02, 2022ZY1-CGZY-02, 2022ZY1-CGZY-03), the Technology Innovation Leading Program of Shaanxi (2023QYPY2-02).
Conflict of interest
Author DZ is employed by Shaanxi Panlong Pharmaceutical Group Limited by Share Ltd.
The remaining 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.1302417/full#supplementary-material
References
Aidoo, O. F., Souza, P. G. C., Silva, R. S., Junior, P. A. S., Picanco, M. C., Heve, W. K., et al. (2023). Modeling climate change impacts on potential global distribution of Tamarixia radiata Waterston (Hymenoptera: Eulophidae). Sci. Total Environ. 864, 160962. doi: 10.1016/j.scitotenv.2022.160962
Akpoti, K., Kabo-bah, A. T., Dossou-Yovo, E. R., Groen, T. A., Zwart, S. J. (2020). Mapping suitability for rice production in inland valley landscapes in Benin and Togo using environmental niche modeling. Sci. Total Environ. 709, 136165. doi: 10.1016/j.scitotenv.2019.136165
Akyol, A., Orucu, O. K., Arslan, E. S., Sarikaya, A. G. (2023). Predicting of the current and future geographical distribution of Laurus nobilis L. under the effects of climate change. Environ. Monit. Assess. 195, 459. doi: 10.1007/s10661-023-11086-z
Amatulli, G., Domisch, S., Tuanmu, M.-N., Parmentier, B., Ranipeta, A., Malczyk, J., et al. (2018). A suite of global, cross-scale topographic variables for environmental and biodiversity modeling. Sci. Data 5 (1), 180040. doi: 10.1038/sdata.2018.40
Araujo, M. B., Luoto, M. (2007). The importance of biotic interactions for modelling species distributions under climate change. Global Ecol. Biogeogr. 16, 743–753. doi: 10.1111/j.1466-8238.2007.00359.x
Araujo, M. B., Pearson, R. G., Thuiller, W., Erhard, M. (2005). Validation of species-climate impact models under climate change. Global Change Biol. 11, 1504–1513. doi: 10.1111/j.1365-2486.2005.01000.x
Arenas-Castro, S., Goncalves, J. F., Moreno, M., Villar, R. (2020). Projected climate changes are expected to decrease the suitability and production of olive varieties in southern Spain. Sci. Total Environ. 709, 20. doi: 10.1016/j.scitotenv.2019.136161
Azareh, A., Rahmati, O., Rafiei-Sardooi, E., Sankey, J. B., Lee, S., Shahabi, H., et al. (2019). Modelling gully-erosion susceptibility in a semi-arid region, Iran: Investigation of applicability of certainty factor and maximum entropy models. Sci. Total Environ. 655, 684–696. doi: 10.1016/j.scitotenv.2018.11.235
Blank, L., Blaustein, L. (2012). Using ecological niche modeling to predict the distributions of two endangered amphibian species in aquatic breeding sites. Hydrobiologia 693, 157–167. doi: 10.1007/s10750-012-1101-5
Briggs, N., Gudmundsson, K., Cetinic, I., D’Asaro, E., Rehm, E., Lee, C., et al. (2018). A multi-method autonomous assessment of primary productivity and export efficiency in the springtime North Atlantic. Biogeosciences 15, 4515–4532. doi: 10.5194/bg-15-4515-2018
Chen, C. J., Chen, H., Zhang, Y., Thomas, H. R., Frank, M. H., He, Y. H., et al. (2020). TBtools: an integrative toolkit developed for interactive analyses of big biological data. J. Citation Rep. 13, 1194–1202. doi: 10.1016/j.molp.2020.06.009
Chen, Y., Lu, S., Mao, L. (2021). Prediction of future changes in suitable distribution area for rare tree species of Dalbergia. J. Zhejiang Univ. 38, 837–845. doi: 10.11833/j.issn.2095-0756.20200522
Cobos, M. E., Peterson, A. T., Barve, N., Osorio-Olvera, L. (2019). Kuenm: an R package for detailed development of ecological niche models using Maxent. Peerj 7, e6281. doi: 10.7717/peerj.6281
Demasi, S., Caser, M., Lonati, M., Cioni, P. L., Pistelli, L., Najar, B., et al. (2018). Latitude and altitude influence secondary metabolite production in peripheral alpine populations of the mediterranean species Lavandula angustifolia Mill. Front. Plant Sci. 9. doi: 10.3389/fpls.2018.00983
Elith, J., Phillips, S. J., Hastie, T., Dudík, M., Chee, Y. E., Yates, C. J. (2011). A statistical explanation of MaxEnt for ecologists. Diversity Distributions 17, 43–57. doi: 10.1111/j.1472-4642.2010.00725.x
Feng, G., Xiong, Y. J., Wei, H. Y., Li, Y., Mao, L. F. (2023). Endemic medicinal plant distribution correlated with stable climate, precipitation, and cultural diversity. Plant Diversity 45 (4), 479–484. doi: 10.1016/j.pld.2022.09.007
Garcia-Valdes, R., Zavala, M. A., Araujo, M. B., Purves, D. W. (2013). Chasing a moving target: projecting climate change-induced shifts in non-equilibrial tree species distributions. J. Ecol. 101, 441–453. doi: 10.1111/1365-2745.12049
Ginestet, C. (2011). ggplot2: elegant graphics for data analysis. J. R. Stat. Soc. Ser. a-Statistics Soc. 174, 245–245. doi: 10.1111/j.1467-985X.2010.00676_9.x
Graham, C. H., Elith, J., Hijmans, R. J., Guisan, A., Peterson, A. T., Loiselle, B. A., et al. (2008). The influence of spatial errors in species occurrence data used in distribution models. J. Appl. Ecol. 45 (1), 239–247. doi: 10.1111/j.1365-2664.2007.01408.x
Guo, Y. L., Li, X., Zhao, Z. F., Nawaz, Z. (2019). Predicting the impacts of climate change, soils and vegetation types on the geographic distribution of Polyporus umbellatus in China. Sci. Total Environ. 648, 1–11. doi: 10.1016/j.scitotenv.2018.07.465
Han, F., Xiao, Z., Li, Q., Liu, J., Han, R., Li, J., et al. (2020). Constrastive study on main genetic characters of different types of Atractylodes lancea. Chin. Wild Plant Resour. 39, 1–8. doi: 10.3969/j.issn.1006-9690.2020.05.001
Heyu, Y., Haiyan, W., Manjie, S., Zhonghui, S., Yajuan, M. (2016). Phenotypic plasticity of Schisandra sphenanthera leaf and the effect of environmental factors on leaf phenotype. Chin. Bull. Bot. 51 (3), 322–334.
Hu, L. L., Zhang, H. Y., Qin, L., Yan, B. Q. (2012). Current distribution of Schisandra chinensis in China and its predicted responses to climate change. J. Appl. Ecol. 23, 2445–2450. doi: 10.13287/j.1001-9332.2012.0338
Huang, S. Q., Zhang, D. D., Li, Y. Z., Fan, H., Liu, Y. Y., Huang, W. L., et al. (2021). Schisandra sphenanthera: a comprehensive review of its botany, phytochemistry, pharmacology, and clinical applications. Am. J. Chin. Med. 49, 1577–1622. doi: 10.1142/s0192415x21500749
Huang, X., Song, F. R., Liu, Z. Q., Liu, S. Y. (2008). Studies on the lignans in extract of the fruits of Schisandra chinensis and Schisandra sphenanthera by high performance liquid chromatography-electrospray ionization mass spectrometry. Acta Chimica Sinica 66 (9), 1059–1066.
Ji, B. Y., Liu, M., Pei, L. X., Yang, L. L. (2021). Ecologically suitable areas for growing Scutellaria baicalensis worldwide: an analysis based on GMPGIS. China J. Chin. Materia Med. 46, 4389–4394. doi: 10.19540/j.cnki.cjcmm.20210625.101
Kad, P., Blau, M. T., Ha, K. J., Zhu, J. (2022). Elevation-dependent temperature response in early Eocene using paleoclimate model experiment. Environ. Res. Lett. 17, (11). doi: 10.1088/1748-9326/ac9c74
Kaky, E., Nolan, V., Alatawi, A., Gilbert, F. (2020). A comparison between Ensemble and MaxEnt species distribution modelling approaches for conservation: a case study with Egyptian medicinal plants. Ecol. Inf. 60. doi: 10.1016/j.ecoinf.2020.101150
Karger, D. N., Conrad, O., Böhner, J., Kawohl, T., Kreft, H., Soria-Auza, R. W., et al. (2017). Climatologies at high resolution for the earth’s land surface areas. Sci. Data 4 (1), 170122. doi: 10.1038/sdata.2017.122
Li, J. J., Fan, G., He, Y. (2020a). 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, Y. Q., Kong, D. X., Fu, Y., Sussman, M. R., Wu, H. (2020b). The effect of developmental and environmental factors on secondary metabolites in medicinal plants. Plant Physiol. Biochem. 148, 80–89. doi: 10.1016/j.plaphy.2020.01.006
Li, F., Li, B., Liu, J., Wei, X., Qiang, T., Mu, X., et al. (2022). Anti-asthmatic fraction screening and mechanisms prediction of Schisandrae Sphenantherae Fructus based on a combined approach. Front. Pharmacol. 13. doi: 10.3389/fphar.2022.902324
Li, G. L., Yang, J., Duan, J. A., Liu, H. B., Zhu, Z. H., Qian, D. W., et al. (2015). Quality analysis and evaluation of Anemarrhena asphodeloides rhizome from different habitats. J. Chin. medicinal materials 38 (6), 1148–1152. doi: 10.13863/j.issn1001-4454.2015.06.009
Liu, B. Y., Gao, X., Ma, J., Jiao, Z. H., Xiao, J. H., Hayat, M. A., et al. (2019). Modeling the present and future distribution of arbovirus vectors Aedes aEgypti and Aedes albopictus under climate change scenarios in mainland China. Sci. Total Environ. 664, 203–214. doi: 10.1016/j.scitotenv.2019.01.301
Love, M. I., Huber, W., Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550. doi: 10.1186/s13059-014-0550-8
Lu, Y., Chen, D. F. (2009). Analysis of Schisandra chinensis and Schisandra sphenanthera. J. Chromatogr. A 1216, 1980–1990. doi: 10.1016/j.chroma.2008.09.070
Luo, Y., Dong, Y. B., Zhu, C., Peng, W. F., Fang, Q. M., Xu, X. L. (2017). Research on suitable distribution of Paris yunnanensis based on remote sensing and GIS. China J. Chin. Materia Med. 42, 4378–4386. doi: 10.19540/j.cnki.cjcmm.2017.0188
Mullin, M., Klutsch, J. G., Cale, J. A., Hussain, A., Zhao, S., Whitehouse, C., et al. (2021). Primary and secondary metabolite profiles of lodgepole pine trees change with elevation, but not with latitude. J. Chem. Ecol. 47, 280–293. doi: 10.1007/s10886-021-01249-y
Pan, J. W., Fan, X., Luo, S. Q., Zhang, Y. Q., Yao, S., Guo, Q. Q., et al. (2020). Predicting the potential distribution of two varieties of Litsea coreana (Leopard-Skin Camphor) in China under climate change. Forests 11, 1159. doi: 10.3390/f11111159
Parveen, S., Kaur, S., Baishya, R., Goel, S. (2022). Predicting the potential suitable habitats of genus Nymphaea in India using MaxEnt modeling. Environ. Monit. Assess. 194, 853. doi: 10.1007/s10661-022-10524-8
Pelletier, J. D., Broxton, P. D., Hazenberg, P., Zeng, X., Troch, P. A., Niu, G., et al. (2016). Global 1-km gridded thickness of soil, regolith, and sedimentary deposit layers (Oak Ridge, Tennessee, USA: ORNL Distributed Active Archive Center) 2016. doi: 10.3334/ORNLDAAC/1304
Peng, L., Guo, Y. (2017). Geographical distribution of Astragali Radix and prediction of its suitable srea in China. J. Sichuan Agric. Univ. 35 (1), 60–68. doi: 10.16036/j.issn.1000-2650.2017.01.009
Phillips, S. J., Anderson, R. P., Dudik, M., Schapire, R. E., Blair, M. E. (2017). Opening the black box: an open-source release of Maxent. Ecography 40, 887–893. doi: 10.1111/ecog.03049
Ramakrishna, A., Ravishankar, G. A. (2011). Influence of abiotic stress signals on secondary metabolites in plants. Plant Signaling Behav. 6, 1720–1731. doi: 10.4161/psb.6.11.17613
Rao, S. Q., Chen, X. Q., Wang, K. H., Zhu, Z. J., Yang, J., Zhu, B. (2021). Effect of short-term high temperature on the accumulation of glucosinolates in Brassica rapa. Plant Physiol. Biochem. 161, 222–233. doi: 10.1016/j.plaphy.2021.02.013
Renault, H., Roussel, V., El Amrani, A., Arzel, M., Renault, D., Bouchereau, A., et al. (2010). The Arabidopsis pop2-1 mutant reveals the involvement of GABA transaminase in salt stress tolerance. BMC Plant Biol. 10. doi: 10.1186/1471-2229-10-20
Robinson, M. D., McCarthy, D. J., Smyth, G. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616
Santos-Hernandez, A. F., Monterroso-Rivas, A. I., Granados-Sanchez, D., Villanueva-Morales, A., Santacruz-Carrillo, M. (2021). Projections for Mexico’s tropical rainforests considering ecological niche and climate change. Forests 12, 119. doi: 10.3390/f12020119
Shan, Z. J., Ye, J. F., Hao, D. C., Xiao, P. G., Chen, Z. D., Lu, A. M. (2022). Distribution patterns and industry planning of commonly used traditional Chinese medicinal plants in China. Plant Diversity 44 (3), 255–261. doi: 10.1016/j.pld.2021.11.003
Shcheglovitova, M., Anderson, R. P. (2013). Estimating optimal complexity for ecological niche models: a jackknife approach for species with small sample sizes. Ecol. Model. 269, 9–17. doi: 10.1016/j.ecolmodel.2013.08.011
Shen, W., Ye, Z., Guo, Q. (2010). Quality evaluation of Schisandra chinensis at different producing areas in northeast China. China J. Chin. Materia Med. 35 (22), 3016–3020.
Song, L., Ding, J. Y., Tang, C., Yin, C. H. (2007). Compositions and biological activities of essential oils of Kadsura longepedunculata and Schisandra sphenanthera. Am. J. Chin. Med. 35, 353–364. doi: 10.1142/s0192415x07004874
Sowndhararajan, K., Deepa, P., Kim, M., Park, S. J., Kim, S. (2018). An overview of neuroprotective and cognitive enhancement properties of lignans from Schisandra chinensis. Biomed. Pharmacother. 97, 958–968. doi: 10.1016/j.biopha.2017.10.145
Sun, S. X., Zhang, Y., Huang, D. Z., Wang, H., Cao, Q., Fan, P. X., et al. (2020). The effect of climate change on the richness distribution pattern of oaks (Quercus L.) in China. Sci. Total Environ. 744, 140786. doi: 10.1016/j.scitotenv.2020.140786
Tang, C. Q., Matsui, T., Ohashi, H., Dong, Y. F., Momohara, A., Herrando-Moraira, S., et al. (2018). Identifying long-term stable refugia for relict plant species in east Asia. Nat. Commun. 9, 4488. doi: 10.1038/s41467-018-07727-4
Tao, X., Sun, Y., Men, X., Xu, Z. (2020). A compound plant extract and its antibacterial and antioxidant properties in vitro and in vivo. 3 Biotech. 10, 532. doi: 10.1007/s13205-020-02529-2
The State Pharmacopoeia Commission of China (2020). Pharmacopoeia of the people’s republic of China, Part I Vol. 255 (Beijing: Chemical Industry Press).
Title, P. O., Bemmels, J. B. (2018). ENVIREM: an expanded set of bioclimatic and topographic variables increases flexibility and improves performance of ecological niche modeling. Ecography 41, 291–307. doi: 10.1111/ecog.02880
Tuanmu, M. N., Jetz, W. (2014). A global 1-km consensus land-cover product for biodiversity and ecosystem modelling. Global Ecol. Biogeogr. 23 (9), 1031–1045. doi: 10.1111/geb.12182
Wang, Z. S., Chen, H. X., Zhang, W. J., Lan, G. S., Zhang, L. K. (2011). Comparative studies on the chemical composition and antioxidant activities of Schisandra chinensis and Schisandra sphenanthera fruits. J. Medicinal Plants Res. 5, 1207–1216. doi: 10.5897/JMPR.9000379
Wang, X. R., Liu, Y., Zhou, S. C., Qin, X. L., Gu, W. (2017). Comparison on the compositions of essential oils from the seed and pulp of Schisandra sphenanthera Rehd. et Wils. J. Essential Oil Bearing Plants 20, 1066–1073. doi: 10.1080/0972060x.2017.1375866
Wang, X. T., Yu, J. H., Li, W., Wang, C. M., Li, H., Ju, W. B., et al. (2018). Characteristics and antioxidant activity of lignans in Schisandra chinensis and Schisandra sphenanthera from different locations. Chem. Biodiversity 15, e1800030. doi: 10.1002/cbdv.201800030
Warren, D. L., Glor, R. E., Turelli, M.. (2010). ENMTools: a toolbox for comparative studies of environmental niche models. Ecography 33 (3), 607–611. doi: 10.1111/j.1600-0587.2009.06142.x
Wei, Y., Zhang, L., Wang, J., Wang, W., Niyati, N., Guo, Y., et al. (2021). Chinese caterpillar fungus (Ophiocordyceps sinensis) in China: current distribution, trading, and futures under climate change and overexploitation. Sci. Total Environ. 755, 142548. doi: 10.1016/j.scitotenv.2020.142548
Wilson, A. M., Jetz, W. (2016). Remotely sensed high-resolution global cloud dynamics for predicting ecosystem and biodiversity distributions. PloS Biol. 14, (3). doi: 10.1371/journal.pbio.1002415
Wu, Y. M., Shen, X. L., Tong, L., Lei, F. W., Mu, X. Y., Zhang, Z. X. (2021). Impact of past and future climate change on the potential distribution of an endangered montane shrub Lonicera oblata and its conservation implications. Forests 12, 125. doi: 10.3390/f12020125
Xu, D. P., Zhuo, Z. H., Wang, R. L., Ye, M., Pu, B. (2019). Modeling the distribution of Zanthoxylum armatum in China with MaxEnt modeling. Global Ecol. Conserv. 19, e00691. doi: 10.1016/j.gecco.2019.e00691
Yan, H. Y., Feng, L., Zhao, Y. F., Feng, L., Zhu, C. P., Qu, Y. F., et al. (2020). Predicting the potential distribution of an invasive species, Erigeron canadensis L., in China with a maximum entropy model. Global Ecol. Conserv. 21, e00822. doi: 10.1016/j.gecco.2019.e00822
Yang, Z., Bai, Y., Alatalo, J. M., Huang, Z., Yang, F., Pu, X., et al. (2021). Spatio-temporal variation in potential habitats for rare and endangered plants and habitat conservation based on the maximum entropy model. Sci. Total Environ. 784, 147080. doi: 10.1016/j.scitotenv.2021.147080
Yang, G. Y., Fan, P., Wang, R. R., Cao, J. L., Xiao, W. L., Yang, L. M., et al. (2010). Dibenzocyclooctadiene lignans from Schisandra lancifolia and their anti-human immunodeficiency virus-1 activities. Chem. Pharm. Bull. 58, 734–737. doi: 10.1248/cpb.58.734
Yang, X. Q., Kushwaha, S. P. S., Saran, S., Xu, J., Roy, P. S. (2013). Maxent modeling for predicting the potential distribution of medicinal plant, Justicia adhatoda L. @ in Lesser Himalayan foothills. Ecol. Eng. 51, 83–87. doi: 10.1016/j.ecoleng.2012.12.004
Ye, X. P., Yu, X. P., Yu, C. Q., Tayibazhaer, A., Xu, F. J., Skidmore, A. K., et al. (2018). Impacts of future climate and land cover changes on threatened mammals in the semi-arid Chinese Altai Mountains. Sci. Total Environ. 612, 775–787. doi: 10.1016/j.scitotenv.2017.08.191
Ye, X. Z., Zhao, G. H., Zhang, M. Z., Cui, X. Y., Fan, H. H., Liu, B. (2020). Distribution pattern of endangered plant Semiliquidambar cathayensis (Hamamelidaceae) in response to climate change after the last interglacial period. Forests 11, 434. doi: 10.3390/f11040434
Yuan, B., Guo, J. P., Ye, M. Z., Zhao, J. F. (2012). Variety distribution pattern and climatic potential productivity of spring maize in northeast China under climate change. Chin. Sci. Bull. 57, 3497–3508. doi: 10.1007/s11434-012-5135-x
Yuan, Y., Tang, X., Jia, Z., Li, C., Ma, J., Zhang, J. (2020). The effects of ecological factors on the main medicinal components of Dendrobium officinale under different cultivation modes. Forests 11, 94. doi: 10.3390/f11010094
Zhang, X. H., Niu, M. Y., da Silva, J. A. T., Zhang, Y. Y., Yuan, Y. F., Jia, Y. X., et al. (2019b). Identification and functional characterization of three new terpene synthase genes involved in chemical defense and abiotic stresses in Santalum album. BMC Plant Biol. 19. doi: 10.1186/s12870-019-1720-3
Zhang, K. L., Sun, L. P., Tao, J. (2020). Impact of climate change on the distribution of Euscaphis japonica(Staphyleaceae) Trees. Forests 11, (5). doi: 10.3390/f11050525
Zhang, Y., Wang, Y., Niu, H. (2017). Spatio-temporal variations in the areas suitable for the cultivation of rice and maize in China under future climate scenarios. Sci. Total Environ. 601-602, 518–531. doi: 10.1016/j.scitotenv.2017.05.232
Zhang, W., Xu, M., Liu, X., Gu, W. (2014). Investigation of wild plant resources of Schisandra in Qinling Mountains. J. Plant Genet. Resour. 15 (2), 236–241. doi: 10.13430/j.cnki.jpgr.2014.02.002
Zhang, K. L., Zhang, Y., Tao, J. (2019a). Predicting the potential distribution of Paeonia veitchii (Paeoniaceae) in China by incorporating climate change into a maxent model. Forests 10, 190. doi: 10.3390/f10020190
Zhao, T., Mao, G. H., Zhang, M., Li, F., Zou, Y., Zhou, Y., et al. (2013). Anti-diabetic effects of polysaccharides from ethanol-insoluble residue of Schisandra chinensis (Turcz.) Baill on alloxan-induced diabetic mice. Chem. Res. Chin. Universities 29, 99–102. doi: 10.1007/s40242-012-2218-9
Zhou, J., Guo, L., Zhang, J., Yang, G., Zhao, M., Huang, L. (2010). Responses of medicinal plant to drought stress and controlled experiment. China J. Chin. Materia Med. 35 (15), 1919–1924.
Zhu, W., Yang, B. X., Komatsu, S., Lu, X. P., Li, X. M., Tian, J. K. (2015). Binary stress induces an increase in indole alkaloid biosynthesis in Catharanthus roseus. Front. Plant Sci. 6. doi: 10.3389/fpls.2015.00582
Keywords: medicinal plant, climate change, maximum entropy modeling, suitable habitats, secondary metabolites
Citation: Shang J, Zhao Q, Yan P, Sun M, Sun H, Liang H, Zhang D, Qian Z and Cui L (2023) Environmental factors influencing potential distribution of Schisandra sphenanthera and its accumulation of medicinal components. Front. Plant Sci. 14:1302417. doi: 10.3389/fpls.2023.1302417
Received: 03 October 2023; Accepted: 28 November 2023;
Published: 12 December 2023.
Edited by:
Zhentian Lei, University of Missouri, United StatesCopyright © 2023 Shang, Zhao, Yan, Sun, Sun, Liang, Zhang, Qian and Cui. 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: Zengqiang Qian, zengqiang.qian@yahoo.com; Langjun Cui, ljcui@snnu.edu.cn
†These authors have contributed equally to this work and share first authorship