Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Factors influencing cultivated ginseng (Panax ginseng C. A. Meyer) bioactive compounds

  • Han Yu,

    Roles Investigation, Methodology, Writing – original draft

    Affiliation College of Agriculture, Jilin Agricultural University, Changchun, Jilin, China

  • Jiaxin Zhao,

    Roles Writing – review & editing

    Affiliations National & Local United Engineering Laboratory for Chinese Herbal Medicine Breeding and Cultivation, Jilin University, Changchun, Jilin, China, School of Life Sciences, Jilin University, Changchun, Jilin, China

  • Jian You,

    Roles Formal analysis, Project administration

    Affiliations National & Local United Engineering Laboratory for Chinese Herbal Medicine Breeding and Cultivation, Jilin University, Changchun, Jilin, China, School of Life Sciences, Jilin University, Changchun, Jilin, China

  • Jiangnan Li,

    Roles Data curation

    Affiliations National & Local United Engineering Laboratory for Chinese Herbal Medicine Breeding and Cultivation, Jilin University, Changchun, Jilin, China, School of Life Sciences, Jilin University, Changchun, Jilin, China

  • Hongyu Ma,

    Roles Data curation, Software

    Affiliation Jilin Provincial Joint Key Laboratory of Changbai Mountain Biocoenosis and Biodiversity, Academy of Science of Changbai Mountain, Yanbian, Jilin, China

  • Xia Chen

    Roles Conceptualization

    chenxiajlu@163.com

    Affiliations National & Local United Engineering Laboratory for Chinese Herbal Medicine Breeding and Cultivation, Jilin University, Changchun, Jilin, China, School of Life Sciences, Jilin University, Changchun, Jilin, China

Abstract

We aimed to investigate the effects of genome, age, and soil factors on cultivated Panax ginseng C. A. Meyer (CPG) compounds under identical climate and agronomic practices. Eight populations of CPG from different years and rhizosphere soils were collected from garden and cropland in the city of Ji’an, China. Inter-simple sequence repeat (ISSR) primers were used to detect genetic diversity and identity, and soil microbial community diversity. Soil enzyme activities and nutrients were also measured. The contents of total ginsenosides (TG), Rg1, Re, Rf, Rd, and ginsenoside extractions of CPG were analyzed by spectrophotometry and HPLC. The relative importance of each factor was analyzed by mathematical methods such as correlation analysis, stepwise line regression, and path analysis. Regression equations of similarity values of HPLC fingerprint (SVHF), richness index of HPLC fingerprint (RIHF) and the TG, Rg1, Re, Rf, and Rd contents with their respective significant correlation factors were obtained. For SVHF, the relative importance is age>microbial community diversity>genetic diversity. For RIHF, the relative importance is age>genetic diversity>microbial community diversity. For TG, Rg1, and Rf contents, the relative importance is age>microbial community diversity. Ginseng age and genetic identity influenced Rd content, and age was more important. Total phosphorus was the only directly negative effect on Re. According to regression equations and path analysis, increasing age and decreasing Shannon (H') could improve the TG, Rg1, and Rf contents, with little effect on SVHF. Adding age, genetic diversity, and decreasing Shannon (H’) increased RIHF. Adding age and genetic identity could also improve Rd content. Appropriate decreases in total phosphorus might increase Re content. These findings are significant for CPG scientific cultivation methods, through which CPG bioactive ingredients could be finely controlled via regulation of genotypes and cultural conditions.

Introduction

For at least 2,000 years, Panax ginseng C. A. Meyer, a perennial herb in the Araliaceae family commonly known as Asian ginseng [1], has been valued as an herbal tonic and stimulant in China [2,3]. P. ginseng is widely cultivated in northeast China, Japan, Russia, and the Korean peninsula [4]. There are two main cultivated types, garden ginseng (GGS) and cropland ginseng (CGS). GGS is grown by traditional cultivation methods by sowing P. ginseng seeds into a garden after deforestation and reclamation. Under purely artificial conditions, their growth usually spans 4–7 years. CGS includes only sowing seeds of P. ginseng into cropland, and its cultivation techniques are the same as GGS. Ji’an is located in southeast Jilin Province, China, and its climate data are shown in S1 Table. The region’s climate and soil permeability are suitable for the growth and development of ginseng.

The major bioactive ingredients of P. ginseng are a group of triterpene saponins known as ginsenosides [5,6,7]. More than 30 ginsenosides have been isolated from ginseng roots and are classified into two main groups, the glycosides of 20(S)-protopanaxadiol (Rb1, Rb2, Rc, Rd, Rg3, and Rh2) and the glycosides of 20(S)-protopanaxatriol (Re, Rf, Rg1, Rg2, Rh1, and R1) [8,9]. Ginsenosides have extensive pharmacological action including neuroprotective [10,11,12], anti-aging [13,14], immunomodulatory [15,16], cardiovascular protective [17,18], anti-tumor [19,20,21], and internal secretion adjustment effects [22,23]. Different ginsenosides can have completely different biological activities and pharmacological effects [24,25,26]. For example, Rg1 has a role in angiogenesis, while Rb1 inhibits the earliest step of angiogenesis [27]. The composition and content of ginsenosides is the most important factor affecting the ginseng medicinal value. However, TG content in different ginseng roots can vary by up to 20% [28]. Assessments of published literature reveals a poor understanding of the factors influencing composition and content of ginsenosides in ginseng roots, including age, genotype, soil factors, cultivation methods, and preservation or extraction methods [8,29]. If we understand the relative contribution of genotype, age, and soil factors to the variations in cultivated Panax ginseng C. A. Meyer (CPG) total ginsenosides (TG) composition, then scientific cultivation methods could be established.

The objective of our study was to find a quantitative relationship between genome, ginseng age, soil factors, and the composition and content of ginsenosides in CPG roots under the same climate and agronomic practices. We focused on Ji’an, where CPG is reputed to be produced at high quality and sold at premium prices. In addition, the experimental results will contribute to establishing best scientific cultivation methods to improve and control the quality and yield of ginseng roots.

Materials and methods

Plant materials

A total of 126 plants, which corresponded to eight cultivated populations of P. ginseng (CGS: 4, GGS: 4), were taken from Taishang in Ji’an, Jilin province, China in 07/2011 (Table 1). We confirmed that permits were obtained from Yisheng Pharmaceutical Company where collecting took place. We also confirmed that the location accessed was not privately owned and the field studies did not involve endangered or protected species. Fresh leaves were collected, dried in plastic bags with silica gel, transported back to the laboratory, and kept at −80°C. At the same time, the soil adhered to the surface of the roots (rhizosphere soil) was collected and put in sterile polyethylene bags, transported back to laboratory, and kept at −20°C. Within one day of root collection, roots were rinsed with tap water to remove soil, blotted dry, and then dried in plastic bags with silica gel. After drying, the whole roots (containing secondary roots and storage roots) of each population were prepared for analysis by grinding to a fine powder with a tissue grinder (KX-11A/B/C, Ji’nan Kexiang Instrument Co., Ltd., China). Powdered samples were stored at room temperature in plastic bags.

thumbnail
Table 1. Cultivated type, age, and sample sizes of cultivated P. ginseng populations.

https://doi.org/10.1371/journal.pone.0223763.t001

DNA extraction

Total genomic DNA was extracted from leaves by using Plant Genomic DNA Isolation Kit (NEP003-1, Beijing Dingguo Changsheng Biotechnology Co., Ltd., China). DNA concentration was then determined by comparing the plant DNA samples with commercial standard lambda DNA on 0.8% (w/v) agarose gel, after which it was adjusted to 5 ng/μl.

ISSR-PCR amplification

ISSR primers used in this study were synthesized by Beijing Dingguo Changsheng Biotechnology Co., Ltd (China), according to the primer set published by the University of British Columbia (UBC). One hundred ISSR primers were initially screened, and twelve that yielded bright and discernible bands, were used for the analysis of all 126 samples (Table 1). Fifteen or sixteen individuals from each population were used for the primer screening, and PCR amplifications were repeated for working primers to check the stability and reproducibility of ISSR fragments. PCR was performed in 25 μl reactions containing 1.75 mM MgCl2, 0.25 mM dNTPs, 1 U Taq DNA polymerase (TaKaRa), 0.2 μM primers and 10 ng DNA templates. PCR amplifications were performed in the Mastercycler Gradient PCR (Eppendorf, Germany) with the following program: initial denaturation at 94°C for 5 min; 40 cycles of 94°C for 50 s, appropriate annealing temperature (see Table 2) for 45 s, 72°C for 1 min; and final synthesis at 72°C for 10 min. A negative control with no DNA added was included in each PCR run. Amplification products were separated with 1.5% agarose gels (1×TAE buffer) at 80 V for 1.5 h, stained with ethidium bromide (0.5 μg/ml), and photographed under UV light using an EC3 Gel Documentation System (UVP, USA).

thumbnail
Table 2. Polymorphisms of inter-simple sequence repeat markers in cultivated P. ginseng populations.

https://doi.org/10.1371/journal.pone.0223763.t002

Extraction of ginsenosides

The extraction method of ginsenosides as based on the protocol by Lim W et al., with some modifications [30]. An accurately weighed sample (100 mg) of each population’s roots was transferred to a 50 ml centrifuge tube. Ginsenosides were extracted in 20 ml of 100% HPLC-grade methanol and placed in a sonicator bath for 15 min at 60°C. The sample tube was centrifuged at 5,625 g for 10 min, and the supernatant was collected. The precipitate was re-extracted two additional times with 20 mL of solvent each time, and the supernatants were combined. The supernatant was reduced to dryness under vacuum with a rotary evaporator at 38°C, and the residue was re-dissolved in 2 mL of 100% methanol. This was dried under a stream of N2 at 38°C and re-dissolved in 500 μl of 70% (v/v) HPLC-grade methanol diluted with HPLC-grade water. Samples were re-filtered and 15 μl of extract was immediately injected in the HPLC system.

HPLC analysis and TG determination

A HP1100 high-performance liquid chromatography (HPLC) system was used (Agilent Technologies Inc., Palo Alto, CA) with gradient elution and a μBondapak C18 reversed phase column (10 μm, 4.6 mm×150 mm) (Waters Inc., Milford, MA). The binary gradient employed the mobile phases: (A) phosphate buffer (10.3 mM KH2PO4 at pH 5.8) and (B) CH3CN with a flow rate of 1.2 ml/min according to the following profile adapted from Lim W et al. [30]: 0–20 min, 84–82% A and 16–18% B; 20–60 min, 82–60% A and 18–40% B, 60–120 min 60%–5% A and 40%–95% B. The UV diode array detector was set at 203 nm. Ginsenoside standards included Rg1, Re, Rf, and Rd (National Institutes for Food and Drug Control, NIFDC). Qualitative identification of ginsenoside peaks was determined by cochromatography (equivalent retention time) with chemically pure standards, and quantification was based on the integration of the peak area compared with a standard curve. Results are reported as percent ginsenoside on a dry weight basis.

The spectrophotometric method was used to determine the TG content (mg/g) present in each population’s roots [31]. Each sample extract (50 μl) was diluted to 0.5 ml methanol and reacted at 60°C for 10 min with 8% vanillin solution (0.5 ml) and 87% sulfuric acid (5 ml). The absorbance of the reaction mixture was read at 544 nm against a blank solution.

Microbial flora analyses

Viable total counts of cultivable bacteria, fungi, actinomycetes, abiogenous Azotobacter, cellulose-decomposing microorganisms, nitrifying bacteria, sulfur bacteria, ammonifying bacteria, and potassium bacteria were determined as colony forming units (CFUs) on agar plates by dilution plate methods. The medium used for the enumeration of bacteria, fungi, actinomycetes, abiogenous Azotobacter, cellulose-decomposing microorganisms, nitrifying bacteria, sulfur bacteria, ammonifying bacteria, and potassium bacteria were beef extract peptone medium, improved Gause’s No.1 medium, Rose bengal medium, Ashby nitrogen free medium, cellulose-Congo red medium, nitrifying bacteria medium, sulfur bacteria medium, peptone ammonifying culture medium, and potassium aluminum silicate agar medium, respectively [32].

Determination of rhizosphere soil

The activities of sucrase, urease, acid phosphatase, catalase and cellulase in P. ginseng rhizosphere soil were determined according to Guan [33]. Chemical analyses (total nitrogen, total phosphorus, total potassium, nitrate nitrogen, ammonium nitrogen, available phosphorus, available potassium, and organic matter) were done according to analysis of soil physical and chemical properties [34].

Data analysis

Amplified bands were scored 1/0 as presence/absence of homologous bands for all samples. The presence/absence data matrix was analyzed using POPGENE version 1.32 [35,36,37] to calculate various genetic diversity parameters, including the percentage of polymorphic loci (Ppl), Shannon’s information index (I) and genetic diversity (h), genetic diversity, gene differentiation coefficient (Gst) and gene flow (Nm), and total genetic diversity (Ht) and within group genetic diversity (Hs). Genetic distance was also generated by POPGENE and a dendrogram was constructed from Nei’s (1978) genetic distance with the unweighted pair-group method of averages (UPGMA) with 1,000 permutations of bootstrapping using MEGA v5.2. SVHFs were computed by the professional software Similarity Evaluation System for Chromatographic Fingerprint of Traditional Chinese Medicine (Version 2004 A), which was developed and recommended by Chinese State Food and Drug Administration. This software was also used to synchronize among different samples [38,39]. DPS 14.10 (data processing system) was employed to compute the correlation of SVHF, RIHF, and the contents of Rg1, Re, Rf, and Rd with age, genetic diversity, genetic identity, Shannon (H’) and soil nutrients, stepwise line regression, and path analysis [40].

Results and discussion

ISSR profile and genetic analysis

The twelve selected ISSR primers generated 1856 clear and repeatable DNA fragments from eight CPG populations. The amplified DNA fragments ranged from 180 to 2,200 bp in size. DNA fragments of the same size were considered as the same band. In total, 232 ISSR bands were detected with repeatability across 126 P. ginsengs from eight cultivated populations. The number of bands per primer varied between 14 (UBC836) and 26 (UBC807), with an average of 23.2 (Table 1). Four of 12 primers revealed ISSR loci with 100% polymorphism at the species level, while other primers detected polymorphic loci from 76.9% (primer UBC815) to 94.7% (primer UBC823), leading to an average of 21.2 polymorphic loci per primer (Table 2).

A high level of genetic variation was detected using ISSR markers, with 91.38% polymorphic loci at the species level. The CGS IV population had the highest diversity (h = 0.1749, I = 0.2595, and Ppl = 49.57), while the CGS I population shown the lowest diversity (h = 0.0938, I = 0.1409, and Ppl = 28.88%) (Table 3). This study revealed that the species-level genetic diversity (Ppl = 91.83%, h = 0.2454, I = 0.3823) in GGS and CGS was higher than that in its cultivated conspecifics (Ppl = 85.42%, h = 0.2294, I = 0.3590) or its cultivated congeneric counterparts, e.g. P. quinquefolius L. (RAPD: Ppl = 45.7%; Allozyme: Ppl = 62.5%) and P. notoginseng (RAPD: Npl = 75.5%), and approximated its wild conspecifics (AFLP: Npl = 94.4%, h = 0.3246) [4,41,42,43,44]. Therefore, genetic diversity in the eight populations (CGS and GGS) selected in this study could represent CPG genetic diversity.

thumbnail
Table 3. Genetic diversity within populations of cultivated P. ginseng based on inter-simple sequence repeat data.

https://doi.org/10.1371/journal.pone.0223763.t003

At the species level, the coefficient of gene differentiation (Gst) was 0.4551, and the limited into population gene flow (Nm) was 0.5987 (Table 4). The estimate of the total genetic diversity (Ht) was 0.2463, and the within group genetic diversity (Hs) was 0.1342, indicating that the total genetic diversity in this species (about 55.5%) was primarily from genetic divergence between horticultural P. ginseng populations. This result indicates that the high genetic diversity in CGS and GGS could be attributed to the dominance of selfing (ranging from 58.14% to 89%) in P. ginseng [45,46]. Therefore, the genetic identity and diversity of CPG populations are relatively stable and the interference by other populations is relatively small.

thumbnail
Table 4. The coefficient of gene differentiation and gene flow at the species level of cultivated P. ginseng.

https://doi.org/10.1371/journal.pone.0223763.t004

Nei’s (1978) genetic distances ranged from 0.0903 (GGS IV vs. CGS III) to 0.2003 (GGS IV vs. CGS II), with an average of 0.1521 (Table 5). Accordingly, the genetic identity ranged from 0.8001 (CGS I vs. CGS II) to 0.9137 (GGS IV vs. CGS III). The genetic identity (from 0.8416 to 0.8997) of GGS was more uniform than CGS (from 0.8001 to 0.9023). The UPGMA cluster analysis clustered all eight cultivated populations into four groups (Fig 1), rather than the eight cultivated populations attached to two cultivated groups (CGS and GGS). In other words, all populations that belonged to the same cultivated type (GGS and CGS) were not clustered together. This is consistent with a randomly chosen P. ginseng seed when sown.

thumbnail
Fig 1. UPGMA dendrogram of garden ginseng (GGS) and cropland ginseng (CGS).

UPGMA dendrogram based on 232 inter-simple sequence repeat (ISSR) loci showing the genetic relationships among eight populations of the two main cultivated P. ginseng (CGS and GGS).

https://doi.org/10.1371/journal.pone.0223763.g001

thumbnail
Table 5. Nei's unbiased measures of genetic identity and genetic distance in garden ginseng (GGS) and cropland ginseng (CGS).

https://doi.org/10.1371/journal.pone.0223763.t005

HPLC fingerprint analysis and ginsenoside content

Standard solutions of Rg1, Re, Rf, and Rd were prepared in 70% (v/v) HPLC-grade methanol diluted with HPLC-grade water at final concentrations of 0.03, 0.06, 0.13, 0.25, 0.50, and 1.00 mg/mL, respectively. Calibration was performed by analyzing the four reference solutions in duplicate at six concentration levels, and then the calibration curves were constructed by plotting the peak areas versus the injection concentrations of each compound.

HPLC fingerprints obtained from eight batches of eight CPG populations and the CGS IV HPLC reference fingerprint are given in Fig 2. The SVHFs versus reference fingerprint are tabulated in Table 6. There were 21 common peaks in all eight batches, and common peaks area accounted for over 52% of the overall peaks area. Common peak area increased with age in the GGS and CGS. Peaks 12, 13, 15, and 20 were identified as Rg1, Re, Rf, and Rd by comparison with the corresponding chemical references chromatogram under the same conditions (Fig 3). Each sample was analyzed in duplicate to determine the mean contents (mg/g) of TG and four selected ginsenosides. The results are shown in Table 6.

thumbnail
Fig 2. HPLC fingerprints of garden ginseng (GGS) and cropland ginseng (CGS).

HPLC fingerprints obtained from batches of eight cultivated P. ginseng populations. The reference fingerprint was defined as the CGS IV HPLC fingerprint. 1–21 correspond to 21 common peaks, while peaks 1–10 were not intense and not marked. Peaks 12, 13, 15, and 20 were identified as Rg1, Re, Rf, and Rd, respectively.

https://doi.org/10.1371/journal.pone.0223763.g002

thumbnail
Fig 3. Chromatogram of mixed standard compounds (Rg1, Re, Rf, and Rd).

https://doi.org/10.1371/journal.pone.0223763.g003

thumbnail
Table 6. HPLC fingerprint analysis and total ginsenoside (TG) content of garden ginseng (GGS) and cropland ginseng (CGS).

https://doi.org/10.1371/journal.pone.0223763.t006

RIHF is calculated by the Monk (1967) index with the formula R = S/N, where S is the number of HPLC peaks for each sample, and N is the number of HPLC peaks for all samples (common peaks were only counted once) [47]. The RIHF indicates the rich degree of chemical components in CPG roots. The computational results are shown Table 6.

The effect of age was not the same for SVHF, RIHF, TG, Rg1, Re, Rf, and Rd. The TG, Rg1, Rf, and Rd contents increased with increasing age in GGS and CGS. Ginseng age was more approximate and SVHF was higher, indicating that the main chemical composition (the peak accounting for over 5% of the total peak area of the peak) was more similar [48,49]. RIHF increased from age I to IV in GGS and CGS, indicating that as ginseng age increased, chemical components were enriched in CPG root. The content of Re also increased with age from II to IV in GGS, but not in CGS, indicating that other factors could affect the content of Re. In general, with the increase of ginseng age, the content and number of CPG root ginsenosides was greater and CGP bioactive value was much better.

Rhizosphere soil microbial community diversity analysis

Fig 4 shows the amount of bacteria (BA), fungi (FU), actinomycetes (AC), abiogenous Azotobacter (AA), cellulose-decomposing microorganisms (CM), nitrifying bacteria (NB), sulfur bacteria (SB), ammonifying bacteria (AB), and potassium bacteria (PB) in P. ginseng rhizosphere soil. The amount of these microorganisms in each age-matched rhizosphere samples was significantly different between CGS and GGS (p<0.01). The microorganisms content in GGS rhizosphere soil was more than five times higher that of the CGS rhizosphere. Diversity and evenness were lower in the CGS than those in GGS, and increased with CGS and GGS age (Table 7). Similar results were also reported by Yong Li et al. and Li Xi-ying et al. [50,51]. Changes in microbial community diversity could be induced by environmental factors, such as overuse of nitrogen, phosphorus fertilizers, and exudates released from roots to their adjacent soil [52,53].

thumbnail
Fig 4. Amounts of rhizosphere soil microorganisms.

Bacteria (BA), Fungi (FU), Actinomycetes (AC), Abiogenous Azotobacter (AA), Cellulose-decomposing Microorganisms (CM), Nitrifying Bacteria (NB), Sulfur Bacteria (SB), Ammonifying Bacteria (AB), and Potassium Bacteria (PB) were counted. Eight samples from each of three replicate plots were averaged. Values are mean±standard deviation.

https://doi.org/10.1371/journal.pone.0223763.g004

thumbnail
Table 7. Culturable microbial community diversity indices for garden ginseng (GGS) and cropland ginseng (CGS).

https://doi.org/10.1371/journal.pone.0223763.t007

Rhizosphere soil enzymatic activities

Soil enzymes play an important role in the material cycle and energy transformation of soil ecological systems. They are also important for catalyzing reactions necessary for the life of microorganisms and plants, decomposition of organic residues, cycling of nutrients, and formation of organic matter and soil structure [37,54]. Soil enzyme activities may be considered early and sensitive indicators to measuring the degree of soil degradation in both natural and agro-ecosystems, and can be an important indicator of soil fertility [37,55,56,57]. Table 8 shows the rhizosphere soil enzymatic activities of eight populations of GGS and CGS. The enzymatic activities of sucrase, acid phosphatase, and cellulase firstly increased from age I to II, and then decreased from age II to III, and finally increased from age III to IV. In contrast, the urease activity firstly decreased from age I to II, and then increased from age II to III, and finaly decreased from age III to IV. Catalase activity increased with CGS and GGS age. The GGS enzymatic activities of sucrase, acid phosphatase catalase and cellulase were higher than their peers, while the urease activity was the opposite. Therefore, the soil fertility in GGS was higher than in CGS.

thumbnail
Table 8. Activities of sucrase, urease, acid phosphatase, catalase, and cellulase in the rhizosphere soil of garden ginseng (GGS) and cropland ginseng (CGS).

https://doi.org/10.1371/journal.pone.0223763.t008

Rhizosphere soil nutrient analysis

Table 9 shows that the GGS rhizosphere soil exhibited the same available P as CGS, and 2–3 fold greater soil total N, nitrate N, ammonium N, available K, and organic matter than CGS. However, total P and total K in GGS were only approximately 70% as much as in CGS.

thumbnail
Table 9. Rhizosphere soil nutrient analysis for eight populations of garden ginseng (GGS) and cropland ginseng (CGS).

https://doi.org/10.1371/journal.pone.0223763.t009

Statistical analysis

The correlation coefficient of SVHF, RIHF, and the content of TG, Rg1, Re, Rf, and Rd with rhizosphere soil enzymatic activities and nutrients was not significant (p>0.05) (Table 10). The correlation coefficient of the content of TG and Rg1 with age and Shannon (H’) was significant (p<0.05) and positive. The correlation coefficient of SVHF, RIHF, and the content of Rf with age, genetic diversity index (h), and Shannon (H') was significant (p<0.05) and positive. The correlation coefficient of Re content with total phosphorus was very significant (p<0.01) and negative. The age and genetic identity were significantly related to the content of Rd (p<0.05). Thus, selection for these significant correlative factors may improve SVHF, RIHFs, TG, and selected four-monomer ginsenoside content (Rg1, Re, Rf, and Rd). The stepwise line regression equations of SVHF, RIHF, TG, and selected four-monomer ginsenoside content as dependent variable (Y) with their own significant correlative factors (Xn) are shown in Table 11. Path analysis results of SVHF, RIHF, TG, and four-monomer ginsenoside content to their own significant correlative factors are shown in Table 12. Ginseng age was the most important influence on SVHF, RIHF, and TG, Rg1, Rf, and Rd contents, but had no effect on Re content. Ginseng age directly affected TG (0.8187), Rg1 (0.9869), and Rf (0.9996), and was higher than their respective Shannon (H') (microbial community diversity index) (0.043, −0.1080, and −0.1019). Shannon (H') had a direct positive effect on TG content, but a direct negative effect on Rg1 and Rf contents. Determination coefficients of TG, Rg1, and Rf contents were 0.7278, 0.8158, and 0.8486, respectively. This means that the theoretical values calculated through the regression equations of TG, Rg1, and Rf contents had considered their respective variability of 72.78%, 81.58%, and 84.86%. The first direct effect to SVHF was age (0.8336), second was Shannon (H') (0.1170), and the third was h (0.0441). The first direct effect to RIHF was age (0.8341), second was h (0.2874), and the third was Shannon (H') (−0.1519). The determination coefficients of SVHF and RIHF were 0.9291 and 0.9224, respectively. This means that the theoretical values of SVHF and RIHF had considered their respective variability of 92.91% and 92.24%. Therefore, ginseng age, h, and Shannon (H') are more approximate and SVHF was higher, indicating that the main chemical composition is more similar. The same ginseng age, h, and Shannon (H’) might produce a similar amount of chemical components in CPG roots. Therefore, increasing age and genetic diversity while reducing microbial community diversity could increase the number of chemical components. For one ginseng population, appropriate selection for age and Shannon (H') could result in increasing RIHF, TG, Rg1, and Rf contents, but had little effect on SVHF. The age (0.5323) direct effect on Rd was higher than genetic identity (0.4511), and its determination coefficient was 0.7330. This indicates that improving age and genetic consistency could increase the content of Rd. Total phosphorus (−0.9249) was a directly negative effect on Re content, and its determination coefficient is 0.8555. Thus adding phosphate fertilizer could decrease the content of Re. This result is the same as Konsler T R. at el. [58].

thumbnail
Table 10. Simple correlation of similarity values of HPLC fingerprint (SVHF), richness values of HPLC fingerprint (RIHF), and the contents of total ginsenosides (TG), Rg1, Re, Rf, and Rd with age, genetic diversity, genetic identify, Shannon (H’), and soil nutrients.

https://doi.org/10.1371/journal.pone.0223763.t010

thumbnail
Table 11. Regression equations of similarity values of HPLC fingerprint (SVHF), richness values of HPLC fingerprint (RIHF), and the contents of total ginsenosides (TG), Rg1, Re, Rf, and Rd with their own significance factors.

https://doi.org/10.1371/journal.pone.0223763.t011

thumbnail
Table 12. Path analyses of similarity values of HPLC fingerprint (SVHF), richness values of HPLC fingerprint (RIHF), and the contents of total ginsenosides (TG), Rg1, Re, Rf, and Rd with their own significance factors.

https://doi.org/10.1371/journal.pone.0223763.t012

During cultivation of P. ginseng, appropriate selection for various factors could improve SVFH, RIHF, and the content of ginsenosides (TG, Rg1, Re, Rf, and Rd). For example, appropriate age increases could improve RIHF and TG, Rg1, Rf, and Rd contents in CPG root. However, rhizosphere soil microbial community diversity increased with age (Tables 7 and 13). This could be due to secretions from the ginseng root causing increases in specific carbon substrates and/or signaling compounds supporting increased rhizosphere microbial community diversity [59]. Because increasing rhizosphere soil microbial community diversity could decrease RIHF and TG, Rg1, and Rf contents according to stepwise line regression equations and path analyses, appropriate management measures can be taken to reduce microbial community diversity (Shannon [H']) while managing CPG. The simple correlation of Shannon (H') with age, genetic diversity, genetic identify, and soil nutrients is shown in Table 13. The moderate correlation (0.5≤|correlation coefficient|<0.8) of Shannon (H') with sucrase, acid phosphatase, catalase, total nitrogen, nitrate nitrogen, and ammonium nitrogen was positive, while the correlation of nitrate nitrogen and ammonium nitrogen with acid Phosphatase, catalase and total nitrogen was positive and significant (p<0.05). The simple correlation of nitrate nitrogen and ammonium nitrogen with RIHF and the content of TG, Rg1, and Rf was not correlated (|correlation coefficient|<0.3) or had low correlation (0.3≤|correlation coefficient|<0.5). Thus, appropriate reduction in the amount of ammonium nitrogen and nitrate nitrogen in fields could reduce Shannon (H') and improve RIHF and TG, Rg1, and Rf contents. This conjecture was consistent with published results that root N is negatively correlated with root Rg1 and the accumulation of TG was severely inhibited when NH4+ content is increased [58,60].

thumbnail
Table 13. Simple correlation of Shannon (H') with age, genetic diversity, genetic identify, and soil nutrients.

https://doi.org/10.1371/journal.pone.0223763.t013

In this study, we defined the quantitative relationship between SVHF, RIHF, and the contents of Rg1, Re, Rf, and Rd and their respective significant correlative factors (age, genetic diversity, genetic identify, Shannon [H’], and soil nutrients). These findings could help progress CPG cultivation methods. The regression coefficients of acquired regression equations were less than 0.999 and remaining path coefficients were also larger (>0.2664) (Tables 11 and 12), indicating that some factors influencing SVHF, RIHF, and TG, Rg1, Re, Rf, and Rd contents were not taken into account. These factors might include light, rainfall, moisture, temperature, uncultivable microbial community, soil physical properties, soil chemical properties, soil trace elements (such as Mn, Me, and Zn), and cultural practices. If these factors influencing the content and constituents of CPG could be controlled, an accurate quantitative relationship between chemical content and factors could be determined by mathematical analysis. These accurate quantitative relationships combined with modern networks and automatic detection technology can establish the best CPG cultivation methods. These methods allow us to control the content and constituents of CPG bioactive ingredients by adjusting the related influencing factors. Because the potential benefits of specific ginsenosides on cancer and diabetes has been published [61,62], CPG cultivation methods enhancing the production of specific monomer ginsenosides and other bioactive ingredients could seriously impact commerce in this medicinal herb and its future role in public health.

Conclusions

In conclusion, we obtained the regression equations of similarity values of HPLC fingerprint (SVHF), richness index of HPLC fingerprint (RIHF) and the TG, Rg1, Re, Rf, and Rd contents with their respective significant correlation factors. SVHF and RIHF were influenced not only by age and microbial community diversity but also genetic diversity. For SVHF, the relative importance is age>microbial community diversity>genetic diversity. For RIHF, the relative importance is age>genetic diversity>microbial community diversity. The factors that influence TG, Rg1, and Rf content were ginseng age and microbial community diversity, by contrast, ginseng age was the main influencing factor. Ginseng age and genetic identity influenced Rd content, and age was more important. Re was influenced only by total phosphorus. Therefore, under the same climate, the relative importance of genes, age, and soil factors were not the same for SVHF, RIHF, and TG Rg1, Re, Rf, and Rd contents in CPG. In general, increasing age and decreasing Shannon (H') could improve RIHF and TG, Rg1 Rf, and Rd contents, but had little effect on SVHF; increasing age and genetic diversity identity could also improve the content of Rd; appropriate decreases in total phosphorus might increase the content of Re. These findings can help progress CPG cultivation methods, which could help achieve customized CPG bioactive ingredients through regulating genotypes and cultural conditions.

Supporting information

S1 Table. Basic climatic information of Ji 'an, Jilin Province, China (according to the data from 1971–2000).

https://doi.org/10.1371/journal.pone.0223763.s001

(DOCX)

Acknowledgments

I would to thank Fangliang Li, Dr. Ming Xing and Dr. Wei Zhao for providing logistical support in the laboratory.

References

  1. 1. Chinese Academy of Science. Editorial Committee of the Flora of China: Flora of China. Volume 54. Beijing: Science Press; 1978. pp. 180.
  2. 2. Wen J, Zimmer EA. Phylogeny and biogeography of PanaxL. (the ginseng genus, Araliaceae): inferences from ITS sequences of nuclear ribosomal DNA. Mol Phylogenet Evol. 1996;6(2): 167–177. pmid:8899721
  3. 3. Lee JW, Kim YC, Jo IH, Seo AY, Lee JH, Kim OT, et al. Development of an ISSR-Derived SCAR Marker in Korean Ginseng Cultivars (Panax ginseng C. A. Meyer). J. Ginseng Res. 2011;35(1): 52–59.
  4. 4. Li S, Li J, Yang XL, Cheng Z, Zhang WJ. Genetic diversity and differentiation of cultivated P. ginseng (Panax ginseng C. A. Meyer) populations in North-east China revealed by inter-simple sequence repeat (ISSR). Genet Resour Crop Evol. 2011;58: 815–824.
  5. 5. Wang CZ, Ni M, Sun S, Li XL, He H, Mehendale SR, et al. Detection of adulteration of notoginseng root extract with other panax species by quantitative HPLC coupled with PCA. J Agric Food Chem. 2009;57(6): 2363–2367. pmid:19256509
  6. 6. Zuo G, Guan T, Chen DL, Li CL, Jiang R, Luo CY, et al. Total saponins of Panax ginseng induces K562 cell differentiation by promoting internalization of the erythropoietin receptor. Am J Chin Med. 2009;37(4): 747–57. pmid:19655412
  7. 7. Taira S, Ikeda R, Yokota N, Osaka I, Sakamoto M, Kato M, et al. Mass spectrometric imaging of ginsenosides localization in Panax ginseng root. Am J Chin Med. 2010;38(3): 485–93. pmid:20503467
  8. 8. Schlag EM, Mclntosh MS. Ginsenoside content and variation among and within American ginseng (Panax quinquefoliusL.) populations. Phytochemistry. 2006;67: 1510–1519. pmid:16839573
  9. 9. Leung KW, Wong AST. Pharmacology of ginsenosides: a literature review. Chinese Medicine. 2010;5(20): 1–7.
  10. 10. Kampen JV, Robertson H, Hagg T, Drobitch R. Neuroprotective actions of the ginseng extract G115 in two rodent models of Parkinson's disease. Experimental neurology. 2003;184(1): 521–529. pmid:14637121
  11. 11. Liao B, Newmark H, Zhou R. Neuroprotective Effects of Ginseng Total Saponin and Ginsenosides Rb1 and Rg1 on Spinal Cord Neurons in Vitro. Experimental neurology. 2002;173(2): 224–234. pmid:11822886
  12. 12. Lim JH, Wen TC, Matsuda S, Tanaka J, Maeda N, Peng H, et al. Protection of ischemic hippocampal neurons by ginsenoside Rb1, a main ingredient of ginseng root. Neuroscience research. 1997;28(3): 191–200. pmid:9237267
  13. 13. Cho J, Park W, Lee S, Ahn W, Lee Y. Ginsenoside-Rb1 from Panax ginseng C.A. Meyer activates estrogen receptor-α and-β, independent of ligand binding. Journal of Clinical Endocrinology & Metabolism. 2004;89(7): 3510–3515.
  14. 14. Cheng Y, Shen LH, Zhang JT. Anti-amnestic and anti-aging effects of ginsenoside Rg1 and Rb1 and its mechanism of action. ACTA pharmacologica sinica. 2005;26(2): 143–149. pmid:15663889
  15. 15. Kenarova B, Neychev H, Hadjiivanova C, Petkov VD. Immunomodulating activity of ginsenoside Rg1 from Panax ginseng. Japanese journal of pharmacology. 1990;54(4): 447–449. pmid:2087006
  16. 16. Liu M, Zhang JT. Studies on the mechanisms of immunoregulatory effects of ginsenoside Rg1 in aged rats. Acta pharmaceutica Sinica. 1995;31(2): 95–100.
  17. 17. Zhou W, Chai H, Lin PH, Lumsden AB, Yao QZ, Chen CY. Ginsenoside Rb1 blocks homocysteine-induced endothelial dysfunction in porcine coronary arteries. Journal of vascular surgery. 2005;41(5): 861–868. pmid:15886672
  18. 18. Zhang QH, Wu CF, Duan L, Yang JY. Protective effects of total saponins from stem and leaf of Panax ginseng against cyclophosphamide-induced genotoxicity and apoptosis in mouse bone marrow cells and peripheral lymphocyte cells. Food and Chemical Toxicology. 2008;46(1): 293–302. pmid:17904265
  19. 19. Liu WK, Xu SX, Che CT. Anti-proliferative effect of ginseng saponins on human prostate cancer cell line. Life sciences. 2000;67(11): 1297–1306. pmid:10972198
  20. 20. Hofseth LJ, Wargovich MJ. Inflammation, cancer, and targets of ginseng. The Journal of nutrition. 2007;137(1): 183–185.
  21. 21. Choi YJ, Lee HJ, Kang DW, Han IH, Choi BK, Cho WH. Ginsenoside Rg3 induces apoptosis in the U87MG human glioblastoma cell line through the MEK signaling pathway and reactive oxygen species. Oncology Reports. 2013;30(3) 1362–1370 pmid:23783960
  22. 22. Ling C, Li Y, Zhu X, Zhang C, Li M. Ginsenosides may reverse the dexamethasone-induced down-regulation of glucocorticoid receptor. General and comparative endocrinology. 2005;140(3): 203–209. pmid:15639148
  23. 23. Nocerino E, Amato M, Izzo AA. The aphrodisiac and adaptogenic properties of ginseng. Fitoterapia. 2000;71: 1–5.
  24. 24. Nah SY, Kim DH, Rhim H. Ginsenosides are any of them candidates for drugs acting on the central nervous system? CNS Drug Rev. 2007;13: 381–404. pmid:18078425
  25. 25. Nah SY, Bhatia KS, Lyles J, Ellinwood EH, Lee TH. Effects of ginseng saponin on acute cocaine-induced alterations in evoked dopamine release and uptake in rat brain nucleus accumbens. Brain Res. 2009;1248: 184–90. pmid:19026615
  26. 26. Chen CF, Chiou WF, Zhang JT. Comparison of the pharmacological effects of Panax ginseng and Panax quinquefolium. Acta Pharmacol Sinica. 2008;29: 1103–8.
  27. 27. Sengupta S, Toh SA, Sellers LA, Skepper JN, Koolwijk P, Leunq HW, et al. Modulating angiogenesis: the yin and the yang in ginseng. Circulation. 2004;110: 1219–25. pmid:15337705
  28. 28. Yue PY, Mak NK, Cheng YK, Leung KW, Ng TB, Fan DT, et al. Pharmacogenomics and the Yin/Yang actions ofginseng: anti-tumor, angiomodulating and steroid-like activities of ginsenosides. Chin. Med. 2007;2(6): 1–21.
  29. 29. Davy CW Lee, Allan SY Lau. Effects of Panax ginsengon Tumor Necrosis Factor-α-Mediated Inflammation: A Mini-Review. Molecules. 2011;16: 2802–2816. pmid:21455094
  30. 30. Lim W, Mudge KW, Vermeylen F. Effects of population, age, and cultivation methods on ginsenoside content of wild American ginseng (Panax quinquefolium). Journal of agricultural and food chemistry. 2005;53(22): 8498–8505. pmid:16248544
  31. 31. Lui JH, Staba EJ. The ginsenosides of various ginseng plants and selected products. Journal of Natural Products. 1980;43(3): 340–346.
  32. 32. Xu GH, Zheng HY. Handbook of Analysis Methods of Soil Microbiology. Beijing: Agricultural Press; 1986. pp. 102.
  33. 33. Guan SY. Methods for Determination of Soil Enzyme Activity. Beijing: Agricultural Press;1986. pp. 234.
  34. 34. Nanjing Institute of Soil Science, CAS. Analysis of soil physical and chemical properties. Shanghai: Shanghai Sci.& Tech. Press; 1981. pp. 134.
  35. 35. Yeh FC, Yang RC, Boyle T; 1999. POPGENE, version 1.32: the user friendly software for population genetic analysis. Molecular Biology and Biotechnology Centre, University of Alberta, Edmonton, AB, Canada.
  36. 36. Francis CY, Yang RC; 2000. Popgene version 1.32. Available from: http://www.seekbio.com/DownloadShow.asp?id=1059
  37. 37. Doelman P, Haanstra L. Short-and long-term effects of heavy metals on phosphatase activity in soils: An ecological dose-response model approach. Biology and Fertility of Soils. 1989;8(3): 235–241.
  38. 38. Tistaert C, Dejaegher B, Heyden YV. Chromatographic separation techniques and data handling methods for herbal fingerprints: a review. Analytica chimica acta. 2011;690(2): 148–161. pmid:21435470
  39. 39. Chen Y, Yan Y, Xie MY, Nie SP, Liu W, Gong xf, et al. Development of a chromatographic fingerprint for the chloroform extracts of Ganoderma lucidum by HPLC and LC–MS. Journal of pharmaceutical and biomedical analysis. 2008;47(3): 469–477. pmid:18337046
  40. 40. Tang QY, Feng MG DPS data processing system: experimental design, statistical analysis and data mining. Science, Beijing. 2007
  41. 41. Bai D, Brandle J, Reeleder R. Genetic diversity in North American ginseng (Panax quinquefolius L.) grown in Ontario detected by RAPD analysis. Genome. 1997;40(1): 111–115. pmid:18464811
  42. 42. Cruse-Sanders JM, Hamrick JL. Genetic diversity in harvested and protected populations of wild American ginseng, Panax quinquefolius L.(Araliaceae). American Journal of Botany. 2004;91(4): 540–548. pmid:21653409
  43. 43. Li S, Li J, Yang XL, Cheng Z, Zhang WJ. Genetic diversity and differentiation of cultivated P. ginseng (Panax ginseng CA Meyer) populations in North-east China revealed by inter-simple sequence repeat (ISSR) markers. Genetic Resources and Crop Evolution. 2011;58(6): 815–824.
  44. 44. Reunova GD, Kats IL, Muzarok TI, Zhuravlev YN. Polymorphism of RAPD, ISSR and AFLP markers of the Panax ginseng CA Meyer (Araliaceae) genome. Russian Journal of Genetics. 2010;46(8): 938–947.
  45. 45. Zhao Y, Gu X, Wu L. Researches on categories, characteristics, and utilization value of cultivated P. ginseng germplasm resources. Chinese Traditional and Herbal Drugs. 2007;38(2): 294.
  46. 46. Culley TM, Wallace LE, Gengler-Nowak KM, Crawford DJ. A comparison of two methods of calculating GST, a genetic measure of population differentiation. American Journal of Botany. 2002;89(3): 460–465. pmid:21665642
  47. 47. Monk CD. Tree species diversity in the eastern deciduous forest with particular reference to north central Florida. Am. Natur. 1967;101: 173–87
  48. 48. Nanping Z, Xinyue X, Ping Z. Study on the establishing of reference fingerprint for the traditional Chinese medicine. Chinese Pharmaceutical Affairs. 2003;6: 347–350.
  49. 49. Qinglei S, Yunliang L, He Z. Study on the Similarity Between HPLC Fingerprints. Chemical Analysis and Meterage. 2006;15(6): 54–55.
  50. 50. Li Y, Ying YX, Zhao DY, Ding WL. Microbial community diversity analysis of Panax ginseng rhizosphere and non-rhizosphere soil using randomly amplified polymorphic DNA method. Open Journal of Genetics. 2012;2: 95–102.
  51. 51. Li XY, Jin HQ, Jia B. Biodiversity of soil microorganisms in the fields of different growing years of Ginseng. Journal of Agricuhural Science Yanbian University. 2011;33(2): 133–136.
  52. 52. Yao HY, Jiao XD, Wu FZ. Effects of continuous cucumber cropping and alternative rotations under protected cultivation on soil microbial community diversity. Plant Soil. 2006;284: 195–203.
  53. 53. Lin XG, Yin R, Zhang HY, Huang JF, Chen RR, Cao ZH. Changes of soil microbiological properties caused by land use changing from rice-wheat rotation to vegetable cultivation. Environmental Geochemistry and Health. 2004;26: 119–128. pmid:15499767
  54. 54. Balota EL, Kanashiro M, Colozzi Filho A, Andrade DS, Dick RP. Soil enzyme activities under long-term tillage and crop rotation systems in subtropical agro-ecosystems. Brazilian Journal of Microbiology. 2004;35(4): 300–306.
  55. 55. Dick RP, Pankhurst C, Doube B M. Soil enzyme activities as integrative indicators of soil health. Biological indicators of soil health. 1997: 121–156.
  56. 56. Piotrowska A, Dlugosz J, Zamorski R, Bogdanowicz P. Changes of enzymatic activity in soil supplemented with microbiological preparation. UGmax®19th World Congress of Soil Science. 2010.
  57. 57. Trasar-Cepeda C, Leiros MC, Seoane S, Gil-Sotres F. Limitations of soil enzymes as indicators of soil pollution. Soil Biology and Biochemistry. 2000;32(13): 1867–1875.
  58. 58. Konsler TR, Zito SW, Shelton JE, Staba EJ. Lime and phosphorus effects on American ginseng: II. Root and leaf ginsenoside content and their relationship. Journal of the American Society for Horticultural Science. 1990;115(4): 575–580.
  59. 59. Broeckling CD, Broz AK, Bergelson J, Manter DK, Vivanco JM. Root exudates regulate soil fungal community composition and diversity. Applied and Environmental Microbiology. 2008;74(3): 738–744. pmid:18083870
  60. 60. Yu KW, Gao WY, Hahn EJ, Paek KY. Effects of macro elements and nitrogen source on adventitious root growth and ginsenoside production in ginseng (Panax ginseng CA Meyer). Journal of Plant Biology. 2001;44(4): 179–184.
  61. 61. Murphy L. Effects of American ginseng on breast cancer and prostate cancer cells. American Ginseng; Production in the 21st Century; Conference Proceedings 2000; Cornell Cooperative Extension of Green County. 2000: 39–45.
  62. 62. Attele AS, Zhou YP. Antidiabetic effects of Panax ginseng berry extract and the identification of an effective component. Diabetes. 2002;51(6): 1851–1858. pmid:12031973