Skip to main content

Metabolome and transcriptome analysis reveals the molecular profiles underlying the ginseng response to rusty root symptoms



Ginseng rusty root symptoms (GRS) is one of the primary diseases of ginseng. This disease leads to a severe decline in the quality of ginseng. It has been shown that the occurrence of GRS is associated with soil environmental degradation, which may involve changes in soil microbiology and physicochemical properties.


In this study, GRS and healthy ginseng (HG) samples were used as experimental materials for comparative analysis of transcriptome and metabolome. Compared with those in HG samples, 949 metabolites and 9451 genes were significantly changed at the metabolic and transcriptional levels in diseased samples. The diseased tissues’ metabolic patterns changed, and the accumulation of various organic acids, alkaloids, alcohols and phenols in diseased tissues increased significantly. There were significant differences in the expression of genes involved in plant hormone signal transduction, phenylpropanoid biosynthesis, the peroxidase pathway, and the plant-pathogen interaction pathway.


The current study involved a comparative metabolome and transcriptome analysis of GRS and HG samples. Based on the findings at the transcriptional and metabolic levels, a mechanism model of the ginseng response to GRS was established. Our results provide new insights into ginseng’s response to GRS, which will reveal the potential molecular mechanisms of this disease in ginseng.

Peer Review reports


Ginseng (Panax ginseng Mayer) is one of the most important foods and herbs because of its high nutritional value and medicinal potential. As a perennial plant, the active ingredients in ginseng accumulate over time [1]. Ginseng has higher requirements for the growing environment and a longer growing period, making it susceptible to various diseases during growth, thus affecting its yield and quality.

Ginseng rusty root symptoms (GRS) is one of the most common diseases in ginseng cultivation and production. It produces reddish-brown spots on the periderm of ginseng roots, and over time, the spots may gradually expand, leading to a decline in the commodity grade and quality of ginseng. Whether GRS is an infectious disease or a non-infectious physiological disease is still controversial, and some scholars in China have tried to answer this question [2, 3]. Previous research on American ginseng (Panax quinquefolius L.) showed that rusty roots might be a defence mechanism against the invasion of certain fungi, resulting in the stimulation of phenolic compound production [4, 5]. Some research results have indicated that fungal infection is the cause of GRS [6,7,8]. However, this disease’s effect on ginseng at the transcriptional and metabolic levels and the response mechanism of ginseng are still unclear.

RNA-seq is a method of transcriptome analysis using deep sequencing technology [9]. This method can detect almost all genes and pathways related to a physiological response with high sensitivity [10]. In recent years, increasing studies have used the Illumina RNA-eq platform based on transcriptome analysis to explore plant responses to abiotic or biological stresses and understand their associated molecular mechanisms [11,12,13,14]. Metabolomics, similar to transcriptomics, is an essential part of systems biology. Metabolomics aims to detect and quantify all metabolites in biological samples [15]. This method is gradually being applied in medicinal plant research and plant pathology as an essential bridge connecting genes, metabolites, and phenotypes. In the research field of plant pathology, this method has been gradually combined with other research methods and applied to research on plant disease diagnosis [16], plant disease resistance [17], plant noninfectious diseases [18, 19], and other fields showing good application prospects.

In this study, GRS and HG tissues from the same ginseng farm were selected as experimental materials. Transcriptome and metabolome analyses of GRS and HG tissues were performed using RNA-seq and ultra-high-performance liquid chromatography-tandem mass spectrometry (UHPLC–MS/MS) to elucidate the mechanism of the ginseng response to this disease. These findings provide not only valuable information for the prevention and treatment of GRS but also important insights for further understanding the molecular mechanism of the ginseng response to external stress.


Quality control of RNA-seq reads

We performed RNA-seq analysis to compare the transcriptome between GRS and HG samples, and mainly analysed mRNA accumulation information obtained from the cDNA libraries. The sequencing of each ginseng sample produced over 100 million reads. After filtering out adaptors and low-quality reads, we obtained high-quality clean reads, accounting for more than 96% of the raw reads (Table S1).

Differentially expressed genes (DEGs) of diseased tissues

The DEGs between the GRS and HG samples are shown in the volcano plot and clustering map (Fig. 1a-b). The results showed that compared with the HG group, the expression of 4570 genes in the GRS group was upregulated, and that of 4881 genes was downregulated.

Fig. 1
figure 1

Expression profiling changes of genes in diseased tissues. Volcano plot indicating upregulated and downregulated genes (a); heat map showing hierarchical clustering of DEGs (b). HG: healthy ginseng; GRS: ginseng rusty root symptoms

To find the gene functions of DEGs in ginseng rusty roots, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. GO is the international standard classification system for gene functions. All of the DEGs can be divided into three categories, including biological process (BP), cellular component (CC), and molecular function (MF).

GO enrichment histograms drawn after classification of the enriched GO terms are shown in Fig. 2a. Based on GO analysis of the transcriptionally upregulated mRNAs, the most significantly enriched BPs were “carbohydrate metabolic process,” “steroid biosynthetic process,” “oxidation-reduction process,” “cell wall organization or biogenesis,” and “response to oxidative stress,” and the most significantly enriched CCs were “cell wall” and “beta-galactosidase complex.” The most significantly enriched MFs were “hydrolase activity,” “oxidoreductase activity,” and “peroxidase activity.” In addition, the histogram drawn according to transcriptionally downregulated mRNAs is shown in Fig. 2b. The results show that the most significantly enriched BPs were “oxidation-reduction process” and “terpenoid biosynthetic process,” and the most significantly enriched CCs were “apoplast” and “cell wall.” The most significantly enriched MFs were “iron ion binding,” “hydrolase activity,” and “peroxidase activity.” All of the enrichment results are shown in Tables S2 and S3.

Fig. 2
figure 2

Functional GO and KEGG pathway classification of DEGs in diseased tissues. a GO analysis of upregulated genes; b GO analysis of downregulated genes; c DEGs enriched the KEGG pathway scatterplot

For KEGG analysis of DEGs, the most significantly involved pathways were “metabolic pathways,” “biosynthesis of secondary metabolites,” and “phenylpropanoid biosynthesis” (Fig. 2c). All enriched KEGG pathways containing DEGs are given in Table S4.

Differentially accumulated metabolites of diseased tissues

UHPLC–MS/MS determined the composition of metabolites in diseased and healthy tissues, and 1101 and 649 compounds were analysed in ginseng root tissues in the positive ion mode and negative ion mode, respectively (Tables S5-S6). The most abundant metabolites extracted included lipids and lipid-like molecules, organic acids and derivatives, alkaloids, and phenylpropanoids.

As shown in Fig. 3a-b for both the positive ion mode and negative ion mode, the metabolites for the two kinds of ginseng tissues were clustered together based on principal component analysis (PCA). This result suggests that there were significant differences in the metabolism of the two samples. We used the supervised partial least squares discriminant analysis (PLS-DA) method to separate the samples based on critical metabolite information. The PLS-DA score plot showed significant separation between the two groups of samples (Fig. 3c-d). Additionally, to evaluate the models’ reliability, we performed the permutation test with 200 iterations. The cross-validation results showed that neither the positive ion model nor the negative ion model was overfitted (R2 = 0.79 and Q2 = -0.90, R2 = 0.74 and Q2 = -0.82, respectively) (Fig. 3e-f).

Fig. 3
figure 3

Metabolomics profiling of diseased and healthy ginseng tissues. PCA plots of the positive ion mode (a) and negative ion mode (b); PLS-DA plots of the positive ion mode (c) and negative ion mode (d); permutation tests for PLS-DA models in the positive ion mode (e) and negative ion mode (f). HG: healthy ginseng; GRS: ginseng rusty root symptoms

We extracted VIP values from the variable importance plots of the PLS-DA models to search for differentially accumulated metabolites in diseased ginseng tissues. Differentially accumulated metabolites were screened based on three criterions: VIP > 1, P < 0.05 and fold change > 2. According to the above criteria, 595 and 344 metabolites were selected from the positive ion mode and negative ion mode, respectively. To better understand the chemical differences between diseased and healthy tissues, we plotted the heatmaps of 12 samples in positive and negative ion modes (Fig. 4a-b). Most of the differentially accumulated metabolites screened in positive ion mode had higher accumulation in diseased tissues (419 upregulated and 176 downregulated). The same phenomenon was observed in negative ion mode (285 upregulated and 59 downregulated). By matching these metabolites with those in the database, we finally identified 669 differentially accumulated metabolites, of which 508 were upregulated and 159 were downregulated (Table 1).

Fig. 4
figure 4

The hierarchical clustering analysis and heatmaps of differentially accumulated metabolites from diseased and healthy ginseng tissues identified in positive (a) and negative ion mode (b); differentially accumulated metabolites enriched the KEGG pathway scatterplot in the diseased tissues (c); changes in the accumulation of metabolites annotated to the KEGG metabolic pathway in ginseng rusty root symptoms tissues (d). HG: healthy ginseng; GRS: ginseng rusty root symptoms

Table 1 Identification results of differential metabolites

The differentially accumulated metabolites in the two tissues were mapped to the KEGG database. These metabolites were mainly enriched in multiple alkaloid biosynthesis pathways, such as “Biosynthesis of alkaloids derived from histidine and purine” and “Biosynthesis of alkaloids derived from ornithine, lysine and nicotinic acid.” In addition, the extremely significant differences in “Stilbenoid, diarylheptanoid and gingerol biosynthesis,” “Phenylpropanoid biosynthesis,” “Biosynthesis of plant hormones,” and “Purine metabolism” are also noteworthy (Fig. 4c). The abundances of differentially accumulated metabolites across 9 metabolic pathways, as well as the significantly enriched differentially accumulated metabolites in KEGG metabolic pathways and those with a compound number ≥ 5 are shown in Fig. 4d.

Response of ginseng to disease

The results of the GO and KEGG enrichment analyses showed the response of GRS tissues. Numerous DEGs in GRS tissues were involved in plant hormone signal transduction, lignin synthesis, plant-pathogen interactions, and oxidative stress compared with those in HG tissues.

We focused on the pathway “plant hormone signal transduction,” involving 135 DEGs (Fig. S1-S2). Almost all genes involved in ABA signal transduction, such as soluble ABA receptors PYL/RCAR, protein phosphatases type-2C (PP2Cs), and SNF1-related protein kinase (SnRK2.3, SnRK2.10), were upregulated in the GRS group (Fig. S1B). The JAR1 and coronatine insensitive 1 (COI1) genes involved in JA signalling were upregulated, and the jasmonate-zim domain (JAZ) gene was downregulated (Fig. S1C). Additionally, differential gene expression of mediator protein nonexpressor of pathogen-related gene 1 (NPR1), TGA factors, and pathogenesis-related gene PR-1 involved in SA signalling was found for GRS (Fig. S1D). The DEGs associated with plant hormone signal transduction are shown in Fig. S2. In addition, we found that methyl jasmonate (MeJA), abscisic acid (ABA), and salicylic acid (SA) were strongly induced in diseased tissues (Fig. S1A).

There were 124 DEGs and 9 differentially accumulated metabolites involved in the “Phenylpropanoid biosynthesis” pathway (Table S4 and Fig. 4d). The detailed locations of differentially accumulated metabolites and DEGs in the pathway are shown in Fig. 5. In GRS tissues, enzymes closely related to lignin and phenol synthesis, such as phenylalanine ammonia-lyase (PAL), O-methyltransferase 1 (OMT1), 4-coumarate: coenzyme A ligase 1 (4CL1), and cinnamoyl coenzyme A reductase (CCR), were detected with a high expression level.

Fig. 5
figure 5

DEGs and differentially accumulated metabolites related to phenylpropanoid biosynthesis. The different color indicates the DEGs and differentially metabolites in GRS compared with HG (Red frame: up-regulation; green frame: down-regulation). PAL, phenylalanine ammonia-lyase; C4H, cinnamate 4-hydroxylase; 4CL, 4-coumarate: CoA ligase; CCR, cinnamoyl-CoA reductase; CAD, cinnamyl alcohol dehydrogenase; OMT1, caffeic acid 3-Omethyltransferase 1; F5H, ferulate-5-hydroxylase

We analysed the “peroxisome” pathway, which is closely related to oxidative stress in plants. As shown in Fig. S3 and Fig. S4, we found that 59 DEGs were involved in the “peroxisome” pathway. Several PEX family genes, such as PEX14, PEX16, and PEX19, were downregulated in the GRS group. PXMP2 and MPV17 were involved in reactive oxygen species (ROS) metabolism and showed obvious downregulated expression. Additionally, catalase (CAT) and superoxide dismutase (SOD) in the antioxidant system were upregulated.

In our enrichment results, 96 DEGs were involved in the “plant-pathogen interaction” pathway (Table S4). The corresponding KEGG pathway is shown in Fig. S5. Several genes, including the respiratory burst oxidase homologue (RBOH) genes RBOHD and RBOHF, flagellin sensing 2 (FLS2), and WRKY family genes (WRKY33 and WRKY22), involved in pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) were differentially expressed. The RIN4, RPS2, and HSP genes related to effector-triggered immunity (ETI) were also differentially expressed. Furthermore, the defence-related genes NHO1 and PR1 were upregulated (Fig. S6).

Correlation analysis of differentially accumulated metabolites and DEGs

We found that many DEGs and differentially accumulated metabolites were enriched in the same KEGG metabolic pathways in both types of tissue. These pathways include “Phenylpropanoid biosynthesis,” “Phenylalanine metabolism,” “Stilbenoid, diarylheptanoid and gingerol biosynthesis,” and others (Fig. 6a). In addition, considering the role of the phenanthrene biosynthesis pathway in disease resistance and the stress response in plants, we constructed a correlation network graph of DEGs and differentially accumulated metabolites based on the Pearson correlation coefficient (PCC). As shown in Fig. 6b, a total of 36 DEGs and 8 differentially accumulated metabolites in the network were highly positively correlated (PCC > 0.8). It is speculated that these genes and metabolites may play a key role in the response of ginseng to stress.

Fig. 6
figure 6

a Statistical results of DEGs and metabolites enriched in the same KEGG pathway; b connection network between DEGs and differentially accumulated metabolites in phenylpropanoid biosynthesis pathway

Real-time quantitative polymerase chain reaction (qPCR) validation of mRNA expression

To verify the accuracy of the RNA-seq results and provide the basis for further research, some genes from the four plant resistance pathways were selected for qRT-PCR analysis. The expression level of these genes was demonstrated by qRT-PCR and transcriptome sequencing. The genes Pg_S1410.4 (UGT72E1), Pg_S1825.18 (CCR1), Pg_S2641.9 (CSD1), Pg_S7602.1 (MPV17), Pg_S4621.10 (FLS2), Pg_S1378.3 (HSP90.1), Pg_S3064.15 (RP-1), and Pg_S0213.26 (JAZ) were analysed by qPCR (Fig. 7a). The expression levels of the corresponding mRNAs obtained by RNA-seq are shown in Fig. 7b. All validation results fully proved the reliability and accuracy of the transcriptome sequencing data.

Fig. 7
figure 7

qRT-PCR validation of significant DEGs. a qRT-PCR, the data were presented as the mean ± SEM (n = 3); **p < 0.05; b RNA-seq. HG: healthyginseng; GRS: ginseng rusty root symptoms


Metabolites and total RNAs were extracted from the same variety of 5-year-old ginseng planted in the same ginseng farm (due to its higher rusty root index). The samples were prepared from healthy and diseased ginseng tissues. The purpose was to ensure that ginseng samples experienced the same field management and climate changes during the growth process [20]. We conducted UHPLC–MS/MS and RNA-seq analysis of GRS and HG tissues. The results showed that in the GRS sample, 939 metabolites were differentially accumulated, and 9451 genes were significantly differentially expressed. Based on previous studies and the results of this metabolomic and transcriptomic analysis, we propose a model to explain ginseng’s response to GRS (Fig. 8).

Fig. 8
figure 8

Response mechanism of ginseng to rusty root symptoms. “?” means not validated; The solid arrows indicate the direct action and the dotted arrows indicate the indirect action. Environmental stress induces oxidative stress, which activates the ETI system and stimulates the synthesis of SA, MeJA, and ABA in ginseng tissues. SA and MeJA induce the accumulation of phenolic, alkaloid, and alcohol metabolites. SA and MeJA also induce lipid and fatty acid metabolism and enhance membrane stability against stress. The production of antioxidant enzymes such as CAT and SOD and the peroxisome regulation affect ROS accumulation in tissues. Meanwhile, DEGs in PTI and ETI in ginseng epidermal tissues may promote PR proteins' production and hypersensitive responses against external stresses. The activation of the phenylpropanoid biosynthesis pathway enhances the production of phenolic metabolites on the one hand. It may increase the accumulation of lignin and promote the cell wall's reinforcement on the other hand. The reinforcement of the cell wall and the enhancement of membrane stability together constitute a physical defense in ginseng tissues

Hormone regulation in diseased tissues

Plants produce a variety of hormones during growth. ABA, SA, and JA play an essential role in mediating plant defence responses to pathogens and abiotic stresses [21, 22]. At the transcriptional level, the ABA receptor PYR/PYL/RCAR, negative regulator PP2C, and positive regulator SnRK2 together comprise the regulatory system [23]. In GRS tissues, the genes that regulate ABA signal transduction were differentially expressed, and ABF transcription factors were upregulated (Fig. S1B). Under JA-stimulated conditions, JA-Ile, synthesized by JAR1, binds to the receptor COI1, leading to the degradation of the repressor protein JAZ, which upregulates the expression level of JA target genes (Fig. S1C) [24, 25]. NPR1 is one of the key regulatory elements of SA-dependent PR gene activation and interacts with GTA factors to regulate PR gene expression [26, 27]. We found differential expression of genes that are by-products of NPR1 (BOP2 and NPR3) and multiple TGA transcription factors genes, which resulted in the upregulation of PR-1 gene expression (Fig. S1D). In addition, three plant hormones and hormone-related signalling molecules (ABA, SA, and MeJA) were detected by UHPLC-MS/MS, and all of them were significantly upregulated in diseased tissues (Fig. S1A). Therefore, the accumulation of hormones and the differential expression of genes related to hormone signal transduction in diseased tissues affected metabolic processes and eventually led to changes in growth patterns to withstand environmental stress.

Immune system in diseased tissues

Plants have two layers of defence against pathogens, PTI and ETI [28]. Here, compared to HG tissues, there were 96 DEGs associated with the plant-pathogen interaction pathway in GRS tissues. RBOH mediates ROS generation, which is strongly related to both PTI and ETI [29, 30]. RBOHD and RBOHF are considered to be critical components of plant defences [31]. FLS2 is a pattern recognition receptor (PRR) that can trigger innate immunity in plants and is also involved in endocytosis [32]. WRKY33 plays a vital role as a transcription factor in plant resistance to pathogens [33]. The downregulation of the above genes, along with NOA1 and MEKK1, may indicate that part of PTI is repressed during GRS.

Interestingly, ETI seems to be activated in GRS samples. RPS2 is a resistance (R) protein in plants, while RIN4 is an essential negative regulator of natural plant immunity and has a regulatory effect on RPS2 [34]. HSP90.1 is required for RPS2 resistance and is rapidly induced in plants in the face of biological stress [35]. SGT1 positively regulates disease resistance produced by many R proteins [36]. In our analysis, RIN4 was downregulated, and PRS2, SGT1, and HSP90.1 were upregulated. Furthermore, NHO1, an essential universal resistance gene in plants, was upregulated [37]. Overall, the immune system of diseased tissues was activated, which indirectly led to hypersensitivity and cell wall reinforcement.

ROS-scavenging system in diseased tissues

Exposure of plants to various biotic and abiotic stress conditions triggers rapid changes in reactive oxygen species (ROS) production and clearance [38]. ROS play an essential role in signalling pathways that regulate inflammatory and defence responses in plants, but their accumulation is often harmful to cells.

Peroxisomes are subcellular organelles with a basic oxidative type of metabolism and may be the primary intracellular ROS production site [39, 40]. KEGG enrichment analysis showed that both PXMP2 and MPV17, which are involved in ROS metabolism, were transcriptionally downregulated (Fig. S2). PXMP2 is a perovskite membrane protein with a channel-forming activity that allows free diffusion of compounds, such as H2O2 [41]. MPV17 is a protein involved in the production of ROS [40]. Several PEX family protein genes were differentially expressed, including PEX13 and PEX14, which are involved in matrix protein import, and PEX16 and PEX19, which are involved in membrane protein import [42, 43]. In addition, CAT and SOD genes were upregulated in diseased tissues. Thus, diseased tissues may respond to oxidative stress through peroxisomes and antioxidant systems [44].

The phenolic content is often considered to be closely related to plant tissues’ total antioxidant capacity [45]. Studies have shown that flavonoids are also one of the main chemical defence substances in plants and can reduce various forms of reactive oxygen species in plant cells [46]. In our results, a significant upregulation of a large number of flavonoids and phenolics was found. Significant enrichment of the flavonoid biosynthetic pathway was also found by KEGG enrichment analysis (Fig. 4c). We also noted that multiple alcohol metabolites showed higher accumulation in diseased tissues. It was shown that the accumulation of alcohol usually plays a role in the fight against oxidative stress. Alcohols can act as intermediates of redox reactions with good scavenging capacity for free radicals and superoxide [47, 48]. It can be speculated that ginseng accumulates more alcohols in diseased tissues to cope with oxidative stress to improve its resistance. Overall, these compounds accumulate in diseased tissues of ginseng and may be jointly involved in the oxidative stress response.

Physical and chemical defence in diseased tissues

Lignin is a highly branched polymer of phenylpropanoid compounds and is the main component of the plant cell wall [49]. Lignin biosynthesis involves a complex genetic network containing many containing enzymes [50, 51]. Several key enzymes, including PAL, 4CL, CCR, OMT, and CAD, were upregulated [52, 53]. In addition, phenols and flavonoids not only play a key role in the plant antioxidant process, but are also essential compounds involved in plant chemical defences [54]. Our analysis showed that both types of data (transcriptomic and metabolomic) were significantly enriched in the phenylpropanoid biosynthesis pathway (Fig. 6a). Considering that the phenylpropanoid biosynthesis pathway is closely related to phenolic production, flavonoid metabolism, and lignin formation, it is speculated that it may be a critical metabolic pathway in ginseng in response to GRS [55]. Therefore, we constructed a correlation network of differentially accumulated metabolites and DEGs in this metabolic pathway based on the correlation analysis results (Fig. 6b). The network includes several upstream regulatory genes of phenolic, lignin, and flavonoid synthesis, such as PAL (Pg_S1239.2, Pg_S0954.1, and Pg_S2256.14), 4CL (Pg_S1367.7), and C4H (Pg_S4060.1), and precursor compounds, such as coniferin, sinapic acid, and coumarin. It is speculated that the upregulation of these genes and the accumulation of metabolites play a crucial role in the ginseng response to GRS.

Lipids are essential membrane components, and changes in their composition may help plants maintain cellular compartmentalization and membrane integrity [56]. In addition, free radicals are metabolized and produced when cells are under stress or in the process of injury. Therefore, reducing lipid peroxidation has an essential role in plant resilience. In diseased tissues of ginseng, there were more lipid metabolites with significantly altered contents. For example, phosphatidylcholine (PE) and phosphatidylethanolamine (PC) are lipids that are present in large amounts in cell membranes. An increase in the PC/PE ratio has been reported to be a marker of the membrane bilayer’s structural stability to prevent membrane degradation [57]. Interestingly, we detected a significant increase in several PC levels and a substantial decrease in PE levels in diseased tissues of ginseng. Therefore, diseased tissues may have a stronger tendency to maintain biofilm integrity and stability than healthy tissues.

The accumulation of alkaloids is a vital defence strategy for plants in response to biotic stresses, and alkaloids are a crucial class of compounds involved in plant chemical defences [58]. In diseased tissues, we found an upregulation of the accumulation of numerous alkaloid metabolites. The KEGG enrichment results of differentially accumulated metabolites also contained numerous alkaloid biosynthetic pathways (Fig. 2c). This result suggests that alkaloid metabolites are essential for the stress response of ginseng. Chemical defences against external stress involving alkaloids may occur in diseased tissues [55].


The current study performed a comparative metabolome and transcriptome analysis of GRS and HG tissues. An analysis of the pathways associated with the plant resistance response revealed that hormone signalling pathways, lignin biosynthesis, ROS regulation by peroxisomes changes during the plant-pathogen response. Additionally, the diseased tissues’ metabolic patterns changed, and the production and accumulation of many secondary metabolites may play an essential role in GRS. Finally, a model was proposed to explain the response of ginseng to GRS (Fig. 8). Our results provide new insights into ginseng’s response to GRS, revealing the potential molecular mechanisms of this disease in ginseng.


Plant material and tissue collection

All ginseng samples (5 years old) were collected from the same ginseng farm in Hunchun city, Jilin Province, China (42.86′N and 130.37′E). In the previous survey, this was a ginseng farm with a high rusty root index [59].

GRS diseased tissues (main root) and HG healthy tissues were collected as samples frozen in liquid nitrogen immediately and stored at − 80 °C before use. Independent biological replicates were prepared, and each replicate included root materials from three or more ginseng plants.

RNA extraction, library preparation, clustering and sequencing

Six complementary DNA (cDNA) libraries, three for GRS tissues, and three for HG tissues, were constructed. TRIzol reagent (Invitrogen, Carlsbad, CA) was used to isolate each sample’s total RNA. RNA degradation and contamination were monitored on 1% agarose gels. RNA purity was checked using a NanoPhotometer® spectrophotometer (IMPLEN, CA, USA). RNA concentration was measured using Qubit® RNA Assay Kit in a Qubit® 2.0 Fluorometer (Life Technologies, CA, USA). RNA integrity was assessed using an RNA Nano 6000 Assay Kit and the Bioanalyzer 2100 system (Agilent Technologies, CA, USA).

A total amount of 3 μg of RNA per sample was used as the input material for RNA sample preparation. Sequencing libraries were generated using an NEBNext® Ultra™ RNA Library Prep Kit for Illumina® (NEB, USA), following the manufacturer’s recommendations, and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from the total RNA using poly-T oligo-attached magnetic beads. Fragmentation was conducted using divalent cations under elevated temperatures in NEBNext First Strand Synthesis Reaction Buffer (5X). First strand cDNA was synthesized using a random hexamer primer and M-MuLV Reverse Transcriptase (RNase H-). Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. The remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of the 3′ ends of DNA fragments, the NEBNext adaptor with a hairpin loop structure was ligated to prepare for hybridization. To preferentially select cDNA fragments that were 150–200 bp in length, the library fragments were purified with the AMPure XP system (Beckman Coulter, Beverly, USA). Then, 3 μl of USER Enzyme (NEB, USA) was used with size-selected, adaptor-ligated cDNA at 37 °C for 15 min followed by 5 min at 95 °C before PCR. PCR was performed with Phusion High-Fidelity DNA polymerase, universal PCR primers, and an index (X) primer. PCR products were purified (AMPure XP system), and library quality was assessed on the Agilent Bioanalyzer 2100 system.

After cluster generation, RNA library preparations were sequenced on an Illumina HiSeq PE150 platform, and 150 bp paired-end reads were generated. The raw RNA-seq data are freely available in the NCBI database under accession no. PRJNA684799.

Transcriptome data analysis

Raw data (raw reads) in fastq format were first processed through in-house Perl scripts. A certain length range was then chosen from clean reads to perform all of the downstream analyses.

The reference genome and gene model annotation files were downloaded from the genome website directly [60]. The reference genome index was built using STAR, and paired-end clean reads were aligned to the reference genome using STAR (v.2.5.1b). STAR used the maximal mappable prefix (MMP) method to generate a precise mapping result for junction reads. HTSeq (v.0.6.0) was used to count the number of reads mapped to each gene. The expected number of fragments per kilobase of transcripts per million mapped reads (FPKM) of each gene was calculated based on the length and mapped read count for that gene.

Differential expression analysis of the two groups was performed using the DESeq2 R package (v.1.10.1). The resulting P values were adjusted using Benjamini and Hochberg’s approach for controlling the false discovery rate.

GO ( enrichment analysis of DEGs was implemented by the GOseq-based Wallenius noncentral hypergeometric distribution, in which gene length bias was corrected [61]. In addition, we performed KEGG enrichment analysis on DEGs using KOBAS 2.0 software [62].

Metabolome analysis

First, we extracted the metabolites. Tissues (100 mg) were individually ground with liquid nitrogen, and the homogenate was resuspended in prechilled 80% methanol and 0.1% formic acid by vortexing. The samples were incubated on ice for 5 min and then were centrifuged at 15,000 rpm and 4 °C for 5 min. Some supernatants were diluted with UHPLC–MS/MS grade water to a final concentration of 53% methanol. The samples were subsequently transferred to a fresh Eppendorf tube and then centrifuged at 15000×g and 4 °C for 10 min.

UHPLC-MS/MS analysis was then performed using a Vanquish UHPLC system (Thermo Fisher, Germany) coupled with an Orbitrap Q Exactive™ HF-X mass spectrometer (Thermo Fisher, Germany). Samples were injected onto a Hypesil Gold column (100 × 2.1 mm, 1.9 μm) using a 17-min linear gradient at a flow rate of 0.2 mL/min. The eluents for the positive polarity mode were eluent A (0.1% FA in water) and eluent B (methanol). The eluents for the negative polarity mode were eluent A (5 mM ammonium acetate, pH 9.0) and eluent B (methanol). The solvent gradient was set as follows: 2% B, 1.5 min; 2–100% B, 12.0 min; 100% B, 14.0 min; 100–2% B, 14.1 min; 2% B, 17 min. A Q Exactive™ HF-X mass spectrometer was operated in positive/negative polarity mode with a spray voltage of 3.2 kV, capillary temperature of 320 °C, sheath gas flow rate of 40 arb, and aux gas flow rate of 10 arb.

The raw data files generated by UHPLC-MS/MS were processed using the Compound Discoverer 3.1 (CD3.1, Thermo Fisher) to perform peak alignment, peak selection, and quantitation for each metabolite. The main parameters were set as follows: retention time tolerance, 0.2 min; actual mass tolerance, 5 ppm; signal intensity tolerance, 30%; signal/noise ratio, 3; and minimum intensity, 100,000. After that, peak intensities were normalized to the total spectral intensity. The normalized data were used to predict the molecular formula based on additive ions, molecular ion peaks, and fragment ions. Then, peaks were matched with the mzCloud, mzVault, and MassList databases to obtain accurate qualitative and relative quantitative results.

PCA and PLS-DA were performed using SIMCA-P 14.1 (Umetrics, Umea, Sweden). The metabolites with a VIP > 1, P-value < 0.05 and fold change > 2 were considered to be differential metabolites. KEGG enrichment analysis of differentially accumulated metabolites was performed using KOBAS 2.0 software [62].

Combined analysis of DEGs and differentially accumulated metabolites

The obtained differentially accumulated metabolites and DEGs were mapped to KEGG pathway for integration. In addition, we performed correlation analysis on transcriptomic and metabolomic data to explore correlations between DEGs and metabolites. The PCC of DEGs and differentially accumulated metabolites was calculated using the COR function in R. The genes and metabolites with a PCC > 0.8 in the KEGG pathway were used to establish the related network that was visualized by Cytoscape software.

Quantitative real-time qPCR

Quantitative PCR was performed on a Real-Time PCR System (Lightcycler 96, Roche, Switzerland) using 2 × RealStar Green Fast Mixture (GenStar, China). For each reaction, 0.5 μl of the forward and reverse primers and 2 μl of the cDNA template were added. All of the primers used in this study are listed in Table S7. The relative gene expression level was calculated according to the 2-ΔΔCt method [63].

Statistical analysis

Statistical analysis was performed with SPSS 20.0 software. All of the data were expressed as the mean ± SEM. P < 0.05 was considered statistically significant. RT-PCR results were analysed using a Student’s t-test. Significance was established at P < 0.05.

Availability of data and materials

All relevant supporting data sets are included in the article and its supplemental files.



Ginseng rusty root symptoms


Healthy ginseng


Liquid chromatography-tandem mass spectrometry


Kyoto Encyclopedia of Genes and Genomes


Differentially expressed genes


Gene Ontology


Biological process


Cellular component


Molecular function


Principal component analysis


Partial least-squares discriminant analysis


Methyl jasmonate


Abscisic acid


Salicylic acid


Protein phosphatases type-2C


SNF1-Related Protein Kinase


Coronatine insensitive 1


Jasmonate-zim domain


Protein non-expressor pr gene 1


Phenylalanine ammonia-lyase




4-coumarate: coenzyme A ligase


Cinnamoyl coenzyme A reductase


Reactive oxygen species




Superoxide dismutase


Respiratory burst oxidase homologues


Flagellin sensing 2


Pathogen-associated molecular patterns


PAMP-triggered immunity


Effector-triggered immunity


Pearson correlation coefficient


Quantitative polymerase chain reaction


Complementary DNA


Maximal mappable prefix


Fragments per kilobase of transcript per million mapped reads


  1. 1.

    Proctor JTA, Bailey WG. Ginseng: industry, botany, and culture. Hortic Rev. 1987;9:187–236.

    Google Scholar 

  2. 2.

    Zhou Y, Yang Z, Gao L, Liu W, Liu R, Zhao J, et al. Changes in element accumulation, phenolic metabolism, and antioxidative enzyme activities in the red-skin roots of Panax ginseng. J Ginseng Res. 2016;41(3):307–15.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Wang Q, Sun H, Xu C, Ma L, Li M, Shao C, et al. Analysis of rhizosphere bacterial and fungal communities associated with rusty root disease of Panax ginseng. Appl Soil Ecol. 2019;138:245–52.

    Article  Google Scholar 

  4. 4.

    Rahman M, Punja ZK. Biochemistry of ginseng root tissues affected by rusty root symptoms. Plant Physiol Biochem. 2005;43(12):1103–14.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Reeleder RD, Hoke SMT, Zhang Y. Rusted root of ginseng (Panax quinquefolius) is caused by a species of Rhexocercosporidium. Phytopathology. 2006;96(11):1243–54.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Lu XH, Zhang XM, Jiao XL, Hao JJ, Zhang XS, Luo Y, et al. Taxonomy of fungal complex causing red-skin root of Panax ginseng in China. J Ginseng Res. 2020;44(3):506–18.

    Article  PubMed  Google Scholar 

  7. 7.

    Farh ME, Kim YJ, Kim YJ, Yang DC. Cylindrocarpon destructans/Ilyonectria radicicola-species complex: causative agent of ginseng root-rot disease and rusty symptoms. J Ginseng Res. 2018;42(1):9–15.

    Article  PubMed  Google Scholar 

  8. 8.

    Liu D, Sun H, Ma H. Deciphering microbiome related to rusty roots of Panax ginseng and evaluation of antagonists against pathogenic Ilyonectria. Front Microbiol. 2019;10:1350.

    Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Wang Z, Gerstein M, Snyder MJNRG. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57–63.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Sun L, Wang J, Song K, Sun Y, Qin Q, Xue YJSR. Transcriptome analysis of rice (Oryza sativa L.) shoots responsive to cadmium stress. Sci Rep. 2019;9(1):1–10.

    Google Scholar 

  11. 11.

    Zhang Y, Gao X, Li J, Gong X, Yang P, Gao J, et al. Comparative analysis of proso millet ( Panicum miliaceum L.) leaf transcriptomes for insight into drought tolerance mechanisms. BMC Plant Biol. 2019;19(1):1–17.

    Article  Google Scholar 

  12. 12.

    Li H, Mo Y, Cui Q, Yang X, Guo Y, Wei C, et al. Transcriptomic and physiological analyses reveal drought adaptation strategies in drought-tolerant and -susceptible watermelon genotypes. Plant Sci. 2019;278:32–43.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Kang X, Wang L, Guo Y, Arifeen MZU, Cai X, Xue Y, et al. A comparative Transcriptomic and proteomic analysis of Hexaploid Wheat's responses to colonization by bacillus velezensis and Gaeumannomyces graminis, both separately and combined. Mol Plant Microbe. 2019;32(10):1336–47.

    CAS  Article  Google Scholar 

  14. 14.

    Bano A, Muqarab RJPB. Plant defence induced by PGPR against Spodoptera litura in tomato (Solanum lycopersicum L.). Plant Biol. 2017;19(3):406–12.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Guo J, Wu Y, Wang G, Wang T, Cao F. Integrated analysis of the transcriptome and metabolome in young and mature leaves of Ginkgo biloba L. Ind Crop Prod. 2020;143:111906.

    CAS  Article  Google Scholar 

  16. 16.

    Cevallos-Cevallos JM, García-Torres R, Etxeberria E, Reyes-De-Corcuera JI. GC-MS analysis of headspace and liquid extracts for Metabolomic differentiation of citrus Huanglongbing and zinc deficiency in leaves of ‘Valencia’ sweet Orange from commercial groves. Phytochem Anal. 2011;22(3):236–46.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Hamzehzarghani H, Paranidharan V, Abu-Nada Y, Kushalappa AC, Mamer O, Somers D. Metabolic profiling to discriminate wheat near isogenic lines, with quantitative trait loci at chromosome 2DL, varying in resistance to fusarium head blight. Can J Plant Sci. 2008;88(4):789–97.

    Article  Google Scholar 

  18. 18.

    Barding GA, Béni S, Fukao T, Bailey-Serres J, Larive CK. Comparison of GC-MS and NMR for metabolite profiling of Rice subjected to submergence stress. J Proteome Res. 2013;12(2):898–909.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Jahangir M, Abdel-Farid IB, Choi YH, Verpoorte R. Metal ion-inducing metabolite accumulation in Brassica rapa. J Plant Physiol. 2008;165(14):1429–37.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Jayakodi M, Lee S, Lee YS, Park H, Kim N, Jang W, et al. Comprehensive analysis of Panax ginseng root transcriptomes. BMC Plant Biol. 2015;15(1):1–12.

    CAS  Article  Google Scholar 

  21. 21.

    Jones JD, Dangl JL. The plant immune system. Nature. 2006;444(7117):323–9.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Bari R, Jones JD. Role of plant hormones in plant defence responses. Plant Mol Biol. 2009;69(4):473–88.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Fujii H, Chinnusamy V, Rodrigues A, Rubio S, Antoni R, Park SY, et al. In vitro reconstitution of an abscisic acid signalling pathway. Nature. 2009;462(7273):660–4.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Chini A, Fonseca S, Fernandez G, Adie B, Chico J, Lorenzo O, et al. The JAZ family of repressors is the missing link in jasmonate signalling. Nature. 2007;448(7154):666–71.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Chico JM, Chini A, Fonseca S, Solano R. JAZ repressors set the rhythm in jasmonate signaling. Curr Opin Plant Biol. 2008;11(5):486–94.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Canet JV, Dobon A, Roig A, Tornero PJPC. Environment. Structure-function analysis of npr1 alleles in Arabidopsis reveals a role for its paralogs in the perception of salicylic acid. Plant Cell Environ. 2010;33(11):1911–22.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Johnson C, Boden E, Arias JJTPC. Salicylic acid and NPR1 induce the recruitment of trans- activating TGA factors to a defense gene promoter in Arabidopsis. Plant Cell. 2003;15(8):1846–58.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Irieda H, Inoue Y, Mori M, Yamada K, Oshikawa Y, Saitoh H, et al. Conserved fungal effector suppresses PAMP-triggered immunity by targeting plant immune kinases. P Natl A Sci. 2019;116(2):496–505.

    CAS  Article  Google Scholar 

  29. 29.

    Marino D, Dunand C, Puppo A. Pauly NJTiPS. A burst of plant NADPH oxidases. Trends Plant Sci. 2012;17(1):9–15.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Zhang J, Shao F, Li Y, Cui H, Chen L, Li H, et al. A pseudomonas syringae effector inactivates MAPKs to suppress PAMP-induced immunity in plants. Cell Host Microbe. 2007;1(3):175–85.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Torres MA, Dangl JL, Jones JDG. Arabidopsis gp91phox homologues AtrbohD and AtrbohF are required for accumulation of reactive oxygen intermediates in the plant defense response. P Natl A Sci. 2002;99(1):517–22.

    CAS  Article  Google Scholar 

  32. 32.

    Robatzek S, Chinchilla D, Boller TJG. Development. Ligand-induced endocytosis of the pattern recognition receptor FLS2 in Arabidopsis. Genes Dev. 2006;20(5):537–42.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Zheng Z, Qamar SA, Chen Z, Mengiste TJPJ. Arabidopsis WRKY33 transcription factor is required for resistance to necrotrophic fungal pathogens. Plant J. 2006;48(4):592–605.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Liu J, Elmore JM, Coaker GJPS. Behavior. Investigating the functions of the RIN4 protein complex during plant innate immune responses. Plant Signal Behav. 2009;4(12):1107–10.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Takahashi A, Casais C, Ichimura K, Shirasu K. HSP90 interacts with RAR1 and SGT1 and is essential for RPS2-mediated disease resistance in Arabidopsis. P Natl A Sci. 2003;100(20):11777–82.

    CAS  Article  Google Scholar 

  36. 36.

    Azevedo C, Betsuyaku S, Peart JR, Takahashi A, Noel LD, Sadanandom A, et al. Role of SGT1 in resistance protein accumulation in plant immunity. EMBO J. 2006;25(9):2007–16.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Kang L, Li J, Zhao T, Xiao F, Tang X, Thilmony R, et al. Interplay of the Arabidopsis nonhost resistance gene NHO1 with bacterial virulence. P Natl A Sci. 2003;100(6):3519–24.

    CAS  Article  Google Scholar 

  38. 38.

    Czarnocka W, Karpinski S. Friend or foe? Reactive oxygen species production, scavenging and signaling in plant response to environmental stresses. Free Radical Bio Med. 2018;122:4–20.

    CAS  Article  Google Scholar 

  39. 39.

    Gill SS, Tuteja N. Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants. Plant Physiol Biochem. 2010;48(12):909–30.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Corpas FJ, Barroso JB, Rio LA. Peroxisomes as a source of reactive oxygen species and nitric oxide signal molecules in plant cells. Trends Plant Sci. 2001;6(4):145–50.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Lismont C, Koster J, Provost S, Baes M, Van Veldhoven PP, Waterham HR, et al. Deciphering the potential involvement of PXMP2 and PEX11B in hydrogen peroxide permeation across the peroxisomal membrane reveals a role for PEX11B in protein sorting. Biochim Biophys Acta Biomembr. 2019;1861(10):182991.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Lee MY, Sumpter R Jr, Zou Z, Sirasanagandla S, Wei Y, Mishra P, et al. Peroxisomal protein PEX13 functions in selective autophagy. EMBO Rep. 2017;18(1):48–60.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Cross LL, Ebeed HT, Baker A. Peroxisome biogenesis, protein targeting mechanisms and PEX gene functions in plants. Biochim Biophys Acta. 2016;1863(5):850–62.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Mittler R. Oxidative stress, antioxidants and stress tolerance. Trends Plant Sci. 2002;7(9):405–10.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Mwamba T, Islam F, Ali B, Lwalaba J, Gill R, Zhang F, et al. Comparative metabolomic responses of low-and high-cadmium accumulating genotypes reveal the cadmium adaptive mechanism in Brassica napus. Chemosphere. 2020;250:126308.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Chen J, Ullah C, Reichelt M, Gershenzon J, Hammerbacher A. Sclerotinia sclerotiorum circumvents flavonoid defenses by catabolizing Flavonol glycosides and Aglycones. Plant Physiol. 2019;180(4):1975–87.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  47. 47.

    Van den Ende W, Valluru R. Sucrose, sucrosyl oligosaccharides, and oxidative stress: scavenging and salvaging? J Exp Bot. 2008;60(1):9–18.

    CAS  Article  PubMed  Google Scholar 

  48. 48.

    KEUNEN E, PESHEV D, VANGRONSVELD J, VAN DEN ENDE W, CUYPERS A. Plant sugars are crucial players in the oxidative challenge during abiotic stress: extending the traditional concept. Plant Cell Environ. 2013;36(7):1242–55.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Moura JCMS, Bonine CAV, Viana JDOF, et al. Abiotic and biotic stresses and changes in the lignin content and composition in plants. J Integr Plant Biol. 2010;52(4):360–76.

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    Rogers LA, Campbell MM. The genetic control of lignin deposition during plant growth and development. New Phytol. 2004;164(1):17–30.

    CAS  Article  PubMed  Google Scholar 

  51. 51.

    Rastogi S, Dwivedi UN. Manipulation of lignin in plants with special reference to O-methyltransferase. Plnat Sci. 2008;174(3):264–77.

    CAS  Article  Google Scholar 

  52. 52.

    Lauvergeat V, Lacomme C, Lacombe E, Lasserre E, Roby D, Grima-Pettenati J. Two cinnamoyl-CoA reductase (CCR) genes from Arabidopsis thaliana are differentially expressed during development and in response to infection with pathogenic bacteria. Phytochemistry. 2001;57(7):1187–95.

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    Anterola AM, Lewis NG. Trends in lignin modification: a comprehensive analysis of the effects of genetic manipulations/mutations on lignification and vascular integrity. Phytochemistry. 2002;61(3):221–94.

    CAS  Article  PubMed  Google Scholar 

  54. 54.

    Gomaa NH, Hassan MO, Fahmy GM, González L, Hammouda O, Atteya AM. Flavonoid profiling and nodulation of some legumes in response to the allelopathic stress of Sonchus oleraceus L. Acta Bot Bras. 2015;29(4):553–60.

    Article  Google Scholar 

  55. 55.

    Zaynab M, Fatima M, Abbas S, Sharif Y, Umair M, Zafar MH, et al. Role of secondary metabolites in plant defense against pathogens. Microb Pathog. 2018;124:198–202.

  56. 56.

    Gigon A, Matos A-R, Laffray D, et al. Effect of drought stress on lipid metabolism in the leaves of Arabidopsis thaliana (ecotype Columbia). Ann Bot-London. 2004;94(3):345–51.

    CAS  Article  Google Scholar 

  57. 57.

    Toumi I, Gargouri M, Nouairi I, Moschou P, Salem-Fnayou AB, Mliki A, et al. Water stress induced changes in the leaf lipid composition of four grapevine genotypes with different drought tolerance. Biol Plant. 2008;52(1):161–4.

    CAS  Article  Google Scholar 

  58. 58.

    Mithöfer A, Maffei ME. General mechanisms of plant defense and plant toxins; 2016. p. 1–22.

    Google Scholar 

  59. 59.

    Bian X, Xiao S, Zhao Y, Xu Y, Yang H, Zhang L. Comparative analysis of rhizosphere soil physiochemical characteristics and microbial communities between rusty and healthy ginseng root. Sci Rep. 2020;10(1):15756.

  60. 60.

    Jayakodi M, Choi B-S, Lee S-C, Kim N-H, Park JY, Jang W, et al. Ginseng genome database: an open-access platform for genomics of Panax ginseng. BMC Plant Biol. 2018;18(1):62.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  61. 61.

    Young MD, Wakefield MJ, Smyth GK, Oshlack AJGB. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2):1–12.

    Article  Google Scholar 

  62. 62.

    Mao X, Cai T, Olyarchuk JG, Wei LJB. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics. 2005;21(19):3787–93.

    CAS  Article  PubMed  Google Scholar 

  63. 63.

    Arocho A, Chen B, Ladanyi M, Pan QJDMP. Validation of the 2-DeltaDeltaCt calculation as an alternate method of data analysis for quantitative PCR of BCR-ABL P210 transcripts. Diagn Mol Pathol. 2006;15(1):56–61.

    CAS  Article  PubMed  Google Scholar 

Download references


Not applicable.


This study was supported in part by grants from the National key R&D program (grant No. 2017YFC1702100) and the Jilin Province Major Science and Technology Special Project (grant No. 20200504003YY). Funders helped design the experiment for this study.

Author information




L.Z. designed the experiment; X.B. and S.X. analyzed the samples; Y.Z., H.Y. and Y.H. analyzed the data; X.B. wrote the first draft of the manuscript. All authors have read and approved the final mauscript.

Corresponding authors

Correspondence to Shengyuan Xiao or Lianxue Zhang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Table S1.

Summary of RNA-seq data and mapped reads for mRNA.

Additional file 2: Table S2.

Gene Ontology results of all the enriched terms for the upregulated mRNAs in GRS tissues.

Additional file 3: Table S3.

Gene Ontology results of all the enriched terms for the downregulated mRNAs in GRS tissues.

Additional file 4: Table S4.

All pathways associated with DEGs based on KEGG enrichment analysis.

Additional file 5: Table S5.

Results of UHPLC–MS/MS analysis in positive ion mode.

Additional file 6: Table S6.

Results of UHPLC–MS/MS analysis in negative ion mode.

Additional file 7: Table S7.

Primer sequences of the genes for qPCR verification.

Additional file 8: Figure S1.

Plant hormone contents and heatmaps of DEGs.

Additional file 9: Figure S2.

The plant hormone signal transduction pathway.

Additional file 10: Figure S3.

Heatmap of DEGs related to peroxisome synthesis.

Additional file 11: Figure S4.

The peroxisome pathway.

Additional file 12: Figure S5.

The plant-pathogen interaction pathway.

Additional file 13: Figure S6.

Heatmap of DEGs related to plant-pathogen interaction.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, 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 changes were made. 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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Bian, X., Zhao, Y., Xiao, S. et al. Metabolome and transcriptome analysis reveals the molecular profiles underlying the ginseng response to rusty root symptoms. BMC Plant Biol 21, 215 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Panax ginseng
  • Rusty root
  • Metabolome
  • Transcriptome
  • Environmental stress