
==== Front
BMC Plant Biol
BMC Plant Biol
BMC Plant Biology
1471-2229
BioMed Central London

5518
10.1186/s12870-024-05518-2
Research
Temperature seasonality and soil phosphorus availability shape ginseng quality via regulating ginsenoside contents
Wu Dehua 12
Xiong Feng 2
Wang Hongyang 2
Liu Siqi 2
Zhu Jitong 2
Zhao Dan 3
Yang Jian 2
Ma Wenqi 2
http://orcid.org/0000-0001-6974-1411
Guo Lanping glp01@126.com

2
Kang Chuanzhi kangchuanzhi1103@163.com

2
1 https://ror.org/02vg7mz57 grid.411847.f 0000 0004 1804 4300 Guangdong Pharmaceutical University, Guangzhou, 510006 China
2 https://ror.org/042pgcv68 grid.410318.f 0000 0004 0632 3409 State Key Laboratory for Quality Ensurance and Sustainable Use of Dao-di Herbs, National Resource Center for Chinese Materia Medica, China Academy of Chinese Medical Sciences, Beijing, 100700 PR China
3 grid.443382.a 0000 0004 1804 268X Guizhou University of Traditional Chinese Medicine, Guiyang, 550025 China
3 9 2024
3 9 2024
2024
24 82416 5 2024
14 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
The accumulation of secondary metabolites in Panax ginseng Meyer (P. ginseng) exhibits significant geographical variation, normally due to environmental factors. The current study aimed at elucidating the key environmental factors modulating the accumulation of secondary metabolites in P. ginseng. Plant and the associated soil samples were collected from ten geographical locations within the latitudinalrange of 27.09°N − 42.39°N and longitudinal range of 99.28°E − 128.19°E. 12 secondary metabolites in P. ginseng toots were measured. And the correlation between secondary metabolites with a series of soil properties and 7 climatic factors were investigated through Pearson’s correlation, mantel test, random forest and pathway analysis. The results revealed that climatic factors were stronger drivers of ginseng secondary metabolite profile than soil nutrients. Specifically, temperature seasonality (TS) and soil available phosphorus (AP) were the most effective environments to have significantly and positively influence on the secondary metabolites of ginseng. This findings contribute to identifying optimal cultivation areas for P. ginseng, and hopefully establishing methods for interfering/shaping microclimate for cultivating high-quality P. ginseng.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12870-024-05518-2.

Keywords

Panax gingseng
Ginsenoside
Environmental factors
Secondary metabolite
http://dx.doi.org/10.13039/501100012166 National Key Research and Development Program of China 2023YFC3503802 2023YFC3503802 2023YFC3503802 2023YFC3503802 2023YFC3503802 2023YFC3503802 2023YFC3503802 2023YFC3503802 2023YFC3503802 2023YFC3503802 Scientific and technological innovation project of China Academy of Chinese Medical SciencesCI2023E002 CI2023E002 CI2023E002 CI2023E002 CI2023E002 CI2023E002 CI2023E002 CI2023E002 CI2023E002 CI2023E002 China Agricultural Research System of MOF and MARACARS-21 CARS-21 CARS-21 CARS-21 CARS-21 CARS-21 CARS-21 CARS-21 CARS-21 CARS-21 issue-copyright-statement© BioMed Central Ltd., part of Springer Nature 2024
==== Body
pmcIntroduction

Panax ginseng Meyer (P. ginseng) of the Araliaceae family is an important economic crop. The roots and rhizome of P. ginseng are known as ginseng, which is widely used in Asia and Europe as a functional food, e.g. in soups and healthy tea, because it boosts vital energy [1]. Ginseng is believed to benefit the spleen and lungs, calm the mind, and improve intelligence. Modern studies show it treats metabolic syndromes, cardiovascular disease, and has antitumor, antifatigue, antidepressant, and anti-inflammatory effects [2, 3]. Ginsenosides was one of the most important group of compounds related to pharmacological effects of ginseng, such as Rb1 [4], Rb2 [5], Rb3 [6], Rd2 [7] etc., making these compounds more attentioned when evaluating the quality of P. ginseng raw products.

Wild P. ginseng, once widespread in Russia, China, and Korea, now only survives in Jilin, Heilongjiang, and parts of Primorsky due to overexploitation and environmental changes [8]. And previous research shows that the quality of P. ginseng varies depending on the place of origin, demonstrating that environment have a significant impact on the distribution and ginsenoside of ginseng populations. The ginsenoside of P. ginseng vary significantly between different plantations, which is closely related to the physicochemical properties of the local soil [9], and the climate factors. For example, contents of Fe, available phosphorus, available potassium, total nitrogen, organic maters, and soil pH could promote the increase of ginsenoside content [10, 11], with Fe particularly proven to upregulate proteins involved in the ginsenoside synthesis pathway [12]. Climatic factors not only affect the migration of ginseng distribution areas but also influence ginsenoside accumulation. For example, low-temperature stimulation can promote the accumulation of ginsenosides [13]. Zhu et al. also found in field experiments that temperature and precipitation affect the accumulation of ginsenosides [14]. The health-promoting efficacy of P. ginseng have been shown to be influenced by environmental factors such as photosynthetically active radiation (PAR), soil water potential [15]. However, research on the relationship between ginseng habitats and its quality has been limited to its natural distribution in northeast China, parts of Russia and Korea [8], and studies on a larger scale are lacking.

In this study, we investigated the climatic and soil factors that potentially affect the ginsenosides content of P. ginseng in a larger geographical scale. The aim of this study was to evaluate: (1) the differences in ginsenosides contents of P. ginseng from different geo-origins, (2) the correlation between ginsenosides contents and climate and soil factors, and (3) the positive mechanisms how environments affect ginsenosides properties of P. ginseng in a raletive larger scale. We conducted correlation analyses between the concentration of ginsenoside compounds, climatic factors and physicochemical properties of the soil. The results will provide guidance for the identification of optimal locations for the cultivation of P. ginseng and the establishment of appropriate soil treatment systems or soil administrating regimes.

Methods

Collection and preparation of plant and soil samples

We collected five-year-old P. ginseng samples and their associated soil samples from 10 locations in three provinces in China. Samples were identified as P. ginseng by Prof. Guo Lanping (National Resource Center for Chinese Materia Medica), and voucher specimens were deposited in the National Resource Center for Chinese Materia Medica, China Academy of Chinese Medical Sciences (Beijing). These locations ranged in elevation from 182 to 2861 m, as shown in Fig. 1 and Table S1. Each P. ginseng plant was excavated together with the surrounding soil 10 cm in diameter and 20 cm deep using a hand shovel. The soil was then carefully removed from the plant. Both the plants and the soil were collected, labelled accordingly, placed on ice, and quickly transported to the laboratory. And 9 plants and soil samples was collected from each location. After the P. ginseng samples were collected, they were thoroughly washed, The parameters of the freeze dryer is acuum gauge 8.0 Pa, trap temperature − 60 ℃, and dried for 72 h, pulverised with a MM400 ball mill, and stored in a refrigerator at 4 °C for later use. And the soil samples were air-dried for subsequent analyses.

Fig. 1 Locations of sample collection

Reagents and chemicals

Methanol (chromatographic grade, UN1230) and acetonitrile (chromatographic grade, UN1648) were purchased from Merck KGaA, Germany. Ultrapure water was produced in-house using a Thermo Barnstead Gen Pure UV/UF system from Thermo Fisher Scientific, USA. The samples were filtered through 0.22 μm hydrophobic PTFE syringe filters (Millipore, USA).

Reference compounds of ginsenosides Rb1 (batch no. 19112451), Rb2 (batch no. 19092591), Rb3 (batch no. 19082903), Rc (batch no. 19111861), Rd (batch no. 19110713), Rd2 (batch no. 19102489), Re (batch no. 19092451), Rf (batch no. 19102798), Rg1 (batch no. 19122514), 20 S-Rg2 (batch no. 191102401), 20 S-Rh1 (batch no. 19101823), and Ro (batch no. 19111720) were all obtained from Shanghai YuanYe Biotechnology Co., Ltd., with a purity of ≥ 98%.

Sample and standard solutions preparation

Sample extraction was conducted following the protocols of Li et al. [16] with minor modifications: Approximately 0.15 g of P. ginseng powder was weighed (with the exact weight recorded as ‘mg’) and homogenized with 1.5 mL of 70% (V/V) methanol solution. The mixture was then subjected to ultrasonic extraction at 250 W, 40 kHz for 60 min. After cooling, add 70% (V/V) methanol to bring the volume back up to 1.5 mL. The mixture was then centrifuged at 13,000 rpm for 10 min. The supernatant was collected and filtered through a 0.22 μm microporous membrane and then diluted five fold with filtered 70% (V/V) methanol as the extracted sample solution.

For the preparation of the reference standard solutions: Accurately weigh an appropriate amount of each reference standard and dissolve it in methanol. Vortex to mix and prepare a 2 mg/mL stock solution. Dilute aliquots of this stock solution with methanol to obtain a range of concentrations for a mixed standard solution, which is then stored at 4 °C for later studies.

UPLC method

Chromatographic separation was performed using a Waters Acquity UPLC I-Class system (Waters Corp., Milford, MA, USA) equipped with a Waters ACQUITY UPLC HSS T3 column (2.1 mm × 100 mm, 1.8 μm). The column temperature was maintained at 40 °C with a flow rate of 0.4 mL/min and an injection volume of 1 µL. The mobile phase consisted of 0.1% formic acid in acetonitrile (A) and 0.1% formic acid in water (B). The elution gradient was set as follows: 0–2.5 min, 28–31% A; 2.5–3 min, 31–35% A; 3–4.5 min, 35–36% A; 4.5–6 min, 36–37% A; 6–7 min, 37–50% A; 7–9 min, 50–98% A; 9–11 min, 98% A; 11–12 min, 98–28% A; 12–15 min, 28% A.

Measurement of soil physiochemical properties

The soil pH was determined using a CO2-free water extraction-potentiometric method [17]. Soil organic matter (SOM) was quantified using the dichromate oxidation method with a heating plate [18]. Alkali-hydrolyzable nitrogen (TN) was measured using the 1 mol/L NaOH alkali dissolution gas diffusion method [19]. Available phosphorus (AP) was determined using a 0.5 mol/L sodium bicarbonate extraction-molybdenum antimony anticolorimetric method [20]. Available potassium (AK) was determined by 1 mol/L ammonium acetate extraction followed by atomic absorption spectrophotometry [21]. Available iron (Fe) was determined by extraction with DTPA solution and atomic absorption spectrophotometry. Available boron (B) was determined by extraction with hot water and azomethine-H colorimetry [22]. Three technical replicates were performed for each measurement.

Acquisition of climate data

The selection of all climatic factors and the statistical approach for the (partial) Mantel tests in this study were based on the methodology of Liu et al. [23]. Climate data for all ten sampling locations, including mean annual temperature (MAT), temperature seasonality (coefficient of variation, TS), maximum temperature of the warmest month (MT), minimum temperature of the coldest month (MIT), annual temperature range (TR), mean annual precipitation (MAP), and precipitation seasonality (PS), were obtained from the WorldClim database (http://www.worldclim.org). Detailed climatic data can be found in Table S1 [24].

Statistical analyses

In this study, chemometric techniques were applied to evaluate the effects of environmental factors on the levels of secondary metabolites of P. ginseng. Our approach included the Kruskal-Wallis H test, principal component analysis (PCA) and hierarchical cluster analysis (HCA), complemented by multiple regression and path analyses, and Mantel tests to understand variable correlations. Due to the significant differences in size between the measurements of the variables, a log transformation was performed when calculating the Pearson correlation coefficients in order to ensure the normality of the data. The Kruskal-Wallis H-test and regression analyses were performed with SPSS. The principal component analyses (PCA) were performed with SIMCA 14.1. Hierarchical cluster analyses (HCA) were performed in R with the package “stats”. Conduct path analysis using AMOS software. For more details, see the Supplementary Material. Data processing was performed with the library “Pandas” [25]. The graphical visualisation of the Spearman correlation coefficients was created with the libraries “Seaborn” and “Matplotlib“ [26, 27]. The Pearson correlation coefficients were calculated and the linear relationships between the variables were evaluated using the “Pandas” library [25]. The significance of the correlation was displayed in a co-occurrence network created with “NetworkX” and “Matplotlib” [28]. For the analyses of correlations between different distance matrices derived from our dataset, Mantel tests used “scipy.spatial.distance” to calculate Euclidean and Bray-Curtis distances with respect to secondary metabolites, soil properties [29] and climatic factors. The Spearman rank correlation coefficients for these matrices were determined using “scipy.stats.spearmanr”, with the partial Mantel test taking into account the influence of a third matrix using “sklearn.metrics.pairwise_distances” [30], with permutation tests used in both cases to check significance. R (v4.3.1) facilitated statistical modelling using random forest models to determine the predictive value of the environmental variables on the response metrics. Here, “randomForest” and “caret” were used for model development and validation [31, 32] and “rfUtilities” and “rfPermute” for significance testing of the predictors [33, 34]. The graphical visualisation was performed in R with “ggplot2” [35].

Results

Ginsenoside differentiation of P. ginseng and associated driving factors

A total of 12 metabolites were identified in P. ginseng from various origins. The current study began with a rigorous quantitative analysis of twelve ginsenoside compounds in these P. ginseng samples, combining these data with soil physicochemical properties. Using Shapiro-Wilk and Levene tests, we assessed the normality and homogeneity of variance in our data sets. As shown in Table S2, deviations from the normal distribution and homogeneity of variance required the use of non-parametric statistical approaches. Consequently, the Kruskal-Wallis H test, supplemented by Dunn’s post-hoc test, was utilized to determine the statistical significance of the differences between the different groups, as shown in Fig. 2.

Fig. 2 Summary of the edaphic properties and Ginsenoside across the sampling groups from different sites in China. JA, Jian (Jilin); CB, Changbai (Jilin); FS, Fusong (Jilin); SY, Shiyan (Hubei); YL, Yulong (Yunnan); pH, soil pH; SOM, soil organic matter; TN, total nitrogen; AP, available phosphorus; AK, available potassium; Fe, iron; and B, boron. Different letters indicate significant differences within the same boxplot (P < 0.05; Kruskal-Wallis H Test, followed by Dunn’s Test for post-hoc comparisons). Bold P-values indicate significant differences (P < 0.05) (n = 9)

Considerable discrepancies were found in the physico-chemical properties of the soil, with the acidic nature being a common feature of all samples. Notably, the pH of the soil in the JA group was identified as the highest. The concentration ranges for total nitrogen (TN), available phosphorus (AP) and available potassium (AK) were determined as 26.47-664.83 mg/kg, 14.30-81.98 mg/kg and 72.10–420.00 mg/kg, respectively, with soil organic matter (SOM) content ranging from 9.96 to 250.41 g/kg. A distinctive elevation in TN, AP and AK concentrations was observed in the CB group, whereas the YL group showed a significant increase in SOM content.

Within the analyzed ginsenoside spectrum, ginsenoside Re proved to be the most prevalent ginsenoside, whereas ginsenoside Rd2 was detected in the lowest concentration, with the CB group having the highest levels of ginsenoside Rb1. The analysis revealed no significant differences between the groups for ginsenosides Rg1, Rf, Rg2 and Rh1. To mitigate possible interference and avoid false-negative implications in the cluster analysis, these four ginsenoside components were excluded from further PCA. This analytical endeavor highlights the extensive environmental variations between sampling sites and facilitates the refinement of variables for subsequent analysis, ensuring a targeted and effective investigation.

Variation of ginsenoside content in P. ginseng cultivated under diverse environmental conditions

PCA is an effective tool for dimensionality reduction in data analysis, facilitating the preservation of original data variability through a transformed coordinate system. In our study of P. ginseng, PCA was strategically used to detect patterns within 8 significantly varying ginsenosides at different sampling sites, resulting in the derivation of 8 principal components that together comprise 100% of the variance, as shown in Table S5. The first 3 components, accounting for 91.400% of the total variance, enabled the construction of a three-dimensional score plot (Fig. 3a), visually dividing the samples into distinct clusters along the first 3 principal components and highlighting the effectiveness of the method in data clustering.

Fig. 3 (a) Unsupervised analysis of ginsenosides by PCA. (b) Dendrograms resulting from HCA of P. ginseng from different locations (k = 3). (c) Co-occurrence map of Pearson correlation coefficients between environmental factors and ginsenosides. The color stands for the strength of the correlation, solid lines indicate a P ≤ 0.05 and dashed lines a P >0.05. (d) Correlation analysis between soil properties, climatic characteristics and secondary metabolites of roots, determined using the Mantel test. The color of the line represents the significance of the differences (Mantel’s P). The pattern of the dot represents the significance of the differences (Spearman’s P). The size of the line and the dot represents the size of the correlation coefficients (Mantel’s r and Spearman’s r)

Despite the usefulness of PCA in revealing data structures, it lacks the functionality of a clustering algorithm. This limitation necessitated the integration of HCA, a method that groups entities based on their similarity or distance, to further stratify the data. Utilizing Ward’s method with Euclidean distance as the similarity metric [31], we aimed to minimize the variance within clusters and thereby improve the homogeneity of samples within clusters. The preliminary PCA-guided clustering suggested a five-cluster model, which, after further examination using the Gap statistic [32], revealed a more refined scheme with three clusters to be optimal, combining superior clustering efficiency with minimized error (Table S6, Fig. 3b). In this model, the YL and JA samples in particular were divided into separate clusters, while the CB and FS samples were combined into a third cluster. In line with the findings of the HCA, the K-Means clustering analysis using the elbow method identified 3 as the optimal number of clusters, which corresponds to the result of the HCA (Fig. S2).

The results of PCA and HCA show that P. ginseng from different countries of origin forms different clusters due to the different content of twelve metabolites, which illustrates the chemical differentiation of P. ginseng from different areas. In addition, the loadings of principal components are crucial for evaluating the contribution of each metabolite to the separation of the clusters. Notably, ginsenoside Rb2 with a loading of -0.541 on the first principal component (PC1) and ginsenoside Rd2 with a loading of 0.902 on the second principal component (PC2) were significant contributors to their respective axes. In addition, ginsenoside Rd with a high loading of -0.776 on the third principal component (PC3) emphasises the differential contribution of these saponins across different geographical origins, suggesting their potential as markers for origin discrimination.

Environmental determinants influencing the total content of twelve ginsenosides

Once the differences in ginsenoside content between the different groups of P. ginseng had been confirmed, further investigations were carried out to explore the possible causes of these differences. Using (partial) Mantel tests, we analyzed the influence of soil physicochemical properties and climatic factors on ginsenoside content (Table S7). A slight positive correlation was found between the ginsenoside content and both soil physicochemical properties of the soil and the climatic factors. However, when controlling for climatic factors, the relationship between ginsenoside content and soil physicochemical properties may no longer be significant. The co-occurrence matrix of Pearson correlation coefficients between environmental factors and ginsenoside content showed significant positive correlations with variables such as Fe, TS, longitude, MT, latitude, B, AP, PS and TR; conversely, MAP, pH, altitude and MIT showed significant negative correlations with various ginsenoside components (Fig. 3c). The analysis showed 66 coefficients with p-values less than 0.05, indicating significant correlations with ginsenoside content in P. ginseng. We focused on the top 50% of these coefficients. Specifically, in Pearson correlations, 15 were climatic, 9 geographic, and 9 edaphic; in Spearman correlations, 15 were climatic, 10 geographic, and 8 edaphic, all significantly correlated with ginsenoside content. Which is consistent with the (partial) Mantel test results. Similar trends were observed in the Spearman’s rank correlation results (Fig. 3d).

The Mantel test results (Table S8) revealed significant positive correlations between AP, Fe and pH in soil physicochemical properties and ginsenoside content, while all other climatic factors except MAT and MIT showed significant positive correlations with ginsenoside content. PCA was employed to reduce the dimensionality of the matrix of ginsenoside content. The extracted PCA1 was subjected to multiple linear regression analysis with environmental factors identified as significantly different by Mantel tests, revealing significant correlations between AP and TS with ginsenoside content PCA1 (Table S9). To further clarify the effect of these two environmental factors on ginsenoside content, scatter plots were generated with ΔAP and ΔTS as x-axes and the Euclidean distance of ginsenoside content as y-axis (Fig. S3), which showed a positive correlation between AP and TS with ginsenoside content.

This analytical section, using (partial) Mantel tests, Mantel tests, and multiple linear regression analysis, confirmed the two most important environmental factors affecting the 12 measured ginsenosides: AP and TS.

The most important environmental factors influencing the individual categories of ginsenosides

To clarify the relationships between various environmental factors and the content of ginsenosides in P. ginseng, Random Forest was used to determine the environmental factors that significantly affect the content of each ginsenoside, where ginsenosides are the dependent variable and environmental factors are the independent variables. The results, as illustrated in Fig. 4, indicate that Rg2 was significantly affected by factors such as pH, AP, PS, and altitude; Rb2 by TS, MT, longitude, latitude, altitude, TR, and PS; Rd by pH and longitude; Rb3 by Fe, AP, longitude, and latitude; Ro by TR, TS, longitude, and latitude; Rb1 by latitude, TS, and TR; and Rd2 by AP, latitude, TR, B, PS, altitude, and TS.

Fig. 4 Random forest mean predictor importance of environmental characteristics (climatic and soil) and geographical factors for ginsenoside content. “%IncMSE” stands for “Percentage Increase in Mean Squared Error” ***P < 0.001; **P < 0.01; *P < 0.05; ns, not significant. The fitting degree (R2) of the random forest model is indicated

To validate the random forest results and further investigate the influence of environmental factors on ginsenoside content, we performed path analysis to calculate both direct and indirect effects of these factors on secondary metabolites. This process involved two main steps: First, we conducted stepwise regression with ginsenosides as the dependent variable and environmental factors as independent variables (Table S10). This helped identify the dominant ecological factors for each ginsenoside component. Next, we calculated the direct and indirect effects of these identified ecological factors on ginsenosides, as shown in Table 1. And the results indicated that AP had a significant direct effect on Rg2, while TN had a significant indirect effect. PS and TS had significant direct and indirect effects on Rb2. TS had a significant direct effect on Rd, whereas pH had a significant indirect effect. AP and MT had significant direct and indirect effects on Rb3. TR and latitude significantly affected Ro both directly and indirectly. AP and MT also had significant direct and indirect effects on Rc. TS and altitude significantly influenced Rb1 both directly and indirectly. B had a significant direct effect on Rd2, while AK and AP had significant indirect effects.

Table 1 Path analysis between environmental factors and secondary metabolites of P. ginseng

Items	Fcators	Direct path coefficients	Indirect path coefficient	
Rg2			Total	→AP	→TN		
	AP	0.288***	0.097		0.097***		
	TN	-0.008	0.338	0.338			
Rb2			Total	→PS	→TS		
	PS	0.148***	0.001		0.123***		
	TS	0.337***	0.001	0.279***			
Rd			Total	→TS	→pH		
	TS	0.301***	-0.091		-0.091***		
	pH	-0.095	0.000	0.000			
Rb3			Total	→AP	→MT		
	AP	0.281***	0.000		0.000***		
	MT	0.259***	0.037	0.037***			
Ro			Total	→TR	→latitude		
	TR	0.135***	0.134		0.134***		
	latitude	0.123**	0.123	0.123**			
Rc			Total	→AP	→MT		
	AP	0.125**	0.018		0.018**		
	MT	0.211***	0.030	0.030***			
Rb1			Total	→Fe	→TS	→altitude	
	Fe	0.009	-0.003		-0.003	0.000	
	TS	0.095**	0.000	0.000**		0.000**	
	altitude	0.130***		0.35***	0.054***		
Rd2			Total	→B	→AK	→AP	
	B	0.316***	0.166		0.086***	0.080***	
	AK	0.015	0.011	0.005		0.006	
	AP	0.017		0.006	0.006		
Note: pH: soil pH; SOM: soil organic matter; TN: total nitrogen; AP: available phosphorus; AK: available potassium; Fe: iron; and B: boron; MAT: Mean Annual Temperature; TS: Temperature Seasonality; MT: Maximum Temperature of Warmest Month; MIT: Minimum Temperature of Coldest Month; TR: Temperature Annual Range; MAP: Mean Annual Precipitation; PS: Precipitation Seasonality. The symbol “→” represents the causal path in path analysis: indicating the indirect effect of the variable following the “→” on the Items variable via the Factors variable

Discussion

Chemotypes of P. ginseng from different habitats

In the current study, soil properties and climatic variables had a significant effect on the concentration of secondary metabolites in P. ginseng, which was similar with a previous study that statses the importance of microenvironment in shaping active compounents of P. ginseng [36]. P. ginseng from different origins were different from each other in secondary metabolites according to PCA and a further HCA and K-mean indicates that they could be clustered as 3 chemotypes. The differences in chemical characteristics could be affected by the environmental differences. Such environmentally induced chemotypic differentiation is also uncovered in Rh. tanguticum, which exhibits significant variations in secondary metabolitesbetween canopy-covered and open habitats [37]. Similar reports exist in the current research on ginseng [38], our study innovatively expands the distribution range of ginseng in its natural habitat, collecting samples from a wider area and covering different climates and altitudes to the best of our ability. This provides scientific reference for predicting high-quality ginseng cultivation sites in the context of global climate change.

The complex interplay of environmental factors correlate with secondary metabolites in P. ginseng

Partial Mantel’s tests was employed to decipher the independent effect of soil properties or climates on secondary metabolite concentrations and their interaction with climatic factors [23]. In the current study climate and soil were found to have joint influences on secondary metabolite content in P. ginseng, with climate likely being the stronger driving factor. This was also supported by the following Pearson’s correlation and Spearman’s rank correlation. Soil properties was also thought to have a media role, in which climate affect secondary metabolites of P. ginseng through soil microbial community and soil properties [39]. Similar results have been reported for different species within the genus Panax, such as P. notoginseng [40] and P. quinquefolius [41], indicating that climatic factors significantly influence ginsenoside content. This emphasises the broad applicability of environmental influences on ginsenoside profiles throughout the genus. The Mantel test, Pearson’s and Spearman’s correlations indicate that both low and high temperatures, such as MT, MIT, and TR, promote ginsenoside accumulation in ginseng. This may be due to ginsenosides’ antioxidant activity protecting ginseng cells from oxidative stress and free radical damage under extreme temperatures [42]. Thus, ginsenosides likely act as a defense mechanism, helping ginseng withstand various stresses, which fit the theory that stressed environments was suitable for producing plant medicines with high active compounds [43]. The negative correlation between pH and ginsenoside content is also an interesting result, indicating that more acidic environments lead to greater ginsenoside accumulation. This could be caused by the upregulated ginsenosides’ synthesis to resist acid stress [44]. However, prvious studies has also revealed continuous ginseng cropping acidifies soil [45], and the allelopathic effect of ginsenosides may also causing regional death of cultivated ginseng and continuous cropping obstacles [46].

Key factors influencing secondary metabolites in P. ginseng

The Mantel test identifies AP and TS as key environmental factors significantly influencing ginsenoside levels. Using PCA for dimensionality reduction, a multiple linear regression analysis reveals their positive impact on PCA1 of the ginsenoside data. Further analysis with path and random forest models shows AP’s direct effect on Rg2 and Rb3 and its indirect effect on Rd2.

Previous studies have shown that phosphorus can affect ginsenoside content in Panax quinquefolius and P. ginseng, which is consistent with the results of the current study [47, 48]. Phosphorus is also a crucial element for plant growth and has been shown to influence the synthesis of various secondary metabolites in different plant genera [49], for example, glycoside content in Stevia rebaudiana Bertoni and flavonoids in maize [50, 51]. This could be due to the fact that phosphorus conditions could influence different gene expressions in the synthesis of regulatory secondary metabolites [52]. In contrast, TS showed significant direct and indirect effects on Rb2, but only direct effects on Rb1. This suggests that TS may be a key factor affecting the quality of P. ginseng. The ginsenoside content of Panax quinquefolius, a close relative of ginseng [53], is also significantly affected by TS [54]. An influence of TS on the secondary metabolite content of Salvia miltiorrhiza and Pilocarpus pennatifolius (Rutaceae) was also found [55, 56]. Previous studies have shown that TS is crucial for the geographical distribution of plants and is essential for plant growth [57]. This indicator shows a predominantly positive correlation with ginsenoside content, suggesting that more significant seasonal temperature variations may be beneficial for the quality of ginseng.

Conclusions

This study employs a combination of methods including PCA, HCA, Mantel tests, path analysis and Random Forest to investigate the correlations between environmental factors and ginsenoside content in cultivated P. ginseng. The results indicate that climatic characteristics have a more significant influence on ginsenoside content than edaphic properties. Of all the factors analyzed, AP and TS exert the greatest influence on ginsenoside content. In particular, TS shows a significant correlation with the primary active compounds of P. ginseng, suggesting that it may be a key determinant in the formation of P. ginseng quality. Our research provides scientific reference for predicting high-quality ginseng cultivation sites.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Acknowledgements

Reviewers are acknowledged for their contribution to the improvement of the manuscript in the revision process.

Author contributions

Dehua Wu: Conceptualization, Methodology, Writing original draft preparation. Feng Xiong: Methodology, original draft preparation. Siqi Liu: Writing review. Jitong Zhu: Writing review. Dan Zhao: Writing review. Jian Yang: Writing review. Chuanzhi Kang: Writing review & editing. Hongyang Wang: Formal analysis. Wenqi Ma: Resources. Lanping Guo: Supervision, Project administration, Writing review & editing. All authors read and approved the final manuscript.

Funding

This work was financially supported by the National Key Research and Development Program of China(2023YFC3503802), Scientific and technological innovation project of China Academy of Chinese Medical Sciences (Grant No. CI2023E002, CI2021A03905), China Agricultural Research System of MOF and MARA (CARS-21).

Data availability

The datasets generated for this study are available on request to the corresponding authors.

Declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Feng Xiong as the co-first-author of this article.
==== Refs
References

1. U.S. Food and Drug Administration. Chromium picolinate; Ginkgo biloba leaf extract; and Ginseng extract. 2000. https://www.cfsanappsexternal.fda.gov/scripts/fdcc/index.cfm?set=GRASNotices&id=36. Accessed 18 Dec 2023.
2. Zhao B Lv C Lu J Natural occurring polysaccharides from Panax ginseng C. A. Meyer: a review of isolation, structures, and bioactivities Int J Biol Macromol 2019 133 324 36 10.1016/j.ijbiomac.2019.03.229 30943421
Zhao B, Lv C, Lu J. Natural occurring polysaccharides from Panax ginseng C. A. Meyer: a review of isolation, structures, and bioactivities. Int J Biol Macromol. 2019;133:324–36. 10.1016/j.ijbiomac.2019.03.22930943421 10.1016/j.ijbiomac.2019.03.229
3. Fan W Huang Y Zheng H Li S Li Z Yuan L Ginsenosides for the treatment of metabolic syndrome and cardiovascular diseases: pharmacology and mechanisms Biomed Pharmacother 2020 132 110915 10.1016/j.biopha.2020.110915 33254433
Fan W, Huang Y, Zheng H, Li S, Li Z, Yuan L, et al. Ginsenosides for the treatment of metabolic syndrome and cardiovascular diseases: pharmacology and mechanisms. Biomed Pharmacother. 2020;132:110915. 10.1016/j.biopha.2020.11091533254433 10.1016/j.biopha.2020.110915
4. Xie W Wang X Xiao T Cao Y Wu Y Yang D Protective effects and network analysis of ginsenoside Rb1 against cerebral ischemia injury: a pharmacological review Front Pharmacol 2021 12 604811 10.3389/fphar.2021.604811 34276353
Xie W, Wang X, Xiao T, Cao Y, Wu Y, Yang D, et al. Protective effects and network analysis of ginsenoside Rb1 against cerebral ischemia injury: a pharmacological review. Front Pharmacol. 2021;12:604811. 10.3389/fphar.2021.60481134276353 10.3389/fphar.2021.604811
5. Miao L Yang Y Li Z Fang Z Zhang Y Han C Ginsenoside Rb2: a review of pharmacokinetics and pharmacological effects J Ginseng Res 2022 46 206 13 10.1016/j.jgr.2021.11.007 35509822
Miao L, Yang Y, Li Z, Fang Z, Zhang Y, Han C. Ginsenoside Rb2: a review of pharmacokinetics and pharmacological effects. J Ginseng Res. 2022;46:206–13. 10.1016/j.jgr.2021.11.00735509822 10.1016/j.jgr.2021.11.007
6. Li W Duan Y Yan X Liu X Fan M Wang Z A mini-review on pharmacological effects of ginsenoside Rb3, a marked saponin from Panax Genus BIOCELL 2022 46 1417 23 10.32604/biocell.2022.017609
Li W, Duan Y, Yan X, Liu X, Fan M, Wang Z. A mini-review on pharmacological effects of ginsenoside Rb3, a marked saponin from Panax Genus. BIOCELL. 2022;46:1417–23. 10.32604/biocell.2022.01760910.32604/biocell.2022.017609
7. Zuo X Li Q Ya F Ma L-J Tian Z Zhao M Ginsenosides Rb2 and Rd2 isolated from Panax notoginseng flowers attenuate platelet function through P2Y 12-mediated cAMP/PKA and PI3K/Akt/Erk1/2 signaling Food Funct 2021 12 5793 805 10.1039/D1FO00531F 34041517
Zuo X, Li Q, Ya F, Ma L-J, Tian Z, Zhao M, et al. Ginsenosides Rb2 and Rd2 isolated from Panax notoginseng flowers attenuate platelet function through P2Y 12-mediated cAMP/PKA and PI3K/Akt/Erk1/2 signaling. Food Funct. 2021;12:5793–805. 10.1039/D1FO00531F34041517 10.1039/D1FO00531F
8. Zhuravlev YN Koren OG Reunova GD Muzarok TI Gorpenchenko TY Kats IL Panax ginseng natural populations: their past, current state and perspectives Acta Pharmacol Sin 2008 29 1127 36 10.1111/j.1745-7254.2008.00866.x 18718182
Zhuravlev YN, Koren OG, Reunova GD, Muzarok TI, Gorpenchenko TY, Kats IL, et al. Panax ginseng natural populations: their past, current state and perspectives. Acta Pharmacol Sin. 2008;29:1127–36. 10.1111/j.1745-7254.2008.00866.x18718182 10.1111/j.1745-7254.2008.00866.x
9. Lan Y Zhang M Han M Yang L Differences in the quality, yield, and soil microecology of ginseng in different planting environments HORTICULTURAE 2023 9 520 10.3390/horticulturae9040520
Lan Y, Zhang M, Han M, Yang L. Differences in the quality, yield, and soil microecology of ginseng in different planting environments. HORTICULTURAE. 2023;9:520. 10.3390/horticulturae904052010.3390/horticulturae9040520
10. Liu C Xia R Tang M Chen X Zhong B Liu X Improved ginseng production under continuous cropping through soil health reinforcement and rhizosphere microbial manipulation with biochar: a field study of Panax ginseng from Northeast China Hortic Res 2022 9 uhac108 10.1093/hr/uhac108 35836471
Liu C, Xia R, Tang M, Chen X, Zhong B, Liu X, et al. Improved ginseng production under continuous cropping through soil health reinforcement and rhizosphere microbial manipulation with biochar: a field study of Panax ginseng from Northeast China. Hortic Res. 2022;9:uhac108. 10.1093/hr/uhac10835836471 10.1093/hr/uhac108
11. Eo J Park K-C Effect of vermicompost application on root growth and ginsenoside content of Panax ginseng J Environ Manage 2019 234 458 63 10.1016/j.jenvman.2018.12.101 30641356
Eo J, Park K-C. Effect of vermicompost application on root growth and ginsenoside content of Panax ginseng. J Environ Manage. 2019;234:458–63. 10.1016/j.jenvman.2018.12.10130641356 10.1016/j.jenvman.2018.12.101
12. Wang Q Jin Q Ma Y Zhang S Zhang L Liu Z Iron toxicity-induced regulation of key secondary metabolic processes associated with the quality and resistance of Panax ginseng and Panax quinquefolius Ecotoxicol Environ Saf 2021 224 112648 10.1016/j.ecoenv.2021.112648 34450425
Wang Q, Jin Q, Ma Y, Zhang S, Zhang L, Liu Z, et al. Iron toxicity-induced regulation of key secondary metabolic processes associated with the quality and resistance of Panax ginseng and Panax quinquefolius. Ecotoxicol Environ Saf. 2021;224:112648. 10.1016/j.ecoenv.2021.11264834450425 10.1016/j.ecoenv.2021.112648
13. Wang S Liang W Yao L Wang J Gao W Effect of temperature on morphology, ginsenosides biosynthesis, functional genes, and transcriptional factors expression in Panax ginseng adventitious roots J Food Biochem 2019 43 e12794 10.1111/jfbc.12794 31353579
Wang S, Liang W, Yao L, Wang J, Gao W. Effect of temperature on morphology, ginsenosides biosynthesis, functional genes, and transcriptional factors expression in Panax ginseng adventitious roots. J Food Biochem. 2019;43:e12794. 10.1111/jfbc.1279431353579 10.1111/jfbc.12794
14. Zhu L Xu L Huang Y Xie C Dou D Xu J Correlations between ecological factors and the chemical compositions of mountainous forest cultivated ginseng J Food Compos Anal 2022 114 104867 10.1016/j.jfca.2022.104867
Zhu L, Xu L, Huang Y, Xie C, Dou D, Xu J. Correlations between ecological factors and the chemical compositions of mountainous forest cultivated ginseng. J Food Compos Anal. 2022;114:104867. 10.1016/j.jfca.2022.10486710.1016/j.jfca.2022.104867
15. Zhang T Han M Yang L Han Z Cheng L Sun Z The effects of environmental factors on ginsenoside biosynthetic enzyme gene expression and saponin abundance Molecules 2019 24 14 10.3390/molecules24010014
Zhang T, Han M, Yang L, Han Z, Cheng L, Sun Z, et al. The effects of environmental factors on ginsenoside biosynthetic enzyme gene expression and saponin abundance. Molecules. 2019;24:14. 10.3390/molecules2401001410.3390/molecules24010014
16. Li S-L Lai S-F Song J-Z Qiao C-F Liu X Zhou Y Decocting-induced chemical transformations and global quality of Du-Shen-Tang, the decoction of ginseng evaluated by UPLC-Q-TOF-MS/MS based chemical profiling approach J Pharm Biomed Anal 2010 53 946 57 10.1016/j.jpba.2010.07.001 20667431
Li S-L, Lai S-F, Song J-Z, Qiao C-F, Liu X, Zhou Y, et al. Decocting-induced chemical transformations and global quality of Du-Shen-Tang, the decoction of ginseng evaluated by UPLC-Q-TOF-MS/MS based chemical profiling approach. J Pharm Biomed Anal. 2010;53:946–57. 10.1016/j.jpba.2010.07.00120667431 10.1016/j.jpba.2010.07.001
17. Bao SD Agricultural and chemistry analysis of soil 2005 Beijing Agric Press Beijing China
Bao SD. Agricultural and chemistry analysis of soil. Beijing: Agric Press Beijing China; 2005.
18. Egli M Sartori G Mirabella A Favilli F Giaccai D Delbos E Effect of north and south exposure on organic matter in high Alpine soils Geoderma 2009 149 124 36 10.1016/j.geoderma.2008.11.027
Egli M, Sartori G, Mirabella A, Favilli F, Giaccai D, Delbos E. Effect of north and south exposure on organic matter in high Alpine soils. Geoderma. 2009;149:124–36. 10.1016/j.geoderma.2008.11.02710.1016/j.geoderma.2008.11.027
19. Bardelli T Gómez-Brandón M Ascher-Jenull J Fornasier F Arfaioli P Francioli D Effects of slope exposure on soil physico-chemical and microbiological properties along an altitudinal climosequence in the Italian Alps Sci Total Environ 2017 575 1041 55 10.1016/j.scitotenv.2016.09.176 27692937
Bardelli T, Gómez-Brandón M, Ascher-Jenull J, Fornasier F, Arfaioli P, Francioli D, et al. Effects of slope exposure on soil physico-chemical and microbiological properties along an altitudinal climosequence in the Italian Alps. Sci Total Environ. 2017;575:1041–55. 10.1016/j.scitotenv.2016.09.17627692937 10.1016/j.scitotenv.2016.09.176
20. Murphy PNC Stevens RJ Christie P Long-term application of animal slurries to grassland alters soil cation balance Soil Use Manag 2005 21 240 4 10.1111/j.1475-2743.2005.tb00130.x
Murphy PNC, Stevens RJ, Christie P. Long-term application of animal slurries to grassland alters soil cation balance. Soil Use Manag. 2005;21:240–4. 10.1111/j.1475-2743.2005.tb00130.x10.1111/j.1475-2743.2005.tb00130.x
21. Liu Y Zhang L Lu J Chen W Wei G Lin Y Topography affects the soil conditions and bacterial communities along a restoration gradient on Loess-Plateau Appl Soil Ecol 2020 150 103471 10.1016/j.apsoil.2019.103471
Liu Y, Zhang L, Lu J, Chen W, Wei G, Lin Y. Topography affects the soil conditions and bacterial communities along a restoration gradient on Loess-Plateau. Appl Soil Ecol. 2020;150:103471. 10.1016/j.apsoil.2019.10347110.1016/j.apsoil.2019.103471
22. Fageria NK Baligar VC Clark RB Micronutrients in crop production Adv Agron 2002 77 185 268 10.1016/S0065-2113(02)77015-6
Fageria NK, Baligar VC, Clark RB. Micronutrients in crop production. Adv Agron. 2002;77:185–268. 10.1016/S0065-2113(02)77015-610.1016/S0065-2113(02)77015-6
23. Liu Y Wang H Peng Z Li D Chen W Jiao S Regulation of root secondary metabolites by partial root-associated microbiotas under the shaping of licorice ecotypic differentiation in northwest China J Integr Plant Biol 2021 63 2093 109 10.1111/jipb.13179 34655272
Liu Y, Wang H, Peng Z, Li D, Chen W, Jiao S, et al. Regulation of root secondary metabolites by partial root-associated microbiotas under the shaping of licorice ecotypic differentiation in northwest China. J Integr Plant Biol. 2021;63:2093–109. 10.1111/jipb.1317934655272 10.1111/jipb.13179
24. Hijmans RJ Cameron SE Parra JL Jones PG Jarvis A Very high resolution interpolated climate surfaces for global land areas Int J Climatol 2005 25 1965 78 10.1002/joc.1276
Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A. Very high resolution interpolated climate surfaces for global land areas. Int J Climatol. 2005;25:1965–78. 10.1002/joc.127610.1002/joc.1276
25. Reback J McKinney W Van Den Bossche J Augspurger T Cloud P Klein A Pandas-dev/pandas: pandas 1.0.5 Zenodo 2020 10.5281/zenodo.3898987
Reback J, McKinney W, Van Den Bossche J, Augspurger T, Cloud P, Klein A, et al. Pandas-dev/pandas: pandas 1.0.5. Zenodo. 2020. 10.5281/zenodo.389898710.5281/zenodo.3898987
26. Waskom ML Seaborn: statistical data visualization J Open Source Softw 2021 6 3021 10.21105/joss.03021
Waskom ML. Seaborn: statistical data visualization. J Open Source Softw. 2021;6:3021. 10.21105/joss.0302110.21105/joss.03021
27. Hunter JD Matplotlib: a 2D graphics environment Comput Sci Eng 2007 9 90 5 10.1109/MCSE.2007.55
Hunter JD. Matplotlib: a 2D graphics environment. Comput Sci Eng. 2007;9:90–5. 10.1109/MCSE.2007.5510.1109/MCSE.2007.55
28. Hagberg A, Swart PJ, Schult DA. Exploring network structure, dynamics, and function using NetworkX. Proceedings of the 7th Python in Science Conference (SciPy 2008). 2008;11–16. https://www.osti.gov/biblio/960616
29. Virtanen P Gommers R Oliphant TE Haberland M Reddy T Cournapeau D SciPy 1.0: fundamental algorithms for scientific computing in Python Nat Methods 2020 17 261 72 10.1038/s41592-019-0686-2 32015543
Virtanen P, Gommers R, Oliphant TE, Haberland M, Reddy T, Cournapeau D, et al. SciPy 1.0: fundamental algorithms for scientific computing in Python. Nat Methods. 2020;17:261–72. 10.1038/s41592-019-0686-232015543 10.1038/s41592-019-0686-2
30. Pedregosa F Varoquaux G Gramfort A Michel V Thirion B Grisel O Scikit-learn: machine learning in Python J Mach Learn Res 2011 12 2825 30 10.5555/1953048.2078195
Pedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, Grisel O, et al. Scikit-learn: machine learning in Python. J Mach Learn Res. 2011;12:2825–30. 10.5555/1953048.207819510.5555/1953048.2078195
31. Svetnik V Liaw A Tong C Culberson JC Sheridan RP Feuston BP Random forest: a classification and regression tool for compound classification and QSAR modeling J Chem Inf Comput Sci 2003 43 1947 58 10.1021/ci034160g 14632445
Svetnik V, Liaw A, Tong C, Culberson JC, Sheridan RP, Feuston BP. Random forest: a classification and regression tool for compound classification and QSAR modeling. J Chem Inf Comput Sci. 2003;43:1947–58. 10.1021/ci034160g14632445 10.1021/ci034160g
32. Kuhn M Building predictive models in R using the caret package J Stat Softw 2008 28 1 26 10.18637/jss.v028.i05 27774042
Kuhn M. Building predictive models in R using the caret package. J Stat Softw. 2008;28:1–26. 10.18637/jss.v028.i0527774042 10.18637/jss.v028.i05
33. Evans JS, Murphy MA. Package ‘rfUtilities.’ R Core Team Vienna Austria. 2019. https://github.com/jeffreyevans/rfUtilities
34. Archer E. rfPermute: estimate permutation p-values for random forest importance metrics. 2016. https://search.r-project.org/CRAN/refmans/rfPermute/html/rfPermute.html
35. Wickham H. ggplot2. WIREs Comput. Stat. 2011;3:180–5. 10.1002/wics.147
36. Yang L Wen K-S Ruan X Zhao Y-X Wei F Wang Q Response of plant secondary metabolites to environmental factors Molecules 2018 23 762 10.3390/molecules23040762 29584636
Yang L, Wen K-S, Ruan X, Zhao Y-X, Wei F, Wang Q. Response of plant secondary metabolites to environmental factors. Molecules. 2018;23:762. 10.3390/molecules2304076229584636 10.3390/molecules23040762
37. Xiong F Nie X Yang L Wang L Li J Zhou G Non-target metabolomics revealed the differences between Rh. Tanguticum plants growing under canopy and open habitats BMC Plant Biol 2021 21 119 10.1186/s12870-021-02897-8 33639841
Xiong F, Nie X, Yang L, Wang L, Li J, Zhou G. Non-target metabolomics revealed the differences between Rh. Tanguticum plants growing under canopy and open habitats. BMC Plant Biol. 2021;21:119. 10.1186/s12870-021-02897-833639841 10.1186/s12870-021-02897-8
38. Kang J Lee S Kang S Kwon HN Park JH Kwon SW NMR-Based metabolomics approach for the differentiation of ginseng (Panax ginseng) roots from different origins Arch Pharm Res 2008 31 330 6 10.1007/s12272-001-1160-2 18409046
Kang J, Lee S, Kang S, Kwon HN, Park JH, Kwon SW, et al. NMR-Based metabolomics approach for the differentiation of ginseng (Panax ginseng) roots from different origins. Arch Pharm Res. 2008;31:330–6. 10.1007/s12272-001-1160-218409046 10.1007/s12272-001-1160-2
39. Puissant J Cecillon L Mills RTE Robroek BJM Gavazov K De Danieli S Seasonal influence of climate manipulation on microbial community structure and function in mountain soils Soil Biol Biochem 2015 80 296 305 10.1016/j.soilbio.2014.10.013
Puissant J, Cecillon L, Mills RTE, Robroek BJM, Gavazov K, De Danieli S, et al. Seasonal influence of climate manipulation on microbial community structure and function in mountain soils. Soil Biol Biochem. 2015;80:296–305. 10.1016/j.soilbio.2014.10.01310.1016/j.soilbio.2014.10.013
40. Liu C Zuo Z Xu F Wang Y Study of the suitable climate factors and geographical origins traceability of Panax notoginseng based on correlation analysis and spectral images combined with machine learning Front Plant Sci 2023 13 1009727 10.3389/fpls.2022.1009727 36825249
Liu C, Zuo Z, Xu F, Wang Y. Study of the suitable climate factors and geographical origins traceability of Panax notoginseng based on correlation analysis and spectral images combined with machine learning. Front Plant Sci. 2023;13:1009727. 10.3389/fpls.2022.100972736825249 10.3389/fpls.2022.1009727
41. Zhang X, Tian L, Zhang G, Huang L. Revealing the impact of the environment on the ginsenosides in American ginseng from three scales: bioclimatic factors, soil physicochemical properties and microbial communities. 2022. 10.21203/rs.3.rs-1860897/v1
42. Yamabe N Kim Y-J Lee S Cho E-J Park S-H Ham J Increase in antioxidant and anticancer effects of ginsenoside re-lysine mixture by Mail lard reaction Food Chem 2013 138 876 83 10.1016/j.foodchem.2012.12.004 23411191
Yamabe N, Kim Y-J, Lee S, Cho E-J, Park S-H, Ham J, et al. Increase in antioxidant and anticancer effects of ginsenoside re-lysine mixture by Mail lard reaction. Food Chem. 2013;138:876–83. 10.1016/j.foodchem.2012.12.00423411191 10.1016/j.foodchem.2012.12.004
43. Guo L-P Zhou L-Y Kang C-Z Wang H-Y Zhang W-J Wang S Strategies for medicinal plants adapting environmental stress andsimulative habitat cultivationof Dao-Di herbs China J Chin Materia Med 2020 45 1969 74 10.19540/j.cnki.cjcmm.20200302.101
Guo L-P, Zhou L-Y, Kang C-Z, Wang H-Y, Zhang W-J, Wang S, et al. Strategies for medicinal plants adapting environmental stress andsimulative habitat cultivationof Dao-Di herbs. China J Chin Materia Med. 2020;45:1969–74. 10.19540/j.cnki.cjcmm.20200302.10110.19540/j.cnki.cjcmm.20200302.101
44. Tewari RK Paek K-Y Salicylic acid-induced nitric oxide and ROS generation stimulate ginsenoside accumulation in Panax ginseng roots J Plant Growth Regul 2011 30 396 404 10.1007/s00344-011-9202-3
Tewari RK, Paek K-Y. Salicylic acid-induced nitric oxide and ROS generation stimulate ginsenoside accumulation in Panax ginseng roots. J Plant Growth Regul. 2011;30:396–404. 10.1007/s00344-011-9202-310.1007/s00344-011-9202-3
45. Chen F Xie Y Jia Q Li S Li S Shen N Effects of the continuous cropping and soilborne diseases of Panax ginseng C. A. Meyer on rhizosphere soil physicochemical properties, enzyme activities, and microbial communities Agronomy 2023 13 210 10.3390/agronomy13010210
Chen F, Xie Y, Jia Q, Li S, Li S, Shen N, et al. Effects of the continuous cropping and soilborne diseases of Panax ginseng C. A. Meyer on rhizosphere soil physicochemical properties, enzyme activities, and microbial communities. Agronomy. 2023;13:210. 10.3390/agronomy1301021010.3390/agronomy13010210
46. Meng X Huang X Li Q Wang E Chen C Application of UPLC-Orbitrap-HRMS targeted metabolomics in screening of allelochemicals and model plants of ginseng J Plant Physiol 2023 285 153996 10.1016/j.jplph.2023.153996 37141674
Meng X, Huang X, Li Q, Wang E, Chen C. Application of UPLC-Orbitrap-HRMS targeted metabolomics in screening of allelochemicals and model plants of ginseng. J Plant Physiol. 2023;285:153996. 10.1016/j.jplph.2023.15399637141674 10.1016/j.jplph.2023.153996
47. Zhang X Zhang G Tian L Huang L Ecological regulation network of quality in American ginseng: insights from macroscopic-mesoscopic-microscopic perspectives Ind Crops Prod 2023 206 117617 10.1016/j.indcrop.2023.117617
Zhang X, Zhang G, Tian L, Huang L. Ecological regulation network of quality in American ginseng: insights from macroscopic-mesoscopic-microscopic perspectives. Ind Crops Prod. 2023;206:117617. 10.1016/j.indcrop.2023.11761710.1016/j.indcrop.2023.117617
48. Sun J Luo H Yu Q Kou B Jiang Y Weng L Optimal NPK fertilizer combination increases Panax ginseng yield and quality and affects diversity and structure of rhizosphere fungal communities Front Microbiol 2022 13 919434 10.3389/fmicb.2022.919434 35801112
Sun J, Luo H, Yu Q, Kou B, Jiang Y, Weng L, et al. Optimal NPK fertilizer combination increases Panax ginseng yield and quality and affects diversity and structure of rhizosphere fungal communities. Front Microbiol. 2022;13:919434. 10.3389/fmicb.2022.91943435801112 10.3389/fmicb.2022.919434
49. Shen J Yuan L Zhang J Li H Bai Z Chen X Phosphorus dynamics: from soil to plant Plant Physiol 2011 156 997 1005 10.1104/pp.111.175232 21571668
Shen J, Yuan L, Zhang J, Li H, Bai Z, Chen X, et al. Phosphorus dynamics: from soil to plant. Plant Physiol. 2011;156:997–1005. 10.1104/pp.111.17523221571668 10.1104/pp.111.175232
50. Luo B Ma P Nie Z Zhang X He X Ding X Metabolite profiling and genome-wide association studies reveal response mechanisms of phosphorus deficiency in maize seedling Plant J 2019 97 947 69 10.1111/tpj.14160 30472798
Luo B, Ma P, Nie Z, Zhang X, He X, Ding X, et al. Metabolite profiling and genome-wide association studies reveal response mechanisms of phosphorus deficiency in maize seedling. Plant J. 2019;97:947–69. 10.1111/tpj.1416030472798 10.1111/tpj.14160
51. Pal PK Kumar R Guleria V Mahajan M Prasad R Pathania V Crop-ecology and nutritional variability influence growth and secondary metabolites of Stevia rebaudiana Bertoni BMC Plant Biol 2015 15 67 10.1186/s12870-015-0457-x 25849326
Pal PK, Kumar R, Guleria V, Mahajan M, Prasad R, Pathania V, et al. Crop-ecology and nutritional variability influence growth and secondary metabolites of Stevia rebaudiana Bertoni. BMC Plant Biol. 2015;15:67. 10.1186/s12870-015-0457-x25849326 10.1186/s12870-015-0457-x
52. Pant B-D Pant P Erban A Huhman D Kopka J Scheible W-R Identification of primary and secondary metabolites with phosphorus status-dependent abundance in Arabidopsis, and of the transcription factor PHR1 as a major regulator of metabolic changes during phosphorus limitation Plant Cell Environ 2015 38 172 87 10.1111/pce.12378 24894834
Pant B-D, Pant P, Erban A, Huhman D, Kopka J, Scheible W-R. Identification of primary and secondary metabolites with phosphorus status-dependent abundance in Arabidopsis, and of the transcription factor PHR1 as a major regulator of metabolic changes during phosphorus limitation. Plant Cell Environ. 2015;38:172–87. 10.1111/pce.1237824894834 10.1111/pce.12378
53. Kim N-H Jayakodi M Lee S-C Choi B-S Jang W Lee J Genome and evolution of the shade-requiring medicinal herb Panax ginseng Plant Biotechnol J 2018 16 1904 17 10.1111/pbi.12926 29604169
Kim N-H, Jayakodi M, Lee S-C, Choi B-S, Jang W, Lee J, et al. Genome and evolution of the shade-requiring medicinal herb Panax ginseng. Plant Biotechnol J. 2018;16:1904–17. 10.1111/pbi.1292629604169 10.1111/pbi.12926
54. Hao M Zhou Y Zhou J Zhang M Yan K Jiang S Cold-induced ginsenosides accumulation is associated with the alteration in DNA methylation and relative gene expression in perennial American ginseng (Panax quinquefolius L.) along with its plant growth and development process J Ginseng Res 2020 44 747 55 10.1016/j.jgr.2019.06.006 32913404
Hao M, Zhou Y, Zhou J, Zhang M, Yan K, Jiang S, et al. Cold-induced ginsenosides accumulation is associated with the alteration in DNA methylation and relative gene expression in perennial American ginseng (Panax quinquefolius L.) along with its plant growth and development process. J Ginseng Res. 2020;44:747–55. 10.1016/j.jgr.2019.06.00632913404 10.1016/j.jgr.2019.06.006
55. Liang H Kong Y Chen W Wang X Jia Z Dai Y The quality of wildSalvia miltiorrhizafrom Dao Di area in China and its correlation with soil parameters and climate factors Phytochem Anal 2021 32 318 25 10.1002/pca.2978 32761717
Liang H, Kong Y, Chen W, Wang X, Jia Z, Dai Y, et al. The quality of wildSalvia miltiorrhizafrom Dao Di area in China and its correlation with soil parameters and climate factors. Phytochem Anal. 2021;32:318–25. 10.1002/pca.297832761717 10.1002/pca.2978
56. Allevato DM Kiyota E Mazzafera P Nixon KC Ecometabolomic analysis of wild populations of Pilocarpus Pennatifolius (rutaceae) using unimodal analyses Front Plant Sci 2019 10 258 10.3389/fpls.2019.00258 30894869
Allevato DM, Kiyota E, Mazzafera P, Nixon KC. Ecometabolomic analysis of wild populations of Pilocarpus Pennatifolius (rutaceae) using unimodal analyses. Front Plant Sci. 2019;10:258. 10.3389/fpls.2019.0025830894869 10.3389/fpls.2019.00258
57. Huang E Chen Y Fang M Zheng Y Yu S Environmental drivers of plant distributions at global and regional scales Glob Ecol Biogeogr 2021 30 697 709 10.1111/geb.13251
Huang E, Chen Y, Fang M, Zheng Y, Yu S. Environmental drivers of plant distributions at global and regional scales. Glob Ecol Biogeogr. 2021;30:697–709. 10.1111/geb.1325110.1111/geb.13251
